Auxiliary-field Quantum Monte Carlo Research Articles

The pseudogap regime in the unitary Fermi gas

We discuss the pseudogap regime in the cold atomic unitary Fermi gas, with a particular emphasis on the auxiliary-field quantum Monte Carlo (AFMC) approach. We discuss possible signatures of the pseudogap, review experimental results, and survey analytic and quantum Monte Carlo techniques before focusing on AFMC calculations in the canonical ensemble. For the latter method, we discuss results for the heat capacity, energy-staggering pairing gap, spin susceptibility, and compare to experiment and other theoretical methods.

Finite-temperature auxiliary-field quantum Monte Carlo: Self-consistent constraint and systematic approach to low temperatures

We describe an approach for many-body calculations with a finite-temperature, grand canonical ensemble formalism using auxiliary-field quantum Monte Carlo (AFQMC) with a self-consistent constraint to control the sign problem. The usual AFQMC formalism of Blankenbecler, Scalapino, and Sugar suffers from the sign problem with most physical Hamiltonians, as is well known. Building on earlier ideas to constrain the paths in auxiliary-field space [Phys. Rev. Lett. \textbf{83}, 2777 (1999)] and incorporating recent developments in zero-temperature, canonical-ensemble methods, we discuss how a self-consistent constraint can be introduced in the finite-temperature, grand-canonical-ensemble framework. This together with several other algorithmic improvements discussed here leads to a more accurate, more efficient, and numerically more stable approach for finite-temperature calculations. We carry out a systematic benchmark study in the two-dimensional repulsive Hubbard model at $1/8$ doping. Temperatures as low as $T=1/80$ (in units of hopping) are reached. The finite-temperature method is exact at very high temperatures, and approaches the result of the zero-temperature constrained-path AFQMC as temperature is lowered. The benchmark shows that systematically accurate results are obtained for thermodynamic properties.

Overcoming the Memory Bottleneck in Auxiliary Field Quantum Monte Carlo Simulations with Interpolative Separable Density Fitting.

We investigate the use of interpolative separable density fitting (ISDF) as a means to reduce the memory bottleneck in auxiliary field quantum Monte Carlo (AFQMC) simulations of real materials in Gaussian basis sets. We find that ISDF can reduce the memory scaling of AFQMC simulations from [Formula: see text] to [Formula: see text]. We test these developments by computing the structural properties of carbon in the diamond phase, comparing to results from existing computational methods and experiment.

Auxiliary-field quantum Monte Carlo calculations of the structural properties of nickel oxide.

Auxiliary-field quantum Monte Carlo (AFQMC) has repeatedly demonstrated itself as one of the most accurate quantum many-body methods, capable of simulating both real and model systems. In this article, we investigate the application of AFQMC to realistic strongly correlated materials in periodic Gaussian basis sets. Using nickel oxide (NiO) as an example, we investigate the importance of finite size effects and basis set errors on the structural properties of the correlated solid. We provide benchmark calculations for NiO and compare our results to both experimental measurements and existing theoretical methods.

Ab initio computations of molecular systems by the auxiliary‐field quantum Monte Carlo method

The auxiliary‐field quantum Monte Carlo (AFQMC) method provides a computational framework for solving the time‐independent Schrödinger equation in atoms, molecules, solids, and a variety of model systems. AFQMC has recently witnessed remarkable growth, especially as a tool for electronic structure computations in real materials. The method has demonstrated excellent accuracy across a variety of correlated electron systems. Taking the form of stochastic evolution in a manifold of nonorthogonal Slater determinants, the method resembles an ensemble of density‐functional theory (DFT) calculations in the presence of fluctuating external potentials. Its computational cost scales as a low‐power of system size, similar to the corresponding independent‐electron calculations. Highly efficient and intrinsically parallel, AFQMC is able to take full advantage of contemporary high‐performance computing platforms and numerical libraries. In this review, we provide a self‐contained introduction to the exact and constrained variants of AFQMC, with emphasis on its applications to the electronic structure of molecular systems. Representative results are presented, and theoretical foundations and implementation details of the method are discussed.This article is categorized under:Electronic Structure Theory > Ab Initio Electronic Structure MethodsStructure and Mechanism > Computational Materials ScienceComputer and Information Science > Computer Algorithms and Programming

Communication: Calculation of interatomic forces and optimization of molecular geometry with auxiliary-field quantum Monte Carlo.

We propose an algorithm for accurate, systematic, and scalable computation of interatomic forces within the auxiliary-field quantum Monte Carlo (AFQMC) method. The algorithm relies on the Hellmann-Feynman theorem and incorporates Pulay corrections in the presence of atomic orbital basis sets. We benchmark the method for small molecules by comparing the computed forces with the derivatives of the AFQMC potential energy surface and by direct comparison with other quantum chemistry methods. We then perform geometry optimizations using the steepest descent algorithm in larger molecules. With realistic basis sets, we obtain equilibrium geometries in agreement, within statistical error bars, with experimental values. The increase in computational cost for computing forces in this approach is only a small prefactor over that of calculating the total energy. This paves the way for a general and efficient approach for geometry optimization and molecular dynamics within AFQMC.

