Extended Hubbard model on fractals: d-Wave superconductivity and competing pairing channels
Pith reviewed 2026-05-15 09:06 UTC · model grok-4.3
The pith
The Sierpiński carpet fractal lattice suppresses the d-wave superconducting dome at half filling while strongly enhancing extended s-wave pairing at other fillings.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using Bogoliubov-de Gennes mean-field theory on the extended Hubbard model with nearest-neighbor attraction, the authors find that the Sierpiński carpet dramatically alters the competition between pairing channels: the predominant d-wave superconducting dome at half filling of the square lattice becomes unstable for the carpet, while at high and low fillings extended s-wave pairing gets strongly enhanced. They attribute this to geometric frustration of sign-changing order parameters by the fractal boundary structure. On the triangular Sierpiński gasket, hybrid s+d+id states show critical temperature enhancement comparable to that previously observed for pure s-wave pairing.
What carries the argument
Geometric frustration of sign-changing superconducting order parameters imposed by the irregular boundary structure of the Sierpiński carpet, within the Bogoliubov-de Gennes mean-field treatment of the extended Hubbard model.
If this is right
- d-wave superconductivity is eliminated at half filling on the carpet while extended s-wave pairing is stabilized at high and low fillings.
- Hybrid s+d+id states on the gasket reach critical temperatures comparable to those of pure s-wave pairing.
- Fractal geometry selectively filters pairing symmetries according to the compatibility of their nodal structure with the lattice boundaries.
- The same geometric mechanism that suppresses d-wave order can be used to tune the dominant channel by changing electron filling.
Where Pith is reading between the lines
- Fractal boundary engineering could be explored as a route to stabilize s-wave or hybrid pairing in artificial nanostructures where d-wave order is normally preferred.
- The selective-filter effect may extend to other fractal geometries or to models with longer-range interactions, offering a design principle for controlling superconducting symmetry.
- Quantitative comparison of mean-field Tc values with fluctuation-corrected methods on the same fractals would test how robust the reported channel competition remains beyond the approximation.
Load-bearing premise
The Bogoliubov-de Gennes mean-field approximation remains quantitatively reliable for capturing pairing competition and critical temperatures on these irregular fractal lattices without significant corrections from fluctuations or exact many-body effects.
What would settle it
Quantum Monte Carlo or exact diagonalization calculations on finite Sierpiński-carpet clusters that find a stable d-wave solution persisting at half filling or that yield lower critical temperatures for extended s-wave pairing than the mean-field prediction.
Figures
read the original abstract
Fractal structures such as the Sierpi\'nski gasket have been predicted to enhance the critical temperature of s-wave superconductivity compared to regular crystals while maintaining macroscopic phase coherence of Cooper pairs. Here we extend this analysis to order parameters with non-trivial symmetry by studying the extended Hubbard model with nearest-neighbor attraction on fractal lattices. Using Bogoliubov-de Gennes mean-field theory, we find that the Sierpi\'nski carpet dramatically alters the competition between pairing channels: the predominant d-wave superconducting dome at half filling of the square lattice becomes unstable for the carpet, while at high and low fillings extended s-wave pairing gets strongly enhanced. We attribute this to geometric frustration of sign-changing order parameters by the fractal boundary structure. On the triangular Sierpi\'nski gasket, hybrid s+d+id states show critical temperature enhancement comparable to that previously observed for pure s-wave pairing. Our results demonstrate that fractal geometry acts as a selective filter for pairing symmetries, with the compatibility between order parameter structure and lattice topology determining which channels are stabilized or suppressed.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript applies Bogoliubov-de Gennes mean-field theory to the extended Hubbard model with nearest-neighbor attraction on Sierpiński carpet and gasket fractals. It reports that the carpet suppresses the d-wave superconducting dome present at half-filling on the square lattice, while strongly enhancing extended s-wave pairing at high and low fillings; this is attributed to geometric frustration of sign-changing order parameters by the fractal boundaries. On the triangular gasket, hybrid s+d+id states exhibit critical-temperature enhancement comparable to prior s-wave results.
Significance. If the mean-field results prove robust, the work demonstrates that fractal geometry can act as a selective filter for pairing symmetries, stabilizing or suppressing channels according to compatibility with lattice topology. This extends earlier s-wave studies on fractals and suggests routes to engineer specific order parameters or higher Tc via boundary-induced frustration.
major comments (2)
- [Methods and Results sections on BdG implementation and carpet at half filling] The central claims on d-wave instability at half-filling and the crossover to extended s-wave rest entirely on self-consistent BdG solutions. No comparisons to fluctuation-corrected approaches (RPA, cluster DMFT) or exact methods on small clusters are provided, despite the carpet having sites with coordination numbers 2–4 where local fluctuations are expected to be large and could stabilize sign-changing order.
- [Results paragraph discussing channel competition on the carpet] The attribution of d-wave suppression to 'geometric frustration by the fractal boundary structure' lacks supporting diagnostics such as site-resolved order-parameter signs, phase windings at corners, or a control calculation on a non-fractal lattice with comparable boundaries.
minor comments (2)
- The abstract and methods should explicitly state the ranges of Hubbard parameters U and V, the fractal iteration depths, and the convergence criteria (energy tolerance, iteration count) used for the self-consistent solutions.
- Figure captions for Tc domes and order-parameter maps should clarify color scales, whether values are averaged or site-resolved, and any smoothing applied.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and for the constructive comments, which help clarify the scope and limitations of our mean-field study. We address each major comment point by point below.
read point-by-point responses
-
Referee: The central claims on d-wave instability at half-filling and the crossover to extended s-wave rest entirely on self-consistent BdG solutions. No comparisons to fluctuation-corrected approaches (RPA, cluster DMFT) or exact methods on small clusters are provided, despite the carpet having sites with coordination numbers 2–4 where local fluctuations are expected to be large and could stabilize sign-changing order.
Authors: We agree that the results rely on self-consistent Bogoliubov-de Gennes mean-field theory. This approach is standard for mapping out pairing instabilities and their competition on complex, non-periodic lattices such as fractals, where the irregular connectivity makes cluster DMFT or RPA implementations computationally prohibitive at the system sizes required to capture the fractal boundaries. While local fluctuations at low-coordination sites (coordination 2–4) are indeed expected to be stronger than on the square lattice, the geometric frustration mechanism we identify arises from the global topology and boundary structure rather than local coordination alone; mean-field therefore still captures the qualitative suppression of sign-changing order. In the revised manuscript we have added a dedicated paragraph in the Methods and Discussion sections explicitly discussing these limitations of the mean-field approximation and outlining how exact diagonalization on small clusters or future fluctuation-corrected calculations could provide further validation. revision: partial
-
Referee: The attribution of d-wave suppression to 'geometric frustration by the fractal boundary structure' lacks supporting diagnostics such as site-resolved order-parameter signs, phase windings at corners, or a control calculation on a non-fractal lattice with comparable boundaries.
Authors: We thank the referee for this helpful suggestion. In the revised Results section we now include site-resolved maps of the superconducting order parameter for both d-wave and extended s-wave channels on the carpet at half filling and at representative dopings. These maps explicitly show the sign changes of the d-wave order parameter being disrupted at the fractal corners and edges, while the extended s-wave order parameter remains uniform. We have also added a control calculation on a square lattice with artificially introduced boundaries that mimic the perimeter of the carpet; this non-fractal control does not exhibit the same d-wave suppression, thereby supporting the role of the self-similar fractal structure in generating the geometric frustration. revision: yes
Circularity Check
No circularity: direct numerical BdG solution on fractals
full rationale
The paper applies standard Bogoliubov-de Gennes mean-field theory to the extended Hubbard model with nearest-neighbor attraction on explicitly constructed fractal lattices (Sierpiński carpet and gasket). Pairing-channel competition and critical temperatures are obtained by direct numerical solution of the self-consistent gap equations on the site-inhomogeneous graphs, with order-parameter symmetries (d-wave, extended s-wave, hybrid s+d+id) defined independently from the lattice geometry. No parameters are fitted to data subsets and then relabeled as predictions; no load-bearing claims reduce to self-citations or prior ansatzes by the same authors; the geometric-frustration interpretation follows from the explicit lattice structure and sign-changing properties of the order parameters. The derivation chain is therefore self-contained in the numerical implementation.
Axiom & Free-Parameter Ledger
free parameters (1)
- Hubbard interaction parameters
axioms (1)
- domain assumption Bogoliubov-de Gennes mean-field theory accurately describes the superconducting order parameter competition on fractal lattices
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Using Bogoliubov-de Gennes mean-field theory, we find that the Sierpiński carpet dramatically alters the competition between pairing channels: the predominant d-wave superconducting dome at half filling of the square lattice becomes unstable for the carpet, while at high and low fillings extended s-wave pairing gets strongly enhanced. We attribute this to geometric frustration of sign-changing order parameters by the fractal boundary structure.
-
IndisputableMonolith/Foundation/BranchSelection.leanbranch_selection unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The extended Hubbard model with both on-site (U) and nearest-neighbor (V) attractive interactions
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]
S. N. Kempkes, M. R. Slot, S. E. Freeney, S. J. Zeven- huizen, D. Vanmaekelbergh, I. Swart, and C. M. Smith, Nature Physics15, 127 (2019)
work page 2019
- [2]
-
[3]
C. Li, X. Zhang, N. Li, Y. Wang, J. Yang, G. Gu, Y. Zhang, S. Hou, L. Peng, K. Wu, and D. Nieckarz, Journal of the American Chemical Society139, 13749 (2017)
work page 2017
-
[4]
F. De Nicola, N. S. P. Purayil, D. Spirito, M. Mis- cuglio, F. Tantussi, A. Tomadin, F. De Angelis, M. Polini, R. Krahne, and V. Pellegrini, ACS Photonics5, 2418 (2018)
work page 2018
- [5]
- [6]
- [7]
-
[8]
A. A. Iliasov, M. I. Katsnelson, and S. Yuan, Physical Review B99, 075402 (2019)
work page 2019
-
[9]
A. A. Iliasov, M. I. Katsnelson, and S. Yuan, Physical Review B102, 075440 (2020)
work page 2020
-
[10]
Q. Yao, X. Yang, A. A. Iliasov, M. I. Katsnelson, and S. Yuan, Physical Review B107, 115424 (2023)
work page 2023
- [11]
-
[12]
Laskin, Physical Review E66, 056108 (2002)
N. Laskin, Physical Review E66, 056108 (2002)
work page 2002
- [13]
-
[14]
E. van Veen, A. Tomadin, M. Polini, M. I. Katsnelson, and S. Yuan, Physical Review B96, 235438 (2017)
work page 2017
-
[15]
F. Salvati, M. I. Katsnelson, and A. A. Bagrov, Com- munications Physics9, 24 (2026)
work page 2026
-
[16]
K. Wang, Y. Liu, and T. Liang, Physica B: Condensed Matter498, 33 (2016)
work page 2016
- [17]
- [18]
-
[19]
M. Brzezi´ nska, A. M. Cook, and T. Neupert, Physical Review B98, 205116 (2018)
work page 2018
- [20]
-
[21]
R. Canyellas, C. Liu, R. Arouca, L. Eek, G. Wang, Y. Yin, D. Guan, Y. Li, S. Wang, H. Zheng, C. Liu, J. Jia, and C. Morais Smith, Nature Physics20, 1421 (2024)
work page 2024
-
[22]
M. N. Ivaki, I. Sahlberg, K. P¨ oyh¨ onen, and T. Ojanen, Communications Physics5, 327 (2022)
work page 2022
-
[23]
X. Yang, W. Zhou, Q. Yao, P. Lv, Y. Wang, and S. Yuan, Physical Review B105, 205433 (2022)
work page 2022
-
[24]
T. Westerhout, E. van Veen, M. I. Katsnelson, and S. Yuan, Physical Review B97, 205434 (2018)
work page 2018
-
[25]
M. St˚ alhammar and C. M. Smith, Physical Review Re- search5, 043043 (2023)
work page 2023
-
[26]
M. Ippoliti, T. Rakovszky, and V. Khemani, Physical Review X12, 011045 (2022)
work page 2022
-
[27]
D. S. Ageev, A. A. Bagrov, and A. A. Iliasov, Physical Review B103, L100302 (2021)
work page 2021
- [28]
- [29]
- [30]
- [31]
-
[32]
K. Dutkiewicz, M. P lodzie´ n, A. Rojo-Franc` as, B. Juli´ a- 11 D´ ıaz, M. Lewenstein, and T. Grass, Phys. Rev. A113, 013320 (2026)
work page 2026
-
[33]
Strong enhancement of superconduc- tivity on finitely ramified fractal lattices,
A. A. Iliasov, R. Canyellas, M. I. Katsnelson, and A. A. Bagrov, “Strong enhancement of superconduc- tivity on finitely ramified fractal lattices,” (2024), arXiv:2310.11497, arXiv:2310.11497 [cond-mat.supr- con]
-
[34]
N. N. Bogoliubov, Nuovo Cimento7, 794 (1958)
work page 1958
-
[35]
P. G. de Gennes,Superconductivity of Metals and Alloys (CRC Press, 2018)
work page 2018
-
[36]
Zhu,Bogoliubov-de Gennes Method and Its Appli- cations, Lecture Notes in Physics, Vol
J.-X. Zhu,Bogoliubov-de Gennes Method and Its Appli- cations, Lecture Notes in Physics, Vol. 924 (Springer, 2016)
work page 2016
- [37]
-
[38]
S. A. Sousa-J´ unior, N. C. Costa, R. R. dos Santos, and T. Paiva, Physical Review B109, 165102 (2024)
work page 2024
-
[39]
D. J. Scalapino, S. R. White, and S. Zhang, Physical Review B47, 7995 (1993)
work page 1993
-
[40]
Lecture 6.5 — RMSProp: Divide the gradient by a running average of its recent magnitude,
T. Tieleman and G. Hinton, “Lecture 6.5 — RMSProp: Divide the gradient by a running average of its recent magnitude,” COURSERA: Neural Networks for Machine Learning (2012)
work page 2012
- [41]
- [42]
-
[43]
E. G. C. P. van Loon and M. I. Katsnelson, Journal of Physics: Conference Series1136, 012006 (2018)
work page 2018
-
[44]
M. Harland, S. Brener, A. I. Lichtenstein, and M. I. Katsnelson, Physical Review B100, 024510 (2019)
work page 2019
- [45]
- [46]
-
[47]
E. G. C. P. van Loon, M. R¨ osner, M. I. Katsnelson, and T. O. Wehling, Physical Review B104, 045134 (2021)
work page 2021
-
[48]
A. Chen, Z.-Q. Wan, A. Sengupta, A. Georges, and C. Roth, arXiv preprint arXiv:2507.10705 (2025), 10.48550/arXiv.2507.10705, arXiv:2507.10705 [cond- mat.str-el]
-
[49]
V. Bashmakov, A. Iliasov, T. Bzduˇ sek, and A. A. Bagrov, arXiv preprint arXiv:2509.09330 (2025), arXiv:2509.09330 [cond-mat.supr-con]
- [50]
-
[51]
Y. Zhang, X. Zhang, Y. Li, S. Zhao, S. Hou, K. Wu, and Y. Wang, Journal of the American Chemical Society 142, 17928 (2020). Appendix A: Local channel mixing induced by bond removal In this Appendix we illustrate, at the level of a single coordination shell, how bond removal modifies the struc- ture of nearest-neighbor pairing channels. The purpose is pure...
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.