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 127 records · Page 7

Lensing and the Warm-hot Intergalactic Medium

The correlation of weak lensing and Cosmic Microwave Anisotropy (CMB) data traces the pressure distribution of the hot, ionized gas and the underlying matter density field. The measured correlation is dominated by baryons residing in halos. Detecting the contribution from unbound gas by measuring the residual cross-correlation after masking all known halos requires a theoretical understanding of this correlation and its dependence with model parameters. Our model assumes that the gas in filaments is well described by a log-normal probability distribution function, with temperatures 10{sup 5–7} K and overdensities ξ ≤ 100. The lensing-Comptonization cross-correlation is dominated by gas with overdensities in the range ξ ≈ [3–33]; the signal is generated at redshifts z ≤ 1. If only 10% of the measured cross-correlation is due to unbound gas, then the most recent measurements set an upper limit of T-bar {sub e}≲10{sup 6} K on the mean temperature of the intergalactic Medium. The amplitude is proportional to the baryon fraction stored in filaments. The lensing-Comptonization power spectrum peaks at a different scale than the gas in halos, making it possible to distinguish both contributions. To trace the distribution of the low-density and low-temperature plasma on cosmological scales, the effect of halos will have to be subtracted from the data, requiring observations with larger signal-to-noise ratios than are currently available.

79 ASTRONOMY AND ASTROPHYSICS↗

Uncertainty in estimates of the number of extraterrestrial civilizations

Estimation of the number N of communicative civilizations by means of Drake's formula involves the combination of several quantities, each of which is to some extent uncertain. The uncertainty in any quantity may be represented by a probability distribution function, even if that quantity is itself a probability. The uncertainty of current estimates of N is derived principally from uncertainty in estimates of the lifetime of advanced civilizations. It is argued that this is due primarily to uncertainty concerning the existence of a "Galactic Federation" which is in turn contingent upon uncertainty about whether the limitations of present-day physics are absolute or (in the event that there exists a yet-undiscovered "hyperphysics") transient . It is further argued that it is advantageous to consider explicitly these underlying assumptions in order to compare the probable numbers of civilizations operating radio beacons, permitting radio leakage, dispatching probes for radio surveillance or dispatching vehicles for manned surveillance.

Peter A Sturrock↗

Using Geothermal Play Types as an Analogue for Estimating Potential Resource Size

Blind geothermal systems are becoming increasingly common as more geothermal fields are developed. Geothermal development is known to have high risk in the early stages of a project development because reservoir characteristics are relatively unknown until wells are drilled. Play types (or occurrence models) categorize potential geothermal fields into groups based on geologic characteristics. To aid in lowering exploration risk, these groups' reservoir characteristics can be used as analogues in new site exploration. The play type schemes used in this paper were Moeck and Beardsmore play types (Moeck et al. 2014) and Brophy occurrence models (Brophy et al. 2011). Operating geothermal fields throughout the world were classified based on their associated play type, and then reservoir characteristics data were catalogued. The distributions of these characteristics were plotted in histograms to develop probability density functions for each individual characteristic. The probability density functions can be used as input analogues in Monte Carlo estimations of resource potential for similar play types in early exploration phases. A spreadsheet model was created to estimate resource potential in undeveloped fields. The user can choose to input their own values for each reservoir characteristic or choose to use the probability distribution functions provided from the selected play type. This paper also addresses the United States Geological Survey's 1978 and 2008 assessment of geothermal resources by comparing their estimated values to reported values from post-site development. Information from the collected data was used in the comparison for thirty developed sites in the United States. No significant trends or suggestions for methodologies could be made by the comparison.

15 GEOTHERMAL ENERGY↗

A measurement of the turbulence-driven density distribution in a non-star-forming molecular cloud

Molecular clouds are supersonically turbulent. This turbulence governs the initial mass function and the star formation rate. In order to understand the details of star formation, it is therefore essential to understand the properties of turbulence, in particular the probability distribution of density in turbulent clouds. We present H{sub 2}CO volume density measurements of a non-star-forming cloud along the line of sight toward W49A. We use these measurements in conjunction with total mass estimates from {sup 13}CO to infer the shape of the density probability distribution function. This method is complementary to measurements of turbulence via the column density distribution and should be applicable to any molecular cloud with detected CO. We show that turbulence in this cloud is probably compressively driven, with a compressive-to-total Mach number ratio b=M{sub C}/M>0.4. We measure the standard deviation of the density distribution, constraining it to the range 1.5 < σ {sub s} < 1.9, assuming that the density is lognormally distributed. This measurement represents an essential input into star formation laws. The method of averaging over different excitation conditions to produce a model of emission from a turbulent cloud is generally applicable to optically thin line observations.

