Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probability Distribution Function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

Charged particle dynamics in the presence of non-Gaussian Lévy electrostatic fluctuations

Full orbit dynamics of charged particles in a 3-dimensional helical magnetic field in the presence of -stable Levy electrostatic fluctuations and linear friction modeling collisional Coulomb drag is studied via Monte Carlo numerical simulations. The Levy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space resulting from intermittent electrostatic turbulence. The probability distribution functions of energy, particle displacements, and Larmor radii are computed and showed to exhibit a transition from exponential decay, in the case of Gaussian fluctuations, to power law decay in the case of Levy fluctuations. The absolute value of the power law decay exponents are linearly proportional to the Levy index. Furthermore, the observed anomalous non-Gaussian statistics of the particles' Larmor radii (resulting from outlier transport events) indicate that, when electrostatic turbulent fluctuations exhibit non-Gaussian Levy statistics, gyro-averaging and guiding centre approximations might face limitations and full particle orbit effects should be taken into account.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Molecular dynamics simulations of concentration-dependent defect production in Fe-Cr and Fe-Cu alloys

Molecular dynamics simulations are conducted to study the effects of alloying elements on the primary damage behaviors in three Fe-based ferritic alloy systems: (1) a Fe-Cr system in which the heat of mixing changes its sign with the Cr concentration; (2) a Fe-Cu system that has a positive heat of mixing; and (3) an ideal but artificial Fe-Cr system that has a zero heat of mixing, which is used as a reference system to investigate solute interstitial formation based on probability. It is found that in these alloys, the solute type and concentration do not have a significant effect on the total number of surviving Frenkel pairs. However, the fraction of solute interstitials has distinct behaviors. In Fe-Cr, the Cr interstitial fraction is much higher than the Cr solute concentration and the Cr interstitial production efficiency decreases with the increasing Cr concentration. By contrast, in Fe-Cu, Cu interstitials are barely produced. In the ideal alloy, the solute interstitial fraction is close to the solute concentration. The defect formation energies in both dilute and concentrated alloys, interstitial binding energies, liquid diffusivities of Fe and solute atoms, and heat of mixing have been calculated for both Fe-Cr and Fe-Cu alloys. Among them, we find that the relative thermodynamic stability between Fe self-interstitials and solute interstitials plays the most important role in the solute interstitial production behaviors. The decrease of Cr interstitial production efficiency with increasing Cr concentration can be explained by the probability distribution functions of solute interstitial formation energy in concentrated alloys.

Zhang, Yaxuan↗

Propagation of input parameter uncertainties in transport models

Here, the many sources of uncertainty in validation studies of plasma turbulence in magnetically confined fusion devices are well-known. In this paper, we investigate how to efficiently transform uncertainties in experimentally derived transport model inputs into model prediction uncertainties, using the quasilinear trapped-gyro-Landau-fluid (TGLF) turbulent transport model. We use the rapidly converging and computationally inexpensive non-intrusive probabilistic collocation method (PCM) to propagate input parameter uncertainty probability distribution functions (PDFs) through TGLF, yielding PDFs of predicted transport fluxes. We observe in many cases that the flux PDFs exhibit significant non-normal features such as strong skewness, even when the input distributions were normal. To illustrate the utility of the PCM approach, we apply this methodology to transport predictions for a DIII-D ITER baseline plasma in which the mix of neutral beam injection (NBI) and electron cyclotron heating (ECH) was varied. The model predictions show clear changes in the parametric dependencies and sensitivities of the turbulence between the two heating mixes. Specifically, when only NBI heating was used, the transport fluxes responded significantly only to the ion temperature gradient scale length. However, when both NBI and ECH were applied, the electron transport channels demonstrate a strong sensitivity to the electron temperature and density gradients not observed in the NBI-only case. Additional context for the PCM approach is provided by comparing its predictions with those obtained via a local flux-matching approach. A new set of validation metrics based on the Wasserstein distance is proposed for PDF-based comparisons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A LAW FOR STAR FORMATION IN GALAXIES