Numerical results on the short-range spin correlation functions in the ground state of the two-dimensional Hubbard model

Optical lattice experiments with ultracold fermion atoms and quantum gas microscopy have recently realized direct measurements of magnetic correlations at the site-resolved level. We calculate the short-range spin-correlation functions in the ground state of the two-dimensional repulsive Hubbard model with the auxiliary-field quantum Monte Carlo (AFQMC) method. The results are numerically exact at half filling where the fermion sign problem is absent. Away from half filling, we employ the constrained path AFQMC approach to eliminate the exponential computational scaling from the sign problem. The constraint employs unrestricted Hartree-Fock trial wave functions with an effective interaction strength $U$, which is optimized self-consistently within AFQMC. Large supercells are studied, with twist averaged boundary conditions as needed, to reach the thermodynamic limit. We find that the nearest-neighbor spin correlation always increases with the interaction strength $U$, contrary to the finite-temperature behavior where a maximum is reached at a finite $U$ value. We also observe a change of sign in the next-nearest-neighbor spin correlation with increasing density, which is a consequence of the buildup of the long-range antiferromagnetic correlation. We expect the results presented in this paper to serve as a benchmark as lower temperatures are reached in ultracold atom experiments.

Response Functions for the Two-Dimensional Ultracold Fermi Gas: Dynamical BCS Theory and Beyond

Response functions are central objects in physics. They provide crucial information about the behavior of physical systems, and they can be directly compared with scattering experiments involving particles like neutrons, or photons. Calculations of such functions starting from the many-body Hamiltonian of a physical system are challenging, and extremely valuable. In this paper we focus on the two-dimensional (2D) ultracold Fermi atomic gas which has been realized experimentally. We present an application of the dynamical BCS theory to obtain response functions for different regimes of interaction strengths in the 2D gas with zero-range attractive interaction. We also discuss auxiliary-field quantum Monte Carlo (AFQMC) methods for the calculation of imaginary-time correlations in these dilute Fermi gas systems. Illustrative results are given and comparisons are made between AFQMC and dynamical BCS theory results to assess the accuracy of the latter.

Chemical Transformations Approaching Chemical Accuracy via Correlated Sampling in Auxiliary-Field Quantum Monte Carlo.

The exact and phaseless variants of auxiliary-field quantum Monte Carlo (AFQMC) have been shown to be capable of producing accurate ground-state energies for a wide variety of systems including those which exhibit substantial electron correlation effects. The computational cost of performing these calculations has to date been relatively high, impeding many important applications of these approaches. Here we present a correlated sampling methodology for AFQMC which relies on error cancellation to dramatically accelerate the calculation of energy differences of relevance to chemical transformations. In particular, we show that our correlated sampling-based AFQMC approach is capable of calculating redox properties, deprotonation free energies, and hydrogen abstraction energies in an efficient manner without sacrificing accuracy. We validate the computational protocol by calculating the ionization potentials and electron affinities of the atoms contained in the G2 test set and then proceed to utilize a composite method, which treats fixed-geometry processes with correlated sampling-based AFQMC and relaxation energies via MP2, to compute the ionization potential, deprotonation free energy, and the O-H bond disocciation energy of methanol, all to within chemical accuracy. We show that the efficiency of correlated sampling relative to uncorrelated calculations increases with system and basis set size and that correlated sampling greatly reduces the required number of random walkers to achieve a target statistical error. This translates to CPU-time speed-up factors of 55, 25, and 24 for the ionization potential of the K atom, the deprotonation of methanol, and hydrogen abstraction from the O-H bond of methanol, respectively. We conclude with a discussion of further efficiency improvements that may open the door to the accurate description of chemical processes in complex systems.

Cluster size convergence of the density matrix embedding theory and its dynamical cluster formulation: A study with an auxiliary-field quantum Monte Carlo solver

We investigate the cluster size convergence of the energy and observables using two forms of density matrix embedding theory (DMET): the original cluster form (CDMET) and a new formulation motivated by the dynamical cluster approximation (DCA-DMET). Both methods are applied to the half-filled one- and two-dimensional Hubbard models using a sign-problem free auxiliary-field quantum Monte Carlo (AFQMC) impurity solver, which allows for the treatment of large impurity clusters of up to 100 sites. While CDMET is more accurate at smaller impurity cluster sizes, DCA- DMET exhibits faster asymptotic convergence towards the thermodynamic limit (TDL). We use our two formulations to produce new accurate estimates for the energy and local moment of the two-dimensional Hubbard model for U/t = 2, 4, 6. These results compare favourably with the best data available in literature, and help resolve earlier uncertainties in the moment for U/t = 2.

