REVIEW 3 major objections 5 minor 45 references
On quantitative linear gravitational relaxation
T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Small smooth perturbations of compact galactic equilibria in the gravitational Vlasov-Poisson system decay at quantified rates: the gravitational force decays as a power of time set by the regularity of the data.
desk verdict A novel and technically demanding proof of quantitative decay for the linearized gravitational Vlasov-Poisson system; the main soft spot is the imported spectral input that needs referee verification. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Two coordinate changes and one decomposition carry the argument. Action-angle variables $(\theta,I)$ convert the transport operator into $\omega(I)\partial_\theta$ with $\omega(I)=1/T(I)$, so pure transport is solved by oscillating phases $e^{-2\pi i m t\omega(I)}$. The new foliation $(R,y,z)=(R,\omega(I),E-\Psi_L(R))$ separates the singular behaviour at the trapping set $E=E_L^{\min}$ from the vacuum boundary, and is invertible when $\partial_L\omega=O(\eta)$. Near the elliptic points, the Birkhoff-Poincaré normal form (a canonical Hamiltonian normal form around a nondegenerate elliptic point) gives explicit action-angle expansions, and the Green's function decomposes as $S_m(\theta(R,I))=E^{-|m|/2}H_m(R,y,z)$ with $H_m$ analytic; this factorisation cancels the singular factor in the Fourier coefficients of the data and yields the quantitative derivative bound (1.32). The full linearised decay is obtained from Stone's formula by proving uniform resolvent estimates, based on a limiting absorption principle and on Plemelj-type estimates split into near-resonant and non-resonant frequencies.
What would settle it
An eigenvalue inside the essential spectrum of the linearised operator for any steady state in the stated family would be a direct obstruction: the Stone-formula representation would then contain a non-decaying oscillatory component, and the gravitational force could not satisfy the claimed $(1+t)^{-b}$ bound uniformly over the full range of $b\le K-1$. A concrete way to look for it is to scan polytropic parameters near the thresholds $\mu=2$, $\nu=1$ or at small but nonzero $\eta$ for resonant modes.
Extended reading notes
Core claim
The paper's central claim is Theorem 1.1. For polytropic steady states of the form (1.5) with $\mu>2$, $\nu>1$ and amplitude $\eta$ sufficiently small, any initial datum $f_0\in C^k(\Omega)$ satisfying the orthogonality condition (1.19) produces a unique solution of the linearised problem (1.13)–(1.14) whose gravitational force obeys $\|\partial_R U_{|\varphi'|f}(t,\cdot)\|_{C^{K-1-b}}\le C\|f_0\|_{C^K(\Omega)}(1+t)^{-b}$ for $0\le b\le K-1$, with $K=\min\{\mu-1,\nu,k\}$. This bound applies to the full linearised flow, not only to the pure transport part: the gravitational back-reaction is shown to contribute a strictly lower-order term, so the transport mechanism dictates the decay. The theorem also implies, for $k\ge 2$, decay of the spatial density in $C^{K-2-b}$, and, for $K\ge 3$, existence of a scattering profile in the natural weighted Hilbert space $H$.
Load-bearing premise
The argument depends on the spectral input that the linearised operator has no embedded eigenvalues, together with the geometric input that the frequency map $\omega(E,L)$ is strictly monotone in $E$ with $\partial_L\omega=O(\eta)$, so that the $(y,z)$-foliation is a true coordinate change; if an embedded eigenvalue existed or the foliation degenerated, the Stone-formula and Plemelj machinery would not deliver the decay estimate.
Editorial extensions
If this is right
- The theorem gives explicit, regularity-determined decay rates: smoother steady states and smoother initial data yield faster algebraic decay, with the exponent limited by $K=\min\{\mu-1,\nu,k\}$.
- For data with $k\ge 2$, the spatial density inherits the decay in $C^{K-2-b}$, so macroscopic observables relax at quantified rates rather than merely oscillating.
- For $K\ge 3$, the linearised solution scatters to a well-defined profile in the weighted space $H$, which is the linear precursor of a relaxation picture for nearby galaxies.
- The result is new even in the pure transport case $\eta=0$; the mechanism that converts regularity into decay is tied to the point-mass potential and to the three-dimensional geometry, so the same rates should not be expected in one-dimensional confining models.
- Quantitative linear decay is an essential stepping stone to nonlinear asymptotic stability, since the same cancellations near trapping control the quadratic terms in a bootstrap argument.
Reading between the lines
- An extension the paper leaves implicit is that the decay mechanism should be stable under small perturbations of the steady state that preserve the single-gap spectrum: the proof's central spectral input is absence of embedded eigenvalues, so nearby equilibria in the same regularity class should satisfy the same bounds with modified constants.
- The factored Green's function identity suggests a sharper heuristic: the effective transfer function is nearly diagonal in the new $(y,z)$ foliation, so a numerical scheme that discretises $(y,z)$ rather than $(E,L)$ should resolve the decay with fewer modes.
- A testable extrapolation is that violating the smallness of $\eta$ or the monotonicity of $\omega$ will destroy or slow the decay; constructing an equilibrium where $\partial_L\omega$ is not small and measuring the gravitational force would probe the sharpness of the foliation condition.
- The method points toward inhomogeneous plasma analogues, where BGK-type equilibria with a trapping region could be treated by the same combination of normal form and resolvent estimates, provided an absence-of-eigenvalues statement is available.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies quantitative decay for the linearized gravitational Vlasov-Poisson system around compactly supported polytropic steady states with a central point mass. Theorem 1.1 claims that, for small steady-state amplitude eta, initial data satisfying the orthogonality condition (1.19), and regularity parameter K = min{mu-1, nu, k}, the gravitational force decays as (1+t)^{-b} for 0 <= b <= K-1, with the decay rate improving with the regularity of the data. The proof combines action-angle variables, a Birkhoff-Poincare normal form, a new two-dimensional foliation of the action support, a quantitative pure-transport decay theorem, Stone's formula, and uniform resolvent bounds with Plemelj-type singular integral estimates. The paper is self-contained from Section 3 onward in the sense of proving the new estimates, but it relies on prior results for the spectral structure of the linearized operator and for regularity of action-angle variables.
Significance. If correct, the result is a significant advance: it gives the first quantitative linear decay theorem around nontrivial inhomogeneous gravitational equilibria with trapping, and it introduces a genuinely new mechanism (the (y,z)-foliation and the normal-form cancellations) that is likely to be reusable in related problems such as massive perturbations of the Schwarzschild solution. The central claim is a genuine theorem, not an identity: no parameter is fitted, and the decay rate is derived rather than imposed. The authors are also honest about the main external inputs: the absence of embedded eigenvalues and the spectral description are quoted from [22,44], and several bootstrap steps are compressed. My assessment is that the main architecture is coherent and the new contributions are substantial, but two load-bearing points need to be made fully explicit before the theorem can be considered established: the exact scope of the imported spectral theorem, and the global invertibility of the foliation change of variables.
major comments (3)
- [Section 2.3, Proposition 2.5 and Eq. (2.29)] Proposition 2.5 is load-bearing: the claims sigma(L)=sigma_ess(L)=sigma(T) with no embedded eigenvalues are used in Stone's formula (2.38), in the definition of sigma_L in (2.39), and in every resolvent bound, in particular Theorems 6.13 and 6.14. The proof of Proposition 2.5 consists of a citation to [22,21,44], but the manuscript does not verify that the theorems of [44] apply verbatim to the weighted space H in (1.12), to the polytropic exponents mu>2, nu>1, and to the small-eta regime used here. This is not merely a cosmetic gap: if an embedded eigenvalue existed, the resolvent bounds (6.48) and (6.57) would acquire point-spectrum contributions and the decay (1.20) would fail. Please either state precisely the theorems being cited and check their hypotheses, or give a self-contained argument for the spectral statement in this setting.
- [Section 3.3, Lemma 3.6 and Remark 3.10] Lemma 3.6 proves only that the Jacobian determinant in (3.38) is nonzero at each point, which gives local invertibility. The global statement that the map (3.36) is a diffeomorphism onto J_R is asserted in Remark 3.10 ('the map ... is a bijection') and is used throughout the paper to justify the change of variables in the Plemelj integrals, e.g. in (5.5), (6.7), and (6.11). The missing step is global injectivity together with a description of the boundary of J_R. For fixed R, injectivity should follow from monotonicity of y=omega(E,L) in L, since d_L y = P_R omega > 0, but the proof is not written. Because all subsequent estimates are formulated in (R,y,z) coordinates, this global fact must be proved explicitly rather than left as an assertion.
- [Section 7.1, Eq. (7.14), and Lemma 7.2] The bound (7.14) is introduced with the phrase 'by the same argument as in Theorem 6.14' and is then used to justify integration by parts in lambda and the passage to the limit in (7.18). This is a load-bearing estimate: it supplies the |lambda|^{-2} decay that produces the (1+t)^{-b} factor for b>=1. It should be proved explicitly, because the objects being estimated, R^{+,lambda,epsilon}_{a,b}[U^+_epsilon] - R^{-,lambda,epsilon}_{a,b}[U^-_epsilon], mix differences of resolvents and differences of Plemelj kernels. In the same lemma, the inductive proof of (7.11) says 'arguing as in (7.8)' and omits the details of the epsilon -> 0 limit for the higher lambda-derivatives; since (7.11) is used to kill boundary terms in the integration by parts, the argument needs to be written out.
minor comments (5)
- [Throughout, esp. Eq. (6.27)] The notation in (6.27)-(6.28), written as omega(I)+lambda+i epsilon/(2 pi m), is ambiguous: it should be omega(I) + (lambda + i epsilon)/(2 pi m), with the parentheses displayed explicitly. The same ambiguity appears in several later formulas.
- [Theorem 5.1, Eq. (5.9)] The proof of Theorem 5.1 passes an infinite sum over m through the integration by parts in y without an explicit summability statement. The estimate (5.9) suggests the sum is finite, but a sentence justifying the interchange would improve readability.
- [Remark 7.4] In the definition of Sigma^{pm}, the indicator chi_{[omega(R,z), omega(R,z)]}(y) appears to denote a degenerate interval and is not used in the displayed trace formula. This should be corrected or deleted.
- [Section 4.3, Proposition 4.6] The symbol C is reused with different meanings in Proposition 4.6, e.g. in 'C_{k,j} C^{|m|} |m|^{k+j}'. The constants are independent of m, but the exposition would be clearer if the various constants were labeled distinctly or if the dependence on k,j was stated once.
- [Abstract and Introduction] The abstract says 'first linear asymptotic stability result around such equilibria', while Theorem 1.1 proves decay of the gravitational force and density, and Corollary 7.3 provides H-scattering only under the extra assumption K>=3. The wording in the abstract is stronger than the proved statements; I suggest qualifying it as decay of macroscopic quantities rather than full asymptotic stability.
Circularity Check
No circularity: the quantitative decay is derived from resolvent and Plemelj estimates, with published spectral inputs that are independent of the target bound.
full rationale
I walked the derivation chain from the linearised system (2.31) through Stone's formula (2.38), the (y,z) foliation (3.36), the product and Plemelj estimates of Sections 4 and 6, and the final integration-by-parts argument in Section 7. The claimed decay rate (1.20) is not an input: the exponent K is the regularity index min{mu-1,nu,k}, and the proof genuinely derives (1+t)^{-b} from derivative counts, boundary vanishing of q, invertibility of F^{+-eps}_lambda, and the resolvent bounds (6.48), (6.57). No parameter is fitted to a subset of data and then renamed a prediction. The main load-bearing import is the purely spectral statement of Proposition 2.5, citing [22, 44] for absence of embedded eigenvalues and the essential-spectrum description. That is a separate, prior spectral fact with stated assumptions that do not include the decay conclusion; it is independent support rather than a circular reduction. The skeptical concerns about whether [44, Theorems 6.4.1 and 6.5.5] apply verbatim to the weighted space H, and about whether Lemma 3.6's local invertibility plus Remark 3.10's bijection assertion give global diffeomorphism properties, are correctness/verification risks, not instances of the paper's derivation reducing to its own inputs. I found no equation or fitted quantity that is equivalent by construction to the output of Theorem 1.1, and no self-definitional step.
Assumptions & free parameters
assumptions (7)
- domain assumption Polytropic steady-state ansatz (1.5): f = eta (E0-E)_+^mu (L-L0)_+^nu with mu>2, nu>1 and phi'<0.
- domain assumption Single-gap condition on kappa (1.9): -2 - (2/3)(M^2/(2L0)) < kappa < 0.
- domain assumption Smallness eta in [0,eta0) with eta0 sufficiently small.
- domain assumption Initial datum f0 orthogonal to the kernel of L, equivalent to zero angular average (1.19), (2.27).
- domain assumption Spectral facts from prior work: sigma(L)=i(-infinity,-lambda_min] union {0} union i[lambda_min,infinity) and no embedded eigenvalues (Proposition 2.5, cited [22,44]).
- standard math Birkhoff-Poincare canonical normal form theorem for planar Hamiltonian systems [35, Theorem A].
- domain assumption Steady-state existence and support properties from [44, Proposition 6.2.6].
Cite this review
Pith. "Pith review of On quantitative linear gravitational relaxation." pith.science (2026). https://pith.science/paper/TRDB2LYF
@misc{pith2026250514856,
author = {Pith},
title = {Pith review of: On quantitative linear gravitational relaxation},
year = {2026},
howpublished = {\url{https://pith.science/paper/TRDB2LYF}},
note = {Machine review of arXiv:2505.14856}
}
read the original abstract
We prove quantitative decay rates for the linearised Vlasov-Poisson system around compactly supported equilibria. More precisely, we prove decay of the gravitational potential induced by the radial dynamics of this system in the presence of a point mass source. Our result can be interpreted as the gravitational version of linear Landau damping in the radial setting and hence the first linear asymptotic stability result around such equilibria. We face fundamental obstacles to decay caused by the presence of stable trapping in the problem. To overcome these issues we introduce several new ideas. We use different tools, including the Birkhoff-Poincar\'e normal form, action-angle type variables, and delicate resolvent bounds to prove a suitable version of the limiting absorption principle and obtain the decay-in-time.
Figures
Reference graph
Works this paper leans on
-
[44]
Straub, C., Pulsating Galaxies , Doctoral thesis, Universit¨ at Bayreuth 2024
work page 2024
-
[1]
A., Remarks on the problems of stability in stellar dynamics
Antonov, V. A., Remarks on the problems of stability in stellar dynamics. Soviet Astronom. AJ 4, 859–867 (1960)
work page 1960
-
[2]
Y., On algebraic damping close to inhomogeneous Vlasov equilibria in multi- dimensional spaces
Barr´ e, J., Yamaguchi, Y. Y., On algebraic damping close to inhomogeneous Vlasov equilibria in multi- dimensional spaces. J. Phys. A: Math. Theor. 46 225501 (2013)
work page 2013
-
[3]
Bedrossian, J., Masmoudi, N., Inviscid damping and the asymptotic stability of planar shear flows in the 2D Euler equations. Publ. Math. Inst. Hautes ´Etudes Sci. 122, 195–300 (2015)
work page 2015
-
[4]
Bedrossian, J., Masmoudi, N., Mouhot, C., Landau Damping: Paraproducts and Gevrey Regularity. Ann. PDE 2, Art. 4, (2016)
work page 2016
-
[5]
Bedrossian, J., Masmoudi, N., Mouhot, C., Landau damping in finite regularity for unconfined systems with screened interactions. Comm. Pure Appl. Math. 71, 537–576 (2018)
work page 2018
-
[6]
Bedrossian, J., Masmoudi, N., Mouhot, C., Linearized wave-damping structure of Vlasov-Poisson in R3. SIAM J. Math. Anal. 54, 4379–4406 (2022)
work page 2022
-
[7]
Binney, J., Tremaine, S., Galactic Dynamics (second edition), Princeton Series in Astrophysics 4, Princeton University Press 2008
work page 2008
Show all 45 references
-
[8]
Chaturvedi, S., Luk, J., Phase mixing for solutions to 1D transport equation in a confining potential. Kinet. Relat. Models 15, 403–416 (2022)
2022
-
[9]
Preprint, arXiv:2409.14626
Chaturvedi, S., Luk, J., Linear and nonlinear phase mixing for the gravitational Vlasov-Poisson system under an external Kepler potential. Preprint, arXiv:2409.14626
-
[10]
M., Savits, T
Constantine, G. M., Savits, T. H., A generalised Fa` a di Bruno formula with applications,Trans. AMS 348 (1996), 503–520
1996
-
[11]
Cvijovi´ c, D., New identities for the partial Bell polynomials, Appl. Math. Lett. 24, 1544–1547, (2011)
2011
-
[12]
Despr´ es, B., Scattering structure and Landau damping for linearized Vlasov equations with inhomo- geneous Boltzmannian states. Ann. Henri Poincar´ e20, 2767–2818 (2019)
2019
-
[13]
Special issue in memory of Walter Craig,J
Faou, E., Horsin, R., Rousset, F., On linear damping around inhomogeneous stationary states of the Vlasov-HMF model. Special issue in memory of Walter Craig,J. Dyn. Diff. Equat. 33,1531-1577 (2021)
2021
-
[14]
Peking Math
Flynn, P., Ouyang, Z., Pausader, B., Widmayer, K., Scattering map for the Vlasov-Poisson system. Peking Math. J. , 6, 365–392, (2023)
2023
-
[15]
Gagnebin, A., Iacobelli, M., Landau damping on the torus for the Vlasov-Poisson system with massless electrons. J. Diff. Equations , 376, 154–203, (2023)
2023
-
[16]
Transport Theory Statist
Glassey, R., Schaeffer, J., Time decay for solutions to the linearized Vlasov equation. Transport Theory Statist. Phys. 23, 411–453 (1994)
1994
-
[17]
T., Rodnianski, I., Landau damping for analytic and Gevrey data
Grenier, E., Nguyen, T. T., Rodnianski, I., Landau damping for analytic and Gevrey data. Math. Res. Lett. 28, 1679–1702 (2021)
2021
-
[18]
Guo, Y., Lin, Z., The existence of stable BGK waves. Comm. Math. Phys. 352, 1121–1152 (2017)
2017
-
[19]
Guo, Y., Rein, G., A non-variational approach to nonlinear stability in stellar dynamics applied to the King model. Comm. Math. Phys. 271, 489–509 (2007)
2007
-
[20]
Hadˇ zi´ c, M., Moreno, M., On absence of embedded eigenvalues and stability of BGK waves.Preprint, arXiv:2412.07025, (2024)
2024 arXiv
-
[21]
Hadˇ zi´ c, M., Rein, G., Straub, C., On the Existence of Linearly Oscillating Galaxies. Arch. Rational Mech. Anal. 243, 611–696 (2022)
2022
-
[22]
Accepted to Arch
Hadˇ zi´ c, M., Rein, G., Schrecker, M., Straub, C., Damping versus oscillations for a gravitational Vlasov- Poisson system. Accepted to Arch. Rational Mech. Anal. , arXiv:2301.07662, (2023)
2023 arXiv
-
[23]
Preprint, arXiv:2405.17153 (2024)
Hadˇ zi´ c, M., Rein, G., Schrecker, M., Straub, C., Quantitative Gravitational Relaxation: pure transport case. Preprint, arXiv:2405.17153 (2024)
2024 arXiv
-
[24]
T., Rousset, F., Asymptotic stability of equilibria for screened Vlasov- Poisson systems via pointwise dispersive estimates
Han-Kwan, D., Nguyen, T. T., Rousset, F., Asymptotic stability of equilibria for screened Vlasov- Poisson systems via pointwise dispersive estimates. Ann. PDE 7, Art. 18, (2021)
2021
-
[25]
T., Rousset, F., On the Linearized Vlasov-Poisson System on the Whole Space Around Stable Homogeneous Equilibria
Han-Kwan, D., Nguyen, T. T., Rousset, F., On the Linearized Vlasov-Poisson System on the Whole Space Around Stable Homogeneous Equilibria. Comm. Math. Phys. 387, 1405–1440 (2021). 68 MAHIR HAD ˇZI ´C AND MATTHEW SCHRECKER
2021
-
[26]
Ionescu, A., Pausader, B., Wang, X., Widmayer, K., Nonlinear Landau damping for the Vlasov-Poisson system in R3: the Poisson equilibrium. Ann. PDE 10, Art. 2, (2024)
2024
-
[27]
Jia, H., Linear inviscid damping near monotone shear flows, SIAM Journal on Mathematical Analysis 52 (2020), 623–652
2020
-
[28]
J., Dynamics of flat galaxies
Kalnajs, A. J., Dynamics of flat galaxies. IV. The integral equation for normal modes in matrix form. Astrophys. J. 212, 637–644 (1977)
1977
-
[29]
Landau, L., On the vibrations of the electronic plasma. Akad. Nauk SSSR. Zhurnal Eksper. Teoret. Fiz. 16, 574–586 (1946)
1946
-
[30]
Kunze, M., A Birman-Schwinger Principle in Galactic Dynamics , Progress in Mathematical Physics 77, Birkh¨ auser Basel 2021
2021
-
[31]
Lemou, M., M´ ehats, F., Rapha¨ el, P., Orbital stability of spherical galactic models.Invent. Math. 187, 145–194 (2012)
2012
-
[32]
Lions, P.-L., Perthame, B., Propagation of moments and regularity for the 3-dimensional Vlasov–Poisson system. Invent. Math. 105, 415–430 (1991)
1991
-
[33]
Lynden-Bell, D., The stability and vibrations of a gas of stars. Mon. Not. R. Astr. Soc. 124, 279–296 (1962)
1962
-
[34]
Monthly Notices of the Royal Astronomical Society, 136, p.101 (1967)
Lynden-Bell, D., Statistical mechanics of violent relaxation in stellar systems. Monthly Notices of the Royal Astronomical Society, 136, p.101 (1967)
1967
-
[35]
Differential Equations 179, 625–646 (2002)
Ma˜ nosas, F., Villadelprat, J., Area-Preserving Normalizations for Centers of Planar Hamiltonian Sys- tems, J. Differential Equations 179, 625–646 (2002)
2002
-
[36]
Moreno, M., Rioseco, P., Van Den Bosch, H., Mixing in anharmonic potential well. J. Math. Phys. 63, 071502 (2022)
2022
-
[37]
Acta Math
Mouhot, C., Villani, C., On Landau damping. Acta Math. 207, 29–201 (2011)
2011
-
[38]
Preprint, arXiv:2412.18620, (2024)
Nguyen, T., Landau damping below survival threshold. Preprint, arXiv:2412.18620, (2024)
2024 arXiv
-
[39]
Pausader, B., Widmayer, K., Stability of a point charge for the Vlasov-Poisson system: the radial case. Commun. Math. Phys. , 385, 1741–1769, (2021)
2021
-
[40]
Pausader, B., Widmayer, K., Yang, J., Stability of a point charge for the repulsive Vlasov-Poisson system J. Eur. Math. Soc. , published online, (2024)
2024
-
[41]
Pfaffelmoser, K., Global classical solutions of the Vlasov–Poisson system in three dimensions for general initial data. J. Differential Equations , 95, 281–303 (1992)
1992
-
[42]
Classical Quantum Gravity 37, 195027 (2020)
Rioseco, P., Sarbach, O., Phase space mixing in an external gravitational central potential. Classical Quantum Gravity 37, 195027 (2020)
2020
-
[43]
Smulevici, J., Small data solutions of the Vlasov-Poisson system and the vector field method. Ann. PDE 2, Art. 11, (2016)
2016
-
[45]
Preprint, arXiv:2501.04175
Weder, R., Galaxy dynamics, gravitational Vlasov-Poisson system, Landau damping, and scattering theory. Preprint, arXiv:2501.04175. Department of Mathematics, University College London, UK Email address : m.hadzic@ucl.ac.uk Department of Mathematics, University of Bath, UK Ema...
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.