Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Reevaluation of the 22 Ne(α,γ) 26 Mg and 22 Ne(α,n) 25 Mg reaction rates

The competing 22 Ne(α, γ) 26 Mg and 22 Ne(α,n) 25 Mg reactions control the production of neutrons for the weak s process in massive and asymptotic giant branch (AGB) stars. In both systems, the ratio between the corresponding reaction rates strongly impacts the total neutron budget and strongly influences the final nucleosynthesis. A number of experimental studies have been performed over recent years which necessitate the reevaluation of the 26 Mg and 22 Ne(α, γ) 25 Mg reactions control the production 26 Mg and 22 Ne(α, n) 25 Mg reaction rates. Evaluations of the reaction rates following the collection of new nuclear data presently show differences of up to a factor of 500, resulting in considerable uncertainty in the resulting nucleosynthesis. To reevaluate the 22 Ne(α, γ) 26 Mg and 22 Ne(α,n) 25 Mg reaction rates using updated nuclear data from a number of sources including updating spin and parity assignments. With updated spin and parity assignments, the levels which can contribute to the reaction rates are identified. The reaction rates are computed using a Monte Carlo method which has been used for previous evaluations of the reaction rates in order to focus solely on the changes due to modified nuclear data. The evaluated 22 Ne(α, γ) 26 Mg reaction rate remains substantially similar to that of Longland et al. but, including recent results from Texas A&M, the 22 Ne(α, n) 25 Mg reaction rate is lower at a range of astrophysically important temperatures. Stellar models computed with newton and mesa predict decreased production of the weak branch s process due to the decreased efficiency of 22 Ne as a neutron source. Using the new reaction rates in the mesa model results in 96 Zr / 94 Zr and 135 Ba / 136 Ba ratios in much better agreement with the measured ratios from presolar SiC grains. The 22 Ne + α reaction rates 22 Ne(α, γ) 26 Mg and 22 Ne(α, n) 25 Mg have been recalculated based on more recent nuclear data. The 22 Ne(α, γ) 26 Mg reaction rate remains substantially unchanged since the previous evaluation but the 22 Ne(α, γ) 25 Mg reaction rate is substantially decreased due to updated nuclear data. This results in significant changes to the nucleosynthesis in the weak branch of the s process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Ab initio calculation of the β-decay spectrum of 6 He

Here, we calculate the β spectrum in the decay of 6 He using quantum Monte Carlo methods with nuclear interactions derived from chiral effective field theory and consistent weak vector and axial currents. We work at second order in the multipole expansion, retaining terms suppressed by O(q 2 /m$^{2}_{π}$), where q denotes low-energy scales such as the reaction’s Q value or the electron energy, and m π is the pion mass. We go beyond the impulse approximation by including the effects of two-body vector and axial currents. We estimate the theoretical error on the spectrum by using four potential models in the Norfolk family of local two- and three-nucleon interactions, which have different cutoffs, fit two-nucleon data up to different energies, and use different observables to determine the couplings in the three-body force. We find the theoretical uncertainty on the β spectrum, normalized by the total rate, to be well below the permille level, and to receive contributions of comparable size from first- and second-order corrections in the multipole expansion. We consider corrections to the β decay spectrum induced by beyond-standard-model charged-current interactions in the standard model effective field theory, with and without sterile neutrinos, and discuss the sensitivity of the next generation of experiments to these interactions.

6 ≤ A ≤ 19↗
DOE OSTI2023

Revealing the strength of three-nucleon interactions with the proposed Einstein Telescope

Nuclear systems, ranging from atomic nuclei to dense matter probed in neutron stars, are governed by strong interactions. Three-nucleon forces have been found to be a crucial ingredient for the reliable description of these systems. Here, in this work, we explore how astrophysical data on neutron stars and their mergers from current and next-generation observatories will enable us to distinguish nuclear Hamiltonians. In particular, we investigate two different nuclear Hamiltonians that have been adjusted to reproduce two-nucleon scattering data and properties of light nuclei, but differ in the three-nucleon interactions among neutrons. We find that no significant constraints can be obtained from current data, but that the proposed Einstein Telescope could provide strong evidence to distinguish among these Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Nuclear three-body short-range correlations in coordinate space

