{"id":"07b2942f-62a7-4cf3-a7f4-cb36335be348","arxiv_id":"2411.10852","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The normalized first-harmonic THz nonlinear response of MgB2 rises monotonically as temperature falls, fitting an overdamped π-band amplitude mode with damping rate 0.55 THz, about four times larger than in NbN.","lead":"In the two-gap superconductor MgB2, terahertz two-dimensional coherent spectroscopy shows that the superconducting amplitude mode is only well defined at the lowest temperatures and is far more strongly damped than in the single-gap superconductor NbN. The result points to interband coupling between the two gaps as the likely cause of the enhanced damping.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Attribution of the normalized FH response to an overdamped π-band amplitude mode is unsupported; the single-band fit cannot exclude non-resonant diamagnetic or quasiparticle contributions.","rationale":"The reader's verdict is CONDITIONAL, and my stress test identifies the same load-bearing weakness: the observed normalized FH monotonic increase is not uniquely attributable to the π-band amplitude mode. The paper's quantitative claim (δ=0.55 THz, about 4.5 times NbN) rests entirely on a single-band damped-oscillator fit with no error bars, no goodness-of-fit, and no subtraction of non-resonant contributions. The natural alternative—a diamagnetic (Meissner) nonlinearity that grows with superfluid density—would produce the same monotonic temperature dependence when the drive is below the gap, as it is here. The authors attempt to rule out quasiparticle heating with the pump-probe control, which is a genuine strength, but they do not separate the diamagnetic term analytically. Additionally, their own cited theory (Ref. 32) says the π-band amplitude mode is suppressed by interband coupling, creating an internal tension with the assumption that the π mode dominates the response at 0.63 THz. A two-band calculation including both paramagnetic and diamagnetic responses, using the stated microscopic parameters, would settle whether the amplitude-mode pole is even required to match Fig. 4. Until such a calculation is provided, the correct verdict remains CONDITIONAL: the experimental observation of the monotonic normalized FH is credible, but the interpretation as an overdamped amplitude mode with interband-coupling-induced linewidth is not yet established.","tokens_in":15459,"tokens_out":8092,"duration_ms":85257,"concrete_test":"Compute the third-order nonlinear response of MgB2 using the two-band formalism of Fiore et al. (Ref. 32), including both paramagnetic (amplitude-mode) and diamagnetic (Meissner) terms with the parameters in the SM, and evaluate the normalized FH and TH temperature dependence for Ω/2π=0.63 THz. If the diamagnetic term alone reproduces the monotonic normalized FH without an amplitude-mode pole, or if the amplitude-mode contribution required to match the data has a different δ or negligible weight, the central claim fails. This calculation is directly comparable to Fig. 4.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (1) assumes the FH response is a single damped harmonic oscillator for the π-band amplitude mode, but the paper provides no calculation separating the amplitude-mode (paramagnetic) contribution from the non-resonant diamagnetic and quasiparticle contributions. In NbN (Ref. 20), amplitude-mode dominance was established at resonance (Ω=2Δ); here Ω=0.63 THz lies below 2Δπ≈1 THz for all measured temperatures, and in this off-resonant regime the diamagnetic nonlinearity, which grows with superfluid density as T decreases, naturally yields a monotonic FH increase. The authors' inference from kF l≈13.5 does not carry over off-resonance. Moreover, their cited two-band theory (Ref. 32) states that the π-band amplitude mode is strongly suppressed when interband coupling is finite, which casts doubt on the assumption that the π mode dominates. Consequently, the fitted δ=0.55 THz, obtained without error bars or a two-band model, is not established as the amplitude-mode damping rate, and the interband-coupling conclusion is not uniquely supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports THz two-dimensional coherent spectroscopy (2DCS) measurements on a MgB2 thin film (Tc = 38 K). Using broadband THz pulses, the authors observe a nonlinear signal peaked near 2Δπ ≈ 1 THz at low temperature. Using narrowband pulses at Ω/2π = 0.63 THz, they observe first-harmonic (FH) and third-harmonic (TH) responses that scale as E^6, consistent with third-order nonlinearity. After normalizing the FH intensity by the sixth power of the transmitted drive field, they find a monotonic increase with decreasing temperature, in contrast to the resonant enhancement previously observed in NbN at Ω = 2Δ. The authors interpret this as evidence for a strongly overdamped π-band amplitude mode (δ = 0.55 THz, about 4.5 times larger than in NbN) and attribute the large damping to interband coupling.","tokens_in":15716,"tokens_out":3282,"duration_ms":35485,"significance":"If the interpretation is correct, the paper would provide the first THz 2DCS study of a multigap superconductor and would demonstrate that interband coupling qualitatively changes the amplitude-mode response, a prediction of recent two-band theories (Ref. [32]). The raw data are of good quality: the nonlinear signal is clearly resolved, the E^6 scaling is verified, and the THz pump-probe controls in the SM show that the narrowband drive does not deplete the superconducting condensate. The paper also carefully addresses the normalization issue that affects the temperature dependence of the nonlinear response. However, the central interpretation rests on an unverified single-band model assumption and does not exclude a non-resonant diamagnetic contribution, which would also produce a monotonic FH increase with decreasing temperature. The quantitative claim about the damping rate therefore requires additional analysis or a more cautious framing.","major_comments":[{"comment":"The assignment of the monotonic normalized FH signal to an overdamped π-band amplitude mode is not uniquely supported because the drive frequency Ω/2π = 0.63 THz lies below 2Δπ ≈ 1 THz for all measured temperatures; in this off-resonant regime, a non-resonant diamagnetic nonlinearity that grows with superfluid density would also produce a monotonic increase with decreasing temperature. The argument from Ref. [20] that the amplitude-mode paramagnetic coupling dominates for kF l ≈ 13.5 was established for resonant conditions in NbN and does not automatically carry over to the off-resonant case here. A quantitative estimate separating the paramagnetic (amplitude-mode) and diamagnetic/quasiparticle contributions to the FH response, or a measurement at a drive frequency matched to 2Δπ, is needed to support the attribution.","section":"§4 (normalized FH signal) and Eq. (1)"},{"comment":"The single-band damped-oscillator model of Eq. (1) is applied to the π band of a two-gap superconductor without a derivation or two-band calculation, and the manuscript itself notes that for realistic MgB2 parameters the π-band amplitude mode is strongly suppressed by interband coupling. This internal tension undermines the premise that the π mode dominates the FH response; the authors should either provide a two-band calculation justifying the dominance or explicitly present the fit as a phenomenological characterization rather than as evidence for the amplitude-mode damping rate.","section":"§4, Eq. (1) and discussion of Ref. [32]"},{"comment":"The reported damping rate δ = 0.55 THz is obtained from a fit with free parameters I0 and δ, and an additional ad hoc factor of 0.81 reduction of 2Δπ at 30 kV/cm, with no reported uncertainties, goodness-of-fit, or justification for fixing δ across field strengths while reducing the gap. Given that this fit is the core quantitative evidence for the claim that the damping is about 4.5 times larger than in NbN, the paper must report error bars and fit residuals, and should test whether the data actually constrain δ or whether a similar-quality fit is possible with very different damping values.","section":"§4, Eq. (1) and fitting procedure"}],"minor_comments":[{"comment":"The phrase 'distinct from the single-gap superconductor NbN' would be more accurate as 'distinct from observations in the single-gap superconductor NbN', since the comparison is to one specific film and measurement.","section":"Abstract and §1"},{"comment":"The word 'intrisic' in 'To obtain the intrisic response' appears to be a typo for 'intrinsic'.","section":"§2, paragraph 2"},{"comment":"In the caption of Fig. 1(a), 'pules' should be 'pulses'.","section":"Fig. 1 caption"},{"comment":"The sentence 'For fields larger than 21 kV/cm (Fig. S3.)' should be rephrased to specify that the long-lived component is observed in the THz pump-probe data shown in Fig. S3, and the reference should include the relevant subfigure.","section":"§2, Fig. S3 reference"},{"comment":"The physical meaning of the normalized quantity INL/EA^6 should be stated explicitly: it is proportional to the magnitude squared of the effective third-order susceptibility, assuming a purely electronic third-order process. This would help the reader understand why the sixth power is used.","section":"§3 and §4 (normalization)"},{"comment":"The orange dashed curve in Fig. 2(b) is described as 2Δπ from a two-band BCS calculation, but the parameters of that calculation are only given in the SM; a brief mention in the main text would improve clarity.","section":"§2, Fig. 2 caption and text"}],"recommendation":"major_revision","confidential_remarks":"The raw data and controls are valuable, but the central interpretation currently leans heavily on an unvalidated single-band fit and the paper's own citation to Ref. [32] undercuts the assumption that the π mode dominates. The authors could either strengthen the theoretical analysis or reframe the paper as reporting a phenomenological observation that is consistent with, but not uniquely evidence for, an overdamped amplitude mode. Given the journal's standards, the present version needs the additional analysis before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on the MgB2 THz 2DCS paper. The headline: the normalized FH temperature dependence is a genuine new experimental result, and the paper is candid about the interpretive step, but that step—from a monotonic FH increase to an overdamped pi-band amplitude mode with delta = 0.55 THz—is load-bearing and not yet supported by the model fit alone.\n\nWhat's new and good: this is the first THz 2DCS applied to a multigap superconductor. The broad-band data show a nonlinear peak at 2Delta_pi at low T, and the narrow-band measurements separate FH and TH responses, both scaling as E^6 as expected for third-order nonlinearities. The field-normalization procedure is thoughtful: the raw FH peak near Tc disappears once you divide by the sixth power of the transmitted field, and you get a monotonic increase with decreasing T. The pump-probe controls showing no condensate depletion at 12-30 kV/cm are solid, and the two-band BCS gap solution tracks the measured 2Delta_pi well. These are careful experiments.\n\nThe soft spots are in the interpretation, not the data. Eq. (1) is a single damped oscillator for the pi-band amplitude mode, but the drive frequency (0.63 THz) only matches 2Delta_pi near Tc. Over most of the temperature range the response is off-resonant, and in that regime the non-resonant diamagnetic or quasiparticle nonlinearity, which grows with superfluid density as T drops, will give a monotonic FH increase with no amplitude mode needed. The kF l argument from the NbN work was established at resonance; it does not transfer automatically off-resonance. There is also a tension with the paper's own cited theory: Ref. [32] says the pi-band amplitude mode is strongly suppressed when interband coupling is finite—the authors acknowledge this and then assume the pi mode dominates because it is closer to the drive. That may be true, but it needs a two-band calculation, not an assertion. Finally, delta = 0.55 THz comes from a fit with no error bars, and the 30 kV/cm data require an ad hoc reduction of 2Delta_pi by 0.81 to keep the model alive.\n\nTo be fair, the paper signals these limitations itself ('it is reasonable to assume,' 'likely indicates'). The raw observation stands; the interband-coupling conclusion is a plausible hypothesis, not an established result. This is exactly what peer review is for. I'd send it to a referee who can push for a quantitative estimate of the diamagnetic contribution and a two-band linewidth calculation, or at minimum error bars on delta.","headline":"The normalized FH result is a real new observation, but the overdamped pi-band amplitude mode with delta = 0.55 THz is a plausible hypothesis that needs more than a single-oscillator fit to support it.","tokens_in":16236,"tokens_out":3585,"would_cite":true,"duration_ms":36545,"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":"The multigap superconductor MgB2's amplitude mode is overdamped, with a fitted damping rate of 0.55 THz that the authors attribute to interband coupling.","keywords":["MgB2","multigap superconductor","amplitude mode","terahertz two-dimensional coherent spectroscopy","first-harmonic nonlinear response","interband coupling","overdamped collective mode","superconducting gap dynamics"],"falsifier":"Measure the normalized first-harmonic intensity with narrow-band drives at several frequencies sweeping Ω through 2Δπ(T) at low temperature and fit each temperature scan with the damped-oscillator model leaving δ free; a resonant peak at Ω=2Δπ with δ near 0.12 THz would refute the overdamped/interband-coupling claim, as would the same monotonic FH behavior in a single-gap superconductor.","tokens_in":15245,"feed_emoji":"⚡","tokens_out":9808,"duration_ms":94866,"temperature":0.7,"pith_summary":"This paper tries to establish that the collective amplitude mode of the superconducting order parameter in the multigap superconductor MgB2 is strongly damped, and that the damping comes from coupling between the two superconducting bands. Using terahertz two-dimensional coherent spectroscopy with both broadband and narrowband driving fields, the authors find a nonlinear response at twice the lower gap energy at low temperature, but the normalized first-harmonic signal shows a monotonic increase on cooling rather than the resonant enhancement seen in single-gap NbN. The data are fit with a damped-oscillator model giving a π-band amplitude-mode damping rate δ=0.55 THz, roughly 4.5 times larger than NbN's. If correct, this shows that interband coupling qualitatively changes collective excitations in multigap superconductors and sets a benchmark for how amplitude modes decay in such systems.","feed_headline":"MgB2's superconducting amplitude mode is overdamped","feed_subtitle":"Absent resonance at twice the gap implies a 0.55 THz damping rate from interband coupling.","key_machinery":"The central object is the superconducting amplitude mode—the collective oscillation of the magnitude of the order parameter—in the lower-gap π band of MgB2. The argument is carried by the first-harmonic (FH) intensity of the THz 2DCS nonlinear signal as a function of temperature, normalized by the sixth power of the transmitted drive field to remove screening effects. The fitting machinery is the damped-oscillator resonance model I(ω,T)=I0 Δπ(T)^2 / ((ω+iδ)^2-(2Δπ(T))^2), used with ω=Ω for the FH signal and ω=2Ω for the third-harmonic signal, where δ is the amplitude-mode damping rate. A resonant enhancement is expected when the drive frequency matches twice the gap; its absence in the normalized FH data is read as overdamping.","core_discovery":"The paper reports that in the multigap superconductor MgB2, terahertz two-dimensional coherent spectroscopy with narrow-band drive at Ω/2π=0.63 THz produces first- and third-harmonic nonlinear signals whose normalized first-harmonic intensity rises monotonically as temperature decreases, with no resonant enhancement at the temperature where Ω=2Δπ(T). This contrasts with NbN, where the normalized FH signal peaks resonantly at that matching condition, and is interpreted as evidence that the π-band amplitude mode in MgB2 is strongly overdamped. Fitting the temperature dependence to a single damped oscillator gives damping rate δ=0.55 THz, about 4.5 times the δ=0.12 THz found in NbN. The paper attributes this overdamping to interband coupling between the π and σ bands, which theory says suppresses the π-band amplitude mode. A well-defined amplitude-mode signal is seen only at the lowest temperatures, in a broadband nonlinear peak at 2Δπ≈1 THz.","pith_inferences":["The paper leaves implicit that a direct test would be to tune the narrow-band drive across a range of frequencies near 2Δπ at fixed low temperature; an underdamped mode would show a Lorentzian peak in the normalized FH intensity as a function of drive frequency, whereas pure overdamping would show only a monotonic increase.","If interband coupling is the cause, MgB2 films with modified σ/π coupling, for example through disorder or doping, should show a systematically different damping rate; measuring δ across such samples would separate interband from intrinsic lifetime effects.","The 2D rephasing and non-rephasing spectra, which the paper notes do not show clean diagonal broadening, might still separate amplitude-mode and quasiparticle contributions if analyzed with a two-band model; that is a natural next step not pursued here."],"forward_implications":["If the π-band amplitude mode is overdamped by interband coupling, the first-harmonic channel of THz 2DCS will not show a clean amplitude-mode resonance in MgB2 at Ω≈2Δπ; searches for the mode must rely on the third harmonic or on lower temperatures where the mode is better defined.","The normalization procedure—dividing raw nonlinear intensity by the sixth power of the transmitted drive field—becomes mandatory for comparing superconductors whose optical conductivity changes strongly with temperature; raw FH peaks near Tc are artifacts of screening.","The fitted δ=0.55 THz predicts that the amplitude-mode line is broadened to a degree that should be observable in other MgB2 nonlinear or Raman experiments, providing a cross-check.","Because theory associates the suppression with interband coupling strength, the result supports the view that materials with weaker interband coupling should exhibit sharper amplitude-mode resonances than MgB2."],"supporting_citations":[{"why":"Supplies the NbN comparison, the amplitude-mode interpretation of the resonant first-harmonic enhancement, the damped-oscillator model used for fitting, and the normalization procedure.","marker":"[20]"},{"why":"Theoretical calculation showing that interband coupling suppresses the π-band amplitude mode in MgB2; this is the basis for attributing the large damping to interband coupling.","marker":"[32]"},{"why":"Earlier THz third-harmonic generation experiment in MgB2 whose temperature dependence the present TH signal is compared against.","marker":"[29]"},{"why":"Earlier pump-probe work on MgB2 that identified quasiparticle nonlinear responses and motivates the low-drive-field regime where the condensate is not depleted.","marker":"[27]"},{"why":"Theory on dirty-limit superconductors identifying conditions under which paramagnetic coupling and amplitude modes dominate the THz nonlinear response.","marker":"[31]"},{"why":"Introduces the THz two-dimensional coherent spectroscopy technique used to isolate the third-order nonlinear signal.","marker":"[2]"},{"why":"Supplies transport parameters used to estimate kF l≈13.5, supporting the argument that the large damping is not impurity scattering.","marker":"[41]"}],"fun_headline_variants":["MgB2 amplitude mode overdamped by interband coupling","THz 2DCS: MgB2 amplitude mode only at lowest temps","Multigap MgB2 shows strongly damped π amplitude mode","MgB2's amplitude mode overdamped: 0.55 THz damping rate","Interband coupling kills MgB2 amplitude mode resonance"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the measured first-harmonic signal is dominated by the π-band amplitude mode and behaves like a single damped oscillator; if ordinary quasiparticle or diamagnetic nonlinearity dominates instead, the fitted 0.55 THz damping rate is not the mode's linewidth and the interband-coupling conclusion collapses.","fun_headline_variants_meta":{"raw":{"variants":["MgB2 amplitude mode overdamped by interband coupling","THz 2DCS: MgB2 amplitude mode only at lowest temps","Multigap MgB2 shows strongly damped π amplitude mode","MgB2's amplitude mode overdamped: 0.55 THz damping rate","Interband coupling kills MgB2 amplitude mode resonance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000557,"raw_usage":{"total_tokens":2654,"prompt_tokens":953,"completion_tokens":1701,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":1608}},"tokens_in":569,"tokens_out":1701,"duration_ms":12553,"temperature":1.0,"reasoning_tokens":1608,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:13:36.108876+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the normalized first-harmonic intensity with narrow-band drives at several frequencies sweeping Ω through 2Δπ(T) at low temperature and fit each temperature scan with the damped-oscillator model leaving δ free; a resonant peak at Ω=2Δπ with δ near 0.12 THz would refute the overdamped/interband-coupling claim, as would the same monotonic FH behavior in a single-gap superconductor.","supporting_citations":[{"cited_title":"Murotani and R","cited_arxiv_id":null,"evidence_quote":"Theory on dirty-limit superconductors identifying conditions under which paramagnetic coupling and amplitude modes dominate the THz nonlinear response."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies transport parameters used to estimate kF l≈13.5, supporting the argument that the large damping is not impurity scattering."}],"review_version":1}