Rashba Spin-Orbit Coupling, Strong Interactions, and the BCS-BEC Crossover in the Ground State of the Two-Dimensional Fermi Gas.

The recent experimental realization of spin-orbit coupled Fermi gases provides a unique opportunity to study the interplay between strong interaction and spin-orbit coupling (SOC) in a tunable, disorder-free system. We present here precision abinitio numerical results on the two-dimensional, unpolarized, uniform Fermi gas with attractive interactions and Rashba SOC. Using the auxiliary-field quantum MonteCarlo method and incorporating recent algorithmic advances, we carry out exact calculations on sufficiently large system sizes to provide accurate results systematically as a function of experimental parameters. We obtain the equation of state, the momentum distributions, the pseudospin correlations, and the pair wave functions. Our results help illuminate the rich pairing structure induced by SOC, and provide benchmarks for theory and guidance to future experimental efforts.

Net-baryon number fluctuations across the chiral phase transition at finite density in strong-coupling lattice QCD

We investigate the net-baryon number fluctuations across the chiral phase transition at finite density in the strong coupling and chiral limit. Mesonic field fluctuations are taken into account by using the auxiliary field Monte-Carlo method. We find that the higher-order cumulant ratios, $S\sigma$ and $\kappa\sigma^2$, show oscillatory behavior around the phase boundary at $\mu/T\gtrsim 0.2$, and there exists the region where the higher-order cumulant ratios are negative. The negative region of $\kappa\sigma^2$ is found to shrink with increasing lattice size. This behavior agrees with the expectations from the scaling analysis.

Solving the fermion sign problem in quantum Monte Carlo simulations by Majorana representation

In this paper, we discover a new quantum Monte Carlo (QMC) method to solve the fermion sign problem in interacting fermion models by employing Majorana representation of complex fermions. We call it "Majorana QMC" (MQMC). Especially, MQMC is fermion sign free in simulating a class of spinless fermion models on bipartite lattices at half filling and with arbitrary range of (unfrustrated) interactions. To the best of our knowledge, MQMC is the first auxiliary field QMC method to solve fermion sign problem in spinless (more generally, odd number of species) fermion models. MQMC simulations can be performed efficiently both at finite and zero temperatures. We believe that MQMC could pave a new avenue to solve fermion sign problem in more generic fermionic models.

Quantum Monte Carlo Calculations in Solids with Downfolded Hamiltonians.

We present a combination of a downfolding many-body approach with auxiliary-field quantum Monte Carlo (AFQMC) calculations for extended systems. Many-body calculations operate on a simpler Hamiltonian which retains material-specific properties. The Hamiltonian is systematically improvable and allows one to dial, in principle, between the simplest model and the original Hamiltonian. As a by-product, pseudopotential errors are essentially eliminated using frozen orbitals constructed adaptively from the solid environment. The computational cost of the many-body calculation is dramatically reduced without sacrificing accuracy. Excellent accuracy is achieved for a range of solids, including semiconductors, ionic insulators, and metals. We apply the method to calculate the equation of state of cubic BN under ultrahigh pressure, and determine the spin gap in NiO, a challenging prototypical material with strong electron correlation effects.

An auxiliary-field quantum Monte Carlo study of the chromium dimer.

The chromium dimer (Cr2) presents an outstanding challenge for many-body electronic structure methods. Its complicated nature of binding, with a formal sextuple bond and an unusual potential energy curve (PEC), is emblematic of the competing tendencies and delicate balance found in many strongly correlated materials. We present an accurate calculation of the PEC and ground state properties of Cr2, using the auxiliary-field quantum Monte Carlo (AFQMC) method. Unconstrained, exact AFQMC calculations are first carried out for a medium-sized but realistic basis set. Elimination of the remaining finite-basis errors and extrapolation to the complete basis set limit are then achieved with a combination of phaseless and exact AFQMC calculations. Final results for the PEC and spectroscopic constants are in excellent agreement with experiment.

Auxiliary field Monte Carlo simulation of strong-coupling lattice QCD for QCD phase diagram

We study the QCD phase diagram in the strong coupling limit with fluctuation effects by using the auxiliary field Monte-Carlo method. We apply the chiral angle fixing technique in order to obtain finite chiral condensate in the chiral limit in finite volume. The behavior of order parameters suggests that chiral phase transition is the second order or crossover at low chemical potential and the first order at high chemical potential. Compared with the mean field results, the hadronic phase is suppressed at low chemical potential, and is extended at high chemical potential as already suggested in the monomer-dimer-polymer simulations. We find that the sign problem originating from the bosonization procedure is weakened by the phase cancellation mechanism; a complex phase from one site tends to be canceled by the nearest neighbor site phase as long as low momentum auxiliary field contributions dominate.

Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method

