Pith. sign in

REVIEW 3 major objections 7 minor 13 references

Acoustic Plasmon Resonance: Breaking the Anderson Stiffness Paradigm in Quasi-Two-Dimensional Superconducting Films

T0 review · 3 major / 7 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A superconducting film of finite thickness supports an acoustic plasmon whose velocity is set by thickness and by the superconducting density, breaking the usual stiffness of the plasma spectrum and coupling the mode to light.

desk verdict Interesting mechanism for acoustic plasmons in superconducting films, but the fitted D0≈70 violates the model's own geometric bounds, so the quantitative agreement with experiment is not supported. read the letter →

arxiv 2608.11610 v1 pith:SJVIBFMT submitted 2026-08-12 cond-mat.supr-con

V. M. Kovalev , A. V. Chaplik This is my paper · ORCID
classification cond-mat.supr-con
keywords acousticplasmonsuperconductingfilmplasmastiffnessCarlson-Goldmanmodetransversechargedistributiondirty-limitconductivitysubgapcollectiveelectromagneticcoupling
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to establish that the plasma spectrum of a superconductor is not necessarily stiff: in a film of finite thickness, the superfluid and normal charge densities occupy different transverse profiles, and their spatial mismatch breaks local neutrality. Using the Poisson equation for these profiles together with kinetic-equation (clean) and dirty-limit conductivities, the authors derive an acoustic branch $\omega_k = s k$ whose velocity is $s = \omega_0 \sqrt{(N_s/N)(M_{so}+k a J_{ss})/(J_{oo}+k a)}$, which for $a \ll d$ reduces to $s = \omega_0 \sqrt{D_0 N_s d/N}$. This mode lives below the superconducting gap, carries a dipole moment, couples to external radiation, and its velocity scales as $\sqrt{d}$ and as $\sqrt{T_c - T}$ near $T_c$, matching recent optical and microwave observations on thin films. If correct, the result reconciles those observations with the basic stiffness constraints of superconductivity and separates the new mode cleanly from the quasi-neutral Carlson-Goldman mode.

What carries the argument

