Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability distribution”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables

We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory (χ EFT) that are provided by bound-state observables in the A = 3 and A = 4 sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the χ EFT expansion at next-to-next-to-leading order. A consistent solution for the 3 H binding energy, the 4 He binding energy and radius, and the 3 H β-decay rate can only be obtained if χ EFT truncation errors are included in the analysis. Here, the β-decay rate is the only one of these that yields a nondegenerate constraint on the 3NF low-energy constants, which makes it crucial for the parameter estimation. We use eigenvector continuation for fast and accurate emulation of no-core shell model calculations of the few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the parameters that quantify the truncation error. We find a χ EFT expansion parameter of Q = 0.33 ± 0.06 for these observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2022

Angular correlation between the fission fragment intrinsic spins

The generation of fission fragment (FF) spins and of their relative orbital angular momentum has been debated for more than six decades, and no consensus has been achieved so far. The interpretation of recent experimental results of Wilson et al., has been challenged by several recent theoretical studies, which are not in agreement with one another. According to the interpretation by Wilson et al., the FF spins emerge very long after scission occurrs. Randrup and Vogt, while agreeing that the FF spins are uncorrelated, conclude on the basis of a phenomenological model that these spins are uncorrelated already before scission. Bulgac et al., in a fully microscopic study demonstrate that the primordial FF spin final values are defined before the emission of prompt neutrons and statistical gammas and are strongly correlated with a relative angle between spins close to 2π/3, a result in full agreement with the present independent analysis. The prompt neutrons and statistical gammas carry a significant amount of angular momentum according to the study of Stetcu et al. which can lead to a decorrelation of the FF bandhead spins of the yrast lines measured by Wilson et al., and which also provides arguments why these measured spins are so different from the primordial FF spins evaluated by Bulgac et al.. Here, I show that the unexpected character of the angular correlation between the primordial FF intrinsic spins, recently evaluated by Bulgac et al., which favors FF intrinsic spins pointing predominantly in opposite directions, can be understood by using simple general phase space arguments. The observation by Wilson et al., that the FF spins are uncorrelated follows from both the results of the microscopic calculation of Bulgac et al., and the present analysis of the full correlated probability distribution of the FF spins together with the relative orbital angular momentum of the primordial FFs. These arguments may apply also to heavy-ion collisions and, since there is no use of specifics of the particle interactions, the present results might apply to atomic and molecular systems as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Statistical analysis for the neutrinoless double-β-decay matrix element of 48 Ca

Neutrinoless double-β-decay (0⁢vββ) nuclear matrix elements (NME) are the object of many theoretical calculation methods, and are very important for analysis and guidance of a large number of experimental efforts. However, there are large discrepancies between the NME values provided by different methods. Here, in this paper, we propose a statistical analysis of the 48 Ca 0v⁢ββ NME using the interacting shell model, emphasizing the range of the NME probable values and their correlations with observables that can be obtained from the existing nuclear data. Based on this statistical analysis with three independent effective Hamiltonians, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (0.45–0.95) at 90% confidence level, and a mean value of 0.68.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Predicting the neutrinoless double-β-decay matrix element of 136 Xe using a statistical approach

Calculation of the nuclear matrix elements (NMEs) for double-β decay is of paramount importance for guiding experiments and for analyzing and interpreting the experimental data, especially for the search of the neutrinoless double β decay mode (0v⁢ββ). However, there are currently still large differences between the NME values calculated by different methods, hence a quantification of their uncertainties is very much required. Here, in this paper, we propose a statistical analysis of 0⁢v⁢ββ NME for the 136 Xe isotope, based on the interacting shell model, but using three independent effective Hamiltonians, emphasizing the range of the NMEs' most probable values and its correlations with observables that can be obtained from the existing nuclear data. Consequently, we propose a common probability distribution function for the 0⁢vββ NME, which has a range of (1.55–2.65) at 90% confidence level, with a mean value of 1.99 and a standard deviation of 0.37.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Production and discovery of neutron-rich isotopes by fragmentation of 198 Pt

Production cross sections were measured for fragments produced by an 85 MeV/u 198 Pt beam incident on a beryllium target. Event-by-event particle identification of A, Z, and q for the reaction products was performed by employing energy loss, time-of-flight, magnetic rigidity, and total kinetic energy measurements. Over 70 nuclei in the Hf-Pt region were identified, including three isotopes first observed in this work: 191,192 Hf and 189 Lu. Due to the existence of multiple charge states between H-like and C-like ions, a new analysis method was introduced, incorporating Monte Carlo calculations of charge state fractions for a given charge state of the projectile residue just after the reaction. For the first time, charge-state probability distribution functions after the reaction have been deduced from experimental data. Furthermore, this study provides insight into how to produce key nuclides near N = 126 and the ability of a fragmentation residue to retain electrons from the primary beam.

190 ≤ A ≤ 219↗
DOE OSTI2025

Relaxation times for disoriented isospin condensates in high energy heavy ion collisions