Quantum Monte Carlo methods are powerful techniques for studying strongly interacting Fermi systems. However, implementing these methods on computers with finite-precision arithmetic requires careful attention to numerical stability. In the auxiliary-field Monte Carlo (AFMC) method, low-temperature or large-model-space calculations require numerically stabilized matrix multiplication. When adapting methods used in the grand-canonical ensemble to the canonical ensemble of fixed particle number, the numerical stabilization increases the number of required floating-point operations for computing observables by a factor of the size of the single-particle model space, and thus can greatly limit the systems that can be studied. We describe an improved method for stabilizing canonical-ensemble calculations in AFMC that exhibits better scaling, and present numerical tests that demonstrate the accuracy and improved performance of the method.

CPMC-Lab: A Matlab package for Constrained Path Monte Carlo calculations

We describe CPMC-Lab, a Matlab program for the constrained-path and phaseless auxiliary-field Monte Carlo methods. These methods have allowed applications ranging from the study of strongly correlated models, such as the Hubbard model, to ab initio calculations in molecules and solids. The present package implements the full ground-state constrained-path Monte Carlo (CPMC) method in Matlab with a graphical interface, using the Hubbard model as an example. The package can perform calculations in finite supercells in any dimensions, under periodic or twist boundary conditions. Importance sampling and all other algorithmic details of a total energy calculation are included and illustrated. This open-source tool allows users to experiment with various model and run parameters and visualize the results. It provides a direct and interactive environment to learn the method and study the code with minimal overhead for setup. Furthermore, the package can be easily generalized for auxiliary-field quantum Monte Carlo (AFQMC) calculations in many other models for correlated electron systems, and can serve as a template for developing a production code for AFQMC total energy calculations in real materials. Several illustrative studies are carried out in one- and two-dimensional lattices on total energy, kinetic energy, potential energy, and charge- and spin-gaps. Program summaryProgram title: CPMC-LabCatalogue identifier: AEUD_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEUD_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 2850No. of bytes in distributed program, including test data, etc.: 24838Distribution format: tar.gzProgramming language:Matlab.Computer: The non-interactive scripts can be executed on any computer capable of running Matlab with all Matlab versions. The GUI requires Matlab R2010b (version 7.11) and above. Operating system: Windows, Mac OS X, Linux.RAM: Variable.Classification: 7.3.External routines:MatlabNature of problem:Obtaining ground state energy of a repulsive Hubbard model in a supercell in any number of dimensions.Solution method:In the Constrained Path Monte Carlo (CPMC) method, the ground state of a many-fermion system is projected from an initial trial wave function by a branching random walk in an overcomplete basis of Slater determinants. Constraining the determinants according to a trial wave function |ΨT〉 removes the exponential decay of the signal-to-noise ratio characteristic of the sign problem. The method is exact if |ΨT〉 is exact.Unusual features:Direct and interactive environment with a Graphical User Interface for beginners to learn and study the Constrained Path Monte Carlo method with minimal overhead for setup.Running time:The sample provided takes a few seconds to run, the batch sample a few minutes.

Imaginary time correlations and the phaseless auxiliary field quantum Monte Carlo

The phaseless Auxiliary Field Quantum Monte Carlo (AFQMC) method provides a well established approximation scheme for accurate calculations of ground state energies of many-fermions systems. Here we address the possibility of calculating imaginary time correlation functions with the phaseless AFQMC. We give a detailed description of the technique and test the quality of the results for static properties and imaginary time correlation functions against exact values for small systems.

Frozen-Orbital and Downfolding Calculations with Auxiliary-Field Quantum Monte Carlo.

We describe the implementation of the frozen-orbital and downfolding approximations in the auxiliary-field quantum Monte Carlo (AFQMC) method. These approaches can provide significant computational savings, compared to fully correlating all of the electrons. While the many-body wave function is never explicit in AFQMC, its random walkers are Slater determinants, whose orbitals may be expressed in terms of any one-particle orbital basis. It is therefore straightforward to partition the full N-particle Hilbert space into active and inactive parts to implement the frozen-orbital method. In the frozen-core approximation, for example, the core electrons can be eliminated in the correlated part of the calculations, greatly increasing the computational efficiency, especially for heavy atoms. Scalar relativistic effects are easily included using the Douglas-Kroll-Hess theory. Using this method, we obtain a way to effectively eliminate the error due to single-projector, norm-conserving pseudopotentials in AFQMC. We also illustrate a generalization of the frozen-orbital approach that downfolds high-energy basis states to a physically relevant low-energy sector, which allows a systematic approach to produce realistic model Hamiltonians to further increase efficiency for extended systems.