79 ASTRONOMY AND ASTROPHYSICS↗

An Observationally Trained Markov Model for MJO Propagation

A Markovian stochastic model is developed for studying the propagation of the Madden-Julian Oscillation (MJO). This model represents the daily changes in real time multivariate MJO (RMM) indices as random functions of their current state and background conditions. The probability distribution function of the RMM changes is obtained using a machine learning algorithm trained to maximize MJO forecast skills using observed daily indices of RMM and different modes of variability. Skillful forecasts are obtained for lead times between 8 and 27 days. Large ensemble simulations by the stochastic model show that with monsoonal changes in the background state, MJO propagation across the Maritime Continent (MC) is most likely to be disrupted in boreal spring and summer when MJO events propagate from favorable conditions over the Indian Ocean to unfavorable ones over the MC, and predictability is higher during spring and summer when MJO activity is away from the MC region.

54 ENVIRONMENTAL SCIENCES↗

A count probability cookbok: Spurious effects and the scaling model

We study the errors brought by finite volume effects and dilution effects on the practical determination of the count probability distribution function P(sub N)(n,l), which is the probability of having N objects in a cell of volume l cubed for a set of average number density n. Dilution effects are particularly revelant to the so-called sparse sampling strategy. This work is mainly done in the framework of the Bailan & Schaeffer scaling model, which assumes that the Q-body correlation functions obey the scaling relation Xi(sub Q)(lambda r(sub l),....lambda r(sub Q) = lambda(exp -(Q-1)gamma) Xi(sub Q)(r(sub 1),....r(sub Q)). We use three synthetic samples as references to perform our analysis: a fractal generated by a Rayleigh-Levy random walk with approximately 3 x 10(exp 4) objects, a sample dominated by a spherical power-law cluster with approximately 3 x 10(exp 4) objects and a cold dark matter (CDM) universe involving approximately 3 x 10(exp 5) matter particles.

Colombi, S.↗

Potential Use of a Bayesian Network for Discriminating Flash Type from Future GOES-R Geostationary Lightning Mapper (GLM) data

Continuous monitoring of the ratio of cloud flashes to ground flashes may provide a better understanding of thunderstorm dynamics, intensification, and evolution, and it may be useful in severe weather warning. The National Lighting Detection Network TM (NLDN) senses ground flashes with exceptional detection efficiency and accuracy over most of the continental United States. A proposed Geostationary Lightning Mapper (GLM) aboard the Geostationary Operational Environmental Satellite (GOES-R) will look at the western hemisphere, and among the lightning data products to be made available will be the fundamental optical flash parameters for both cloud and ground flashes: radiance, area, duration, number of optical groups, and number of optical events. Previous studies have demonstrated that the optical flash parameter statistics of ground and cloud lightning, which are observable from space, are significantly different. This study investigates a Bayesian network methodology for discriminating lightning flash type (ground or cloud) using the lightning optical data and ancillary GOES-R data. A Directed Acyclic Graph (DAG) is set up with lightning as a "root" and data observed by GLM as the "leaves." This allows for a direct calculation of the joint probability distribution function for the lighting type and radiance, area, etc. Initially, the conditional probabilities that will be required can be estimated from the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) together with NLDN data. Directly manipulating the joint distribution will yield the conditional probability that a lightning flash is a ground flash given the evidence, which consists of the observed lightning optical data [and possibly cloud data retrieved from the GOES-R Advanced Baseline Imager (ABI) in a more mature Bayesian network configuration]. Later, actual GLM and NLDN data can be used to refine the estimates of the conditional probabilities used in the model; i.e., the Bayesian network is a learning network. Methods for efficient calculation of the conditional probabilities (e.g., an algorithm using junction trees), finding data conflicts, goodness of fit, and dealing with missing data will also be addressed.

Solakiewiz, Richard↗

Vector wind and vector wind shear models 0 to 27 km altitude for Cape Kennedy, Florida, and Vandenberg AFB, California

