The recently discovered nickelate superconductor
La3Ni2O7 has a high transition temperature near 80 K under pressure, providing an additional avenue for exploring unconventional superconductivity. Here, with state-of-the-art tensor-network methods, we study a bilayer ${t\text{\ensuremath{-}}J\text{\ensuremath{-}}J}_{\ensuremath{\perp}}$ model for
La3Ni2O7 and find a robust
s-wave superconductive (SC) order mediated by interlayer magnetic couplings. Large-scale density matrix renormalization group calculations find algebraic pairing correlations with Luttinger parameter ${K}_{\mathrm{SC}}\ensuremath{\lesssim}1$. Infinite projected entangled-pair state method obtains a nonzero SC order directly in the thermodynamic limit, and estimates a strong pairing strength ${\overline{\mathrm{\ensuremath{\Delta}}}}_{z}\ensuremath{\sim}\mathcal{O}(0.1)$. Tangent-space tensor renormalization group simulations elucidate the temperature evolution of SC pairing and further determine a high SC temperature ${T}_{c}^{*}/J\ensuremath{\sim}\mathcal{O}(0.1)$. Because of the intriguing orbital selective behaviors and strong Hund's rule coupling in the compound, ${t\text{\ensuremath{-}}J\text{\ensuremath{-}}J}_{\ensuremath{\perp}}$ model has strong interlayer spin exchange (while negligible interlayer hopping), which greatly enhances the SC pairing in the bilayer system. Such a magnetically mediated pairing has also been observed recently in the optical lattice of ultracold atoms. Our accurate and comprehensive tensor-network calculations reveal a robust SC order in the bilayer ${t\text{\ensuremath{-}}J\text{\ensuremath{-}}J}_{\ensuremath{\perp}}$ model and shed light on the pairing mechanism of the high-
Tc nickelate superconductor.
Motivated by the bilayer cuprate superconductors and nickelate superconductor La
3Ni
2O
7, we investigate the evolution from intralayer to interlayer superconductivity based on a bilayer two-leg
t-
J-
J⊥ model, where
t is the in-plane electron hopping,
J is the in-plane spin interaction, and
J⊥ is the inter-plane spin interaction. By means of the density matrix renormalization group calculations, we obtain the quantum phase diagram of the system by tuning
J⊥ in a large doping range
δ=1/8−1/2. We find that a large
J⊥ can always drive an interlayer superconductivity by coupling the two layers in both the Luther-Emery liquid and Luttinger liquid states. By coupling two Luther-Emery liquid states, the in-plane superconductivity evolves to inter-plane superconductivity either through an intermediate charge density wave (CDW) phase or directly, depending on doping ratio. This emergent CDW phase, which exists over a finite doping range, appears to develop from the CDW state of the two-leg ladder at
δ=1/4. By coupling two Luttinger liquids, the in-plane Luttinger liquids show a transition to the inter-plane superconducting phase at large
J⊥, as reported in previous literature. Interestingly, in the intermediate
J⊥ regime we find that while the in-plane Luttinger-liquid features remain stable, the inter-plane superconductivity can develop an enhanced quasi-long-range order with the power exponent
KSCzz∼1. At last, we show that the interlayer superconductivity is also stable by coupling the bilayer three-leg
t-
J ladders by a strong
J⊥ interaction, from both the Luther-Emery liquid and Luttinger-liquid states.
Abstract Ever since the discovery of high-temperature superconductivity in cuprates, gaining microscopic insights into the nature of pairing in strongly correlated systems has remained one of the greatest challenges in modern condensed matter physics. Following recent experiments reporting superconductivity in the bilayer nickelate La3Ni2O7 (LNO) with remarkably high critical temperatures of Tc = 80 K, it has been argued that the low-energy physics of LNO can be described by the strongly correlated, mixed-dimensional bilayer t–J model. Here we investigate this bilayer system and utilize density matrix renormalization group techniques to establish a thorough understanding of the model and the magnetically induced pairing through comparison to the perturbative limit of dominating inter-layer spin couplings. In particular, this allows us to explain appearing finite-size effects, firmly establishing the existence of long-range superconducting order in the thermodynamic limit. By analyzing binding energies, we predict a BEC–BCS crossover as a function of the Hamiltonian parameters. We find large binding energies of the order of the inter-layer coupling that suggest strikingly high critical temperatures of the Berezinskii–Kosterlitz–Thouless transition, raising the question of whether (mixD) bilayer superconductors possibly facilitate critical temperatures above room temperature.
Competing orders represent a central challenge in understanding strongly correlated systems. In this work, we employ projector quantum Monte Carlo simulations to study a sign-problem-free bilayer extended Hubbard model. In this model, a charge stripe phase, characterized by a peak at (0±2πδ,0) is induced by highly anisotropic interlayer antiferromagnetic spin-exchange coupling Jz, and strongly suppressed by the decreasing spin-exchange anisotropy (i.e., introducing the spin-flip term J⊥); in contrast, the introduction of J⊥ favors the emergence of interlayer pairing superconductivity. We further demonstrate that the anisotropy of the interlayer spin-exchange directly governs the competition between these two phases, while the on-site interaction U plays a complex role in tuning both the charge stripe and superconductivity. Furthermore, against the background of charge stripes we discover the spin stripes, characterized by twice the period of charge stripes, indicating the profound internal connection between them. Our work identifies the key factors driving charge stripe formation, highlights the sensitivity of both the charge stripe and superconducting phases to interaction parameters, and thereby provides valuable insights into competing orders in strongly correlated systems.
Abstract In this article, we review some recent theoretical developments on potential high-temperature superconductors and unconventional metallic states that can arise from doping a spin-one Mott insulator in the d 8 valence. These studies are particularly relevant—though not limited—to the recently discovered bilayer nickelate superconductor La 3 Ni 2 O 7 . We focus on a ferromagnetic Kondo lattice model with mobile electrons in the d x 2 − y 2 orbital coupled to the localized spin moments in d z 2 orbital through a large Hund’s coupling J H . In the large J H limit, the model reduces to the type II t − J model with a mixture of spin-half singlon states and spin-one doublon states. We summarize density matrix renormalization group results on the Luther-Emery liquid in one dimensional chain and two-leg ladder. Then we mainly focus on bilayer square lattice and show that a large inter-layer coupling J ⊥ of d z 2 orbital can induce strong inter-layer pairing of d x 2 − y 2 orbital. In the strong J ⊥ limit, a kinetic-energy driven high T c superconductivity is demonstrated in an ideal model with only a single hopping term. Furthermore, the model predicts a symmetric pseudogap metal—dubbed ‘second Fermi liquid (FL)’—in the underdoped regime, yielding a phase diagram analogous to that of hole-doped cuprates. The bilayer Kondo model therefore presents a promising platform for both realizing higher- T c superconductors and exploring non-FL physics. We also comment on the possible limitations of the current models for the bilayer nickelate material and point out some future directions.
Motivated by the recent discovery of superconductivity in the bilayer nickelate La
3Ni
2O
7 (LNO) under pressure, we study a mixed-dimensional (mixD) bilayer
t∥-
J∥-
J⊥ model, which has been proposed as an effective low-energy description of LNO. Using neural quantum states (NQS), and in particular Gutzwiller-projected Hidden Fermion Pfaffian State, we access the ground-state properties on large lattices up to
8×8×2 sites. We show that this model exhibits superconductivity across a wide range of dopings and couplings, and analyze the pairing behavior in detail. We identify a crossover from tightly bound, Bose-Einstein-condensed interlayer pairs at strong interlayer exchange to more spatially extended Bardeen-Cooper-Schrieffer-like pairs as the interlayer exchange is decreased. Furthermore, upon tuning the intralayer exchange, we observe a sharp transition from interlayer
s-wave pairing to intralayer
d-wave pairing, consistent with a first-order change in the pairing symmetry. We verify that our simulations are accurate by comparing with matrix product state simulations on coupled ladders. Our results represent the first simulation of a fermionic multi-orbital system with NQS, and provide the first evidence for superconductivity in two-dimensonal bilayers using high-precision numerics. These findings provide insight into superconductivity in bilayer nickelates and cold atom quantum simulation platforms.