Field Quantum Monte Carlo Method Research Articles

Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids.

We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and, thus, reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and the corresponding supercell calculations numerically verify our implementation. We calculate the equation of state for diamond and the correlation energies for a range of prototypical solid materials. A down-sampling technique along with natural orbitals accelerates the convergence with respect to the number of orbitals and crystal momentum points. We illustrate the competitiveness of our implementation in accuracy and computational cost for dense crystal momentum point meshes compared to a well-established quantum-chemistry approach, the coupled-cluster ansatz including singles, doubles, and perturbative triple particle-hole excitation operators.

Self-consistent optimization of the trial wave function within the constrained path auxiliary field quantum Monte Carlo method using mixed estimators

We propose a scheme to implement the self-consistent optimization of the trial wave function within the constrained path auxiliary field quantum Monte Carlo (CP-AFQMC) method in the framework of natural orbitals. In this scheme, a new trial wave function in the form of a Slater determinant is constructed from the CP-AFQMC results by diagonalizing the mixed estimator of the one-body reduced density matrix. We compare two ways (from real and mixed estimators in CP-AFQMC) to calculate the one-body reduced density matrix in the self-consistent process and study the ground state of a doped two-dimensional Hubbard model to test the accuracy of the two schemes. By comparing the local density, occupancy, and ground state energy we find the scheme in which a one-body reduced density matrix is calculated from the mixed estimator is computationally more efficient and provides more accurate results with less fluctuation. The local densities from the mixed estimator scheme agree well with the numerically exact values. This scheme provides a useful tool for studying strongly correlated electron systems.

Higher Order Auxiliary Field Quantum Monte Carlo Methods

The auxiliary field quantum Monte Carlo (AFQMC) method has been a workhorse in the field of strongly correlated electrons for a long time and has found its most recent implementation in the ALF package (alf.physik.uni-wuerzburg.de). The utilization of the Trotter decomposition to decouple the interaction from the non-interacting Hamiltonian makes this method inherently second order in terms of the imaginary time slice. We show that due to the use of the Hubbard-Stratonovich transformation (HST) a semigroup structure on the time evolution is imposed that necessitates the introduction of a new family of complex-hermitian splitting methods for the purpose of reaching higher order. We will give examples of these new methods and study their efficiency, as well as perform comparisons with other established second and higher order methods in the realm of the AFQMC method.

Metal-insulator and magnetic phase diagram of Ca2RuO4 from auxiliary field quantum Monte Carlo and dynamical mean field theory

Layered perovskite ruthenium oxides exhibit a striking series of metal-insulator and magnetic-nonmagnetic phase transitions easily tuned by temperature, pressure, epitaxy, and nonlinear drive. In this work, we combine results from two complementary state of the art many-body methods, Auxiliary Field Quantum Monte Carlo and Dynamical Mean Field Theory, to determine the low-temperature phase diagram of Ca$_2$RuO$_4$. Both methods predict a low temperature, pressure-driven metal-insulator transition accompanied by a ferromagnetic-antiferromagnetic transition. The properties of the ferromagnetic state vary non-monotonically with pressure and are dominated by the ruthenium $d_{xy}$ orbital, while the properties of the antiferromagnetic state are dominated by the $d_{xz}$ and $d_{yz}$ orbitals. Differences of detail in the predictions of the two methods are analyzed. This work is theoretically important as it presents the first application of the Auxiliary Field Quantum Monte Carlo method to an orbitally-degenerate system with both Mott and Hunds physics, and provides an important comparison of the Dynamical Mean Field and Auxiliary Field Quantum Monte Carlo methods.

Ab Initio Finite Temperature Auxiliary Field Quantum Monte Carlo.

We present an ab initio auxiliary field quantum Monte Carlo method for studying the electronic structure of molecules, solids, and model Hamiltonians at finite temperature. The algorithm marries the ab initio phaseless auxiliary field quantum Monte Carlo algorithm known to produce high accuracy ground state energies of molecules and solids with its finite temperature variant, long used by condensed matter physicists for studying model Hamiltonian phase diagrams, to yield a phaseless, ab initio finite temperature method. We demonstrate that the method produces internal energies within chemical accuracy of exact diagonalization results across a wide range of temperatures for H2O (STO-3G), C2 (STO-6G), the one-dimensional hydrogen chain (STO-6G), and the multiorbital Hubbard model. Our method effectively controls the phase problem through importance sampling, often even without invoking the phaseless approximation, down to temperatures at which the systems studied approach their ground states and may therefore be viewed as exact over wide temperature ranges. This technique embodies a versatile tool for studying the finite temperature phase diagrams of a plethora of systems whose properties cannot be captured by a Hubbard U term alone. Our results moreover illustrate that the severity of the phase problem for model Hamiltonians far exceeds that for many molecules at all of the temperatures studied.

Majorana Positivity and the Fermion Sign Problem of Quantum MonteCarlo Simulations.

The sign problem is a major obstacle in quantum MonteCarlo simulations for many-body fermion systems. We examine this problem with a new perspective based on the Majorana reflection positivity and Majorana Kramers positivity. Two sufficient conditions are proven for the absence of the fermion sign problem. Our proof provides a unified description for all the interacting lattice fermion models previously known to be free of the sign problem based on the auxiliary field quantum MonteCarlo method. It also allows us to identify a number of new sign-problem-free interacting fermion models including, but not limited to, lattice fermion models with repulsive interactions but without particle-hole symmetry, and interacting topological insulators with spin-flip terms.

Imaginary time density-density correlations for two-dimensional electron gases at high density.

We evaluate imaginary time density-density correlation functions for two-dimensional homogeneous electron gases of up to 42 particles in the continuum using the phaseless auxiliary field quantum Monte Carlo method. We use periodic boundary conditions and up to 300 plane waves as basis set elements. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products, and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems. We discuss the numerical stabilization techniques and the computational complexity of the methodology and we present the limitations related to the size of the systems on a quantitative basis. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of medium-sized homogeneous fermion systems.

Quantum Monte Carlo study of the superconducting Kosterlitz-Thouless transition of the attractive Hubbard model on a triangular lattice

We study the superconducting Kosterlitz-Thouless transition of the attractive Hubbard model on a two-dimensional triangular lattice using auxiliary field quantum Monte Carlo method for system sizes up to $12\times 12$ sites. Combining three methods to analyze the numerical data, we find, for the attractive interaction of $U=-4t$, that the transition temperature stays almost constant within the band filling range of $1.0 < n < 1.4$, while it is found to be much lower in the $n<1$ region.

Numerical Study for the Ground State of Multi-Orbital Hubbard Models

Ground state properties of multi-orbital Hubbard models are investigated by the auxiliary field quantum Monte Carlo method. A Monte Carlo technique generalized to the multi-orbital systems is introduced and examined in detail. The algorithm contains non-trivial cases where the negative sign problem does not exist. We investigate one-dimensional systems with doubly degenerate orbitals by this new technique. Properties of the Mott insulating state are quantitatively clarified as the strongly correlated insulator, where the charge gap amplitude is much larger than the spin gap. The insulator-metal transitions driven by the chemical potential shows a universality class with the correlation length exponent $\nu=1/2$, which is consistent with the scaling arguments. Increasing level split between two orbitals drives crossover from the Mott insulator with high spin state to the band insulator with low spin state, where the spin gap amplitude increases and becomes closer to the charge gap. Experimental relevance of our results especially to Haldane materials is discussed.