Nonpertubative many-body theory for the two-dimensional Hubbard model at low temperature: From weak to strong coupling regimes

In theoretical studies of two-dimensional (2D) systems, the Mermin-Wagner theorem prevents continuous symmetry breaking at any finite temperature, thus forbidding a Landau phase transition at a critical temperature Tc T c. The difficulty arises when many-body theoretical studies predict a Landau phase transition at finite temperatures, which contradicts the Mermin-Wagner theorem and is termed a pseudo phase transition. To tackle this problem, we systematically develop a symmetrization scheme, defined as averaging physical quantities over all symmetry-breaking states, thus ensuring that it preserves the Mermin-Wagner theorem. We apply the symmetrization scheme to the GW-covariance calculation for the 2D repulsive Hubbard model at half-filling in the intermediate-to-strong coupling regime and at low temperatures, obtaining the one-body Green’s function and spin-spin correlation function, and benchmark them against Determinant Quantum Monte Carlo (DQMC) with good agreement. The spin-spin correlation functions are approached within the covariance theory, a general method for calculating two-body correlation functions from a one-particle starting point, such as the GW formalism used here, which ensures the preservation of the fundamental fluctuation-dissipation relation (FDR) and Ward-Takahashi identities (WTI). With the FDR and WTI satisfied, we conjecture that the χ -sum rule, a fundamental relation from the Pauli exclusion principle, can be used to probe the reliability of many-body methods, and demonstrate this by comparing the GW-covariance and mean-field-covariance approaches. This work provides a novel framework to investigate the strong-coupling and doped regime of the 2D Hubbard model, which is believed to be applicable to real high- Tc T c cuprate superconductors.