Dynamics of superconducting pairs in the two-dimensional Hubbard model
Pith reviewed 2026-05-17 02:21 UTC · model grok-4.3
The pith
In the two-dimensional Hubbard model, d-wave superconducting pairing is driven exclusively by low-frequency processes at the superexchange scale, not by high-frequency processes at the interaction strength U.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For all values of U and δ, pair-forming processes are confined to frequencies set by the superexchange interaction and followed by pair-breaking processes, ruling out both pair-forming and pair-breaking processes on the scale of U. This suggests that at high frequencies, the effect of U is eliminated by the d-wave pairing, and that at small frequencies, U generates the superexchange interaction that leads to low-frequency pair-forming processes providing the net contribution to pairing.
What carries the argument
The frequency structure of d-wave pairing correlations that separates pair-forming processes at superexchange frequencies from subsequent pair-breaking processes.
If this is right
- The net contribution to pairing arises only from low-frequency processes generated by superexchange.
- High-frequency effects from the on-site repulsion cancel due to the d-wave symmetry.
- The separation of forming and breaking processes is independent of specific values of U and doping δ.
- Pairing is mediated by magnetic superexchange rather than direct high-energy repulsion.
Where Pith is reading between the lines
- This result suggests that effective low-energy models focusing on superexchange could capture the essential pairing physics in related systems.
- Similar frequency separations might be observable in experimental spectra of cuprate materials or other strongly correlated models.
- Future work could test whether the pattern persists with longer-range interactions or in three-dimensional versions of the model.
Load-bearing premise
The cellular dynamical mean-field theory with the chosen cluster size and solver accurately captures the full frequency structure of d-wave pairing correlations in the thermodynamic limit.
What would settle it
An exact calculation on larger lattices that shows significant pair-forming or pair-breaking contributions at frequencies around the interaction strength U would disprove the claimed confinement to superexchange scales.
Figures
read the original abstract
The frequency structure of the superconducting correlations in cuprates gives insights on the pairing mechanism. Here we present an exhaustive study of this problem in the two-dimensional Hubbard model with cellular dynamical mean-field theory. To this end, we systematically quantify the dependence on doping $\delta$ and interaction strength $U$ of the superconducting gap, of the frequency scales where $d$-wave pairing occurs, and of their relative contribution to pairing. For all values of $U$ and $\delta$, we find pair-forming processes confined to frequencies set by the superexchange interaction and followed by pair-breaking processes, ruling out both pair-forming and pair-breaking processes on the scale of $U$. This suggests that at high frequencies, the effect of $U$ is eliminated by the $d$-wave paring, and that at small frequencies, $U$ generates the superexchange interaction that leads to low-frequency pair-forming processes providing the net contribution to pairing.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports an exhaustive cellular dynamical mean-field theory (cDMFT) study of the frequency structure of d-wave superconducting correlations in the two-dimensional Hubbard model. It systematically analyzes the dependence on interaction strength U and doping δ of the superconducting gap, the frequency scales of pair formation and breaking, and their relative contributions. The key result is that, for all U and δ, pair-forming processes occur at frequencies set by the superexchange J = 4t²/U, followed by pair-breaking processes, with no contributions on the bare U scale. This leads to the conclusion that d-wave pairing eliminates high-frequency U effects while low-frequency superexchange drives the net pairing.
Significance. If the central finding holds in the thermodynamic limit, this work would provide compelling numerical evidence that the pairing mechanism in the Hubbard model is dominated by superexchange interactions at low frequencies, with the on-site repulsion U playing no direct role in either pair formation or breaking at high frequencies due to d-wave symmetry. The comprehensive scan over parameters strengthens the generality of the claim. The use of cDMFT allows access to real-frequency information, which is a valuable contribution to the field.
major comments (1)
- The headline claim (abstract) that pair-forming and pair-breaking processes on the U scale are ruled out for all U and δ depends on the cDMFT frequency resolution of the pairing vertex or susceptibility. Because cDMFT replaces the lattice self-energy with a finite-cluster (typically 2×2) impurity problem, any momentum-dependent charge fluctuations that could feed U-scale weight into the low-frequency channel are filtered by construction. Without an explicit cluster-size extrapolation or cross-check against a method retaining full momentum resolution (e.g., DQMC), the observed clean J-scale / U-scale separation may be an artifact of the truncation rather than a property of the thermodynamic-limit model.
minor comments (2)
- Clarify in the methods section the precise cluster size, impurity solver, and post-processing procedure used to extract the frequency-dependent pairing quantities; this is needed to assess the claimed resolution at both J and U scales.
- In figures showing frequency dependence, label the characteristic scales J = 4t²/U and U explicitly on the axes to facilitate direct visual comparison across different parameter values.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive report. The major comment raises an important point about the possible influence of the finite cluster size in our cDMFT calculations. We address this below and will revise the manuscript accordingly to strengthen the discussion of methodological limitations.
read point-by-point responses
-
Referee: The headline claim (abstract) that pair-forming and pair-breaking processes on the U scale are ruled out for all U and δ depends on the cDMFT frequency resolution of the pairing vertex or susceptibility. Because cDMFT replaces the lattice self-energy with a finite-cluster (typically 2×2) impurity problem, any momentum-dependent charge fluctuations that could feed U-scale weight into the low-frequency channel are filtered by construction. Without an explicit cluster-size extrapolation or cross-check against a method retaining full momentum resolution (e.g., DQMC), the observed clean J-scale / U-scale separation may be an artifact of the truncation rather than a property of the thermodynamic-limit model.
Authors: We agree that the 2×2 cDMFT setup approximates the lattice and inherently limits the treatment of long-range momentum-dependent fluctuations that could, in principle, mix high-frequency U-scale weight into lower frequencies. However, the d-wave pairing symmetry itself provides a robust mechanism for suppressing on-site U contributions at high frequencies, since the pair wavefunction changes sign between neighboring sites and vanishes on-site; this symmetry argument is captured already at the cluster level. The superexchange scale J emerges from virtual hopping processes within the cluster and consistently sets the pair-forming window across our broad parameter scan in U and δ. While we have not performed a systematic cluster-size extrapolation in this work, the uniformity of the J/U separation for all studied parameters makes an artifact less likely. We will revise the manuscript to include an expanded discussion of these limitations, referencing supporting evidence from DQMC and larger-cluster studies in the literature that report compatible low-frequency pairing dominance. revision: partial
Circularity Check
No circularity: direct numerical output from cDMFT equations
full rationale
The paper reports frequency scales of d-wave pairing correlations obtained by solving the cellular DMFT equations on finite clusters for the 2D Hubbard model. The central claim—that pair-forming weight is confined to frequencies ~J=4t²/U while pair-breaking occurs at higher scales, with no U-scale processes—is an output of the impurity solver and cluster self-consistency loop rather than a quantity fitted or defined to reproduce a target result. No self-definitional loops, fitted inputs renamed as predictions, or load-bearing self-citations that reduce the result to prior unverified claims by the same authors appear in the derivation chain. The work is self-contained against external benchmarks in the sense that the reported spectra follow from the stated approximation and can be cross-checked by independent implementations of the same method.
Axiom & Free-Parameter Ledger
free parameters (2)
- interaction strength U
- doping δ
axioms (2)
- domain assumption Cellular DMFT with finite cluster provides a controlled approximation to the 2D Hubbard model dynamics
- standard math Superexchange J is generated by virtual processes at scale t²/U
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
For all values of U and δ, we find pair-forming processes confined to frequencies set by the superexchange interaction and followed by pair-breaking processes, ruling out both pair-forming and pair-breaking processes on the scale of U.
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the long-lived pair-forming processes occurring in this frequency interval are associated with short-range spin fluctuations... This frequency scale is of the order of J/2
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
M. R. Norman, The challenge of unconventional super- conductivity, Science332, 196 (2011)
work page 2011
-
[2]
D. J. Scalapino, A common thread: The pairing inter- action for unconventional superconductors, Rev. Mod. Phys.84, 1383 (2012)
work page 2012
- [3]
-
[4]
E. van Heumen, E. Muhlethaler, A. B. Kuzmenko, H. Eisaki, W. Meevasana, M. Greven, and D. van der 6 Marel, Optical determination of the relation between the electron-boson coupling function and the critical tem- perature in high-T c cuprates, Phys. Rev. B79, 184512 (2009)
work page 2009
-
[5]
J. P. Carbotte, T. Timusk, and J. Hwang, Bosons in high- temperature superconductors: an experimental survey, Reports on Progress in Physics74, 066501 (2011)
work page 2011
-
[6]
S. D. Conte, C. Giannetti, G. Coslovich, F. Cilento, D. Bossini, T. Abebaw, F. Banfi, G. Ferrini, H. Eisaki, M. Greven, A. Damascelli, D. van der Marel, and F. Parmigiani, Disentangling the Electronic and Phononic Glue in a High-T c Superconductor, Science 335, 1600 (2012)
work page 2012
-
[7]
F. Cilento, S. D. Conte, G. Coslovich, F. Banfi, G. Fer- rini, H. Eisaki, M. Greven, A. Damascelli, D. v. d. Marel, F. Parmigiani, and C. Giannetti, In search for the pairing glue in cuprates by non-equilibrium optical spectroscopy, Journal of Physics: Conference Series449, 012003 (2013)
work page 2013
-
[8]
S. Dal Conte, L. Vidmar, D. Goleˇ z, M. Mierzejewski, G. Soavi, S. Peli, F. Banfi, G. Ferrini, R. Comin, B. M. Ludbrook, L. Chauviere, N. D. Zhigadlo, H. Eisaki, M. Greven, S. Lupi, A. Damascelli, D. Brida, M. Capone, J. Bonˇ ca, G. Cerullo, and C. Giannetti, Snapshots of the retarded interaction of charge carriers with ultrafast fluc- tuations in cuprate...
work page 2015
-
[9]
P. W. Anderson, Is there glue in cuprate superconduc- tors?, Science316, 1705 (2007)
work page 2007
- [10]
- [11]
-
[12]
S. Acharya, C. Weber, E. Plekhanov, D. Pashov, A. Tara- phder, and M. Van Schilfgaarde, Metal-Insulator Transi- tion in Copper Oxides Induced by Apex Displacements, Phys. Rev. X8, 021038 (2018)
work page 2018
-
[13]
B. Bacq-Labreuil, B. Lacasse, A.-M. S. Tremblay, D. S´ en´ echal, and K. Haule, Toward an Ab Initio The- ory of High-Temperature Superconductors: A Study of Multilayer Cuprates, Phys. Rev. X15, 021071 (2025)
work page 2025
-
[14]
Z.-H. Cui, J. Yang, J. T¨ olle, H.-Z. Ye, S. Yuan, H. Zhai, G. Park, R. Kim, X. Zhang, L. Lin, T. C. Berkelbach, and G. K.-L. Chan, Ab initio quantum many-body descrip- tion of superconducting trends in the cuprates, Nature Communications16, 1845 (2025)
work page 2025
- [15]
-
[16]
C. Stahl and M. Eckstein, Noise correlations in time- and angle-resolved photoemission spectroscopy, Phys. Rev. B 99, 241111 (2019)
work page 2019
-
[17]
T. P. Devereaux, M. Claassen, X.-X. Huang, M. Zaletel, J. E. Moore, D. Morr, F. Mahmood, P. Abbamonte, and Z.-X. Shen, Angle-resolved pair photoemission theory for correlated electrons, Phys. Rev. B108, 165134 (2023)
work page 2023
-
[18]
F. Boschini, M. Zonno, and A. Damascelli, Time-resolved arpes studies of quantum materials, Rev. Mod. Phys.96, 015003 (2024)
work page 2024
-
[19]
P. W. Anderson, The resonating valence bond state in La 2CuO4 and superconductivity, Science235, 1196 (1987)
work page 1987
-
[20]
A.-M. S. Tremblay, Strongly correlated superconductiv- ity, inEmergent Phenomena in Correlated Matter Mod- eling and Simulation, Vol. 3, edited by E. Pavarini, E. Koch, and U. Schollw¨ ock (Verlag des Forschungszen- trum, J¨ ulich, 2013) Chap. 10
work page 2013
-
[21]
G. Kotliar and J. Liu, Superexchange mechanism and d-wave superconductivity, Phys. Rev. B38, 5142 (1988)
work page 1988
- [22]
-
[23]
A.-M. S. Tremblay, B. Kyung, and D. S´ en´ echal, Pseudo- gap and high-temperature superconductivity from weak to strong coupling. Towards a quantitative theory, Low Temp. Phys.32, 424 (2006)
work page 2006
-
[24]
M. Qin, T. Sch¨ afer, S. Andergassen, P. Corboz, and E. Gull, The Hubbard Model: A Computational Per- spective, Annual Review of Condensed Matter Physics 13, 275 (2022)
work page 2022
-
[25]
G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. Parcollet, and C. A. Marianetti, Electronic struc- ture calculations with dynamical mean-field theory, Rev. Mod. Phys.78, 865 (2006)
work page 2006
-
[26]
A. Georges, G. Kotliar, W. Krauth, and M. J. Rozen- berg, Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions, Rev. Mod. Phys.68, 13 (1996)
work page 1996
-
[27]
T. A. Maier, D. Poilblanc, and D. J. Scalapino, Dynamics of the Pairing Interaction in the Hubbard andt−JMod- els of High-Temperature Superconductors, Phys. Rev. Lett.100, 237001 (2008)
work page 2008
- [28]
-
[29]
M. Civelli, Evolution of the dynamical pairing across the phase diagram of a strongly correlated high-temperature superconductor, Phys. Rev. Lett.103, 136402 (2009)
work page 2009
-
[30]
D. S´ en´ echal, A. G. R. Day, V. Bouliane, and A.- M. S. Tremblay, Resilience ofd-wave superconductivity to nearest-neighbor repulsion, Phys. Rev. B87, 075123 (2013)
work page 2013
-
[31]
E. Gull and A. J. Millis, Pairing glue in the two- dimensional Hubbard model, Phys. Rev. B90, 041110 (2014)
work page 2014
-
[32]
A. Reymbaut, M. Charlebois, M. F. Asiani, L. Fratino, P. S´ emon, G. Sordi, and A.-M. S. Tremblay, Antagonis- tic effects of nearest-neighbor repulsion on the supercon- ducting pairing dynamics in the doped Mott insulator regime, Phys. Rev. B94, 155146 (2016)
work page 2016
-
[33]
X. Dong, L. Del Re, A. Toschi, and E. Gull, Mechanism of superconductivity in the Hubbard model at interme- diate interaction strength, Proceedings of the National Academy of Science119, e2205048119 (2022)
work page 2022
-
[34]
X. Dong, E. Gull, and A. J. Millis, Quantifying the role of antiferromagnetic fluctuations in the superconductivity of the doped Hubbard model, Nature Physics18, 1293 (2022)
work page 2022
-
[35]
K. Haule, Quantum Monte Carlo impurity solver for clus- ter dynamical mean-field theory and electronic structure calculations with adjustable cluster base, Phys. Rev. B 75, 155113 (2007)
work page 2007
- [36]
-
[37]
E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. Troyer, and P. Werner, Continuous-time Monte Carlo methods for quantum impurity models, Rev. Mod. Phys. 83, 349 (2011)
work page 2011
-
[38]
P. S´ emon, C.-H. Yee, K. Haule, and A.-M. S. Trem- blay, Lazy skip-lists: An algorithm for fast hybridization- expansion quantum Monte Carlo, Phys. Rev. B90, 075149 (2014)
work page 2014
-
[39]
P. S´ emon, G. Sordi, and A.-M. S. Tremblay, Ergodicity of the hybridization-expansion Monte Carlo algorithm for broken-symmetry states, Phys. Rev. B89, 165113 (2014)
work page 2014
-
[40]
D. Bergeron and A.-M. S. Tremblay, Algorithms for opti- mized maximum entropy and diagnostic tools for analytic continuation, Phys. Rev. E94, 023303 (2016)
work page 2016
-
[41]
A. Reymbaut, D. Bergeron, and A.-M. S. Tremblay, Max- imum entropy analytic continuation for spectral func- tions with nonpositive spectral weight, Phys. Rev. B92, 060509 (2015)
work page 2015
- [42]
-
[43]
K. Haule and G. Kotliar, Strongly correlated super- conductivity: A plaquette dynamical mean-field theory study, Phys. Rev. B76, 104509 (2007)
work page 2007
- [44]
-
[45]
L. Fratino, P. S´ emon, G. Sordi, and A.-M. S. Tremblay, An organizing principle for two-dimensional strongly cor- related superconductivity, Sci. Rep.6, 22715 (2016)
work page 2016
-
[46]
This neglects Kosterlitz-Thouless physics
- [47]
-
[48]
A. Paramekanti, M. Randeria, and N. Trivedi, High-T c superconductors: A variational theory of the supercon- ducting state, Phys. Rev. B70, 054504 (2004)
work page 2004
-
[49]
S. S. Kancharla, B. Kyung, D. S´ en´ echal, M. Civelli, M. Capone, G. Kotliar, and A.-M. S. Tremblay, Anoma- lous superconductivity and its competition with antifer- romagnetism in doped Mott insulators, Phys. Rev. B77, 184516 (2008)
work page 2008
-
[50]
M. Civelli, M. Capone, A. Georges, K. Haule, O. Parcol- let, T. D. Stanescu, and G. Kotliar, Nodal-Antinodal Di- chotomy and the Two Gaps of a Superconducting Doped Mott Insulator, Phys. Rev. Lett.100, 046402 (2008)
work page 2008
-
[51]
E. Gull, O. Parcollet, and A. J. Millis, Superconductiv- ity and the Pseudogap in the Two-Dimensional Hubbard Model, Phys. Rev. Lett.110, 216405 (2013)
work page 2013
- [52]
-
[53]
J. A. Sobota, Y. He, and Z.-X. Shen, Angle-resolved photoemission studies of quantum materials, Rev. Mod. Phys.93, 025006 (2021)
work page 2021
- [54]
- [55]
-
[56]
G. Sordi, P. S´ emon, K. Haule, and A.-M. S. Tremblay, Pseudogap temperature as a Widom line in doped Mott insulators, Sci. Rep.2, 547 (2012). Appendix Figure 5 shows a zoom in of Fig. 3 on the low frequencies. 0.00 0.05 δ 0.0 0.1 0.2 0.3 0.4 0.5 ω (a) U = 5.2 0.00 0.05 0.10 δ 0.0 0.1 0.2 0.3 0.4 0.5 ω (b) U = 6.2 0.00 0.05 0.10 0.15 δ 0.0 0.1 0.2 0.3 0.4...
work page 2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.