Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability distributions”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Disentangling the physics of the attractive Hubbard model as a fully interacting model of fermions via the accessible and symmetry-resolved entanglement entropies

The complicated ways in which electrons interact in many-body systems such as molecules and materials have long been viewed through the lens of local electron correlation and associated correlation functions. However, quantum information science has demonstrated that more global diagnostics of quantum states like the entanglement entropy can provide a complementary and clarifying lens on electronic behavior. One particularly useful measure that can be used to distinguish between quantum and classical sources of entanglement is the accessible entanglement, the entanglement available as a quantum resource for systems subject to conservation laws, such as fixed particle number, due to superselection rules. In this work, we introduce an algorithm and demonstrate how to compute accessible and symmetry-resolved entanglements for interacting fermion systems. This is accomplished by combining an incremental version of the swap algorithm with a recursive auxiliary field quantum Monte Carlo algorithm recently developed by the authors. We apply these tools to study the pairing and charge density waves exhibited in the paradigmatic attractive Hubbard model via entanglement. We find that the particle and spin symmetry-resolved entanglements and their related full probability distribution functions show very clear—and unique—signatures of the underlying electronic behavior even when those features are less pronounced in conventional correlation functions. Altogether, this work provides a systematic means of characterizing the entanglement within quantum systems that can grant a deeper understanding of the complicated electronic behavior that underlies quantum phase transitions and crossovers in many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Bipartite mutual information in classical many-body dynamics

Information theoretic measures have helped to sharpen our understanding of many-body quantum states. As perhaps the most well-known example, the entanglement entropy (or more generally, the bipartite mutual information) has become a powerful tool for characterizing the dynamical growth of quantum correlations. By contrast, although computable, the bipartite mutual information (MI) is almost never explored in classical many particle systems; this owes in part to the fact that computing the MI requires keeping track of the evolution of the full probability distribution, a feat which is rarely done (or thought to be needed) in classical many-body simulations. Here, we utilize the MI to analyze the spreading of information in 1D elementary cellular automata (CA). Broadly speaking, we find that the behavior of the MI in these dynamical systems exhibits a few different types of scaling that roughly correspond to known CA universality classes. Of particular note is that we observe a set of automata for which the MI converges parametrically slowly to its thermodynamic value. We develop a microscopic understanding of this behavior by analyzing a two-species model of annihilating particles moving in opposite directions. Furthermore, our work suggests the possibility that information theoretic tools such as the MI might enable a more fine-grained characterization of classical many-body states and dynamics.

Cellular automata↗
DOE OSTI2020

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Bayesian inference of nuclear symmetry energy from measured and imagined neutron skin thickness in Sn 116 , 118 , 120 , 122 , 124 , 130 , 132 , Pb 208 , and Ca 48

The neutron skin thickness Δr np in heavy nuclei has been known as one of the most sensitive terrestrial probes of the nuclear symmetry energy E sym (ρ) around $\frac{2}{3}$ of the saturation density ρ 0 of nuclear matter. Existing neutron skin data mostly from hadronic observables suffer from large uncertainties and their extraction from experiments are often strongly model dependent. While waiting eagerly for the promised model-independent and high-precision neutron skin data for 208 Pb and 48 Ca from the parity-violating electron scattering experiments (PREX-II and CREX at JLab as well as MREX at MESA), within the Bayesian statistical framework using the Skyrme-Hartree-Fock model we infer the posterior probability distribution functions (PDFs) of the slope parameter L of the nuclear symmetry energy at ρ 0 from imagined Δr np ( 208 Pb)=0.15, 0.20, and 0.30 fm with a 1σ error bar of 0.02, 0.04, and 0.06 fm, respectively, as well as Δr np ( 48 Ca)=0.12, 0.15, and 0.25 fm, with different 1σ error bar of 0.01 and 0.02 fm, respectively. The results are compared with the PDFs of L inferred using the same approach from the available Δr np data for 116, 118, 120, 122, 124, 130, 132 Sn from hadronic probes. They are also compared with results from a recent Bayesian analysis of the radius and tidal deformability data of canonical neutron stars from GW170817 and NICER. The neutron skin data for Sn isotopes gives L = 45.5 $^{+ 26.5}_{-21.6}$ MeV surrounding its mean value or L = 53 . 4 $^{+ 18.6}_{ -29.5}$ MeV surrounding its maximum a posteriori value, respectively, with the latter smaller than but consistent with the L = 66 $^{+ 12}_{-20}$ MeV from the neutron star data within their 68% confidence intervals. We found that Δr np = 0.17 –0.18 fm in 208 Pb with an error bar of about 0.02 fm leads to a PDF of L compatible with that from analyzing the Sn data. To provide additionally useful information on L extracted from the Δr np of Sn isotopes, the experimental error bar of Δr np in 208 Pb should be at least smaller than 0.06 fm aimed by some current experiments. In addition, the Δr np ( 48 Ca) needs to be larger than 0.15 fm but smaller than 0.25 fm to be compatible with the Sn and/or neutron star results. To further improve our current knowledge about L and distinguish its PDFs in the examples considered, even higher precisions leading to significantly less than ±20 MeV error bars for L at 68% confidence level are necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Bayesian uncertainty quantification for nuclear matter incompressibility

Within a Bayesian statistical framework using the standard Skyrme-Hartree-Fock model, the maximum a posteriori (MAP) values and uncertainties of nuclear matter incompressibility and isovector interaction parameters are inferred from the experimental data of giant resonances and neutron-skin thicknesses of typical heavy nuclei. Furthermore, with the uncertainties of the isovector interaction parameters constrained by the data of the isovector giant dipole resonance and the neutron-skin thickness, we have obtained K 0 = 223$^{+7}_{–8}$ MeV at 68% confidence level using the data of the isoscalar giant monopole resonance in 208 Pb measured at the Research Center for Nuclear Physics (RCNP), Japan, and at the Texas A&M University, USA. Although the corresponding 120 Sn data gives a MAP value for K 0 about 5 MeV smaller than the 208 Pb data, there are significant overlaps in their posterior probability distribution functions.

190 ≤ A ≤ 219↗
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↗