{"id":"b29cc620-cded-4e32-be6d-e57b5d2c8a34","arxiv_id":"2607.16055","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Terahertz pulses can excite the Higgs mode in superconductors by statically shifting the order parameter—a displacive mechanism that explains the nonlinear first-harmonic enhancement and its measured phase in NbN.","lead":"A theoretical and experimental study shows that intense terahertz pulses excite the Higgs mode in superconductors through a static, non-resonant shift of the superconducting order parameter—the same 'displacive' mechanism that launches coherent phonons in opaque crystals. This explains enhanced nonlinear signals and long-lived pump-probe responses, and gives nonlinear terahertz spectroscopy access to energy relaxation.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Factorization in Eq. (7) is validated only by a self-consistency check at moderate disorder; the experimental phase prediction uses stronger disorder where Anderson-theorem assumptions may fail, and η is an uncomputed input.","rationale":"The reader's weakest assumption identifies exactly the factorization in Eq. (7) and the status of η as the most fragile points. My independent read agrees: the paper's central claim depends on a factorization that is verified only by a self-referential numerical decomposition (Fig. 3d) at moderate disorder (γ/2Δ≈0.5), not at the stronger disorder used for the experimental phase prediction (γ/2Δ≈1.5 in Fig. 4a). The analytical derivation of Eq. (7) assumes Anderson's theorem and a uniform order parameter, while the numerical model (BdG on disordered lattices) can support inhomogeneous pairing; if this mismatch is significant, the 1/ω pole and the resulting phase shift are not robust. Additionally, the paper explicitly states that η is not computed microscopically, so the 1/η divergence and the decay time are inputs, not predictions. These concerns do not refute the plausibility of the displacive mechanism, but they do justify the reader's Conditional verdict: the quantitative claims about phase and relaxation await stronger validation. I recommend no change to the verdict.","tokens_in":18074,"tokens_out":5489,"duration_ms":50003,"concrete_test":"At the parameters of Fig. 4a (V0/t=0.5, c=1, n=0.875, γ/2Δ≈1.5), compute the full BdG order parameter Δ_i for each disorder realization and record its spatial standard deviation. Then separately compute j_NL^Higgs(Ω) from the full calculation and the product [δΔ(0)·(∂χ/∂Δ)] at Ω=2Δ, using δΔ(0) and ∂χ/∂Δ obtained from the same full calculation. If the product deviates from j_NL^Higgs by more than ~20% at this disorder level, or if the spatial variance of Δ_i is comparable to the mean, the Anderson-theorem factorization is not valid in the regime used for the experimental prediction. Alternatively, replace the hand-set η with a microscopically computed inelastic self-energy and check whether the predicted phase shift and decay time survive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central reduction (Eqs. 7–10) factors the nonlinear FH current into ∂χ/∂Δ · δΔ(ω), with δΔ(ω) ∝ (2i Imχ(Ω)/(ω+iη)) H(ω). This relies on (i) Anderson-theorem pairing between exact time-reversed disorder eigenstates and a spatially uniform Δ, and (ii) a sharp 1/(ω+iη) pole whose residue is fixed by Imχ(Ω). Neither (i) nor (ii) is independently established. Fig. 3d tests Eq. (9) by extracting δΔ(0) and ∂χ/∂Δ from the same numerical calculation and multiplying them; this is a consistency check, not an ab initio verification, and it is shown only for γ/2Δ≈0.5. The experimental phase simulation in Fig. 4a uses γ/2Δ=1.5, where the factorization's accuracy is undocumented and where BdG on a disordered lattice can develop spatial Δ fluctuations outside Anderson's theorem. If such inhomogeneous pairing contributes, the pole structure in Eq. (7) changes, and the predicted phase shift of ∂χ/∂Δ loses its quantitative meaning. The paper explicitly concedes that η is not computed ('inelastic scattering processes ... not explicitly computed'), so the 1/η divergence is an input parameter, not a prediction; the 1/η scaling in Fig. 6 merely varies this input. Thus the quantitative content of the headline claim—the phase shift and decay time—is not yet secured.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the dominant nonlinear first-harmonic current in a superconductor driven by a nearly monochromatic THz pulse arises from a displacive mechanism: quasiparticle absorption above 2Δ activates a zero-frequency (static) shift δΔ of the superconducting order parameter, which then modulates the linear optical susceptibility and generates the nonlinear signal j_NL(Ω) ≈ (∂χ/∂Δ) δΔ(0) A(Ω). This factorization is derived analytically from a Lehmann representation of the relevant response functions (Eqs. 5–10), with δΔ(ω) ∼ [2i Imχ(Ω)/(ω+iη)] H(ω) A(Ω) A(ω−Ω). The authors support the derivation with BdG/RPA numerical simulations on disordered lattices (Fig. 3), an analytical Zimmermann-type calculation, and experimental measurements of the temperature-dependent phase of the nonlinear first harmonic in s-wave NbN (Fig. 4). The paper further argues that the same mechanism explains the long-lived pump-probe signal, whose decay is set by the inelastic scattering rate η, and the absence of a sharp FH resonance in d-wave superconductors because low-energy nodal excitations smooth out Imχ(Ω).","tokens_in":18461,"tokens_out":7612,"duration_ms":77068,"significance":"If the central factorization Eq. (9) holds, the paper provides a unifying and testable explanation for several nonlinear THz experiments in superconductors: the FH enhancement at Ω=2Δ in s-wave NbN, the accompanying phase shift of the FH signal, the long-lived 1/η pump-probe response, and the smooth d-wave FH response. The analytical Lehmann construction and the exact-disorder BdG/RPA approach are well matched to the problem, and the paper makes concrete falsifiable predictions (e.g., the absence of a sharp FH resonance in d-wave systems). The analogy to displacive excitation of coherent phonons is conceptually valuable and clearly drawn. These are substantial strengths. The main weaknesses are that the quantitative validation of the factorization is performed in a narrower disorder range than that used for the experimental phase comparison, and that the key dissipation rate η is an input parameter rather than a calculated quantity. Both issues are acknowledged in the text, but they limit the strength of the quantitative claims made in the abstract and conclusions.","major_comments":[{"comment":"The central factorization Eq. (9) is directly tested numerically only for the disorder parameters of Fig. 1b, γ/(2Δ)=0.5, and the analytical Zimmermann comparison in Fig. 3c uses γ/(2Δ)≈0.7. The text states that the agreement improves for γ/(2Δ)<1. However, the phase simulation for the experimental comparison in Fig. 4a uses γ/(2Δ)=1.5, i.e., outside the validated regime and in a parameter range where the Anderson-theorem/uniform-Δ assumptions behind Eqs. (5)–(7) are most questionable. Since the phase shift is the paper's main experimental validation, the authors should either extend the numerical test of Eq. (9) to γ/(2Δ)≈1.5 or provide a separate analytical estimate of the corrections from inhomogeneous pairing at strong disorder. Without this, the experimental agreement in Fig. 4b does not establish the mechanism under the experimental conditions.","section":"FH in s-wave superconductors; Fig. 3 and Fig. 4"},{"comment":"The paper explicitly states that η is not computed: 'any additional source of damping, like the eta parameter appearing in Eq. (2), can only be provided by inelastic scattering processes, which are not explicitly computed in the present treatment.' Equations (10) and (14) nevertheless make the quantitative predictions δΔ(0)∼1/η and a decay time 1/η. Equation (11) is a formal power-absorption identity, not a microscopic calculation of η. Thus the absolute magnitude of the enhancement and the relaxation time are parameterized by an input, not derived. The scaling form is a valid prediction, but the claim that the work 'provides and quantifies a specific microscopic mechanism' overstates what is established. This limitation should be stated in the abstract and conclusions, and the text should distinguish more carefully between the predicted 1/η scaling and the unknown prefactor set by the i","section":"Model; Displacive driving of the Higgs mode; Eq. (11) and Eq. (14)"},{"comment":"The numerical check in Fig. 3d extracts δΔ(0) and ∂χ/∂Δ from the same BdG/RPA calculation that produces j_NL^{Higgs}. Agreement with Eq. (9) is then partly a self-consistency test, because the identity is built from the same response functions. The more independent check is the Zimmermann analytical product in Fig. 3c, but that curve uses phenomenological γ and η values chosen to be 'consistent with numerical simulations.' The authors should explicitly acknowledge the partial circularity and preferably provide at least one direct evaluation of Eq. (9) where the two factors are obtained from independent calculations (e.g., a clean-limit analytical estimate or a separate numerical protocol).","section":"Fig. 3d; validation of Eq. (9)"},{"comment":"The experimental phase data are presented without error bars, statistical information, or a description of the extraction procedure in the main text; the reader is referred to Supplementary [44]. Because this is the only experimental validation of the central mechanism, the main text should at least report uncertainties on φ_exp and state how many samples/measurements are included. The statement 'we report the measured value ... of φ_exp' is too thin for a quantitative claim, especially since the theoretical curve depends on the phenomenological parameters γ/(2Δ)=1.5 and η/(2Δ)=0.05.","section":"Fingerprints of displacive Higgs excitation; Fig. 4"}],"minor_comments":[{"comment":"The labels 'N' and 'A' in the caption are not defined; please spell out 'numerical' and 'analytical' in the caption. Also define 'H' in Fig. 1 caption (Higgs fluctuations).","section":"Fig. 3 caption"},{"comment":"The notation 'iΩ12...' is introduced without a definition. A one-sentence definition of iΩij... would help the reader follow the frequency assignments and the analytical continuations in Eqs. (5)–(8).","section":"Eq. (4) and surrounding text"},{"comment":"Most derivations, including the Lehmann representation (Eq. 5), the Zimmermann approximation, and the pump-probe formula (Eq. 12), are delegated to Supplementary [44] with no equation numbers. Please add specific equation references (e.g., 'SM Eq. (S10)') so that the claims can be checked.","section":"General"},{"comment":"The phrase 'non-resonant static displacement' may be misleading: the displacement is at zero frequency, but it is activated by resonant absorption above 2Δ. Consider rewording to 'zero-frequency (rectified) displacement driven by quasiparticle absorption.'","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper is a strong candidate for publication after revision. The central mechanism is physically compelling and supported by a combination of analytical and numerical work. The main risk is overclaiming quantitative validation: the experimental phase comparison uses a disorder strength where the factorization is not tested, and the η-dependent predictions are parameterized by an uncomputed inelastic rate. I would ask the authors to either extend the numerical check to the experimental regime or clearly state the limitation in the abstract and conclusions, and to provide more detail on the experimental phase extraction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead this one. The core claim is that the enhanced first-harmonic signal at 2Delta in s-wave NbN is not resonant Higgs excitation but a static order-parameter shift driven by quasiparticle absorption, which then modulates the linear conductivity: j_NL = (d chi/d Delta) delta Delta A. The paper derives that factorization from a BCS/RPA calculation with disorder treated exactly, checks it numerically in Fig. 3d and against a Zimmermann approximation in Fig. 3c, and shows the same 1/eta scaling in d-wave. That is genuinely new—previous work attributed the FH enhancement to resonant processes, and this gives a concrete microscopic picture plus a phase prediction that appears to match.\n\nWhat I like: the derivation is analytic, not a fit. The decomposition is physical and testable. The paper is honest about what it cannot compute—it states in the text that inelastic scattering is not explicitly computed, so eta is an input. The d-wave prediction of no sharp FH resonance is a useful discriminant.\n\nWhere I'd push: the numerical validation of the factorization is partly a consistency check—delta Delta(0) and d chi/d Delta are extracted from the same BdG/RPA calculation that produced j_NL, so Fig. 3d confirms the decomposition but not the ingredients independently. More importantly, the experimental phase data in Fig. 4b appear without error bars or a methods paragraph, and the simulation of that data uses gamma/(2Delta)=1.5, while the factorization check in Fig. 3d is shown at gamma/2Delta~0.5. The stress-test note is right that the accuracy of the factorization at stronger disorder, where Anderson-theorem assumptions about homogeneous pairing are less secure, is undocumented. If inhomogeneous pairing contributes there, the pole structure in Eq. (7) changes and the phase shift prediction loses quantitative grounding. Also, eta is never computed from microscopics; the 1/eta divergence and the decay time 1/eta in pump-probe are input dependent, not predictions. That doesn't kill the paper—the qualitative mechanism is plausible and well-supported—but it means the quantitative content of the headline claim (phase shift magnitude, decay time) is not yet nailed.\n\nBottom line: this deserves a serious referee. I'd send it out, with attention on the experimental methods and on whether the factorization survives at the disorder level used in the phase simulation. I'd cite it if I work in this area, and I'd bring it to the reading group—it's a clear, useful mechanism with testable consequences.\n\nRecommendation: send to referees, but the experimental phase measurement needs error bars and the eta identification needs either a microscopic estimate or a clear statement that it's a phenomenological parameter.","headline":"The displacive-Higgs mechanism is a real, analytically grounded new result; the uncomputed inelastic rate and the thin experimental phase data are the soft spots.","tokens_in":18988,"tokens_out":3271,"would_cite":true,"duration_ms":28918,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.25.Gz"],"model":"deepseek-v4-flash","headline":"This paper argues that intense terahertz light above the gap displaces the superconducting order parameter, and that static shift drives the nonlinear first-harmonic signal.","keywords":["Higgs mode","displacive excitation","nonlinear first harmonic","terahertz spectroscopy","superconducting order parameter","disorder","pump-probe","energy relaxation"],"falsifier":"Measure, in the same s-wave NbN film, the pump-probe decay time and the phase of the nonlinear first harmonic while increasing disorder (for example by ion irradiation) at fixed pump frequency: the model predicts the decay time tracks the energy-relaxation rate 1/η and the phase jump stays at 2Δ(T)=Ω; if the decay time instead tracks the momentum-relaxation rate, or the phase jump broadens with disorder, the factorization in Eq. (7) fails.","tokens_in":17941,"feed_emoji":"⚡","tokens_out":7903,"duration_ms":73360,"temperature":0.7,"pith_summary":"The paper identifies a mechanism by which intense terahertz light drives the Higgs (amplitude) mode in superconductors: not by resonantly ringing it at twice the gap, but by absorbing quasiparticles above the gap, which shifts the static order parameter. That static displacement then modulates the linear optical susceptibility and generates a long-lived, rectified nonlinear signal. The authors derive a closed expression for this displacive current and validate it numerically for disordered s-wave and d-wave models and experimentally through the temperature-dependent phase of the nonlinear first harmonic in NbN. If correct, the mechanism unifies several observed terahertz effects — first-harmonic enhancement at 2Δ, a phase jump at that temperature, and slow pump-probe decay — and explains why d-wave superconductors show no sharp first-harmonic resonance. The result matters because it shows nonlinear spectroscopy can measure how light moves a collective order parameter, not just how it excites oscillations.","feed_headline":"THz light bends the superconducting gap to boost nonlinear signal","feed_subtitle":"Static gap shift explains the 2Δ peak, phase jump, and slow pump-probe tail.","key_machinery":"The central object is the static order-parameter displacement δΔ(ω), the zero-frequency Higgs fluctuation induced by difference-frequency mixing of two pump photons. The key identity factors the nonlinear current as a parametric modulation of the linear susceptibility: j_NL(ω+Ω) = (∂χ(Ω)/∂Δ) δΔ(ω) A(Ω), with δΔ(ω) ≈ [2i Im χ(Ω)/(ω+iη)] H(ω) A(Ω) A(ω−Ω). The 1/(ω+iη) pole arises because the pairing operator is diagonal in the quasiparticle basis under Anderson-theorem pairing, so absorbed energy is converted into a static deformation rather than into decaying oscillations. The Higgs propagator H(ω) is finite at ω=0, so the enhancement comes from the divergent susceptibility, not from resonant","core_discovery":"The central claim is that the dominant nonlinear first-harmonic current in a superconductor driven by a nearly monochromatic terahertz pulse has displacive form: j_NL(ω+Ω) ≈ (∂χ(Ω)/∂Δ) δΔ(ω) A(Ω), where δΔ(ω) ≈ [2i Im χ(Ω)/(ω+iη)] H(ω) A(Ω) A(ω−Ω). Here χ is the linear optical susceptibility, Δ the superconducting order parameter, H the Higgs-mode propagator, and η an inelastic scattering rate. The mechanism is activated when the pump frequency exceeds the gap: quasiparticle absorption (Im χ ≠ 0) drives a static shift of the order parameter, which then modulates the susceptibility to emit a rectified signal. The authors show that this process dominates over resonant Higgs driving at 2Ω=2Δ, t","pith_inferences":["This suggests the same absorption-to-static-deformation logic should apply to other collective orders whose order-parameter operator has a finite diagonal projection on quasiparticle eigenstates, making terahertz nonlinear spectroscopy a generic probe of how light moves order parameters.","If the 1/η pole is robust, the pump-probe decay time offers a direct measurement of the energy-relaxation rate, potentially comparable with inelastic rates inferred from transport or other spectroscopies in the same materials.","Adding spatial gap fluctuations beyond Anderson's theorem would test whether the sharp phase jump survives in more strongly disordered films; if it does not, the mechanism may be limited to regimes where pairing is locally translationally symmetric."],"forward_implications":["In dirty s-wave superconductors, the first-harmonic peak at Ω=2Δ is displacive, not resonant: the resonant 2Ω channel is subleading (scaling as 1/√η) compared with the 1/η static displacement.","The phase of the nonlinear first harmonic tracks ∂χ/∂Δ, so it shifts by roughly π/2 as temperature moves 2Δ(T) through the pump frequency; this matches the measured NbN phase behavior.","Pump-probe experiments should show a non-oscillating, rectified signal decaying as exp(−ητ), with η the energy-relaxation rate — distinguishing energy dissipation from the momentum relaxation probed by linear response.","In d-wave superconductors, nodal quasiparticles smooth the absorption edge, removing the sharp onset of δΔ and the ∂χ/∂Δ divergence, so no strong first-harmonic resonance at 2Δ is expected, while the 1/η scaling of the static gap shift remains.","The displacive mechanism assigns long-lived pump-probe signals, previously attributed generically to condensate depletion, to a specific static modulation of the order-parameter amplitude."],"fun_headline_variants":["Displacive gap shift drives Higgs-mode boost in superconductors","THz pulse displaces gap to ignite Higgs nonlinearity","Static gap shift: hidden trigger for THz nonlinear response","Superconductor's gap bend boosts terahertz signal via Higgs"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the static gap shift comes entirely from quasiparticle absorption in a disordered superconductor with Anderson-theorem pairing, and that all inelastic damping can be rolled into one adjustable rate η; if inelastic vertex corrections or spatial gap inhomogeneity alter the 1/(ω+iη) pole, the predicted phase shift and decay time lose their quantitative grounding.","fun_headline_variants_meta":{"raw":{"variants":["Displacive gap shift drives Higgs-mode boost in superconductors","THz pulse displaces gap to ignite Higgs nonlinearity","Static gap shift: hidden trigger for THz nonlinear response","Superconductor's gap bend boosts terahertz signal via Higgs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1233,"prompt_tokens":761,"completion_tokens":472,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":402}},"tokens_in":505,"tokens_out":472,"duration_ms":4680,"temperature":1.0,"reasoning_tokens":402,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T21:29:25.333056+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, in the same s-wave NbN film, the pump-probe decay time and the phase of the nonlinear first harmonic while increasing disorder (for example by ion irradiation) at fixed pump frequency: the model predicts the decay time tracks the energy-relaxation rate 1/η and the phase jump stays at 2Δ(T)=Ω; if the decay time instead tracks the momentum-relaxation rate, or the phase jump broadens with disorder, the factorization in Eq. (7) fails.","supporting_citations":[],"review_version":1}