The central object is the set of transverse charge distribution functions $A_\alpha(z)$, normalized so $\int A_\alpha(z)\,dz = 1$, and their geometric form factors $J_{\alpha\beta}(k) = \int\!\int dz\,dz'\, A_\alpha(z)\,e^{-k|z-z'|}\,A_\beta(z')$. These enter a dispersion relation built from the 2D Coulomb potential and polarization operators $\Pi_s$, $\Pi_n$, and $\Pi_0$; in the clean limit the response comes from the kinetic equation, and in the dirty low-frequency limit from the standard dirty-limit conductivities. The crucial simplification is that for $a \ll d$ the combination $(M_{so} + k a J_{ss})/(J_{oo} + k a)$ becomes $D_0 k d$ with $D_0$ a fitted dimensionless constant, which turns what would be a gapped plasma oscillation into a linearly dispersing acoustic mode with velocity $s = \omega_0 \sqrt{D_0 N_s d/N}$. Damping is governed by the analogous structure factor $S(k) \approx S_0$.

What would settle it

Measure the acoustic mode in films of fixed areal density but different thickness: the theory requires $s \propto \sqrt{d}$, while a Carlson-Goldman-like mode is independent of $d$, and the mode must vanish in the strictly two-dimensional limit $d \to 0$. A second check is to resolve the dispersion $\omega_k = s k$ and compare the temperature dependence of the slope with the predicted $N_s(T)/N$; if the measured $s$ does not vanish as $\sqrt{T_c - T}$ near $T_c$, the mechanism is not the one described.

Watch

Extended reading notes

Core claim

The central claim is that a finite-thickness superconducting film supports a subgap collective mode with linear dispersion $\omega_k = s k$, where the velocity $s$ is controlled by the superconducting density, the film thickness, and the temperature. The mechanism is geometric: the normalized transverse charge distributions $A_s(z)$, $A_n(z)$, and $A_o(z)$ for superfluid, normal, and deep-lying electrons are not identical, so the long-range Coulomb interaction couples to a net local charge in an oscillation that would be neutral in the Carlson-Goldman picture. Substituting the low-frequency polarization operators into the full dispersion relation and taking $a \ll d$ gives $\omega_k = \omega_0 \sqrt{D_0 N_s k d/N}$, i.e., an acoustic branch whose velocity scales as $\sqrt{d}$ and, near $T_c$, as $\sqrt{T_c - T}$. The normal component is overdamped and supplies a damping rate $\Gamma = (\omega_0^2 \tau/2)(N_n/N) S_0$, which grows with temperature and shows a Hebel-Slichter-like peak. The paper states that the transverse charge mismatch is what gives the mode its dipole moment and strong coupling to external electromagnetic radiation.

Load-bearing premise

The load-bearing premise is that the superfluid and normal charge distributions really are arranged differently across the film thickness; the paper assumes this mismatch rather than deriving it from a model, and the constant $D_0$ that fixes the velocity is fitted to the measured value.

Editorial extensions

If this is right

  • Near $T_c$, the acoustic velocity obeys $s \propto \sqrt{T_c - T}$ in both clean and dirty limits, with an extra suppression by $\sqrt{\tau T_c}$ in a disordered film, giving an optical handle on the condensate fraction.
  • Because the velocity scales as $\sqrt{d}$ at fixed areal density, film thickness is a control parameter for the subgap resonance frequency, and strictly two-dimensional films should show no such mode.
  • If the bulk three-dimensional electron density is fixed instead, the frequency scales as $\omega_k \propto d$, a prediction the paper says is already checked for 5, 8, and 12 nm films.
  • The mode's damping contains a Hebel-Slichter-like coherence peak near $T_c$, so linewidth measurements can probe the same effect seen in nuclear magnetic resonance.
  • Unlike Carlson-Goldman modes, this acoustic plasmon has a dipole moment and couples to external radiation, so far-field optical and microwave spectroscopy can observe it directly.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • An extension the authors leave implicit is that $D_0$ encodes the unmodeled transverse profiles; deliberately engineering the film interfaces or disorder (capping layers, asymmetric substrates, irradiation) should change the acoustic velocity, and sample-to-sample variation of $D_0$ would expose what the fitted constant is hiding.
  • By analogy, any two-component system whose charge carriers have different transverse confining wavefunctions could show the same mechanism—for example, coupled electron-hole bilayers or superconductor/normal-metal hybrids—so the paper's argument may generalize beyond superconducting films.
  • The polarization dependence of the optical coupling should carry direct information about $A_s(z)$ and $A_n(z)$; measuring the oscillator strength versus incidence angle could reconstruct the mismatch and test the assumed geometry without a specific microscopic model.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 7 minor

Summary. The manuscript proposes a microscopic theory for an acoustic plasmon mode in quasi-two-dimensional superconducting films. The authors model a superconducting film of finite thickness with transverse charge distribution functions for superconducting, normal, and deep-lying electrons, and derive a general dispersion relation from linearized kinetic theory or Mattis-Bardeen conductivities. In the long-wavelength limit they obtain an acoustic branch with linear dispersion, velocity scaling as the square root of the film thickness and of (Tc - T) near the critical temperature, and claim quantitative agreement with the experiments of Andreeva et al. [7]. The paper explicitly states that the microscopic origin of the transverse profiles is an open question, and the structure-factor constants D0 and S0 are fitted to the experimental mode velocity and damping.

Significance. If the proposed mechanism were sound, it would challenge the Anderson-Higgs stiffness picture and explain the recent observation of an EM-active subgap acoustic mode in superconducting films. The paper offers a concrete analytical framework and correctly identifies a physically interesting geometric mechanism (transverse charge separation). However, the central quantitative prediction rests on a fitted constant D0 that violates a fundamental geometric bound for normalized transverse profiles. The claimed agreement with experiment is therefore not supported by the model, and the significance of the result as stated cannot be established. The paper also honestly discloses the ad hoc nature of the transverse profiles, which further undercuts the predictive content.

major comments (3)
  1. [Collective mode, Eqs. (8), (17), (18)] The constant D0 ≈ 70 used to reproduce the experimental mode velocity is inadmissible for the normalized transverse profiles defined in the manuscript. Expanding J_αβ(k) = ∫∫ A_α(z) A_β(z') e^{-k|z-z'|} dz dz' for small k gives J_αβ = 1 - k L_αβ + O(k²) with L_αβ = ∫∫ A_α(z) A_β(z') |z-z'| dz dz'. Since the A_α are nonnegative and normalized to unity on a support of width d, one has 0 ≤ L_αβ ≤ d. Consequently (M_so + ka J_ss)/(J_oo + ka) = k(2L_so - L_ss - L_oo + a) + O(k²), so the coefficient D0 in Eq. (17) satisfies D0 ≤ 2 + a/d. With d = 8 nm and a ≈ 0.3 nm this yields D0 ≤ 2.04, not 70. Thus the acoustic velocity in Eq. (18) cannot be achieved with any physical choice of A_s, A_n, A_o consistent with the stated normalization. The quantitative agreement with experiment is enforced by an unphysical fitting constant.
  2. [Model, Eq. (2) and final paragraph before Conclusions] The transverse charge distribution functions A_s(z), A_n(z), and A_o(z) are postulated rather than derived, and the manuscript itself states that their microscopic origin remains an open question. Because the entire acoustic branch in Eq. (18) depends on the structure factor through D0, the central claim is conditional on an unmodeled geometric input. This would be acceptable if the resulting D0 were within the physical bound derived above, but it is not. The mechanism as formulated therefore does not explain the observed high-velocity acoustic mode; it only demonstrates that an acoustic dispersion would follow from an arbitrarily large, physically impossible structure factor.
  3. [Collective mode and numerical analysis, Eq. (18) and following paragraph] The claim of quantitative agreement with experiment is compromised by circularity: D0 is estimated from the experimental zero-temperature velocity and S0 is set by the measured damping at Tc. The subsequent temperature dependence is then carried by the superconducting density Ns(T) and normal density Nn(T), which are not independently fitted but also are not verified against the experimental data in the paper. Thus the agreement in Fig. 1 is not a predictive test of the theory but a consequence of fitting the two structure-factor constants. This circularity is compounded by the inadmissibility of D0 noted above.
