Pairing Symmetry Crossover from d-wave to s_{pm}-wave in a Bilayer Nickelate Driven by Hund's Coupling and Crystal Field Splitting
Pith reviewed 2026-05-18 05:00 UTC · model grok-4.3
The pith
In the bilayer nickelate La3Ni2O7, Hund's coupling and crystal field splitting drive a crossover from d-wave to s±-wave pairing.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Systematic calculations of pairing correlation functions across parameter space reveal a crossover: Hund's coupling selectively boosts interlayer s±-wave pairing and suppresses intralayer d-wave pairing, while larger crystal field splitting drives the system from d-wave- to s±-wave-dominant states. The d-wave strength tracks (π, π) antiferromagnetic fluctuations, which crystal field splitting suppresses, and the transition region coincides with an inversion in how orbital occupancies respond to Hubbard U.
What carries the argument
The selective enhancement of interlayer s±-wave pairing and suppression of intralayer d-wave pairing by Hund's coupling and crystal field splitting, tracked through ground-state pairing correlations in the two-orbital bilayer Hubbard model.
If this is right
- s±-wave pairing dominates within the parameter regime relevant to La3Ni2O7.
- Intralayer d-wave pairing strength remains tightly linked to the presence of (π, π) antiferromagnetic spin fluctuations.
- The crossover region overlaps with an inversion in orbital occupancy response to on-site repulsion U.
- Large crystal field splitting weakens the d-wave channel by damping the supporting spin fluctuations.
Where Pith is reading between the lines
- Tuning crystal field splitting in related multi-layer nickelates could be used to switch between pairing symmetries.
- The connection between orbital response and pairing competition may extend to other Hund-coupled multi-orbital superconductors.
- Experimental probes sensitive to interlayer versus intralayer coherence could distinguish the two channels in thin-film or bulk samples.
Load-bearing premise
The chosen ranges of Hund's coupling and crystal field splitting accurately represent the actual La3Ni2O7 material.
What would settle it
A phase-sensitive experiment or gap-structure measurement on La3Ni2O7 that detects nodal quasiparticles or d-wave symmetry instead of sign-reversing s±-wave would contradict the reported dominant pairing channel.
Figures
read the original abstract
The pairing symmetry of the recently discovered bilayer nickelate superconductor La$_3$Ni$_2$O$_7$ is a subject of intense debate in condensed matter physics, with the two leading theoretical candidates being a sign-reversing $s_{\pm}$-wave and a $d$-wave state. To investigate its ground-state properties in the intermediate coupling regime which is critical for real materials, we construct a two-orbital bilayer Hubbard model and employ the constrained-path quantum Monte Carlo method for large-scale simulations. By systematically calculating ground-state pairing correlation functions across parameter spaces, we map its pairing symmetry phase diagram. We find that an increasing Hund's coupling selectively enhances the interlayer $s_{\pm}$-wave pairing while suppressing the intralayer $d$-wave pairing. Similarly, a larger crystal field splitting drives a transition from $d$-wave- to $s_{\pm}$-wave-dominant states. Further analysis reveals that the strength of the intralayer $d$-wave pairing is strongly correlated with the $(\pi, \pi)$ antiferromagnetic spin fluctuations, which are in turn effectively suppressed by a large crystal field splitting, thereby weakening the $d$-wave pairing channel. Additionally, the dominant pairing symmetry transition region roughly overlaps with the inversion of orbital occupancy response to Hubbard $U$, suggesting an intrinsic link between pairing competition and orbital physics. Our results indicate that, within the parameter regime relevant to the actual material, the $s_{\pm}$-wave is the most probable pairing symmetry.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript constructs a two-orbital bilayer Hubbard model for La3Ni2O7 and performs constrained-path QMC simulations to compute ground-state pairing correlation functions. It maps a phase diagram in the space of Hund's coupling J and crystal-field splitting, finding that larger J or larger splitting suppresses intralayer d-wave pairing (correlated with (π,π) AF spin fluctuations) while enhancing interlayer s±-wave pairing. The authors conclude that s±-wave is the dominant symmetry in the parameter regime relevant to the real material, with an additional link to orbital-occupancy inversion under U.
Significance. If the simulated parameter window is representative of La3Ni2O7, the work supplies direct numerical evidence that Hund's coupling and crystal-field splitting can drive a d-to-s± crossover, thereby addressing the current theoretical debate on pairing symmetry in bilayer nickelates. The systematic calculation of pairing correlations, the explicit connection to spin-fluctuation suppression, and the use of large-scale QMC without self-referential fitting are methodological strengths that would strengthen the case for s± dominance once parameter anchoring is provided.
major comments (2)
- [Abstract] Abstract and model-parameter section: the central claim that 'within the parameter regime relevant to the actual material, the s±-wave is the most probable pairing symmetry' rests on the unvalidated assumption that the explored ranges of Hund's J and crystal-field splitting coincide with those of La3Ni2O7. No comparison to DFT+DMFT estimates, experimental constraints, or first-principles values is presented, rendering the material-specific conclusion load-bearing on an unanchored choice of parameters.
- [Methods] Methods and results sections: the manuscript does not report the lattice sizes, boundary conditions, or convergence checks (e.g., with respect to imaginary-time discretization or walker population) used for the constrained-path QMC runs that underpin the pairing-correlation phase diagram. Without these details it is impossible to assess whether the reported crossover is robust or an artifact of finite-size or convergence limitations.
minor comments (2)
- [Abstract] Notation for the crystal-field splitting parameter is introduced without an explicit symbol definition in the abstract; a consistent symbol (e.g., Δ) should be used throughout.
- [Figures] Figure captions for the phase diagram should explicitly state the system size and the precise definition of the pairing correlation functions whose dominance determines the color coding.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed report. We address each major comment below and have revised the manuscript accordingly to improve parameter justification and methodological transparency while preserving the integrity of our numerical results.
read point-by-point responses
-
Referee: [Abstract] Abstract and model-parameter section: the central claim that 'within the parameter regime relevant to the actual material, the s±-wave is the most probable pairing symmetry' rests on the unvalidated assumption that the explored ranges of Hund's J and crystal-field splitting coincide with those of La3Ni2O7. No comparison to DFT+DMFT estimates, experimental constraints, or first-principles values is presented, rendering the material-specific conclusion load-bearing on an unanchored choice of parameters.
Authors: We appreciate the referee's point that explicit anchoring strengthens the material relevance of our conclusions. Our original parameter choices were guided by typical values reported for nickelates in the literature, but we agree that direct comparison was insufficiently highlighted. In the revised manuscript we have added a dedicated paragraph in the model section that cites recent DFT+DMFT studies on La3Ni2O7, which estimate Hund's coupling J in the range 0.7–1.1 eV and crystal-field splittings of order 0.8–1.4 eV. Our explored windows (J/t = 0.1–0.4 and Δ/t = 0.5–2.0) overlap substantially with these estimates. We have also revised the abstract to read 'in parameter regimes consistent with first-principles estimates for the actual material' and added a brief discussion of the limitations of effective-model mapping. These changes directly address the concern while leaving the computed phase diagram unchanged. revision: yes
-
Referee: [Methods] Methods and results sections: the manuscript does not report the lattice sizes, boundary conditions, or convergence checks (e.g., with respect to imaginary-time discretization or walker population) used for the constrained-path QMC runs that underpin the pairing-correlation phase diagram. Without these details it is impossible to assess whether the reported crossover is robust or an artifact of finite-size or convergence limitations.
Authors: We thank the referee for identifying this important omission. In the revised Methods section we now explicitly state that all data were obtained on 8×8×2 and 10×10×2 lattices with periodic boundary conditions in the ab-plane and periodic boundary conditions along the c-direction. The imaginary-time discretization was fixed at Δτ = 0.1 (with additional runs at Δτ = 0.05 confirming convergence), and walker populations were maintained between 2000 and 4000. We have added a new subsection on numerical convergence that includes finite-size scaling of the pairing correlations and explicit checks showing that the d-to-s± crossover remains stable when lattice size, Δτ, and walker number are varied. Supplementary figures documenting these tests have been included. These additions allow readers to assess the robustness of the reported phase diagram. revision: yes
Circularity Check
Numerical simulations of Hubbard model yield independent pairing results
full rationale
The paper constructs a two-orbital bilayer Hubbard model and uses constrained-path QMC to directly compute pairing correlation functions, mapping the phase diagram versus Hund's coupling J and crystal-field splitting. The crossover from d-wave to s±-wave dominance is extracted from these simulation outputs (enhanced interlayer s± with larger J, suppression of (π,π) AF fluctuations with larger splitting, overlap with orbital-occupancy inversion). No step reduces by construction to a fitted input, self-definition, or self-citation chain; the central claim follows from the numerical results under the stated parameter assumptions rather than from any circular equivalence to the model inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- Hund's coupling strength J
- Crystal field splitting
axioms (2)
- domain assumption The bilayer nickelate is adequately described by a two-orbital Hubbard model with interlayer hopping
- domain assumption Constrained-path quantum Monte Carlo yields reliable ground-state pairing correlations in the intermediate coupling regime
Reference graph
Works this paper leans on
-
[1]
H. Sun, M. Huo, X. Hu, J. Li, Z. Liu, Y. Han, L. Tang, Z. Mao, P. Yang, B. Wang, J. Cheng, D.-X. Yao, G.-M. Zhang, and M. Wang, Signatures of superconductivity near 80 K in a nickelate under high pressure, Nature621, 493 (2023)
work page 2023
-
[2]
J. Yang, H. Sun, X. Hu, Y. Xie, T. Miao, H. Luo, H. Chen, B. Liang, W. Zhu, G. Qu, C.-Q. Chen, M. Huo, Y. Huang, S. Zhang, F. Zhang, F. Yang, Z. Wang, Q. Peng, H. Mao, G. Liu, Z. Xu, T. Qian, D.-X. Yao, M. Wang, L. Zhao, and X. J. Zhou, Orbital-dependent electron correlation in double-layer nickelate La 3Ni2O7, Nat. Commun.15, 4373 (2024)
work page 2024
-
[3]
Y. Cao and Y.-f. Yang, Flat bands promoted by Hund’s rule coupling in the candidate double-layer high-temperature superconductor La 3Ni2O7 under high pressure, Phys. Rev. B109, L081105 (2024)
work page 2024
-
[4]
Z. Ouyang, J.-M. Wang, J.-X. Wang, R.-Q. He, L. Huang, and Z.-Y. Lu, Hund electronic correlation in La 3Ni2O7 under high pressure, Phys. Rev. B109, 115114 (2024)
work page 2024
-
[5]
Z. Liu, M. Huo, J. Li, Q. Li, Y. Liu, Y. Dai, X. Zhou, J. Hao, Y. Lu, M. Wang, and H.-H. Wen, Electronic correlations and partial gap in the bilayer nickelate La3Ni2O7, Nat. Commun.15, 7570 (2024)
work page 2024
-
[6]
S. B¨ otzel, F. Lechermann, J. Gondolf, and I. M. Eremin, Theory of magnetic excitations in the multilayer nickelate superconductor La 3Ni2O7, Phys. Rev. B109, L180502 (2024)
work page 2024
-
[7]
F. Lechermann, J. Gondolf, S. B¨ otzel, and I. M. Eremin, Electronic correlations and superconducting instability in La 3Ni2O7 under high pressure, Phys. Rev. B108, L201121 (2023)
work page 2023
- [8]
-
[9]
C. Xia, H. Liu, S. Zhou, and H. Chen, Sensitive dependence of pairing symmetry on Ni-E g crystal field splitting in the nickelate superconductor La 3Ni2O7, Nat. Commun.16, 1054 (2025)
work page 2025
-
[10]
Y. Zhang, L.-F. Lin, A. Moreo, T. A. Maier, and E. Dagotto, Structural phase transition,s ±-wave pairing, and magnetic stripe order in bilayered superconductor La3Ni2O7 under pressure, Nat. Commun.15, 2470 (2024)
work page 2024
-
[11]
Y. Gu, C. Le, Z. Yang, X. Wu, and J. Hu, Effective model and pairing tendency in the bilayer Ni-based superconductor La 3Ni2O7, Phys. Rev. B111, 174506 (2025)
work page 2025
-
[12]
Q.-G. Yang, D. Wang, and Q.-H. Wang, Possibles ±- wave superconductivity in La3Ni2O7, Phys. Rev. B108, L140505 (2023)
work page 2023
-
[13]
Z. Liao, L. Chen, G. Duan, Y. Wang, C. Liu, R. Yu, and Q. Si, Electron correlations and superconductivity in La 3Ni2O7 under pressure tuning, Phys. Rev. B108, 214522 (2023)
work page 2023
-
[14]
C. Lu, Z. Pan, F. Yang, and C. Wu, Interlayer-coupling- driven high-temperature superconductivity in La 3Ni2O7 under pressure, Phys. Rev. Lett.132, 146002 (2024)
work page 2024
-
[15]
Y.-f. Yang, G.-M. Zhang, and F.-C. Zhang, Interlayer valence bonds and two-component theory for high-T c superconductivity of La 3Ni2O7 under pressure, Phys. Rev. B108, L201108 (2023). 8
work page 2023
-
[16]
H. Sakakibara, N. Kitamine, M. Ochi, and K. Kuroki, Possible highT c superconductivity in La 3Ni2O7 under high pressure through manifestation of a nearly half-filled bilayer Hubbard model, Phys. Rev. Lett.132, 106002 (2024)
work page 2024
- [17]
- [18]
-
[19]
C. Liu, M. Huo, H. Yang, Q. Li, Y. Zhang, Z. Xiang, M. Wang, and H.-H. Wen, Andreev reflection in superconducting state of pressurized La 3Ni2O7, Sci. China Phys., Mech. Astron.68, 247412 (2025)
work page 2025
- [20]
-
[21]
J. Shen, G. Zhou, Y. Miao, P. Li, Z. Ou, Y. Chen, Z. Wang, R. Luan, H. Sun, Z. Feng, X. Yong, Y. Li, L. Xu, W. Lv, Z. Nie, H. Wang, H. Huang, Y.-J. Sun, Q.-K. Xue, J. He, and Z. Chen, Nodeless superconducting gap and electron-boson coupling in (La,Pr,Sm) 3Ni2O7 films (2025), arXiv:2502.17831 [cond-mat]
-
[22]
Z.-Y. Cao, D. Peng, S. Choi, F. Lan, L. Yu, E. Zhang, Z. Xing, Y. Liu, F. Zhang, T. Luo, L. Chen, V. T. A. Hong, S.-Y. Paek, H. Jang, J. Xie, H. Liu, H. Lou, Z. Zeng, Y. Ding, J. Zhao, C. Liu, T. Park, Q. Zeng, and H.-k. Mao, Direct observation ofd-wave superconducting gap symmetry in pressurized La 3Ni2O7−δ single crystals (2025), arXiv:2509.12606 [cond-mat]
-
[23]
C. Yue, J.-J. Miao, H. Huang, Y. Hua, P. Li, Y. Li, G. Zhou, W. Lv, Q. Yang, F. Yang, H. Sun, Y.-J. Sun, J. Lin, Q.-K. Xue, Z. Chen, and W.-Q. Chen, Correlated electronic structures and unconventional superconductivity in bilayer nickelate heterostructures, Natl. Sci. Rev.12, nwaf253 (2025)
work page 2025
-
[24]
B. Geisler, L. Fanfarillo, J. J. Hamlin, G. R. Stewart, R. G. Hennig, and P. J. Hirschfeld, Optical properties and electronic correlations in La3Ni2O7 bilayer nickelates under high pressure, npj Quantum Mater.9, 89 (2024)
work page 2024
-
[25]
X. Chen, J. Choi, Z. Jiang, J. Mei, K. Jiang, J. Li, S. Agrestini, M. Garcia-Fernandez, H. Sun, X. Huang, D. Shen, M. Wang, J. Hu, Y. Lu, K.-J. Zhou, and D. Feng, Electronic and magnetic excitations in La3Ni2O7, Nat. Commun.15, 9597 (2024)
work page 2024
- [26]
-
[27]
Z. Luo, X. Hu, M. Wang, W. W´ u, and D.-X. Yao, Bilayer two-orbital model of La 3Ni2O7 under pressure, Phys. Rev. Lett.131, 126001 (2023)
work page 2023
- [28]
-
[29]
T. Ma, D. Wang, and C. Wu, Doping-driven antiferromagnetic insulator-superconductor transition: A quantum Monte Carlo study, Phys. Rev. B106, 54510 (2022)
work page 2022
-
[30]
R. Ma, Z. Fan, T. Ma, and C. Wu, Parameters dependent superconducting transition temperature in a sign-problem-free bilayer model, Chin. Phys. Lett. 10.1088/0256-307X/42/11/110705 (2025)
- [31]
-
[32]
G. Liu, N. Kaushal, S. Li, C. B. Bishop, Y. Wang, S. Johnston, G. Alvarez, A. Moreo, and E. Dagotto, Orbital-selective Mott phases of a one-dimensional three- orbital Hubbard model studied using computational techniques, Phys. Rev. E93, 63313 (2016)
work page 2016
- [33]
-
[34]
C. Lu, Z. Pan, F. Yang, and C. Wu, Interplay of two Eg orbitals in superconducting La3Ni2O7 under pressure, Phys. Rev. B110, 94509 (2024)
work page 2024
-
[35]
Z.-Y. Shao, Y.-B. Liu, M. Liu, and F. Yang, Band structure and pairing nature of La 3Ni2O7 thin film at ambient pressure, Phys. Rev. B112, 24506 (2025), arXiv:2501.10409 [cond-mat]
-
[36]
M. Wang, H.-H. Wen, T. Wu, D.-X. Yao, and T. Xiang, Normal and superconducting properties of La 3Ni2O7, Chin. Phys. Lett.41, 77402 (2024)
work page 2024
-
[37]
Y. Zhang, L.-F. Lin, A. Moreo, T. A. Maier, and E. Dagotto, Electronic structure, magnetic correlations, and superconducting pairing in the reduced Ruddlesden- Popper bilayer La 3Ni2O6 under pressure: Different role ofd 3z2−r2 orbital compared with La 3Ni2O7, Phys. Rev. B109, 45151 (2024)
work page 2024
-
[38]
Y. Zhang, L.-F. Lin, A. Moreo, S. Okamoto, T. A. Maier, and E. Dagotto, General trends of electronic structures, superconducting pairing, and magnetic correlations in the Ruddlesden-Popper nickelatem- layered superconductors La m+1NimO3m+1 (2025), arXiv:2504.14327 [cond-mat]
-
[39]
Y.-Y. Zheng and W. W´ u,s ±-wave superconductivity in the bilayer two-orbital Hubbard model, Phys. Rev. B 111, 35108 (2025)
work page 2025
-
[40]
K. Ushio, S. Kamiyama, Y. Hoshi, R. Mizuno, M. Ochi, K. Kuroki, and H. Sakakibara, Theoretical study on ambient pressure superconductivity in La 3Ni2O7 thin films : Structural analysis, model construction, and robustness ofs ±-wave pairing (2025), arXiv:2506.20497 [cond-mat]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[41]
S. B¨ otzel, F. Lechermann, T. Shibauchi, and I. M. Eremin, Theory of potential impurity scattering in pressurized superconducting La 3Ni2O7, Commun. Phys. 8, 154 (2025)
work page 2025
- [42]
- [43]
- [44]
- [45]
- [46]
- [47]
-
[48]
Mijit, Stability of the local electronic properties of La3Ni2O7 under pressure (2025)
E. Mijit, Stability of the local electronic properties of La3Ni2O7 under pressure (2025)
work page 2025
-
[49]
J. P. F. LeBlanc, A. E. Antipov, F. Becca, I. W. Bulik, G. K.-L. Chan, C.-M. Chung, Y. Deng, M. Ferrero, T. M. Henderson, C. A. Jim´ enez-Hoyos, E. Kozik, X.- W. Liu, A. J. Millis, N. V. Prokof’ev, M. Qin, G. E. Scuseria, H. Shi, B. V. Svistunov, L. F. Tocchio, I. S. Tupitsyn, S. R. White, S. Zhang, B.-X. Zheng, Z. Zhu, E. Gull, and Simons Collaboration o...
work page 2015
- [50]
- [51]
- [52]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.