Projector quantum Monte Carlo with matrix product states
We marry tensor network states (TNS) and projector quantum Monte Carlo (PMC) to overcome the high computational scaling of TNS and the sign problem of PMC. Using TNS as trial wave functions provides a route to systematically improve the sign structure and to eliminate the bias in fixed-node and constrained-path PMC. As a specific example, we describe phaseless auxiliary-field quantum Monte Carlo with matrix product states (MPS-AFQMC). MPS-AFQMC improves significantly on the density-matrix renormalization group (DMRG) ground-state energy. For the J_1 − J_2 model on two-dimensional square lattices, we observe with MPS-AFQMC an order of magnitude reduction in the error for all couplings, compared to DMRG. The improvement is independent of walker bond dimension, and we therefore use bond dimension 1 for the walkers. The computational cost of MPS-AFQMC is then quadratic in the bond dimension of the trial wave function, which is lower than the cubic scaling of DMRG. The error due to the constrained-path bias is proportional to the variational error of the trial wave function. We show that for the J_1 − J_2 model on two-dimensional square lattices, a linear extrapolation of the MPS-AFQMC energy with the discarded weight from the DMRG calculation allows us to remove the constrained-path bias. Extensions to other tensor networks are briefly discussed.
© 2014 American Physical Society. Received 12 March 2014; revised manuscript received 15 April 2014; published 3 July 2014. S.W. received a Ph.D. fellowship and B.V. a postdoctoral fellowship from the Research Foundation Flanders (FWO Vlaanderen). The work was supported by the US National Science Foundation through Grants No. OCI-1265278 and No. CHE-1265277.
Accepted Version - 1403.3125.pdf
Published - PhysRevB.90.045104.pdf