SEARCH · Engineering Papers

Results for “Monte Carlo method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Quantum advantage for differential equation analysis

Quantum algorithms for differential equation solving, data processing, and machine learning potentially offer an exponential speedup over all known classical algorithms. However, there also exist obstacles to obtaining this potential speedup in useful problem instances. The essential obstacle for quantum differential equation solving is that outputting useful information may require difficult postprocessing, and the essential obstacle for quantum data processing and machine learning is that inputting the data is a difficult task just by itself. In this study, we demonstrate that, when combined, these difficulties solve one another. We show how the output of quantum differential equation solving can serve as the input for quantum data processing and machine learning, allowing dynamical analysis in terms of principal components, power spectra, and wavelet decompositions. To illustrate this, we consider continuous-time Markov processes on epidemiological and social networks. These quantum algorithms provide an exponential advantage over existing classical Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Charge density wave order on a $\textit{π}$-flux square lattice

The effect of electron-phonon coupling on Dirac fermions was recently explored numerically on a honeycomb lattice, leading to precise quantitative values for the finite temperature and quantum critical points. In this paper, we use the unbiased determinant quantum Monte Carlo method to study the Holstein model on a half-filled staggered-flux square lattice and compare the results with the honeycomb lattice geometry, presenting results for a range of phonon frequencies 0.1 ≤ $\omega$ ≤ 2.0. In this work, we find that the interactions give rise to charge density wave order, but only above a finite coupling strength $\gamma_{\text{crit}}$. The transition temperature is evaluated and presented in a $T_{c^{–λ}}$ phase diagram. An accompanying mean-field theory calculation also predicts the existence of a quantum phase transition, but at a substantially smaller coupling strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Witnessing entanglement in quantum magnets using neutron scattering

We demonstrate how quantum entanglement can be directly witnessed in the quasi-1D Heisenberg antiferromagnet KCuF 3 . We apply three entanglement witnesses—one tangle, two tangle, and quantum Fisher information—to its inelastic neutron spectrum and compare with spectra simulated by finite-temperature density matrix renormalization group (DMRG) and classical Monte Carlo methods. We find that each witness provides direct access to entanglement. Of these, quantum Fisher information is the most robust experimentally and indicates the presence of at least bipartite entanglement up to at least 50 K, corresponding to around 10% of the spinon zone-boundary energy. We apply quantum Fisher information to higher spin-S Heisenberg chains and show theoretically that the witnessable entanglement gets suppressed to lower temperatures as the quantum number increases. Finally, we outline how these results can be applied to higher dimensional quantum materials to witness and quantify entanglement.

1-dimensional systems↗

DOE OSTI2024

Magnetic, charge, and bond order in the two-dimensional Su-Schrieffer-Heeger-Holstein model

Most nonperturbative numerical studies of electron-phonon interactions focus on model Hamiltonians where the electrons interact with a phonon branch via a single type of microscopic mechanism. Two commonly explored couplings in this context are the Holstein and Su-Schrieffer-Heeger (SSH) interactions, which describe phonons modulating the on-site energy and intersite electron hopping, respectively. Many materials, however, have multiple phonon branches that can each interact with electronic degrees of freedom in different ways. We present here a determinant quantum Monte Carlo study of the half-filled two-dimensional (bond) SSH-Holstein Hamiltonian, where electrons couple to different phonon branches via either the Holstein or SSH mechanism. As a result, we map the model's phase diagram and determine the nature of the transitions between charge-density wave, bond-order wave, and antiferromagnetic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Indirect determination of neutron-capture cross sections for Sm isotopes

Neutron-capture cross sections in the energy range of 0.01–10 MeV for Sm 146 , 148 , 150 , 152 nuclei were calculated using the γ -ray strength and level density functions extracted from the Oslo-type experiments. The uncertainties in the cross-section values were determined using a Monte Carlo method. For the Sm 148 , 150 , 152 isotopes, the calculated cross sections are in a good agreement with the existing experimental data and for the Sm 146 nucleus, an experimental ( n , γ ) cross section is reported for the first time. The results are compared with the ENDF, EAF, and TENDL evaluations. Maxwellian-averaged cross sections were also calculated using the same input γ SF and LD functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Charge conservation and higher moments of charge fluctuations

Higher moments of distributions of net charge and baryon number in heavy-ion collisions have been proposed as signals of fundamental QCD phase transitions. In order to better understand background processes for these observables, models are presented which enable one to gauge the effects of local charge conservation, decays of resonances and clusters, Bose symmetrization, and volume fluctuations. Monte Carlo methods for generating samplings of particles consistent with local charge conservation are presented and are followed by a review of simple analytic models involving a single type of charge with a constant experimental efficiency. The main model consists of thermal emission superimposed onto a simple parametrization of collective flow, known as a blast wave, with emission being consistent with individual canonical ensembles. The spatial extent of local charge conservation is parameterized by the size and extent over which charge is conserved. Here, the sensitivity of third- and fourth-order moments, skewness and kurtosis, to these parameters, and to beam energy and baryon density is explored. Comparisons with STAR data show that a significant part of the observed non-Poissonian fluctuations in net-proton fluctuations are explained by charge and baryon-number conservation, but that measurements of the STAR collaboration for fluctuations of net electric charge significantly differ from expectations of the models presented here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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 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 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 OSTI2021

Effect of slip on vortex shedding from a circular cylinder in a gas flow

Most studies of vortex shedding from a circular cylinder in a gas flow have explicitly or implicitly assumed that the no-slip condition applies on the cylinder surface. To investigate the effect of slip, vortex shedding is simulated using molecular gas dynamics (the direct simulation Monte Carlo method) and computational fluid dynamics (the incompressible Navier-Stokes equations with a slip boundary condition). A Reynolds number of 100, a Mach number of 0.3, and a corresponding Knudsen number of 0.0048 are examined. For these conditions, compressibility effects are small, and periodic laminar vortex shedding is obtained. Slip on the cylinder is varied using combinations of diffuse and specular molecular reflections with accommodation coefficients from zero (maximum slip) to unity (minimum slip). Although unrealistic, bounce-back molecular reflections are also examined because they approximate the no-slip boundary condition (zero slip). The results from both methods are in reasonable agreement. The shedding frequency increases slightly as the accommodation coefficient is decreased, and shedding ceases at low accommodation coefficients (large slip). The streamwise and transverse forces decrease as the accommodation coefficient is decreased. Based on the good agreement between the two methods, computational fluid dynamics is used to determine the critical accommodation coefficient below which vortex shedding ceases for Reynolds numbers of 60–100 at a Mach number of 0.3. Finally, conditions to observe the effect of slip on vortex shedding appear to be experimentally realizable, although challenging.

42 ENGINEERING↗

DOE OSTI2020

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

Here we present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Ab Initio Nuclear Thermodynamics

Here, we propose a new Monte Carlo method called the pinhole trace algorithm for ab initio calculations of the thermodynamics of nuclear systems. For typical simulations of interest, the computational speedup relative to conventional grand-canonical ensemble calculations can be as large as a factor of one thousand. Using a leading-order effective interaction that reproduces the properties of many atomic nuclei and neutron matter to a few percent accuracy, we determine the location of the critical point and the liquid-vapor coexistence line for symmetric nuclear matter with equal numbers of protons and neutrons. We also present the first ab initio study of the density and temperature dependence of nuclear clustering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