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 793 records · Page 44

Field ionization of Rydberg atoms in a single-cycle pulse

Here, we study the ionization of Rydberg atoms in a single-cycle-pulse electric field based on both classical and quantum calculations. The ionization-probability curve exhibits a “ripple” structure as a function of the pulse duration and the field amplitude. These ripple structures are found to be dependent on the angular distribution of the initial state. A large electron-emission asymmetry is observed, and the ionized electron is almost completely emitted to one side of the atom except when the pulse length is roughly one Rydberg period. In both the long-pulse and the short-pulse regimes, larger electron energy can be expected from the ionization of lower-lying Rydberg states, matching the observation in a recent experiment [S. Li and R. R. Jones, Phys. Rev. Lett. 112, 143006 (2014)]. This trend is closely related to the electron-emission asymmetry associated with the field-direction change in a single-cycle pulse. The possible implications of the different energy transfer in a single-cycle pulse from that in a multicycle pulse are also discussed briefly for the strong-field ionization of the ground atomic state.

74 ATOMIC AND MOLECULAR PHYSICS↗

Basis light-front quantization for a chiral nucleon-pion Lagrangian

We present the first application of the basis light-front quantization method to a simple chiral model of the nucleon-pion system as a relativistic bound state for the physical proton. In this work, the light-front mass-squared matrix of the nucleon-pion system is obtained within a truncated basis. The mass and the corresponding light-front wave function (LFWF) of the proton are computed by numerical diagonalization of the resulting mass-squared matrix. With the boost-invariant LFWF, we calculate the probability density distribution of the pion's longitudinal momentum fraction and the Dirac form factor of the proton.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scrambling and quantum chaos indicators from long-time properties of operator distributions

Scrambling is a key concept in the analysis of nonequilibrium properties of quantum many-body systems. Most studies focus on its characterization via out-of-time-ordered correlation (OTOC) functions, particularly through the early-time decay of the OTOC. However, scrambling is a complex process which involves operator spreading and operator entanglement, and a full characterization requires one to access more refined information on the operator dynamics at several timescales. In this work we analyze operator scrambling by expanding the target operator in a complete basis and studying the structure of the expansion coefficients treated as a coarse-grained probability distribution in the space of operators. Here we study different features of this distribution, such as its mean, variance, and participation ratio, for the Ising model with longitudinal and transverse fields, kicked collective spin models, and random circuit models. We show that the long-time properties of the operator distribution display common features across these cases and discuss how these properties can be used as a proxy for the onset of quantum chaos. Finally, we discuss the connection with OTOCs and analyze the cost of probing the operator distribution experimentally using these correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

How to Obtain the Redshift Distribution from Probabilistic Redshift Estimates

Abstract A reliable estimate of the redshift distribution n ( z ) is crucial for using weak gravitational lensing and large-scale structures of galaxy catalogs to study cosmology. Spectroscopic redshifts for the dim and numerous galaxies of next-generation weak-lensing surveys are expected to be unavailable, making photometric redshift (photo- z ) probability density functions (PDFs) the next best alternative for comprehensively encapsulating the nontrivial systematics affecting photo- z point estimation. The established stacked estimator of n ( z ) avoids reducing photo- z PDFs to point estimates but yields a systematically biased estimate of n ( z ) that worsens with a decreasing signal-to-noise ratio, the very regime where photo- z PDFs are most necessary. We introduce Cosmological Hierarchical Inference with Probabilistic Photometric Redshifts ( CHIPPR ), a statistically rigorous probabilistic graphical model of redshift-dependent photometry that correctly propagates the redshift uncertainty information beyond the best-fit estimator of n ( z ) produced by traditional procedures and is provably the only self-consistent way to recover n ( z ) from photo- z PDFs. We present the chippr prototype code, noting that the mathematically justifiable approach incurs computational cost. The CHIPPR approach is applicable to any one-point statistic of any random variable, provided the prior probability density used to produce the posteriors is explicitly known; if the prior is implicit, as may be the case for popular photo- z techniques, then the resulting posterior PDFs cannot be used for scientific inference. We therefore recommend that the photo- z community focus on developing methodologies that enable the recovery of photo- z likelihoods with support over all redshifts, either directly or via a known prior probability density.

79 ASTRONOMY AND ASTROPHYSICS↗

Microgrid energy scheduling under uncertain extreme weather: Adaptation from parallelized reinforcement learning agents

