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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Model for scattering with proliferating resonances: Many coupled square wells

We present a multichannel model for elastic interactions, composed of an arbitrary number of coupled finite square-well potentials, and derive semianalytic solutions for its scattering behavior. Despite the model's simplicity, it is flexible enough to include many coupled short-ranged resonances in the vicinity of the collision threshold, as is necessary to describe ongoing experiments in ultracold molecules and lanthanide atoms. Here, we also introduce a simple but physically realistic statistical ensemble for parameters in this model. We compute the resulting probability distributions of nearest-neighbor resonance spacings and analyze them by fitting to the Brody distribution. We quantify the ability of alternative distribution functions, for resonance spacing and resonance number variance, to describe the crossover regime. The analysis demonstrates that the multichannel square-well model with the chosen ensemble of parameters naturally captures the crossover from integrable to chaotic scattering as a function of closed-channel coupling strength.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Southern Photometric Quasar Catalog from the Dark Energy Survey Data Release 2

We present a catalog of 1.4 million photometrically selected quasar candidates in the southern hemisphere over the ~5000 deg 2 Dark Energy Survey (DES) wide survey area. We combine optical photometry from the DES second data release (DR2) with available near-infrared (NIR) and the all-sky unWISE mid-infrared photometry in the selection. We build models of quasars, galaxies, and stars with multivariate skew-t distributions in the multidimensional space of relative fluxes as functions of redshift (or color for stars) and magnitude. Our selection algorithm assigns probabilities for quasars, galaxies, and stars and simultaneously calculates photometric redshifts (photo-z) for quasar and galaxy candidates. Benchmarking on spectroscopically confirmed objects, we successfully classify (with photometry) 94.7% of quasars, 99.3% of galaxies, and 96.3% of stars when all IR bands (NIR YJHK and WISE W1W2) are available. The classification and photo-z regression success rates decrease when fewer bands are available. Our quasar (galaxy) photo-z quality, defined as the fraction of objects with the difference between the photo-z z p and the spectroscopic redshift z s , |Δz| ≡ |z s - z p |/(1 + z s ) ≤ 0.1, is 92.2% (98.1%) when all IR bands are available, decreasing to 72.2% (90.0%) using optical DES data only. Our photometric quasar catalog achieves an estimated completeness of 89% and purity of 79% at r < 21.5 (0.68 million quasar candidates), with reduced completeness and purity at 21.5 < r ≲ 24. Among the 1.4 million quasar candidates, 87,857 have existing spectra, and 84,978 (96.7%) of them are spectroscopically confirmed quasars. Finally, we provide quasar, galaxy, and star probabilities for all (0.69 billion) photometric sources in the DES DR2 coadded photometric catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Experimental observation of electron bounce resonance through electron energy distribution measurement in a finite size inductively coupled plasma

The electron bounce resonance was experimentally investigated in a low pressure planar inductively coupled plasma. The electron energy probability functions (EEPFs) were measured at different chamber heights and the energy diffusion coefficients were calculated by the kinetic model. It is found that the EEPFs begin to flatten at the first electron bounce resonance condition, and the plateau shifts to a higher electron energy as the chamber height increases. The plateau which indicates strong electron heating corresponds not only to the electron bounce resonance condition but also to the peaks of the first component of the energy diffusion coefficients. As a result, the plateau formation in the EEPFs is mainly due to the electron bounce resonance in a finite inductive discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Minimum Detectable Quantity Calculation for Radiation Portal Monitors

The minimum detectable quantity of a radiation portal monitor is the smallest amount of radioactive material that can be detected passing through the monitor (with a specified detection probability). The minimum detectable quantity is a function of many factors, including isotope, velocity, distance between pillars, height of portal and vehicle, source distribution and position, background radiation, background suppression, detector volume and positions, detector efficiency, decision metrics and algorithms, the presence of NORM (naturally occurring radioactive material), and alarm thresholds. Experimentally testing a radiation portal monitor to failure for all combinations of these factors is extremely time-consuming and often cost prohibitive. This document outlines a straightforward method for quickly estimating the minimum detectable quantity for moving sources over a large variety of conditions with a minimal number of static measurements. The document also describes a proof of concept software application that the International Atomic Energy Agency has designed based on the described method that can be used by Member States to estimate a monitor’s minimum detectable quantity.

