Auxiliary Field Quantum Monte Carlo Simulations of Real Strongly Correlated Materials

Wirawan Purwanto, Shiwei Zhang, Henry Krakauer

Shuai Zhang, Fionn D. Malone, Miguel A. Morales

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.

Shiwei Zhang

Fionn D. Malone, Shuai Zhang, Miguel A. Morales

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.

S. Hochkeppel, Thomas C. Lang, C. Brünger +2 more

Hao Shi, Shiwei Zhang

The auxiliary-field quantum Monte Carlo (AFQMC) method is a general numerical method for correlated many-electron systems, which is being increasingly applied in lattice models, atoms, molecules, and solids. Here, we introduce the theory and algorithm of the method specialized for real materials and present several recent developments. We give a systematic exposition of the key steps of AFQMC, closely tracking the framework of a modern software library we are developing. The building of a Monte Carlo Hamiltonian, projecting to the ground state, sampling two-body operators, phaseless approximation, and measuring ground state properties are discussed in detail. An advanced implementation for multi-determinant trial wave functions is described, which dramatically speeds up the algorithm and reduces the memory cost. We propose a self-consistent constraint for real materials, and discuss two flavors for its realization, either by coupling the AFQMC calculation to an effective independent-electron calculation or via the natural orbitals of the computed one-body density matrix.