We study the galactic-scale triggering of star formation. We find that the largest mass scale not stabilized by rotation, a well-defined quantity in a rotating system and with clear dynamical meaning, strongly correlates with the star formation rate in a wide range of galaxies. We find that this relation can be understood in terms of self-regulation toward marginal Toomre stability and the amount of turbulence allowed to sustain the system in this self-regulated quasi-stationary state. We test such an interpretation by computing the predicted star formation rates for a galactic interstellar medium characterized by a lognormal probability distribution function and find good agreement with the observed relation.

79 ASTRONOMY AND ASTROPHYSICS↗

INTERSTELLAR SONIC AND ALFVENIC MACH NUMBERS AND THE TSALLIS DISTRIBUTION

In an effort to characterize the Mach numbers of interstellar medium (ISM) magnetohydrodynamic (MHD) turbulence, we study the probability distribution functions (PDFs) of spatial increments of density, velocity, and magnetic field for 14 ideal isothermal MHD simulations at a resolution of 512{sup 3}. In particular, we fit the PDFs using the Tsallis function and study the dependency of the fit parameters on the compressibility and magnetization of the gas. We find that the Tsallis function fits PDFs of MHD turbulence well, with fit parameters showing sensitivities to the sonic and Alfven Mach numbers. For three-dimensional density, column density, and Position-Position-Velocity data, we find that the amplitude and width of the PDFs show a dependency on the sonic Mach number. We also find that the width of the PDF is sensitive to the global Alfvenic Mach number especially in cases where the sonic number is high. These dependencies are also found for mock observational cases, where cloud-like boundary conditions, smoothing, and noise are introduced. The ability of Tsallis statistics to characterize the sonic and Alfvenic Mach numbers of simulated ISM turbulence points to it being a useful tool in the analysis of the observed ISM, especially when used simultaneously with other statistical techniques.

79 ASTRONOMY AND ASTROPHYSICS↗

THE OCCURRENCE RATE OF EARTH ANALOG PLANETS ORBITING SUN-LIKE STARS

Kepler is a space telescope that searches Sun-like stars for planets. Its major goal is to determine {eta}{sub Earth}, the fraction of Sun-like stars that have planets like Earth. When a planet 'transits' or moves in front of a star, Kepler can measure the concomitant dimming of the starlight. From analysis of the first four months of those measurements for over 150,000 stars, Kepler's Science Team has determined sizes, surface temperatures, orbit sizes, and periods for over a thousand new planet candidates. In this paper, we characterize the period probability distribution function of the super-Earth and Neptune planet candidates with periods up to 132 days, and find three distinct period regimes. For candidates with periods below 3 days, the density increases sharply with increasing period; for periods between 3 and 30 days, the density rises more gradually with increasing period, and for periods longer than 30 days, the density drops gradually with increasing period. We estimate that 1%-3% of stars like the Sun are expected to have Earth analog planets, based on the Kepler data release of 2011 February. This estimate of {eta}{sub Earth} is based on extrapolation from a fiducial subsample of the Kepler planet candidates that we chose to be nominally 'complete' (i.e., no missed detections) to the realm of the Earth-like planets, by means of simple power-law models. The accuracy of the extrapolation will improve as more data from the Kepler mission are folded in. Accurate knowledge of {eta}{sub Earth} is essential for the planning of future missions that will image and take spectra of Earth-like planets. Our result that Earths are relatively scarce means that a substantial effort will be needed to identify suitable target stars prior to these future missions.

79 ASTRONOMY AND ASTROPHYSICS↗

MAGNETIZATION OF CLOUD CORES AND ENVELOPES AND OTHER OBSERVATIONAL CONSEQUENCES OF RECONNECTION DIFFUSION