Here, we study the effects of three-nucleon short-range correlations on nuclear coordinate-space densities. For this purpose, novel three-body densities are calculated for ground-state nuclei using the auxiliary-field diffusion Monte Carlo method. The results are analyzed in terms of the generalized contact formalism, extended to include three-body correlations, revealing the universal behavior of nucleon triplets at short distances. We identify the quantum numbers of such correlated triplets and extract scaling factors of triplet abundances that can be compared to upcoming inclusive electron-scattering data. We also show that the dynamics of these triplets is sensitive to three-body forces, and, therefore, the short-range part of three-body force models could be tested against appropriate exclusive electron-scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Quantum Monte Carlo calculations of dark matter scattering off light nuclei

We compute the matrix elements for elastic scattering of dark matter (DM) particles off light nuclei ( 2 H , 3 H , 3 He , 4 He , and 6 Li ) using quantum Monte Carlo methods. We focus on scalar-mediated DM-nucleus interactions and use scalar currents obtained to next-to-leading order in chiral effective theory. The nuclear ground states are obtained from a phenomenological nuclear Hamiltonian that includes the Argonne v 18 two-body interaction and the three-body Urbana IX interaction. Within this approach, we study the impact of one- and two-body currents and discuss the size of nuclear uncertainties, including for the first time two-body effects in A = 4 and A = 6 systems. Lastly, our results provide the nuclear structure input needed to assess the sensitivity of future experimental searches of (light) dark matter using light nuclei, such as 3 He and 4 He .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Doppler broadening in Mg 20 ( β p γ ) Ne 19 decay

The $^{15}$$\textbf{O}$$\Big({α,γ}\Big)$$^{19}$$\textbf{Ne}$ bottleneck reaction in Type I x-ray bursts is the most important thermonuclear reaction rate to constrain experimentally, to improve the accuracy of burst light-curve simulations. A proposed technique to determine the thermonuclear rate of this reaction employs the $^{20}$$\textbf{Mg}$$\Big({βpα}\Big)$$^{15}$$\textbf{O}$ decay sequence. The key $^{15}$$\textbf{O}$$\Big({α,γ}\Big)$$^{19}$$\textbf{Ne}$ resonance at an excitation of 4.03 MeV is now known to be fed in $^{20}$$\textbf{Mg}$$\Big({βpγ}\Big)$$^{19}$$\textbf{Ne}$ however, the energies of the protons feeding the 4.03 MeV state are unknown. Knowledge of the proton energies will facilitate future $^{20}$$\textbf{Mg}$$\Big({βpα}\Big)$$^{15}$$\textbf{O}$ measurements. To determine the energy of the proton transition feeding the 4.03 MeV state in $^{19}$$\textbf{Ne}$. A fast beam of $^{20}$$\textbf{Mg}$ was implanted into a plastic scintillator, which was used to detect $β$ particles. 16 high purity germanium detectors were used to detect $γ$ rays emitted following $βp$ decay. A Monte Carlo method was used to simulate the Doppler broadening of $^{19}$$\textbf{Ne}γ$-ray lines and compare to the experimental data. The center of mass energy between the proton and $^{19}$$\textbf{Ne}$, feeding the 4.03 MeV state, is measured to be $1.21^{+0.25}_{-0.22}$ $\textbf{MeV}$, corresponding to a $^{20}$$\textbf{Na}$ excitation energy of $7.44^{+0.25}_{-0.22}$$\textbf{MeV}$. Absolute feeding intensities and $γ$-decay branching ratios of $^{19}$$\textbf{Ne}$ states were determined including the 1615 keV state, which has not been observed before in this decay. A new $γ$ decay branch from the 1536 keV state in $^{19}$$\textbf{Ne}$ to the ground state is reported. The lifetime of the 1507 keV state in $^{19}$$\textbf{Ne}$ is measured to be $4.3^{+1.3}_{-1.1}$ ps resolving discrepancies in the literature. Conflicting $^{20}$$\textbf{Mg}$$\Big({βp}\Big)$ decay schemes in published literature are clarified.The utility of this Doppler broadening technique to provide information on $β$-delayed nucleon emission and excited-state lifetimes has been further demonstrated. In particular, knowledge of the proton energies feeding the 4.03 MeV $^{19}$$\textbf{Ne}$ state in $^{20}$$\textbf{Mg}$ $β$ decay will facilitate future measurements of the $α$-particle branching ratio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Full transport model of GW170817-like disk produces a blue kilonova

