{"id":"b144fdeb-79a5-439a-b05e-2eb52cb08701","arxiv_id":"2602.05607","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A gauge-invariant Raman susceptibility is derived for arbitrary multicomponent BdG superconductors, with point-group selection rules and a UTe2 prediction of sharp in-gap peaks from intraband relative modes.","lead":"This paper builds a general microscopic theory of Raman scattering from the internal collective vibrations of multicomponent superconductors. It supplies a computable formula, complete symmetry tables, and predicts sharp in-gap Raman peaks for the candidate superconductor UTe2.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"UTe2 peaks rest on hand-picked order-parameter amplitudes; without a self-consistent gap-equation check, the in-gap Raman resonances in Fig. 4(b) may be artifacts.","rationale":"The reader identifies the missing gap equation and interaction strengths as the weakest assumption. I agree this is the most load-bearing gap, but with an important nuance: the gap equation (A20) can be inverted to define U_i for any nonzero Δ vector, so the issue is not that no interaction exists. The real open questions are whether the inferred U_i are positive, whether the chosen Δ is a self-consistent solution of the nonlinear gap equation, and whether the Gaussian saddle point is stable. The paper's Fig. 4(a) plots absolute eigenvalues, which hides any negative eigenvalues and therefore does not demonstrate stability. The central formula Eq. (24) and the selection-rule derivation appear structurally sound, and the UTe2 section does verify the Ward-Takahashi zero mode, which is independent supporting evidence for the formalism. Thus the concern is addressable and does not overturn the framework; it makes the numerical prediction conditional, matching the reader's verdict.","tokens_in":34913,"tokens_out":19651,"duration_ms":219000,"concrete_test":"From the BdG solution with Δx=Δy=0.04, Δz=0.01, compute U_i = -2β Δ_i / Σ_k Σ_m f(E_m(k))<u_m|U21 Δhat†_i|u_m> using Eq. (A20). (1) Check U_i>0. (2) Insert these U_i into the full nonlinear gap equation and verify that the chosen Δ is a solution. (3) Diagonalize U_eff^{-1}(0,ω=0) without taking absolute values; verify that all nonzero eigenvalues are positive and the only zero is the gauge mode. If any check fails, repeat the Raman calculation with a self-consistently determined Δ; if the in-gap peaks disappear or move outside the quasiparticle gap, the Sec IV prediction is an artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The UTe2 application in Sec IV is the paper's main falsifiable prediction, but it is not tied to a specified microscopic interaction. The amplitudes in Eq. (43) (Δx=Δy=0.04, Δz=0.01) are chosen by hand; no interaction strengths U_i are given and the gap equation (A20) is not solved. The collective-mode kernel U_eff^{-1} in Eq. (A38) is only the correct RPA fluctuation kernel if the expansion in Eq. (13) is about a saddle point of the action (A16); otherwise the linear terms in δΔ do not vanish and the Gaussian integration leading to Eq. (24) is not valid for that state. It is true that for any nonzero Δ vector the gap equation can be inverted to define U_i, so the logical gap is narrower than 'no specified interaction exists.' The missing checks are: (i) the inferred U_i are positive, as assumed in Eq. (3); (ii) the Hessian at ω=0 has no negative eigenvalues; and (iii) the full nonlinear gap equation for those U_i actually returns the chosen Δ. Without these checks, the in-gap peaks in Fig. 4(b) at ω≈0.063 and 0.071 cannot be distinguished from artifacts of an unstable or non-stationary expansion point. The paper only plots |λ_i(ω)| in Fig. 4(a), which folds any negative eigenvalues, so stability is not visible.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Yamazaki and Morimoto present a microscopic functional-integral theory of Raman scattering from collective modes in multicomponent superconductors. For a general BdG Hamiltonian with separable pairing interactions and a scalar potential representing the long-range Coulomb interaction, they integrate out fermions, expand to quadratic order in order-parameter fluctuations, and derive a gauge-invariant Raman susceptibility, Eq. (24)/(A43), expressed through explicit BdG kernel functions. They also derive a group-theoretical selection rule (Γ_{Δ†_α}⊗Γ_γ⊗Γ_{Δ_β} ∋ A_tot) and provide classification tables for all point groups. The formalism is applied to a tight-binding model of UTe2 with an Au odd-parity three-component pairing state; the calculation finds two sharp in-gap Raman resonances below the quasiparticle continuum, attributed to intraband relative modes between Au components, and confirms numerically the Ward–Takahashi zero mode and the equality χ_RR(0,ω)=π_RR(0,ω).","tokens_in":35284,"tokens_out":17560,"duration_ms":182391,"significance":"If the central formula is correct, it is a useful and practical advance: it extends earlier singlet-based formulations to arbitrary BdG Hamiltonians with spin-triplet pairing, gives explicit kernel expressions that can be evaluated numerically from BdG eigenvalues and eigenvectors, and provides a symmetry-based framework for identifying Raman-active Leggett, Bardasis–Schrieffer, clapping, and relative modes. The numerical verification of the Ward–Takahashi identity and the equality of screened and unscreened Raman response at q=0 are strengths. The UTe2 application is a falsifiable prediction, but its credibility currently rests on hand-picked mean-field amplitudes; this needs to be fixed before the prediction can be accepted. The group-theoretical tables also need qualification for multidimensional IRs. With these revisions, the paper would be a valuable contribution.","major_comments":[{"comment":"The UTe2 calculation is not tied to a specified microscopic interaction. The amplitudes in Eq. (43) are chosen by hand, and no Ui are given; the inverse fluctuation kernel U_eff^{-1} in Eq. (A38) is the correct RPA kernel only if the expansion in Eq. (13) is around a saddle point of the action (A16), i.e., if the mean-field gap equation (A20) is satisfied. The gap equation can be inverted to define Ui for any nonzero Δ, so the issue is not that no interaction exists; the missing checks are (i) Ui>0 as assumed in Eq. (3), (ii) no negative eigenvalues of the Hessian at ω=0, and (iii) the full nonlinear gap equation returns the chosen Δ. Fig. 4(a) plots |λ_i(ω)|, so the sign of the eigenvalues is not visible. Without these checks, the sharp in-gap peaks at ω≈0.063 and 0.071 in Fig. 4(b) cannot be distinguished from artifacts of an unstable or non-stationary expansion point. The same applies","section":"Sec. IV.A, Eq. (43), Appendix A 2"},{"comment":"The derivation of the selection rule assumes each gap function transforms as a one-dimensional IR of G (Eq. (27)), but Tables I and II list pairs involving multidimensional IRs (E, T, Eg, Tg). For a component of a multidimensional IR, a generic linear combination does not preserve the full point group, so Eq. (28) does not hold, and the character sum in Eq. (35) is not the correct criterion. The footnote [79] saying higher-dimensional IRs are regarded as 1D IRs of a subgroup does not resolve this, because the tables label full-group IRs and do not specify the subgroup or which component is used. Please either restrict the classification to 1D IRs or derive the correct little-group treatment for multidimensional order parameters. As written, the claim of a classification for 'all point groups' is not fully supported.","section":"Sec. III, Eqs. (27)–(36), Tables I/II"}],"minor_comments":[{"comment":"Fig. 4(a) plots |λ_i(ω)|, so the sign of the inverse-fluctuation eigenvalues is not shown. If the stability checks requested above are added, please plot Re λ_i or otherwise display the sign; otherwise unstable modes are not distinguishable from stable ones.","section":"Sec. IV.B, Fig. 4(a)"},{"comment":"Please specify the convention for the BdG eigenvalues E_m(k) (signed eigenvalues vs. positive quasiparticle energies) and the temperature/filling used in the numerical evaluation of f(E_m). The occupation factors in Eq. (A38) are convention-dependent and the numerics in Fig. 4 are not otherwise reproducible.","section":"Appendix A, Eq. (A38)"},{"comment":"The selection rule is derived in the Ginzburg–Landau regime βΔ≪1, but the UTe2 application is at low temperature. Since the vanishing of symmetry-forbidden vertices is exact, the criterion is presumably general, but the text should state that the linearization in Eq. (33) is only a convenience for the GL analysis and is not required for the symmetry statement.","section":"Sec. III.B"},{"comment":"There are several typos and minor inconsistencies: 'Euqulid' in Appendix C, 'fequency' in the Fig. 5 caption, 'obatin' in Appendix A, and the statement 'ω≈0.083≈E_gap' (Sec. IV.B) is inconsistent with E_gap≈0.074 quoted earlier. Please correct.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper's formal core is solid and the UTe2 application is potentially interesting, but the application is not yet tied to a self-consistent mean-field solution and the group-theory classification overreaches for multidimensional IRs. Both issues are fixable within the manuscript's scope. I would be willing to review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the central formula, Eq. (24)/(A43), is a real step forward. It gives a computable, gauge-invariant Raman response for generic BdG Hamiltonians with separable pairing, and it explicitly covers spin-triplet and time-reversal-symmetry-breaking states, which the earlier spin-singlet formulations didn't. The appendices are detailed and careful, and the numerics confirm the Ward–Takahashi zero mode and the equality χ_RR(0,ω)=π_RR(0,ω). The group-theory classification for all 32 point groups, built on the 'second-order Lifshitz-invariant' idea, is a genuinely useful reference table. That part deserves a serious referee.\n\nWhat's less solid is the UTe2 section. The amplitudes in Eq. (43) are chosen by hand, and no interaction strengths U_i are given. The authors never solve the linearized gap equation (A26) or the full gap equation (A20) to check that their chosen Δ is a stationary point, nor do they check that the inferred U_i are positive or that the RPA Hessian has no negative eigenvalues. Fig. 4(a) plots |λ_i(ω)|, so an instability would be hidden by the absolute value. Without those checks, the in-gap peaks at ω≈0.063 and 0.071 are suggestive, but not established. This is an addressable deficiency—for any chosen Δ one can invert the gap equation to define U_i, so the logical gap is narrow—but the missing positivity and stability checks are exactly what a referee should request.\n\nThe only other quibble: the selection rule is derived in the Ginzburg–Landau regime (βΔ≪1) but the classification is presented as general. In practice the symmetry argument likely goes through without the GL approximation, but the paper should state the exact conditions under which the product rule is proved.\n\nWho is this for? People working on Raman spectroscopy of unconventional superconductors, especially UTe2 and other triplet candidates. It deserves peer review, not desk rejection. My recommendation: send it out, with instructions to either solve the gap equation for specified U_i or explicitly show that the chosen state is a stable saddle point, and to add a short discussion of how far the GL-based selection rule extends.","headline":"A careful, useful extension of Raman susceptibility formalism to spin-triplet and TRSB superconductors, with solid derivations and handy group-theory tables; the UTe2 application needs self-consistent checks before its in-gap peaks can be trusted.","tokens_in":35772,"tokens_out":2515,"would_cite":true,"duration_ms":28235,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A gauge-invariant Raman formula for any multicomponent superconductor is derived, and it predicts sharp in-gap Raman peaks in UTe2 from intraband relative modes, not a Leggett mode.","keywords":["Raman response","collective modes","multicomponent superconductors","Leggett mode","Bardasis-Schrieffer mode","UTe2","gauge invariance","selection rules"],"falsifier":"For the UTe2 model, solve the coupled gap equations (A20) with explicit interaction strengths U_i and check whether Δ_x = Δ_y = 0.04, Δ_z = 0.01 is a stationary point. If it is not, recompute Im χ_RR at the self-consistent amplitudes: the in-gap peaks at ω≈0.063 and 0.071 should disappear or shift if they were artifacts of an inconsistent saddle point. Experimentally, polarization-resolved Raman on clean UTe2 below T_c should either show A_g-symmetric in-gap peaks near the predicted frequencies or not.","tokens_in":1668,"feed_emoji":"🔬","tokens_out":2451,"duration_ms":82238,"temperature":0.7,"pith_summary":"The paper tries to establish a general, gauge-invariant microscopic theory of Raman scattering from collective modes in multicomponent superconductors, covering spin-singlet and triplet pairing, multiband systems, and time-reversal-symmetry-breaking states. The centerpiece is a closed formula for the Raman susceptibility that is directly computable from any BdG Hamiltonian with separable pairing and Gaussian fluctuations, together with a symmetry selection rule and classification tables for all point groups. If correct, it turns Raman spectroscopy into a systematic probe of Leggett, Bardasis-Schrieffer, and clapping modes, and it predicts that fully gapped odd-parity superconductors such as UTe2 can show sharp in-gap Raman peaks from intraband relative modes between pairing components. A sympathetic reader would care because the result is a practical computational recipe and a symmetry dictionary, not just a formal exercise.","feed_headline":"Sharp in-gap Raman peaks predicted for UTe2","feed_subtitle":"A gauge-invariant susceptibility turns any BdG Hamiltonian into a Raman spectrum, exposing Leggett, Bardasis-Schrieffer, and clapping modes.","key_machinery":"The load-bearing object is the gauge-invariant inverse fluctuation kernel U^{-1}_{eff,ϕ}(q,ω), whose inverse U_{eff,ϕ} is the RPA-resummed propagator for the coupled system of order-parameter fluctuations and the Coulomb scalar field. Poles of U_{eff,ϕ} give the collective-mode spectrum, and its coupling vector Q_{R,ϕ} to the Raman source generates the second term of χ_RR. A Hubbard-Stratonovich scalar field implements the long-range Coulomb interaction and the Anderson-Higgs mechanism; the Ward-Takahashi identity makes the overall phase mode a pure gauge zero mode and cancels Coulomb screening at q→0. The selection rule Γ_{Δ†_α} ⊗ Γ_γ ⊗ Γ_{Δ_β} ∋ A_tot is derived in the Ginzburg-Landau regi","core_discovery":"Central claim: For any BdG Hamiltonian with separable pairing and Gaussian fluctuations, the Raman susceptibility can be written χ_RR = 1/4 Φ_RR − 1/8 Q^T_{R,ϕ} U_{eff,ϕ} Q_{R,ϕ}, where the first term is the bare quasiparticle bubble and the second term collects collective-mode contributions via an RPA effective interaction that includes Coulomb through a scalar field. This expression is gauge-invariant, computable from BdG eigenvalues/eigenvectors, and applies to singlet/triplet, single/multiband, TRS/TRSB. The authors also derive the selection rule Γ_Δ†_α ⊗ Γ_γ ⊗ Γ_Δ_β ∋ A_tot and tabulate all point groups. For UTe2 with fully gapped Au pairing, they find sharp in-gap Raman peaks at ω≈0.06","pith_inferences":["[Editorial inference] The same symmetry tables can be used in reverse: a measured polarization-dependent in-gap Raman peak would pin down the symmetry of the two pairing components involved, offering an order-parameter diagnostic for candidate spin-triplet superconductors beyond UTe2.","[Editorial inference] Because the formula cleanly separates quasiparticle and collective contributions and is basis-independent, it should transfer to superconductors coupled to phonons or magnons by enlarging the fluctuation sector—a route the paper flags as future work but does not develop.","[Editorial inference] The prediction of A_g-symmetric in-gap peaks whose intensity tracks spin-orbit coupling provides a concrete experimental target: polarization-resolved Raman on clean UTe2 crystals below T_c. Observing such peaks would also indirectly test whether the assumed Au pairing amplitudes are the correct mean-field state."],"forward_implications":["With a BdG Hamiltonian as input, the full Raman spectrum—quasiparticle background plus collective-mode peaks—follows from explicit kernel expressions without further model-building, enabling computational screening of candidate materials.","The selection-rule tables imply that a given Raman polarization geometry excites only specific pairs of pairing channels; for example, in D4h the B1g geometry can reach collective modes built from (A1g,B1g) or (A1u,B1u) order-parameter pairs.","In UTe2, the two sharp in-gap peaks at ω≈0.063 and 0.071 are predicted to be intraband relative modes among the three Au components, distinguishing this state from conventional two-band Leggett scenarios.","The Raman-active relative Higgs mode at ω≈0.044 is symmetry-allowed but nearly invisible for the standard vertex; increasing the spin-orbit component R_z in the Raman vertex enhances its intensity, so polarization and resonance conditions can tune which mode is observed.","For nodal B3u pairing the same formalism yields only a broad continuum, confirming that a finite quasiparticle gap is required for sharp collective-mode Raman peaks."],"fun_headline_variants":["Raman theory exposes hidden collective modes in superconductors","Gauge-invariant Raman response for any superconductor","Selection rules for Raman-active modes in multiband superconductors","UTe2 in-gap Raman resonances from multicomponent pairing","Theory: Raman reveals Leggett and Bardasis-Schrieffer modes"],"cache_read_input_tokens":36992,"weakest_assumption_plain":"The UTe2 calculation assumes the manually chosen gap sizes Δ_x = Δ_y = 0.04 and Δ_z = 0.01 are a genuine stationary point of some specified pairing interaction, but the paper never gives the interaction strengths or solves the gap equation that would confirm this.","fun_headline_variants_meta":{"raw":{"variants":["Raman theory exposes hidden collective modes in superconductors","Gauge-invariant Raman response for any superconductor","Selection rules for Raman-active modes in multiband superconductors","UTe2 in-gap Raman resonances from multicomponent pairing","Theory: Raman reveals Leggett and Bardasis-Schrieffer modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00081,"raw_usage":{"total_tokens":3436,"prompt_tokens":832,"completion_tokens":2604,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":2528}},"tokens_in":576,"tokens_out":2604,"duration_ms":18934,"temperature":1.0,"reasoning_tokens":2528,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T04:10:56.350867+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the UTe2 model, solve the coupled gap equations (A20) with explicit interaction strengths U_i and check whether Δ_x = Δ_y = 0.04, Δ_z = 0.01 is a stationary point. If it is not, recompute Im χ_RR at the self-consistent amplitudes: the in-gap peaks at ω≈0.063 and 0.071 should disappear or shift if they were artifacts of an inconsistent saddle point. Experimentally, polarization-resolved Raman on clean UTe2 below T_c should either show A_g-symmetric in-gap peaks near the predicted frequencies or not.","supporting_citations":[],"review_version":1}