{"id":"605b60fa-1ea1-422a-9a96-37f68c92c346","arxiv_id":"2506.15531","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper derives closed-form mass-dependent transport coefficients for massive RTA gases and proposes that the discontinuity across the correlator cut defines an effective lightcone velocity.","lead":"Using the relaxation time approximation of kinetic theory, this paper computes exact transport coefficients for a massive relativistic gas, including thermoelectric response. It also shows how the analytic structure of retarded correlators can be read as a causal lightcone, with a cut encoding the maximum and typical particle speeds.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (59), the closed-form shear viscosity, is dimensionally inconsistent as printed: it is missing a factor 1/T_0, so at m=0 it gives η = 4τ_R T_0^5/(5π^2), not the paper's own cross-check value η = D_sh(ε_0+P_0) = 4τ_R T_0^4/(5π^2).","rationale":"The reader's verdict identifies the cut-indeformability argument (Section IV B) as the weakest assumption; I agree that this is a genuine gap, and I would also flag the threshold-dependence of v_cut in Eq. (90): the toy model in Eq. (B16) makes v_cut(ε) explicit, and no physical principle fixes ε. However, on reading the transport sections closely I found a more concrete, checkable problem that the reader's weakest_assumption did not catch: Eq. (59), a headline closed-form formula, is dimensionally wrong as printed. The integral ∫_m^∞ dx x^{-1}(x^2−m^2)^{5/2} e^{−x/T_0} scales as T_0^5, so with the prefactor τ_R/(30π^2) the right-hand side has the dimension of pressure, not of viscosity; at m=0 it gives (4/5)τ_R T_0^5/π^2, whereas the paper's own identifications η = D_sh(ε_0+P_0), with D_sh = τ_R/5 (40) and ε_0+P_0 = 4T_0^4/π^2 (19)–(20), give (4/5)τ_R T_0^4/π^2. The factor of T_0 is present in the analogous Eq. (71) (χD = τ_R/(6π^2 T_0)∫...), so the omission in (59) is almost certainly a transcription error rather than a conceptual failure; the cross-checks cited (Kubo (58), D_sh identification, [46]) are consistent with the corrected normalization. I therefore do not recommend a harsher verdict than the reader's: CONDITIONAL is right, since the paper should (i) correct Eq. (59), and (ii) either prove the indeformability (e.g., by a monodromy analysis around ω = k − i/τ_R) or soften the causal interpretation. I also credit the solid parts: the transport coefficients D and D_sh follow from ratios of manifestly positive integrals (A10), (A21), so positivity for all m and the correction of the [31] instability are secure; the thermoelectric relations (71)–(79) pass dimensional and massless checks; and the numerical verifications of the discontinuity profile (Fig. 5) and truncated-cut Fourier transforms (Fig. 8) are genuine supporting evidence. The decisive next step is the trivial dimensional/cross-check test above, which any referee can run in minutes.","tokens_in":22329,"tokens_out":60144,"duration_ms":559989,"concrete_test":"Set T_0 as a symbolic parameter and evaluate the paper's cross-check at m=0. From (59) compute η = 4τ_R T_0^5/(5π^2); from (39)–(40) together with (19)–(20) compute η = D_sh(ε_0+P_0) = 4τ_R T_0^4/(5π^2). If replacing the prefactor in (59) by τ_R/(30π^2 T_0) makes the two agree for all T_0, the printed formula is confirmed misprinted and the corrected one is the shear viscosity. A purely numerical variant: evaluate η from (59) and from D_sh(ε_0+P_0) at T_0 = 2 with m/T_0 fixed (say 1); the printed version fails by a factor of 2, the corrected version passes. The same check settles the claimed agreement with [46].","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III B 3, Eq. (59), is printed as η = τ_R/(30π^2) ∫_m^∞ dx x^{-1}(x^2−m^2)^{5/2} e^{−x/T_0}. This is dimensionally inconsistent: the integral scales as T_0^5 and τ_R as T_0^{-1}, so the right-hand side has dimension [E]^4 (pressure), while viscosity scales as [E]^3. At m=0 the printed formula evaluates to η = 4τ_R T_0^5/(5π^2). The paper itself requires η = D_sh(ε_0+P_0) (stated in Section III B 1 and said to be independently confirmed in Section III B 3) and agreement with the Chapman–Enskog results of [46]. Using the paper's Eqs. (19)–(20) and (40), at m=0 one has ε_0+P_0 = 4T_0^4/π^2 and D_sh = τ_R/5, giving η = D_sh(ε_0+P_0) = 4τ_R T_0^4/(5π^2): one power of T_0 smaller than Eq. (59) as printed. The same integral appears with the 1/T_0 factor present in Eq. (71) (τ_R/(6π^2 T_0)∫... = χD, dimensionally consistent), indicating Eq. (59) lost a 1/T_0 in transcription. Because Eq. (59) is one of the five closed-form transport coefficients named in the central claim, the claim 'complete analytic expressions for all first-order transport coefficients' is not true as printed for η at T_0 ≠ 1; any reader implementing (59) obtains viscosities off by a factor T_0. The positivity/stability physics is unaffected: D and D_sh are ratios of positive integrals (A10), (A21), and the correction of [31] stands. Secondary but real: the reader's flagged gap in Section IV B remains — indeformability is argued only via the 2D-hole property of deformed-contour representations, and the paper itself concedes 'This may be inconsistent with the general properties of the causal correlators' — and v_cut in Eq. (90) depends on an arbitrary threshold ε (explicitly so in Eq. (B16)). These affect the interpretive causal claims; Eq. (59) is a stated result that is simply wrong as printed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies massive relativistic particles in the relaxation time approximation of kinetic theory at finite temperature and density. It derives closed-form expressions for the first-order transport coefficients—charge and shear diffusion constants D and D_sh, bulk viscosity ζ, shear viscosity η, and DC thermoelectric conductivity σ_Q—as functions of m/T0, using the variational/Kubo approach to retarded correlators. It also analyzes the analytic structure of the correlators, identifies a cut at Im ω = -1/τ_R between ω = ±k - i/τ_R, studies its discontinuity profile as a function of mass, and proposes a criterion (disc G = ε) that extracts an 'effective lightcone velocity' v_cut. The paper further claims to correct an instability reported in Ref. [31], showing positive diffusion coefficients for all masses in d>1, and discusses the indeformability of the cut at finite mass.","tokens_in":22865,"tokens_out":6800,"duration_ms":70682,"significance":"If the results hold, the paper provides a useful exact reference: all first-order transport coefficients of the massive RTA gas in closed form, with cross-checks against Kubo formulas, Ward identities, Chapman–Enskog results, and positivity/stability in arbitrary dimensions. The paper contains no fitted parameters in the transport coefficients and gives explicit analytic formulas (e.g., (32), (39), (50), (77)) that can be implemented directly. These strengths are offset by two caveats: the printed shear viscosity (59) is dimensionally inconsistent, and the claimed indeformability/uniqueness of the cut in Section IV B is argued heuristically rather than proven. The v_cut criterion also depends on an arbitrary threshold ε. The central positivity result is robust, but the 'complete analytic expressions' claim is not true as printed until Eq. (59) is corrected.","major_comments":[{"comment":"Equation (59) as printed is dimensionally inconsistent: the right-hand side has energy dimension four (τ_R times an integral scaling as T0^5), whereas the shear viscosity η has energy dimension three. The missing factor is 1/T0 in front of the integral, as is visible by comparing with Eq. (71) and with the required cross-check η = D_sh(ε0 + P0). At m=0, Eq. (59) gives η = 4τ_R T0^5/(5π^2) instead of 4τ_R T0^4/(5π^2). Since Eq. (59) is one of the five closed-form transport coefficients advertised in the abstract, the central claim is not true as printed; any reader implementing (59) obtains viscosities off by a factor of T0. The positivity/stability conclusion is unaffected because D and D_sh are ratios of positive integrals, but the formula must be corrected before the manuscript can be accepted.","section":"III B 3, Eq. (59)"},{"comment":"The claim that the cut is undeformable and unique is not established. The argument that deforming the angular contour removes a two-dimensional region and 'may be inconsistent' with causal correlators does not prove that the real-u, real-z contour is forced by causality; the sentence 'this implies a unique choice of the cut' overstates the result. As the authors themselves note, other analytic continuations might be physically allowed, and the discussion of the pole merging with the cut at k_* depends on this choice. The paper should either provide a proof from analyticity/causality of the retarded correlator or explicitly present the straight cut as a preferred convention rather than as a uniqueness statement.","section":"IV B, Eq. (94)"},{"comment":"The definition of v_cut via disc G = ε depends on an arbitrarily chosen threshold ε, and no physical principle selects a value of ε. The numerical demonstrations use ε = 10^-3 without showing that v_cut is stable or meaningful under changes of ε, and the same ε-dependence enters the toy model expression (B16). As it stands, the 'extraction of physical information from the discontinuity profile' is convention-dependent, which weakens the claim that v_cut encodes an effective lightcone. The authors should quantify the ε-dependence and either justify a canonical choice or soften the interpretation.","section":"IV A, Eq. (90)"}],"minor_comments":[{"comment":"The symbol m is used both for the dimensionful mass (in the integration limits and in quantities such as ε0 and P0) and for the dimensionless ratio m/T0 (in the Bessel functions and the Meijer-G expressions). This is confusing; a separate notation such as m̂ = m/T0 would improve readability.","section":"II A, Eqs. (17)–(20)"},{"comment":"The abstract advertises 'finite density', but the thermodynamic expressions in Section II A are derived 'in the absence of chemical potential'. The authors should clarify whether finite density means nonzero n0 at zero chemical potential or whether µ0 ≠ 0 results are included, since the thermoelectric section later uses n0 and s0.","section":"Abstract and Section II A"},{"comment":"The caption writes the cut endpoints as ω = ±ck − i/τ_R, while the text and Eq. (95) use ±k − i/τ_R with c implicitly set to unity. The speed of light c should be defined or omitted to avoid confusion.","section":"Figure 1 caption"},{"comment":"The statement that the cutoff wavevector k_* (where the diffusive pole reaches the cut) increases with mass and spatial dimension is made without presenting the corresponding formula or numerical evidence, making it difficult to verify.","section":"IV B, k_* discussion"},{"comment":"In the toy-model definition of v_cut, the result depends on ε through both log ε and ε itself; using a different normalization of the discontinuity profile would change v_cut. It would be helpful to state explicitly that v_cut is a threshold-dependent quantity and to give its behavior for representative ε values.","section":"Appendix B, Eq. (B16)"}],"recommendation":"major_revision","confidential_remarks":"The dimensional typo in Eq. (59) is likely fixable, but it affects one of the five advertised closed-form transport coefficients and must be corrected. The indeformability and v_cut claims should be reframed as a preferred-convention analysis or supported by a proof. I believe the paper is publishable after major revision; the positivity result and the analytic transport formulas are valuable contributions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the complete set of mass-dependent first-order RTA transport coefficients in d dimensions, including the thermoelectric coefficients, and the explicit demonstration that D and D_sh stay positive, which corrects the instability claim in [31]. The derivations are cross-checked against Kubo formulas, Ward identities, and Chapman-Enskog, so the underlying math seems right.\n\nBut there is a concrete error in the printed Eq. (59): the shear viscosity formula is dimensionally inconsistent. The integral scales as T0^5, τ_R as 1/T0, so the right-hand side has dimension energy^4, not energy^3. At m=0 it evaluates to η = 4τ_R T0^5/(5π^2), while the paper's own relation η = D_sh(ε0+P0) gives 4τ_R T0^4/(5π^2). The analogous integral in Eq. (71) has the missing 1/T0 factor and is dimensionally consistent. This looks like a transcription error, but as printed the central claim \"complete analytic expressions for all first-order transport coefficients\" is false for η, and anyone implementing (59) gets viscosities off by a factor T0.\n\nThe other soft spots are interpretive rather than mathematical. The \"indeformable cut\" argument in Sec. IV B relies on keeping u in [0,1] and z in [-1,1] real; the paper does not prove causality forces that contour, and it concedes other choices may be inconsistent with causal correlators. The effective lightcone velocity vcut depends on an arbitrary threshold ε, as shown in Eq. (B16). These weaken the causal claims but do not undermine the transport coefficients.\n\nBottom line: the transport part is a solid and useful contribution for people working on RTA kinetic theory and heavy-ion phenomenology. The cut interpretation is suggestive and nicely illustrated but not yet rigorous. The paper deserves a serious referee; the main fixes are a corrected Eq. (59) and a more careful framing of the indeformability and vcut claims.","headline":"Solid analytic transport results for massive RTA, but Eq. (59) is dimensionally off by a factor T0 and the causal-cut interpretation is heuristic.","tokens_in":23417,"tokens_out":3583,"would_cite":true,"duration_ms":35193,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.60.-k","05.20.Dd"],"model":"deepseek-v4-flash","headline":"All first-order transport coefficients of a massive RTA gas are now known in closed form and remain positive for every mass.","keywords":["relaxation time approximation","massive kinetic theory","transport coefficients","thermoelectric coefficients","retarded correlators","branch cuts","effective lightcone velocity","causality"],"falsifier":"Evaluate the exact diffusion constants, for instance Eq. (39) in $d=3$ or Eq. (A22) in $d=2$, at large $m/T_0$: if either $D$ or $D_{\\rm sh}$ turns negative at any mass, the positivity claim fails. Separately, deform the angular integration contour $\\xi_c$ in Eq. (94) and test whether the resulting correlator satisfies the usual causal analyticity conditions; a causal deformed-contour correlator with a different cut would refute the claimed indeformability.","tokens_in":22124,"feed_emoji":"","tokens_out":6412,"duration_ms":62395,"temperature":0.7,"pith_summary":"The paper claims that the first-order transport coefficients of a massive gas in the relaxation time approximation of kinetic theory—charge and shear diffusion, shear and bulk viscosity, and the thermoelectric DC conductivity—can be written in closed analytic form as functions of $m/T_0$ in any spatial dimension $d>1$. If the claim is right, first-order hydrodynamics of this model is known exactly, with no free parameters beyond mass, temperature, chemical potential, and relaxation time. The results also correct an earlier report of a diffusive instability: the coefficients are positive for every mass, with controlled light- and heavy-mass asymptotics. The same calculation maps the non-hydrodynamic structure of retarded correlators: finite mass replaces the freely deformable branch cut of the massless case with a unique, undeformable cut, and the discontinuity profile along that cut yields an effective lightcone velocity.","feed_headline":"Massive gas transport coefficients computed exactly; none turn negative","feed_subtitle":"The paper gives closed-form transport coefficients for the massive RTA gas and shows all stay positive, correcting a claimed instability.","key_machinery":"The load-bearing object is the retarded two-point correlator written as a double integral over energy, $u\\in[0,1]$, and angle, $z\\in[-1,1]$, as in Eq. (94), carrying the Maxwell-J\\\"uttner weight $\\exp(-m/(T_0\\sqrt{1-u^2}))$. The mass enters through the particle velocity $\\xi=\\sqrt{1-m^2/x^2}$, coupling the momentum integral to the angular integral so that the denominator $1-i\\omega\\tau_R+ik\\tau_R z u$ has a continuum of zeros, generating a cut instead of a movable logarithmic branch point. The argument then extracts transport coefficients from the diffusive pole of this integral in the hydrodynamic limit, computes the discontinuity profile by the residue-style procedure of Appendix B, and defines $v_{\\rm cut}$ as the threshold $\\epsilon$ below which the discontinuity is exponentially suppressed.","core_discovery":"On the paper's own terms, the central discovery is that mass does not destabilize RTA hydrodynamics but reshapes its analytic structure. The diffusive-pole integral equations are solved in closed form: $D/\\tau_R$ is a Meijer $G$-function combination in Eq. (32), $D_{\\rm sh}/\\tau_R$ in Eq. (39), the bulk viscosity in Eq. (50), the shear viscosity in Eq. (59), and the thermoelectric DC conductivity $\\sigma_Q$ in Eq. (77); all remain positive for all masses and reduce correctly in the $m\\to 0$ and $m\\to\\infty$ limits. In the complex frequency plane, finite mass turns the logarithmic branch cut of the massless theory into a straight cut between $\\omega=\\pm k-i/\\tau_R$ that cannot be deformed, because the energy and angular integrations are coupled: deforming either contour removes a two-dimensional region from the domain of the correlator. The endpoints of this cut mark the maximal propagation speed, while the shape of the discontinuity along it defines a threshold criterion $v_{\\rm cut}$ (Eq. (90)) for the effective ballistic lightcone, which is distinct from the mean or rms particle speed and from the speed of sound.","pith_inferences":["If the indeformability argument carries over to RTA with energy-dependent relaxation time or to weakly coupled massive field theories, one expects the same two-dimensional non-analytic region rather than a condensate of branch points; this is a testable extension of Section IV B, not established here.","The $v_{\\rm cut}$ criterion suggests an operational definition: in any system whose retarded correlator has a discontinuity tail, an effective lightcone speed can be extracted by thresholding the normalized discontinuity at a fixed small $\\epsilon$; applying this to other channels could reveal whether $v_{\\rm cut}$ is universal or channel-dependent.","The exact positivity of all first-order coefficients for $d>1$ makes it plausible that higher-order RTA transport coefficients remain well-behaved in the massive theory, although the paper does not compute them beyond finding $\\kappa=0$ and $\\tau_\\pi=\\tau_R$ independent of mass."],"forward_implications":["For any $d\\ge 2$, charge and shear diffusion constants are non-negative for all $m/T_0$; the earlier claimed instability in [31] is an artefact of extrapolating the small-mass expansion, and the exact expressions give $D/\\tau_R\\to 1/m$ at large mass.","The sound channel is consistent with first-order hydrodynamics: the leading pole speed equals the equation-of-state speed of sound, and the attenuation $\\Gamma_s$ is non-negative, with bulk viscosity matching the Kubo formula and Chapman-Enskog results.","At finite mass the correlator's cut between $\\pm k-i/\\tau_R$ is undeformable, so the hydrodynamic pole eventually meets the cut at a finite wavenumber $k_*$ that grows with mass and dimension; in the massless theory the cut can be deformed to avoid the pole, so no such cutoff is forced.","The discontinuity profile of the cut carries causal information: its endpoints give the maximum speed, namely the speed of light, while the criterion $v_{\\rm cut}$ gives an effective lightcone that moves to smaller velocities as mass increases.","The thermoelectric transport coefficients are parameterized by a single microscopic quantity, the DC conductivity $\\sigma_Q$, with all other response fixed by thermodynamics and Ward identities."],"supporting_citations":[{"why":"Supplies the baseline thermodynamics and sound-speed result for the massive gas that the paper extends to all first-order transport coefficients.","marker":"[30]"},{"why":"The earlier claim of a mass-induced diffusive instability that the paper's exact coefficients are designed to correct.","marker":"[31]"},{"why":"Provides the variational correlator method, conventions, and massless RTA results, including thermoelectric coefficients, that the paper generalizes to finite mass.","marker":"[24]"},{"why":"Gives the massless RTA diffusion constant $D=\\tau_R/3$ and the pole-branch-cut merger that motivates the finite-mass cutoff analysis.","marker":"[6]"},{"why":"Establishes the ballistic interpretation of the RTA branch cut as a continuum of poles from particles moving at different angles, which the paper extends to the massive velocity distribution.","marker":"[10]"},{"why":"The concurrent proposal of an undeformable branch cut as a condensation of branch points, which the paper partially agrees with but reinterprets as a unique non-deformable cut.","marker":"[32]"},{"why":"Supplies the hydrodynamic dictionary used to identify $D_{\\rm sh}$ with $\\eta/(\\varepsilon_0+P_0)$ and $\\Gamma_s$ with the combination of shear and bulk viscosity.","marker":"[45]"},{"why":"Provides Chapman-Enskog results for shear and bulk viscosity against which the paper cross-checks its Kubo-formula extractions.","marker":"[46]"}],"fun_headline_variants":["Exact transport coefficients for massive RTA gas, all positive","Massive gas: closed-form transport, no negative coefficients","Causal structure from correlators: massive gas transport exact","Lightcone speed from massive gas correlator cut"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's unique straight cut rests on the assumption that the only physically admissible analytic continuation keeps both integration variables real, $u\\in[0,1]$ and $z\\in[-1,1]$; the paper notes that deformed contours may be inconsistent with causal correlators but does not prove that causality forces this particular contour.","fun_headline_variants_meta":{"raw":{"variants":["Exact transport coefficients for massive RTA gas, all positive","Massive gas: closed-form transport, no negative coefficients","Causal structure from correlators: massive gas transport exact","Lightcone speed from massive gas correlator cut"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000597,"raw_usage":{"total_tokens":2760,"prompt_tokens":879,"completion_tokens":1881,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":1815}},"tokens_in":495,"tokens_out":1881,"duration_ms":15747,"temperature":1.0,"reasoning_tokens":1815,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:33:20.158631+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the exact diffusion constants, for instance Eq. (39) in $d=3$ or Eq. (A22) in $d=2$, at large $m/T_0$: if either $D$ or $D_{\\rm sh}$ turns negative at any mass, the positivity claim fails. Separately, deform the angular integration contour $\\xi_c$ in Eq. (94) and test whether the resulting correlator satisfies the usual causal analyticity conditions; a causal deformed-contour correlator with a different cut would refute the claimed indeformability.","supporting_citations":[{"cited_title":"Quasi-particle hydrodynamics with momentum-dependent relaxation time","cited_arxiv_id":"2504.11572","evidence_quote":"The earlier claim of a mass-induced diffusive instability that the paper's exact coefficients are designed to correct."},{"cited_title":"the sound","cited_arxiv_id":null,"evidence_quote":"Gives the massless RTA diffusion constant $D=\\tau_R/3$ and the pole-branch-cut merger that motivates the finite-mass cutoff analysis."},{"cited_title":"Example 1","cited_arxiv_id":null,"evidence_quote":"Establishes the ballistic interpretation of the RTA branch cut as a continuum of poles from particles moving at different angles, which the paper extends to the massive velocity distribution."}],"review_version":2}