arxiv: 2605.10387 · v2 · submitted 2026-05-11 · ❄️ cond-mat.supr-con · cond-mat.str-el

Recognition: 2 theorem links

· Lean Theorem

Layer-antisymmetric pair-phase resonance at the bonding-antibonding splitting in the AA-stacked bilayer attractive Hubbard model

Yogeshwar Prasad

Authors on Pith no claims yet

Pith reviewed 2026-05-13 00:44 UTC · model grok-4.3

classification ❄️ cond-mat.supr-con cond-mat.str-el

keywords bilayer superconductorpair phase modeattractive Hubbard modelcollective excitationsGaussian fluctuationsbonding-antibonding splittingAA-stacked bilayerLeggett mode

0 comments

The pith

In an AA-stacked bilayer superconductor the layer-antisymmetric pair phase supports an in-gap collective resonance at twice the interlayer hopping.

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

This paper establishes that the layer-antisymmetric pair-phase fluctuations in the AA-stacked attractive Hubbard honeycomb bilayer produce a collective pole inside the superconducting gap at frequency 2t_h. This frequency equals the bonding-antibonding splitting of the single-particle bands, and the coincidence arises algebraically because the antisymmetric phase bubble equals the static symmetric bubble at that specific frequency. A reader would care because this sets the resonance energy by the bare hybridization parameter instead of the interaction strength that usually governs Josephson modes. The finding is exact within Gaussian theory at any doping and is backed by numerical Bogoliubov-de Gennes checks.

Core claim

Working at the Gaussian fluctuation level, the antisymmetric pair-phase channel hosts an in-gap collective pole at 2t_h, the bonding-antibonding band splitting. At this frequency the antisymmetric phase bubble reduces pointwise in momentum space to the static symmetric phase bubble that enforces the in-phase Goldstone pole. The resonance scale is therefore fixed by single-particle hybridization rather than by the interaction-driven Josephson coupling. The diagonal kernel zero is exact at any chemical potential, while the full amplitude-phase pole coincides at half filling.

What carries the argument

The pointwise momentum-space reduction of the antisymmetric phase bubble to the static symmetric phase bubble at frequency 2t_h.

If this is right

The resonance energy scale is determined by the single-particle interlayer hopping rather than the Josephson coupling.
The antisymmetric phase-channel kernel zero is exact within Gaussian theory at any chemical potential.
The full coupled amplitude-phase pole coincides with the kernel zero at half filling and tracks it closely away from half filling.
A layer-imbalance drive overlaps with the pair-phase sector at the Gaussian level.
The excitation is Raman-forbidden by inversion symmetry, suggesting layer-odd detection methods.

Where Pith is reading between the lines

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

If the algebraic reduction survives beyond Gaussian fluctuations, the resonance would remain pinned at 2t_h in more complete treatments.
Layer-bias drives in cold-atom realizations could reveal a response feature near the sub-kilohertz scale for typical optical lattice parameters.
Similar resonances might appear in other bilayer systems where single-particle hybridization dominates over interaction effects.

Load-bearing premise

The exact algebraic reduction of the antisymmetric bubble holds only at the Gaussian fluctuation level, and higher-order corrections could shift the resonance away from 2t_h.

What would settle it

A direct evaluation of the antisymmetric phase bubble at frequency equal to 2t_h that fails to match the static symmetric phase bubble pointwise in momentum space would disprove the algebraic mechanism.

Figures

Figures reproduced from arXiv: 2605.10387 by Yogeshwar Prasad.

