Delocalized Davydov D1 Ansatz for the Holstein polaron

An efficient, yet very accurate trial wave function, constructed from projecting the well-known Davydov D1 Ansatz onto momentum eigenstates, is employed to study the ground state properties of the generalized Holstein Hamiltonian with simultaneous diagonal and off-diagonal coupling. Ground-state energies have been obtained with a precision matching that of the computationally much more demanding density-matrix renormalization group method. The delocalized D1 Ansatz lowers the ground-state energies at the Brillouin zone boundary significantly compared with the Toyozawa and Global-Local Ansätze in the weak coupling regime, while considerable improvement is demonstrated to have been achieved over the entire Brillouin zone in the strong coupling regime. Unique solutions are obtained with the delocalized D1 for different initial conditions when the transfer integral is 20 times the phonon frequency at the zone center, implying the absence of any self-trapping discontinuity. The scaled correlation variance is found to fit satisfactorily well with the predictions of the perturbation theories.