Microgrids are useful solutions for integrating renewable energy resources and providing seamless green electricity to minimize carbon footprint. In recent years, extreme weather events happened often worldwide and caused significant economic and societal losses. Such events bring uncertainties to the microgrid energy scheduling problems and increase the challenges of microgrid operation. Traditional optimization approaches suffer from the inaccuracy of the uncertain microgrid model and the unseen events. Existing reinforcement learning (RL) - based approaches are also hampered by the limited generalization and the increasing computational burden when stochastic formulations are required to accommodate the uncertainties. This paper proposes a new parallelized reinforcement learning (PRL) method based on the probabilistic events to handle the microgrid energy uncertainties. Specifically, several local learning agents are employed to interact with pertinent microgrid environments in a distributed manner and report outcomes to the global agent, which will optimize microgrid energy resources online during extreme events. The stochastic microgrid energy optimization problem is reformulated to include all possible scenarios with probabilities. The advantage estimate functions of learning agents are designed with a backward sweep to transfer the outcomes to the value function updating process. Two simulation studies, stochastic optimization and online testing, are performed to compare with several existing RL approaches. Results substantiate that the proposed PRL method can achieve up to 20% optimization performance improvement with 4 and 28 times less computation cost than Q-learning with experience replay and multi-agent Q-learning approaches, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Influence of conformational changes in complex molecules on photon statistics of single molecule fluorescence

A single complex molecule with conformational changes (conformations 0 and 2) is considered. When such a molecule is irradiated by cw-laser light it can randomly change the intensity or polarization of its fluorescence due to jumps from one conformation to another. In fact, the molecule manifests itself either like the 0-type or the 2-type emitter. An expression for the matrix s{sub {alpha}{beta}}(t) called the start-stop correlator (waiting time distribution) in which {alpha}=0,2 and {beta}=0,2 is derived. An expression for the matrix p{sub {alpha}{beta}}(t) called the full correlator is derived as well. It determines the density of the probability of finding an event of {alpha} type and an event of {beta} type separated by time interval t. A relation between matrices s{sub {alpha}{beta}}(t) and p{sub {alpha}{beta}}(t) is found. A mathematical expression for the distribution w{sub N}(T) of events measured in time interval T is derived. It is expressed solely via the matrix s{sub {alpha}{beta}}(t). Numerical calculations of the event distribution function for various rates of intra- and interconformational jumps are carried out with the help of the formula for w{sub N}(T) and by the Monte Carlo method. Both methods of the calculation yield identical distributions. Fluctuating fluorescence intensity I(t) for a bin time of 5 ms is calculated for slow and fast interconformational jumps. A relation is found between the autocorrelation function g{sup (2)}(t) of fluorescence measurable in experiments and the matrix p{sub {alpha}{beta}}(t) calculated theoretically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Collective properties of neutron-star X-ray binary populations of galaxies. II. Pre-low-mass X-ray binary properties, formation rates, and constraints

We continue our exploration of the collective properties of neutron-star X-ray binaries in the stellar fields (i.e., outside globular clusters) of normal galaxies. In Paper I of this series, we considered high-mass X-ray binaries (HMXBs). In this paper (Paper II), we consider low-mass X-ray binaries (LMXBs), whose evolutionary scenario is very different from that of HMXBs. We consider the evolution of primordial binaries up to the stage where the neutron star just formed in the supernova explosion of the primary is in a binary with its low-mass, unevolved companion, and this binary has circularized tidally, producing what we call a pre-low-mass X-ray binary (pre-LMXB). We study the constraints on the formation of such pre-LMXBs in detail (since these are low-probability events), and calculate their collective properties and formation rates. To this end, we first consider the changes in the binary parameters in the various steps involved, viz., the common-envelope phase, the supernova, and the tidal evolution. This naturally leads to a clarification of the constraints. We then describe our calculation of the evolution of the distributions of primordial binary parameters into those of pre-LMXB parameters, following the standard evolutionary scenario for individual binaries. We display the latter as both bivariate and monovariate distributions, discuss their essential properties, and indicate the influences of some essential factors on these. Finally, we calculate the formation rate of these pre-LMXBs. The results of this paper will be used in a subsequent one to compute the expected X-ray luminosity function of LMXBs.

79 ASTRONOMY AND ASTROPHYSICS↗