Fluctuations between charged and neutral kaons measured by the ALICE Collaboration in Pb-Pb collisions at the CERN Large Hadron Collider (LHC) exceed conventional explanations. Previously it was shown that if the scalar condensate is accompanied by an electrically neutral isospin-1 field then the combination can produce large equilibrium fluctuations where ⟨$\overline{𝑢}$ ⁢𝑢⟩≠⟨ $\overline{𝑑}$𝑑⟩. Hadronizing strange and antistrange quarks might then strongly fluctuate between charged ($\overline{𝑢}$⁢𝑠 or $\overline{𝑠}$⁢𝑢) and neutral (𝑑⁢$\overline{𝑠}$⁢ or 𝑠⁢$\overline{𝑑}$) kaons. Here, we estimate the times for the condensates to achieve their equilibrium probability distributions within causal volumes in high energy heavy ion collisions. This is achieved by modeling the temperature dependence of the condensates, mesonic collective excitations, decay rates of the associated fields, and employing the Langevin and Fokker-Planck equations. Within this model, we find that the equilibration times are short compared with the expansion time, making disoriented isospin condensates a viable explanation for the anomalous fluctuations observed at the LHC.

Quark-gluon plasma↗
DOE OSTI2020

Forecasting cosmological constraints from the weak lensing magnification of type Ia supernovae measured by the Nancy Grace Roman Space Telescope

We report that the weak lensing magnification of type Ia supernovae (SNe Ia) is sensitive to the clustering of matter and provides an independent cosmological probe complementary to SN Ia distance measurements. The Nancy Grace Roman Space Telescope is uniquely sensitive to this measurement as it can discover high redshift SNe Ia and measure them with high precision. We present a methodology for reconstructing the probability distribution of the weak lensing magnification μ of SNe Ia, p(μ), from observational data, and using it to constrain cosmological parameters. We find that the reconstructed p(μ) can be fitted accurately by a stretched Gaussian distribution and used to measure the variance of μ, ξ μ , which can be compared to theoretical predictions in a likelihood analysis. Applying our methodology to a set of realistically simulated SNe Ia expected from the Roman Space Telescope, we find that using the weak lensing magnification of the SNe Ia constrains a combination of matter density Ω m and matter clustering amplitude σ 8 . SN Ia distances alone lead to a better than 1% measurement of Ωm. The combination of SN Ia weak lensing magnification and distance measurements result in a ~10 % measurement on σ 8 . The SNe Ia from Roman will be powerful in constraining the cosmological model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Bayesian inference of fine features of the nuclear equation of state from future neutron star radius measurements to 0.1 km accuracy

To more precisely constrain the equation of state (EOS) of supradense neutron-rich nuclear matter, future high-precision x-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a metamodel EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility K 0 and skewness J 0 of symmetric nuclear matter (SNM) as well as the slope L, curvature K sym , and skewness J sym characterizing the density dependence of nuclear symmetry energy E sym ⁡(ρ), respectively, from mean values of NS radii consistent with existing observations and an expected accuracy Δ⁢R ranging from about 1.0 to 0.1 km. Here, we found that (1) the Δ⁢R has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller Δ⁢R reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density E sym ⁡(ρ), (3) a double-peak feature of the PDF(K sym ) corresponding to the strong K sym – J sym and K sym – L anticorrelations is revealed when Δ⁢R is less than about 0.2 km, and the locations of the two peaks are sensitive to the maximum value of J sym reflecting the stiffness of E sym ⁡(ρ) above about 3 times the saturation density ρ 0 of SNM, and (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2–3)⁢ρ 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Impact of baryonic feedback on HSC-Y1 weak lensing non-Gaussian statistics

Baryonic feedback is a major systematic in weak lensing cosmology. Its most studied effect is the suppression of the lensing power spectrum, a second-order statistic, on small scales. Motivated by the growing interest in statistics beyond the second order, we investigate the effect of baryons on lensing non-Gaussian statistics and the resulting biases in the matter clustering amplitude S 8 = σ 8 Ω m / 0.3 . We focus on the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) data which, with its high source number density, closely resembles those expected from the upcoming Euclid and Rubin Legacy Survey of Space and Time. We study four non-Gaussian statistics of convergence maps—peak counts, minimum counts, the probability distribution function, and the scattering transform coefficients—in addition to the usual power spectrum. We first estimate the biases in S 8 using mock observations built from the IllustrisTNG and BAHAMAS hydrodynamical simulations and theoretical models built from dark-matter-only simulations. We find up to 1 σ bias in S 8 when the smallest scales (2 arcmin) and the highest feedback level are considered. We then analyze the HSC-Y1 data and compare the S 8 obtained for each statistic with different smoothing scales or scale cuts. As we expect that baryons mostly affect the small scales, comparing the results obtained from including and excluding small scales can indicate the level of impact from baryons. With HSC data, we find only minor ( ≤ 0.5 σ ) differences in S 8 for all statistics, even when considering very small scales (2 arcmin). Our results suggest that the effect of baryons is insignificant at the level of HSC-Y1 down to 2 arcmin for all statistics examined here, or it is canceled by other scale-dependent systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Maximally entangled gluons for any x