minor comments (7)
  1. [Collective mode, Eq. (13)] The 'after algebraic simplification' step leading from Eq. (13) to Eq. (14) is not shown; a derivation or a supplementary note should be provided to allow verification of the approximations used.
  2. [Eqs. (13)-(18)] The notation ω0 is k-dependent (ω0² = 2π e² k N/m), which is not stated explicitly. As written, Eq. (18) can be misread as implying ω_k ∝ sqrt(k); the authors should clarify that the k-dependence of ω0 cancels the sqrt(k) to yield a linear dispersion.
  3. [Collective mode, sentence after Eq. (18)] The statement 'our theory predicts ω_k ∝ d' is only true at fixed k and with N = N_V d; please specify these conditions explicitly.
  4. [Numerical analysis and Fig. 1] The damping Γ in Eq. (18) is independent of k because S0 is a constant, but the experimental data may exhibit k dependence; the paper does not discuss whether this is consistent with the observations.
  5. [Table I] The dirty-limit entry T ≫ Δ for ilde N_s/N is written as 'π Tτ_i Δ²/T²' in a way that is ambiguous; based on Eq. (20) it should presumably be π τ_i Δ²/(2T). Please correct the typesetting.
  6. [References] Reference [10] is 'private communications'; for a quantitative claim about experimental verification, an archival reference or a public data source should be provided.
  7. [Throughout] There are several typographical and grammatical issues, including 'full a full cubic equation' in the paragraph after Eq. (14) and 'Beyond explanation recent measurements' in the Conclusions.

Circularity Check

2 steps flagged · score 6.0 of 10

The acoustic dispersion and absolute frequency/damping in Eq. (18) are set by the fitted constants D0 and S0, so the central quantitative 'prediction' is partially circular; only the scaling laws are independent.

  1. fitted input called prediction [Sec. 'Collective mode', Eqs. (17)-(18)]
    "Using experimental data [7], we estimate a∼0.3nm≪d=8nm, which yields the long-wavelength behavior of the structure factors: Mso +kaJss Joo +ka ≈D0kd, S(k)≈S0, (17) Consequently, the frequency exhibits an acoustic-type dispersion law, ωk =sk=ω0 r D0 Ns N kd, Γ = ω2 0τ 2 Nn N S0. (18)"

    The acoustic dispersion is not derived from the transverse profiles Aα(z): instead, Eq. (17) simply asserts that the geometric structure-factor combination is D0 kd. That linear-in-k replacement is exactly what produces ωk = s k in Eq. (18). Thus the 'prediction' of an acoustic mode is the same statement as the assumed long-wavelength form of the structure factors; the linear dispersion is inserted as an ansatz, not obtained from the microscopic model. The remaining independent content is the scaling with Ns/N and d, which holds for any constant D0.

  2. fitted input called prediction [Sec. 'Collective mode', numerical estimates, p.4]
    "These data give the estimations ω k(T= 0)≈113GHz, what corresponds to s(T= 0)≈1.4·10 9cm/s for D0∼70. The damping Γ(Tc)≈7.4GHz for S0∼7.4·10 −3."

    The constants D0 and S0 are not predicted from the theory; they are fixed by the experimental zero-temperature velocity and the measured damping at Tc. D0≈70 is chosen so that Eq. (18) reproduces s(T=0)≈1.4×10^9 cm/s, and S0≈7.4×10^-3 is chosen to reproduce Γ(Tc)≈7.4 GHz. Consequently, the absolute acoustic frequency and damping quoted as theoretical results are fitted inputs that are then presented as agreement with experiment. Only normalized temperature and thickness dependences are genuine predictions.