Recent observational results for magnetic fields in molecular clouds reviewed by Crutcher seem to be inconsistent with the predictions of the ambipolar diffusion theory of star formation. These include the measured decrease in mass to flux ratio between envelopes and cores, the failure to detect any self-gravitating magnetically subcritical clouds, the determination of the flat probability distribution function (PDF) of the total magnetic field strengths implying that there are many clouds with very weak magnetic fields, and the observed scaling B{proportional_to}{rho}{sup 2/3} that implies gravitational contraction with weak magnetic fields. We consider the problem of magnetic field evolution in turbulent molecular clouds and discuss the process of magnetic field diffusion mediated by magnetic reconnection. For this process that we termed 'reconnection diffusion', we provide a simple physical model and explain that this process is inevitable in view of the present-day understanding of MHD turbulence. We address the issue of the expected magnetization of cores and envelopes in the process of star formation and show that reconnection diffusion provides an efficient removal of magnetic flux that depends only on the properties of MHD turbulence in the core and the envelope. We show that as the amplitude of turbulence as well as the scale of turbulent motions decrease from the envelope to the core of the cloud, the diffusion of the magnetic field is faster in the envelope. As a result, the magnetic flux trapped during the collapse in the envelope is being released faster than the flux trapped in the core, resulting in much weaker fields in envelopes than in cores, as observed. We provide simple semi-analytical model calculations which support this conclusion and qualitatively agree with the observational results. Magnetic reconnection is also consistent with the lack of subcritical self-gravitating clouds, with the observed flat PDF of field strengths, and with the scaling of field strength with density. In addition, we demonstrate that the reconnection diffusion process can account for the empirical Larson relations and list a few other implications of the reconnection diffusion concept. We argue that magnetic reconnection provides a solution to the magnetic flux problem of star formation that agrees better with observations than the long-standing ambipolar diffusion paradigm. Due to the illustrative nature of our simplified model we do not seek quantitative agreement, but discuss the complementary nature of our approach to the three-dimensional MHD numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

IMPROVING THE PRECISION OF TIME-DELAY COSMOGRAPHY WITH OBSERVATIONS OF GALAXIES ALONG THE LINE OF SIGHT

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line-of-sight convergence {kappa}{sub ext}. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with {kappa}{sub ext} obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of {kappa}{sub ext} from {sigma}{sub {kappa}} {approx} 0.06 to {approx}0.04 for very overdense LOSs like that of the system B1608+656. For more common LOSs, {sigma}{sub {kappa}} is reduced to {approx}<0.03, corresponding to an uncertainty of {approx}< 3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain {kappa}{sub ext} almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable {kappa}{sub ext} estimator gives {sigma}{sub {kappa}} = 0.047 down from 0.065 using only galaxy counts. Although deeper multiband observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

79 ASTRONOMY AND ASTROPHYSICS↗

MODELING THE POLLUTION OF PRISTINE GAS IN THE EARLY UNIVERSE

We conduct a comprehensive theoretical and numerical investigation of the pollution of pristine gas in turbulent flows, designed to provide useful new tools for modeling the evolution of the first generation of stars. The properties of such Population III (Pop III) stars are thought to be very different than those of later stellar generations, because cooling is dramatically different in gas with a metallicity below a critical value Z{sub c}, which lies between ∼10{sup –6} and ∼10{sup –3} Z{sub ☉}. The critical value is much smaller than the typical overall average metallicity, , and therefore the mixing efficiency of the pristine gas in the interstellar medium plays a crucial role in determining the transition from Pop III to normal star formation. The small critical value, Z{sub c}, corresponds to the far left tail of the probability distribution function (PDF) of the metal abundance. Based on closure models for the PDF formulation of turbulent mixing, we derive evolution equations for the fraction of gas, P, lying below Z{sub c}, in statistically homogeneous compressible turbulence. Our simulation data show that the evolution of the pristine fraction P can be well approximated by a generalized 'self-convolution' model, which predicts that P-dot = - (n/τ{sub con}) P (1-P{sup 1/n}), where n is a measure of the locality of the mixing or PDF convolution events and the convolution timescale τ{sub con} is determined by the rate at which turbulence stretches the pollutants. Carrying out a suite of numerical simulations with turbulent Mach numbers ranging from M = 0.9 to 6.2, we are able to provide accurate fits to n and τ{sub con} as a function of M, Z{sub c}/(Z), and the length scale, L{sub p}, at which pollutants are added to the flow. For pristine fractions above P = 0.9, mixing occurs only in the regions surrounding blobs of pollutants, such that n = 1. For smaller values of P, n is larger as the mixing process becomes more global. We show how these results can be used to construct one-zone models for the evolution of Pop III stars in a single high-redshift galaxy, as well as subgrid models for tracking the evolution of the first stars in large cosmological numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