New Measurements of the Deuteron-to-Proton 𝐹 2 Structure-Function Ratio

Nucleon structure functions, as measured in lepton-nucleon scattering, have historically provided a critical observable in the study of partonic dynamics within the nucleon. However, at very large parton momenta, it is both experimentally and theoretically challenging to extract parton distributions due to the probable onset of nonperturbative contributions and the unavailability of high-precision data at critical kinematics. Extraction of the neutron structure and the d quark distribution have been further challenging because of the necessity of applying nuclear corrections when utilizing scattering data from a deuteron target to extract the free neutron structure. However, a program of experiments has been carried out recently at the energy-upgraded Jefferson Lab electron accelerator aimed at significantly reducing the nuclear correction uncertainties on the d quark distribution function at large partonic momentum. This allows leveraging the vast body of deuterium data covering a large kinematic range to be utilized for d quark parton distribution function extraction. In this Letter, we present new data from experiment E12-10-002, carried out in Jefferson Lab Experimental Hall C, on the deuteron to proton cross section ratio at large Bjorken 𝑥. These results significantly improve the precision of existing data and provide a first look at the expected impact on quark distributions extracted from parton distribution function fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗

Evaluation of a Class of Simple and Effective Uncertainty Methods for Sparse Samples of Random Variables and Functions

When very few samples of a random quantity are available from a source distribution of unknown shape, it is usually not possible to accurately infer the exact distribution from which the data samples come. Under-estimation of important quantities such as response variance and failure probabilities can result. For many engineering purposes, including design and risk analysis, we attempt to avoid under-estimation with a strategy to conservatively estimate (bound) these types of quantities -- without being overly conservative -- when only a few samples of a random quantity are available from model predictions or replicate experiments. This report examines a class of related sparse-data uncertainty representation and inference approaches that are relatively simple, inexpensive, and effective. Tradeoffs between the methods' conservatism, reliability, and risk versus number of data samples (cost) are quantified with multi-attribute metrics use d to assess method performance for conservative estimation of two representative quantities: central 95% of response; and 10 -4 probability of exceeding a response threshold in a tail of the distribution. Each method's performance is characterized with 10,000 random trials on a large number of diverse and challenging distributions. The best method and number of samples to use in a given circumstance depends on the uncertainty quantity to be estimated, the PDF character, and the desired reliability of bounding the true value. On the basis of this large data base and study, a strategy is proposed for selecting the method and number of samples for attaining reasonable credibility levels in bounding these types of quantities when sparse samples of random variables or functions are available from experiments or simulations.

97 MATHEMATICS AND COMPUTING↗

Pulse height model for deuterated scintillation detectors

This work includes the theoretical calculation of pulse height distribution and elaborate measurement at Ohio University using deuterated organic scintillation detector, which allows us to thoroughly investigate the light output spectrum induced by fast neutrons and convincingly compensate the pulse deficiency due to application of threshold to minimize noises either from electronics or environmental background. Specifically, a complete and self-contained analytic model of pulse height distribution for EJ-315 is created. In the model, considering impacts of detector geometry on light output and neutron leakage, a non-uniform angular distribution is also taken into account. An energy range of incident neutron from 1 MeV to 10 MeV is investigated. In return, based on the Peierls-formula, the completed calculation exhibits quantitative probability of fast neutron slowing down history due to stochastic collisions with carbon or deuterium nuclei. The quantitative collision probability conjoining with light output function, which was measured and corrected in view of detector response to several gamma rays, for continuous neutron flux generated by the tandem Van De Graaff Generator gives the light output for various collisions. Furthermore, the individual collision structure of the light pulse spectrum is analyzed and compared to the corresponding pulse height distribution from experimental measurement. Lastly, the resulting light pulse height distribution generated in this model offers intriguing advantages for identifying or characterizing nuclear materials and estimating potentials of various fast recoil detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Electronic properties of electron and hole in type-II semiconductor nano-heterostructures

In this project, we record the orbitals of electron and hole in type-II (CdTe/CdSe/CdTe/CdSe) semiconductor nanocrystal using effective mass approximation. In type-II the band edges of both valance and conduction band are higher than that of shell. So the electron and hole get confined in different layers of the hetero-structure. The energy eigen values and eigen functions are calculated by solving Schrodinger equation using finite difference matrix method. Based on this we investigate the effect of shell thickness and well width on energy and probability distribution of ground state (1s) and few excited states (1p,1d,etc). Our results predict that, type-II quantum dots have significant importance in photovoltaic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-E-T-227: Could the Alpha/Beta Ratio Change in a Strong Magnetic Field?

