| The two-dimensional (2D) Hubbard model is solved in a new way. The zero-frequency components of the Hubbard interaction are separated out and solved to infinite order to avoid the traditional difficulties associated with perturbation theory in the strong-coupling limit. Specifically the zero-frequency spin interaction is decoupled by the introduction of Stratonovich-Hubbard fields, while the zero-frequency charge interaction is treated by "slave" Hartree-Fock. This approach leads to the formation of large fluctuating spin moments, and as a consequence the one-electron spectral density exhibits a pseudogap. A Metropolis Monte Carlo method is used to perform the configuration averaging over Stratonovich-Hubbard fields. Within single Monte Carlo configurations the non-zero frequency components of the Hubbard interaction are treated by a high-order perturbation method, which incorporates the exchange of dynamical particle-hole fluctuations. Since a single configuration has neither translational symmetry nor spin rotational symmetry, these fluctuations do not have a well-defined crystal momentum or total spin. When the dynamical particle-hole fluctuations are included in the calculation, localized electron pairing results at elevated temperatures. Superconductivity presumably onsets when these localized pairs develop coherence, but this is not studied in the present work. The nature of the pairing interaction and the resulting pair eigenfunctions is studied. |
