Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational Monte Carlo”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Static Self-Energy and Effective Mass of the Homogeneous Electron Gas from Quantum Monte Carlo Calculations

We discuss the methodology of quantum Monte Carlo calculations of the effective mass based on the static self-energy Σ(k,0). We then use variational Monte Carlo calculations of Σ(k,0) of the homogeneous electron gas at various densities to obtain results very close to perturbative G 0 W 0 calculations for values of the density parameter 1 ≤ r s ≤ 10. Here, the obtained values for the effective mass are close to diagrammatic Monte Carlo results and disagree with previous quantum Monte Carlo calculations based on a heuristic mapping of excitation energies to those of an ideal gas.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Method-independent cusps for atomic orbitals in quantum Monte Carlo

Here, we present an approach for augmenting Gaussian atomic orbitals with correct nuclear cusps. Like the atomic orbital basis set itself and unlike previous cusp corrections, this approach is independent of the many-body method used to prepare wave functions for quantum Monte Carlo. Once the basis set and molecular geometry are specified, the cusp-corrected atomic orbitals are uniquely specified, regardless of which density functionals, quantum chemistry methods, or subsequent variational Monte Carlo optimizations are employed. We analyze the statistical improvement offered by these cusps in a number of molecules and find them to offer similar advantages as molecular-orbital-based approaches while remaining independent of the choice of many-body method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo

Abstract We present a discussion of recent progress in excited‐state‐specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited‐state‐specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited‐state‐specific complete active space self‐consistent field theory, selected configuration interaction, and state‐specific variance minimization. These demonstrations show that combining excited‐state‐specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi‐reference perturbation theory, and that it can offer new clarity in cases where existing high‐level methods do not agree. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

Otis, Leon↗
NASA NTRS1986

Cometary atmospheres: Modeling the spatial distribution of observed neutral radicals

Progress during the second year of a program of research on the modeling of the spatial distributions of cometary radicals is discussed herein in several major areas. New scale length laws for cometary C2 and CN were determined which explain that the previously-held apparent drop of the C2/CN ratio for large heliocentric distances does not exist and that there is no systematic variation. Monte Carlo particle trajectory model (MCPTM) analysis of sunward and anti-sunward brightness profiles of cometary C2 was completed. This analysis implies a lifetime of 31,000 seconds for the C2 parent and an ejection speed for C2 of approximately 0.5 km/sec upon dissociation from the parent. A systematic reanalysis of published C3 and OH data was begun. Preliminary results find a heliocentric distance dependence for C3 scale lengths with a much larger variation than for C2 and CN. Scale lengths for OH are generally somewhat larger than currently accepted values. The MCPTM was updated to include the coma temperature. Finally, the collaborative effort with the University of Arizona programs has yielded some preliminary CCD images of Comet P/Halley.

Combi, Michael R.↗
DOE OSTI2026

Dirac Node Pinning from Dzyaloshinskii-Moriya Interactions in a Kagome Spin Liquid

Recent experiments on the Kagome spin liquid candidate material YCu 3 ⁢(OH) 6 ⁢Br 2 ⁢[Br 1−𝑦 ⁢(OH) 𝑦 ] suggest the presence of Dirac fermionic spinons near the magnetization plateau at 1/9. Theories suggest that the spinons are charge neutral spin-1/2 excitations, in a 2⁢𝜋/3 flux, which triples the unit cell. Generally a gap is expected, and there is no symmetry protection for the Dirac nodes in this system. The question arises as to what causes the nodes and stabilizes them. In this work, we propose a node-creation and node-pinning mechanism driven by Dzyaloshinskii-Moriya (DM) interactions. Employing Gutzwiller-projected variational Monte Carlo calculations, we demonstrate that DM interactions induce a band closing phase transition in the spinon spectrum. There is a change in the Chern number when the bands are inverted. Together with the DM-generated internal gauge flux, the coupling to the spinon orbital magnetization counteracts the band reopening. Furthermore, this interplay energetically pins the Dirac nodes over a range of parameters, resulting in a pinning mechanism distinct from the usual one from symmetry protection.

Kagome lattice↗
DOE OSTI2024

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