full rationale

The algebra from the general dispersion relation (6) through Eq. (15) is internally consistent and not circular. The circularity enters at the decisive step, Eq. (17), where the geometric factor (Mso+kaJss)/(Joo+ka) is replaced by D0 kd without computing D0 from the normalized transverse charge distributions. The value D0≈70 is then fixed by the experimental s(T=0)≈1.4×10^9 cm/s, and S0≈7.4×10^-3 is fixed by Γ(Tc)≈7.4 GHz. Therefore the absolute frequency and damping in Eq. (18) reduce by construction to the very experimental values they claim to explain. The genuinely independent content is the predicted scaling: ωk ∝ sqrt(Ns/N), s ∝ sqrt(d), ωk ∝ d when the bulk 3D density is fixed, and the sqrt(Tc−T) law near Tc; none of these require the specific fitted values of D0 or S0. The admitted 'open question' of the microscopic origin of Aα(z) is a model input rather than a circularity, and the self-citation [12] is peripheral, not load-bearing. Separately, the fitted D0≈70 appears inconsistent with the geometric bound D0≤2+O(a/d) for normalized profiles on a film of thickness d, but the circularity finding rests on the fitted-input structure itself. Overall, the central quantitative prediction is partially circular, while the scaling laws retain independent content: score 6.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The central claim depends on two fitted constants (D0 and S0) and on several assumptions, most importantly the existence of transverse charge-profile mismatch, whose microscopic origin the paper explicitly leaves open. These are loaded into the model rather than derived from the underlying Hamiltonian, so the ledger measures the gap between the stated 'microscopic' ambition and the actual inputs.

free parameters (2)
  • D0 = ~70
    Dimensionless constant controlling the long-wavelength structure factor (Mso+kaJss)/(Joo+ka) approximately D0 kd; matched to the experimental zero-temperature velocity s(0) approximately 1.4e9 cm/s.
  • S0 = ~7.4e-3
    Dimensionless constant controlling the damping structure factor S(k) approximately S0; matched to the experimental damping at Tc, Gamma(Tc) approximately 7.4 GHz.
assumptions (5)
  • domain assumption Continuity equation holds for each component (superfluid and normal) individually
    The paper states this assumption in the Model section after Eq. (2), justified by suppression of pair-breaking in the adiabatic limit omega << Delta; it is needed to derive the dispersion relation.
  • domain assumption Relaxation-time approximation with energy-dependent rate tau_p^{-1} = tau^{-1}|zeta|/epsilon
    Assumed for the collision operator for point-like impurities; used to linearize the kinetic equation and obtain Eq. (3).
  • ad hoc to paper The low-frequency, long-wavelength structure factor (Mso + kaJss)/(Joo+ka) scales as D0 kd with constant D0
    This scaling is introduced after Eq. (16) as an estimate 'using experimental data'; it is central to obtaining linear dispersion, and D0 is fitted, not derived.
  • ad hoc to paper There exists a transverse spatial mismatch between normal and superfluid charge density profiles
    The entire mechanism rests on this mismatch, yet the paper states its microscopic origin remains an open question; the functions A_alpha(z) are never derived.
  • standard math Mattis-Bardeen conductivities describe the dirty limit
    Standard theory of anomalous skin effect in superconductors, used to express polarization operators via conductivities.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Acoustic Plasmon Resonance: Breaking the Anderson Stiffness Paradigm in Quasi-Two-Dimensional Superconducting Films." pith.science (2026). https://pith.science/paper/SJVIBFMT

