REVIEW 3 major objections 4 minor 40 references
Arbitrary Polynomial Decay Rates of Neutral, Collisionless Plasmas
T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read A neutral plasma can be made to decay at any polynomial rate.
desk verdict The inductive decay machinery is the real contribution; the existence step in Theorem 1.3 leans on an unverified scattering-map application that is probably repairable but must be fixed before the result stands. 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
The central object is the translated distribution $g^\alpha(t,x,v)=f^\alpha(t,x+vt,v)$ and its spatial-moment coefficients $F^{\alpha,\ell}(t,v)=\sum_{|\beta|=\ell}\frac{1}{\beta!}\int(-y)^\beta D_v^\beta g^\alpha(t,y,v)\,dy$. A Taylor expansion of $g^\alpha$ in the velocity argument turns the charge density into a sum of $t^{-\ell}\rho_\ell(t,x)$, whose limits $\rho_{\ell,\infty}$ are the asymptotic profiles. The field-decay assumption (A) keeps the spatial support of $g^\alpha$ bounded, which controls the Taylor remainder and the derivative norms $G^k_v$ and $G^k_{x,v}$; with those bounds, an induction over $\ell$ shows that each vanishing moment improves the decay by one power of $t$ and determines the next profile.
What would settle it
Numerically integrate the Vlasov-Poisson system for a two-species neutral plasma whose limiting spatial distributions are Gaussians with equal means but different variances, as described in Remark 1.4; the paper predicts $\|E(t)\|_{\infty}\sim t^{-3}$, so observing the dispersive $t^{-2}$ rate would refute Theorem 1.3. A separate check is whether the compactly supported limiting states built with Gegenbauer polynomials satisfy the smallness hypothesis of the cited scattering-map theorem; if they do not, the existence step is unsupported.
Extended reading notes
Core claim
The central claim is Theorem 1.3: for every $m\in\mathbb{N}_0$ there exist nontrivial solutions in $C^{m+1}((0,\infty)\times\mathbb{R}^6)$ of the multispecies Vlasov-Poisson system with zero total net charge such that $\|\rho(t)\|_{\infty}\sim t^{-m-3}$ and $\|E(t)\|_{\infty}\sim t^{-m-2}$; for $m\ge 1$ all derivatives $\nabla_x^k E$ up to order $m$ decay like $t^{-m-3}$, and $f^\alpha(t,x+vt,v)\to f^\alpha_\infty(x,v)$ in $L^\infty$ at rate $t^{-m}$. The mechanism is charge cancellation in the limiting spatial averages: the functions $F^{\alpha,\ell}_\infty$ and the associated densities $\rho_{\ell,\infty}$ form the coefficients of a polyhomogeneous expansion of $\rho$ and $E$, and each vanishing coefficient buys one extra power of time decay. Theorem 1.4 propagates this by induction: assuming $\rho_{\ell,\infty}\equiv 0$ for $\ell=0,\ldots,n$ yields the next-order convergence at rate $t^{-1}$ and uniform bounds on all relevant derivatives of the translated distribution. The paper states this is the first construction of nontrivial decay faster than the dispersive rates for the Vlasov-Poisson system in $\mathbb{R}^3$.
Load-bearing premise
The argument stands on the electric field decaying at least as fast as $t^{-5/3}$ and on the scattering-map theorem supplying a true solution from each constructed asymptotic state; if either fails, the existence claim in Theorem 1.3 collapses.
Editorial extensions
If this is right
- For every $m\ge 0$ there are nontrivial neutral solutions whose charge density and electric field decay exactly like $t^{-m-3}$ and $t^{-m-2}$, respectively, so no single universal polynomial rate governs neutral plasmas.
- For $m\ge 1$ the particle distribution scatters linearly at rate $t^{-m}$, while the $m=0$ case scatters in the modified sense with a logarithmic correction; these are the only two scattering behaviors that occur.
- If $\rho_{0,\infty},\ldots,\rho_{n,\infty}$ all vanish, then $\|\rho(t)\|_{\infty}\lesssim t^{-n-4}$ and $\|E(t)\|_{\infty}\lesssim t^{-n-3}$, with the next profile $\rho_{n+1,\infty}$ determining the leading term.
- The same induction yields bounds on all derivatives of the electric field and of the translated distribution, so the asymptotic profiles of $E$, $\nabla E$, $\rho$, and their derivatives are all identified.
- The construction uses compactly supported data and, as the paper notes, the method extends to higher dimensions and to the small-data setting where assumption (A) is known to hold.
Reading between the lines
- The mechanism suggests a practical diagnostic: measuring the first nonvanishing limiting moment of the charge density of a neutral plasma would directly predict the long-time decay exponent of its field and density.
- The two-species Gaussian example in the paper could be simulated numerically with compactly supported approximants; matching the predicted $t^{-3}$ field decay would confirm that the cancellation mechanism does not depend on the special Gegenbauer construction used in the proof.
- If assumption (A) holds with $p$ only slightly above $5/3$, the support-growth estimates become marginal; whether the same decay hierarchy survives weaker field decay is left open by the paper and is a natural stress test.
- An analogous hierarchy of polynomial rates may hold for relativistic Vlasov-Maxwell or screened Coulomb systems whenever a scattering map with suitable asymptotic states exists, as the authors suggest in a remark; that extension is not proved here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the multispecies Vlasov-Poisson system in R^3 under a neutrality assumption M=0 and an electric-field decay assumption (A). The main analytic result, Theorem 1.4, shows that if the first n+1 asymptotic moment densities ρ_{ℓ,∞} vanish, then the charge density, electric field, and their derivatives admit sharp asymptotic profiles with rates t^{-n-4}, t^{-n-3}, etc., and the distribution functions scatter linearly. The companion result, Theorem 1.3, claims that for every m in N0 there exist nontrivial solutions realizing the decay rates ‖ρ(t)‖∞ ∼ t^{-m-3} and ‖E(t)‖∞ ∼ t^{-m-2}. Sections 2 and 4 develop a lengthy induction based on Taylor expansion of the translated distribution functions and on estimates for G^k_v and G^k_{x,v}; the proof of Theorem 1.4 is carried out in Section 3. The existence part of Theorem 1.3 is delegated to an external scattering map theorem [14] at the end of Section 3, and this delegation is not verified against the hypotheses used elsewhere in the paper.
Significance. If the claimed results hold, they provide, for the first time, a countably infinite family of sharp polynomial decay rates for neutral plasmas in R^3, interpolating between the standard dispersive rates and the much faster rates associated with phase mixing. The internal derivation is genuinely parameter-free: the decay rate is forced by the first nonvanishing moment of the asymptotic state, and no fitted parameters appear. The inductive scheme in Sections 2–4 is a substantial and coherent piece of analysis, and the paper is careful to identify which estimates are preliminary and which are sharp. The main weakness is the existence step for Theorem 1.3, which relies on an unverified application of the inverse scattering map of [14]; because the paper's own Theorem 1.4 requires compactly supported initial data and moment conditions on the actual evolved solution, this gap is load-bearing. The result is likely repairable, but the manuscript as written does not close the existence argument.
major comments (3)
- [Section 3, final paragraph (proof of Theorem 1.3)] The existence assertion for every m in N0 rests on the sentence invoking [14, Theorem 1.1(ii) and Remark 1.2(5)], but the hypotheses of that scattering map are not stated and are not verified for the constructed states f^α_{m,∞}=φ^α_m(x)ψ_m(v). In particular, the paper does not check the required smallness condition of [14] for these compactly supported profiles, and it does not check that the initial data produced by the scattering map are compactly supported. This matters for two reasons: Theorem 1.4 is stated only for initial data in C^{n+2}_c(R^6), and its proof repeatedly uses compact support (Lemma 2.1 and the uniform spatial support bound (17)); moreover, compact support of the asymptotic state f_∞ does not imply compact support of the initial data obtained from an inverse scattering construction. The existence claim of Theorem 1.3 is therefore not proved as written. The gap is likely repairable either by stating and verifying the hypotheses of [14] (for instance, by a small-amplitude rescaling if the relevant norm is homogeneous) or by an independent constructive argument, but the present text does not supply such a verification.
- [Section 3, final paragraph and Theorem 1.4] The application of the scattering map must also ensure that the conditions ρ_{ℓ,∞}≡0 for ℓ≤n hold for the actual evolved solution, not merely for the prescribed profile f^α_{m,∞}. Theorem 1.4 assumes these conditions for the solution and later identifies F^{α,ℓ}_∞ with moments of f^α_∞ through equation (11). If [14] supplies only L∞ scattering of f^α, or scattering in a norm that does not control the derivatives appearing in the definitions (9) and (11), then the moment limits of the evolved solution might differ from the moments of the constructed f^α_{m,∞}. The compact support of the prescribed profile does not by itself guarantee the needed vanishing of the evolved solution's moments. The proof should either quote a derivative-scattering statement from [14]/[38] that covers these moments or justify the moment limits directly.
- [Section 3, proof of Theorem 1.3, m=0 case] For m=0 the construction chooses nonnegative ψ^α_0 satisfying Σ_α q_α ψ^α_0(v)=η(v), where η is an arbitrary C^1_c function with zero integral. Such a representation requires an explicit feasibility argument when the charges q_α have mixed signs, since the right-hand side must be decomposable into nonnegative compactly supported pieces with the prescribed charges. A short argument, for example choosing large common nonnegative profiles and then adding small signed corrections, would make the construction complete; as written the constraint is asserted without proof.
minor comments (4)
- [Section 3, proof of Theorem 1.3] In the definition of Φ_m(x), the expression '/BD_{[-1,1]}' appears to be a typographical error for the indicator function 1_{[-1,1]}; please correct the notation.
- [Reference list] Reference [22] contains the typo 'Arsigmave for rational mechanics and analysis'; this should read 'Archive for Rational Mechanics and Analysis'.
- [Remark 1.4] In the first bullet, the informal phrase 'the distributions don't overlap much at all' is imprecise; the vanishing of ρ_{0,∞} in that case follows from the normalization ∫(f^+ - f^-) dx = 0 regardless of the means, so the wording could be clarified.
- [Theorem 1.3, m=0 case] The statement uses the notation f^α_∞ ∈ C^m_c(R^6) with m=0; this is acceptable but slightly nonstandard, and it may be clearer to write C^0_c or to state the regularity separately for the m=0 case.
Circularity Check
No circularity found: decay rates are consequences of prescribed asymptotic moments and an independent induction; the unverified [14] application is a correctness gap, not a circular reduction.
full rationale
The central derivation is not circular: Theorem 1.3's rates are not fitted parameters but are forced by the first nonvanishing moment of a prescribed asymptotic state. The paper constructs f^α_{m,∞} with ∫ x^β μ_m(x) dx = 0 for |β| < m and nonzero for β = (m,0,0), so ρ_{ℓ,∞} ≡ 0 for ℓ < m and ρ_{m,∞} ≠ 0. Theorem 1.4 then proves, by induction and Taylor expansion, that any solution satisfying these moment conditions has ρ ~ t^{-m-3} and E ~ t^{-m-2}. No equation in the proof assumes the target rate; assumption (A) is an input hypothesis, not a fitted value. The use of [31] is a self-citation, but [31] is an independent prior theorem proving the base modified-scattering estimates and does not rely on the present conclusions. The scattering-map invocation of [14] and [38] is external. A genuine gap exists: the paper does not verify [14]'s smallness hypothesis for the constructed states, nor that the resulting Cauchy data are compactly supported, both of which are needed for Theorem 1.4; however, this is a hypothesis-checking omission, not a definitional equivalence or a fitted input renamed as a prediction, so it does not constitute circularity under the hard rules.
Assumptions & free parameters
assumptions (4)
- domain assumption The electric field decays at a rate ||E(t)||_infty <= C t^{-p} for some p>5/3 (assumption (A)).
- domain assumption The scattering map theorem of [14] applies to the constructed asymptotic states f^alpha_{m,infty}.
- domain assumption Global existence of smooth solutions for compactly supported initial data.
- domain assumption Neutrality M=0 and, in Theorem 1.4, the vanishing of moments rho_{ell,infty} identically 0 for ell=0..n.
Cite this review
Pith. "Pith review of Arbitrary Polynomial Decay Rates of Neutral, Collisionless Plasmas." pith.science (2026). https://pith.science/paper/5BSNY3MG
@misc{pith2026250202678,
author = {Pith},
title = {Pith review of: Arbitrary Polynomial Decay Rates of Neutral, Collisionless Plasmas},
year = {2026},
howpublished = {\url{https://pith.science/paper/5BSNY3MG}},
note = {Machine review of arXiv:2502.02678}
}
abstract
A multispecies, collisionless plasma is modeled by the Vlasov-Poisson system. Assuming the plasma is neutral and the electric field decays with sufficient rapidity as $t \to\infty$, we show that solutions can be constructed with arbitrarily fast, polynomial rates of decay, depending upon the properties of the limiting spatial average and its derivatives. In doing so, we establish, for the first time, a countably infinite number of asymptotic profiles for the charge density, electric field, and their derivatives, each of which is necessarily realized by a sufficiently smooth solution and exceeds the established dispersive decay rates. Finally, in each case we establish a linear $L^\infty$ scattering result for every particle distribution function, namely we show that they converge as $t \to \infty$ along the transported spatial characteristics at increasingly faster rates.
Reference graph
Works this paper leans on
-
[31]
Pankavich, S., Asymptotic Dynamics of Dispersive, Col lisionless Plasmas. Comm. Math. Phys. 2022 391: 455-493
work page 2022
-
[14]
et al., Scattering M ap for the Vlasov-Poisson System, Peking Math J
Flynn, P., Ouyang, Z., Pausader, B. et al., Scattering M ap for the Vlasov-Poisson System, Peking Math J. 2023 6: 365-392
work page 2023
-
[38]
Schlue, V. and Taylor, M., Inverse modified scattering a nd polyhomogeneous expansions for the Vlasov–Poisson syst em. arXiv:2404.15885 2024
arXiv 2024
-
[1]
and Stegun, I., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables
Abramowitz, M. and Stegun, I., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Dover (9th Ed) 1964
work page 1964
-
[2]
Bardos, C. and Degond, P., Global existence for the Vlaso v-Poisson equation in 3 space variables with small initial d ata. Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire1985, 2(2): 101-118
-
[3]
Advances in Differential Equations 1998, 3: 271-292
Batt, J., Kunze, M., and Rein, G., On the asymptotic behav ior of a one-dimensional, monocharged plasma and a rescalin g method. Advances in Differential Equations 1998, 3: 271-292
work page 1998
-
[4]
Communications on Pure and Applied Mathemat ics 2017, 71(3): 537-576
Bedrossian, J., Masmoudi, N., and Mouhot, C., Landau Dam ping in Finite Regularity for Unconfined Systems with Screen ed Interactions. Communications on Pure and Applied Mathemat ics 2017, 71(3): 537-576
work page 2017
-
[5]
Ben-Artzi, J., Calogero, S., and Pankavich, S., Arbitra rily large solutions of the Vlasov-Poisson system. SIAM J. M ath. Anal. 2018, 50(4): 4311-4326
work page 2018
Show all 40 references
-
[6]
Kinetic and Related Models 2024, 17(1): 29-51
Ben-Artzi, J., Morisse, B., and Pankavich, S., Asymptot ic Growth and Decay of Two-dimensional Symmetric Plasmas. Kinetic and Related Models 2024, 17(1): 29-51
2024
-
[7]
and Pankavich, S., Modified Scattering of S olutions to the Relativistic Vlasov-Maxwell System Inside the Light Cone, arXiv:2306.11725
Ben-Artzi, J. and Pankavich, S., Modified Scattering of S olutions to the Relativistic Vlasov-Maxwell System Inside the Light Cone, arXiv:2306.11725
-
[8]
arXiv:2208.08360
Bigorgne, L., Global existence and modified scattering f or the solutions to the Vlasov-Maxwell system with a small di stri- bution function. arXiv:2208.08360
-
[9]
arXiv:2312 .12214
Bigorgne, L., Scattering map for the Vlasov-Maxwell sys tem around source-free electromagnetic fields. arXiv:2312 .12214
-
[10]
and Ruiz R
Bigorgne, L. and Ruiz R. V., Late-time asymptotics of sm all data solutions for the Vlasov-Poisson system. arXiv:24 04.05812 2024
2024
-
[11]
Castella F., Propagation of space moments in the Vlasov -Poisson equation and further results. Ann. de l’Inst. Henr i Poincare-Anal. non lineaire 1999, 16: 503-533
1999
-
[12]
and Zhang, X., Global Existence to the Vlasov-P oisson System and Propagation of Moments Without Assumptio n of Finite Kinetic Energy
Chen, Z. and Zhang, X., Global Existence to the Vlasov-P oisson System and Propagation of Moments Without Assumptio n of Finite Kinetic Energy. Communications in Mathematical P hysics 2016, 343: 851-879
2016
-
[13]
and Kwon, S., Modified scattering for the Vlasov -Poisson system, Nonlinearity 2016 29: 2755
Choi, S. and Kwon, S., Modified scattering for the Vlasov -Poisson system, Nonlinearity 2016 29: 2755
2016
-
[15]
SIAM : 1996
Glassey, R., The Cauchy Problem in Kinetic Theory. SIAM : 1996
1996
-
[16]
Glassey, R., Pankavich, S., and Schaeffer, J., Decay in T ime for a One-Dimensional, Two Component Plasma. Math. Meth Appl. Sci. 2008, 31:2115-2132
2008
-
[17]
Kinetic and Related Models 2009, 2: 465-488
Glassey, R., Pankavich, S., and Schaeffer, J., On long-t ime behavior of monocharged and neutral plasma in one and one -half dimensions. Kinetic and Related Models 2009, 2: 465-488
2009
-
[18]
Quarterly of Applied Mathematics 2010, 68: 135-141
Glassey, R., Pankavich, S., and Schaeffer, J., Time Deca y for Solutions to the One-dimensional Equations of Plasma Dynamics. Quarterly of Applied Mathematics 2010, 68: 135-141
2010
-
[19]
Differential & Integral Equations 2010 23: 61-77
Glassey, R., Pankavich, S., and Schaeffer, J., Large Tim e Behavior of the Relativistic Vlasov-Maxwell System in Low Space Dimension. Differential & Integral Equations 2010 23: 61-77
2010
-
[20]
Hormander, L., The Analysis of Linear Partial Different ial Operators I: Distribution Theory and Fourier Analysis, 2nd ed., Springer-Verlag, Berlin, 1990
1990
-
[21]
Horst, E., Symmetric plasmas and their decay. Comm. Mat h. Phys. 1990, 126: 613-633
1990
-
[22]
Arsigmave for rational mech anics and analysis 2011, 200: 313-360
Hwang, H., Rendall, A., and Velazquez, J., Optimal grad ient estimates and asymptotic behaviour for the Vlasov-Poi sson system with small initial data. Arsigmave for rational mech anics and analysis 2011, 200: 313-360
2011
-
[23]
International Mathematics Research Notices 2022, 2022: 8865-8889
Ionescu, A., Pausader, B., W ang, X., Widmayer, K., On th e asymptotic behavior of solutions to the Vlasov-Poisson sy stem. International Mathematics Research Notices 2022, 2022: 8865-8889
2022
-
[24]
Annals of PDE 2024, 10(1), Article 2
Ionescu, A., Pausader, B., W ang, X., Widmayer, K., Nonl inear Landau Damping for the Vlasov-Poisson System in R3: The Poisson Equilibrium. Annals of PDE 2024, 10(1), Article 2
2024
-
[25]
D., Classical Electrodynamics
Jackson, J. D., Classical Electrodynamics. Wiley: 1999
1999
-
[26]
and Perthame, B., Propogation of moments an d regularity for the three dimensional Vlasov-Poisson syst em
Lions, P.L. and Perthame, B., Propogation of moments an d regularity for the three dimensional Vlasov-Poisson syst em. Invent. Math. 1991, 105: 415-430
1991
-
[27]
and Villani, C., On Landau damping
Mouhot, C. and Villani, C., On Landau damping. Acta Math 2011, 207, 29-201. 34
2011
-
[28]
Pallard, C., Space Moments of the Vlasov–Poisson Syste m: Propagation and Regularity. SIAM J. Math. Anal. 2014 46(3): 1754-1770
2014
-
[29]
Kinetic & Related Models 2011, 4(2): 549-567
Pallard, C., Growth estimates and uniform decay for a co llisionless plasma. Kinetic & Related Models 2011, 4(2): 549-567
2011
-
[30]
Pankavich, S., Exact Large Time Behavior of Sphericall y-Symmetric Plasmas. SIAM J. Math. Anal. 2021 53(4): 4474-4512
2021
-
[32]
Pankavich, S., Scattering and Asymptotic Behavior of S olutions to the Vlasov-Poisson system in High Dimension. SI AM J. Math. Anal. 2023 55(5): 4727-4750
2023
-
[33]
Pfaffelmoser, K., Global classical solution of the Vlas ov-Poisson system in three dimensions for general initial d ata. J. Diff. Eq. 1992, 95(2): 281-303
1992
-
[34]
Rein G., Collisionless Kinetic Equations from Astroph ysics - The Vlasov-Poisson System, in Handbook of Differenti al Equations, Evolutionary Equations, 2007, Eds. C. M. Dafermos and E. Feireisl, Elsevier: 383-479
2007
-
[35]
PDE 1991, 16(8-9): 1313-1335
Schaeffer, J., Global existence of smooth solutions to t he Vlasov-Poisson system in three dimensions, Comm. PDE 1991, 16(8-9): 1313-1335
1991
-
[36]
Schaeffer, J., Asymptotic growth bounds for the Vlasov- Poisson system. Math. Meth. Appl. Sci. 2011, 34: 262-277
2011
-
[37]
Schaeffer, J., Large-time behavior of a one-dimensiona l monocharged plasma. Diff. and Int. Equations 2007, 20(3): 277-292
2007
-
[39]
Smulevici, J., Small data solutions of the Vlasov-Pois son system and the vector field method. Ann. PDE 2016 2(2)
2016
-
[40]
Methods Appl
Yang, D., Growth estimates and uniform decay for the Vla sov-Poisson system, Math. Methods Appl. Sci. 2017, 40: 4906-4916. Department of Applied Mathematics and Statistics, Colorado School of Mines, Golden, CO 80401. Email address : gmattingly@mines.edu Email address : pankavi...
2017
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.