Purpose: Magnetic resonance imaging (MRI) is being integrated into radiotherapy delivery for MRI-guided radiotherapy. The presence of a strong magnetic field from a MRI machine during radiotherapy delivery presents a new challenge since the trajectories of electrons liberated by ionizing radiation in patients are strongly dependent on the applied magnetic field. The purpose of this work is to explore the potential effect of a strong magnetic field on the α/β ratio, an important radiobiological parameter in radiotherapy. Methods: Based on the theory of dual radiation action, the α/β ratio can be expressed by an integral of the product of two microdosimetry quantities γ(x) and t(x), where γ(x) is the probability that two energy transfers, a distance x apart, results in a lesion, and t(x) is the proximity function, which is the energy-weighted point-pair distribution of distances between energy transfer points in a track. The quantity t(x) depends on the applied magnetic field. An analytical approach has been used to derive a formula that can be used to calculate the α/β ratio in an extremely strong magnetic field. Results: The α/β ratio has been evaluated in the special case when the applied magnetic field approaches infinity, which gives the upper limit of the potential change of the α/β ratio due to the presence of a strong magnetic field. For V79 Chinese hamster cells it has been shown that the α/β ratio could be increased by 2.90 times for Pd-103, 2.97 times for I-125 and about 2.3 times for Co-60 sources when the applied magnetic field approaches infinity. Conclusion: It has been shown theoretically that the α/β ratio can change in a strong magnetic field, and there could be up to a nearly three-fold increase in the α/β ratio, depending on the strength of the applied magnetic field, the cell type and the radiation used.

60 APPLIED LIFE SCIENCES↗

Risk-based System Upgrade Planning for AOT-IC [Capstone Project]

The Accelerator Operations and Technology Instrumentation and Controls Group (AOT IC) at the Los Alamos Neutron Science Center (LANSCE) has established a comprehensive risk-based system upgrade planning strategy which has been utilized as a tool to prioritize system upgrade projects for the past several years. A challenge with the current system is that the group often lacks the data to quantify the probability and impact of system failures, so the group often relies instead on qualitative risk assessments as well as evaluations of a potential upgrade project’s alignment with group strategy and vision to prioritize projects. A proposed framework to enable quantitative risk probability and impact assessments has been developed and tested on three pilot systems chosen to broadly represent the types and conditions of equipment used by the group. The proposed framework incorporates availability data gathered from the LANSCE control room logbooks and the LANSCE work control system, as well as a system health evaluation which is conducted from a variety of resources to assess the probability of equipment failure. The impact of system failure is viewed from the perspective of impact on mission and schedule, where spares status, system documentation status, and the functional distribution of deployed systems are used to quantify these impacts.

43 PARTICLE ACCELERATORS↗

A new subgrid-scale representation of hydrometeor fields using a multivariate PDF

Abstract. The subgrid-scale representation of hydrometeor fields is important for calculating microphysical process rates. In order to represent subgrid-scale variability, the Cloud Layers Unified By Binormals (CLUBB) parameterization uses a multivariate probability density function (PDF). In addition to vertical velocity, temperature, and moisture fields, the PDF includes hydrometeor fields. Previously, hydrometeor fields were assumed to follow a multivariate single lognormal distribution. Now, in order to better represent the distribution of hydrometeors, two new multivariate PDFs are formulated and introduced.The new PDFs represent hydrometeors using either a delta-lognormal or a delta-double-lognormal shape. The two new PDF distributions, plus the previous single lognormal shape, are compared to histograms of data taken from large-eddy simulations (LESs) of a precipitating cumulus case, a drizzling stratocumulus case, and a deep convective case. Finally, the warm microphysical process rates produced by the different hydrometeor PDFs are compared to the same process rates produced by the LES.

97 MATHEMATICS AND COMPUTING↗

Branch-and-Bound algorithm applied to uncertainty quantification of a Boiling Water Reactor Station Blackout