Blessinger, Christopher S.↗
DOE OSTI2017

A practical method to test the validity of the standard Gumbel distribution in logit-based multinomial choice models of travel behavior

Most multinomial choice models (e.g., the multinomial logit model) adopted in practice assume an extreme-value Gumbel distribution for the random components (error terms) of utility functions. This distributional assumption offers a closed-form likelihood expression when the utility maximization principle is applied to model choice behaviors. As a result, model coefficients can be easily estimated using the standard maximum likelihood estimation method. However, maximum likelihood estimators are consistent and efficient only if distributional assumptions on the random error terms are valid. It is therefore critical to test the validity of underlying distributional assumptions on the error terms that form the basis of parameter estimation and policy evaluation. In this paper, a practical yet statistically rigorous method is proposed to test the validity of the distributional assumption on the random components of utility functions in both the multinomial logit (MNL) model and multiple discrete-continuous extreme value (MDCEV) model. Based on a semi-nonparametric approach, a closed-form likelihood function that nests the MNL or MDCEV model being tested is derived. The proposed method allows traditional likelihood ratio tests to be used to test violations of the standard Gumbel distribution assumption. Simulation experiments are conducted to demonstrate that the proposed test yields acceptable Type-I and Type-II error probabilities at commonly available sample sizes. The test is then applied to three real-world discrete and discrete-continuous choice models. For all three models, the proposed test rejects the validity of the standard Gumbel distribution in most utility functions, calling for the development of robust choice models that overcome adverse effects of violations of distributional assumptions on the error terms in random utility functions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2016

MCNP6 Modification for the Calculation of Rossi and Feynman Distributions

The Rossi distribution is given by the number of captured neutron pairs as a function of the time interval between the two capture events. This distribution contains a random part, which is given by neutron captures caused by neutrons of different fission chains, and a correlated part, which is given by neutron captures caused by neutrons from the same fission chain. The correlated part exponentially vanishes with the {Delta}t time interval following the decay of prompt neutrons. The Rossi distribution peaks when the time interval between two capture events is small since this condition is satisfied for neutrons generated from the same fission chain. In the Poisson probability distribution, the variance is equal to the mean and the next event probability is independent of the elapsed time since the previous event. Neutron captures events do not follow the Poisson probability distribution because they are correlated events via fission chains. Consequently, the Feynman distribution Y is defined as the variance to mean ratio minus one, of the neutron counts array as a function of a fixed time interval. The MCNP6 source code has been modified to allow the calculation of the Rossi and Feynman distributions. The source code modification takes into account the time when a neutron is emitted by the external neutron source and gives a precision in the timestamp lower than a fraction of nanosecond. The modified MCNP6 version has been used to model KUKA facility driven by a californium neutron source. The prompt neutron decay constants obtained by the Rossi and Feynman methods are in agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Stochastic quasi-steady state approximations for asymptotic solutions of the chemical master equation

In this paper, we propose two methods to carry out the quasi-steady state approximation in stochastic models of enzyme catalytic regulation, based on WKB asymptotics of the chemical master equation or of the corresponding partial differential equation for the generating function. The first of the methods we propose involves the development of multiscale generalisation of a WKB approximation of the solution of the master equation, where the separation of time scales is made explicit which allows us to apply the quasi-steady state approximation in a straightforward manner. To the lowest order, the multi-scale WKB method provides a quasi-steady state, Gaussian approximation of the probability distribution. The second method is based on the Hamilton-Jacobi representation of the stochastic process where, as predicted by large deviation theory, the solution of the partial differential equation for the corresponding characteristic function is given in terms of an effective action functional. The optimal transition paths between two states are then given by those paths that maximise the effective action. Such paths are the solutions of the Hamilton equations for the Hamiltonian associated to the effective action functional. The quasi-steady state approximation is applied to the Hamilton equations thus providing an approximation to the optimal transition paths and the transition time between two states. Using this approximation we predict that, unlike the mean-field quasi-steady approximation result, the rate of enzyme catalysis depends explicitly on the initial number of enzyme molecules. The accuracy and validity of our approximated results as well as that of our predictions regarding the behaviour of the stochastic enzyme catalytic models are verified by direct simulation of the stochastic model using Gillespie stochastic simulation algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Machine learning approach to pattern recognition in nuclear dynamics from the ab initio symmetry-adapted no-core shell model

