Search NASA⌕ Search

SEARCH · Search NASA

Results for “Quantum 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Visualizing the impact of quenched disorder on 2D electron Wigner solids

Electron Wigner solids (WSs) provide an ideal system for understanding the competing effects of electron–electron and electron–disorder interactions, a central unsolved problem in condensed matter physics. Progress in this topic has been limited by a lack of single-defect-resolved experimental measurements as well as accurate theoretical tools to enable realistic experiment/theory comparison. Here, in this study, we overcome these limitations by combining atomically resolved scanning tunnelling microscopy (STM) with neural-quantum-state quantum Monte Carlo (NQS-QMC) simulation of disordered 2D electron WSs to discover new disorder-induced physical regimes of correlated electron behaviour. STM was used to image the electron density (n e )-dependent evolution of electron WSs in gate-tunable bilayer MoSe 2 (BL-MoSe 2 ) devices with varying long-range (n LR ) and short-range (n SR ) disorder densities. These images were compared with NQS-QMC simulations using realistic disorder maps extracted from experiment, thus allowing the roles of different disorder types to be disentangled. We identify two distinct physical regimes for disordered electron WSs that depend on n SR . For n SR ≲ n e , the WS behaviour is dominated by long-range disorder and features extensive mixed solid–liquid phases, a new type of local re-entrant melting/crystallization and prominent Friedel oscillations. By contrast, when n SR ≫ n e , these features are suppressed and a more robust amorphous WS phase emerges that persists to higher n e , highlighting the importance of short-range disorder in this regime. Our work establishes a powerful framework for studying disordered quantum solids through a combined experimental–theoretical approach.