The 2017 detection of the inspiral and merger of two neutron stars in gravitational waves and gamma rays was accompanied by a quickly reddening transient. Such a transient was predicted to occur following a rapid neutron capture (r-process) nucleosynthesis event, which synthesizes neutron-rich, radioactive nuclei and can take place in both dynamical ejecta and in the wind driven off the accretion torus formed after a neutron star merger. In this study, we present the first three-dimensional general relativistic, full transport neutrino radiation magnetohydrodynamics simulations of the black hole-accretion disk-wind system produced by the GW170817 merger. We show that the small but non-negligible optical depths lead to neutrino transport globally coupling the disk electron fraction, which we capture by solving the transport equation with a Monte Carlo method. The resulting absorption drives up the electron fraction in a structured, continuous outflow, with electron fraction as high as Y e ~ 0.4 in the extreme polar region. Finally, we show via nuclear reaction network and radiative transfer calculations that nucleosynthesis in the disk wind will produce a blue kilonova.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Lattice QED in an external magnetic field: Evidence for dynamical chiral symmetry breaking

We simulate quantum electrodynamics (QED) in a strong constant homogeneous external magnetic field on a Euclidean space-time lattice using the rational hybrid Monte Carlo method, developed for simulating lattice quantum chromodynamics (QCD). Our primary goal is to measure the chiral condensate in the limit when the input electron mass m is zero. We observe a nonzero value, indicating that the external magnetic field catalyzes chiral symmetry breaking as predicted by approximate truncated Schwinger-Dyson methods. Such behavior is associated with dominance by the lowest Landau level which causes the effective dimensional reduction from 3 + 1 dimensions to 1 + 1 dimensions for charged particles (electrons and positrons) where the attractive forces of QED can produce chiral symmetry breaking with a dynamical electron mass and associated chiral condensate. Since our lattice simulations use bare (lattice) parameters, while the Schwinger-Dyson analyses work with renormalized quantities, direct numerical comparison will require renormalization of our lattice results. Published by the American Physical Society 2024

Kogut, J. B. (ORCID:000000030151548X)↗
DOE OSTI2015

Double-regge exchange limit for the γp→ K⁺K⁻p reaction

We apply the generalized Veneziano model (B₅ model) in the double-Regge exchange limit to the γp→K⁺K⁻p reaction. Four different cases defined by the possible combinations of the signature factors of leading Regge exchanges ((K * ,a₂/f₂), (K * ,ρ/ω), (K * ₂,a₂/f₂), and (K * ₂,ρ/ω)) have been simulated through the Monte Carlo method. Suitable event candidates for the double-Regge exchange high-energy limit were selected employing Van Hove plots as a better alternative to kinematical cuts in the K⁺K⁻p Dalitz plot. In this way we predict and analyze the double-Regge contribution to the K⁺K⁻p Dalitz plot, which constitutes one of the major backgrounds in the search for strangeonia, hybrids and exotics using γp→K⁺K⁻p reaction. We expect that data currently under analysis, and that to come in the future, will allow verification of the double-Regge behavior and a better assessment of this component of the amplitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Phenomenology of large scale structure in scalar-tensor theories: Joint prior covariance of w DE , Σ, and μ in Horndeski theories

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, w DE , and the phenomenology of large scale structure. They will allow us to constrain deviations from the Λ cold dark matter predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions Σ and μ, respectively. In this work, we derive the theoretical prior for the joint covariance of w DE , Σ and μ, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. Here, we identify several interesting features and trends in the distribution functions of w DE , Σ and μ, as well as in their covariances; we confirm the high degree of correlation between Σ and μ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly w DE , Σ and μ in a nonparametric way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Microemulsions in the driven Widom-Rowlinson lattice gas

An investigation of the two-dimensional Widom-Rowlinson lattice gas under an applied drive uncovered a remarkable nonequilibrium steady state in which uniform stripes (reminiscent of an equilibrium lamellar phase) form perpendicular to the drive direction [R. Dickman and R. K. P. Zia, Phys. Rev. E 97, 062126 (2018)]. Here we study this model at low particle densities in two and three dimensions, where we find a disordered phase with a characteristic length scale (a “microemulsion”) along the drive direction. Further we develop a continuum theory of this disordered phase to derive a coarse-grained field-theoretic action for the nonequilibrium dynamics. The action has the form of two coupled driven diffusive systems with different characteristic velocities, generated by an interplay between the particle repulsion and the drive. We then show how fluctuation corrections in the field theory may generate the characteristic features of the microemulsion phase, including a peak in the static structure factor corresponding to the characteristic length scale. This work lays the foundation for understanding the stripe phenomenon more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Solving the nuclear pairing model with neural network quantum states