A novel machine learning approach is used to provide further insight into atomic nuclei and to detect orderly patterns amid a vast data of large-scale calculations. The method utilizes a neural network that is trained on ab initio results from the symmetry-adapted no-core shell model (SA-NCSM) for light nuclei. We show that the SA-NCSM, which expands ab initio applications up to medium-mass nuclei by using dominant symmetries of nuclear dynamics, can reach heavier nuclei when coupled with the machine learning approach. In particular, we find that a neural network trained on probability amplitudes for s- and p-shell nuclear wave functions not only predicts dominant configurations for heavier nuclei but in addition, when tested for the 20 Ne ground state, accurately reproduces the probability distribution. The non-negligible configurations predicted by the network provide an important input to the SA-NCSM for reducing ultralarge model spaces to manageable sizes that can be, in turn, utilized in SA-NCSM calculations to obtain accurate observables. The neural network is capable of describing nuclear deformation and is used to track the shape evolution along the 20-42 Mg isotopic chain, suggesting a shape coexistence that is more pronounced toward the very neutron-rich isotopes. We provide first descriptions of the structure and deformation of 24 Si and 40 Mg of interest to x-ray burst nucleosynthesis, and even of the extremely heavy nuclei such as 166,168 Er and 236 U, that build on first-principles considerations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Lifetimes in Te 124 : Examining critical-point symmetry in the Te nuclei

The Doppler-shift attenuation method following inelastic neutron scattering was used to determine the lifetimes of nuclear levels to 3.3-MeV excitation in 124 Te. Level energies and spins, γ -ray energies and branching ratios, and multipole-mixing ratios were deduced from measured γ-ray angular distributions at incident neutron energies of 2.40 and 3.30 MeV, γ-ray excitation functions, and γγ coincidence measurements. The newly obtained reduced transition probabilities and level energies for 124 Te were compared to critical-point symmetry model predictions. The E(5) and β 4 potential critical-point symmetries were also investigated in 122 Te and 126 Te.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

An actuarial approach to retrofit savings in buildings

An actuarial method has been developed for determining energy savings from retrofits from energy use data for a number of buildings. This method should be contrasted with the traditional method of using pre- and post-retrofit data on the same building. This method supports the U.S. Department of Energy Building Performance Database of real building performance data and related tools that enable engineering and financial practitioners to evaluate retrofits. The actuarial approach derives, from the database, probability density functions (PDFs) for energy savings from retrofits by creating peer groups for the user’s pre post buildings. From the energy use distribution of the two groups, the savings PDF is derived. This provides the basis for engineering analysis as well as financial risk analysis leading to investment decisions. Several technical issues are addressed: The savings PDF is obtained from the pre- and post-PDF through a convolution. Smoothing using kernel density estimation is applied to make the PDF more realistic. The low data density problem can be mitigated through a neighborhood methodology. Correlations between pre and post buildings are addressed to improve the savings PDF. Sample size effects are addressed through the Kolmogorov--Smirnov tests and quantile-quantile plots.

Subbarao, Krishnappa↗
DOE OSTI2015

X{sub max}{sup μ} vs. N{sup μ} from extensive air showers as estimator for the mass of primary UHECR's. Application for the Pierre Auger Observatory