TURBULENCE-INDUCED RELATIVE VELOCITY OF DUST PARTICLES. I. IDENTICAL PARTICLES

We study the relative velocity of inertial particles suspended in turbulent flows and discuss implications for dust particle collisions in protoplanetary disks. We simulate a weakly compressible turbulent flow, evolving 14 particle species with friction timescale, τ{sub p}, covering the entire range of scales in the flow. The particle Stokes numbers, St, measuring the ratio of τ{sub p} to the Kolmogorov timescale, are in the range 0.1 ∼< St ∼< 800. Using simulation results, we show that the model by Pan and Padoan gives satisfactory predictions for the rms relative velocity between identical particles. The probability distribution function (PDF) of the relative velocity is found to be highly non-Gaussian. The PDF tails are well described by a 4/3 stretched exponential function for particles with τ{sub p} ≅ 1-2 T{sub L}, where T{sub L} is the Lagrangian correlation timescale, consistent with a prediction based on PP10. The PDF approaches Gaussian only for very large particles with τ{sub p} ∼> 54 T{sub L}. We split particle pairs at given distances into two types with low and high relative speeds, referred to as continuous and caustic types, respectively, and compute their contributions to the collision kernel. Although amplified by the effect of clustering, the continuous contribution vanishes in the limit of infinitesimal particle distance, where the caustic contribution dominates. The caustic kernel per unit cross section rises rapidly as St increases toward ≅ 1, reaches a maximum at τ{sub p} ≅ 2 T{sub L}, and decreases as τ{sub p}{sup -1/2} for τ{sub p} >> T{sub L}.

79 ASTRONOMY AND ASTROPHYSICS↗

A large catalog of accurate distances to molecular clouds from PS1 photometry

Distance measurements to molecular clouds are important but are often made separately for each cloud of interest, employing very different data and techniques. We present a large, homogeneous catalog of distances to molecular clouds, most of which are of unprecedented accuracy. We determine distances using optical photometry of stars along lines of sight toward these clouds, obtained from PanSTARRS-1. We simultaneously infer the reddenings and distances to these stars, tracking the full probability distribution function using a technique presented in Green et al. We fit these star-by-star measurements using a simple dust screen model to find the distance to each cloud. We thus estimate the distances to almost all of the clouds in the Magnani et al. catalog, as well as many other well-studied clouds, including Orion, Perseus, Taurus, Cepheus, Polaris, California, and Monoceros R2, avoiding only the inner Galaxy. Typical statistical uncertainties in the distances are 5%, though the systematic uncertainty stemming from the quality of our stellar models is about 10%. The resulting catalog is the largest catalog of accurate, directly measured distances to molecular clouds. Our distance estimates are generally consistent with available distance estimates from the literature, though in some cases the literature estimates are off by a factor of more than two.

79 ASTRONOMY AND ASTROPHYSICS↗

A NOVEL TECHNIQUE TO IMPROVE PHOTOMETRY IN CONFUSED IMAGES USING GRAPHS AND BAYESIAN PRIORS

