pith. machine review for the scientific record. sign in

arxiv: 2601.11055 · v2 · submitted 2026-01-16 · ❄️ cond-mat.str-el · cond-mat.supr-con

Recognition: 2 theorem links

· Lean Theorem

Surface Functional Renormalization Group for Layered Quantum Materials

Authors on Pith no claims yet

Pith reviewed 2026-05-16 13:57 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.supr-con
keywords functional renormalization groupHubbard modelinterlayer couplingd-wave superconductivityspin-density wavespin-bond ordersurface stateslayered materials
0
0 comments X

The pith

Extending functional renormalization group to surfaces splits the superconducting phase with a narrow incommensurate magnetic window at moderate interlayer hopping.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a surface-adapted version of the functional renormalization group that incorporates interlayer hopping into the flow equations for a semi-infinite stack of square lattices. It applies this method to a Hubbard model confined to the top layer with alternating interlayer couplings. Across most of the phase diagram the surface retains the familiar two-dimensional Hubbard physics, dominated by antiferromagnetic, d-wave superconducting, and ferromagnetic correlations. At intermediate interlayer hopping strengths and moderate interaction values, however, the superconducting regime splits into two separate pockets separated by a narrow strip of incommensurate spin-density-wave and spin-bond order.

Core claim

The central claim is that an SSH-like interlayer coupling added to the surface Hubbard model leaves the two-dimensional physics intact for wide parameter ranges but, at intermediate interlayer hopping, inserts a small region of incommensurate spin-density-wave and spin-bond order that divides the superconducting state into two regimes, thereby opening a possible route to chiral spin-bond order.

What carries the argument

The surface functional renormalization group flow, obtained by extending the standard two-dimensional fRG equations to include interlayer hopping terms while keeping the interaction vertex on the surface layer.

If this is right

  • Antiferromagnetic, d-wave superconducting, and ferromagnetic correlations remain the dominant instabilities for most values of interlayer hopping.
  • The superconducting state at intermediate interaction strengths is interrupted only in a limited window of interlayer couplings.
  • Incommensurate spin-density-wave and spin-bond order appear together in that narrow window.
  • The spin-bond order inside the window can in principle take a chiral form.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could be used to study interfaces in other van-der-Waals heterostructures where only the top layer carries strong correlations.
  • The incommensurate window might be located experimentally by tuning interlayer distance in a layered material or by applying pressure.
  • If the chiral spin-bond order is realized, it would break time-reversal symmetry and could produce detectable edge currents or anomalous Hall signals.

Load-bearing premise

The truncation and discretization used in the functional renormalization group flow remain accurate once interlayer hopping is added to the surface Hubbard model.

What would settle it

A calculation with finer momentum discretization or a higher-order truncation that either eliminates the narrow incommensurate region or moves its location substantially would falsify the reported splitting of the superconducting dome.

Figures

Figures reproduced from arXiv: 2601.11055 by Dante M. Kennes, Lennart Klebl.