We study the possibility of primary mass estimation for Ultra High Energy Cosmic Rays (UHECR's) using the X{sub max}{sup μ} (the height where the number of muons produced on the core of Extensive Air Showers (EAS) is maximum) and the number N{sup μ} of muons detected on ground. We use the 2D distribution - X{sub max}{sup μ} against N{sup μ} in order to find its sensitivity to the mass of the primary particle. For that, we construct a 2D Probability Function Prob(p,Fe | X{sub max}{sup μ}, N{sup μ}) which estimates the probability that a certain point from the plane (X{sub max}{sup μ}, N{sup μ}) corresponds to a shower induced by a proton, respectively an iron nucleus. To test the procedure, we analyze a set of simulated EAS induced by protons and iron nuclei at energies of 10{sup 19}eV and 20° zenith angle with CORSIKA. Using the Bayesian approach and taking into account the geometry of the infill detectors from the Pierre Auger Observatory, we observe an improvement in the accuracy of the primary mass reconstruction in comparison with the results obtained using only the X{sub max}{sup μ} distributions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Visualisation and outlier detection for probability density function ensembles

Abstract Exploratory data analysis (EDA) for functional data—data objects where observations are entire functions—is a difficult problem that has seen significant attention in recent literature. This surge in interest is motivated by the ubiquitous nature of functional data, which are prevalent in applications across fields such as meteorology, biology, medicine and engineering. Empirical probability density functions (PDFs) can be viewed as constrained functional data objects that must integrate to one and be nonnegative. They show up in contexts such as yearly income distributions, zooplankton size structure in oceanography and in connectivity patterns in the brain, among others. While PDF data are certainly common in modern research, little attention has been given to EDA specifically for PDFs. In this paper, we extend several methods for EDA on functional data for PDFs and compare them on simulated data that exhibit different types of variation, designed to mimic that seen in real‐world applications. We then use our new methods to perform EDA on the breakthrough curves observed in gas transport simulations for underground fracture networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Collisionless “thermalization” in the sheath of an argon discharge

We performed kinetic Vlasov simulations of the plasma-wall transition for a low-pressure argon discharge without external magnetic fields, using the same plasma parameters as in the experiments of Claire et al. [Phys. Plasmas 13, 062103 (2006)]. Experimentally, it was found that the ion velocity distribution function is highly asymmetric in the presheath, but, surprisingly, becomes again close to Maxwellian inside the sheath. Here, we show that this “thermalization” can be explained by purely collisionless effects that are akin to the velocity bunching phenomenon observed in charged particles beams. Such collisionless thermalization is also observed in the presheath region close to the sheath entrance, although it is much weaker there and in practice probably swamped by collisional processes (standard or enhanced by instabilities)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A NEW METHOD FOR DERIVING THE STELLAR BIRTH FUNCTION OF RESOLVED STELLAR POPULATIONS

We present a new method for deriving the stellar birth function (SBF) of resolved stellar populations. The SBF (stars born per unit mass, time, and metallicity) is the combination of the initial mass function (IMF), the star formation history (SFH), and the metallicity distribution function (MDF). The framework of our analysis is that of Poisson Point Processes (PPPs), a class of statistical models suitable when dealing with points (stars) in a multidimensional space (the measurement space of multiple photometric bands). The theory of PPPs easily accommodates the modeling of measurement errors as well as that of incompleteness. Our method avoids binning stars in the color–magnitude diagram and uses the whole likelihood function for each data point; combining the individual likelihoods allows the computation of the posterior probability for the population's SBF. Within the proposed framework it is possible to include nuisance parameters, such as distance and extinction, by specifying their prior distributions and marginalizing over them. The aim of this paper is to assess the validity of this new approach under a range of assumptions, using only simulated data. Forthcoming work will show applications to real data. Although it has a broad scope of possible applications, we have developed this method to study multi-band Hubble Space Telescope observations of the Milky Way Bulge. Therefore we will focus on simulations with characteristics similar to those of the Galactic Bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

An all-sky sample of intermediate-mass star-forming regions

We present an all-sky sample of 984 candidate intermediate-mass Galactic star-forming regions that are color selected from the Infrared Astronomical Satellite (IRAS) Point Source Catalog and morphologically classify each object using mid-infrared Wide-field Infrared Survey Explorer (WISE) images. Of the 984 candidates, 616 are probable star-forming regions (62.6%), 128 are filamentary structures (13.0%), 39 are point-like objects of unknown nature (4.0%), and 201 are galaxies (20.4%). We conduct a study of four of these regions, IRAS 00259+5625, IRAS 00420+5530, IRAS 01080+5717, and IRAS 05380+2020, at Galactic latitudes |b| > 5° using optical spectroscopy from the Wyoming Infrared Observatory, along with near-infrared photometry from the Two-Micron All Sky Survey, to investigate their stellar content. New optical spectra, color-magnitude diagrams, and color-color diagrams reveal their extinctions, spectrophotometric distances, and the presence of small stellar clusters containing 20-78 M {sub ☉} of stars. These low-mass diffuse star clusters contain ∼65-250 stars for a typical initial mass function, including one or more mid-B stars as their most massive constituents. Using infrared spectral energy distributions we identify young stellar objects near each region and assign probable masses and evolutionary stages to the protostars. The total infrared luminosity lies in the range 190-960 L {sub ☉}, consistent with the sum of the luminosities of the individually identified young stellar objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