The techniques are presented to derive several statistical wind models. The techniques are from the properties of the multivariate normal probability function. Assuming that the winds can be considered as bivariate normally distributed, then (1) the wind components and conditional wind components are univariate normally distributed, (2) the wind speed is Rayleigh distributed, (3) the conditional distribution of wind speed given a wind direction is Rayleigh distributed, and (4) the frequency of wind direction can be derived. All of these distributions are derived from the 5-sample parameter of wind for the bivariate normal distribution. By further assuming that the winds at two altitudes are quadravariate normally distributed, then the vector wind shear is bivariate normally distributed and the modulus of the vector wind shear is Rayleigh distributed. The conditional probability of wind component shears given a wind component is normally distributed. Examples of these and other properties of the multivariate normal probability distribution function as applied to Cape Kennedy, Florida, and Vandenberg AFB, California, wind data samples are given. A technique to develop a synthetic vector wind profile model of interest to aerospace vehicle applications is presented.

Smith, O. E.↗

SURVEY DESIGN FOR SPECTRAL ENERGY DISTRIBUTION FITTING: A FISHER MATRIX APPROACH

The spectral energy distribution (SED) of a galaxy contains information on the galaxy's physical properties, and multi-wavelength observations are needed in order to measure these properties via SED fitting. In planning these surveys, optimization of the resources is essential. The Fisher Matrix (FM) formalism can be used to quickly determine the best possible experimental setup to achieve the desired constraints on the SED-fitting parameters. However, because it relies on the assumption of a Gaussian likelihood function, it is in general less accurate than other slower techniques that reconstruct the probability distribution function (PDF) from the direct comparison between models and data. We compare the uncertainties on SED-fitting parameters predicted by the FM to the ones obtained using the more thorough PDF-fitting techniques. We use both simulated spectra and real data, and consider a large variety of target galaxies differing in redshift, mass, age, star formation history, dust content, and wavelength coverage. We find that the uncertainties reported by the two methods agree within a factor of two in the vast majority ({approx}90%) of cases. If the age determination is uncertain, the top-hat prior in age used in PDF fitting to prevent each galaxy from being older than the universe needs to be incorporated in the FM, at least approximately, before the two methods can be properly compared. We conclude that the FM is a useful tool for astronomical survey design.

79 ASTRONOMY AND ASTROPHYSICS↗

How pervasive is the Hirshfeld partitioning?

One can partition the molecular density into its atomic contributions by minimizing the divergence of the atom-in-molecule densities from their corresponding reference pro-atomic densities, subject to the constraint that the sum of the atom-in-molecule densities is the total molecular density. We expose conditions on the divergence measure that are necessary, and sufficient, to recover the popular Hirshfeld partitioning. Specifically, among all local measures of the divergence between two probability distribution functions, the Hirshfeld partitioning is obtained only for f-divergences.

74 ATOMIC AND MOLECULAR PHYSICS↗

INVESTIGATION OF THE ERRORS IN SLOAN DIGITAL SKY SURVEY PROPER-MOTION MEASUREMENTS USING SAMPLES OF QUASARS

We investigate in detail the probability distribution function (pdf) of the proper-motion measurement errors in the SDSS+USNO-B proper-motion catalog of Munn et al. using clean quasar samples. The pdf of the errors is well represented by a Gaussian core with extended wings, plus a very small fraction (<0.1%) of 'outliers'. We find that while formally the pdf could be well fit by a five-parameter fitting function, for many purposes it is also adequate to represent the pdf with a one-parameter approximation to this function. We apply this pdf to the calculation of the confidence intervals on the true proper motion for an SDSS+USNO-B proper-motion measurement, and discuss several scientific applications of the SDSS proper-motion catalog. Our results have various applications in studies of the galactic structure and stellar kinematics. Specifically, they are crucial for searching hyper-velocity stars in the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

THE DEPENDENCE OF PRESTELLAR CORE MASS DISTRIBUTIONS ON THE STRUCTURE OF THE PARENTAL CLOUD

