Spectral evolution with doping of an antiferromagnetic Mott state

Since the discovery of half-filled cuprate to be a Mott insulator, the excitation spectra above the chemical potential for the unoccupied states has attracted much research attention. There were many theoretical works using different numerical techniques to study this problem, but many have reached different conclusions. One of the reasons is the lack of very detailed high-resolution experimental results for the theories to be compared with. Recently, the scanning tunneling spectroscopy [P. Cai et al., Nat. Phys. 12, 1047 (2016); C. Ye et al., Nat. Commun. 4, 1365 (2013)] on lightly doped Mott insulator with an antiferromagnetic order found the presence of in-gap states with energy of order half an eV above the chemical potential. The measured spectral properties with doping are not quite consistent with earlier theoretical works. Although the experiment has disorder and localization effect, but for the energy scale we will study here, a model without disorder is sufficed to illustrate the underlying physics. We perform a diagonalization method on top of the variational Monte Carlo calculation to study the evolution of antiferromagnetic Mott state with doped hole concentration in the Hubbard model. Our results found in-gap states that behave similarly with ones reported by STS. These in-gap states acquire a substantial amount of dynamical spectral weight transferred from the upper Hubbard band. The in-gap states move toward chemical potential with increasing spectral weight as doping increases. Our result also provides information about the energy scale of these in-gap states in relation with the Coulomb coupling strength U.