{"id":"889c880c-abee-4760-862d-42e0b04322be","arxiv_id":"1908.08559","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"In spin- and mass-imbalanced Fermi superfluids, Anderson-Bogolyubov sound modes acquire Landau damping at zero temperature once the Fermi-surface mismatch exceeds a threshold set by the mass ratio and the superfluid gap.","lead":"This paper studies how sound-like collective oscillations in a superfluid Fermi gas die out when the two spin species have different masses or densities. It derives a simple inequality that predicts when these Anderson-Bogolyubov modes are damped even at absolute zero, a regime relevant for ultracold atom experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (19) double-counts spin branches, but the central damping threshold (21) survives.","rationale":"Reader's weakest assumption — the unproved dominance of σ=σ′ — is the right place to look, but a direct analysis of Eq. (13) shows the other terms cannot produce a |q_m|/|q| term: σ≠σ′ denominators are finite at q=0 for Δ>0, and the σ=σ′ linear frequency term is odd in the angle and integrates to zero. The central threshold (21) therefore rests on solid ground and is also confirmed by the numerical activation point in Figs. 5 and 6. The concrete defect I found is a factor-of-two (or spin-asymmetry) overcount in Eq. (19), which affects the quoted analytical damping coefficient but not the existence condition. The paper should correct this and soften the 'full agreement' statement, but the conditional verdict stands as is; no further verdict movement is warranted.","tokens_in":13676,"tokens_out":21042,"duration_ms":207932,"concrete_test":"Set r=1 (ζ=0), choose μ, Δ, h with h>Δ, and evaluate Eq. (16) directly using δ(E^+_ε) with the two roots ε_{1,2}=μ±√(h^2−Δ^2), discarding σ=−. Compare the resulting γ^{-1} with Eq. (19), which sums both σ at each root. If the direct value is half of Eq. (19), the double-counting is confirmed. Then repeat the Fig. 5 activation scan with the corrected single-spin sum to verify that the activation value of h is unchanged.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing step in deriving the threshold is the claim that only σ=σ′ terms in Eq. (13) generate the |q_m|/|q| nonlocality. This claim is actually justified: for σ≠σ′ the energy denominator E^{σ,σ′}_{k,q} is gapped by 2E_k > 0 for Δ > 0, and the σ=σ′ linear-iq_m contribution is odd in cosθ and vanishes after angular integration. The genuine weakness is in Eq. (19), which sums ∑_σ L^{-1}_σ(ε_i) at each root ε_i of the squared equation. For Δ > 0, E^σ_ε = 0 and E^{−σ}_ε = 0 cannot both hold at the same ε, and for h−ζμ ≥ 0 with ε_i ≥ 0 one finds h−ζε_i > 0, so only σ=+ contributes to δ(E^σ_ε). Summing over both spins at each root therefore overcounts the damping coefficient, by a factor of 2 in the mass-balanced (ζ=0) limit. This error does not affect the existence condition (21) — the non-negativity of ε_i is independent of the prefactor — but it invalidates the literal analytical expression for γ^{-1} and weakens the paper's claim of full quantitative agreement for damping rates.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Anderson-Bogolyubov (Goldstone) phonon in a two-component Fermi superfluid with spin and mass imbalance, working in the Gaussian pair fluctuation (GPF) approximation. The authors isolate the temporally nonlocal |q_m|/|q| contribution to the inverse pair fluctuation propagator and relate it to Landau damping of the collective mode. Their central result is Eq. (21), a necessary condition for damping to be active at T=0, expressed in terms of the mismatch parameter h−ζμ, the mass-imbalance parameter ζ, the gap Δ, and the chemical potential μ. They complement the analytic derivation with a numerical solution of the complex pole of the GPF propagator for mass ratios r=1 and r=6.67, finding that damping activates at the predicted value of h in both a continuous and a first-order transition regime.","tokens_in":13988,"tokens_out":12085,"duration_ms":113719,"significance":"If the threshold condition (21) holds, the paper provides a simple, falsifiable criterion for zero-temperature Landau damping of the Anderson-Bogolyubov mode in imbalanced Fermi superfluids. This is a useful analytical result, especially since the numerics indicate that damping sets in for h substantially below the superfluid-normal transition value hc. The derivation is largely transparent, the numerical check is consistent with the threshold, and no free parameters are fitted apart from the UV cutoff. The central claim is defensible, but the analytical expression for the damping coefficient γ^{-1} in Eq. (19) contains a spin-summation error that overcounts the contribution; this does not affect the threshold but invalidates the literal rate formula and weakens the paper's quantitative claims.","major_comments":[{"comment":"Equation (19) sums Σ_σ L^{-1}_σ(ε_i) at each root ε_i of the squared equation, but after squaring, only one spin branch satisfies the original equation E^σ_ε=0. For h−ζμ≥0, the case analyzed in Sec. III, both roots correspond to σ=+; the σ=− delta function vanishes identically. Summing over both spins therefore overcounts γ^{-1}, giving a factor of 2 in the mass-balanced limit ζ=0. This does not affect the threshold condition (21), which depends only on ε_i≥0, but it invalidates the literal expression for the damping coefficient and weakens the claim in Sec. V of 'full agreement' with numerically computed damping rates. The authors should replace the σ-sum by the single branch σ=sign(h−ζμ) and re-examine the quantitative comparison.","section":"Eq. (19)"},{"comment":"The concluding statement that analytical predictions are 'in full agreement with damping rates obtained numerically' is not supported by the presented analysis. The numerical study (Figs. 4–6) compares the activation value of h, not the magnitude of γ^{-1} from Eq. (19). Since Eq. (19) is erroneous (see the previous comment), the claim should be revised to state agreement only for the threshold condition, or the comparison should be repeated with the corrected expression.","section":"Sec. V"}],"minor_comments":[{"comment":"The text contains the typo 'Fiq. 1', which should read 'Fig. 1'.","section":"Sec. III, after Eq. (21)"},{"comment":"The function f^{σ,σ'}_{k,q} is used in Eq. (13) before it is defined in Appendix A. Please define it in the main text or add an explicit pointer to the appendix.","section":"Sec. III, Eq. (13)"},{"comment":"The numerical parameters include a hard UV cutoff Λ=10. A brief statement about the sensitivity of the threshold condition and of the reported damping rates to the cutoff would be helpful.","section":"Sec. IV"},{"comment":"The extremely small damping rates (of order 10^{-10}–10^{-11}) are quoted without an estimate of numerical precision. A short note on how the analytic continuation and root-finding errors were controlled would be useful.","section":"Sec. IV, Figs. 3 and 4"}],"recommendation":"major_revision","confidential_remarks":"The paper's central result, the zero-temperature damping threshold Eq. (21), appears sound and is confirmed by the numerics. The error in Eq. (19) is a genuine but localized mistake in the spin summation of the damping coefficient; it is fixable and does not invalidate the main claim. The authors should be asked to correct Eq. (19), revise the 'full agreement' statement in Sec. V, and add the missing justifications and numerical-precision notes listed in the minor comments."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, you should know this paper gives a clean necessary condition for zero-temperature Landau damping of the Anderson-Bogolyubov mode in spin- and mass-imbalanced Fermi mixtures, Eq. (21), and it survives scrutiny. The genuinely new piece is the closed condition in terms of the imbalance parameters, plus the observation that for r > 3.01 the damping region unavoidably reaches a continuous quantum critical point. The derivation is careful: they extract the |q_m|/|q| nonlocality from the Gaussian pair fluctuation propagator, and the claim that only σ=σ' terms matter is actually justified—the σ≠σ' denominators are gapped by Δ, and the linear-iq_m piece is odd in cosθ and integrates to zero.\n\nThe main soft spot is Eq. (19). At each root ε_i of E^σ_ε = 0, summing over σ counts both spin branches, but for Δ > 0 only one branch can vanish at a given ε_i (for h−ζμ ≥ 0, only σ=+ contributes). So the literal expression for γ^{-1}_q is overcounted, roughly by a factor of two in the mass-balanced limit. That does not change the existence condition (21), which only needs ε_i ≥ 0, but it does break the paper's advertised 'full quantitative agreement' for the damping rate itself. The numerical comparison in Figs. 4–6 mostly checks the activation threshold, so the error is not glaring, but a referee should ask for the prefactor to be fixed and the comparison redone.\n\nThe numerical test is also a consistency check, not an independent benchmark—same GPF propagator, no free parameters besides the cutoff. That is acceptable for an analytic derivation, but a direct quantitative check against Matera-Wagner (Ref. 60) in the r=1 limit would strengthen the prefactor claim. No code or data is shipped, which makes the prefactor question hard to resolve from the text alone.\n\nWho gets value: cold-atom theorists working on imbalanced superfluids, and people interested in the fate of Goldstone modes near quantum critical points. The paper deserves a serious referee—the central criterion is new, the derivation is mostly sound, and the flaw is localized and fixable. I recommend sending it to review with a request to correct Eq. (19) and clarify the numerical comparison.","headline":"Clean analytic threshold for T=0 Landau damping in mass-imbalanced Fermi superfluids, with a localized prefactor error in the damping-rate formula that does not sink the central result.","tokens_in":14449,"tokens_out":2058,"would_cite":true,"duration_ms":20833,"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":"Fermi-surface mismatch damps the superfluid sound mode even at zero temperature, with a threshold set by Eq. (21).","keywords":["Anderson-Bogolyubov mode","Landau damping","Fermi mixture","mass imbalance","spin imbalance","pair fluctuation propagator","gradient expansion","zero-temperature damping"],"falsifier":"Take a mass-imbalanced parameter set used in the paper (e.g. $r=6.67$, $\\mu=0.1$, $\\Delta$ from the mean-field gap equation) and compute the complex pole of the full pair-fluctuation propagator without dropping any terms, at $h$ just below $\\sqrt{\\mu^2+\\Delta^2}$. If the damping rate is not zero (or not orders of magnitude smaller than just above the threshold), or if any $\\sigma\\neq\\sigma'$ term contributes to the $|q_m|/|q|$ coefficient at the same order, the central claim fails.","tokens_in":13468,"feed_emoji":"⚛️","tokens_out":11369,"duration_ms":110527,"temperature":0.7,"pith_summary":"At low temperature the usual Landau-damping channel for the Anderson-Bogolyubov sound mode—phonon scattering off thermally excited fermionic quasiparticles—is exponentially suppressed because the quasiparticle spectrum is gapped. This paper argues that spin and mass imbalance open a second channel that survives at T=0. It derives the necessary condition, Eq. (21), for a nonzero damping rate, and explains it as the occupation of the upper quasiparticle branch caused by Fermi-surface mismatch. If correct, ultracold two-species Fermi mixtures with sufficient imbalance should show damped phonons at temperatures where thermal damping is frozen out. The analytical threshold agrees with numerically computed complex poles of the pair-fluctuation propagator.","feed_headline":"Mass-mismatched Fermi superfluids damp their sound at T=0","feed_subtitle":"A pair-fluctuation analysis pins the threshold imbalance at which the Goldstone phonon gains a nonzero damping rate.","key_machinery":"The central object is the inverse Gaussian pair-fluctuation propagator $F^{-1}_q$ of the two-component Fermi gas, whose determinant's complex root gives the collective mode frequency and damping rate. The load-bearing step is the low-momentum, low-frequency expansion of its matrix elements, in particular the temporally nonlocal term\n$$\n-\\frac{|q_m|}{|q|}\\int d\\varepsilon\\, c_d D_d(\\varepsilon) $u^{2}$_\\varepsilon $v^{2}$_\\varepsilon \\sum_\\$\\sigma$ \\frac{f'(E^\\sigma_\\varepsilon)}{|a_\\$\\sigma$(\\varepsilon)|} = \\frac{|q_m|}{\\gamma_q},\n$$\nwhich is the same structure that produces Landau damping of spin fluctuations in itinerant magnets. At $T\\to 0$ the Fermi function derivative becomes a delta function, so $\\gamma_q^{-1}$ counts the density of zeros of $E^\\sigma_\\varepsilon$; demanding that those zeros occur at nonnegative kinetic energy $\\varepsilon$ turns the analysis into the threshold inequality (21). The retained $\\sigma=\\sigma'$ contribution is what makes the pole of the Gaussian action complex, and hence what damps the mode.","core_discovery":"The paper's central claim is that the gapless sound mode guaranteed by spontaneous breaking of the continuous U(1) symmetry—the Anderson-Bogolyubov phonon—is Landau-damped at T=0 once the Fermi surfaces of the two species are sufficiently mismatched. Within Gaussian pair-fluctuation theory, the authors expand the inverse propagator at small momentum $q$ and frequency $q_m$ and isolate a nonlocal contribution proportional to $|q_m|/|q|$. At zero temperature the coefficient of this term is nonzero precisely when the quasiparticle dispersion $E^\\sigma_\\varepsilon$ has a zero at a physically allowed kinetic energy $\\varepsilon \\ge 0$. This condition is Eq. (21): for $h-\\zeta\\mu\\ge 0$, damping requires $h-\\zeta\\mu\\ge \\Delta\\sqrt{1-\\zeta^2}$ when $\\mu\\ge \\zeta\\Delta/\\sqrt{1-\\zeta^2}$, and $h\\ge \\sqrt{\\mu^2+\\Delta^2}$ otherwise, where $\\zeta=(r-1)/(r+1)$ with $r=m_-/m_+$ measures mass imbalance and $h$ measures spin (chemical-potential) imbalance. Numerically locating the complex pole of the propagator confirms that damping switches on at the predicted $h$ for mass ratios $r=6.67$ and $r=3.47$, while the balanced case $r=1$ remains undamped.","pith_inferences":["Editorial inference: the derivation is at Gaussian (RPA) level; in low dimensions, where long-range order is reduced to algebraic order by fluctuations, the pole structure and the threshold could be modified, so Eq. (21) should be re-examined beyond Gaussian order there.","Editorial inference: because the threshold comes only from the requirement $\\varepsilon_i\\ge 0$, the same condition likely applies in two dimensions, while the magnitude of the damping rate changes through the density of states $D_d(\\varepsilon)$.","Editorial inference: the same nonlocal mechanism should also damp the gapped amplitude mode; deriving the analogous condition by expanding the off-diagonal matrix element $M_{1,2}$ is a direct extension the paper leaves open.","Editorial inference: a cold-atom experiment with a $^6$Li-$^{40}$K mixture could test the sharp onset by measuring the phonon lifetime as the polarization $h$ is swept through $\\sqrt{\\mu^2+\\Delta^2}$."],"forward_implications":["At T=0 a superfluid with enough Fermi-surface mismatch has a damped Anderson-Bogolyubov phonon; the damping rate is nonzero and increases as $h$ moves further beyond the threshold.","The activation of damping occurs well below the critical field for the superfluid-normal transition, so a large part of the superfluid phase diagram features damped sound modes (for the $r=6.67$ parameters used here, onset near $h\\simeq 1.59$ versus $h_c\\simeq 1.97$).","The zero-temperature damping is insensitive to the order of the superfluid-normal transition: numerical checks with a first-order transition ($r=3.47$) and a continuous one ($r=6.67$) both follow the same threshold condition.","For $r>3.01$, in the continuous phase-transition region, the analytical condition is automatically fulfilled for $h>\\mu\\ge 0$, so Landau damping is unavoidable in the proximity of the quantum critical point.","For equal masses ($r=1$) the threshold $h\\ge \\Delta$ is not reached for the studied parameters, recovering the undamped zero-temperature phonon of the balanced case."],"supporting_citations":[{"why":"Supplies the low-momentum, low-frequency gradient expansion of the pair-fluctuation propagator that the paper extends with temporally nonlocal terms.","marker":"[50]"},{"why":"Companion gradient-expansion analysis of the same propagator; its expansion without the $|q_m|/|q|$ term is the baseline that must be supplemented.","marker":"[51]"},{"why":"Gives the mass-balanced limit ($r=1$) in which the damping threshold is $h\\ge \\Delta$; the paper generalizes this to arbitrary mass ratio.","marker":"[60]"},{"why":"Supplies the mean-field phase diagram and the statement that for $r>3.01$ the condition is fulfilled near the quantum critical point, making damping unavoidable there.","marker":"[22]"},{"why":"Provides the Gaussian pair-fluctuation formalism and the mass ratio $r=6.67$ of a $^6$Li-$^{40}$K mixture used in the numerical checks.","marker":"[48]"},{"why":"Supplies the analytic-continuation procedure for locating the complex pole of the pair-fluctuation propagator.","marker":"[54]"},{"why":"Introduces the $|q_m|/|q|$ nonlocal term and its interpretation as Landau damping, the structural template for the paper's damping analysis.","marker":"[65]"},{"why":"Shows that in the gapped case Landau damping is exponentially suppressed at low T in the balanced case, framing the imbalance-activated damping as new.","marker":"[58]"}],"fun_headline_variants":["Mismatched Fermi seas damp superfluid phonons at T=0","Zero-temp Landau damping from Fermi mismatch in superfluids","Mass and spin imbalance kill superfluid sound at T=0","Fermi mismatch triggers T=0 damping of Goldstone phonon"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that among all the scattering processes in the pair-fluctuation bubble, only the one with equal quasiparticle branch indices ($\\sigma=\\sigma'$) contributes the nonlocal term that causes damping at leading order; if any other process contributed at the same order, the damping rate and the threshold condition could be different.","fun_headline_variants_meta":{"raw":{"variants":["Mismatched Fermi seas damp superfluid phonons at T=0","Zero-temp Landau damping from Fermi mismatch in superfluids","Mass and spin imbalance kill superfluid sound at T=0","Fermi mismatch triggers T=0 damping of Goldstone phonon"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000539,"raw_usage":{"total_tokens":2578,"prompt_tokens":930,"completion_tokens":1648,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":546,"completion_tokens_details":{"reasoning_tokens":1573}},"tokens_in":546,"tokens_out":1648,"duration_ms":12587,"temperature":1.0,"reasoning_tokens":1573,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:37:05.624794+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a mass-imbalanced parameter set used in the paper (e.g. $r=6.67$, $\\mu=0.1$, $\\Delta$ from the mean-field gap equation) and compute the complex pole of the full pair-fluctuation propagator without dropping any terms, at $h$ just below $\\sqrt{\\mu^2+\\Delta^2}$. If the damping rate is not zero (or not orders of magnitude smaller than just above the threshold), or if any $\\sigma\\neq\\sigma'$ term contributes to the $|q_m|/|q|$ coefficient at the same order, the central claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Companion gradient-expansion analysis of the same propagator; its expansion without the $|q_m|/|q|$ term is the baseline that must be supplemented."},{"cited_title":"\\ Matera \\ and\\ author M","cited_arxiv_id":null,"evidence_quote":"Gives the mass-balanced limit ($r=1$) in which the damping threshold is $h\\ge \\Delta$; the paper generalizes this to arbitrary mass ratio."},{"cited_title":"\\ Zdybel \\ and\\ author P","cited_arxiv_id":null,"evidence_quote":"Supplies the mean-field phase diagram and the statement that for $r>3.01$ the condition is fulfilled near the quantum critical point, making damping unavoidable there."},{"cited_title":"\\ Iskin \\ and\\ author C","cited_arxiv_id":null,"evidence_quote":"Provides the Gaussian pair-fluctuation formalism and the mass ratio $r=6.67$ of a $^6$Li-$^{40}$K mixture used in the numerical checks."},{"cited_title":"\\ Kurkjian \\ and\\ author J","cited_arxiv_id":null,"evidence_quote":"Shows that in the gapped case Landau damping is exponentially suppressed at low T in the balanced case, framing the imbalance-activated damping as new."}],"review_version":1}