Solving 2D and 3D Lattice Models of Correlated Fermions—Combining Matrix Product States with Mean-Field Theory

Correlated electron states are at the root of many important phenomena including unconventional superconductivity (USC), where electron pairing arises from repulsive interactions. Computing the properties of correlated electrons, such as the critical temperature T_{c} for the onset of USC, efficiently and reliably from the microscopic physics with quantitative methods remains a major challenge for almost all models and materials. In this theoretical work, we combine matrix product states (MPS) with static mean field (MF) to provide a solution to this challenge for quasi-one-dimensional (Q1D) systems: two- and three-dimensional materials comprised of weakly coupled correlated 1D fermions. This MPS+MF framework for the ground state and thermal equilibrium properties of Q1D fermions is developed and validated for attractive Hubbard systems first, and further enhanced via analytical field theory. We then deploy it to compute T_{c} for superconductivity in 3D arrays of weakly coupled, doped, and repulsive Hubbard ladders. The MPS+MF framework thus enables the quantitative study of USC and high-T_{c} superconductivity—and potentially many more correlated phases—in fermionic Q1D systems based directly on their microscopic parameters, in ways inaccessible to previous methods. This approach further allows one to treat competing macroscopic orders, such as superconducting and insulating ones, on an equal footing. Benchmarks of the framework using auxiliary-field quantum Monte Carlo techniques show that the overestimation of, e.g., T_{c} due to its mean-field component, is near constant in microscopic parameters. These features of the MPS+MF approach to correlated fermions open up the possibility of designing deliberately optimized Q1D superconductors, from experiments in ultracold gases to synthesizing new materials.