**Figure 2.** Figure 2: FIG. 2. Noninteracting band structure and density of states [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Layer-antisymmetric pair-phase response in the [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Mean-field pairing scale and half-filled pseu [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Layer-imbalance–phase cross-susceptibility [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Layer-antisymmetric pair-phase response in the [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗

read the original abstract

The relative phase between the two pair condensates of a bilayer s-wave superconductor is a collective degree of freedom distinct from the usual in-phase Anderson-Bogoliubov mode. Working at the Gaussian fluctuation level for the AA-stacked attractive-Hubbard honeycomb bilayer, we show analytically that the layer-antisymmetric pair-phase channel hosts an in-gap collective pole at twice the single-particle interlayer hopping, $2t_h$, precisely the bonding-antibonding band splitting. The mechanism is algebraic: at this frequency, the antisymmetric phase bubble reduces pointwise in momentum space to the static symmetric phase bubble that enforces the in-phase Goldstone pole. The resulting resonance scale is therefore fixed by the single-particle hybridization, rather than by the interaction-driven Josephson coupling that controls the canonical Leggett mode. The identity is verified numerically by direct Bogoliubov-de Gennes calculations. The diagonal antisymmetric phase-channel kernel zero is exact within Gaussian theory at any chemical potential; the full coupled amplitude-phase pole coincides with it at half filling and tracks it closely away from half filling. The excitation is Raman-forbidden by inversion, which motivates layer-odd probes. We find that a layer-imbalance drive has finite Gaussian-level overlap with the pair-phase sector, suggesting a possible cold-atom layer-bias response feature near the sub-kilohertz scale for typical optical-lattice parameters.

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.

Desk Editor's Note private letter to a colleague

The paper pins an antisymmetric pair-phase resonance exactly at the single-particle bonding-antibonding splitting 2t_h via an algebraic identity in the Gaussian bubble for this bilayer Hubbard model.

read the letter

The central observation is that the layer-antisymmetric phase channel develops a collective pole right at twice the interlayer hopping because the antisymmetric bubble at that frequency reduces pointwise to the static symmetric bubble that already enforces the in-phase Goldstone mode. This fixes the scale by the single-particle hybridization rather than by any interaction strength or Josephson coupling, which distinguishes it from the usual Leggett mode discussion. The identity is shown analytically for the AA-stacked attractive Hubbard honeycomb bilayer and checked numerically with BdG calculations; the diagonal kernel zero holds exactly at any chemical potential, while the full coupled amplitude-phase pole tracks it closely and coincides at half filling. That algebraic step is clean and does not rely on fitted parameters or external data. The work is useful for anyone tracking collective modes in layered superconductors or looking for cold-atom signatures, since the resonance sits in the sub-kilohertz range for typical lattice parameters and has finite overlap with a layer-imbalance drive. The main limitation is that the entire analysis stays inside Gaussian fluctuations. Nothing in the manuscript addresses whether higher-order corrections or finite-size effects would shift or damp the pole, and the numerical verification is stated without reported system sizes, convergence checks, or scans away from the half-filled point. The Raman-forbidden character and the suggestion of layer-odd probes follow directly from the symmetry but inherit the same approximation. Readers working on multilayer superconductivity or optical-lattice realizations of Josephson physics will get the most from it. The claim is internally consistent and the supporting steps are reproducible within the stated framework, so the paper merits a serious referee rather than a desk rejection.

Referee Report

0 major / 3 minor

Summary. The paper analytically demonstrates within Gaussian fluctuation theory for the AA-stacked bilayer attractive Hubbard model on the honeycomb lattice that the layer-antisymmetric pair-phase bubble reduces pointwise in momentum space to the static symmetric phase bubble at frequency ω=2t_h (the bonding-antibonding splitting). This places an in-gap collective pole in the antisymmetric channel. The reduction is exact for the diagonal kernel at arbitrary chemical potential; the full coupled amplitude-phase pole coincides exactly at half filling and tracks the resonance closely away from it. Numerical verification is provided via Bogoliubov-de Gennes calculations, and the mode is Raman-forbidden but potentially accessible via layer-imbalance drives.

Significance. If the Gaussian-level identity holds, the work identifies a collective mode in bilayer superconductors whose scale is fixed by single-particle hybridization t_h rather than Josephson coupling, distinguishing it from the Leggett mode. The algebraic, parameter-free character of the bubble reduction at any μ is a notable strength, as is the direct BdG verification of the full pole. This provides a falsifiable prediction and motivates layer-odd probes in cold-atom systems. The result is internally consistent within its stated Gaussian scope; higher-order fluctuation effects lie outside the manuscript's claims.

minor comments (3)

[Abstract] The explicit algebraic steps demonstrating the pointwise bubble reduction are not shown in the abstract or summary text, which would allow independent verification of the identity without re-deriving the full Gaussian kernel.
The BdG numerical verification is stated but provides no details on lattice size, momentum discretization, extraction of the pole position, or quantitative deviation from 2t_h away from half filling.
No parameter scans or error estimates are included to illustrate robustness of the tracking behavior or sensitivity to filling.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive and accurate summary of the manuscript, as well as for the recommendation of minor revision. The assessment correctly identifies the algebraic pointwise reduction of the antisymmetric bubble and its implications for the resonance scale being set by single-particle hybridization rather than Josephson coupling.

Circularity Check

0 steps flagged

No significant circularity: algebraic identity is self-contained

full rationale

The paper derives an exact pointwise reduction of the antisymmetric phase bubble to the static symmetric phase bubble at ω=2t_h directly from the Gaussian fluctuation equations, with the resonance scale fixed by the input single-particle hopping t_h. This holds analytically at any μ for the diagonal kernel and is verified numerically via BdG without fitted parameters or load-bearing self-citations. The derivation chain is internal to the bubble formalism and does not reduce any prediction to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claim rests on the Gaussian fluctuation approximation applied to the attractive Hubbard model; no new entities are introduced and the only input scale is the single-particle hopping.

axioms (2)

domain assumption Gaussian fluctuation level suffices for the pair-phase channel
All analytic results are derived under the quadratic fluctuation approximation around the mean-field solution.
domain assumption AA-stacked attractive Hubbard model on honeycomb bilayer accurately captures the physics
The model choice and stacking geometry are taken as given for the derivation.

pith-pipeline@v0.9.0 · 5552 in / 1481 out tokens · 67547 ms · 2026-05-13T00:44:41.431117+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

at this frequency, the antisymmetric phase bubble reduces pointwise in momentum space to the static symmetric phase bubble... (Ea,s + Eb,s)² − (2th)² = 2 N_ph_ab,s
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The resulting resonance scale is therefore fixed by the single-particle hybridization, rather than by the interaction-driven Josephson coupling

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

47 extracted references · 47 canonical work pages

[1]

Each complex fluctuation field may in turn be de- composed into amplitude and phase components

This basis separates the layer-symmetric saddle from the layer-odd collective channel of interest. Each complex fluctuation field may in turn be de- composed into amplitude and phase components. Ac- cordingly, the Gaussian sector is organized into four channels: symmetric phase, antisymmetric phase, sym- metric amplitude, and antisymmetric amplitude. At h...

work page
[2]

The kinematic feature that drives the central identity is the rigid parity splitting ξa,s(k)−ξ b,s(k) = 2th,(6) independent of momentum, sublattice branch, and chem- ical potential. The block decomposition follows from AA layer parity, the layer-symmetric uniform saddle (∆a = ∆ b), and singlet pairing on a real saddle, and is therefore valid at any chemic...

work page
[3]

Micnas, J

R. Micnas, J. Ranninger, and S. Robaszkiewicz, Rev. Mod. Phys.62, 113 (1990)

work page 1990
[4]

Toschi, M

A. Toschi, M. Capone, and C. Castellani, Phys. Rev. B 72, 235118 (2005)

work page 2005
[5]

R. T. Scalettar, E. Y. Loh, J. E. Gubernatis, A. Moreo, S. R. White, D. J. Scalapino, R. L. Sugar, and E. Dagotto, Phys. Rev. Lett.62, 1407 (1989)

work page 1989
[6]

Paiva, Y

T. Paiva, Y. L. Loh, M. Randeria, R. T. Scalettar, and N. Trivedi, Phys. Rev. Lett.107, 086401 (2011)

work page 2011
[7]

R. A. Fontenele, N. C. Costa, R. R. dos Santos, and T. Paiva, Phys. Rev. B105, 184502 (2022)

work page 2022
[8]

Zhao and A

E. Zhao and A. Paramekanti, Phys. Rev. Lett.97, 230404 (2006)

work page 2006
[9]

N. C. Costa, M. V. Ara´ ujo, J. P. Lima, T. Paiva, R. R. dos Santos, and R. T. Scalettar, Phys. Rev. B97, 085123 (2018)

work page 2018
[10]

C. N. Yang, Phys. Rev. Lett.63, 2144 (1989)

work page 1989
[11]

Zhang, Phys

S.-C. Zhang, Phys. Rev. Lett.65, 120 (1990)

work page 1990
[12]

Bloch, J

I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008)

work page 2008
[13]

Prasad, A

Y. Prasad, A. Medhi, and V. B. Shenoy, Phys. Rev. A 89, 043605 (2014)

work page 2014
[14]

Prasad, Phys

Y. Prasad, Phys. Rev. B106, 184506 (2022)

work page 2022
[15]

M. Gall, C. F. Chan, N. Wurz, and M. K¨ ohl, Phys. Rev. Lett.124, 010403 (2020)

work page 2020
[16]

R. A. Fontenele, N. C. Costa, T. Paiva, and R. R. dos Santos, Phys. Rev. A110, 053315 (2024), arXiv:2408.17405 [cond-mat.supr-con]

work page arXiv 2024
[17]

Grubiˇ si´ c-Cabo, J

A. Grubiˇ si´ c-Cabo, J. C. Kotsakidis, Y. Yin, A. Tadich, M. Haldon, S. Solari, J. Riley, E. Huwald, K. M. Daniels, R. L. Myers-Ward, M. T. Edmonds, N. Med- hekar, D. K. Gaskill, and M. S. Fuhrer, Front. Nanotech- nol.5, 1333127 (2023)

work page 2023
[18]

A. L. Rakhmanov, A. V. Rozhkov, A. O. Sboychakov, and F. Nori, Phys. Rev. Lett.109, 206801 (2012)

work page 2012
[19]

Ghadimi, C

R. Ghadimi, C. Mondal, S. Kim, and B.-J. Yang, Phys. Rev. Lett.133, 196603 (2024)

work page 2024
[20]

Mondal, R

C. Mondal, R. Ghadimi, and B.-J. Yang, Phys. Rev. B 113, L081101 (2026)

work page 2026
[21]

A. J. Leggett, Prog. Theor. Phys.36, 901 (1966)

work page 1966
[22]

S. G. Sharapov, V. P. Gusynin, and H. Beck, Eur. Phys. J. B30, 45 (2002)

work page 2002
[23]

Blumberg, A

G. Blumberg, A. Mialitsin, B. S. Dennis, M. V. Klein, N. D. Zhigadlo, and J. Karpinski, Phys. Rev. Lett.99, 227002 (2007)

work page 2007
[24]

Cea and L

T. Cea and L. Benfatto, Phys. Rev. B94, 064512 (2016)

work page 2016
[25]

J. J. Cuozzo, W. Yu, P. Davids,et al., Nat. Phys.20, 1118 (2024)

work page 2024
[26]

N. A. Hackner and P. M. R. Brydon, Physical Review B 108, L220505 (2023)

work page 2023
[27]

B. A. Levitan, Y. Oreg, E. Berg, M. S. Rudner, and I. Iorsh, Physical Review Research6, 043170 (2024)

work page 2024
[28]

P. W. Anderson, Phys. Rev.112, 1900 (1958)

work page 1900
[29]

N. N. Bogoliubov, Sov. Phys. JETP7, 41 (1958)

work page 1958
[30]

See Supplemental Material at [URL inserted by pub- lisher] for detailed lattice and mean-field conventions, the parity-resolved derivation of the pointwise susceptibility identity, the full-basis 8×8 Bogoliubov–de Gennes ver- ification, finite-doping and full coupled amplitude–phase pole analysis, finite-temperature reference scales, finite- momentum disp...

work page
[31]

Shiba, Prog

H. Shiba, Prog. Theor. Phys.48, 2171 (1972)

work page 1972
[32]

Tinkham,Introduction to Superconductivity, 2nd ed

M. Tinkham,Introduction to Superconductivity, 2nd ed. (Dover, 2004) coherence factors discussed in§3.9

work page 2004
[33]

Matsunaga, N

R. Matsunaga, N. Tsuji, H. Fujita, A. Sugioka, K. Makise, Y. Uzawa, H. Terai, Z. Wang, H. Aoki, and R. Shimano, Science345, 1145 (2014)

work page 2014
[34]

Giorgianni, T

F. Giorgianni, T. Cea, C. Vicario, C. P. Hauri, W. K. Withanage, X. Xi, and L. Benfatto, Nat. Commun.10, 11 2235 (2019)

work page 2019
[35]

M. Gall, N. Wurz, J. Samland, C. F. Chan, and M. K¨ ohl, Nature589, 397 (2021)

work page 2021
[36]

Pekker and C

D. Pekker and C. M. Varma, Annu. Rev. Condens. Mat- ter Phys.6, 269 (2015)

work page 2015
[37]

Shimano and N

R. Shimano and N. Tsuji, Annu. Rev. Condens. Matter Phys.11, 103 (2020)

work page 2020
[38]

Tsuji and H

N. Tsuji and H. Aoki, Phys. Rev. B92, 064508 (2015)

work page 2015
[39]

V. L. Berezinskii, Sov. Phys. JETP34, 610 (1972)

work page 1972
[40]

J. M. Kosterlitz and D. J. Thouless, J. Phys. C: Solid State Phys.6, 1181 (1973)

work page 1973
[41]

N. D. Mermin and H. Wagner, Phys. Rev. Lett.17, 1133 (1966)

work page 1966
[42]

D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett.39, 1201 (1977)

work page 1977
[43]

Ambegaokar, B

V. Ambegaokar, B. I. Halperin, D. R. Nelson, and E. D. Siggia, Phys. Rev. B21, 1806 (1980). 12 Supplemental Material: Layer-antisymmetric pair-phase resonance at the bonding–antibonding splitting in the AA-stacked bilayer attractive Hubbard model Yogeshwar Prasad This Supplemental Material provides technical details supporting the main text. It summarizes...

work page 1980
[44]

Each term inP †,un − carries one bonding and one antibonding creation operator, so the contributions ofH th cancel pairwise and one finds [H th , P †,un − ] = 0

The interlayer hopping in the same basis reads Hth =t h X i,σ c† iaσciaσ −c † ibσcibσ ,(A17) where bonding (b) sits at energy−t h and antibonding (a) at +t h. Each term inP †,un − carries one bonding and one antibonding creation operator, so the contributions ofH th cancel pairwise and one finds [H th , P †,un − ] = 0. Thus the interlayer hopping does not...

work page
[45]

BdG blocks and Bogoliubov coefficients At half filling, the BdG Hamiltonian is block-diagonal in the layer-parity labelκ∈ {b, a} ≡ {−,+}and in the sublattice-band indexs∈ {+,−}. The block dispersions are ξκ,s(k) =κ t h +s|f(k)|,(A22) and the quasiparticle energies are Eκ,s(k) = q ξ2κ,s(k) + ∆2 0.(A23) The corresponding Bogoliubov coefficients satisfy u2 κ...

work page
[46]

Antisymmetric phase bubble in the parity-resolved basis In thephysicallayer⊗Nambu space, the symmetric and antisymmetric phase vertices are Γ+ϕ =1 layer ⊗τ Nambu y ,Γ −ϕ =σ layer z ⊗τ Nambu y .(A26) Equation (A26) is written in the physical-layer representation; after rotating to the bonding/antibonding parity basis used in the block calculation,σ layer z...

work page
[47]

Bare-band identities atω= 2t h At half filling, the parity-resolved algebra simplifies becauseξ a,s −ξ b,s = 2th is independent of bothkands. Using (Ea −E b)2 =E 2 a +E 2 b −2E aEb, E 2 =ξ 2 + ∆2 0,(A31) one finds (Ea,s −E b,s)2 = 4t2 h + 2ξa,sξb,s + 2∆2 0 −2E a,sEb,s.(A32) It follows that (Ea,s +E b,s)2 −(2t h)2 = 2 Ea,sEb,s +ξ a,sξb,s + ∆2 0 = 2N ph ab,...

work page