Interlayer Pairing in Bilayer Nickelates

Abstract: The discovery of $T_c\sim 80$~K superconductivity in pressurized La$3$Ni$_2$O$_7$ has launched a new platform to study high-temperature superconductivity. Using non-perturbative dynamic cluster approximation quantum Monte Carlo calculations, we characterize the magnetic and superconducting pairing behavior of a realistic bilayer two-orbital Hubbard-Hund model of this system that describes the relevant Ni $e_g$ states with physically relevant interaction strengths. We find a leading $s^\pm$ superconducting instability in this model and show that this state primarily arises from interlayer pairing in the $d{3z^2-r^2}$ orbital that is driven by strong interlayer spin-fluctuations in that orbital. These results provide non-perturbative evidence supporting the picture that a simple single-orbital bilayer Hubbard model for the Ni $d_{3z^2-r^2}$ orbital provides an excellent low-energy effective description of the superconducting behavior of La$_3$Ni$_2$O$_7$.