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

1/4 is the new 1/2 when topology is intertwined with Mottness

In non-interacting systems, bands from non-trivial topology emerge strictly at half-filling and exhibit either the quantum anomalous Hall or spin Hall effects. Here we show using determinantal quantum Monte Carlo and an exactly solvable strongly interacting model that these topological states now shift to quarter filling. A topological Mott insulator is the underlying cause. The peak in the spin susceptibility is consistent with a possible ferromagnetic state at T = 0. The onset of such magnetism would convert the quantum spin Hall to a quantum anomalous Hall effect. While such a symmetry-broken phase typically is accompanied by a gap, we find that the interaction strength must exceed a critical value for this to occur. Hence, we predict that topology can obtain in a gapless phase but only in the presence of interactions in dispersive bands. These results explain the recent quarter-filled quantum anomalous Hall effects seen in moiré systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Broken Luttinger theorem in the two-dimensional Fermi-Hubbard model

One of the fundamental questions about high-temperature cuprate superconductors is the size of the Fermi surface underlying the superconducting state. By analyzing the single-particle spectral function for the Fermi-Hubbard model as a function of repulsion U and chemical potential μ, we find that the Fermi surface in the normal state undergoes a transition from a large Fermi surface matching the Luttinger volume as expected in a Fermi liquid, to a Fermi surface that encloses fewer electrons that we dub the “Luttinger breaking” phase, as the Mott insulator is approached. This transition into a non-Fermi-liquid phase that violates the Luttinger count occurs at a critical density in the absence of any other broken symmetry. We obtain the Fermi-surface contour from the spectral weight A k (ω = 0) and from an analysis of the singularities of the Green's function Re G k (E = 0), calculated using determinantal quantum Monte Carlo and analytic continuation methods. We discuss our numerical results in connection with experiments on Hall measurements, scanning tunneling spectroscopy, and angle-resolved photoemission spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Large-𝑁 SU(4) Schwinger boson theory for coupled-dimer antiferromagnets

Here, we develop a systematic large-𝑁 expansion based on the Schwinger boson representation of SU(4) coherent states of dimers for the paradigmatic spin-1/2 bilayer square lattice Heisenberg antiferromagnet. This system exhibits a quantum phase transition between a quantum paramagnetic state and a Néel order state, driven by the coupling constant 𝑔 = 𝐽′/𝐽, which is defined as the ratio between the interdimer 𝐽′ and intradimer 𝐽 exchange interactions. We demonstrate that this approach accurately describes static and dynamic properties on both sides of the quantum phase transition. The critical coupling constant 𝑔 𝑐 ≈ 0.42 and the dynamic spin structure factor reproduce quantum Monte Carlo results with high precision. Notably, the 1/𝑁 corrections reveal the longitudinal mode of the magnetically ordered phase along with the overdamping caused by its decay into the two-magnon continuum. The present large-𝑁 SU(4) Schwinger boson theory can be extended to more general cases of quantum paramagnets that undergo a quantum phase transition into magnetically ordered states.

Schwinger boson method↗
DOE OSTI2020

Dirty bosons on the Cayley tree: Bose-Einstein condensation versus ergodicity breaking

Building on large-scale quantum Monte Carlo simulations, we investigate the zero-temperature phase diagram of hard-core bosons in a random potential on site-centered Cayley trees with branching number K=2. In order to follow how the Bose-Einstein condensate (BEC) is affected by the disorder, we focus on both the zero-momentum density, probing the quantum coherence, and the one-body density matrix (1BDM) whose largest eigenvalue monitors the off-diagonal long-range order. We further study its associated eigenstate which brings useful information about the real-space properties of this leading eigenmode. Upon increasing randomness, we find that the system undergoes a quantum phase transition at finite disorder strength between a long-range ordered BEC state, fully ergodic at large scale, and a new disordered Bose glass regime showing conventional localization for the coherence fraction while the 1BDM displays a nontrivial algebraic vanishing BEC density together with a nonergodic occupation in real space. These peculiar properties can be analytically captured by a simple phenomenological description on the Cayley tree which provides a physical picture of the Bose glass regime.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fluctuating intertwined stripes in the strange metal regime of the Hubbard model

Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here, in this paper, we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

High Accuracy Transition Metal Effective Cores for the Many-Body Diffusion Monte Carlo Method

Practical applications of the real-space diffusion Monte Carlo (DMC) method require the removal of core electrons, where currently localization approximations of semilocal potentials are generally used in the projector. Accurate calculations of complex solids and large molecules demand minimizing the impact of approximated atomic cores. Prior works have shown that the errors from such approximations can be sizable in both finite and periodic systems. In this work, we show that a class of differential pseudopotentials, known as pseudo-Hamiltonians, can be constructed for the 3d transition metal atoms, entirely removing the need for any localization scheme in the DMC projector. As a proof of principle, we demonstrate the approach for the case of Co. In order to minimize errors in the pseudo-Hamiltonian at the many-body level, we generalize the recently proposed correlation-consistent pseudopotential generation scheme to successively close semilocal representations of the differential potentials. Our generation scheme successfully produces potentials tailored specifically for real space projector quantum Monte Carlo methods with low error at the many-body level, i.e., with many-body scattering properties very close to relativistic all-electron results. In particular, we show that the agreement with respect to atomic and molecular quantities reach chemical accuracy in many cases-on par with the most accurate semilocal pseudopotentials available. Further, our pseudo-Hamiltonian generation scheme utilizes standard quantum chemistry codes designed only to work with semilocal pseudopotentials, enabling straightforward generation of pseudo-Hamiltonians for additional elements in future works.

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