We present a new technique for overcoming confusion noise in deep far-infrared Herschel space telescope images making use of prior information from shorter λ < 2 μm wavelengths. For the deepest images obtained by Herschel, the flux limit due to source confusion is about a factor of three brighter than the flux limit due to instrumental noise and (smooth) sky background. We have investigated the possibility of de-confusing simulated Herschel PACS 160 μm images by using strong Bayesian priors on the positions and weak priors on the flux of sources. We find the blended sources and group them together and simultaneously fit their fluxes. We derive the posterior probability distribution function of fluxes subject to these priors through Monte Carlo Markov Chain (MCMC) sampling by fitting the image. Assuming we can predict the FIR flux of sources based on the ultraviolet-optical part of their SEDs to within an order of magnitude, the simulations show that we can obtain reliable fluxes and uncertainties at least a factor of three fainter than the confusion noise limit of 3σ {sub c} = 2.7 mJy in our simulated PACS-160 image. This technique could in principle be used to mitigate the effects of source confusion in any situation where one has prior information of positions and plausible fluxes of blended sources. For Herschel, application of this technique will improve our ability to constrain the dust content in normal galaxies at high redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

Rethinking a mysterious molecular cloud

I present high-resolution column density maps of two molecular clouds (MCs) having strikingly different star formation rates. To better understand the unusual, massive G216-2.5, an MC with no massive star formation, the distribution of its molecular gas is compared to that of the Rosette MC. Far-infrared data from Herschel are used to derive N(H{sub 2}) maps of each cloud and are combined with I{sub CO} data to determine the CO-to-H{sub 2} ratio, X{sub CO}. In addition, the probability distribution functions (PDFs) and cumulative mass fractions of the clouds are compared. For G216-2.5, 〈N(H{sub 2})〉=7.8×10{sup 20} cm{sup −2} and 〈X{sub CO}〉=2.2×10{sup 20} cm{sup −2} (K km s{sup −1}){sup −1}; for the Rosette, 〈N(H{sub 2})〉=1.8×10{sup 21} cm{sup −2} and 〈X{sub CO}〉=2.8×10{sup 20} cm{sup −2} (K km s{sup −1}){sup −1}. The PDFs of both clouds are log-normal for extinctions below ∼2 mag and both show departures from log-normality at high extinctions. Although it is the less-massive cloud, the Rosette has a higher fraction of its mass in the form of dense gas and contains 1389 M{sub ⊙} of gas above the so-called extinction threshold for star formation, A{sub V}=7.3 mag. The G216-2.5 cloud has 874 M{sub ⊙} of dense gas above this threshold.

79 ASTRONOMY AND ASTROPHYSICS↗

Lens statistics with gravitationally lensed yet morphologically regular images

GRAMORs are GRAvitationally highly magnified yet MORphologically regular images. An example of this phenomenon was discovered in the cluster MACS J1149.5+2223 in 2009. We investigate the lens statistics of GRAMORs in detail. Assuming a NFW profile for a sample of clusters, we calculate the expected number and redshift distribution of GRAMORs using parameters from COSMOS data for the number density of the background galaxy. A model with a cluster placed at z = 0.544 based on WMAP5 cosmology predicts the redshift of a GRAMOR at z≃1.49 which is close to the observed z = 1.4906. These results show that the expected number of GRAMORs is about two per cluster in the most likely case, and thus a large number of GRAMORs would be observed in a systematic survey. The probability distribution function of source redshift for GRAMORs depends strongly on dark energy and may be useful for constraining the nature of dark energy.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics of intrinsic rotation in flux-driven ITG turbulence

Global, heat flux-driven ITG gyrokinetic simulations which manifest the formation of macroscopic, mean toroidal flow profiles with peak thermal Mach number 0.05, are reported. Both a particle-in-cell (XGC1p) and a semi-Lagrangian (GYSELA) approach are utilized without a priori assumptions of scale separation between turbulence and mean fields. Flux-driven ITG simulations with different edge flow boundary conditions show in both approaches the development of net unidirectional intrinsic rotation in the co-current direction. Intrinsic torque is shown to scale approximately linearly with the inverse scale length of the ion temperature gradient. External momentum input is shown to effectively cancel the intrinsic rotation profile, thus confirming the existence of a local residual stress and intrinsic torque. Fluctuation intensity, intrinsic torque and mean flow are demonstrated to develop inwards from the boundary. The measured correlations between residual stress and two fluctuation spectrum symmetry breakers, namely E × B shear and intensity gradient, are similar. Avalanches of (positive) heat flux, which propagate either outwards or inwards, are correlated with avalanches of (negative) parallel momentum flux, so that outward transport of heat and inward transport of parallel momentum are correlated and mediated by avalanches. The probability distribution functions of the outward heat flux and the inward momentum flux show strong structural similarity.

