{"id":"98ddf70c-e012-406b-b72f-79e3fb939593","arxiv_id":"2607.06390","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"Charge-imbalance potential relaxation produces a retarded vortex viscosity kernel (1−iωτ_E)⁻¹ that yields a negative effective mass at low frequencies and a reactive pinning-like response at high frequencies.","lead":"The paper derives a vortex equation of motion showing that charge-imbalance potential relaxation causes retardation effects in Abrikosov vortex dynamics at microwave frequencies. This matters for designing superconducting resonators, detectors, and qubits where vortex-induced kinetic inductance affects device performance.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"The derivation of the retarded kernel (1−iωτ_E)⁻¹ relies on u²=ξ²/l²_E≫1 at two distinct steps in Appendix A, and the authors transparently acknowledge this limits quantitative predictions to a regime many experiments do not satisfy.","rationale":"The reader correctly identified the most load-bearing concern: the derivation of the retarded kernel (1−iωτ_E)⁻¹ depends critically on u²≫1 at two stages in Appendix A (dropping ∇|ψ|² terms in A3→A4, and Bessel asymptotics in A14→24), and the authors themselves flag that many experimental systems violate this hierarchy. The paper is internally consistent within its stated domain (gapless TDGL, u²=12 or 5.79), and the negative effective mass result is consistent with prior work [38]. No circular reasoning or algebraic errors were found in checking the derivation. The CONDITIONAL verdict is appropriate: the theoretical result is sound within its domain, but its quantitative applicability to experimentally relevant superconductors (where l_E≳ξ) remains unestablished. The reader's characterization of the concern is accurate and well-targeted. I agree with the verdict and see no reason to adjust it.","tokens_in":15900,"tokens_out":3325,"duration_ms":240922,"concrete_test":"Numerically evaluate the full expression (A14) — before applying the Bessel asymptotic simplification — for u²=1 (i.e., l_E=ξ) and u²=0.1 (l_E≫ξ), across a range of ωτ_E from 0.01 to 100. Fit the resulting frequency-dependent viscosity to the form C/(1−iωτ_eff) and check whether the kernel retains the same functional form with modified parameters, or whether qualitatively new frequency dependence (e.g., additional poles, non-Lorentzian structure) appears. If the kernel structure changes qualitatively for u²≲1, the central claim that the retardation takes the specific form (1−iωτ_E)⁻¹ does not extend beyond the gapless TDGL regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result — the retarded viscosity kernel (1−iωτ_E)⁻¹ in Eq. (24) — is derived through two approximations in Appendix A, both requiring ξ²≫l²_E (i.e., u²≫1): (1) The reduction from Eq. (A3) to the diffusion equation (A4) drops terms containing ∇|ψ|², justified solely by the short-l_E assumption. These terms couple the charge-imbalance potential μ to spatial variations of the order parameter amplitude, which become important when l_E is comparable to or larger than ξ. (2) The evaluation of integral (A14) uses large-argument Bessel asymptotics K₀(z)≈K₁(z)≈(π/2z)^{1/2}e^{-z} for z=ξγ/l_E≫1, which simplifies the ratio in (A14) to unity and yields exactly α_Ohm/(1−iωτ_E). For general u², the full expression (A14) gives a more complex frequency-dependent kernel that does not factorize as (1−iωτ_E)⁻¹. The authors acknowledge in Section IV that 'in many experimentally relevant superconducting systems one can expect the opposite hierarchy l_E≫ξ,' which would both invalidate the derivation and potentially enhance the effect. This is not an internal inconsistency — the derivation is correct within its stated domain (gapless TDGL with u²=12 or 5.79) — but it means the quantitative predictions (Eqs. 34, 36) cannot be directly applied to the experimentally interesting regime without further work. The concern is load-bearing because the kernel structure itself, not just its amplitude, may change when u²~1.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript derives a frequency-dependent vortex equation of motion within the time-dependent Ginzburg-Landau (TDGL) framework for gapless superconductors. The central result is Eq. (27)-(28), in which the Ohmic component of the vortex viscosity acquires a retarded kernel (1−iωτ_E)^{-1} arising from the dynamics of the charge-imbalance potential μ. The authors show that this kernel produces a negative effective mass m_eff = −η_Ohm τ_E at low frequencies (ωτ_E ≪ 1) and a reactive, pinning-like response at high frequencies (ωτ_E ≫ 1). The vortex equation is then used to derive corrections to the microwave impedance and kinetic inductance of superconducting circuits (Eqs. 34, 36). The derivation proceeds via a collective-coordinate perturbation expansion around a static vortex, using the translational zero-mode solvability condition (Eq. 26), with the key technical steps relegated to Appendix A.","tokens_in":16236,"tokens_out":1310,"duration_ms":291919,"significance":"The paper addresses a well-defined problem in vortex dynamics — the intrinsic retardation of the Ohmic viscosity channel — and provides a self-contained derivation from standard TDGL equations. The connection between the charge-imbalance potential relaxation and the (1−iωτ_E)^{-1} kernel is physically transparent and the analogy with the two-frequency-dependent penetration depth (Eq. 5) is instructive. The prediction of a negative effective mass consistent with prior results (Ref. [38]) and the renormalization of the depinning frequency at high frequencies are concrete, falsifiable results. The derivation is internally consistent within its stated domain (gapless TDGL with u² ≫ 1).","major_comments":[{"comment":"§IV, Eqs. (34), (36), and (37): The quantitative predictions for impedance and kinetic inductance depend on the kernel (1−iωτ_E)^{-1}, whose derivation in Appendix A requires ξ² ≫ l_E² (i.e., u² ≫ 1) at two distinct steps: (1) the reduction of Eq. (A3) to the diffusion equation (A4) by dropping terms containing ∇|ψ|², and (2) the evaluation of the integral (A14) using large-argument Bessel asymptotics K₀(z)≈K₁(z)≈(π/2z)^{1/2}e^{-z} for z = ξγ/l_E ≫ 1, which forces the ratio in (A14) to unity. The authors transparently acknowledge (Section IV, final paragraphs) that 'in many experimentally relevant superconducting systems one can expect the opposite hierarchy l_E ≫ ξ,' which would invalidate both steps. This is load-bearing because for general u² ~ 1, the full expression (A14) yields a more complex frequency-dependent kernel that does not factorize as (1−iωτ_E)^{-1}, and the regime where ","section":null}],"minor_comments":[{"comment":"Eq. (5): The heuristic argument relating η_Ohm(ω) to λ²_eff(ω) is appealing but the logical flow is slightly circular — the proportionality η_Ohm ∝ l²_E/ξ² is stated as a known result and then used to motivate the (1−iωτ_E)^{-1} factor before the TDGL derivation. A brief clarifying sentence noting that this is a physical motivation, not a derivation, would help the reader.","section":null},{"comment":"Eq. (29): The viscosity expressions contain α_Rel and α_Ohm, which are stated as numerical constants (α_Rel ≈ 0.279, α_Ohm ≈ 0.159 for u² = 12). It would help to state explicitly whether these specific numerical values are used in the subsequent impedance formulas or whether the general α_Ohm ≈ 2u^{-2} scaling is employed.","section":null},{"comment":"Eq. (30): The depinning frequency for a vortex pinned by a cavity is given as ω_k τ_E = 2αu²/(3f'(0)²a²), but the notation f'(0) is not immediately defined at the point of use (it is stated as f'(0) ≈ 0.583/ξ only parenthetically). A forward reference or brief definition would improve readability.","section":null},{"comment":"§IV, Eq. (34): The simplifying approximation j̃_{s,∞} = j̃_tr is noted as a limitation. The authors should briefly comment on the magnitude of the error introduced by neglecting the normal component of the transport current at finite frequency, particularly for the high-frequency regime ωτ_E ≫ 1.","section":null},{"comment":"Appendix A, Eq. (A2): The transition from Eq. (A1) to (A2) involves using the continuity equation (9) together with Eq. (11), but the algebraic steps are compressed. A brief intermediate expression would aid verification.","section":null},{"comment":"References: The paper by Kogan and Nakagawa (Ref. [43]) is cited for the two-fluid penetration depth, but the specific result λ²_eff(ω) = λ²(1−iωτ_E)^{-1} should be more precisely attributed or derived, as this is a standard result but the specific form with τ_E may differ from the conventional two-fluid expression.","section":null},{"comment":"Typo in §IV: 'can also affect affect the dynamics' — duplicated 'affect'.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The stress-test concern about the u² ≫ 1 assumption is valid and well-articulated by the authors themselves. The key question for the editor is whether the paper's contribution is primarily the derivation within the gapless TDGL domain (which is correct and self-contained) or the applicability to experiments (which is limited). I lean toward the former: the derivation is sound within its stated domain, the limitations are transparently acknowledged, and the extension to u² ~ 1 is explicitly identified as future work. This warrants minor revision rather than major revision, provided the authors add a brief discussion of what changes qualitatively when u² ~ 1 (e.g., whether the kernel structure itself changes or only its amplitude)."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful reading and for identifying the precise technical steps in Appendix A where the u² ≫ 1 assumption enters. The referee's assessment of our central results and their physical significance is appreciated. Below we address the major comment point by point.","responses":[{"response":"The referee is correct that the factorized kernel (1−iωτ_E)^{-1} relies on the u² ≫ 1 limit at the two steps identified, and we agree that this should be stated more explicitly in the main text rather than being implicit from the Appendix. We note the following points in response. First, the manuscript already presents the full expression (A14) before taking the large-u limit, so the general-case result is available to the reader. The ratio of Bessel functions in (A14) is the exact kernel within the Bardeen–Stephen approximation for the vortex core, and it reduces to (1−iωτ_E)^{-1} only when ξ ≫ l_E. For u² ~ 1, this ratio deviates from the simple Lorentzian form, and the referee is right that the factorization no longer holds. Second, we agree that the manuscript should more clearly delineate which results are exact within the stated TDGL framework and which additionally require u² ≫ 1. Specifically, Eqs. (34), (36), and (37) inherit the u² ≫ 1 assumption from Eq. (24), and this dependency should be flagged at the point where those equations are introduced, not only in the final paragraphs of Section IV. Third, regarding the experimentally relevant regime l_E ≫ ξ: we already note in the final paragraphs of Section IV that in this regime the spatial distributions of E and μ extend beyond the vortex core, which can enhance both η_Ohm and τ_E, and that a quantitative treatment requires solving for μ and E without the short-l_E approximation. We believe the physical mechanism — retardation of the Ohmic viscosity channel due to charge-imbalance relaxation — is robust beyond the strict u² ≫ 1 limit, since it originates from the finite relaxation time τ_E of the gauge-invariant potential μ, which is a general feature of the TDGL model. However, the specific Lorentzian form of","revision_made":"no","referee_comment":"§IV, Eqs. (34), (36), (37): The quantitative predictions depend on the kernel (1−iωτ_E)^{-1}, whose derivation requires ξ² ≫ l_E² at two steps in Appendix A: (1) dropping ∇|ψ|² terms in going from (A3) to (A4), and (2) using large-argument Bessel asymptotics in (A14). The authors acknowledge that experimentally l_E ≫ ξ may hold, invalidating both steps. For general u² ~ 1, the full (A14) yields a more complex kernel that does not factorize as (1−iωτ_E)^{-1}."}],"tokens_in":15501,"tokens_out":1218,"duration_ms":78492,"standing_objections":["The referee's comment appears to be truncated mid-sentence ('and the regime where'), so we cannot fully determine whether there is an additional specific request beyond what we have addressed above. If the referee intended to ask for a quantitative treatment of the u² ~ 1 or l_E ≫ ξ regime, we note that this would require numerical solution of Eq. (A9) with the full Bessel-function ratio in (A14) and, for l_E ≫ ξ, abandoning the Bardeen–Stephen core approximation altogether. This is a substantial extension beyond the scope of the present work, which is explicitly framed within the gapless TDGL model with u² ≫ 1."]},"desk_editor":{"model":"glm-5.2","letter":"The main thing to know: this paper derives the retarded vortex viscosity kernel (1−iωτ_E)⁻¹ from TDGL charge-imbalance dynamics, replacing the phenomenological vortex mass with a microscopically grounded retardation mechanism. The derivation is legitimate and self-contained. Starting from standard TDGL equations, the authors use a collective-coordinate perturbation expansion around a static vortex, apply the solvability condition from the translational zero mode, and arrive at the equation of motion η(ω)ṽ = F̃_L with η(ω) = η_Rel + η_Ohm/(1−iωτ_E) − η_tot ω_k/(iω). The low-frequency expansion gives a negative effective mass m_eff = −η_Ohm τ_E, consistent with Kupriyanov and Likharev's 1975 numerical result. The high-frequency limit gives a reactive pinning-like term that renormalizes the depinning frequency. Both limits are physically sensible and the algebra in Appendix A checks out. The connection to kinetic inductance corrections (Eqs. 34, 36) is straightforward once you have the equation of motion. The two-fluid argument in the introduction (Eq. 5) is a nice physical motivation — the same (1−iωτ_E)⁻¹ factor appears in the effective penetration depth, and the paper makes the connection to vortex viscosity explicit. The soft spot is real but not fatal. The derivation of the kernel (1−iωτ_E)⁻¹ in Appendix A requires u² = ξ²/l_E² ≫ 1 at two steps: dropping ∇|ψ|² terms to get the diffusion equation (A4) from (A3), and using large-argument Bessel asymptotics to simplify (A14) to (A15). The authors are transparent about this — they note in Section IV that many experimentally relevant superconductors have l_E ≫ ξ, which would both invalidate the derivation and potentially enhance the effect. The concern is load-bearing in the sense that the kernel structure itself, not just its amplitude, may change when u² ~ 1. The full expression (A14) before the asymptotic simplification is more complex and does not obviously factorize as (1−iωτ_E)⁻¹. So the quantitative predictions for kinetic inductance and impedance cannot be directly applied to the experimentally interesting regime without further work. That said, the derivation is correct within its stated domain (gapless TDGL with u² = 12 or 5.79), and the qualitative physics — retardation from charge-imbalance relaxation producing a negative mass at low frequencies and a pinning-like response at high frequencies — is likely robust even if the exact kernel form changes. This paper is for theorists and experimentalists working on superconducting microwave devices who need a microscopic picture for frequency-dependent vortex dynamics. It deserves a serious referee who can check the Appendix A algebra and assess whether the u² ≫ 1 limitation undermines the paper's claims for its target audience.","headline":"Clean TDGL derivation of retarded vortex viscosity kernel (1−iωτ_E)⁻¹, limited to the l_E≪ξ regime where many experiments sit in the opposite limit","tokens_in":16981,"tokens_out":713,"would_cite":true,"duration_ms":86964,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Charge imbalance gives vortices a memory, reshaping superconducting circuits","keywords":[],"falsifier":"Measure the microwave impedance of a superconducting film with controlled vortex density as a function of frequency and temperature near T_c. If the retarded kernel (1 - i*omega*tau_E)^{-1} correctly describes the vortex dynamics, the impedance should show a crossover from dissipative (flux-flow) to reactive (pinning-like) behavior at omega ~ 1/tau_E, with the crossover frequency scaling as (1 - T/T_c). Absence of this crossover, or a crossover frequency that does not diverge near T_c, would contradict the predicted mechanism.","tokens_in":16200,"feed_emoji":"🌀","tokens_out":1454,"duration_ms":115317,"temperature":0.7,"pith_summary":"The paper derives a modified equation of motion for Abrikosov vortices in a type-II superconductor driven by microwave currents, starting from the time-dependent Ginzburg-Landau (TDGL) equations. The key physical object is the charge-imbalance potential — a gauge-invariant scalar potential that arises inside and around a moving vortex core as supercurrent converts to normal current. The authors show that this potential does not relax instantaneously; its finite relaxation time tau_E introduces a retarded kernel (1 - i*omega*tau_E)^{-1} into the Ohmic part of the vortex viscosity. This retardation has two distinct consequences depending on frequency. At low frequencies (omega*tau_E << 1), expanding the kernel produces a negative effective vortex mass m_eff = -eta_Ohm * tau_E, consistent with earlier TDGL results. At high frequencies (omega*tau_E >> 1), the retarded Ohmic contribution becomes purely reactive and acts like an additional pinning force, renormalizing the depinning frequency while the remaining viscosity becomes purely relaxational. The authors then connect this modified vortex dynamics to the complex impedance and kinetic inductance of superconducting circuits, deriving explicit expressions for the vortex contribution to impedance in both frequency regimes. The correction to kinetic inductance from the charge-imbalance mechanism is proportional to the ratio of Ohmic to total viscosity and to the vortex density, and it diverges near T_c because tau_E diverges there.","feed_headline":"Charge imbalance gives vortices a memory, reshaping superconducting circuits","feed_subtitle":"A retarded kernel in vortex viscosity, born from charge-imbalance relaxation, creates negative mass at low frequencies and pinning-like rig1","key_machinery":"TDGL equations (Eqs. 6-7); collective-coordinate perturbation theory for vortex displacement (Eq. 12); solvability condition from translational zero mode (Eq. 26); diffusion equation for charge-imbalance potential mu (Eq. A5) derived under assumption xi^2 >> l_E^2; Bardeen-Stephen step-function approximation for vortex core; modified Bessel function K_1 matching at core boundary (Eq. A12-A13); Bessel asymptotics yielding the retarded kernel (Eq. A15); two-fluid model relation lambda_eff^2(omega) = lambda^2 / (1 - i*omega*tau_E)","core_discovery":"The central result is the vortex equation of motion, Eq. (27)-(28): eta(omega) * v_tilde = F_tilde_L, where the frequency-dependent viscosity eta(omega) = eta_Rel + eta_Ohm / (1 - i*omega*tau_E) - eta_tot * omega_k / (i*omega). The retarded factor (1 - i*omega*tau_E)^{-1} multiplying the Ohmic viscosity is derived from the dynamics of the charge-imbalance potential mu, which obeys a diffusion-type equation (A5) with relaxation time tau_E = l_E^2 / D. This single kernel produces both a negative effective mass at low frequencies and a pinning-like reactive response at high frequencies, unifying two previously separate phenomena under one mechanism: the finite relaxation time of the charge-imb1","pith_inferences":["The regime where the authors predict the largest enhancement of the kinetic-inductance correction (l_E >> xi) is precisely the regime where their derivation breaks down (they assume xi >> l_E). An experimental test in materials with large l_E/xi ratio would simultaneously probe the physics and the limits of the TDGL gapless approximation.","If the charge-imbalance relaxation time tau_E depends on the order-parameter amplitude between vortices (as the authors hint for fields near H_c2), then in a dense vortex lattice the retardation effects could appear at even lower frequencies than the isolated-vortex estimate suggests, making them relevant for devices operating in moderate magnetic fields.","The negative effective mass m_eff = -eta_Ohm * tau_E implies that near T_c, the vortex response has an anti-inertial character: the vortex accelerates opposite to the applied force direction at low frequencies, which could produce unusual phase shifts in resonator S-parameters that are distinguishable from conventional pinning.","A direct experimental signature would be measuring the frequency-dependent crossover in the vortex impedance: below 1/tau_E the response is dissipative with a negative-mass correction, above 1/tau_E it becomes reactive with a renormalized depinning frequency. Mapping this crossover as a function of temperature (since tau_E ~ 1/(1 - T/T_c)) would test the predicted divergence."],"forward_implications":["Superconducting microwave resonators operating near T_c should exhibit a frequency-dependent kinetic inductance correction from vortex motion that cannot be captured by a simple mass-plus-viscosity model; the full retarded kernel must be used.","In the high-frequency regime omega*tau_E >> 1, trapped vortices contribute a reactive pinning-like term to impedance even without extrinsic pinning centers, renormalizing the depinning frequency upward by a factor involving eta_Ohm / (eta_Rel * tau_E).","The divergence of tau_E near T_c means the crossover frequency 1/tau_E shifts downward, making retardation effects accessible at lower microwave frequencies for devices operated close to the superconducting transition.","For multiband or nematic superconductors, the same retardation kernel (1 - i*omega*tau)^{-1} should appear with a modified relaxation time tau, potentially offering an additional tuning knob for kinetic inductance in complex order-parameter systems."],"fun_headline_variants":["Charge imbalance unifies negative mass and pinning in vortex dynamics","Retarded vortex dynamics from charge imbalance alters kinetic inductance","Frequency-dependent vortex viscosity from charge-imbalance relaxation","Charge-imbalance relaxation drives high-frequency vortex memory","Charge imbalance creates negative mass and pinning in superconductors"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The derivation requires the charge-imbalance length l_E to be much smaller than the coherence length xi, which enters when simplifying the equation for the charge-imbalance potential to a diffusion equation and when using Bessel function asymptotics. The authors note that many experimentally relevant superconductors likely have the opposite hierarchy (l_E >> xi), and the regime where the predicted effects are strongest is exactly where the derivation is not valid.","fun_headline_variants_meta":{"raw":{"variants":["Charge imbalance unifies negative mass and pinning in vortex dynamics","Retarded vortex dynamics from charge imbalance alters kinetic inductance","Frequency-dependent vortex viscosity from charge-imbalance relaxation","Charge-imbalance relaxation drives high-frequency vortex memory","Charge imbalance creates negative mass and pinning in superconductors"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1410,"prompt_tokens":497,"completion_tokens":913,"prompt_tokens_details":null},"tokens_in":497,"tokens_out":913,"duration_ms":40754,"temperature":1.0,"reasoning_tokens":887,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-08T07:16:45.809578+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"Measure the microwave impedance of a superconducting film with controlled vortex density as a function of frequency and temperature near T_c. If the retarded kernel (1 - i*omega*tau_E)^{-1} correctly describes the vortex dynamics, the impedance should show a crossover from dissipative (flux-flow) to reactive (pinning-like) behavior at omega ~ 1/tau_E, with the crossover frequency scaling as (1 - T/T_c). Absence of this crossover, or a crossover frequency that does not diverge near T_c, would contradict the predicted mechanism.","supporting_citations":[],"review_version":1}