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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Neural-network quantum states for the nuclear many-body problem

A long-standing goal of nuclear theory is to explain how the structure and dynamics of atomic nuclei and neutron-star matter emerge from the underlying interactions among protons and neutrons. Achieving this goal requires solving the nuclear quantum many-body problem with high accuracy across a wide range of length scales and density regimes. In this review, we discuss how artificial neural network representations of the nuclear many-body wave function have significantly extended the capabilities of continuum quantum Monte Carlo methods. In particular, neural network quantum states enable calculations of larger systems than were previously accessible and provide a flexible framework for capturing phenomena that challenge conventional approaches, including the emergence of nuclear clusters and superfluid phases in dense matter. We highlight recent applications to finite nuclei, infinite nuclear and neutron matter, and dynamical processes relevant to lepton-nucleus and nucleus-nucleus scattering. We also discuss conceptual and methodological connections with condensed matter physics, emphasizing developments in neural network quantum states that bridge strongly correlated systems across disciplines. Together, these developments demonstrate how neural-network methods open new avenues toward unified and accurate descriptions of nuclear structure, matter, and reactions.

Lovato, Alessandro [Argonne; TIFPA-INFN, Trento; V↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Langevin simulations of the half-filled cubic Holstein model

Over the past several years, reliable quantum Monte Carlo results for the charge density wave transition temperature T cdw of the half-filled two-dimensional Holstein model in square and honeycomb lattices have become available for the first time. Here, exploiting the further development of numerical methodology, here we present results in three dimensions, which are made possible through the use of Langevin evolution of the quantum phonon degrees of freedom. In addition to determining T cdw from the scaling of the charge correlations, we also examine the nature of charge order at general wave vectors for different temperatures, couplings, and phonon frequencies, and the behavior of the spectral function and specific heat.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Thermodynamic Modeling of Intrinsic Defects in MnBi₂Te₄

This repository contains the computational data supporting the manuscript titled “The critical role of intrinsic defects and many-body interactions on the stability of MnBi₂Te₄.” It includes: 1. DFT data generated using VASP, used for training and benchmarking electronic structure models. 2. Quantum Monte Carlo (QMC) data produced with QMCPACK, used to apply many-body corrections and validate the electronic and magnetic properties of MnBi₂Te₄. 3. Relevant scripts used to run, analyze, and process the calculations, enabling reproducibility and transparency of the workflows.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase diagram of the Su-Schrieffer-Heeger-Hubbard model on a square lattice

The Hubbard and Su-Schrieffer-Heeger (SSH) Hamiltonians are iconic models for understanding the qualitative effects of electron-electron and electron-phonon interactions, respectively. In the two-dimensional square lattice Hubbard model at half filling, the on-site Coulomb repulsion U between up and down electrons induces antiferromagnetic (AFM) order and a Mott insulating phase. On the other hand, for the SSH model, there is an AFM phase when the electron-phonon coupling λ is less than a critical value λ c and a bond order wave when λ > λ c . In this Letter, we perform numerical studies on the square lattice optical Su-Schrieffer-Heeger-Hubbard Hamiltonian, which combines both interactions. We use the determinant quantum Monte Carlo method which does not suffer from the fermionic sign problem at half filling. We map out the phase diagram and find that it exhibits a direct first-order transition between an antiferromagnetic phase and a bond-ordered wave as λ increases. Here, the AFM phase is characterized by two different regions. At smaller λ the behavior is similar to that of the pure Hubbard model; the other region, while maintaining long-range AFM order, exhibits larger kinetic energies and double occupancy, i.e., larger quantum fluctuations, similar to the AFM phase found in the pure SSH model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Charge density wave and superconductivity in the disordered Holstein model

The interplay between electron-electron correlations and disorder has been a central theme of condensed matter physics over the past several decades, with particular interest in the possibility that interactions might cause delocalization of an Anderson insulator into a metallic state, and the disrupting effects of randomness on magnetic order and the Mott phase. In this work, we extend this physics to explore electron-phonon interactions and show, via exact quantum Monte Carlo simulations, that the suppression of the charge density wave correlations in the half-filled Holstein model by disorder can stabilize a superconducting phase. Our simulations thus capture qualitatively the suppression of charge ordered phases and emergent superconductivity recently seen experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

π phase shift across stripes in a charge density wave system

Many strongly correlated materials are characterized by deeply intertwined charge and spin order. Besides their high superconducting transition temperatures, one of the central features of these complex patterns in cuprates is a phase shift which occurs across lines of decreased hole density. That is, when doped away from their antiferromagnetic phase, the additional charge is not distributed uniformly, but rather in “stripes.” The sublattices preferentially occupied by up and down spin are reversed across these stripes, a phenomenon referred to as a “π phase shift.” Many of the spin-charge patterns, including the π phase shift, are reproduced by density matrix renormalization group and quantum Monte Carlo calculations of simplified tight binding (repulsive Hubbard) models. In this paper we demonstrate that this sublattice reversal is generic by considering the corresponding phenomenon in the attractive Hubbard Hamiltonian, where a charge density wave phase forms at half filling. Here, we introduce charge stripes via an appropriate local chemical potential; measurements of charge correlation across the resulting lines of lowered density reveal a clear π phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantum membrane phases in synthetic lattices of cold molecules or Rydberg atoms

Here, we calculate properties of dipolar interacting ultracold molecules or Rydberg atoms in a semisynthetic three-dimensional configuration—one synthetic dimension plus a two-dimensional real-space optical lattice or periodic microtrap array—using the stochastic Green's function quantum Monte Carlo method. Through a calculation of thermodynamic quantities and appropriate correlation functions, along with their finite-size scalings, we show that there is a second-order transition to a low-temperature phase in which two-dimensional “sheets” form in the synthetic dimension of internal rotational or electronic states of the molecules or Rydberg atoms, respectively. Simulations for different values of the interaction V, which acts between atoms or molecules that are adjacent both in real and synthetic space, allow us to compute a phase diagram. We find a finite-temperature transition at sufficiently large V as well as a quantum phase transition—a critical value V c below which the transition temperature vanishes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Modeling Unconventional Superconductivity at the Crossover between Strong and Weak Electronic Interactions

High-temperature superconductivity emerges in many different quantum materials, often in regions of the phase diagram where the electronic kinetic energy is comparable to the electron-electron repulsion. Describing such intermediate-coupling regimes has proven challenging as standard perturbative approaches are inapplicable. In this study, we employ quantum Monte Carlo methods to solve a multiband Hubbard model that does not suffer from the sign problem and in which only repulsive interband interactions are present. In contrast to previous sign-problem-free studies, we treat magnetic, superconducting, and charge degrees of freedom on an equal footing. We find an antiferromagnetic dome accompanied by a metal-to-insulator crossover line in the intermediate-coupling regime, with a smaller superconducting dome appearing in the metallic region. Across the antiferromagnetic dome, the magnetic fluctuations change from overdamped in the metallic region to propagating in the insulating region. Our findings shed new light on the intertwining between superconductivity, magnetism, and charge correlations in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Z 2 topological order and first-order quantum phase transitions in systems with combinatorial gauge symmetry

We study a generalization of the two-dimensional transverse-field Ising model, combining both ferromagnetic and antiferromagnetic two-body interactions, that hosts exact global and local Z 2 gauge symmetries. Using exact diagonalization and stochastic series expansion quantum Monte Carlo methods, we confirm the existence of the topological phase in line with previous theoretical predictions. Our simulation results show that the transition between the confined topological phase and the deconfined paramagnetic phase is of first order, in contrast to the conventional Z 2 lattice gauge model in which the transition maps onto that of the standard Ising model and is continuous. We further generalize the model by replacing the transverse field on the gauge spins with a ferromagnetic X X interaction while keeping the local gauge symmetry intact. We find that the Z 2 topological phase remains stable, while the paramagnetic phase is replaced by a ferromagnetic phase. The topological-ferromagnetic quantum phase transition is also of first order. For both models, we discuss the low-energy spinon and vison excitations of the topological phase and their avoided level crossings associated with the first-order quantum phase transitions.

2-dimensional systems↗
DOE OSTI2023

Magnetic and singlet phases in the three-dimensional periodic Anderson Model

Heavy-fermion materials are compounds in which localized f orbitals hybridize with delocalized d ones, leading to quasiparticles with large renormalized masses. The presence of strongly correlated f electrons at the Fermi level may also lead to long-range order, such as magnetism, or unconventional superconductivity. From a theoretical point of view, the “standard model” for heavy-fermion compounds is the periodic Anderson model (PAM). Despite being extensively scrutinized, its thermodynamic properties in three-dimensional (3D) lattices have not been carefully addressed by unbiased methodologies. Here we investigate the 3D PAM employing state-of-the-art finite-temperature auxiliary-field quantum Monte Carlo simulations. Furthermore, we present the behavior of the kinetic energy, the entropy, the specific heat, and the double occupancy as functions of the temperature and the hybridization strength. From these quantities, and by the analysis of the spin-spin correlation functions, we investigate the occurrence of magnetic phase transitions at finite temperatures, and we determine the phase diagram of the model, including the behavior of the Néel temperature as a function of the external parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

How correlations change the magnetic structure factor of the kagome Hubbard model

The kagome Hubbard model (KHM) is a paradigmatic example of a frustrated two-dimensional model. While its strongly correlated regime, described by a Heisenberg model, is of topical interest due to its enigmatic prospective spin-liquid ground state, the weakly and moderately correlated regimes remain largely unexplored. Motivated by the rapidly growing number of metallic kagome materials (e.g., Mn 3 Sn, Fe 3 Sn 2 , FeSn, Co 3 Sn 2 S 2 , Gd 3 Ru 4 Al 12 , and AV 3 Sb 5 with A = K, Rb, Cs), here we study the respective regimes of the KHM by means of three complementary numerical methods: the dynamical mean-field theory, the dynamical vertex approximation, and determinant quantum Monte Carlo. In contrast to the archetypal square lattice, we find no tendencies toward magnetic ordering, as magnetic correlations remain short-range. Nevertheless, the magnetic correlations undergo a remarkable crossover as the system approaches the metal-to-insulator transition. The Mott transition itself does not affect the magnetic correlations. Our equal-time and dynamical structure factors can be used as a reference for inelastic neutron scattering experiments on the growing family of metallic kagome materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Data Files for P. Mai et al., "Fluctuating charge-density-wave correlations in the three-band Hubbard model" (2024)

The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model’s fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.

cuprates↗
DOE OSTI2024

Ab initio electroweak corrections to superallowed β decays and their impact on V ud

Radiative corrections are essential for an accurate determination of V ud from superallowed β decays. In view of recent progress in the single-nucleon sector, the uncertainty is dominated by the theoretical description of nucleus-dependent effects, limiting the precision that can currently be achieved for V ud . In this work, we provide a detailed account of the electroweak corrections to superallowed β decays in effective field theory (EFT), including the power counting, potential and ultrasoft contributions, and factorization in the decay rate. Here, we present a first numerical evaluation of the dominant corrections in light nuclei based on quantum Monte Carlo methods, confirming the expectations from the EFT power counting. Finally, we discuss strategies how to extract from data the low-energy constants that parametrize short-distance contributions and whose values are not predicted by the EFT. Combined with advances in ab initio nuclear-structure calculations, this EFT framework allows one to systematically address the dominant uncertainty in V ud , as illustrated in detail for the 14 O → 14 N transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Quantitative assessment of the universal thermopower in the Hubbard model

As primarily an electronic observable, the room-temperature thermopower S in cuprates provides possibilities for a quantitative assessment of the Hubbard model. Using determinant quantum Monte Carlo, we demonstrate agreement between Hubbard model calculations and experimentally measured room-temperature S across multiple cuprate families, both qualitatively in terms of the doping dependence and quantitatively in terms of magnitude. We observe an upturn in S with decreasing temperatures, which possesses a slope comparable to that observed experimentally in cuprates. From our calculations, the doping at which S changes sign occurs in close proximity to a vanishing temperature dependence of the chemical potential at fixed density. Our results emphasize the importance of interaction effects in the systematic assessment of the thermopower S in cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantum magnetism on small-world networks

While classical spin systems in random networks have been intensively studied, much less is known about quantum magnets in random graphs. In this work, we investigate interacting quantum spins on small-world networks, building on mean-field theory and extensive quantum Monte Carlo simulations. Starting from one-dimensional (1D) rings, we consider two situations: All-to-all interacting and long-range interactions randomly added. The effective infinite dimension of the lattice leads to a magnetic ordering at finite temperature Tc with mean-field criticality. Nevertheless, in contrast to the classical case, we find two distinct power-law behaviors for Tc versus the average strength of the extra couplings. This is controlled by a competition between a characteristic length scale of the random graph and the thermal correlation length of the underlying 1D system, thus challenging mean-field theories. We also investigate the fate of a gapped 1D spin chain against the small-world effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