SchrödingerNet: A Universal Neural Network Solver for the Schrödinger Equation

Recent advances in machine learning have facilitated numerically accurate solution of the electronic Schrödinger equation (SE) by integrating various neural network (NN)-based wave function ansatzes with variational Monte Carlo methods. Nevertheless, such NN-based methods are all based on the Born–Oppenheimer approximation (BOA) and require computationally expensive training for each nuclear configuration. In this work, we propose a novel NN architecture, SchrödingerNet, to solve the full electronic-nuclear SE by defining a loss function designed to equalize local energies across the system. This approach is based on a translationally, rotationally and permutationally symmetry-adapted total wave function ansatz that includes both nuclear and electronic coordinates. Furthermore, this strategy not only allows for an efficient and accurate generation of a continuous potential energy surface at any geometry within the well-sampled nuclear configuration space, but also incorporates non-BOA corrections, through a single training process. Comparison with benchmarks of atomic and small molecular systems demonstrates its accuracy and efficiency.

Chemical calculations↗
DOE OSTI2022

Stripe correlations in the two-dimensional Hubbard-Holstein model

Several state-of-the-art numerical methods have observed static or fluctuating spin and charge stripes in doped two-dimensional Hubbard models, suggesting that these orders play a significant role in shaping the cuprate phase diagram. Many experiments, however, also indicate that the cuprates have strong electron-phonon (e-ph) coupling, and it is unclear how this interaction influences stripe correlations. We study static and fluctuating stripe orders in the doped single-band Hubbard-Holstein model using zero temperature variational Monte Carlo and finite temperature determinant quantum Monte Carlo. We find that the lattice couples more strongly with the charge component of the stripes, leading to an enhancement or suppression of stripe correlations, depending on model parameters like the next-nearest-neighbor hopping t' or phonon energy Ω. Our results help elucidate how the e-ph interaction can tip the delicate balance between stripe and superconducting correlations in the Hubbard-Holstein model with implications for our understanding of the high-T c cuprates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
DOE OSTI2026

Three-nucleon lepton-number-violating potentials in chiral effective field theory and their matrix elements in light nuclei

Here, we derive the three-nucleon neutrinoless double-𝛽 decay potential in a Δ-full chiral effective field theory through next-to-next-to-next-to leading order in Weinberg's power counting. The matrix elements of the resulting operators are computed in light nuclei using variational Monte Carlo with wave functions constructed from the Norfolk family of nuclear interactions. We find that three-nucleon corrections induce a modest quenching of the total nuclear matrix elements. We discuss model dependencies and the potential impact of these corrections on the sensitivity of experimental programs to probe lepton number violating parameters. These results provide a benchmark for many-body methods capable of reaching heavier nuclei of experimental interest.

