{"id":"3b1a374b-ae90-405e-a39f-7bb7ea597aa1","arxiv_id":"2601.19474","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The quasi-normal spectrum of ultrarelativistic Fokker-Planck kinetic theory consists of an exact diffusive hydrodynamic mode, continuous ballistic bands, and a hydrogenic discrete tower in three dimensions.","lead":"This paper calculates the full set of decay modes for a relativistic gas of particles that diffuse in momentum through many small-angle collisions. It shows that the slowest mode is always a simple diffusion mode, while faster modes form a continuous band that can travel ballistically, unlike in nonrelativistic gases.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"3D continuous-spectrum result (Eq. 40) rests on an unjustified extension of Weyl's lemma to a singular, non-self-adjoint potential; if the essential spectrum differs, the claimed 3D ballistic band is unproven.","rationale":"After a careful reading, the 1D analysis (Sec. IV) is internally consistent and appears rigorous: the bound-state matching yields the exact hydrodynamic mode (27), and the piecewise-oscillatory/hybrid solutions account for the continuous branches for both imaginary and real k. The zero-wavenumber 3D result (Sec. V.A) is just hydrogen, and the 3D hydrodynamic eigenfunction (39) is verified directly. The inner-product typo in Eq. (2) is real but evidently a missing 1/f_eq in the printed measure; Eq. (7) and the subsequent unitary transform read correctly once this is fixed, so it is not load-bearing. The one place where the strongest claim outruns its proof is Sec. V.C: the Weyl characterization is invoked for a potential with a Coulomb singularity and a non-self-adjoint perturbation, and the cited lemma does not cover this case. Since Eq. (40) is the only basis for the claimed 3D ballistic continuum, this missing justification is the single most load-bearing concern. It does not invalidate the rest of the paper, but it warrants a conditional acceptance pending a proof or a counterexample.","tokens_in":12647,"tokens_out":20247,"duration_ms":224522,"concrete_test":"Prove that the Coulomb term -β/(2p) is a relatively compact perturbation of -1/2Δ in L²(R³) (e.g., by verifying the resolvent-compactness or form-bound criterion), then compute σ_ess of the limit operators -1/2Δ + β²/8 + χ cosθ for each direction θ; if the relative compactness holds and the non-self-adjoint perturbation has the expected limit-operator essential spectrum, Eq. (40) follows. Equivalently, construct explicit Weyl sequences localized near |p|=R in direction θ, showing ||(H - E)ψ_{R,θ}||→0 for every E in the claimed set and that no such sequence exists for E outside it. If the Coulomb term cannot be dropped, identify a point outside the rectangle with a Weyl sequence, which would refute Eq. (40).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that relativistic kinematics produces a continuous ballistic band in 3D depends entirely on Eq. (40), derived in Sec. V.C. The argument cites a Weyl-type lemma ([30, Lemma 6.17]) that is stated for Schrödinger operators with bounded potentials. The potential in Eq. (37), Vχ=β²/8 - β/(2p)+χ p1/p, is not bounded: the Coulomb term -β/(2p) is singular at p=0, and for real k (the physically relevant case for the ballistic band) χ=ik/(2ν) is purely imaginary, so the operator is non-self-adjoint. The paper does not prove that the essential spectrum is unaffected by the singular local term, nor that the Weyl construction survives non-self-adjointness. Kato's theorem cited in Sec. V.D only gives an inclusion (spectrum lies within |χ| of unperturbed spectrum), not the equality in Eq. (40). If the actual essential spectrum differs (e.g., is larger, or contains residual spectrum), the claimed rectangle ω∈k[-1,1]-i[β²ν/4,∞) is not the spectrum, and the 3D ballistic-band conclusion fails. The 1D spectrum and the zero-wavenumber hydrogenic levels are unaffected, so the paper's other results stand; this is a correctness risk isolated to the most novel three-dimensional claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the well-known ground-state (unitary) equivalence between Fokker-Planck operators and Schrödinger operators to compute the quasi-normal-mode spectrum of ultrarelativistic kinetic theory with momentum-space diffusion. In one spatial dimension the effective potential is a Dirac-delta-plus-step problem; the author derives an exact hydrodynamic diffusive mode ω = -i k^2/(β^2ν) and continuous nonhydrodynamic branches. In three dimensions the k=0 problem is a hydrogen-like Coulomb Hamiltonian with discrete levels and a continuum; the hydrodynamic mode is verified exactly for all k, and the author claims a finite-k continuous ballistic band ω ∈ k[-1,1] - i[β^2ν/4,∞). The remainder of the discrete spectrum is constrained by a bounded-perturbation argument. The paper is analytic and has no fitted parameters.","tokens_in":12969,"tokens_out":15628,"duration_ms":168976,"significance":"If the spectral claims are correct, the paper provides an unusually clean exact mapping between relativistic kinetic theory and solvable quantum mechanics, giving nontrivial predictions about transient dynamics: exact Fick-type diffusion for the hydrodynamic mode and a ballistic continuum in the nonhydrodynamic sector that is absent in the Newtonian regime. The 1D analysis is self-contained and the k=0 3D hydrogenic result is a strong benchmark. The main added value over existing RTA kinetic computations is the closed-form spectral description for a Fokker-Planck collision kernel, and the explicit demonstration of relativistic kinematic effects on the nonhydrodynamic sector. The principal vulnerability is the 3D finite-k continuous-spectrum derivation, which is the only support for the paper's most novel three-dimensional claim.","major_comments":[{"comment":"The derivation of the 3D continuous spectrum cites [30, Lemma 6.17], but that Weyl-type argument is stated for Schrödinger operators with bounded potentials. The potential in Eq. (37) contains the singular Coulomb term -β/(2p); moreover, for real k (exactly the case needed for the ballistic band), χ = ik/(2ν) is purely imaginary, so the Hamiltonian is non-self-adjoint. The text asserts that repeating the direction-dependent construction 'reconstructs the full continuous spectrum', but no proof is given for the converse inclusion or for the absence of residual/other essential spectrum in the non-self-adjoint case. Since Eq. (40) is the basis of the central claim of a 3D ballistic continuum, this gap is load-bearing. Please either supply a rigorous reference or a self-contained proof for the essential spectrum of singular non-self-adjoint Schrödinger operators with direction-dependent limi","section":"§V.C, Eq. (40)"},{"comment":"For k∈iR, Eq. (30) states a continuum with iω ≥ β²ν/4 + |k|, while Eq. (35) gives the lower edge β²ν/4 - |ik|. The two inequalities differ by the sign of the wavenumber term. I understand that Eq. (30) describes the purely oscillatory sector and Eq. (35) includes the hybrid bands of §IV.E, but this is not explained. The present wording makes the 1D spectrum look internally inconsistent. Please add a sentence clarifying that (30) is a subset statement and that (35) is the full nonhydrodynamic continuum for imaginary wavenumbers.","section":"§IV.D–IV.F, Eqs. (30), (35)"},{"comment":"The paper calls the continuous branches 'quasi-normal modes' even though the associated wavefunctions are not L^2 (they remain oscillatory at infinity or decay only on one half-line). This is a standard convention in the relativistic kinetic/QNM literature, but it should be stated explicitly because the Hilbert-space spectrum of the transformed operator would not contain these generalized eigenfunctions. A precise definition of the spectral notion used in Eqs. (36) and (40) would help readers distinguish exact eigenmodes from approximate/generalized modes.","section":"§IV, §V"}],"minor_comments":[{"comment":"The D=1 delta-potential limit is clear, but the distributional meaning of ∂v/∂p for the massless velocity field v = sgn(p) is only implicit. A one-sentence remark that Eq. (16) is a distributional identity would improve rigor.","section":"§III.B, Eq. (15)"},{"comment":"The text immediately after Eq. (29) says 'for any (s1,s2)∈R² there exists a choice of coefficients.' It may be helpful to give the explicit matching-coefficient construction, since the existence statement is essential for the completeness of the continuous band.","section":"§IV.D, Eq. (29)"},{"comment":"The phrase 'this result remains valid even for complex χ' is asserted without qualification. Since the case χ∈iR is non-self-adjoint, this is precisely the point that needs proof; please at least flag it as an assumption here if a full proof is deferred.","section":"§V.C, text before Eq. (40)"},{"comment":"The phrase 'vertical rectangle with upper corners (±k,-igap)' is imprecise: for real k the set in Eq. (40) is an infinite downward half-strip, not a finite rectangle. The caption should say 'half-strip' or 'infinite vertical strip'.","section":"§V.E, caption of Fig. 3"},{"comment":"The property IΦ=0 iff Φ∝feq is stated as a general physical requirement, but its proof in Appendix A relies on the absence of additional conserved currents and on PT symmetry. This is reasonable but might be marked more clearly as an assumption about the model class.","section":"§II.B, Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well written and the 1D and k=0 3D results are likely solid. The main risk is isolated to Sec. V.C: the 3D finite-k continuous spectrum—the paper's most novel claim—rests on an unproven extension of a bounded-potential Weyl argument to a singular, non-self-adjoint Hamiltonian. If the author can supply a convincing proof or a directly applicable reference for this spectral statement, I would be inclined to accept; without it, the central 3D conclusion remains unsubstantiated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does two things cleanly. First, it gives the complete quasi-normal spectrum for ultrarelativistic Fokker-Planck kinetics in one spatial dimension: the exact hydrodynamic mode ω=-i k²/(β²ν) plus the two continua Re ω = ±k, Im ω ≤ -β²ν/4, derived by explicit matching conditions. That part is careful and, as far as I can tell, correct. Second, it shows that the zero-wavenumber 3D problem maps to a hydrogenic Schrödinger equation, recovering the known spectrum, and proves that the hydrodynamic mode stays exactly Fickian at all real k. The diffusion coefficient is the same in one and three dimensions, which is a nice nontrivial consistency. No fitting, no numerology; everything follows analytically from the fixed inputs.\n\nThe soft spot is the 3D continuous spectrum at finite wavenumber, Eq. (40). The argument in Sec. V.C invokes a Weyl lemma stated for bounded potentials, but the potential here has a Coulomb singularity at p=0 and, for real k, the χp1/p term is purely imaginary, so the operator is not self-adjoint. The paper does not show that the lemma extends to this setting. Kato's bounded-perturbation bound gives an inclusion, not the equality claimed in Eq. (40). So the ballistic band Re ω ∈ [-k,k], Im ω ∈ (-∞, -β²ν/4] is plausible—likely true physically—but it is not proven by the text as written. This is the most novel result of the paper, so that gap matters. The 1D spectrum and the k=0 hydrogenic part are unaffected.\n\nMinor issues: the two illustrative 3D numerical modes have no code or data; that's fine for a theory paper but worth noting. Also, the claim that for imaginary k the 3D spectrum reproduces the 1D result exactly is asserted, not demonstrated.\n\nWho is this for? People working on spectral approaches to hydrodynamization, heavy-ion transport models, and anyone interested in non-self-adjoint Fokker-Planck operators. It is clearly written and the physical interpretation of the ballistic band (velocity saturation protecting high-energy transport) is compelling. The author is in control of the background literature; the citation pattern is appropriate, and the self-citations are standard and relevant. This deserves a serious referee. The referee should push on the essential-spectrum argument and ask for a correct proof or an explicit caveat. If that is fixed, I'd take the paper.","headline":"The 1D spectrum and exact diffusive mode are solid and new; the 3D ballistic continuum is plausible but the Weyl argument doesn't carry the weight.","tokens_in":13448,"tokens_out":3883,"would_cite":true,"duration_ms":40251,"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":"Ultrarelativistic Fokker-Planck kinetic theory has quasi-normal modes governed by Schrödinger potentials: exact Fick diffusion plus a ballistic continuum.","keywords":["quasi-normal modes","Fokker-Planck kinetic theory","relativistic kinetic theory","Schrödinger equivalence","hydrodynamic mode","ballistic continuum","hydrogenic spectrum","momentum-space diffusion"],"falsifier":"Compute the essential spectrum of the operator −1/2 Δ + β²/8 − β/(2p) + χ p₁/p on L²(R³) for complex χ (the problem stated in Section V.C), and check whether it equals ∪_θ (χ cosθ + [β²/8,∞)). Any deviation — a missing segment of the band or an additional continuous component — would refute equation (40). A second check: verify numerically whether the exact eigenfunction exp(−βp/2 − 2χp₁/β) remains square-integrable for |Re χ|>β²/4, which would falsify the bound |Im k|<β²ν/2.","tokens_in":12525,"feed_emoji":"⚛️","tokens_out":5943,"duration_ms":63335,"temperature":0.7,"pith_summary":"This paper tries to establish that the quasi-normal modes of an ultrarelativistic gas undergoing Fokker-Planck momentum diffusion are exactly the spectral states of two familiar quantum problems. In one spatial dimension the collision operator becomes a Dirac-delta Schrödinger Hamiltonian; in three dimensions it becomes a Coulomb (hydrogen-like) Hamiltonian. From this map the paper derives the full spectrum: a unique hydrodynamic mode with exact Fick dispersion ω=−i k²/(β²ν) at every real wavenumber, and a non-hydrodynamic continuum that propagates ballistically, with real frequency spanning ±k in one dimension and within [−k,k] in three. The result matters because it shows that in relativistic kinetics diffusion and ballistic transport coexist, and it gives a rare kinetic model where the hydrodynamic mode is exactly diffusive at all wavelengths.","feed_headline":"Relativistic gas: exact Fick diffusion plus ballistic modes","feed_subtitle":"A quantum-map trick turns the kinetic problem into delta and Coulomb potentials, predicting coexisting diffusion and ballistic transport.","key_machinery":"The load-bearing tool is the unitary equivalence between a detailed-balanced Fokker-Planck generator and a Schrödinger operator in momentum space: writing Φ = sqrt((2π)^D f_eq) ψ turns the kinetic equation into −1/2 ∂²ψ + V_χ(p)ψ = Eψ, with E = iω/(2ν). The effective potential encodes both the equilibrium state and the wavenumber perturbation; finite wavenumber k acts as a linear perturbation χ v_1 that breaks rotational symmetry. This map converts the quasi-normal-mode problem into a quantum spectral problem whose discrete and continuous spectra can be read off from known delta and Coulomb Hamiltonians.","core_discovery":"The central claim is that the linearized Vlasov-Fokker-Planck equation, after rescaling and a unitary ground-state transform, is equivalent to the time-independent Schrödinger equation with effective potential V_χ(p)=β²v²/8 − (β/4)∂_j v_j + χ v_1, where χ = ik/(2ν). For massless particles this collapses to a delta potential in one dimension and a 1/r Coulomb potential in three dimensions. Consequently the zero-wavenumber spectrum in 3D is hydrogenic: iω_N=(β²ν/4)(1−1/N²) with fourfold N=2 degeneracy, plus a continuum above β²ν/4. At finite wavenumber the hydrodynamic mode remains exactly ω=−i k²/(β²ν) for all real k, existing only for |Im k|<β²ν/2 and merging into the continuum beyond that s","pith_inferences":["Extension: the exact-diffusion result suggests that in this class of theories a hydrodynamic description is not an approximation at long wavelengths but an exact sector of the linear response, which may place strict constraints on resummed hydrodynamic frameworks.","Extension: since the ballistic continuum arises from large-momentum velocity conservation, one could test the prediction in weakly coupled gluon plasmas by looking for propagating, damped non-hydrodynamic signals in early-time far-from-equilibrium dynamics.","Extension: determining whether the discrete hydrogenic levels persist under the complex dipolar perturbation for |Im k| below β²ν/2, or convert into resonances, is a calculable problem that would sharpen the analyticity claims made near χ=0.","Extension: the 1D result that the hydrodynamic mode exists only for |Im k|<β²ν/2 may be interpreted as a spectral criterion for when a diffusive effective field theory can be matched to the underlying kinetic theory."],"forward_implications":["In one and three spatial dimensions, the hydrodynamic mode of ultrarelativistic Fokker-Planck kinetics is exactly diffusive at all real wavenumbers, with diffusion coefficient D=(β²ν)⁻¹; no higher-order corrections appear.","Non-hydrodynamic modes are not purely damped: they form continuous branches with real frequencies within ±k (1D: exactly ±k), showing ballistic transport survives stochastic momentum diffusion in the relativistic regime.","The diffusive mode exists only inside the strip |Im k|<β²ν/2; outside it merges into the continuum, so the theory never violates the causality bound Im ω ≤ |Im k|.","At zero wavenumber in 3D, the relaxation spectrum is the hydrogen spectrum, giving an infinite tower of discrete modes (slowest non-hydrodynamic rate 3β²ν/16) plus a continuum above β²ν/4.","Because the 3D perturbation is dipolar, the zero-wavenumber degeneracies are lifted at finite k in the same way a Stark field splits hydrogen levels."],"fun_headline_variants":["Quantum map gives exact diffusion and ballistic modes","Fokker-Planck gas maps to delta and Coulomb potentials","Relativistic kinetics: hydrogenic spectrum and exact Fick law","Exact Fick diffusion plus ballistic modes from a quantum trick"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The 3D continuous spectrum is derived in Section V.C by applying a Weyl-type lemma stated for bounded potentials to a Coulomb potential with a 1/p singularity and a complex perturbation, an extension the paper does not justify; if the true essential spectrum differs, the claimed ballistic band ω∈k[−1,1]−i[β²ν/4,∞) may be wrong.","fun_headline_variants_meta":{"raw":{"variants":["Quantum map gives exact diffusion and ballistic modes","Fokker-Planck gas maps to delta and Coulomb potentials","Relativistic kinetics: hydrogenic spectrum and exact Fick law","Exact Fick diffusion plus ballistic modes from a quantum trick"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000301,"raw_usage":{"total_tokens":1549,"prompt_tokens":699,"completion_tokens":850,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":443,"completion_tokens_details":{"reasoning_tokens":781}},"tokens_in":443,"tokens_out":850,"duration_ms":9188,"temperature":1.0,"reasoning_tokens":781,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T07:36:32.519506+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the essential spectrum of the operator −1/2 Δ + β²/8 − β/(2p) + χ p₁/p on L²(R³) for complex χ (the problem stated in Section V.C), and check whether it equals ∪_θ (χ cosθ + [β²/8,∞)). Any deviation — a missing segment of the band or an additional continuous component — would refute equation (40). A second check: verify numerically whether the exact eigenfunction exp(−βp/2 − 2χp₁/β) remains square-integrable for |Re χ|>β²/4, which would falsify the bound |Im k|<β²ν/2.","supporting_citations":[],"review_version":1}