Individual quarks and gluons at small x inside an unpolarized hadron can be regarded as Bell states in which qubits in the spin and orbital angular momentum spaces are maximally entangled. Using the machinery of quantum information science, we generalize this observation to all values 0 < x <1 and describe gluons (but not quarks) as maximally entangled states between a qubit and a qudit. We introduce the conditional probability distribution P⁡(l z |s z ) of a gluon’s orbital angular momentum l z given its helicity s z . Restricting to the three states l z =0,±1, which constitute a qutrit, we explicitly compute P as a function of x.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Constraints on dark photon dark matter from Lyman- α forest simulations and an ultrahigh signal-to-noise quasar spectrum

The ultralight dark photon is a well-motivated, hypothetical dark matter candidate. In a dilute plasma, they can resonantly convert into photons, and heat up the intergalactic medium between galaxies. In this work, we explore the dark photon dark matter parameter space by comparing synthetic Lyman- α forest data from cosmological hydrodynamical simulations to observational data from VLT/UVES of the quasar HE0940-1050 ( z em = 3.09 ). We use a novel flux normalization technique that targets underdense gas, reshaping the flux probability distribution. Not only do we place robust constraints on the kinetic mixing parameter of dark photon dark matter, but notably our findings suggest that this model can still reconcile simulated and observed Doppler parameter distributions of z ∼ 0 Lyman- α lines, as seen by HST/COS. This work opens new pathways for the use of the Lyman- α forest to explore new physics, and can be extended to other scenarios such as primordial black hole evaporation, dark matter decay, and annihilation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Mean-field control for efficient mixing of energy loads

In this work, we pose an engineering challenge of controlling an ensemble of energy devices via coordinated, implementation-light, and randomized on/off switching as a problem in nonequilibrium statistical mechanics. We show that mean-field control with nonlinear feedback on the cumulative consumption, assumed available to the aggregator via direct physical measurements of the energy flow, allows the ensemble to recover from its use in the demand response regime, i.e., transition to a statistical steady state, significantly faster than in the case of the fixed feedback. Moreover when the nonlinearity is sufficiently strong, one observes the phenomenon of “super-relaxation,” where the total instantaneous energy consumption of the ensemble transitions to the steady state much faster than the underlying probability distribution of the devices over their state space, while also leaving almost no devices outside of the comfort zone.

25 ENERGY STORAGE↗
DOE OSTI2020

Weak scaling of the contact distance between two fluctuating interfaces with system size

A pair of flat parallel surfaces, each freely diffusing along the direction of their separation, will eventually come into contact. If the shapes of these surfaces also fluctuate, then contact will occur when their centers-of-mass remain separated by a nonzero distance ℓ. An example of such a situation is the motion of interfaces between two phases at conditions of thermodynamic coexistence, and in particular the annihilation of domain wall pairs under periodic boundary conditions. Here we present a general approach to calculate the probability distribution of the contact distance ℓ and determine how its most likely value ℓ* depends on the surfaces' lateral size L. Using the Edward-Wilkinson equation as a model for interfaces, we demonstrate that ℓ* scales weakly with system size, i.e., the dependence of ℓ* on L for both (1+1)- and (2+1)-dimensional interfaces is such that lim L→∞ (ℓ*/L) = 0. In particular, for (2+1)-dimensional interfaces ℓ* is an algebraic function of logL, a result that is confirmed by computer simulations of slab-shaped domains formed under periodic boundary conditions. Overall, this weak scaling implies that such domains remain topologically intact until ℓ becomes very small compared to the lateral size of the interface, contradicting expectations from equilibrium thermodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Using kernel-based statistical distance to study the dynamics of charged particle beams in particle-based simulation codes

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. Here, the focus is on kernel-based methods such as maximum mean discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics

Quantum computers can be used to simulate nonlinear non-Hamiltonian classical dynamics on phase space by using the generalized Koopman–von Neumann formulation of classical mechanics. The Koopman–von Neumann formulation implies that the conservation of the probability distribution function on phase space, as expressed by the Liouville equation, can be recast as an equivalent Schrödinger equation on Hilbert space with a Hermitian Hamiltonian operator and a unitary propagator. This Schrödinger equation is linear in the momenta because it derives from a constrained Hamiltonian system with twice the classical phase-space dimension. A quantum computer with finite resources can be used to simulate a finite-dimensional approximation of this unitary evolution operator. Quantum simulation of classical dynamics is exponentially more efficient than a deterministic Eulerian discretization of the Liouville equation if the Koopman–von Neumann Hamiltonian is sparse. Utilizing quantum walk techniques for state preparation and amplitude estimation for the calculation of observables leads to a quadratic improvement over classical probabilistic Monte Carlo algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