Ge, Zhehao [University of California, Berkeley, CA↗
DOE OSTI2025

Kolmogorov-Arnold wavefunctions

Here, this work investigates Kolmogorov-Arnold network-based (KAN) wave-function Ansätz as viable representations for quantum Monte Carlo simulations. Through systematic analysis of one-dimensional model systems, we evaluate their computational efficiency and representational power against established methods. Our numerical experiments suggest some efficient training methods and we explore how the computational cost scales with desired precision, particle number, and system parameters. Roughly speaking, KANs seem to be 10 times cheaper computationally than other neural-network-based Ansätz . We also introduce a novel approach for handling strong short-range potentials—a persistent challenge for many numerical techniques—which generalizes efficiently to higher-dimensional, physically relevant systems with short-ranged strong potentials common in atomic and nuclear physics.

1-dimensional systems↗
DOE OSTI2022

DeePKS + ABACUS as a Bridge between Expensive Quantum Mechanical Models and Machine Learning Potentials

Recently, the development of machine learning (ML) potentials has made it possible to perform large-scale and long-time molecular simulations with the accuracy of quantum mechanical (QM) models. However, for different levels of QM methods, such as density functional theory (DFT) at the meta-GGA level and/or with exact exchange, quantum Monte Carlo, etc., generating a sufficient amount of data for training an ML potential has remained computationally challenging due to their high cost. In this work, we demonstrate that this issue can be largely alleviated with Deep Kohn–Sham (DeePKS), an ML-based DFT model. DeePKS employs a computationally efficient neural network-based functional model to construct a correction term added upon a cheap DFT model. Upon training, DeePKS offers closely matched energies and forces compared with high-level QM method, but the number of training data required is orders of magnitude less than that required for training a reliable ML potential. As such, DeePKS can serve as a bridge between expensive QM models and ML potentials: one can generate a decent amount of high-accuracy QM data to train a DeePKS model and then use the DeePKS model to label a much larger amount of configurations to train an ML potential. Further, this scheme for periodic systems is implemented in a DFT package ABACUS, which is open source and ready for use in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Characterizing GPU Energy Usage in Exascale-Ready Portable Science Applications

We characterize the GPU energy usage of two widely adopted exascale-ready applications representing two classes of particle and mesh solvers: (i) QMCPACK, a quantum Monte Carlo package, and (ii) AMReX-Castro, an adaptive mesh astrophysical code. We analyze power, temperature, utilization, and energy traces from double-/single (mixed)-precision benchmarks on NVIDIA’s A100 and H100 and AMD’s MI250X GPUs using queries in NVML and rocm_smi_lib, respectively. We explore application-specific metrics to provide insights on energy vs. performance trade-offs. Our results suggest that mixed-precision energy savings range between 6–25% on QMCPACK and 45% on AMReX-Castro. Also, we found gaps in the AMD tooling used on Frontier GPUs that need to be understood, while query resolutions on NVML have little variability between 1 ms-1 s. Overall, application level knowledge is crucial to define energy-cost/science-benefit opportunities for the codesign of future supercomputer architectures in the post-Moore era.

Godoy, William [ORNL] (ORCID:0000000225905178)↗
DOE OSTI2024

Software stewardship and advancement of a high-performance computing scientific application: QMCPACK

Here, we provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions own runners, and NVIDIA and AMD GPUs used in pre-exascale systems, (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Docker containers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the stewardship and advancement of community HPC codes to enable scientific discovery at scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗
DOE OSTI2020

Accurate Atomic Correlation and Total Energies for Correlation Consistent Effective Core Potentials

Very recently, we introduced a set of correlation consistent effective core potentials (ccECPs) constructed within full many-body approaches. By employing significantly more accurate correlated approaches, we were able to reach a new level of accuracy for the resulting effective core Hamiltonians. We also strived for simplicity of use and easy transferability into a variety of electronic structure methods in quantum chemistry and condensed matter physics. Here, as a reference for future use, we present exact or nearly exact total energy calculations for these ccECPs. The calculations cover H–Kr elements and are based on the state-of-the-art configuration interaction (CI), coupled-cluster (CC), and quantum Monte Carlo (QMC) calculations with systematically eliminated/improved errors. In particular, we carry out full CI/CCSD(T)/CCSDT(Q) calculations with cc-pVnZ with up to n = 6 basis sets and we estimate the complete basis set limits. Using combinations of these approaches, we achieved an accuracy of ≈1–10 mHa for K–Zn atoms and ≈0.1–0.3 mHa for all other elements—within about 1% or better of the ccECP total correlation energies. We also estimate the corresponding kinetic energies within the feasible limit of full CI calculations. In order to provide data for QMC calculations, we include fixed-node diffusion Monte Carlo energies for each element that give quantitative insights into the fixed-node biases for single-reference trial wave functions. The results offer a clear benchmark for future high-accuracy calculations in a broad variety of correlated wave function methods such as CI and CC as well is in stochastic approaches such as real space sampling QMC.

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

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

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

Evolution of spin excitations from bulk to monolayer FeSe

Abstract In ultrathin films of FeSe grown on SrTiO 3 (FeSe/STO), the superconducting transition temperature T c is increased by almost an order of magnitude, raising questions on the pairing mechanism. As in other superconductors, antiferromagnetic spin fluctuations have been proposed to mediate SC making it essential to study the evolution of the spin dynamics of FeSe from the bulk to the ultrathin limit. Here, we investigate the spin excitations in bulk and monolayer FeSe/STO using resonant inelastic x-ray scattering (RIXS) and quantum Monte Carlo (QMC) calculations. Despite the absence of long-range magnetic order, bulk FeSe displays dispersive magnetic excitations reminiscent of other Fe-pnictides. Conversely, the spin excitations in FeSe/STO are gapped, dispersionless, and significantly hardened relative to its bulk counterpart. By comparing our RIXS results with simulations of a bilayer Hubbard model, we connect the evolution of the spin excitations to the Fermiology of the two systems revealing a remarkable reconfiguration of spin excitations in FeSe/STO, essential to understand the role of spin fluctuations in the pairing mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Robust charge-density-wave correlations in the electron-doped single-band Hubbard model

There is growing evidence that the hole-doped single-band Hubbard and t–J models do not have a superconducting ground state reflective of the high-temperature cuprate superconductors but instead have striped spin- and charge-ordered ground states. Nevertheless, it is proposed that these models may still provide an effective low-energy model for electron-doped materials. Here we study the finite temperature spin and charge correlations in the electron-doped Hubbard model using quantum Monte Carlo dynamical cluster approximation calculations and contrast their behavior with those found on the hole-doped side of the phase diagram. We find evidence for a charge modulation with both checkerboard and unidirectional components decoupled from any spin-density modulations. These correlations are inconsistent with a weak-coupling description based on Fermi surface nesting, and their doping dependence agrees qualitatively with resonant inelastic x-ray scattering measurements. Our results provide evidence that the single-band Hubbard model describes the electron-doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Haldane topological spin-1 chains in a planar metal-organic framework

Haldane topological materials contain unique antiferromagnetic chains with symmetry-protected energy gaps. Such materials have potential applications in spintronics and future quantum computers. Haldane topological solids typically consist of spin-1 chains embedded in extended three-dimensional (3D) crystal structures. Here, we demonstrate that [Ni(μ-4,4'-bipyridine)(μ-oxalate)] n (NiBO) instead adopts a two-dimensional (2D) metal-organic framework (MOF) structure of Ni 2+ spin-1 chains weakly linked by 4,4'-bipyridine. NiBO exhibits Haldane topological properties with a gap between the singlet ground state and the triplet excited state. The latter is split by weak axial and rhombic anisotropies. Several experimental probes, including single-crystal X-ray diffraction, variable-temperature powder neutron diffraction (VT-PND), VT inelastic neutron scattering (VT-INS), DC susceptibility and specific heat measurements, high-field electron spin resonance, and unbiased quantum Monte Carlo simulations, provide a detailed, comprehensive characterization of NiBO. Vibrational (also known as phonon) properties of NiBO have been probed by INS and density-functional theory (DFT) calculations, indicating the absence of phonons near magnetic excitations in NiBO, suppressing spin-phonon coupling. The work here demonstrates that NiBO is indeed a rare 2D-MOF Haldane topological material.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗
DOE OSTI2020

The range of non-Kitaev terms and fractional particles in α-RuCl 3

Significant efforts have focused on the magnetic excitations of relativistic Mott insulators, predicted to realize the Kitaev quantum spin liquid (QSL). This exactly solvable model involves a highly entangled state resulting from bond-dependent Ising interactions that produce excitations which are non-local in terms of spin flips. A key challenge in real materials is identifying the relative size of the non-Kitaev terms and their role in the emergence or suppression of fractional excitations. Here, we identify the energy and temperature boundaries of non-Kitaev interactions by direct comparison of the Raman susceptibility of α-RuCl 3 with quantum Monte Carlo (QMC) results for the Kitaev QSL. Moreover, we further confirm the fractional nature of the magnetic excitations, which is given by creating a pair of fermionic quasiparticles. Interestingly, this fermionic response remains valid in the non-Kitaev range. Our results and focus on the use of the Raman susceptibility provide a stringent new test for future theoretical and experimental studies of QSLs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Orbital structure of the effective pairing interaction in the high-temperature superconducting cuprates

The nature of the effective interaction responsible for pairing in the high-temperature superconducting cuprates remains unsettled. This question has been studied extensively using the simplified single-band Hubbard model, which does not explicitly consider the orbital degrees of freedom of the relevant CuO 2 planes. Here, we use a dynamical cluster quantum Monte Carlo approximation to study the orbital structure of the pairing interaction in the three-band Hubbard model, which treats the orbital degrees of freedom explicitly. We find that the interaction predominately acts between neighboring copper orbitals, but with significant additional weight appearing on the surrounding bonding molecular oxygen orbitals. By explicitly comparing these results to those from the simpler single-band Hubbard model, our study provides strong support for the single-band framework for describing superconductivity in the cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The role of electron correlations in the electronic structure of putative Chern magnet TbMn6Sn6

Abstract A member of the RMn 6 Sn 6 rare-earth family materials, TbMn 6 Sn 6 , recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard U (DFT + U ) calculations to examine the electronic structure of TbMn 6 Sn 6 . To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration–interaction (CI) methods. We find that DFT + U and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT + U with Mn magnetic moments adjusted to the experiment predict the Dirac crossing in bulk to be close to the Fermi level, within ~120 meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Emergence of antiferromagnetic correlations and Kondolike features in a model for infinite layer nickelates

We report a determinant quantum Monte Carlo study of a two-band model, inspired by infinite-layer nickelates, focusing on the influence of interlayer hybridization between $3d_{x^{2}-y^{2}}$ orbitals derived from Ni (or Ni and O) in one layer and rare-earth (R) 5d orbitals in the other layer, hereafter the Ni and R layers, respectively. For a filling with one electron shared between the two layers on average, interlayer hybridization leads to “self-doped" holes in the Ni layer and the absence of antiferromagnetic ordering, but rather the appearance of spin-density and charge-density stripe-like states. As the interlayer hybridization increases, both the Ni and R layers develop antiferromagnetic correlations, even though either layer individually remains away from half-filling. For hybridization within an intermediate range, roughly comparable to the intralayer nearest-neighbor hopping t Ni , the model develops signatures of Kondo-like physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Interlayer pairing in bilayer nickelates

The discovery of T c ~ 80 K superconductivity in pressurized La 3 Ni 2 O 7 has launched a new platform to study high-temperature superconductivity. Using non-perturbative dynamic cluster approximation quantum Monte Carlo calculations, we characterize the magnetic and superconducting pairing behavior of a realistic bilayer two-orbital Hubbard-Hund model of this system that describes the relevant Ni e g states with physically relevant interaction strengths. We find a leading s ± superconducting instability in this model at a temperature T ~ 100 K close to the experimentally observed T c . Analyzing the orbital and spatial structure of the effective pairing interaction giving rise to this state reveals that the interaction predominantly acts between local interlayer pairs of the $d_{3z^2-r^2}$ orbital. By correlating the strength of the interaction with that of the magnetic spin fluctuations we show that it is driven by strong interlayer spin-fluctuations arising from the $d_{3z^2-r^2}$ orbital. These results provide first-time non-perturbative evidence supporting the picture that a simple single-orbital bilayer Hubbard model for the Ni $d_{3z^2-r^2}$ orbital provides an excellent low-energy effective description of the superconducting behavior of La 3 Ni 2 O 7 .

Maier, Thomas A. [Oak Ridge National Laboratory (O↗
DOE OSTI2020

Many-body factorization and position–momentum equivalence of nuclear short-range correlations

While mean-field approximations, such as the nuclear shell model, provide a good description of many bulk nuclear properties, they fail to capture the important effects of nucleon–nucleon correlations such as the short-distance and high-momentum components of the nuclear many-body wave function1. Here, we study these components using the effective pair-based generalized contact formalism and ab initio quantum Monte Carlo calculations of nuclei from deuteron to 40 Ca. We observe a universal factorization of the many-body nuclear wave function at short distance into a strongly interacting pair and a weakly interacting residual system. The residual system distribution is consistent with that of an uncorrelated system, showing that short-distance correlation effects are predominantly embedded in two-body correlations. Spin- and isospin-dependent ‘nuclear contact terms’ are extracted in both coordinate and momentum space for different realistic nuclear potentials. The contact coefficient ratio between two different nuclei shows very little dependence on the nuclear interaction model. These findings thus allow extending the application of mean-field approximations to short-range correlated pair formation by showing that the relative abundance of short-range pairs in the nucleus is a long-range (that is, mean field) quantity that is insensitive to the short-distance nature of the nuclear force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