@misc{pith2026260811610,
  author       = {Pith},
  title        = {Pith review of: Acoustic Plasmon Resonance: Breaking the Anderson Stiffness Paradigm in Quasi-Two-Dimensional Superconducting Films},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SJVIBFMT}},
  note         = {Machine review of arXiv:2608.11610}
}
read the original abstract

Recent experiments on superconducting films have revealed an acoustic plasmon mode that depends critically on the superconducting transition, directly challenging the long-standing Anderson-Higgs paradigm regarding the stiffness of the plasma spectrum in superconductors. In this Letter, we provide a microscopic theoretical framework that explains this behavior and establishes the physical conditions under which classical Anderson-Higgs constraints are bypassed. We demonstrate that in films of finite thickness, the transverse redistribution of normal and superfluid charge densities enables a unique coupling mechanism to electromagnetic radiation - a feature fundamentally absent in the conventional Carlson-Goldman scenario. Our theory predicts an acoustic mode whose dispersion, temperature scaling, and dependence on film thickness are in remarkable agreement with recent experimental observations. By delineating the regime of this acoustic response, we reconcile the observed electromagnetic activity of collective excitations with the fundamental principles of superconductivity.

Figures

Figures reproduced from arXiv: 2608.11610 by the authors.

Figure 1
Figure 1. The frequency decreases with increasing temper [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 1
Figure 1. FIG. 1. Acoustic mode properties as a function of temperature in the dirty regime. (a) Temperature dependence of the [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 12 canonical work pages

  1. [7]

    S. N. Artemenko and A. F. Volkov, Collective excitations in superconductors, Sov. Phys. Usp.22, 295 (1979) [Usp. Fiz. Nauk 128, 3 (1979)]. 6

    work page 1979

  2. [1]

    stiff- ness

    = k2Nn mω2 (−iωτ), where we introduced the normal density asN n = 2N R∞ 0 dζ(ϵ/|ζ|)(−n′ 0). In the dirty limit,l≪ξ 0, the kinetic equation ap- proach is inapplicable. This regime requires a description based on the conductivities of the normal and supercon- ducting fractions under strong scattering. In this regime, we apply the Mattis-Bardeen approach [9]...

    work page 2025

  3. [2]

    P. W. Anderson, Random-Phase Approximation in the Theory of Superconductivity, Phys. Rev.112, 1900, (1958)

    work page 1958

  4. [3]

    Ohashi and S

    Y. Ohashi and S. Takada, On the plasma oscillations in superconductivity, J. Phys. Soc. Jpn.67, 551 (1998)

    work page 1998

  5. [4]

    R. V. Carlson and A. M. Goldman, Superconducting Order-ParameterFluctuationsbelowT c, Phys.Rev.Lett. 31, 880 (1973)

    work page 1973

  6. [5]

    R. V. Carlson and A. M. Goldman, Propagating order- parameter collective modes in superconducting films, Phys. Rev. Lett.34, 11 (1975)

    work page 1975

  7. [6]

    Schmid and G

    A. Schmid and G. Schön, Collective oscillations in a dirty superconductor, Phys. Rev. Lett.34, 941 (1975)

    work page 1975

  8. [8]

    S. A. Andreeva, K. R. Dzhikirba, M. V. Shibalov et.al., Observation of plasma excitations of superconducting electrons in NbN thin films, Phys. Rev. B113, L100503 (2026)

    work page 2026

Show all 13 references
  1. [9]

    A. G. Aronov, Yu. M. Gal’perin, V. L. Gurevich and V. I. Kozub, The Boltzmann-equation description of transport in superconductors, Advances in Physics, 30:4, 539-592 (1981)

    1981

  2. [10]

    D. C. Mattis and J. Bardeen, Theory of the Anoma- lous Skin Effect in Normal and Superconducting Metals, Phys. Rev.111, 412 (1958)

    1958

  3. [11]

    I. V. Kukushkin,private communications

  4. [12]

    L. C. Hebel and C. P. Slichter, Phys. Rev.113, 1504 (1959)

    1959

  5. [13]

    121, 360 (2025)

    R.Z.Vitlina, L.I.Magarill and A.V.Chaplik, JETP Lett. 121, 360 (2025)

    2025

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.