Quenching of the quantum p -strength in light exotic nuclei

One-nucleon spectroscopic overlaps, their strength or spectroscopic factors (SFs), and nucleon removal cross sections for light nuclei in the mass range A≤12 were evaluated using the fully correlated Quantum Monte Carlo (QMC) wave functions (WFs). Harder (AV18+UX) and softer (NV2+3) bare interactions were used providing consistent results. The SFs were also taken from simple Shell Model (SM) calculations. This structure information was incorporated in the standard three-body Faddeev/Alt-Grassberger-Sandhas reaction formalism to evaluate the (p,pN) cross sections. The results further our understanding of the quenching of the quantum p-shell strength obtained from structure and reactions. We have found the ratios of the total QMC sums of SFs to the SM ones to be uniform and ∼ 3/4. The corresponding ratios of the sums Below Particle Threshold (BPT) of SFs, as well as of total cross sections, deviate strikingly from the uniform trend in some special cases. We find these ratios to be close to unit for 9 Li → 8 Li + n and 9 C → 8 B + p. In contrast, they are strongly reduced for 11 C → 10 C + n and 11 B → 10 Be + p mirror transitions, resulting from the quenching of the strength for low-lying 2 + final state transitions. The QMC theoretical cross section BPT for 11 C(p,pn) is about two times smaller than the experimental data.

(p,pN) reactions↗
DOE OSTI2025

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Composite Qdrift-product formulas for quantum and classical simulations in real and imaginary time

Recent study has shown that it can be advantageous to implement a composite channel that partitions the Hamiltonian H for a given simulation problem into subsets A and B such that H = A + B , where the terms in A are simulated with a Trotter-Suzuki channel and the B terms are randomly sampled via the Qdrift algorithm. Here we extend Qdrift and composite product formulas to imaginary time, formulating candidate classical algorithms for quantum Monte Carlo calculations. We upper bound the induced Schatten- 1 → 1 norm on both imaginary-time Qdrift and composite channels. Another recent result demonstrated that simulations of lattice Hamiltonians containing geometrically local interactions can be improved using a Lieb-Robinson argument to decompose H into subsets that contain only terms supported on that subset of the lattice. Here, we provide a quantum algorithm by unifying this result with the composite approach into “local composite channels” and we upper bound the diamond distance. We provide exact numerical simulations of algorithmic cost by counting the number of gates of the form e − i H j t and e − H j β to meet a certain error tolerance ε . In doing so, we optimize the partitioning into sets A and B using gradient boosted tree models from machine learning. These numerical studies are important given that product formulas have been historically known to outperform analytic upper bounds. We show constant factor advantages for a variety of interesting Hamiltonians, the maximum of which is a ≈ 20 -fold speedup that occurs in the simulation of Jellium. Published by the American Physical Society 2024

Pocrnic, Matthew (ORCID:0000000203089376)↗
DOE OSTI2024

Particle-Hole Asymmetric Ferromagnetism and Spin Textures in the Triangular Hubbard-Hofstadter Model

In a lattice model subject to a perpendicular magnetic field, when the lattice constant is comparable to the magnetic length, one enters the “Hofstadter regime,” where continuum Landau levels become fractal magnetic Bloch bands. Strong mixing between bands alters the nature of the resulting quantum phases compared to the continuum limit; lattice potential, magnetic field, and Coulomb interaction must be treated on equal footing. Using determinant quantum Monte Carlo and density matrix renormalization group techniques, we study this regime numerically in the context of the Hubbard-Hofstadter model on a triangular lattice. In the field-filling phase diagram, we find a broad wedge-shaped region of ferromagnetic ground states for filling factor ν ≤ 1 , bounded below by filling factor ν = 1 and bounded above by half filling the lowest Hofstadter subband. We observe signatures of SU(2) quantum Hall ferromagnetism at filling factors ν = 1 and ν = 3 . The phases near ν = 1 are particle-hole asymmetric, and we observe a rapid decrease in ground-state spin polarization consistent with the formation of skyrmions only on the electron doped side. At large fields, above the ferromagnetic wedge, we observe a low-spin metallic region with spin correlations peaked at small momenta. We argue that the phenomenology of this region likely results from exchange interaction mixing fractal Hofstadter subbands. The phase diagram derived beyond the continuum limit points to a rich landscape to explore interaction effects in magnetic Bloch bands. Published by the American Physical Society 2024

Ding, Jixun K.↗
DOE OSTI2020

Improved Fast Randomized Iteration Approach to Full Configuration Interaction