Figure 1
Figure 1. Figure 1: FIG. 1. Three-dimensional Hubbard-SSH model. Each of the [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Diagrammatic representation of the flow equation of [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Critical scale Λ [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Magnetic susceptibilities [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Presence of the iSDW state (red) and the [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
read the original abstract

We present an extension to the two-dimensional functional renormalization group to efficiently treat interactions on the surface or at interfaces of three-dimensional systems. As an application, we consider a semi-infinite stack of two-dimensional square lattices, including a Hubbard interaction on the surface layer and an alternating interlayer coupling. We investigate how strongly correlated states of the decoupled two-dimensional Hubbard model on the surface evolve under inclusion of such an SSH-like interlayer coupling. For large parts of the phase diagram as a function of the interlayer hopping parameters, the physics of the two-dimensional system prevails, with antiferromagnetic, superconducting $d$-wave, and ferromagnetic correlations taking center stage. However, for intermediate interlayer couplings the superconducting state at intermediate interaction strengths separates into two regimes by a small region of incommensurate spin-density-wave and spin-bond order, enabling the potential realization of chiral spin-bond order.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper introduces an extension of the two-dimensional functional renormalization group (fRG) to surfaces and interfaces of three-dimensional layered systems. Applied to a semi-infinite stack of square lattices with on-site Hubbard repulsion confined to the surface layer and an alternating (SSH-like) interlayer hopping, the work maps the evolution of the 2D Hubbard-model phases (antiferromagnetism, d-wave superconductivity, ferromagnetism) as interlayer coupling is turned on. The central claim is that, for intermediate interlayer hoppings and interaction strengths, the d-wave superconducting regime is interrupted by a narrow window of incommensurate spin-density-wave plus spin-bond order, potentially allowing chiral spin-bond order.

Significance. If the reported narrow incommensurate window is robust, the result would provide a concrete route to stabilize chiral spin-bond order at surfaces of layered materials and would demonstrate how weak interlayer coupling can qualitatively reorganize the 2D Hubbard phase diagram. The methodological extension itself is a useful technical step for treating surface correlations in quasi-2D systems.

major comments (2)
  1. [Numerical implementation and results section] The headline result—a small incommensurate SDW/spin-bond region that splits the d-wave superconducting dome for intermediate t_perp—rests on the accuracy of the momentum and frequency discretization of the extended fRG flow. No convergence tests with respect to patch number, grid density, or frequency mesh are reported once finite interlayer hopping is introduced, even though the location of incommensurate instabilities is known to be sensitive to nesting-vector resolution and to the opening of new interlayer scattering channels.
  2. [Method and truncation discussion] The truncation to the two-particle vertex (standard in fRG) is retained without additional justification or error estimates when the interlayer term couples momenta along the stacking direction. Because the claimed window is narrow, it is essential to demonstrate that the truncation does not artificially stabilize or shift the incommensurate instabilities relative to the t_perp = 0 limit.
minor comments (1)
  1. [Abstract and introduction] The abstract and introduction would benefit from a brief statement of the discretization parameters (number of patches, frequency points) used for the surface fRG flow.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We appreciate the emphasis on numerical convergence and the validity of the truncation scheme. Below we address the major comments point by point. We will incorporate revisions to strengthen the presentation of our results.

read point-by-point responses
  1. Referee: [Numerical implementation and results section] The headline result—a small incommensurate SDW/spin-bond region that splits the d-wave superconducting dome for intermediate t_perp—rests on the accuracy of the momentum and frequency discretization of the extended fRG flow. No convergence tests with respect to patch number, grid density, or frequency mesh are reported once finite interlayer hopping is introduced, even though the location of incommensurate instabilities is known to be sensitive to nesting-vector resolution and to the opening of new interlayer scattering channels.

    Authors: We agree that additional convergence tests are important to establish the robustness of the narrow incommensurate window. In our previous work on the pure 2D Hubbard model, we have verified convergence with respect to the number of patches (typically 24-48) and frequency discretization. For the extended model with finite t_perp, the discretization scheme is identical, with the interlayer hopping incorporated into the propagator. To address this concern, we will include in the revised manuscript a dedicated paragraph or appendix showing results for different patch numbers (e.g., 16, 24, 32) at representative intermediate t_perp values, confirming that the position of the incommensurate instability remains stable within the resolution. revision: yes

  2. Referee: [Method and truncation discussion] The truncation to the two-particle vertex (standard in fRG) is retained without additional justification or error estimates when the interlayer term couples momenta along the stacking direction. Because the claimed window is narrow, it is essential to demonstrate that the truncation does not artificially stabilize or shift the incommensurate instabilities relative to the t_perp = 0 limit.

    Authors: The truncation to the two-particle vertex is the standard approximation in fRG applications to the Hubbard model, justified by the hierarchy of scales in the flow and the fact that higher-order vertices are generated but often remain small in the relevant regimes. The interlayer hopping enters primarily through the single-particle dispersion and the loop integrals, without altering the vertex truncation itself. We note that for small t_perp the perturbation is weak, and the instabilities evolve continuously from the 2D case. However, we acknowledge the referee's point regarding the narrowness of the window and the lack of explicit error estimates. In the revision, we will add a discussion in the methods section elaborating on the truncation and its expected validity for the interlayer coupling, including a qualitative argument why higher vertices are unlikely to shift the narrow region significantly. Full error estimates would require a different methodological approach beyond the scope of this work. revision: partial

Circularity Check

0 steps flagged

No circularity: standard fRG flow applied to extended surface geometry

full rationale

The paper defines an extension of the two-dimensional functional renormalization group to semi-infinite layered systems by adding an alternating interlayer hopping term to the surface Hubbard model and integrates the standard flow equations numerically. The central claim of a narrow incommensurate SDW/spin-bond window splitting the d-wave superconducting regime is obtained as an output of that integration rather than being presupposed by any self-definition, fitted-parameter renaming, or load-bearing self-citation. No equation or result reduces to its own input by construction; the truncation and discretization choices are methodological assumptions whose accuracy is external to the derivation itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard fRG truncation scheme being sufficient for the surface-plus-interlayer model; no new entities are postulated and the interlayer hopping parameters are scanned rather than fitted.

axioms (1)
  • domain assumption The functional renormalization group flow equations with standard truncation approximate the many-body physics of the Hubbard model on the surface layer
    Invoked throughout the application section of the abstract

pith-pipeline@v0.9.0 · 5444 in / 1211 out tokens · 51965 ms · 2026-05-16T13:57:50.553169+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

70 extracted references · 70 canonical work pages · 3 internal anchors

  1. [1]

    Hubbard, Electron correlations in narrow energy bands, Proceedings of the Royal Society of London

    J. Hubbard, Electron correlations in narrow energy bands, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences276, 238 (1963)

  2. [2]

    P. W. Anderson, The Resonating Valence Bond State in La2CuO4 and Superconductivity, Science235, 1196 (1987)

  3. [3]

    V. J. Emery, Theory of high-t c superconductivity in ox- ides, Phys. Rev. Lett.58, 2794 (1987). 7

  4. [4]

    D. P. Arovas, E. Berg, S. A. Kivelson, and S. Raghu, The hubbard model, Annual Review of Condensed Mat- ter Physics13, 239 (2022)

  5. [5]

    Qin, C.-M

    M. Qin, C.-M. Chung, H. Shi, E. Vitali, C. Hubig, U. Schollw¨ ock, S. R. White, and S. Zhang (Simons Col- laboration on the Many-Electron Problem), Absence of superconductivity in the pure two-dimensional hubbard model, Phys. Rev. X10, 031016 (2020)

  6. [6]

    Sch¨ afer, N

    T. Sch¨ afer, N. Wentzell, F.ˇSimkovic, Y.-Y. He, C. Hille, M. Klett, C. J. Eckhardt, B. Arzhang, V. Harkov, F. m. c.-M. Le R´ egent, A. Kirsch, Y. Wang, A. J. Kim, E. Kozik, E. A. Stepanov, A. Kauch, S. Andergassen, P. Hansmann, D. Rohe, Y. M. Vilk, J. P. F. LeBlanc, S. Zhang, A.-M. S. Tremblay, M. Ferrero, O. Parcollet, and A. Georges, Tracking the foot...

  7. [7]

    M. Qin, T. Sch¨ afer, S. Andergassen, P. Corboz, and E. Gull, The hubbard model: A computational perspec- tive, Annual Review of Condensed Matter Physics13, 275 (2022)

  8. [8]

    Exact evolution equation for the effective potential

    C. Wetterich, Exact evolution equation for the effective potential, Physics Letters B301, 90 (1993), 1710.05815

  9. [9]

    T. R. Morris, The exact renormalization group and ap- proximate solutions, International Journal of Modern Physics A09, 2411 (1994)

  10. [10]

    Salmhofer and C

    M. Salmhofer and C. Honerkamp, Fermionic renormal- ization group flows: Technique and theory, Progress of Theoretical Physics105, 1 (2001)

  11. [11]

    Metzner, M

    W. Metzner, M. Salmhofer, C. Honerkamp, V. Meden, and K. Sch¨ onhammer, Functional renormalization group approach to correlated fermion systems, Reviews of Mod- ern Physics84, 299 (2012)

  12. [12]

    Wang, Y.-Y

    W.-S. Wang, Y.-Y. Xiang, Q.-H. Wang, F. Wang, F. Yang, and D.-H. Lee, Functional renormalization group and variational monte carlo studies of the elec- tronic instabilities in graphene near 1 4 doping, Phys. Rev. B85, 035414 (2012)

  13. [13]

    Platt, W

    C. Platt, W. Hanke, and R. Thomale, Functional renor- malization group for multi-orbital fermi surface instabil- ities, Advances in Physics62, 453 (2013)

  14. [14]

    D. M. Kennes, J. Lischner, and C. Karrasch, Strong cor- relations andd+idsuperconductivity in twisted bilayer graphene, Phys. Rev. B98, 241407 (2018)

  15. [15]

    Classen, C

    L. Classen, C. Honerkamp, and M. M. Scherer, Com- peting phases of interacting electrons on triangular lat- tices in moir´ e heterostructures, Phys. Rev. B99, 195120 (2019)

  16. [16]

    Dupuis, L

    N. Dupuis, L. Canet, A. Eichhorn, W. Metzner, J. M. Pawlowski, M. Tissier, and N. Wschebor, The nonper- turbative functional renormalization group and its appli- cations, Physics Reports910, 1 (2021)

  17. [17]

    M. L. Kiesel, C. Platt, and R. Thomale, Unconventional fermi surface instabilities in the kagome hubbard model, Phys. Rev. Lett.110, 126405 (2013)

  18. [18]

    S´ anchez de la Pe˜ na,Competing orders in honeycomb Hubbard models with nonlocal Coulomb interactions: A functional renormalization group approach, Ph.D

    D. S´ anchez de la Pe˜ na,Competing orders in honeycomb Hubbard models with nonlocal Coulomb interactions: A functional renormalization group approach, Ph.D. thesis, RWTH Aachen University, Aachen (2018)

  19. [19]

    Klett, T

    M. Klett, T. Schwemmer, S. Wolf, X. Wu, D. Riegler, A. Dittmaier, D. Di Sante, G. Li, W. Hanke, S. Rachel, and R. Thomale, From highT c to lowT c: Multiorbital effects in transition metal oxides, Phys. Rev. B104, L100502 (2021)

  20. [20]

    Klebl, Q

    L. Klebl, Q. Xu, A. Fischer, L. Xian, M. Claassen, A. Ru- bio, and D. M. Kennes, Moir´ e engineering of spin–orbit coupling in twisted platinum diselenide, Electronic Struc- ture4, 014004 (2022)

  21. [21]

    W. Qin, C. Huang, T. Wolf, N. Wei, I. Blinov, and A. H. MacDonald, Functional renormalization group study of superconductivity in rhombohedral trilayer graphene, Phys. Rev. Lett.130, 146001 (2023)

  22. [22]

    Klebl, A

    L. Klebl, A. Fischer, L. Classen, M. M. Scherer, and D. M. Kennes, Competition of density waves and super- conductivity in twisted tungsten diselenide, Phys. Rev. Res.5, L012034 (2023)

  23. [23]

    Beyer, J

    J. Beyer, J. B. Profe, L. Klebl, T. Schwemmer, D. M. Kennes, R. Thomale, C. Honerkamp, and S. Rachel, Rashba spin-orbit coupling in the square-lattice hub- bard model: A truncated-unity functional renormaliza- tion group study, Phys. Rev. B107, 125115 (2023)

  24. [24]

    Y. He, K. Yang, J. B. Profe, E. J. Bergholtz, and D. M. Kennes, Superconductivity of repulsive spinless fermions with sublattice potentials, Phys. Rev. Res.5, L012009 (2023)

  25. [25]

    J. B. Profe, L. Klebl, F. Grandi, H. Hohmann, M. D¨ urrnagel, T. Schwemmer, R. Thomale, and D. M. Kennes, Kagome hubbard model from a functional renor- malization group perspective, Phys. Rev. Res.6, 043078 (2024)

  26. [26]

    J. B. Profe, L. C. Rhodes, M. D¨ urrnagel, R. Bisset, C. A. Marques, S. Chi, T. Schwemmer, R. Thomale, D. M. Kennes, C. A. Hooley, and P. Wahl, Magic angle of sr 2ruo4: Optimizing correlation-driven superconduc- tivity, Phys. Rev. Res.6, 043057 (2024)

  27. [27]

    Fischer, L

    A. Fischer, L. Klebl, V. Cr´ epel, S. Ryee, A. Rubio, L. Xian, T. O. Wehling, A. Georges, D. M. Kennes, and A. J. Millis, Theory of intervalley-coherent afm order and topological superconductivity in twse2, Phys. Rev. X15, 041055 (2025)

  28. [28]

    J. Beck, J. Bodky, M. D¨ urrnagel, R. Thomale, J. Ingham, L. Klebl, and H. Hohmann, Kekul´ e Order from Diffuse Nesting near Higher-Order Van Hove Points, Physical Review Letters136, 106503 (2026)

  29. [29]

    D¨ urrnagel, H

    M. D¨ urrnagel, H. Hohmann, A. Maity, J. Seufert, M. Klett, L. Klebl, and R. Thomale, Altermagnetic phase transition in a lieb metal, Phys. Rev. Lett.135, 036502 (2025)

  30. [30]

    Y. Guo, J. Cenker, A. Fischer, D. Mu˜ noz-Segovia, J. Pack, L. Holtzman, L. Klebl, K. Watanabe, T. Taniguchi, K. Barmak, J. Hone, A. Rubio, D. M. Kennes, A. J. Millis, A. Pasupathy, and C. R. Dean, An- gle evolution of the superconducting phase diagram in twisted bilayer wse2, arXiv:2512.06265, accepted in Na- ture (2025)

  31. [31]

    Ehrlich and C

    J. Ehrlich and C. Honerkamp, Functional renormaliza- tion group for fermion lattice models in three dimensions: Application to the hubbard model on the cubic lattice, Phys. Rev. B102, 195108 (2020)

  32. [32]

    J. B. Profe, S. Beck, D. M. Kennes, A. Georges, and O. Gingras, Competition between d-wave superconduc- tivity and magnetism in uniaxially strained Sr2RuO4, npj Quantum Materials9, 53 (2024)

  33. [33]

    C. Bigi, M. D¨ urrnagel, L. Klebl, A. Consiglio, G. Pokharel, M. Zonno, F. Bertran, P. Le F` evre, T. Jaouen, H. C. Tchouekem, P. Turban, A. De Vita, J. A. Miwa, J. W. Wells, D. Oh, R. Comin, R. Thomale, 8 I. Zeljkovic, B. R. Ortiz, S. D. Wilson, G. Sangiovanni, F. Mazzola, and D. Di Sante, Pomeranchuk instability from electronic correlations in CsTi3Bi5 ...

  34. [34]

    P. A. Lee, From high temperature superconductivity to quantum spin liquid: Progress in strong correla- tion physics, Reports on Progress in Physics71, 012501 (2007)

  35. [35]

    Rachel, Interacting topological insulators: A review, Reports on Progress in Physics81, 116501 (2018)

    S. Rachel, Interacting topological insulators: A review, Reports on Progress in Physics81, 116501 (2018)

  36. [36]

    W. P. Su, J. R. Schrieffer, and A. J. Heeger, Solitons in polyacetylene, Phys. Rev. Lett.42, 1698 (1979)

  37. [37]

    Honerkamp and M

    C. Honerkamp and M. Salmhofer, Temperature-flow renormalization group and the competition between su- perconductivity and ferromagnetism, Phys. Rev. B64, 184516 (2001)

  38. [38]

    Kalkstein and P

    D. Kalkstein and P. Soven, A green’s function theory of surface states, Surface Science26, 85 (1971)

  39. [39]

    M. P. L. Sancho, J. M. L. Sancho, J. M. L. Sancho, and J. Rubio, Highly convergent schemes for the calculation of bulk and surface green functions, Journal of Physics F: Metal Physics15, 851 (1985)

  40. [40]

    Husemann and M

    C. Husemann and M. Salmhofer, Efficient parametriza- tion of the vertex function, Ω scheme, and thet, t ′ hub- bard model at van hove filling, Phys. Rev. B79, 195125 (2009)

  41. [41]

    Lichtenstein, D

    J. Lichtenstein, D. S´ anchez de la Pe˜ na, D. Rohe, E. Di Napoli, C. Honerkamp, and S. Maier, High-performance functional Renormalization Group calculations for in- teracting fermions, Computer Physics Communications 213, 100 (2017)

  42. [42]

    Honerkamp, Efficient vertex parametrization for the constrained functional renormalization group for effective low-energy interactions in multiband systems, Phys

    C. Honerkamp, Efficient vertex parametrization for the constrained functional renormalization group for effective low-energy interactions in multiband systems, Phys. Rev. B98, 155132 (2018)

  43. [43]

    Honerkamp, H

    C. Honerkamp, H. Shinaoka, F. F. Assaad, and P. Werner, Limitations of constrained random phase approximation downfolding, Phys. Rev. B98, 235151 (2018)

  44. [44]

    C. J. Eckhardt, C. Honerkamp, K. Held, and A. Kauch, Truncated unity parquet solver, Phys. Rev. B101, 155104 (2020)

  45. [45]

    Hille, F

    C. Hille, F. B. Kugler, C. J. Eckhardt, Y.-Y. He, A. Kauch, C. Honerkamp, A. Toschi, and S. Andergassen, Quantitative functional renormalization group descrip- tion of the two-dimensional hubbard model, Phys. Rev. Research2, 033372 (2020)

  46. [46]

    Klebl, D

    L. Klebl, D. M. Kennes, and C. Honerkamp, Functional renormalization group for a large moir´ e unit cell, Physical Review B102, 085109 (2020)

  47. [47]

    J. B. Profe, C. Honerkamp, S. Achilles, and D. M. Kennes, Electronic instabilities in penrose quasicrystals: Competition, coexistence, and collaboration of order, Phys. Rev. Res.3, 023180 (2021)

  48. [48]

    J. B. Profe and D. M. Kennes, TU 2FRG: A scalable approach for truncated unity functional renormalization group in generic fermionic models, The European Phys- ical Journal B95, 60 (2022)

  49. [49]

    Honerkamp, M

    C. Honerkamp, M. Salmhofer, N. Furukawa, and T. M. Rice, Breakdown of the landau-fermi liquid in two di- mensions due to umklapp scattering, Phys. Rev. B63, 035109 (2001)

  50. [50]

    Honerkamp, Electron-doping versus hole-doping in the 2D t-t’ Hubbard model, The European Physical Jour- nal B - Condensed Matter and Complex Systems21, 81 (2001)

    C. Honerkamp, Electron-doping versus hole-doping in the 2D t-t’ Hubbard model, The European Physical Jour- nal B - Condensed Matter and Complex Systems21, 81 (2001)

  51. [51]

    Beyer, J

    J. Beyer, J. B. Hauck, and L. Klebl, Reference results for the momentum space functional renormalization group, Eur. Phys. J. B95, 65 (2022)

  52. [52]

    Karrasch,The functional renormalization group for zero-dimensional quantum systems in and out of equilib- rium, Ph.D

    C. Karrasch,The functional renormalization group for zero-dimensional quantum systems in and out of equilib- rium, Ph.D. thesis, RWTH Aachen University, Aachen (2010)

  53. [53]

    Q. Wu, S. Zhang, H.-F. Song, M. Troyer, and A. A. Soluyanov, Wanniertools: An open-source software pack- age for novel topological materials, Computer Physics Communications224, 405 (2018)

  54. [54]

    [53] not only provides surface Green’s functions, but at the same time local bulk Green’s functions

    In fact, the implementation provided in Ref. [53] not only provides surface Green’s functions, but at the same time local bulk Green’s functions. The bulk Green’s function can be used in the same manner as the surface Green’s function without further modifications to the FRG equa- tions to study intralayer interactions

  55. [55]

    For details on the numerical methods see Ref. [51]. The open-source implementation of the Hubbard-SSH model FRG [68] makes use of the publicly available divERGe library [69, 70]

  56. [56]

    Laubach, C

    M. Laubach, C. Platt, R. Thomale, T. Neupert, and S. Rachel, Density wave instabilities and surface state evolution in interacting weyl semimetals, Phys. Rev. B 94, 241102 (2016)

  57. [57]

    N. P. Armitage, E. J. Mele, and A. Vishwanath, Weyl and dirac semimetals in three-dimensional solids, Rev. Mod. Phys.90, 015001 (2018)

  58. [58]

    Xiong, X

    F. Xiong, X. Han, and C. Honerkamp, Spin susceptibili- ties in magnetic type-i and type-ii weyl semimetals, Phys. Rev. B104, 115151 (2021)

  59. [59]

    Kundu, H

    R. Kundu, H. A. Fertig, and A. Kundu, Broken symmetry and competing orders in weyl semimetal interfaces, Phys. Rev. B107, L041402 (2023)

  60. [60]

    Liu, P.-J

    Q. Liu, P.-J. Guo, X.-Y. Yue, Z.-K. Yi, Q.-X. Dong, H. Liang, D.-D. Wu, Y. Sun, Q.-J. Li, W.-L. Zhu, T.- L. Xia, X.-F. Sun, and Y.-Y. Wang, Observation of Sur- face Superconductivity in a 3D Dirac Material, Advanced Functional Materials32, 2208616 (2022)

  61. [61]

    Neupert, S

    T. Neupert, S. Rachel, R. Thomale, and M. Greiter, In- teracting surface states of three-dimensional topological insulators, Phys. Rev. Lett.115, 017001 (2015)

  62. [62]

    Lundgren, H

    R. Lundgren, H. Yerzhakov, and J. Maciejko, Nematic order on the surface of a three-dimensional topological insulator, Phys. Rev. B96, 235140 (2017)

  63. [63]

    Mæland, M

    K. Mæland, M. Bahari, and B. Trauzettel, Phonon- mediated intrinsic topological superconductivity in fermi arcs, Phys. Rev. B112, 104507 (2025)

  64. [64]

    Mechanism for Nodal Topological Superconductivity on PtBi$_2$ Surface

    K. Mæland, G. Sangiovanni, and B. Trauzettel, Mecha- nism for Nodal Topological Superconductivity on PtBi 2 Surface, arXiv:2512.09994 (2025)

  65. [65]

    Gariglio, A

    S. Gariglio, A. Fˆ ete, and J.-M. Triscone, Electron con- finement at the LaAlO3/SrTiO3 interface, Journal of Physics: Condensed Matter27, 283201 (2015)

  66. [66]

    M. S. Scheurer and J. Schmalian, Topological supercon- ductivity and unconventional pairing in oxide interfaces, Nature Communications6, 6005 (2015)

  67. [67]

    P. Th¨ ornig, JURECA: Data Centric and Booster Mod- ules implementing the Modular Supercomputing Archi- tecture at J¨ ulich Supercomputing Centre, Journal of large-scale research facilities JLSRF7, A182 (2021)

  68. [68]

    Klebl, Surface FRG for a cubic lattice,https:// 9 github.com/alcubierre-drive/surface-frg(2025)

    L. Klebl, Surface FRG for a cubic lattice,https:// 9 github.com/alcubierre-drive/surface-frg(2025)

  69. [69]

    J. B. Profe, D. M. Kennes, and L. Klebl, divERGe imple- ments various Exact Renormalization Group examples, SciPost Phys. Codebases , 26 (2024)

  70. [70]

    J. B. Profe, D. M. Kennes, and L. Klebl, Codebase release 0.5 for divERGe, SciPost Phys. Codebases , 26 (2024). 10 Appendix A: Robustness of thet ′ =−0.25iSDW To clarify that the presence of the iSDW state att ′ =−0.25 is not a numerical artifact, we perform truncated unity FRG simulations (via the divERGe library [69, 70]) for various momentum resolution...