Estimation of properties of low-lying excited states of Hubbard models: A multiconfigurational symmetrized projector quantum Monte Carlo approach

Abstract

We present in detail the recently developed multi-configurational symmetrized projector quantum Monte Carlo (MSPQMC) method for excited states of the Hubbard model. We describe the implementation of the Monte Carlo method for a multi-configurational trial wavefunction. We give a detailed discussion of issues related to the symmetry of the projection procedure which validates our Monte Carlo procedure for excited states and leads naturally to the idea of symmetrized sampling for correlation functions, developed earlier in the context of ground state simulations. It also leads to three possible averaging schemes. We have analyzed the errors incurred in these various averaging procedures and discuss and detail the preferred averaging procedure for correlations that do not have the full symmetry of the Hamiltonian. We study the energies and correlation functions of the low-lying excited states of the half-filled Hubbard model in 1-D. We have used this technique to study the pair-binding energies of two holes in $4n$ and $4n+2$ systems, which compare well the Bethe ansatz data of Fye, Martins and Scalettar. We have also studied small clusters amenable to exact diagonalization studies in 2-D and have reproduced their energies and correlation functions by the MSPQMC method. We identify two ways in which a multiconfigurational trial wavefunction can lead to a negative sign problem. We observe that this effect is not severe in 1-D and tends to vanish with increasing system size. We also note that this does not enhance the severity of the sign problem in two dimensions.