{"id":"6c7daf5f-4d3a-4a0f-98dd-44e56a2890ab","arxiv_id":"2508.21127","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a single-flavor flat band with inversion symmetry, a local attraction between opposite-parity orbitals yields exact superconducting ground states, including topological pairing.","lead":"A new class of exactly solvable models shows that a single flat band with inversion symmetry can host superconductivity when two opposite-parity orbitals attract each other. The models produce two competing, particle-hole conjugate superconducting states and tie the superfluid stiffness to the band's quantum geometry.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Singular pairing form factor makes the exact ground state ill-defined without a regularization that is promised but not provided","rationale":"The reader identified exactly this as the weakest assumption: the pairing form factor f(k)=u_A*/u_B(-k) is singular where u_B(-k)=0, and the text promises but does not supply a complete proof that this is harmless. This is load-bearing because the exact ground states are constructed from η†, and if η† is not a well-defined operator on the Fock space, the statement that (η†)^N|vac> is a zero-energy ground state loses its meaning. The paper does provide a Weierstrass compactification for the numerical examples, which smooths f, but the central claim is general and not explicitly limited to compactified models. The SM's algebraic commutator proof is elegant and likely correct distributionally, but it does not address normalizability of the many-body states or finiteness of derived quantities like the two-particle mass. The concern does not amount to a demonstrated inconsistency: the construction may be fully salvageable by consistently working on the compactified lattice and taking the thermodynamic limit there. Therefore I do not change the reader's CONDITIONAL verdict; the paper should be accepted conditional on providing the missing regularization proof or explicitly restricting the exactness claim to the lattice-regularized models.","tokens_in":40758,"tokens_out":10519,"duration_ms":111248,"concrete_test":"Using the compactified lattice model defined by Eq. 13 (or a finite momentum-grid approximation of the continuum model Eq. 12 with a small-k cutoff Λ), compute the two-particle norm N_2 = ∑_k |f(k)|^2 and the two-particle inverse mass from SM Eq. 115 as functions of system size L or cutoff Λ. If N_2 diverges logarithmically with L (or m_2p^{-1} diverges as Λ→0), the exact ground state is not normalizable in the thermodynamic limit and the 'no divergence' claim fails without compactification. Conversely, if the compactified results converge and the continuum limit is finite after removing the cutoff in the correct order, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central exact-solvability claim rests on the pairing kernel f(k)=u_A*(k)/u_B(-k) defining a legitimate many-body operator η†. In the continuum two-band model (Eq. 12 with M odd), f ~ k^{-M}, so the two-particle state η†|vac> has norm ∫|f|^2, which diverges in both IR (k→0) and UV (k→∞) for M≥1 in 2D. Thus (η†)^N|vac> is not a normalizable vector in the Fock space of the thermodynamic limit without a regularization. The main text explicitly acknowledges singularities and promises 'we will show that this does not lead to any divergence in any physical observables', but no complete proof appears in the main text or the SM excerpt. The SM commutator argument (Sec. II.B) formally uses F(k)=-F(-k) and holds as a distributional identity, but it does not establish that η† is a well-defined operator or that the zero-energy states are legitimate. The lattice compactification via Eq. 13 replaces 1/k^M by [k0 ζ(k)]^M, which is smooth and vanishes at the former pole, curing the divergence; however, the general claim 'for any pair of orbitals...' is not restricted to such compactified models. Moreover, SM Eq. 115 for the two-particle mass contains ∫ |u_A|^2/|u_B|^2, which diverges for the continuum wavefunction of Eq. 12, suggesting some observables do diverge unless a cutoff is imposed. This missing regularization proof is the load-bearing weak point.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a class of exactly solvable models for superconducting ground states in a single isolated flat band with inversion symmetry. For two orbitals A,B of opposite parity, a local attractive density-density interaction plus a carefully chosen single-particle counterterm projects onto a positive semi-definite quartic Hamiltonian in the flat-band subspace. The central claim is that the pairing operator η† = ∫(u_A*/u_B(-k)) γ†_k γ†_{-k} commutes with the projected inter-orbital hopping terms in the required sense, so that (η†)^N|vac> are zero-energy ground states. The authors construct a BCS coherent state from these number-sector ground states, analyze the Majorana Chern number, derive an analytical variational upper bound for the superfluid stiffness, and supplement it with a numerical many-body-bootstrap lower bound. They illustrate the construction with a two-band continuum model of Chern number -M, a lattice compactification of that model, and a multi-orbital effective model related to rhombohedral graphene.","tokens_in":41156,"tokens_out":9254,"duration_ms":110403,"significance":"If the central construction is fully valid, this is a significant contribution to strong-coupling flat-band superconductivity: it gives a rare family of exact superconducting ground states whose pairing structure is controlled by the flat-band quantum geometry, with topological order inherited from the parent band and with two particle-hole conjugate pairing modes that may have different Majorana Chern numbers. The paper contains several concrete strengths: the positive-semidefinite plus commutation proof is explicit; the two-particle sector is solved exactly; the variational stiffness bound is derived from a gauge-response calculation, and the bootstrap lower bound provides a genuine numerical check. The construction is model-building rather than circular: the solvability is engineered by design, and no parameters are fitted. However, the validity of the main theorem is currently conditioned on a regularization of the singular pairing form factor f(k)=u_A*/u_B(-k), which the manuscript promises but does not supply. That issue is load-bearing for the exact-ground-state claim in the continuum formulation and for the general statement 'for any pair of orbitals.'","major_comments":[{"comment":"The pairing kernel f(k)=u_A*(k)/u_B(-k) is singular when u_B(-k)=0. In the primary continuum example, Eq. (12) with M odd has f(k) ~ k^{-M} for A=1,B=2, so for M=1 the norm of η†|vac> is ∫ d^2k |f|^2, which diverges logarithmically in both the infrared and the ultraviolet. Thus (η†)^N|vac> is not a normalizable vector in the thermodynamic-limit Fock space without an additional prescription. The main text explicitly promises to show that this does not cause divergences in physical observables, but the provided main text and SM excerpt do not contain such a proof. The commutator argument in SM Sec. II.B holds as a distributional identity, but it does not establish that η† is a well-defined operator. SM Eq. (115) also contains the denominator ∫ |u_A|^2/|u_B|^2, which diverges for the wavefunction of Eq. (12), so the assertion that physical quantities are automatically finite cannot be taken","section":"Main text, Eq. (5) and Eq. (12); SM Sec. II.B and Eq. (115)"},{"comment":"The claim 'for any pair of orbitals A,B with opposite parities' is broader than what is established. The single-band construction requires F(k)=u_A*(k)/u_B(-k) to define a legitimate many-body operator with the antisymmetry F(k)=-F(-k). When u_B has zeros, the commutator identity [η†, P S_R P]=0 is only formal unless the singularities are regulated. In the two-band model Eq. (12), the two dual choices A=1,B=2 and A=2,B=1 give form factors k^{-M} and k^M respectively; the latter is regular while the former is not. The particle-hole dual pair are therefore not on equal footing in the continuum model. The authors should either extend the regularization proof to cover the singular dual or explicitly restrict the exact-solvability theorem to models where f(k) is a bounded, or at least locally integrable, function with controlled singularities, with all physical observables shown to be finite.","section":"Main text, Eq. (5) and Eq. (12); SM Sec. III.B"},{"comment":"The variational stiffness upper bound is derived using the parity structure of ∂ ln(u_A* u_B) and δ ln f, and the statement that δf=0 is optimal at A=0. This derivation assumes that all momentum-space integrals are absolutely convergent and that integration by parts has no boundary terms. When f has poles, the intermediate expressions (e.g., the two-particle mass in SM Eq. (115)) are not manifestly finite, so the bound D_s,var in the main text has the same regularization problem. The numerical plots in Fig. 2 use the compactified lattice model, but the analytical formula Eq. (11) is presented as a general result. The manuscript should state the domain of validity of Eq. (11) and supply the corresponding regularization for the continuum case.","section":"SM Sec. IV.E and main text Eq. (11)"}],"minor_comments":[{"comment":"The phrase 'flavorless orbitals' is unclear; presumably 'flavor-less' or 'same-flavor' is intended. Please clarify.","section":"Abstract"},{"comment":"The notation V/V in the counterterm expression is confusing. Use a distinct symbol for the number of unit cells or the BZ volume, e.g., N_c or Ω, to avoid division of identical letters.","section":"SM Eq. (68)"},{"comment":"The parameter z controls the filling through ν(Re z) = ⟨1⟩_z, but the domain of Re z (e.g., whether z→±∞ corresponds to the dilute limits) is stated only in words. It would help to define the limiting fillings explicitly.","section":"Main text Eq. (6) and Eq. (8)"},{"comment":"The statement that the maximum stiffness is achieved when the 'hot spots' of quantum geometry are 'marginally filled' is empirical. Please define 'marginal filling' precisely and state the numerical criterion used in the figure.","section":"Fig. 2 and Fig. 3"},{"comment":"The Schur-complement truncation to the A orbitals is acknowledged to be uncontrolled. The Chern-number predictions for rhombohedral graphene should be stated as properties of the effective model, not as quantitative predictions for the original multi-orbital Hamiltonian.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a strong model-building contribution with a clear construction, an exact two-particle solution, and a genuine numerical bootstrap lower bound. The central concern is the singular pairing form factor: the main text promises a proof that singularities do not affect physical observables, but that proof is absent from the provided main text and SM excerpt. This is a fixable but load-bearing issue: the authors should either provide a rigorous regularization and demonstrate finiteness of the relevant observables, or restrict the theorem to the compactified/lattice class where f is well defined. If that is done, I would be supportive of publication in this journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core construction is solid and worth taking seriously. The paper goes cleanly beyond the earlier QGN models by allowing non-uniform, momentum-dependent pairing in a single flavorless flat band. The PSD-plus-commutation logic is clear: the projected Hamiltonian is a positive sum of operators that move past the pairing operator and annihilate the vacuum, so zero-energy ground states in each number sector are exact. The particle-hole dual pairing modes, with the same angular momentum but possibly different Chern numbers, are a nice twist and are well motivated by the two-dome phenomenology of rhombohedral graphene. I also appreciate that the stiffness upper bound is derived analytically and the lower bound comes from a bootstrap calculation, including explicit finite-size scaling. That is real evidence, not hand-waving.\n\nNow the soft spots, in proportion. The biggest one is exactly the stress-test concern: the pairing form factor f(k) = u_A*(k)/u_B(-k) is singular where u_B(-k)=0, and the main text promises that this does not lead to divergences in any physical observable, but I could not find a complete proof in the main text or the SM. In the continuum model of Eq. (12), f ~ 1/k^M, so the two-particle state has divergent norm for M >= 1, and SM Eq. (115) contains an integral of |u_A|^2/|u_B|^2 that also diverges. That is not a minor detail: without a regularization, (eta^dagger)^N |vac> is not a normalizable state in the thermodynamic limit, so calling it an exact ground state is not rigorous. The lattice compactification via Eq. (13) fixes this by making the form factor smooth and vanishing at the former poles, and the numerical results in Figs. 2 and 3 are for that regularized model. So the concrete results for those models likely stand. But the general claim in the abstract and the beginning of the 'Exact solvability' section is over-broad as stated. The authors should either restrict the general statement to regularized models or supply a rigorous distributional treatment showing that all physical observables are finite.\n\nThe other soft spots are minor. The rhombohedral-graphene connection uses a Schur-complement approximation that is flagged as uncontrolled, and the bootstrap lower bound lacks details on uncertainties and reproducibility. Both are fixable in revision.\n\nThis is a serious paper deserving refereeing. I would accept it after the regularization issue is addressed, not before.","headline":"A genuinely new exact-solvability construction for flat-band superconductivity, with a real gap in the treatment of the singular pairing form factor that needs to be fixed before the general claim is accepted.","tokens_in":774,"tokens_out":812,"would_cite":true,"duration_ms":35135,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A local attraction between opposite-parity orbitals makes flat-band superconductivity exactly solvable.","keywords":["flat-band superconductivity","exactly solvable models","inversion symmetry","quantum geometry","topological superconductivity","particle-hole duality","superfluid stiffness","rhombohedral graphene"],"falsifier":"Compute the expectation value ⟨(η)^N H (η†)^N⟩ on a finite lattice using the Weierstrass compactified model at a filling where u_B(-k) has a zero; if the energy density does not vanish as system size grows, or if the BCS state norm diverges, the exact-solvability claim fails for that regularization.","tokens_in":40646,"feed_emoji":"⚛","tokens_out":6720,"duration_ms":64731,"temperature":0.7,"pith_summary":"This paper claims that a single flavor-polarized flat band with inversion symmetry can host superconductivity that is exactly solvable: if two local orbitals have opposite inversion parities, a local attractive interaction between them, after projection onto the flat band, has zero-energy ground states of the form (η†)^N|vac>, with η† built from the band's Bloch wavefunctions. This matters because superconductivity in rhombohedral graphene and twisted MoTe2 appears in flavor-polarized, nearly flat bands, where standard weak-coupling BCS reasoning does not apply. The construction yields exact statements about the superconducting state: its stiffness is proportional to the attraction strength, it inherits the parent band's topology, and each model has two nearly degenerate pairing modes related by particle-hole duality that may have different Chern numbers. These features reproduce, at a qualitative level, the two distinct superconducting domes observed in rhombohedral graphene.","feed_headline":"Local attraction makes flat-band superconductivity exactly solvable","feed_subtitle":"Exact zero-energy states appear in any inversion-symmetric flat band; stiffness scales with interaction strength.","key_machinery":"The load-bearing object is the ratio-symmetric pairing form factor f(k)=u_A*(k)/u_B(-k), which inversion symmetry makes antisymmetric, f(k)=-f(-k). Around this form factor the paper constructs the local two-orbital attraction and its counterterm so that the projected Hamiltonian is a sum of positive semi-definite operators S_R† S_R with S_R=c†_RA c_RB; the commutation relation [S_R, η†]=0 in the flat-band subspace is what makes (η†)^N|vac> exact zero-energy ground states. The second key object is the particle-hole/time-reversal duality (C and T) that maps the A↔B exchanged Hamiltonian to the original one, relating fillings ν↔1-ν and giving two nearly degenerate SC modes. For explicit models","core_discovery":"The central claim is that for any pair of orbitals A,B with opposite inversion parities (p_A p_B = -1) in an isolated flat band, the projected Hamiltonian H_AB = -V Σ_R n_RA n_RB plus single-particle counterterms is exactly solvable. It is positive semi-definite, and the pairing operator η† = ∫ dk [u_A*(k)/u_B(-k)] γ†_k γ†_-k commutes with each projected local term in the sense that lets the Hamiltonian act on (η†)^N|vac> and then annihilate the vacuum, so these states have zero energy and are ground states. The ground states form a BCS family |z> with arbitrary phase, breaking U(1); in two dimensions their Majorana Chern number is C_maj = C_BdG + C_FB + 2C_<, and for two-orbital models C_Bd","pith_inferences":["The construction should extend to systems with approximate or emergent inversion symmetry, such as the SO(2)-symmetric low-energy models of rhombohedral graphene or twisted MoTe2, since only the nesting relation p_α u_α(k)=u_α(-k)e^{iξ(k)} is needed; the authors hint at this but do not carry it out.","The predicted Chern numbers for L-layer rhombohedral graphene (C_maj=1 and C_maj=2L-3 for the two domes) could be tested by thermal Hall or edge-transport measurements if the two domes are indeed the particle-hole conjugate pair.","The exact states could serve as benchmarks for numerical methods on strongly coupled flat-band superconductors, for example by comparing the bootstrap lower bound with the variational upper bound at finite sizes.","A finite-size scaling study near the singular points of the form factor, where u_B(-k)=0, would sharpen the thermodynamic-limit argument; the paper's lattice compactification makes this a concrete numerical test."],"forward_implications":["Any inversion-symmetric isolated flat band with a local attraction between opposite-parity orbitals has exact zero-energy superconducting ground states at all even fillings, independent of band topology.","When the parent flat band is topological, the superconducting state is topological too (C_maj = C_FB in the two-orbital lower-band case), implying chiral Majorana edge modes.","Each model contains two nearly degenerate SC modes related by particle-hole duality, with D_s(ν)=D_s(1-ν), so the two modes appear as two SC domes tunable by charge-transfer gap and density.","The superfluid stiffness is bounded above and below by the interaction scale V times a quantum-geometric factor; in the dilute limit the upper bound coincides with the exact two-particle bound-state mass.","The exact BCS ground states have off-diagonal long-range order and arbitrary phase, so the model exhibits genuine U(1) symmetry breaking despite being strongly coupled."],"supporting_citations":[{"why":"Experimental report of chiral superconductivity in rhombohedral graphene that motivates the single-flavor flat-band setup.","marker":"[1]"},{"why":"Second experimental report of superconductivity in rhombohedral graphene, used together with [1] to motivate the two-dome phenomenology.","marker":"[2]"},{"why":"Classic pseudopotential models whose exact-solvability strategy this work extends to superconducting flat bands.","marker":"[27, 28]"},{"why":"Result that flat-band superfluid stiffness is proportional to interaction strength, which the present bounds generalize.","marker":"[51]"},{"why":"Prior quantum geometric nesting construction that this paper generalizes by allowing momentum-dependent order parameters.","marker":"[67]"},{"why":"Quantum many-body bootstrap method used to obtain a rigorous lower bound on the superfluid stiffness.","marker":"[74]"},{"why":"Weierstrass-function compactification used to regularize the continuum flat-band model on a hexagonal lattice.","marker":"[92]"},{"why":"BHZ model whose flattened version is the paper's two-band example with Chern number -M.","marker":"[95]"},{"why":"Effective two-band model of rhombohedral graphene whose quantum geometry the M=5 example mimics.","marker":"[100]"}],"fun_headline_variants":["Exact flat-band superconductors from local attraction","Inversion symmetry enables exact flat-band pairing","Superfluid stiffness tied to interaction in flat bands","Chiral flat-band superconductors: exact ground states"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The pairing form factor f(k)=u_A*(k)/u_B(-k) can be singular where u_B(-k)=0; the argument that these singularities cause no divergence in physical observables must hold, for example through the lattice compactification, for the exact ground states to be well defined in the thermodynamic limit.","fun_headline_variants_meta":{"raw":{"variants":["Exact flat-band superconductors from local attraction","Inversion symmetry enables exact flat-band pairing","Superfluid stiffness tied to interaction in flat bands","Chiral flat-band superconductors: exact ground states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00019,"raw_usage":{"total_tokens":1202,"prompt_tokens":796,"completion_tokens":406,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":347}},"tokens_in":540,"tokens_out":406,"duration_ms":4438,"temperature":1.0,"reasoning_tokens":347,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T14:34:30.867891+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the expectation value ⟨(η)^N H (η†)^N⟩ on a finite lattice using the Weierstrass compactified model at a filling where u_B(-k) has a zero; if the energy density does not vanish as system size grows, or if the BCS state norm diverges, the exact-solvability claim fails for that regularization.","supporting_citations":[{"cited_title":"Bootstrapping Flat-band Superconductors: Rigorous Lower Bounds on Superfluid Stiffness","cited_arxiv_id":"2506.18969","evidence_quote":"Quantum many-body bootstrap method used to obtain a rigorous lower bound on the superfluid stiffness."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Weierstrass-function compactification used to regularize the continuum flat-band model on a hexagonal lattice."},{"cited_title":"Herzog-Arbeitman, Y","cited_arxiv_id":null,"evidence_quote":"Effective two-band model of rhombohedral graphene whose quantum geometry the M=5 example mimics."}],"review_version":1}