In this work, we present three modifications to our recently introduced fast randomized iteration method for full configuration interaction (FCI-FRI) and investigate their effects on the method’s performance for Ne, H 2 O, and N 2 . The initiator approximation, originally developed for full configuration interaction quantum Monte Carlo, significantly reduces statistical error in FCI-FRI when few samples are used in compression operations, enabling its application to larger chemical systems. The semistochastic extension, which involves exactly preserving a fixed subset of elements in each compression, improves statistical efficiency in some cases but reduces it in others. We also developed a new approach to sampling excitations that yields consistent improvements in statistical efficiency and reductions in computational cost. We discuss possible strategies based on our findings for improving the performance of stochastic quantum chemistry methods more generally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2025

Quantum mechanical dataset of 836k neutral closed-shell molecules with up to 5 heavy atoms from C, N, O, F, Si, P, S, Cl, Br

Abstract We introduce the Vector-QM24 (VQM24) dataset comprehensively covering all possible neutral closed-shell small organic and inorganic molecules with up to five heavy (p-block) atoms: C, N, O, F, Si, P, S, Cl, Br. All valid stoichiometries, Lewis-rule-consistent graphs, and stable conformers (identified via GFN2-xTB) were enumerated combinatorially, yielding 577k conformational isomers spanning 258k constitutional isomers and 5,599 unique stoichiometries. DFT (ωB97X-D3/cc-pVDZ) optimizations were performed for all, and diffusion quantum Monte Carlo (DMC@PBE0(ccECP/cc-pVQZ)) energies are provided for 10,793 lowest-energy conformers with up to 4 heavy atoms. VQM24 includes structures, vibrational modes, rotational constants, thermodynamic properties (Gibbs free energies, enthalpies, ZPVEs, entropies, heat capacities), and electronic properties such as atomization, electron interaction, exchange-correlation, dispersion energies, multipole moments (dipole to hexadecapole), alchemical potentials, Mulliken charges, and wavefunctions. Machine learning models of atomization energies on this dataset reveal significantly higher complexity than QM9, with none achieving chemical accuracy. VQM24 offers a rigorous, high-fidelity benchmark for evaluating quantum machine learning models.

Science & Technology - Other Topics↗
DOE OSTI2021

Monolayer CrCl 3 as an Ideal Test Bed for the Universality Classes of 2D Magnetism

The monolayer halides CrX 3 (X=Cl, Br, I) attract significant attention for realizing 2D magnets with genuine long-range order (LRO), challenging the Mermin-Wagner theorem. In this work, we show that monolayer CrCl 3 has the unique benefit of exhibiting tunable magnetic anisotropy upon applying a compressive strain. This opens the possibility to use CrCl 3 for producing and studying both ferromagnetic and antiferromagnetic 2D Ising-type LRO as well as the Berezinskii-Kosterlitz-Thouless (BKT) regime of 2D magnetism with quasi-LRO. Using state-of-the-art density functional theory, we explain how realistic compressive strain could be used to tune the monolayer's magnetic properties so that it could exhibit any of these phases. Building on large-scale quantum Monte Carlo simulations, we compute the phase diagram of strained CrCl 3 , as well as the magnon spectrum with spin-wave theory. Our results highlight the eminent suitability of monolayer CrCl 3 to achieve very high BKT transition temperatures, around 50 K, due to their singular dependence on the weak easy-plane anisotropy of the material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Emergent Glassy Behavior in a Kagome Rydberg Atom Array

Here, we present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal emergent glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids—as well as the crossover towards a trivial paramagnetic phase—are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

DC Hall coefficient of the strongly correlated Hubbard model

The Hall coefficient is related to the effective carrier density and Fermi surface topology in non-interacting and weakly interacting systems. In strongly correlated systems, the relation between the Hall coefficient and single-particle properties is less clear. Clarifying this relation would give insight into the nature of transport in strongly correlated materials that lack well-formed quasiparticles. In this work, we investigate the DC Hall coefficient of the Hubbard model using determinant quantum Monte Carlo in conjunction with a recently developed expansion of magneto-transport coefficients in terms of thermodynamic susceptibilities. At leading order in the expansion, we observe a change of sign in the Hall coefficient as a function of temperature and interaction strength, which we relate to a change in the topology of the apparent Fermi surface. We also combine our Hall coefficient results with optical conductivity values to evaluate the Hall angle, as well as effective mobility and effective mass based on Drude theory of metals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-level coupled-cluster energetics by merging moment expansions with selected configuration interaction

Inspired by our earlier semi-stochastic work aimed at converging high-level coupled-cluster (CC) energetics, we propose a novel form of the CC(P; Q) theory in which the stochastic Quantum Monte Carlo propagations, used to identify dominant higher-than-doubly excited determinants, are replaced by the selected configuration interaction (CI) approach using the perturbative selection made iteratively (CIPSI) algorithm. The advantages of the resulting CIPSI-driven CC(P; Q) methodology are illustrated by a few molecular examples, including the dissociation of F 2 and the automerization of cyclobutadiene, where we recover the electronic energies corresponding to the CC calculations with a full treatment of singles, doubles, and triples based on the information extracted from compact CI wave functions originating from relatively inexpensive Hamiltonian diagonalizations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