{"id":"d11eb8bc-4abf-45a6-b505-b86602084c8d","arxiv_id":"2504.14591","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For a massive Maxwell-Boltzmann gas in the relaxation-time approximation, finite mass removes the propagating sound mode, leaving purely imaginary modes, while shear and diffusion modes stay close to the massless results.","lead":"This paper computes how a gas of massive relativistic particles responds to small disturbances, using a simplified kinetic theory with a single collision time. It finds that a finite particle mass can stop sound-like collective modes from propagating, and it derives viscosity coefficients for both light and heavy quasiparticles.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The sound-channel conclusion in Eq. (67) depends on taking x→0 before k→0; the singular roots (ω1 ∼ k^2/x^6, ω3 ∼ x^6) do not rule out propagating sound for finite x at larger k.","rationale":"The reader's weakest assumption identifies the same load-bearing issue: the small-mass expansion and the hydrodynamic limit are taken in a fixed order, and Eq. (67) shows singular x-dependence that makes the two limits non-commuting. My stress-test confirms and sharpens this: the massless limit at fixed kτ gives propagating sound, while the kτ→0 limit at fixed x gives purely imaginary modes, so the headline statement is order-dependent. The paper has genuine supporting evidence — the massless limit reproduces the known results of Ref. [27] in the shear and charge sectors, the Ward identities are checked, and the analytic expressions are explicit — so this is not a rejection-level flaw. It is, however, exactly the condition that makes the central claim unverified: the paper does not solve the unexpanded pole equation, and its own final paragraph suggests that a numerical treatment is the natural next step. If the proposed numerical test shows Reω = 0 at all hydrodynamic k for finite x, the claim is confirmed; if it shows propagating sound above an x-dependent threshold, the abstract and conclusions must be qualified. The reader's CONDITIONAL verdict is therefore appropriate, and no verdict change is needed.","tokens_in":26138,"tokens_out":14361,"duration_ms":132327,"concrete_test":"Solve the full sound-channel pole equation D_sound(ω,k,x)=0 for G00,00, using the unexpanded denominator obtained from the Appendix B integrals, for x = 0.05, 0.1, 0.2 and kτeq in [1e-4, 1e-1]. Track the two smallest nonzero poles: if their real parts are zero below an x-dependent crossover scale but become nonzero above it, the purely-imaginary claim is an order-of-limits artifact; if Reω = 0 at all hydrodynamic k for every finite x, the claim survives. Also check the trend as kτeq is scaled like x^2 and x^3 to identify the crossover.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, that finite mass converts the propagating sound mode into a purely imaginary mode, is extracted from the small-x expansion of the sound denominator made before the small-k hydrodynamic limit. Eq. (67) gives ω1 ≈ 8ik^2τeq/[x^6(γE−ln2)] and ω3 ≈ ix^6(γE−ln2)/(24τeq) + 8ik^2τeq/[x^6(ln2−γE)]. Both roots diverge as x→0 at fixed k, so the small-x expansion is non-uniform and these roots cannot by themselves establish the behavior of the full theory. The denominator near ω=0 contains x-dependent terms that shift the zero which in the massless case produces Reω = c_s k; the real part vanishes only when kτ is below an x-dependent crossover scale, while for larger kτ the same roots can acquire real parts and become propagating. The abstract states the claim without this restriction, so as written the headline result is an artifact of a particular order of limits rather than a demonstrated property of massive RTA kinetic theory. The paper's own closing suggestion, to perform the integrals numerically at fixed x, is precisely the check that would resolve this.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the linearized relaxation-time-approximation (RTA) Boltzmann equation for a massive Maxwell-Boltzmann gas with constant mass m. Using a variational approach, it derives retarded charge and energy-momentum correlators at arbitrary (ω,k), then expands in small and large x = m/T. It reports diffusion constants, shear and bulk viscosities, third-order transport coefficients, and hydrodynamic mode spectra. The main novel claim is that in the sound channel a finite mass removes the propagating sound pole and replaces it with purely imaginary modes.","tokens_in":26354,"tokens_out":4342,"duration_ms":39721,"significance":"If the central claim survives scrutiny, the paper is a useful first analytic treatment of massive RTA correlators: it provides explicit expressions, satisfies the relevant Ward identities, and connects to known massless results in limiting cases. The transport-coefficient extraction against the third-order hydrodynamic ansatz of [56] is a strong external check. However, the headline sound-channel conclusion depends critically on an assumption about commuting the x→0 and k→0 limits, and that assumption is not tested; the paper's own closing sentence identifies the numerical evaluation that would settle the issue.","major_comments":[{"comment":"The conclusion that finite mass converts propagating sound into purely imaginary modes is extracted from the small-x expansion of the sound denominator performed before the small-k hydrodynamic limit. The roots in Eq. (67), ω1 ~ 8ik²τ/(x⁶(γ_E−ln2)) and ω3 ~ ix⁶(γ_E−ln2)/(24τ) + 8ik²τ/(x⁶(ln2−γ_E)), are non-uniform in x: they diverge as x→0 at fixed k, so the two limits do not commute. The paper does not demonstrate that these roots describe poles of the full correlator for small finite x, and the numerical confirmation in Fig. 5 uses the same x-expanded expression, so it is not an independent check. A fixed-x numerical evaluation of the original momentum integrals, as suggested in the conclusion, or a controlled double-scaling analysis is required before the abstract's claim can be supported.","section":"Section IV, Eq. (67)"},{"comment":"The bulk viscosity formula in Eq. (63) and the large-x third-order transport coefficients in Eq. (57) are presented after 'some computations' without derivation. These are load-bearing quantitative results of the transport section; without intermediate steps or an appendix reproducing the calculation, the results are not verifiable. The authors should provide the derivation or a reproducible algebraic outline.","section":"Section IV, Eqs. (63) and (57)"},{"comment":"The expansion of (1/v) ln[(ω−kv+i/τ)/(ω+kv+i/τ)] is an expansion in x/|p|, and its correction terms contain denominators such as ((1−iτeqω)²+k²τeq²) that vanish in the hydrodynamic limit. Using this expansion before taking k→0 can generate spurious poles of order 1/x⁶, which is the same non-commutativity issue identified in the first comment. The paper should state the intended ordering of limits and justify it explicitly.","section":"Section III, Eq. (23)"}],"minor_comments":[{"comment":"The equilibrium distribution in Eq. (2), feq = exp[(gαβ p^α u^β + μ)/T], appears to have the opposite sign from Eq. (9), feq = exp[−(p0−μ0)/T]. For the chosen metric signature and u=(1,0,0,0), these expressions are inconsistent; please clarify the intended sign convention.","section":"Section II, Eqs. (2) and (9)"},{"comment":"The Introduction states 'in the weak coupling regime (τeqT→∞)' and 'in the strong coupling regime (τeqT→∞)' with the same arrow; the strong-coupling limit should presumably be τeqT→0.","section":"Introduction"},{"comment":"The sentence 'the first term on the right-hand side does not contribute' refers to the electromagnetic force, but the mass-gradient term M∂^α M also vanishes for constant mass; this should be stated explicitly.","section":"Section II, after Eq. (4)"},{"comment":"The notation τ and τeq is used interchangeably in these equations; please standardize the relaxation-time symbol throughout Section IV.","section":"Section IV, Eqs. (52) and (55)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid technical contribution to kinetic-theory correlators, but the main physical claim about sound-channel modes requires the fixed-x numerical check that the authors themselves suggest. The derivation is self-contained and the transport matching to the hydrodynamic ansatz is a good feature. I see no citation or novelty concerns, only the need to close the order-of-limits gap before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a real calculation: they extend Romatschke's massless RTA retarded-correlator program to fixed-mass Maxwell-Boltzmann particles, giving explicit correlators, mode spectra, and transport coefficients through third order in both small and large m/T. Second, the advertised result—that mass kills the propagating sound mode—is probably an artifact of expanding in x=m/T before taking the hydrodynamic limit, and as written the paper does not demonstrate it.\n\nWhat is genuinely useful: the machinery is self-contained and the massless limit checks out in the shear and charge sectors, which gives confidence in the formalism. The bulk-viscosity formula, the small- and large-x diffusion constants, and the third-order coefficient combinations are new explicit outputs. The paper is also honest about its own limitations; the closing suggestion to do the integrals numerically at fixed x is exactly the right check.\n\nThe soft spot is load-bearing. Eq. (67) gives sound-channel poles with terms like ω1 ~ k^2/(x^6) and ω3 ~ x^6; both blow up as x→0 at fixed k. That non-uniformity means the small-x expansion of the denominator cannot by itself tell you what happens to the true sound poles at small finite x. The stress-test note is right: for any fixed x, there should be an x-dependent crossover scale in kτ below which the real part vanishes and above which the massive roots can acquire real parts. The abstract and conclusion state the stronger claim unconditionally. Until the authors evaluate the correlators at fixed x (numerically is fine) and show the poles have no real part, the central result is conditional.\n\nMinor soft spots: the bulk viscosity expression (Eq. 63) and the large-x third-order coefficients (Eq. 57) are quoted without derivation or a reproducibility aid; that matters because some of those results—e.g., vanishing of one third-order combination at large x—are surprising and need checking. The paper says it follows the strategy of x expansion, and for a calculation this long, shipping the algebra or code would remove a lot of doubt. The matching to the hydrodynamic ansatz from [56] is appropriate, but the uniqueness of that matching is not discussed.\n\nBottom line: this deserves a serious referee, not a desk reject, but the referee should send it back for a fixed-mass numerical check of the sound poles. If that check supports the purely-imaginary claim, this becomes a useful paper. If not, the substantial parts—the correlators and transport coefficients—may still stand as a reference for massive RTA. I would not cite it until the sound-mode question is settled.","headline":"A serious analytic extension of RTA kinetic theory to massive particles, but the headline sound-mode claim likely depends on an untested order of limits.","tokens_in":26894,"tokens_out":3360,"would_cite":false,"duration_ms":29615,"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":"Finite particle mass turns the propagating sound mode of a massive relativistic gas into a purely imaginary, dissipative mode.","keywords":["relativistic kinetic theory","relaxation time approximation","retarded Green's functions","quasiparticles","hydrodynamic modes","sound channel","branch cut","bulk viscosity"],"falsifier":"Evaluate the sound-channel pole condition $D_{\\rm sound}(\\omega,k,x)=0$ at fixed finite $x$ (say $x=0.1$) directly in small $k$ without first expanding in $x$; if the lowest pole has a real part proportional to $k$, the claimed conversion of propagating sound into a purely imaginary mode fails. The same test can be done by numerically extracting the spectral function peak of $G_{00,00}$ at finite $x$.","tokens_in":25900,"feed_emoji":"🔊","tokens_out":5332,"duration_ms":48520,"temperature":0.7,"pith_summary":"This paper claims that giving the particles in a relaxation-time-approximation Boltzmann gas a finite rest mass qualitatively changes its collective dynamics. Working with a constant mass profile and a Maxwell-Boltzmann equilibrium distribution, the authors derive the retarded two-point correlation functions for charge and energy-momentum currents and study their poles and branch cuts at small and large mass-to-temperature ratio $x=m/T$. The central result is that in the sound channel the usual propagating sound mode disappears: finite mass converts it into a purely imaginary, dissipative mode, while the shear channel keeps modes matching the massless results asymptotically. The paper also extracts transport coefficients, including bulk and shear viscosity and third-order corrections, as expansions in $x$. If correct, this gives a concrete microscopic model in which the presence of a mass scale alone reshapes the hydrodynamic spectrum.","feed_headline":"Mass turns sound into a purely imaginary mode","feed_subtitle":"A massive Boltzmann gas loses its propagating sound mode, while shear modes match massless results.","key_machinery":"The load-bearing object is the set of retarded two-point functions obtained from the linearized Boltzmann equation by the variational method: the response of the induced current and stress tensor to weak metric and gauge perturbations. Poles of these correlators define the collective modes, and logarithmic terms of the form $\\ln\\!\\left(\\frac{\\omega-kv+i/\\tau_{\\rm eq}}{\\omega+kv+i/\\tau_{\\rm eq}}\\right)$ integrated over the massive Maxwell-Boltzmann distribution generate both the branch cuts and the hydrodynamic denominators. The analysis proceeds by expanding those integrals in $x=m/T$ before taking the small-$(\\omega,k)$ limit, which produces explicit dispersion relations and the transport coefficients.","core_discovery":"On the paper's own terms, the discovery is that for a massive relativistic gas in the relaxation-time approximation the retarded energy-momentum correlator in the sound channel has only purely imaginary poles at small momentum. Expanding around $x=m/T$, the three lowest poles are $\\omega_1 \\simeq \\frac{8 i k^2 \\tau_{\\rm eq}}{x^6}(\\gamma_E-\\ln 2)$, $\\omega_2 \\simeq -\\frac{i}{\\tau_{\\rm eq}} - i k^2 \\tau_{\\rm eq}\\left(\\frac{4}{15}-\\frac{x^2}{45}\\right)$, and $\\omega_3 \\simeq \\frac{i x^6(\\gamma_E-\\ln 2)}{24\\tau_{\\rm eq}}$ plus higher-order terms; none has a real part linear in $k$, so the mode that would be propagating sound in the massless limit becomes dissipative. The same expansion yields shear-channel modes that reduce to the massless spectrum, a logarithmic branch cut between $\\omega=k$ and $\\omega=-k$ in the weak-coupling limit, and hydrodynamic poles above this cut in the strong-coupling regime below a critical $k\\tau_{\\rm eq}$. Transport coefficients are then computed by matching the low-frequency correlators to the hydrodynamic Kubo relations.","pith_inferences":["The paper expands in $x$ before taking the hydrodynamic limit; if the two limits do not commute, the true poles at finite small $x$ could retain a real part proportional to $k$, restoring propagating sound. The singular coefficients $1/x^6$ and $x^6$ in Eq. (67) make this a concrete alternative to check.","One could settle the question numerically by solving the full pole condition for $x$ around $0.1$ without the $x$-expansion; a nonzero $\\operatorname{Re}\\omega \\propto k$ at fixed small $x$ would contradict the claim.","If the claim holds, quasiparticle models with a running mass, such as those used near a chiral critical point, would exhibit a qualitative change in hydrodynamic mode structure at the mass scale, with consequences for how long hydrodynamics remains valid."],"forward_implications":["In a massive Maxwell-Boltzmann gas in the relaxation-time approximation, sound-channel hydrodynamic modes do not propagate: all three lowest poles are purely imaginary, so a density perturbation decays without oscillating.","Shear-channel modes asymptotically approach their massless counterparts, so the mass changes only coefficients, not the qualitative structure, in that channel.","The logarithmic branch cut between $\\omega=k$ and $\\omega=-k$ persists for massive particles, and hydrodynamic poles appear above the cut in the strong-coupling regime only below a critical $k\\tau_{\\rm eq}$.","The computed transport coefficients, including a bulk viscosity that starts at $O(x^4)$ for small $x$, provide concrete massive corrections to standard massless relaxation-time-approximation results."],"supporting_citations":[{"why":"supplies the massless relaxation-time-approximation correlation functions and branch-cut/pole analysis that this paper extends to massive particles.","marker":"[27]"},{"why":"introduces the Bhatnagar-Gross-Krook relaxation-time collision term used as the kinetic model.","marker":"[50]"},{"why":"provides the relativistic Anderson-Witting form of the Boltzmann equation adopted here.","marker":"[51]"},{"why":"supplies the massive Maxwell-Boltzmann thermodynamics and Bessel-function integrals used in the equation of state.","marker":"[54]"},{"why":"gives the third-order hydrodynamic expansion of the shear correlator used to read off higher-order transport coefficients.","marker":"[56]"},{"why":"provides the Kubo-type formula relating bulk viscosity to the trace-of-stress correlator.","marker":"[57]"}],"fun_headline_variants":["Mass flips sound to a purely imaginary pole","Mass damps sound to a purely imaginary mode","Mass turns sound non-propagating in a hot gas","Mass turns sound into a damped imaginary pole"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the expansion in $x=m/T$ can be performed before the small-momentum hydrodynamic limit and that these two limits commute.","fun_headline_variants_meta":{"raw":{"variants":["Mass flips sound to a purely imaginary pole","Mass damps sound to a purely imaginary mode","Mass turns sound non-propagating in a hot gas","Mass turns sound into a damped imaginary pole"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000936,"raw_usage":{"total_tokens":4046,"prompt_tokens":1033,"completion_tokens":3013,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":649,"completion_tokens_details":{"reasoning_tokens":2952}},"tokens_in":649,"tokens_out":3013,"duration_ms":18457,"temperature":1.0,"reasoning_tokens":2952,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:45:42.149788+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the sound-channel pole condition $D_{\\rm sound}(\\omega,k,x)=0$ at fixed finite $x$ (say $x=0.1$) directly in small $k$ without first expanding in $x$; if the lowest pole has a real part proportional to $k$, the claimed conversion of propagating sound into a purely imaginary mode fails. The same test can be done by numerically extracting the spectral function peak of $G_{00,00}$ at finite $x$.","supporting_citations":[{"cited_title":"A Model for Collision Processes in Gases. 1. Small Amplitude Processes in Charged and Neutral One-Component Systems,","cited_arxiv_id":null,"evidence_quote":"supplies the massive Maxwell-Boltzmann thermodynamics and Bessel-function integrals used in the equation of state."},{"cited_title":"CHIRALLY INVARIANT TRANSPORT EQUATIONS FOR QUARK MATTER","cited_arxiv_id":"hep-ph/9505407","evidence_quote":"gives the third-order hydrodynamic expansion of the shear correlator used to read off higher-order transport coefficients."}],"review_version":1}