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

Abstract

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.