d-Wave pair density wave superconductivity in a two-orbital model
Pith reviewed 2026-05-10 17:00 UTC · model grok-4.3
The pith
Two-orbital models on the square lattice exhibit instability to incommensurate d-wave pair density wave superconductivity driven by interband pairing.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For minimal interorbital t-J or t-V interactions in two-orbital models, a random phase approximation analysis uncovers regimes of instability towards incommensurate d_xy pair density wave superconductivity driven by interband pairing. At strong coupling, the derived effective hard-core Cooper pair Hamiltonian, treated with a bosonic Gutzwiller ansatz, yields a period-2 PDW over a wide range of fillings and a checkerboard CDW at quarter-filling. The orbital content and multiband nature of the Fermi surface play a key role in stabilizing these interband PDW states.
What carries the argument
Interband pairing interactions that stabilize incommensurate d_xy pair density wave order in two-orbital lattice models.
If this is right
- PDW superconductivity emerges as the leading instability for specific fillings and interaction strengths.
- This modulated order competes with uniform d_xy pairing and with magnetic and charge density wave states.
- The strong-coupling effective Hamiltonian supports a period-2 PDW state across a broad range of fillings.
- A checkerboard charge density wave stabilizes at quarter filling.
- The results apply to correlated materials with quasi-1D bands and to atomic fermions in p-orbitals.
Where Pith is reading between the lines
- Orbital-selective pairing could be tuned to favor PDW order in other multi-orbital systems not studied here.
- Experiments with ultracold atoms loaded into p-orbitals on optical lattices might directly realize the predicted period-2 PDW.
- Varying the orbital hybridization strength may switch the ground state between PDW and uniform superconductivity.
Load-bearing premise
The random phase approximation remains reliable for identifying the leading instability when only minimal interorbital t-J or t-V interactions are present, and the bosonic Gutzwiller ansatz accurately captures the ground state of the derived hard-core pair model.
What would settle it
A calculation or measurement showing no instability to incommensurate d-PDW in the RPA for these models, or finding no period-2 PDW in the strong-coupling limit over wide fillings, would falsify the central claims.
Figures
read the original abstract
Motivated by exploring superconductivity in multi-orbital systems, we study two orbital models of spinful fermions representing ($p_x,p_y$) or ($d_{xz}, d_{yz})$ orbitals on the square lattice. For minimal interorbital $t$-$J$ or $t$-$V$ on-site interactions, a random phase approximation uncovers regimes of instability towards incommensurate $d_{xy}$ pair density wave ($d$-PDW) superconductivity with driven by interband pairing. We study the competition of PDW order with uniform nodal $d_{xy}$ pairing states and magnetic and charge density wave (CDW) instabilities. At strong coupling, we derive an effective hard-core Cooper pair Hamiltonian which we study using a bosonic Gutzwiller ansatz to reveal a period-$2$ PDW over a wide range of fillings as well as a checkerboard CDW at quarter-filling. Our results apply to correlated multi-orbital materials with quasi-1D bands, Hubbard models on the square-octagon lattice, and atomic fermions in $p$-orbitals. Our work highlights the role of the orbital content and multiband Fermi surfaces in stabilizing interband PDW states.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies two-orbital models of spinful fermions on the square lattice (representing (p_x,p_y) or (d_xz,d_yz) orbitals) with minimal interorbital t-J or t-V on-site interactions. RPA calculations identify regimes of instability to incommensurate d_xy pair-density-wave (d-PDW) superconductivity driven by interband pairing, and compare this to uniform nodal d_xy SC, magnetic, and CDW channels. At strong coupling an effective hard-core Cooper-pair Hamiltonian is derived and solved via bosonic Gutzwiller ansatz, yielding a period-2 PDW over a wide filling range plus checkerboard CDW at quarter filling. The results are positioned as relevant to quasi-1D multi-orbital materials, square-octagon Hubbard models, and p-orbital atomic fermions.
Significance. If the reported instability ordering and strong-coupling ground states hold, the work usefully illustrates how orbital content and multiband Fermi surfaces can stabilize interband PDW states. The explicit mapping from the microscopic model to a hard-core pair Hamiltonian is a concrete strength that enables further study and falsifiable predictions for filling-dependent PDW periodicity.
major comments (2)
- [§3] §3 (RPA susceptibility analysis): the identification of incommensurate d_xy PDW as the leading instability rests on the divergence of the RPA pairing bubble driven by interband terms, yet the manuscript provides no estimate or argument for the size of omitted vertex corrections or self-energy feedback. Because the interactions are on-site t-J/t-V (strong-coupling projections), these corrections can reorder PDW versus uniform d_xy SC and CDW channels; this directly undermines the claimed regimes of PDW instability.
- [§4] §4 (strong-coupling mapping and Gutzwiller ansatz): the period-2 PDW result over wide fillings is obtained from the bosonic Gutzwiller treatment of the derived hard-core pair Hamiltonian. No benchmark against exact diagonalization on small clusters or against alternative variational states is reported, so the quantitative reliability of the filling range and the competition with checkerboard CDW at quarter filling cannot be assessed.
minor comments (2)
- [Figure 3] Notation for the incommensurate wave-vector Q in the PDW susceptibility should be defined explicitly in the text and figure captions rather than only in the appendix.
- [§3] The abstract states that the PDW is 'driven by interband pairing'; the corresponding sentence in §3 should cite the explicit interorbital matrix element that produces this channel.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. We address each major comment below and have incorporated revisions to clarify the limitations of our methods while maintaining the core findings of the manuscript.
read point-by-point responses
-
Referee: §3 (RPA susceptibility analysis): the identification of incommensurate d_xy PDW as the leading instability rests on the divergence of the RPA pairing bubble driven by interband terms, yet the manuscript provides no estimate or argument for the size of omitted vertex corrections or self-energy feedback. Because the interactions are on-site t-J/t-V (strong-coupling projections), these corrections can reorder PDW versus uniform d_xy SC and CDW channels; this directly undermines the claimed regimes of PDW instability.
Authors: We agree that vertex corrections and self-energy feedback are not estimated in our RPA calculations, and this is a known limitation of the RPA approach, especially for strong-coupling interactions like t-J and t-V. In the revised manuscript, we have added a paragraph in Section 3 discussing this point, noting that while such corrections could potentially reorder the instabilities, the interband pairing mechanism identified here is robust and is corroborated by the independent strong-coupling analysis in Section 4. We argue that the PDW instability remains a viable candidate in the parameter regimes studied. revision: partial
-
Referee: §4 (strong-coupling mapping and Gutzwiller ansatz): the period-2 PDW result over wide fillings is obtained from the bosonic Gutzwiller treatment of the derived hard-core pair Hamiltonian. No benchmark against exact diagonalization on small clusters or against alternative variational states is reported, so the quantitative reliability of the filling range and the competition with checkerboard CDW at quarter filling cannot be assessed.
Authors: We acknowledge the absence of benchmarks against exact diagonalization or other variational methods, which limits the quantitative assessment of the results. The bosonic Gutzwiller ansatz provides a variational upper bound to the ground-state energy and is known to accurately describe the superfluid and density-wave phases in hard-core boson models on lattices. In the revision, we have expanded the discussion in Section 4 to include a more detailed justification of the method, its expected accuracy for the PDW state, and a note on the competition with CDW at quarter filling. We believe this addresses the concern while the qualitative features, such as the wide filling range for period-2 PDW, are supported by the consistency with RPA. revision: partial
Circularity Check
Derivation chain self-contained with no circular reductions
full rationale
The paper starts from a microscopic two-orbital model with specified interorbital t-J or t-V interactions, applies RPA linear response to extract pairing susceptibilities and identify leading instabilities (including incommensurate d_xy PDW), then performs an explicit strong-coupling projection to obtain an effective hard-core Cooper-pair Hamiltonian whose ground state is analyzed via bosonic Gutzwiller ansatz. None of these steps reduce by the paper's own equations to a fitted parameter renamed as a prediction, a self-citation load-bearing uniqueness claim, or an ansatz smuggled in from prior work; the reported PDW regimes and period-2 order emerge directly from the response functions and the derived effective model without tautological equivalence to the inputs.
Axiom & Free-Parameter Ledger
free parameters (1)
- interorbital interaction strengths (J, V)
axioms (2)
- domain assumption Random phase approximation captures the leading superconducting instability
- domain assumption Strong-coupling projection to hard-core Cooper pairs is valid
Reference graph
Works this paper leans on
-
[1]
D. P. Arovas, E. Berg, S. A. Kivelson, and S. Raghu, The hubbard model, Annual Review of Condensed Mat- ter Physics13, 239 (2022)
work page 2022
-
[2]
Dagotto, Correlated electrons in high-temperature su- perconductors, Rev
E. Dagotto, Correlated electrons in high-temperature su- perconductors, Rev. Mod. Phys.66, 763 (1994)
work page 1994
-
[3]
P. A. Lee, N. Nagaosa, and X.-G. Wen, Doping a mott in- sulator: Physics of high-temperature superconductivity, Rev. Mod. Phys.78, 17 (2006)
work page 2006
-
[4]
P. W. Anderson, P. A. Lee, M. Randeria, T. M. Rice, N. Trivedi, and F. C. Zhang, The physics behind high- temperature superconducting cuprates: the plain vanilla version of rvb, Journal of Physics: Condensed Matter16, R755 (2004)
work page 2004
-
[5]
D. J. Scalapino, Superconductivity and spin fluctuations (1999), cond-mat/9908287
work page internal anchor Pith review arXiv 1999
- [6]
-
[7]
A. V. Chubukov, D. V. Efremov, and I. Eremin, Mag- netism, superconductivity, and pairing symmetry in iron- based superconductors, Phys. Rev. B78, 134512 (2008)
work page 2008
-
[8]
A. Chubukov, Pairing mechanism in fe-based supercon- ductors, Annual Review of Condensed Matter Physics3, 57 (2012)
work page 2012
-
[9]
V. J. Emery, Theory of high-t c superconductivity in ox- ides, Phys. Rev. Lett.58, 2794 (1987)
work page 1987
-
[10]
V. J. Emery and G. Reiter, Mechanism for high- temperature superconductivity, Phys. Rev. B38, 4547 (1988)
work page 1988
-
[11]
C. M. Varma, Non-fermi-liquid states and pairing insta- bility of a general model of copper oxide metals, Phys. Rev. B55, 14554 (1997)
work page 1997
-
[12]
H. Watanabe, T. Shirakawa, K. Seki, H. Sakakibara, T. Kotani, H. Ikeda, and S. Yunoki, Unified description of cuprate superconductors using a four-bandd−pmodel, Phys. Rev. Res.3, 033157 (2021)
work page 2021
-
[13]
P. Mai, G. Balduzzi, S. Johnston, and T. A. Maier, Or- bital structure of the effective pairing interaction in the high-temperature superconducting cuprates, npj Quan- tum Materials6, 26 (2021)
work page 2021
-
[14]
H. Miyahara, R. Arita, and H. Ikeda, Development of a two-particle self-consistent method for multiorbital sys- tems and its application to unconventional superconduc- tors, Phys. Rev. B87, 045113 (2013)
work page 2013
-
[15]
O. Gingras, N. Allaglo, R. Nourafkan, M. Cˆ ot´ e, and A.- M. S. Tremblay, Superconductivity in correlated multi- orbital systems with spin-orbit coupling: Coexistence of even- and odd-frequency pairing, and the case of sr2ruo4, Phys. Rev. B106, 064513 (2022)
work page 2022
-
[16]
A. C. Yuan, E. Berg, and S. A. Kivelson, Multiband mean-field theory of thed+igsuperconductivity scenario in sr2ruo4, Phys. Rev. B108, 014502 (2023)
work page 2023
-
[17]
C.-Y. Moon, Effects of orbital selective dynamical cor- relation on the spin susceptibility and superconduct- ing symmetries in sr 2ruo4, Phys. Rev. Res.5, L022058 (2023)
work page 2023
-
[18]
H. Suzuki, L. Wang, J. Bertinshaw, H. U. R. Strand, S. K¨ aser, M. Krautloher, Z. Yang, N. Wentzell, O. Par- collet, F. Jerzembeck, N. Kikugawa, A. P. Mackenzie, A. Georges, P. Hansmann, H. Gretarsson, and B. Keimer, Distinct spin and orbital dynamics in Sr2RuO4, Nature Communications14, 7042 (2023)
work page 2023
-
[19]
Q. Si, R. Yu, and E. Abrahams, High-temperature super- conductivity in iron pnictides and chalcogenides, Nature Reviews Materials1, 16017 (2016)
work page 2016
-
[20]
Dai, Antiferromagnetic order and spin dynamics in iron-based superconductors, Rev
P. Dai, Antiferromagnetic order and spin dynamics in iron-based superconductors, Rev. Mod. Phys.87, 855 (2015)
work page 2015
- [21]
-
[22]
J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. Wu, G. F. Chen, X. Dai, J. L. Luo, Z. Fang, and N. L. Wang, Competing orders and spin-density-wave instability in La(O1xFx)FeAs, Europhysics Letters83, 27006 (2008)
work page 2008
-
[23]
Fawcett, Spin-density-wave antiferromagnetism in chromium, Rev
E. Fawcett, Spin-density-wave antiferromagnetism in chromium, Rev. Mod. Phys.60, 209 (1988)
work page 1988
-
[24]
D. Li, K. Lee, B. Y. Wang, M. Osada, S. Crossley, H. R. Lee, Y. Cui, Y. Hikita, and H. Y. Hwang, Supercon- ductivity in an infinite-layer nickelate, Nature572, 624 (2019)
work page 2019
-
[25]
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
-
[26]
H. Sakakibara, H. Usui, K. Suzuki, T. Kotani, H. Aoki, and K. Kuroki, Model construction and a possibility of cupratelike pairing in a newd 9 nickelate superconductor (Nd,Sr)nio 2, Phys. Rev. Lett.125, 077003 (2020)
work page 2020
- [27]
-
[28]
P. Adhikary, S. Bandyopadhyay, T. Das, I. Dasgupta, and T. Saha-Dasgupta, Orbital-selective superconductivity in a two-band model of infinite-layer nickelates, Phys. Rev. B102, 100501 (2020)
work page 2020
-
[29]
X. Wu, D. Di Sante, T. Schwemmer, W. Hanke, H. Y. Hwang, S. Raghu, and R. Thomale, Robustd x2−y2-wave superconductivity of infinite-layer nickelates, Phys. Rev. B101, 060504 (2020)
work page 2020
-
[30]
Z. Wang, G.-M. Zhang, Y.-f. Yang, and F.-C. Zhang, Dis- tinct pairing symmetries of superconductivity in infinite- layer nickelates, Phys. Rev. B102, 220501 (2020)
work page 2020
-
[31]
P. Werner and S. Hoshino, Nickelate superconductors: Multiorbital nature and spin freezing, Phys. Rev. B101, 041104 (2020)
work page 2020
-
[32]
Y.-H. Zhang and A. Vishwanath, Type-iit−jmodel in superconducting nickelate nd1−xsrxnio2, Phys. Rev. Res. 2, 023112 (2020)
work page 2020
-
[33]
M. H. Fischer, S. Wu, M. Lawler, A. Paramekanti, and E.-A. Kim, Nematic and spin-charge orders driven by hole-doping a charge-transfer insulator, New Journal of Physics16, 093057 (2014)
work page 2014
-
[34]
S.-H. Baek, D. V. Efremov, J. M. Ok, J. S. Kim, J. van den Brink, and B. B¨ uchner, Orbital-driven ne- maticity in FeSe, Nature Materials14, 210 (2015)
work page 2015
-
[35]
J. K. Glasbrenner, I. I. Mazin, H. O. Jeschke, P. J. Hirschfeld, R. M. Fernandes, and R. Valenti, Effect of magnetic frustration on nematicity and superconductiv- ity in iron chalcogenides, Nature Physics11, 953 (2015). 8
work page 2015
-
[37]
Q. Wang, Y. Shen, B. Pan, Y. Hao, M. Ma, F. Zhou, P. Steffens, K. Schmalzl, T. R. Forrest, M. Abdel-Hafiez, X. Chen, D. A. Chareev, A. N. Vasiliev, P. Bourges, Y. Sidis, H. Cao, and J. Zhao, Strong interplay between stripe spin fluctuations, nematicity and superconductiv- ity in FeSe, Nature Materials15, 159 (2016)
work page 2016
-
[38]
E. Berg, E. Fradkin, E.-A. Kim, S. A. Kivelson, V. Oganesyan, J. M. Tranquada, and S. C. Zhang, Dy- namical layer decoupling in a stripe-ordered high-T c su- perconductor, Phys. Rev. Lett.99, 127003 (2007)
work page 2007
-
[39]
D. F. Agterberg, J. C. S. Davis, S. D. Edkins, E. Fradkin, D. J. Van Harlingen, S. A. Kivelson, P. A. Lee, L. Radz- ihovsky, J. M. Tranquada, and Y. Wang, The Physics of Pair-Density Waves: Cuprate Superconductors and Be- yond, Annual Review of Condensed Matter Physics11, 231 (2020)
work page 2020
-
[40]
A. I. Larkin and Y. N. Ovchinnikov, Inhomogeneous state of superconductors, Soviet Physics JETP20, 762 (1965)
work page 1965
-
[41]
P. Fulde and R. A. Ferrell, Superconductivity in a strong spin-exchange field, Phys. Rev.135, A550 (1964)
work page 1964
-
[42]
E. Berg, E. Fradkin, and S. A. Kivelson, Pair-density- wave correlations in the kondo-heisenberg model, Phys. Rev. Lett.105, 146403 (2010)
work page 2010
-
[43]
P. Nikoli´ c, A. A. Burkov, and A. Paramekanti, Finite mo- mentum pairing instability of band insulators with mul- tiple bands, Phys. Rev. B81, 012504 (2010)
work page 2010
-
[44]
H.-C. Jiang and T. P. Devereaux, Pair density wave and superconductivity in a kinetically frustrated doped emery model on a square lattice (2023), arXiv:2309.11786 [cond- mat.str-el]
-
[45]
Y.-M. Wu, P. A. Nosov, A. A. Patel, and S. Raghu, Pair density wave order from electron repulsion, Phys. Rev. Lett.130, 026001 (2023)
work page 2023
- [46]
- [47]
-
[48]
R. Soto-Garrido and E. Fradkin, Pair-density-wave su- perconducting states and electronic liquid-crystal phases, Phys. Rev. B89, 165126 (2014)
work page 2014
-
[49]
H.-X. Wang, Y.-J. Hu, W. Huang, and H. Yao, A ”negative” route to pair density wave order (2025), arXiv:2512.06100 [cond-mat.supr-con]
-
[50]
G. Jiang and Y. Barlas, Pair density waves from local band geometry, Physical Review Letters131, 10.1103/physrevlett.131.016002 (2023)
- [51]
-
[52]
E. O. Lamponen, S. K. P¨ ontys, and P. T¨ orm¨ a, Supercon- ductivity and pair density waves from nearest-neighbor interactions in frustrated lattice geometries, Phys. Rev. B112, 144514 (2025)
work page 2025
- [53]
-
[54]
M. Yao, Y. Wang, D. Wang, J.-X. Yin, and Q.-H. Wang, Self-consistent theory of pair density waves in kagome superconductors, Physical Review B111, 10.1103/phys- revb.111.094505 (2025)
-
[55]
H. Oh and Y.-H. Zhang, Pair-density-wave supercon- ductivity and anderson’s theorem in bilayer nickelates (2025), arXiv:2512.15023 [cond-mat.str-el]
-
[56]
Y.-M. Wu, Z. Wu, and H. Yao, Pair-density-wave and chiral superconductivity in twisted bilayer transition metal dichalcogenides, Physical Review Letters130, 10.1103/physrevlett.130.126001 (2023)
-
[57]
K. Wang and K. Levin, Kekul´ e superconductiv- ity in twisted magic angle bilayer graphene (2026), arXiv:2510.06451 [cond-mat.supr-con]
-
[58]
S. A. Murshed and B. Roy, Nodal pair density waves from a quarter-metal in crystalline graphene multilayers, Physical Review B112, 10.1103/wy3f-hgr9 (2025)
- [59]
- [60]
-
[61]
A. Bose, S. Vadnais, and A. Paramekanti, Altermag- netism and superconductivity in a multiorbital tj model, Physical Review B110, 10.1103/physrevb.110.205120 (2024)
-
[62]
H.-Y. Liu, D. Wang, Z. Wang, and Q.-H. Wang, Genuine pair density wave order on the kagome lattice (2026), arXiv:2604.03531 [cond-mat.supr-con]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[63]
M. Takahashi, Half-filled Hubbard model at low temper- ature, Journal of Physics C: Solid State Physics10, 1289 (1977). 9 Appendix A: Weak Coupling
work page 1977
-
[64]
Charge/Spin Susceptibility We start by defining the generalized particle-hole bubble: χ0,αβγδ abcd (r, r′, τ) =T τ c† aα(r, τ)cbβ(r, τ)c† cγ(r′,0)c dδ(r′,0) c (A1) χ0,αβγδ abcd (Q, τ) = X k −1 N Gbc(k, τ)G da(k+Q,−τ)δ βγ δαδ (A2) χ0,αβγδ abcd (Q, iΩn) = −1 N β X k,nm X iωn [Ubn(k)U ∗ cn(k)Udm(k+Q)U ∗ am(k+Q)]G βγ n (k, iωn)Gαδ m (k+Q, iω n +iΩ n)δβγ δαδ (...
-
[65]
Singlet Pairing Susceptibility From the sign and overall form of the interactions considered, we derive the spin singlet orbital triplet pairing susceptibility. We define the spin singlet pairing operator and the pairing susceptibility: ∆†(Q) = 1√ 2 X k (X † k↑Y † −k+Q↓ −X † k↓Y † −k+Q↑) (A7) χp(Q, iΩn) = Z β 0 eiΩnτ dτ ∆(Q, τ)∆ †(Q,0) =1 2 Z β 0 eiΩnτ dτ...
-
[66]
RPA effective interactions The interactions we are interested in are the following: HJ =J X i ⃗SX · ⃗SY −V X i nX i nY i (A8) We decompose these interactions in the pairing, spin and charge channels. U c abcd = 0 0 0−V 0 0 0 0 0 0 0 0 −V0 0 0 (A9) U z abcd = 0 0 0J 0 0 0 0 0 0 0 0 J0 0 0 (A10) (A11) The multiorbital RPA equation reads:...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.