In this work, we present a variational Monte Carlo method that solves the nuclear many-body problem in the occupation number formalism exploiting an artificial neural network representation of the groundstate wave function. A memory-efficient version of the stochastic reconfiguration algorithm is developed to train the network by minimizing the expectation value of the Hamiltonian. We benchmark this approach against widely used nuclear many-body methods by solving a model used to describe pairing in nuclei for different types of interaction and different values of the interaction strength. Despite its polynomial computational cost, our method outperforms coupled-cluster and provides energies that are in excellent agreement with the numerically-exact full configuration interaction values.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Application of the variational autoencoder to detect the critical points of the anisotropic Ising model

We generalize the previous study on the application of variational autoencoders to the two-dimensional Ising model to a system with anisotropy. Due to the self-duality property of the system, the critical points can be located exactly for the entire range of anisotropic coupling. This presents an excellent test bed for the validity of using a variational autoencoder to characterize an anisotropic classical model. Furthermore, we reproduce the phase diagram for a wide range of anisotropic couplings and temperatures via a variational autoencoder without the explicit construction of an order parameter. Considering that the partition function of ($d$ + 1)-dimensional anisotropic models can be mapped to that of the $d$-dimensional quantum spin models, the present study provides numerical evidence that a variational autoencoder can be applied to analyze quantum systems via the quantum Monte Carlo method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Magnetic phase transition in coupled spin-lattice systems: A replica-exchange Wang-Landau study

Coupled, dynamical spin-lattice models provide a unique test ground for simulations investigating the finite-temperature magnetic properties of materials under the direct influence of the lattice vibrations. These models are constructed by combining a coordinate-dependent interatomic potential with a Heisenberg-like spin Hamiltonian, facilitating the treatment of both the atomic coordinates and the spins as explicit phase variables. Using a model parameterized for bcc iron, we study the magnetic phase transition in these complex systems via the recently introduced, massively parallel replica-exchange Wang-Landau Monte Carlo method. Comparison with the results obtained from rigid lattice (spin-only) simulations shows that the transition temperature as well as the amplitude of the peak in the specific heat curve is marginally affected by the lattice vibrations. Moreover, we report the results were found to be sensitive to the particular choice of interatomic potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Quantum Monte Carlo with variable spins: Fixed-phase and fixed-node approximations

We study several aspects of the recently introduced fixed-phase spinor diffusion Monte Carlo method, in particular, its relation to the fixed-node method and its potential use as a general approach for electronic structure calculations. We illustrate constructions of spinor-based wave functions with the full space-spin symmetry without assigning up or down spin labels to particular electrons, effectively “complexifying” even ordinary real-valued wave functions for Hamiltonians without spin terms. Interestingly, with proper choice of the simulation parameters and spin variables, such fixed-phase calculations enable one to reach also the fixed-node limit. The fixed-phase approximation has several desirable properties when compared to the fixed-node approximation. The fixed-phase solution provides a straightforward interpretation as the lowest bosonic state in a given effective potential generated by the many-body approximate phase, whereas nodal boundary conditions are defined through less intuitive and complicated hypersurfaces with one dimension less than the original configuration space. In addition, the divergences of the local energy and drift at real wave function nodes are smoothed out to lower dimensionality when the wave function is complexified, thus decreasing the variation of sampled quantities and eliminating artificial nodal domain issues that can occur in the fixed-node formalism. We illustrate some of these properties on calculations of selected first-row systems that recover the fixed-node results with quantitatively similar levels of the corresponding biases. At the same time, the fixed-phase approach opens new possibilities for more general trial wave functions with further opportunities for increasing accuracy in practical calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Hypothesis testing of scientific Monte Carlo calculations

The steadily increasing size of scientific Monte Carlo simulations and the desire for robust, correct, and reproducible results necessitates rigorous testing procedures for scientific simulations in order to detect numerical problems and programming bugs. However, the testing paradigms developed for deterministic algorithms have proven to be ill suited for stochastic algorithms. In this paper we demonstrate explicitly how the technique of statistical hypothesis testing, which is in wide use in other fields of science, can be used to devise automatic and reliable tests for Monte Carlo methods, and we show that these tests are able to detect some of the common problems encountered in stochastic scientific simulations. We argue that hypothesis testing should become part of the standard testing toolkit for scientific simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