Physics↗

Constraining dark matter microphysics with the annihilation signal from subhalos

In the cold dark matter scenario, galactic dark matter halos are populated with a large number of smaller subhalos. Previous work has shown that dark matter annihilations in subhalos can generate a distinctive, non-Poisson signal in the gamma-ray photon counts probability distribution function (PDF). Here we show that the gamma-ray PDF also carries information about the velocity dependence of the dark matter annihilation cross section. After calculating the PDF assuming s-wave and Sommerfeld-enhanced annihilation, we perform a mock data analysis to illustrate how current and future observations can constrain the microphysics of the dark matter annihilation. We find that, with current Fermi data, and assuming a dark matter annihilation cross section roughly at the limit of current bounds from annihilation in dwarf spheroidal galaxies, one can potentially distinguish the non-Poissonian fluctuations expected from dark matter annihilation in subhalos from Poisson sources, as well as from dark matter models with an incorrect velocity-dependence. Here we explore how robust these results are to assumptions about the modeling of the galactic gamma-ray background, but further work is needed to determine the impact of realistic astrophysical source populations on our results. We also point out a four-parameter degeneracy between the velocity dependence of the dark matter annihilation, the minimum subhalo mass, the power law index of the subhalo mass function, and the normalization of the dark matter signal. This degeneracy can be broken with priors from N-body simulations or from observational constraints on the subhalo mass function.

79 ASTRONOMY AND ASTROPHYSICS↗

CMB-PAInT: An inpainting tool for the cosmic microwave background

Abstract The presence of astrophysical emissions in microwave observations forces us to perform component separation to extract the Cosmic Microwave Background (CMB) signal. However, even in the most optimistic cases, there are still strongly contaminated regions, such as the Galactic plane or those with emission from extragalactic point sources, which require the use of a mask. Since many CMB analyses, especially the ones working in harmonic space, need the whole sky map, it is crucial to develop a reliable inpainting algorithm that replaces the values of the excluded pixels by others statistically compatible with the rest of the sky. This is especially important when working withQandUsky maps in order to obtainE- andB-mode maps which are free fromE-to-Bleakage. In this work we study a method based on Gaussian Constrained Realizations (GCR), that can deal with both intensity and polarization. Several tests have been performed to asses the validation of the method, including the study of the one-dimensional probability distribution function (1-PDF),E- andB-mode map reconstruction, and power spectra estimation. We have considered two scenarios for the input simulation: one case with only CMB signal and a second one including also Planck PR4 semi-realistic noise. Even if we are limited to low resolution maps, N side = 64 ifT,QandUare considered, we believe that this is a useful approach to be applied to future missions such as LiteBIRD, where the target are the largest scales.

Astronomy & Astrophysics↗

A Conditional Autoencoder for Galaxy Photometric Parameter Estimation

Astronomical photometric surveys routinely image billions of galaxies, and traditionally infer the parameters of a parametric model for each galaxy. This approach has served us well, but the computational expense of deriving a full posterior probability distribution function is a challenge for increasingly ambitious surveys. In this paper, we use deep learning methods to characterize galaxy images, training a conditional autoencoder on mock data. The autoencoder can reconstruct and denoise galaxy images via a latent space engineered to include semantically meaningful parameters, such as brightness, location, size, and shape. Our model recovers galaxy fluxes and shapes on mock data with a lower variance than the Hyper Suprime-Cam photometry pipeline, and returns reasonable answers even for inputs outside the range of its training data. When applied to data in the training range, the regression errors on all extracted parameters are nearly unbiased with a variance near the Cramr-Rao bound.

79 ASTRONOMY AND ASTROPHYSICS↗