Chambers-Wall, Graham [Washington University, St. ↗
DOE OSTI2022

Machine learning predictions for local electronic properties of disordered correlated electron systems

We present a scalable machine learning (ML) model to predict local electronic properties such as on-site electron number and double occupation for disordered correlated electron systems. Our approach is based on the locality principle, or the nearsightedness nature, of many-electron systems, which means local electronic properties depend mainly on the immediate environment. A ML model is developed to encode this complex dependence of local quantities on the neighborhood. We demonstrate our approach using the square-lattice Anderson-Hubbard model, which is a paradigmatic system for studying the interplay between Mott transition and Anderson localization. We develop a lattice descriptor based on the group-theoretical method to represent the on-site random potentials within a finite region. The resultant feature variables are used as input to a multilayer fully connected neural network, which is trained from data sets of variational Monte Carlo (VMC) simulations on small systems. We show that the ML predictions agree reasonably well with the VMC data. Our work underscores the promising potential of ML methods for multiscale modeling of correlated electron systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing spin-isospin excitations in proton-rich nuclei via the C 11 ( p , n ) N 11 reaction

Tracking the evolution of nuclear properties away from stability serves as a valuable test for nuclear models. In the present work, the (p,n) charge-exchange reaction was used to test the extraction of β- Gamow-Teller transition strengths, B(GT), from proton-rich unstable isotopes, and the resulting B(GT) values were compared to shell-model and ab-initio calculations. The 11 C(p,n) 11 N reaction was measured in inverse kinematics at 95 MeV/u at the National Superconducting Cyclotron Laboratory (NSCL). The B(GT) values to the $\frac{1^-}{2}$ state at 0.73 MeV and the $\frac{3^-}{2}$ state at 2.86 MeV in 11 N were determined to be 0.18(1) stat (3) sys and 0.18(1) stat (4) sys , respectively. These results are consistent with shell-model calculations using the wbp interaction after introducing a phenomenological quenching factor and with ab initio variational Monte Carlo calculations using the NV2 + 3Ia * NN and 3N interactions without any scaling. Additionally, this result is consistent with the B(GT) values extracted from mirror 11 B(n,p) and 11 B(t, 3 He) reactions. In conclusion, this experiment demonstrates the feasibility of using the (p,n) probe in inverse kinematics to extract B(GT) from proton-rich nuclei, although improved background suppression will be important in future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Densities and momentum distributions in A ≤ 12 nuclei from chiral effective field theory interactions

Current and future electron and neutrino scattering experiments will be greatly aided by a better understanding of the role played by short-range correlations in nuclei. Two-body physics, including nucleon-nucleon correlations and two-body electroweak currents, is required to explain the body of experimental data for both static and dynamical nuclear properties. In this work, we focus on examining nucleon-nucleon correlations from a chiral effective field theory perspective and provide a comprehensive set of new variational Monte Carlo calculations of one- and two-body densities and momentum distributions based on the Norfolk many-body nuclear Hamiltonians for A ≤ 12 systems. Here, online access to detailed tables and figures is available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Longitudinal form Factors of A ≤ 10 Nuclei in a Chiral Effective Field Theory Approach

Here, in this work, we present the elastic electron scattering longitudinal form factors of A ≤ 10 nuclei computed in a variational Monte Carlo approach. We employ the Norfolk family of local chiral interactions and a consistent electromagnetic charge operator. Our calculations are compared both to data and past theoretical evaluations. This work represents, to our knowledge, the first accurate calculation of longitudinal form factors using manybody methods based on interacting nucleon degrees of freedom in the 7 ≤ A ≤ 10 mass range. Finally, we identify 9 Be and 10 B as candidate targets for renewed experimental interest, as they exhibit the potential to provide more stringent constraints on the theoretical models.

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

Distilling the Essential Elements of Nuclear Binding via Neural-Network Quantum States

To distill the essential elements of nuclear binding, we seek the simplest Hamiltonian capable of modeling atomic nuclei with percent-level accuracy. A critical aspect of this endeavor consists of accurately solving the quantum many-body problem without incurring an exponential computing cost with the number of nucleons. Furthermore, we address this challenge by leveraging a variational Monte Carlo method based on a highly expressive neural-network quantum state ansatz. In addition to computing binding energies and charge radii of nuclei with up to 𝐴 = 20 nucleons, by evaluating their magnetic moments, we demonstrate that neural-network quantum states are able to correctly capture the self-emerging nuclear shell structure. To this end, we introduce a novel computational protocol based on adding an external magnetic field to the nuclear Hamiltonian, which allows the neural network to learn the preferred polarization of the nucleus within the given magnetic field.

Binding energy & masses↗
DOE OSTI2026

Nuclear Responses with Neural-Network Quantum States

We introduce a variational Monte Carlo framework that combines neural-network quantum states with the Lorentz integral transform technique to compute the dynamical properties of self-bound quantum many-body systems in continuous Hilbert spaces. While broadly applicable to various quantum systems, including atoms and molecules, in this initial application we focus on the photoabsorption cross section of light nuclei, where benchmarks against numerically exact techniques are available. Our accurate theoretical predictions are complemented by robust uncertainty quantification, enabling meaningful comparisons with experiments. Here, we demonstrate that a relatively simple nuclear Hamiltonian—based on a leading-order pionless EFT expansion and known to accurately reproduce ground-state energies of nuclei with 𝐴 ≤ 40—also provides a reliable description of the photoabsorption cross section.

Ab initio calculations↗