REVIEW 3 major objections 5 minor 5 cited by
Ultrarelativistic Fokker-Planck kinetic theory has quasi-normal modes governed by Schrödinger potentials: exact Fick diffusion plus a ballistic continuum.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 07:36 UTC pith:I5N5TOGZ
load-bearing objection 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. the 3 major comments →
The quasi-normal modes of relativistic Fokker-Planck kinetic theory
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
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
What carries the argument
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.
Load-bearing premise
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.
What would settle it
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.
If this is right
- 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.
Where Pith is reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [§V.C, Eq. (40)] 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
- [§IV.D–IV.F, Eqs. (30), (35)] 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.
- [§IV, §V] 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.
minor comments (5)
- [§III.B, Eq. (15)] 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.
- [§IV.D, Eq. (29)] 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.
- [§V.C, text before Eq. (40)] 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.
- [§V.E, caption of Fig. 3] 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'.
- [§II.B, Eq. (3)] 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.
Circularity Check
No significant circularity: the central spectral results are derived from the stated Fokker-Planck operator and standard external mathematical results; minor self-citations are not load-bearing.
full rationale
The paper's spectral derivation is self-contained and non-circular. The unitary Fokker-Planck/Schrödinger correspondence is derived in Sec. III.A (Eqs. 8-11), with the potential following from the stated FP operator (5) and dispersion relation; it is not assumed from the target spectrum. The 1D spectrum (Sec. IV) is obtained by solving the resulting ODE with matching conditions, yielding the hydrodynamic dispersion (27) and the continua (35)-(36) directly from the eigenvalue equation. The 3D k=0 result (38) is the known hydrogenic spectrum for the explicitly derived Coulomb potential, and the finite-k hydrodynamic mode is verified as an exact eigenfunction (39). No parameter is fitted to a subset of data and then renamed a prediction: β and ν are fixed model inputs, and all dispersion relations follow analytically. The self-citations (e.g., [27], [29], [33]) occur in background, in the entropy/Massieu-current discussion of Appendix A, and in stability interpretation; the FP self-adjointness used in the main text is verified explicitly in Eq. (7), so those citations are not load-bearing. The only substantive concern is a correctness/rigor gap in Sec. V.C: the Weyl-type lemma [30, Lemma 6.17] is quoted for bounded potentials, while the potential (37) contains a Coulomb singularity and, for complex χ, is non-self-adjoint; Eq. (40) may therefore lack proof. That is a mathematical-risk issue, however, not a circularity: the result is not assumed as an input or obtained by re-labeling a fit. The benchmark k=0 hydrogenic spectrum and the independent 1D ODE solution support the non-circular character of the paper's main derivations.
Axiom & Free-Parameter Ledger
axioms (4)
- domain assumption The Fokker-Planck collision operator (5) with momentum-space diffusion and detailed balance is the relevant collision model.
- domain assumption The linearized collision operator satisfies self-adjointness and positive semi-definiteness (3), with kernel span{f_eq}.
- standard math The unitary equivalence between Fokker-Planck operators and Schrödinger Hamiltonians, with inner product weight 1/f_eq.
- standard math Essential spectrum of the 3D Schrödinger operator can be characterized by the large-momentum limit of the potential (Weyl-type argument).
read the original abstract
Employing the well-known unitary equivalence between Fokker-Planck operators and Schr\"odinger Hamiltonians, we compute the quasi-normal-mode spectrum of ultrarelativistic kinetic theories with momentum-space diffusion. We show that the collision operator reduces to a Dirac-delta Schr\"odinger problem in one spatial dimension, and to a Coulomb Schr\"odinger operator with hydrogenic spectrum in three dimensions. Finite spatial wavenumber appears as a perturbation of the associated quantum potential. The hydrodynamic mode is found to obey exact Fick-type diffusion at all real wavenumbers, whereas relativistic kinematics generically produces a continuous ballistic band in the non-hydrodynamic sector, a feature absent in the Newtonian regime.
Figures
Forward citations
Cited by 5 Pith papers
-
The Lorentzian geometry of relaxation
Causality forces purely relaxational dispersion relations to follow spacelike trajectories on the Lorentzian {iω,ik} plane, producing universal bounds on diffusivity, viscosity, time-dilation deviations, and hydrodyna...
-
How acausal equations emerge from causal dynamics
A causal microscopic kinetic model can reproduce arbitrary rest-frame stable dissipative dispersion relations at real k through suitable initialization, providing a counterexample to claims that micro-causality alone ...
-
Causal UV completions of relativistic hydrodynamics
Any stand-alone dissipative hydrodynamic EFT is acausal; causal UV completions exist that preserve the exact IR hydrodynamic mode for subluminal sound speed and finite gradient-expansion radius.
-
Lorentz-boosted diffusion: initial value formulation and exact solutions
Lorentz-boosted diffusion becomes a well-posed initial-value problem on a band-limited (Paley-Wiener) function space, with an exact closed-form Shannon-Whittaker Green function.
-
Causal UV completions of relativistic hydrodynamics
Any standalone hydrodynamic EFT is acausal and requires UV completions with transient modes to restore causality.
Reference graph
Works this paper leans on
-
[1]
Dudy´ nski, Journal of Statistical Physics57, 199 (1989)
M. Dudy´ nski, Journal of Statistical Physics57, 199 (1989)
1989
-
[2]
S. Grozdanov, P. K. Kovtun, A. O. Starinets, and P. Tadi´ c, Phys. Rev. Lett.122, 251601 (2019), arXiv:1904.01018 [hep-th]
Pith/arXiv arXiv 2019
-
[3]
L. Gavassino, Phys. Rev. D110, 094012 (2024), arXiv:2408.14316 [nucl-th]
Pith/arXiv arXiv 2024
-
[4]
L. Gavassino, Astrophys. Space Sci.370, 4 (2025), arXiv:2412.00275 [astro-ph.HE]
Pith/arXiv arXiv 2025
-
[5]
P. K. Kovtun and A. O. Starinets, Phys. Rev. D72, 086009 (2005)
2005
-
[6]
G. S. Denicol, J. Noronha, H. Niemi, and D. H. Rischke, Phys. Rev. D83, 074019 (2011)
2011
-
[7]
P. Romatschke, Eur. Phys. J. C76, 352 (2016), arXiv:1512.02641 [hep-th]
Pith/arXiv arXiv 2016
-
[8]
A. Kurkela and U. A. Wiedemann, Eur. Phys. J. C79, 776 (2019), arXiv:1712.04376 [hep-ph]
Pith/arXiv arXiv 2019
-
[9]
G. D. Moore, JHEP05, 084 (2018), arXiv:1803.00736 [hep-ph]
Pith/arXiv arXiv 2018
-
[10]
L. Gavassino, Phys. Rev. Res.6, L042043 (2024), arXiv:2404.12327 [nucl-th]
Pith/arXiv arXiv 2024
-
[11]
M. Bajec, S. Grozdanov, and A. Soloviev, JHEP08, 065 (2024), arXiv:2403.17769 [hep-th]
Pith/arXiv arXiv 2024
-
[12]
R. Brants, Phys. Rev. D110, 116027 (2024), arXiv:2409.09022 [hep-th]. 3 In relativistic thermodynamics, the energy current density is defined asWµ=−Kν T µν , whereK ν is the timelike Killing vector that defines the velocity of the environment [33]. In our case, the environment is the bath, which is at rest in Minkowski space, soK ν = (1,0,0, ...). 12
Pith/arXiv arXiv 2024
-
[13]
G. S. Rocha, I. Danhoni, K. Ingles, G. S. Denicol, and J. Noronha, Phys. Rev. D110, 076003 (2024), arXiv:2404.04679 [nucl-th]
Pith/arXiv arXiv 2024
-
[14]
Dudy´ nski and M
M. Dudy´ nski and M. L. Ekiel-Jezewska, Communications in Mathematical Physics102, 17 (1985)
1985
-
[15]
L. Gavassino, Phys. Rev. D111, 103011 (2025), arXiv:2502.08740 [nucl-th]
Pith/arXiv arXiv 2025
-
[16]
A. H. Mueller, Phys. Lett. B475, 220 (2000), arXiv:hep-ph/9909388
Pith/arXiv arXiv 2000
-
[17]
J.-P. Blaizot, J. Liao, and L. McLerran, Nucl. Phys. A920, 58 (2013), arXiv:1305.2119 [hep-ph]
Pith/arXiv arXiv 2013
-
[18]
Risken,The Fokker–Planck Equation: Methods of Solution and Applications, 2nd ed., Springer Series in Synergetics, Vol
H. Risken,The Fokker–Planck Equation: Methods of Solution and Applications, 2nd ed., Springer Series in Synergetics, Vol. 18 (Springer-Verlag, Berlin, Heidelberg, 1996)
1996
-
[19]
Villani,Hypocoercivity, Memoirs of the American Mathematical Society, Vol
C. Villani,Hypocoercivity, Memoirs of the American Mathematical Society, Vol. 202 (American Mathematical Society, Providence, RI, 2009)
2009
-
[20]
Helffer and F
B. Helffer and F. Nier,Hypoelliptic Estimates and Spectral Theory for Fokker–Planck Operators, Lecture Notes in Mathe- matics, Vol. 1862 (Springer-Verlag, Berlin, Heidelberg, 2005)
2005
-
[21]
G. A. Pavliotis,Stochastic Processes and Applications: Diffusion Processes, the Fokker–Planck and Langevin Equations, Texts in Applied Mathematics, Vol. 60 (Springer, New York, 2014)
2014
-
[22]
Debbasch, K
F. Debbasch, K. Mallick, and J.-P. Rivet, Journal of Statistical Physics88, 945 (1997)
1997
-
[23]
J. Dunkel and P. H ¨anggi, Phys. Rept.471, 1 (2009), arXiv:0812.1996 [cond-mat.stat-mech]
Pith/arXiv arXiv 2009
-
[24]
B. U. Felderhof, Phys. Rev. E86, 061103 (2012)
2012
-
[25]
Pitaevskii and E
L. Pitaevskii and E. Lifshitz,Physical Kinetics, v. 10 (Elsevier Science, 2012)
2012
-
[26]
Cercignani and G
C. Cercignani and G. M. Kremer,The relativistic Boltzmann equation: theory and applications(2002)
2002
-
[27]
Gavassino, (2026), arXiv:2601.03081 [gr-qc]
L. Gavassino, (2026), arXiv:2601.03081 [gr-qc]
Pith/arXiv arXiv 2026
-
[28]
M. P. Heller, A. Serantes, M. Spali´ nski, and B. Withers, Phys. Rev. Lett.130, 261601 (2023), arXiv:2212.07434 [hep-th]
Pith/arXiv arXiv 2023
-
[29]
L. Gavassino, Phys. Lett. B840, 137854 (2023), arXiv:2301.06651 [hep-th]
Pith/arXiv arXiv 2023
-
[30]
Teschl,Mathematical Methods in Quantum Mechanics With Applications to Schr ¨odinger Operators”, Graduate Studies in Mathematics, Vol
G. Teschl,Mathematical Methods in Quantum Mechanics With Applications to Schr ¨odinger Operators”, Graduate Studies in Mathematics, Vol. 99 (American Mathematical Society, Providence, RI, 2009)
2009
-
[31]
Kato, Progress of Theoretical Physics4, 514 (1949)
T. Kato, Progress of Theoretical Physics4, 514 (1949)
1949
-
[32]
Kato,Perturbation Theory for Linear Operators, reprint of the 1980 edition ed., Classics in Mathematics (Springer, 1995)
T. Kato,Perturbation Theory for Linear Operators, reprint of the 1980 edition ed., Classics in Mathematics (Springer, 1995)
1980
-
[33]
Gavassino, Classical and Quantum Gravity38, 21LT02 (2021), arXiv:2104.09142 [gr-qc]
L. Gavassino, Classical and Quantum Gravity38, 21LT02 (2021), arXiv:2104.09142 [gr-qc]
Pith/arXiv arXiv 2021
-
[34]
Landau and E
L. Landau and E. Lifshitz,Statistical Physics, v. 5 (Elsevier Science, 2013)
2013
-
[35]
L. Gavassino, M. Antonelli, and B. Haskell, Phys. Rev. Lett.128, 010606 (2022), arXiv:2105.14621 [gr-qc]
Pith/arXiv arXiv 2022
-
[36]
G. Soares Rocha, L. Gavassino, and N. Mullins, Phys. Rev. D110, 016020 (2024), arXiv:2405.10878 [nucl-th]
Pith/arXiv arXiv 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.