Evaluation of the impacts of uncertainty and sensitivity in modeling presents a significant set of challenges in particular to high fidelity modeling. Computational costs and validation of models creates a need for cost effective decision making with regards to experiment design. Experiments designed to validate computation models can be used to reduce uncertainty in the physical model. In some cases, large uncertainty in a particular aspect of the model may or may not have a large impact on the final results. For example, modeling of a relief valve may result in large uncertainty, however, the actual effects on final peak clad temperature in a reactor transient may be small and the large uncertainty with respect to valve modeling may be considered acceptable. Additionally, the ability to determine the adequacy of a model and the validation supporting it should be considered within a risk informed framework. Low fidelity modeling with large uncertainty may be considered adequate if the uncertainty is considered acceptable with respect to risk. In other words, models that are used to evaluate the probability of failure should be evaluated more rigorously with the intent of increasing safety margin. Probabilistic risk assessment (PRA) techniques have traditionally been used to identify accident conditions and transients. Traditional classical event tree methods utilize analysts’ knowledge and experience to identify the important timing of events in coordination with thermal-hydraulic modeling. These methods lack the capability to evaluate complex dynamic systems. In these systems, time and energy scales associated with transient events may vary as a function of transition times and energies to arrive at a different physical state. Dynamic PRA (DPRA) methods provide a more rigorous analysis of complex dynamic systems. Unfortunately DPRA methods introduce issues associated with combinatorial explosion of states. This study presents a methodology to address combinatorial explosion using a Branch-and-Bound algorithm applied to Dynamic Event Trees (DET), which utilize LENDIT (L – Length, E – Energy, N – Number, D – Distribution, I – Information, and T – Time) as well as a set theory to describe system, state, resource, and response (S2R2) sets to create bounding functions for the DET. The optimization of the DET in identifying high probability failure branches is extended to create a Phenomenological Identification and Ranking Table (PIRT) methodology to evaluate modeling parameters important to safety of those failure branches that have a high probability of failure. The PIRT can then be used as a tool to identify and evaluate the need for experimental validation of models that have the potential to reduce risk. Finally, in order to demonstrate this methodology, a Boiling Water Reactor (BWR) Station Blackout (SBO) case study is presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

2D k -th nearest neighbour statistics: a highly informative probe of galaxy clustering

ABSTRACT Beyond standard summary statistics are necessary to summarize the rich information on non-linear scales in the era of precision galaxy clustering measurements. For the first time, we introduce the 2D k-th nearest neighbour (kNN) statistics as a summary statistic for discrete galaxy fields. This is a direct generalization of the standard 1D kNN by disentangling the projected galaxy distribution from the redshift-space distortion signature along the line-of-sight. We further introduce two different flavours of 2D kNNs that trace different aspects of the galaxy field: the standard flavour which tabulates the distances between galaxies and random query points, and a ‘DD’ flavour that tabulates the distances between galaxies and galaxies. We showcase the 2D kNNs’ strong constraining power both through theoretical arguments and by testing on realistic galaxy mocks. Theoretically, we show that 2D kNNs are computationally efficient and directly generate other statistics such as the popular two-point correlation function (2PCF), voids probability function, and counts-in-cell statistics. In a more practical test, we apply the 2D kNN statistics to simulated galaxy mocks that fold in a large range of observational realism and recover parameters of the underlying extended halo occupation distribution (HOD) model that includes velocity bias and galaxy assembly bias. We find unbiased and significantly tighter constraints on all aspects of the HOD model with the 2D kNNs, both compared to the standard 1D kNN, and the classical redshift-space 2PCF.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of dilute aqueous NaCl solution on caffeine aggregation

The effect of salt concentration on association properties of caffeine molecule was investigated by employing molecular dynamics simulations in isothermal-isobaric ensemble of eight caffeine molecules in pure water and three different salt (NaCl) concentrations, at 300 K temperature and 1 atm pressure. The concentration of caffeine was taken almost at the solubility limit. With increasing salt concentration, we observe enhancement of first peak height and appearance of a second peak in the caffeine-caffeine distribution function. Furthermore, our calculated solvent accessible area values and cluster structure analyses suggest formation of higher order caffeine cluster on addition of salt. The calculated hydrogen bond properties reveal that there is a modest decrease in the average number of water-caffeine hydrogen bonds on addition of NaCl salt. Also observed are: (i) decrease in probability of salt contact ion pair as well as decrease in the solvent separated ion pair formation with increasing salt concentration, (ii) a modest second shell collapse in the water structure, and (iii) dehydration of hydrophobic atomic sites of caffeine on addition of NaCl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