The mass distribution of prestellar cores is obtained for clouds with arbitrary internal mass distributions using a selection criterion based on the thermal and turbulent Jeans mass and applied hierarchically from small to large scales. We have checked this methodology by comparing our results for a log-normal density probability distribution function with the theoretical core mass function (CMF) derived by Hennebelle and Chabrier, namely a power law at large scales and a log-normal cutoff at low scales, but our method can be applied to any mass distributions representing a star-forming cloud. This methodology enables us to connect the parental cloud structure with the mass distribution of the cores and their spatial distribution, providing an efficient tool for investigating the physical properties of the molecular clouds that give rise to the prestellar core distributions observed. Simulated fractional Brownian motion (fBm) clouds with the Hurst exponent close to the value H = 1/3 give the best agreement with the theoretical CMF derived by Hennebelle and Chabrier and Chabrier's system initial mass function. Likewise, the spatial distribution of the cores derived from our methodology shows a surface density of companions compatible with those observed in Trapezium and Ophiucus star-forming regions. This method also allows us to analyze the properties of the mass distribution of cores for different realizations. We found that the variations in the number of cores formed in different realizations of fBm clouds (with the same Hurst exponent) are much larger than the expected root N statistical fluctuations, increasing with H.

79 ASTRONOMY AND ASTROPHYSICS↗

FROM DUST TO PLANETESIMALS: AN IMPROVED MODEL FOR COLLISIONAL GROWTH IN PROTOPLANETARY DISKS

Planet formation occurs within the gas- and dust-rich environments of protoplanetary disks. Observations of these objects show that the growth of primordial submicron-sized particles into larger aggregates occurs at the earliest evolutionary stages of the disks. However, theoretical models of particle growth that use the Smoluchowski equation to describe collisional coagulation and fragmentation have so far failed to produce large particles while maintaining a significant population of small grains. This has generally been attributed to the existence of two barriers impeding growth due to bouncing and fragmentation of colliding particles. In this paper, we demonstrate that the importance of these barriers has been artificially inflated through the use of simplified models that do not take into account the stochastic nature of the particle motions within the gas disk. We present a new approach in which the relative velocities between two particles are described by a probability distribution function that models both deterministic motion (from the vertical settling, radial drift, and azimuthal drift) and stochastic motion (from Brownian motion and turbulence). Taking both into account can give quite different results to what has been considered recently in other studies. We demonstrate the vital effect of two 'ingredients' for particle growth: the proper implementation of a velocity distribution function that overcomes the bouncing barrier and, in combination with mass transfer in high-mass-ratio collisions, boosts the growth of larger particles beyond the fragmentation barrier. A robust result of our simulations is the emergence of two particle populations (small and large), potentially explaining simultaneously a number of longstanding problems in protoplanetary disks, including planetesimal formation close to the central star, the presence of millimeter- to centimeter-sized particles far out in the disk, and the persistence of {mu}m-sized grains for millions of years.

79 ASTRONOMY AND ASTROPHYSICS↗

STAR FORMATION IN TURBULENT MOLECULAR CLOUDS WITH COLLIDING FLOW

Using self-gravitational hydrodynamical numerical simulations, we investigated the evolution of high-density turbulent molecular clouds swept by a colliding flow. The interaction of shock waves due to turbulence produces networks of thin filamentary clouds with a sub-parsec width. The colliding flow accumulates the filamentary clouds into a sheet cloud and promotes active star formation for initially high-density clouds. Clouds with a colliding flow exhibit a finer filamentary network than clouds without a colliding flow. The probability distribution functions (PDFs) for the density and column density can be fitted by lognormal functions for clouds without colliding flow. When the initial turbulence is weak, the column density PDF has a power-law wing at high column densities. The colliding flow considerably deforms the PDF, such that the PDF exhibits a double peak. The stellar mass distributions reproduced here are consistent with the classical initial mass function with a power-law index of –1.35 when the initial clouds have a high density. The distribution of stellar velocities agrees with the gas velocity distribution, which can be fitted by Gaussian functions for clouds without colliding flow. For clouds with colliding flow, the velocity dispersion of gas tends to be larger than the stellar velocity dispersion. The signatures of colliding flows and turbulence appear in channel maps reconstructed from the simulation data. Clouds without colliding flow exhibit a cloud-scale velocity shear due to the turbulence. In contrast, clouds with colliding flow show a prominent anti-correlated distribution of thin filaments between the different velocity channels, suggesting collisions between the filamentary clouds.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-variate joint PDF for non-Gaussianities: exact formulation and generic approximations

