{"id":"d7db04db-059d-498b-b3d6-01054376a812","arxiv_id":"1908.06125","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In two-band BCS superconductors, post-quench gap oscillations decay as t^-3/2 instead of t^-1/2, with beating from two distinct asymptotic gap values.","lead":"This paper calculates how the superconducting gap evolves after a sudden change of pairing interactions in a two-band superconductor. It finds the gap oscillations damp as t^-3/2, faster than the t^-1/2 decay of single-band superconductors, and show beating between two gap frequencies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Steady-state pseudospin ansatz (Eq. 32) is unproven; the t^-3/2 law rests on its branch-point structure, and numerical validation covers only η=0.8, r=0, weak quenches.","rationale":"After reading the paper in detail, I find the reader's weakest-assumption analysis accurate: the entire route to the central t^-3/2 prediction goes through the steady-state pseudospin ansatz (32). The derivation is internally consistent: the Laplace solution (18)-(25) is a legitimate linearization, the final-value self-consistency reproduces the exact single-band Δ∞, and the branch-point algebra in Appendix B is carefully executed. There is no contradiction with current consensus; the claim that multi-band dephasing accelerates relaxation is plausible and supported by the numerics in the tested regime. The concern is one of unproven generality rather than proven error. The numerical validation is narrow: one value of η, r=0 with a single r=0.5 check, and only weak quenches. Because the small-ε behavior of the steady-state S^x determines the order of the branch point, a failure of the ansatz at other parameters would directly change the exponent, not just prefactors. This is precisely the condition that must hold for the central claim to be fully true. The paper's own statements in Sec. IV C ('our analysis is restricted to weak quenches') and Sec. IV B ('we focus on the case r=0') show the authors are aware of the restrictions, while the abstract and conclusions state the result more broadly. The proposed test (varying η, r, and quench strength and checking both the pointwise steady state and the fitted exponent) would settle whether the concern lands. Given the strong agreement in the tested regime and the recovery of the single-band exact limit, a conditional acceptance remains the right verdict.","tokens_in":23754,"tokens_out":19020,"duration_ms":169802,"concrete_test":"Run the numerical integration of Eqs. (11)-(12) for η=0.5 and η=0.2 (with r=0 and r=0.5), for a weak quench (v_i=0.18→v_f=0.2) and a stronger phase-II quench (v_i=0.15→v_f=0.2), using a fine ε grid and a long-time cutoff. Extract the steady-state S^x_{α,∞}(ε) by time-averaging over several periods after transients and compare pointwise to Eq. (32); also fit the late-time envelope of δ_α(t) over, say, one decade in t to t^{-α}. If α departs from 3/2 by more than 0.1 or the ansatz mismatch exceeds 10% (L2 norm), the central power-law claim is not generic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—the t^-3/2 Landau-damping tail in Eqs. (53)—is derived from the branch-point structure of Γ∞_α(s), built on the ansatz Eq. (32) for the steady-state pseudospin orientation S^x_{α,∞}/Δ_{α,∞} = (Δ̃_{α,f}/Δ̃_{ᾱ,f})/(2√(ε²+Δ²_{α,f})). Substituting this into Eq. (22) gives Φ∞_α(s) = (Δ̃_{α,f}/Δ̃_{ᾱ,f}) Υ(Δ̃_{α,f}, s/(2Δ_{α,∞})) (Eq. 33); the asymptotic analysis in Appendix B then shows Im[zδ(z)] ~ √(y-1) near the branch points, converting the single-band 1/√(y-1) singularity into t^-3/2 decay. The √-behavior depends on the small-ε value of the actual S^x_{α,∞}(ε): if that distribution vanished at ε=0, the branch point would be of higher order and the exponent would change. No derivation of Eq. (32) from the equations of motion or from conserved quantities is given; the authors justify it a posteriori by numerics at η=0.8, r=0 (plus one r=0.5 trace, Fig. 4), for weak quenches 5% and 10% (Fig. 8), and by the single-band η→1 limit. The manuscript contains explicit caveats: after Eq. (40) it notes that agreement with the exact single-band solution does not necessarily imply the distribution function is exact, and Sec. IV C restricts the analysis to weak quenches (phase II) and Sec. IV B to r=0. The abstract's 'wide range of quench amplitudes' overstates this restricted validity. If the ansatz fails for other η, r, or quench strengths, the exponent and prefactors of Eqs. (53) are not generic.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes the collisionless post-quench dynamics of two-band fully gapped s-wave BCS superconductors. The authors generalize the Volkov-Kogan Laplace-space method by linearizing the pseudospin equations around the self-consistently determined non-equilibrium steady state rather than the final equilibrium state. They obtain closed-form expressions for the gap deviations in Laplace space, derive long-time asymptotic formulas for Δ1(t) and Δ2(t), and claim that multi-band coupling changes the algebraic decay from t^-1/2 to t^-3/2 and produces beating at frequencies 2Δ1,∞ and 2|Δ2,∞|. The analytic results are compared with Runge-Kutta solutions for η=0.8, r=0 and weak quenches, and the single-band limit is recovered exactly.","tokens_in":24212,"tokens_out":9428,"duration_ms":89063,"significance":"The claimed t^-3/2 Landau-damping tail is a clear, falsifiable prediction that distinguishes two-band from single-band BCS dynamics, and the steady-state expansion is a useful technique for non-integrable quench problems. The manuscript is carefully written: the Laplace-space chain from Eqs. (15) to (25) is logical, the final-value-theorem determination of Δα,∞ is self-consistent, the single-band limit reproduces the exact Lax-vector result, and the analytic formulas in Appendix D give full prefactors. The main concern is that the central power law rests on the unproven steady-state pseudospin ansatz (32), with numerical support only in a limited parameter region.","major_comments":[{"comment":"The steady-state pseudospin distribution S^x_{α,∞}/Δα,∞ is introduced as an ansatz and is not derived from the equations of motion or from conserved quantities. This distribution enters Φ^∞_α(s) via Eq. (33), and the small-ε behavior of Φ^∞ near the branch points is what converts the single-band 1/√(y-1) singularity into the √(y-1) singularity that yields t^-3/2. The numerical justification in Fig. 8 covers r=0, η=0.8 and two weak quenches, which is not enough to establish the generic small-ε form of S^x_{α,∞}. Please either provide a derivation of Eq. (32) or extract the steady-state pseudospin configuration from long-time numerics for a range of η, r, and quench amplitudes to demonstrate that the assumed √ε branch-point structure is generic.","section":"Sec. IV C, Eq. (32)"},{"comment":"As printed, the asymptotic expressions for Im[1/D(y)] and Im[Υ(·,y)/D(y)] near y→1 contain 1/√(y-1), which would produce a t^-1/2 tail if used in Eq. (49), in contradiction with the text's statement in Sec. IV D that Im[zδ(z)] ~ √(y-1). Direct expansion of D(y) from Eq. (51) around y=1 gives D = (1/κ)Υ(Δ̃2,f,κ) + O(√(y-1)), hence Im[1/D] ∝ √(y-1). Please correct the appendix and check the prefactor signs in the expansions, since this is the load-bearing step for Eqs. (53).","section":"Appendix B, Eq. (B2)"},{"comment":"The claim that the method applies to a 'wide range of quench amplitudes' and that the t^-3/2 damping is 'independent on the quench amplitude' is not supported by the manuscript's own scope, since Sec. IV C restricts the analysis to weak quenches in phase II and Sec. IV B sets r=0, with the r≠0 behavior inferred from a single trace in Fig. 4. Please restate the result as a weak-quench, r=0 finding, with the r dependence as an open question, or provide additional numerical evidence for the claimed generality.","section":"Abstract and Sec. V"}],"minor_comments":[{"comment":"The statement that quenches of r 'are expected to not lead to qualitative changes' is not demonstrated, since the analytic derivation assumes r=0; this should be explicitly marked as an expectation rather than a result.","section":"Sec. II B"},{"comment":"The use of cosh(2Δ1,∞zt) in Eq. (49) is potentially confusing because the integration runs along z=iy; writing the factor as cos(2Δ1,∞yt) would make the inversion transparent.","section":"Eq. (49)"},{"comment":"Please state in the caption which gap component is plotted in Fig. 4B and confirm that the quench parameters are the same as in Fig. 3.","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is publishable in principle, but the main claim currently hinges on the unproven ansatz (32) and the printed Appendix B expansions appear to have the wrong branch-point scaling. I recommend major revision rather than rejection because the single-band limit and the numerical comparisons support the overall framework."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline result is genuine: for two-band s-wave superconductors, the long-time gap oscillation tail is t^-3/2, not the single-band t^-1/2. The mechanism is clean—interband coupling splits the Laplace-space branch point, turning a 1/sqrt singularity into sqrt, and they confirm it numerically with the full prefactors and beating pattern. The methodological step is also worth something: expanding around a self-consistently determined non-equilibrium steady state, rather than the final equilibrium state, and showing it recovers the exact single-band solution. The derivation is explicit and the appendices give all prefactors. Credit where due—this is reproducible, honest work.\n\nThe soft spot is exactly what the stress-test flags. The steady-state pseudospin distribution, Eq. (32), is an unproven ansatz. It satisfies the gap constraint and reduces correctly in the eta->1 single-band limit, but it is not derived from the equations of motion. The entire t^-3/2 branch-point structure depends on the small-epsilon behavior of that ansatz. If the true distribution vanishes differently at epsilon=0, the exponent changes. Numerical validation covers eta=0.8, r=0, two weak quenches (vi/vf=0.95 and 0.9), and one r=0.5 trace. That is enough to make the claim plausible for weak quenches in that regime, but not enough for the abstract's \"wide range of quench amplitudes.\" The paper itself restricts to phase II, so the abstract oversells.\n\nIs this fatal? No. The central physics—two-band dephasing gives faster collisionless relaxation and beating at two gap frequencies—is robust. The precise t^-3/2 exponent is conditional on the ansatz, but the numerics strongly support it for the tested regime. I would want the ansatz tested over a broader eta and r range, or better, derived from a controlled approximation, before treating t^-3/2 as a universal two-band law. For phase II, r=0, it holds up.\n\nThis paper is for anyone working on pump-probe or quench dynamics in multiband superconductors; it gives a concrete observable that distinguishes multi-band pairing. It deserves a serious referee, and I would accept it for peer review with the expectation that the authors either extend the numerics or temper the claim in the abstract.","headline":"A solid analytic result—two-band quench gap tails decay as t^-3/2 instead of t^-1/2—with a real caveat: the key steady-state ansatz is unproven and numerically tested only in a narrow regime.","tokens_in":24757,"tokens_out":1841,"would_cite":true,"duration_ms":19653,"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":"This paper claims that after a quench of the pairing interaction, a two-band BCS superconductor's gap oscillations damp as $t^{-3/2}$ rather than the single-band $t^{-1/2}$, and beat at two frequencies set by the two asymptotic gap values.","keywords":["two-band superconductors","post-quench dynamics","BCS gap oscillations","collisionless damping","power-law decay","pseudospin dynamics","Laplace-space perturbation theory","beating in gap oscillations"],"falsifier":"Solve the full time-dependent two-band BCS equations for a case with unequal densities of states and non-zero intra-band pairing, for example $\\eta=0.5$ and $r=1$, over times $t \\gtrsim 50/\\Delta_{1,\\infty}$, and fit the envelope of $\\Delta_\\alpha(t)-\\Delta_{\\alpha,\\infty}$: a $t^{-3/2}$ decay confirms the claim, whereas a $t^{-1/2}$ or different power law refutes it. An experimental counterpart would be time-resolved terahertz spectroscopy of a two-band superconductor film, where the asymptotic envelope of the coherent gap oscillation should fall as $t^{-3/2}$ rather than $t^{-1/2}$.","tokens_in":23483,"feed_emoji":"⚡","tokens_out":16091,"duration_ms":138454,"temperature":0.7,"pith_summary":"This paper asks what happens to the superconducting gap immediately after the pairing interaction is suddenly changed ('quenched') in a two-band superconductor, before any coupling to the lattice or quasiparticle scattering has had time to act. The paper claims that in the collisionless regime the gap oscillations relax with a power-law envelope $t^{-3/2}$, significantly faster than the well-known $t^{-1/2}$ tail of a single-band BCS superconductor, and that they beat because the two bands settle at two different asymptotic gap values $\\Delta_{1,\\infty}$ and $\\Delta_{2,\\infty}$. This matters because two-band superconductivity is realized in MgB$_2$, iron-based superconductors, and SrTiO$_3$, so the decay law and the beating pattern are concrete signatures of multi-band pairing in ultrafast pump-probe experiments.","feed_headline":"Two-band superconductors damp gap oscillations as t^-3/2, not t^-1/2","feed_subtitle":"Inter-band dephasing accelerates collisionless relaxation and produces two-frequency beating in materials like MgB2 and iron pnictides.","key_machinery":"The argument is carried by the pseudospin representation of the BCS model, in which each Cooper pair is a spin precessing in a self-consistent magnetic field set by the gap, and by a Laplace-space perturbation theory that expands the equations of motion around the long-time non-equilibrium steady state rather than around equilibrium. The central analytic object is the function $\\Upsilon(\\Delta,x)=v_f\\frac{\\sqrt{x^2+1/\\Delta^2}\\,\\arccos(\\sqrt{x^2+1/\\Delta^2})}{\\sqrt{1-(1+x^2)/\\Delta^2}}$, whose branch cuts in the complex-frequency plane control the long-time decay. In the single-band case the Laplace-space response has a $1/\\sqrt{\\epsilon}$ branch point at $s=2i\\Delta_\\infty$, which produces the $t^{-1/2}$ tail; in the two-band case, inter-band coupling converts this to a $\\sqrt{\\epsilon}$ singularity at two branch points, $s=2i\\Delta_{1,\\infty}$ and $s=2i|\\Delta_{2,\\infty}|$, producing the faster $t^{-3/2}$ tail. To close the equations, the paper introduces an ansatz for the steady-state pseudospin distribution, $S^x_{\\alpha,\\infty}/\\Delta_{\\alpha,\\infty} = (\\tilde\\Delta_{\\alpha,f}/\\tilde\\Delta_{\\bar\\alpha,f})/(2\\sqrt{\\varepsilon^2+\\Delta_{\\alpha,f}^2})$, which satisfies the gap constraint and is used with the Laplace final-value theorem to determine $\\Delta_{\\alpha,\\infty}$ self-consistently.","core_discovery":"The central claim is that the long-time gap dynamics after a weak interaction quench is not a sum of two independent single-band responses. For a fully gapped $s$-wave two-band superconductor, the deviation of each gap from its asymptotic value has the form $\\delta_\\alpha(t) \\sim A_\\alpha \\sin(2\\Delta_{1,\\infty}t+\\pi/4)/(\\Delta_{1,\\infty}t)^{3/2} + B_\\alpha \\sin(2|\\Delta_{2,\\infty}|t-\\pi/4)/(|\\Delta_{2,\\infty}|t)^{3/2}$, with analogous terms in the companion band (Eqs. 53a, 53b). The paper argues that the exponent change from $t^{-1/2}$ to $t^{-3/2}$ is generic for two-band systems and independent of quench amplitude within the damped-oscillation regime: inter-band coupling splits the single branch point of the Laplace-space response into two square-root branch points, which shifts the asymptotic decay by one power of time. The same analysis reproduces the exact single-band solution in the $\\eta\\to 1$ limit and matches numerical solutions of the two-band equations for weak quenches with $\\eta=0.8$.","pith_inferences":["A testable extension: terahertz pump-probe studies of a two-band film, e.g., MgB$_2$ or an iron-based superconductor, could extract the envelope exponent from the delay trace; the steeper $t^{-3/2}$ tail would distinguish two-band pairing from the single-band $t^{-1/2}$ response seen in simple metallic films.","In a cold-atom Fermi gas with two components, crossing a pairing resonance should show the same branch-point splitting, providing a clean platform to vary the density-of-states ratio and quench amplitude continuously.","The steady-state pseudospin distribution behind the ansatz of Eq. (32) could be extracted from the same numerical solutions at long times; checking it directly would show whether the $t^{-3/2}$ exponent survives for generic intra-band couplings or only near the pure inter-band limit."],"forward_implications":["Multi-band gap dynamics cannot be decomposed into independent single-band responses: the faster decay is generated by inter-band dephasing of pairs from the two bands.","A clear spectral fingerprint is two oscillation frequencies $2\\Delta_{1,\\infty}$ and $2|\\Delta_{2,\\infty}|$; when the asymptotic gaps are close, the gap trace shows beating at their difference.","The same branch-point splitting argument implies that superconductors or superfluids with more than two bands should also exhibit the $t^{-3/2}$ collisionless decay.","Because the self-consistent Laplace method reproduces the exact single-band solution, it can be applied to quenches in non-integrable fully gapped states such as $s+is$ or $s+id$ pairing, where no exact solution exists.","Within the damped-oscillation regime the exponent is independent of quench amplitude, so observing $t^{-3/2}$ rather than $t^{-1/2}$ would identify multi-band pairing even without knowing the exact quench strength."],"supporting_citations":[{"why":"Provides the Laplace-space perturbative method and the single-band $t^{-1/2}$ collisionless-damping baseline that this paper generalizes to two bands.","marker":"7"},{"why":"Gives the exact single-band BCS solution that the new self-consistent expansion is required to reproduce.","marker":"10"},{"why":"Derives the exact single-band long-time asymptotics and phase boundaries used to locate the damped-oscillation regime.","marker":"11"},{"why":"Classifies the three non-equilibrium phases and shows the single-band $t^{-1/2}$ tail that the two-band result is compared against.","marker":"12"},{"why":"Establishes the single-band non-equilibrium phase diagram and reports a $t^{-3/2}$ tail in the BEC regime by a different mechanism.","marker":"32"},{"why":"Provides the exact single-band dynamical-phase analysis that the new Laplace method reduces to in the $\\eta\\to 1$ limit.","marker":"34"},{"why":"Reports earlier numerical two-band gap dynamics driven by laser pulses, including two-frequency beating and a $t^{-1/2}$ tail that the quench result distinguishes from $t^{-3/2}$.","marker":"39"},{"why":"Numerically studies two-band gap dynamics and the coupling between amplitude and relative-phase modes, motivating the multi-band regime.","marker":"48"},{"why":"Identifies exactly solvable special cases of two-band BCS in which the dynamics reduces to the single-band case, the contrast to the generic model studied here.","marker":"49"}],"fun_headline_variants":["Two-band gaps decay as t^-3/2, not t^-1/2","Inter-band dephasing speeds gap decay to t^-3/2","Two-band superconductors: gap oscillations die faster","t^-3/2 gap decay from inter-band dephasing","Two-band gap beating and t^-3/2 relaxation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire derivation rests on an assumed formula for how the two bands share pairing amplitude after the system has relaxed: the long-time steady-state spin distribution is taken to have the same energy dependence as the final equilibrium distribution, with the two gaps rescaled by the ratio of their asymptotic values; if this distribution is not what generic two-band systems actually reach, the predicted $t^{-3/2}$ exponent would be replaced by something else.","fun_headline_variants_meta":{"raw":{"variants":["Two-band gaps decay as t^-3/2, not t^-1/2","Inter-band dephasing speeds gap decay to t^-3/2","Two-band superconductors: gap oscillations die faster","t^-3/2 gap decay from inter-band dephasing","Two-band gap beating and t^-3/2 relaxation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000207,"raw_usage":{"total_tokens":1498,"prompt_tokens":1139,"completion_tokens":359,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":755,"completion_tokens_details":{"reasoning_tokens":269}},"tokens_in":755,"tokens_out":359,"duration_ms":3950,"temperature":1.0,"reasoning_tokens":269,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:54:55.619836+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full time-dependent two-band BCS equations for a case with unequal densities of states and non-zero intra-band pairing, for example $\\eta=0.5$ and $r=1$, over times $t \\gtrsim 50/\\Delta_{1,\\infty}$, and fit the envelope of $\\Delta_\\alpha(t)-\\Delta_{\\alpha,\\infty}$: a $t^{-3/2}$ decay confirms the claim, whereas a $t^{-1/2}$ or different power law refutes it. An experimental counterpart would be time-resolved terahertz spectroscopy of a two-band superconductor film, where the asymptotic envelope of the coherent gap oscillation should fall as $t^{-3/2}$ rather than $t^{-1/2}$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Laplace-space perturbative method and the single-band $t^{-1/2}$ collisionless-damping baseline that this paper generalizes to two bands."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the exact single-band BCS solution that the new self-consistent expansion is required to reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the single-band non-equilibrium phase diagram and reports a $t^{-3/2}$ tail in the BEC regime by a different mechanism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the exact single-band dynamical-phase analysis that the new Laplace method reduces to in the $\\eta\\to 1$ limit."},{"cited_title":"Akbari , author A","cited_arxiv_id":null,"evidence_quote":"Reports earlier numerical two-band gap dynamics driven by laser pulses, including two-frequency beating and a $t^{-1/2}$ tail that the quench result distinguishes from $t^{-3/2}$."},{"cited_title":"Krull , author N","cited_arxiv_id":null,"evidence_quote":"Numerically studies two-band gap dynamics and the coupling between amplitude and relative-phase modes, motivating the multi-band regime."},{"cited_title":"Dzero , author M","cited_arxiv_id":null,"evidence_quote":"Identifies exactly solvable special cases of two-band BCS in which the dynamics reduces to the single-band case, the contrast to the generic model studied here."}],"review_version":1}