Superconductivity and inhomogeneous charge-ordered state in the two-Dimensional Hubbard model -- Off-Diagonal Wave Function Monte Carlo Studies of Hubbard Model IV

We investigate the ground-state properties of the two-dimensional Hubbard model, based on the off-diagonal wave function variational Monte Carlo method. We use an optimized wave function that is improved from an initial one-body wave function by multiplying by multiple correlation operators that are given by the form of exp($-S$)-type where $S$ is a suitable operator. We examine the inhomogeneous ground state at near 1/8 doping in the strongly correlated region where the on-site Coulomb interaction is larger than the bandwidth. The $d$-wave superconductivity with a spatially oscillating gap function can coexist with the charge ordering in the ground state without magnetic ordering, and superconducting condensation energy increases by this coexistence. We also show the $d$-wave pair correlation function as a function of lattice sites. The correlation function indicates that the long-range superconducting order indeed exists, and also that the strong electron correlation suppresses the pair correlation function.