We provide an exact expression for the multi-variate joint probability distribution function of non-Gaussian fields primordially arising from local transformations of a Gaussian field. This kind of non-Gaussianity is generated in many models of inflation. We apply our expression to the non-Gaussianity estimation from Cosmic Microwave Background maps and the halo mass function where we obtain analytical expressions. We also provide analytic approximations and their range of validity. For the Cosmic Microwave Background we give a fast way to compute the PDF which is valid up to more than 7σ for f{sub NL} values (both true and sampled) not ruled out by current observations, which consists of expressing the PDF as a combination of bispectrum and trispectrum of the temperature maps. The resulting expression is valid for any kind of non-Gaussianity and is not limited to the local type. The above results may serve as the basis for a fully Bayesian analysis of the non-Gaussianity parameter.

79 ASTRONOMY AND ASTROPHYSICS↗

The Kullback-Leibler divergence as an estimator of the statistical properties of CMB maps

The identification of unsubtracted foreground residuals in the cosmic microwave background maps on large scales is of crucial importance for the analysis of polarization signals. These residuals add a non-Gaussian contribution to the data. We propose the Kullback-Leibler (KL) divergence as an effective, non-parametric test on the one-point probability distribution function of the data. With motivation in information theory, the KL divergence takes into account the entire range of the distribution and is highly non-local. We demonstrate its use by analyzing the large scales of the Planck 2013 SMICA temperature fluctuation map and find it consistent with the expected distribution at a level of 6%. Comparing the results to those obtained using the more popular Kolmogorov-Smirnov test, we find the two methods to be in general agreement.

79 ASTRONOMY AND ASTROPHYSICS↗

Late time solution for interacting scalar in accelerating spaces

We consider stochastic inflation in an interacting scalar field in spatially homogeneous accelerating space-times with a constant principal slow roll parameter ε. We show that, if the scalar potential is scale invariant (which is the case when scalar contains quartic self-interaction and couples non-minimally to gravity), the late-time solution on accelerating FLRW spaces can be described by a probability distribution function (PDF) ρ which is a function of φ/H only, where φ=φ( x-vector ) is the scalar field and H=H(t) denotes the Hubble parameter. We give explicit late-time solutions for ρarrow ρ{sub ∞}(φ/H), and thereby find the order ε corrections to the Starobinsky-Yokoyama result. This PDF can then be used to calculate e.g. various n-point functions of the (self-interacting) scalar field, which are valid at late times in arbitrary accelerating space-times with ε= constant.

79 ASTRONOMY AND ASTROPHYSICS↗

The PAndAS view of the andromeda satellite system. II. Detailed properties of 23 M31 dwarf spheroidal galaxies

We present a comprehensive analysis of the structural properties and luminosities of the 23 dwarf spheroidal galaxies that fall within the footprint of the Pan-Andromeda Archaeological Survey (PAndAS). These dwarf galaxies represent the large majority of Andromeda’s known satellite dwarf galaxies and cover a wide range in luminosity (−11.6≲M{sub V}≲−5.8 or 10{sup 4.2}≲L≲10{sup 6.5} L{sub ⊙}) and surface brightness (25.1≲μ{sub 0}≲29.3 mag arcsec{sup −2}). We confirm most previous measurements, but we find And XIX to be significantly larger than before (r{sub h}=3065{sub −935}{sup +1065} pc, M{sub V}=−10.1{sub −0.4}{sup +0.8}) and cannot derive parameters for And XXVII as it is likely not a bound stellar system. We also significantly revise downward the luminosities of And XV and And XVI, which are now M{sub V}∼−7.5 or L∼10{sup 5} L{sub ⊙}. Finally, we provide the first detailed analysis of Cas II/And XXX, a fairly faint system (M{sub V}=−8.0{sub −0.3}{sup +0.4}) of typical size (r{sub h}=270±50 pc), located in close proximity to the two bright elliptical dwarf galaxies NGC 147 and NGC 185. Combined with the set of homogeneous distances published in an earlier contribution, our analysis dutifully tracks all relevant sources of uncertainty in the determination of the properties of the dwarf galaxies from the PAndAS photometric catalog. We further publish the posterior probability distribution functions of all the parameters we fit for in the form of MCMC chains available online; these inputs should be used in any analysis that aims to remain truthful to the data and properly account for covariance between parameters.

79 ASTRONOMY AND ASTROPHYSICS↗