Nonlinear stability of the one dimensional screened Vlasov Poisson equation

Dongyi Wei

REVIEW 1 major objections 4 minor 20 references

Nonlinear stability of the one dimensional screened Vlasov Poisson equation

T0 review · 1 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper proves that small Gevrey-2 initial data for the one-dimensional screened Vlasov-Poisson (Vlasov-Yukawa) equation lead to global solutions whose density derivatives decay at exactly the free-streaming rate, $(t+1)^{-n-1}$.

desk verdict Proves the open 1D Vlasov-Yukawa stability problem with Gevrey-2 data; technically strong, but the existence passage needs a spelled-out compactness argument. read the letter →

arxiv 2411.13798 v1 pith:AUZZFWFP submitted 2024-11-21 math.AP

Dongyi Wei This is my paper · ORCID

classification math.AP MSC 35Q8335B40

keywords screenedVlasov-PoissonVlasov-YukawaGevreyregularitynonlinearstabilitydecayestimatessmalldataonedimensionFaàdiBrunoformula

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper proves nonlinear asymptotic stability near vacuum for the one-dimensional screened Vlasov-Poisson (Vlasov-Yukawa) system, for both attractive and repulsive forces. For smooth initial data that are small in a Gevrey-2 sense, the solution exists globally in time and every spatial derivative of the density decays like $(t+1)^{-n-1}$, exactly the rate that free streaming would produce. This settles the one-dimensional case of a stability problem that had remained open, with only analytic small-data long-time stability known previously. The proof is an iterative construction in which the key step is a uniform estimate for the density built from Faà di Bruno expansions and carefully chosen Gevrey weight functions.

What carries the argument

The argument is carried by an iteration scheme that starts with zero density and solves the linear Vlasov equation at each step, combined with an auxiliary boundary value problem for characteristics $X(s;x,x_0,t)$ defined by $\frac{d}{ds}X = V$, $\frac{d}{ds}V = -q(\partial_x\varphi)(X,s)$, $X(t)=x$, $X(0)-V(0)=x_0$. This yields the representation $\rho^*(x,t)=\int_{\mathbb{R}} \tilde f_0(x_0, w_0(x,x_0,t))\,\partial_{x_0} w(x,x_0,t)\,dx_0$. The core technical engine is a set of combinatorial inequalities, principally (2.5) in Lemma 2.1, that bound Faà di Bruno sums by the Gevrey weights $\varphi_n(t)=e^{(n-2)\sqrt{t}/(n+\sqrt{t})}$ with explicit constants, together with the auxiliary function $\gamma(t)=0.01\ln(t+1)+0.99t+1$, which controls the characteristic estimates through comparison arguments. These ingredients let every derivative estimate close in a bootstrap, yielding the uniform decay for the density.

What would settle it

A concrete check would be to evaluate the sum on the left of (2.5) numerically for moderate $n$ (say $n=3,4,5$) and $t \in [0,10]$ and see whether it ever exceeds $n!\varphi_n(t)e^{0.15}$; if it does, the stated constants cannot close the bootstrap and Theorem 1.1 would have to be weakened. Alternatively, a high-resolution numerical solution of the 1D screened Vlasov-Poisson equation with initial data saturating (1.2) could test directly whether $(t+1)^{n+1}|\partial_x^n\rho(x,t)|$ stays uniformly bounded.

Watch

Extended reading notes

Core claim

Theorem 1.1 states: if the initial data $f_0$ is smooth, $f_0 \in W^{4,1}(\mathbb{R}\times\mathbb{R})$, and $\|(\partial_x+\partial_v)^{n+1} f_0\|_{L^1} \le (n!)^2/10^4$ for all $n \ge 0$, then the solution is global in time and the density satisfies $|\partial_x^n \rho(x,t)| \le 3^n (n!)^2 (t+1)^{-n-1}/10^3$ for all $n \ge 0$, $t \ge 0$. In other words, under Gevrey-2 smallness of the transformed initial datum, the nonlinear plasma behaves asymptotically like the free-streaming equation: the screening potential does not alter the decay rate of any density derivative. The paper also notes that by time reversibility the same statement holds for $t \to -\infty$ when the condition is imposed with $\partial_x - \partial_v$ instead of $\partial_x + \partial_v$, and that modified scattering can be proven along the lines of [15].

Load-bearing premise

The proof rests on the Gevrey-2 smallness condition (1.2) and, at its technical core, on the combinatorial inequality (2.5) in Lemma 2.1, which must hold for all $n$ and $t$ with the stated constants; every principal estimate in Lemmas 3.4, 3.5, and 3.7 depends on it, so a single failure of that inequality would invalidate the main theorem.

Editorial extensions

If this is right

Global existence holds for all $q=\pm 1$, attractive or repulsive, under the Gevrey-2 smallness condition.
Each density derivative decays at the free-streaming rate: $|\partial_x^n \rho(x,t)| \le 3^n (n!)^2 (t+1)^{-n-1}/10^3$ for every $n\ge 0$.
Time reversibility gives the same decay as $t\to -\infty$ when the smallness is imposed with $\partial_x-\partial_v$.
The same method can be pushed to prove modified scattering, analogous to the result cited in [15].
The iteration sequence converges to a $C^1$ solution, so the decay estimates hold for the actual nonlinear solution, not merely for approximate densities.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A natural next test is whether the Gevrey-2 smallness can be relaxed to a finite-regularity condition with all the derivative bounds replaced by a weighted Sobolev norm; the present proof needs all orders, so the rate may degrade if only finitely many derivatives are controlled.
The characteristic representation used here could plausibly extend to other screened kinetic models, such as Vlasov equations with Yukawa kernels on bounded or curved domains, where the same comparison-function method may apply.
If the theorem is correct, the decay rate $(t+1)^{-n-1}$ is optimal in the sense that it equals the linear free-streaming rate, so nonlinear screening causes no slowdown of density dispersion in one dimension; a numerical check of the uniform boundedness of $(t+1)^{n+1}|\partial_x^n\rho|$ would be a direct test.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

This paper proves a genuinely new result: asymptotic stability for the 1D screened Vlasov-Poisson (Vlasov-Yukawa) system near vacuum with Gevrey-2 small data. The decay rate (t+1)^{-n-1} for ∂_x^n ρ matches free streaming exactly, and this closes the gap in the literature, where d≥2 or analytic data in 1D was the previous state of the art. The core of the paper is a well-executed bootstrap on densities and characteristics, using carefully chosen Gevrey weights and Faà di Bruno estimates. I found the main lemmas (2.1, 3.2, 3.5, 3.7) detailed and internally consistent, and the bootstrap closes with a comfortable margin: the assumption uses a factor 1500 and the conclusion gives 3000. The smallness assumption in (1.2) is used honestly, with no hidden circularity.

The soft spot is the existence step in the proof of Theorem 1.1. After obtaining uniform bounds on the iterates, the paper says, in one sentence, that (1.8) plus 'the same arguments as in [1,6]' gives a Cauchy sequence converging to a C^1 solution. That is the only load-bearing step that is not actually shown. The uniform Gevrey bounds on ρ(k) and on the characteristic derivatives are strong enough that a standard compactness argument (Arzelà-Ascoli plus passage to the limit in the characteristics) should work, and I do not think the claim is false. But it is an omitted proof, and the referenced papers are finite-regularity settings, not a precise template for this Gevrey framework. I would want a referee to require this to be written out or replaced by a precise reference. It is the difference between a complete proof and a proof sketch at the critical point.

There are also two mechanical typos. The statement of (1.10) gives γ^{-n} while the proved (3.24) gives γ^{-(n+1)}. In Lemma 3.7, the proof cites Lemma 3.4 where it should cite Lemma 3.6, and the first term in the definition of F(s) should be ∂_1^2 φ under ∂_x^n, not ∂_1 φ. Neither affects the estimates.

Overall: worth sending to a serious referee. The result is significant, the technical core looks sound, and the gaps are fixable. After the author expands the convergence argument and corrects the typos, I would expect this to be accepted in a strong analysis journal. I would cite it if I worked in kinetic theory.

Referee Report

1 major / 4 minor

Summary. The paper studies the one-dimensional screened Vlasov-Poisson (Vlasov-Yukawa) system ∂_t f + v∂_x f - q∂_x φ ∂_v f = 0, (1-∂_x^2)φ = ρ, with d=1. It claims that for smooth initial data satisfying the Gevrey-2 smallness condition ‖(∂x+∂v)^{n+1} f0‖_{L^1} ≤ (n!)^2/10^4 for all n ≥ 0, the solution exists globally in time and the density satisfies |∂_x^n ρ(x,t)| ≤ 3^n (n!)^2 (t+1)^{-n-1}/10^3, i.e., the free-streaming decay rate. The proof is based on an iteration scheme, a characteristic representation of the density, and a bootstrap argument with Gevrey weights φ_n(t) and a time weight γ(t), closing with a factor-of-two margin. The main body consists of detailed estimates for Faà di Bruno sums, the characteristic ODEs, and the density representation.

Significance. The result, if fully established, would resolve an open problem stated in the introduction: the nonlinear asymptotic stability of the 1D screened Vlasov-Poisson system near vacuum with Gevrey-2 small data, with the exact decay rate of free streaming. The paper's strongest feature is its explicit, self-contained derivation of the bootstrap estimates: Lemma 2.1 gives a sharp combinatorial bound, Lemmas 3.2–3.7 provide the characteristic and density estimates, and all constants are explicit (1500 to 3000 margin). The verification of the decay rate against the free-streaming solution is a good falsifiable check. The principal weakness is that the existence part of Theorem 1.1 is not proved in the manuscript; the transition from uniform bounds to a C^1 solution is imported from [1,6] without a Cauchy or compactness argument.

major comments (1)

[Sec. 1, Proof of Theorem 1.1] The proof of global existence is not self-contained. After establishing the uniform Gevrey bounds (1.8), the paper states: 'Using (1.8) and the same arguments as in [1,6], we can prove that (f^(k), ρ^(k), φ^(k)) is a Cauchy sequence converging to a C^1 solution'. However, no estimate for the difference of consecutive iterates is given, and the references [1,6] address different settings (finite-regularity VP and 2D Vlasov-Yukawa dispersion). The estimates in Sections 2–3 are pointwise Gevrey bounds on ρ^(k) and on characteristic derivatives; they do not by themselves yield strong compactness in time or a limit that satisfies the nonlinear equation. Since global existence is part of Theorem 1.1, this is a load-bearing gap. The revision should include a complete argument, for example a contraction estimate in a weighted Gevrey norm for (ρ^(k+1)-ρ^(k)), or an explicit equicontinuity/compactness argument with the limit passage in the nonlinear term.

minor comments (4)

[Sec. 1, Eq. (1.10)] The displayed bound has γ(t)^n in the denominator, while the corresponding statement (3.24) in Lemma 3.7 has γ(t)^{n+1}. Since the proof of Proposition 1.1 needs the γ^{n+1} version, the exponent in (1.10) should be corrected.
[Sec. 3, Proof of Lemma 3.7] The text invokes 'Lemma 3.4 (assuming (3.21))' but the estimate needed is Lemma 3.6; please correct the citation. In addition, the displayed definition of F(s) contains ∂_x^n[(∂_1 φ)(X,s)∂_{x0}X]; the subsequent use and equation (3.3) require ∂_x^n[(∂_1^2 φ)(X,s)∂_{x0}X].
[Sec. 1, Remarks after Theorem 1.1] The remark 'Similar to [15], we can prove the modified scattering' is a promise without proof. If it remains in the final version, it should be clearly labeled as a remark/conjecture or supported by an argument, since it is not part of the theorem.
[Sec. 1, Introduction] The reference [12] (Iacobelli–Rossi–Widmayer) is cited but not discussed. Since that paper appears to address the same stability question for the screened Vlasov-Poisson equation, a sentence clarifying the distinction (e.g., regularity class, dimension, or decay rates) would help the reader.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the bootstrap estimates are self-contained and the decay law is verified, not fitted.

full rationale

The paper's central claim is derived by a standard iteration/bootstrap argument. Proposition 1.1 assumes a uniform Gevrey bound (1.7) on the input density rho^(n-1), proves a strictly stronger bound on the output density rho^(n), and this closes by induction starting from rho^(0)=0. The Gevrey weight phi_n and the auxiliary function gamma are chosen to satisfy the inequalities of Lemmas 2.1, 2.4 and 2.5, not to reproduce the target decay; the target power (t+1)^(-n-1) is independently checked against the free-streaming calculation following Theorem 1.1. The Faa di Bruno coefficient estimates and characteristic estimates in Section 3 are proved in the paper with explicit constants. No fitted parameter is renamed as a prediction, and no ansatz is imported through a citation: [1] and [6] are used only as references for a standard Cauchy-sequence compactness argument, not as a substitute for the Gevrey bootstrap, and the author's self-citation [19] is a general pointer to Vlasov-Maxwell iteration and is not load-bearing. The brief reference to the Cauchy-convergence argument of [1,6] is an omitted proof detail, a completeness gap rather than a circular reduction.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The proof introduces no new particles, forces, or dimensions. The weights φ_n(t) and γ(t) are test functions chosen to satisfy the estimates. The main imported ingredients are standard tools (Faà di Bruno, maximum principle) and the iteration-convergence framework from prior work.

free parameters (4)

Gevrey weight φ_n(t) = φ_n(t)=exp((n-2)√t/(n+√t)) for n≥2, φ_0=φ_1=1
Chosen by hand to satisfy the algebraic inequalities (2.3), (2.5), (2.8), and Lemma 2.4. It is not derived from physical considerations.
Time weight γ(t) = γ(t)=0.01 ln(t+1)+0.99t+1
Chosen so γ(0)=γ'(0)=1, γ''<0, and 0.99<γ/(t+1)≤1. The constants 0.01 and 0.99 are ad hoc.
Bootstrap constants = 1500 in (1.7), 3000 in Proposition 1.1
Set by hand to give a factor-of-2 improvement for the bootstrap. The exact values are not forced by the problem.
Smallness constant = 10^-4 in (1.2)
Quantitative smallness threshold for the Gevrey-2 data; larger values would break the bounds in (3.20) and (3.17).

assumptions (4)

standard math Faà di Bruno formula
Used repeatedly (Eq. 2.1) to expand derivatives of composed functions.
standard math Maximum principle for (1-Δ)φ=ρ
Used in Section 3 to bound ‖∂_x^n φ‖ by ‖∂_x^n ρ‖.
domain assumption Cauchy-sequence convergence of the iteration as in [1,6]
The proof of Theorem 1.1 asserts convergence of iterates by reference to [1,6] without re-proving it in the 1D Gevrey-2 setting.
standard math L^1 change-of-variables identity
Used to translate the initial condition (1.2) into the f̃0 frame; the Jacobian is 1.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Nonlinear stability of the one dimensional screened Vlasov Poisson equation." pith.science (2026). https://pith.science/paper/AUZZFWFP

@misc{pith2026241113798,
  author       = {Pith},
  title        = {Pith review of: Nonlinear stability of the one dimensional screened Vlasov Poisson equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AUZZFWFP}},
  note         = {Machine review of arXiv:2411.13798}
}

abstract

We study the asymptotic behavior of small data solutions to the screened Vlasov Poisson(i.e. Vlasov-Yukawa) equation on ${\mathbb R}\times{\mathbb R}$ near vacuum. We show that for initial data small in Gevrey-2 regularity, the derivative of the density of order $n$ decays like $(t+1)^{-n-1}$.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 17 canonical work pages

[1]

Bardos, P

C. Bardos, P. Degond, Global existence for the Vlasov–Poisson equation in 3 space variables with small initial data , Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire 2(2):101–118, 1985

work page 1985
[2]

Bedrossian, A brief introduction to the mathematics of Landau damping , arXiv:2211.13707

J. Bedrossian, A brief introduction to the mathematics of Landau damping , arXiv:2211.13707

arXiv
[3]

Bedrossian, N

J. Bedrossian, N. Masmoudi, C. Mouhot, Landau damping: paraproducts and Gevrey regularity , Ann. PDE, 2, pp. Art. 4, 71, 2016

work page 2016
[4]

Bedrossian, N

J. Bedrossian, N. Masmoudi, C. Mouhot, Landau damping in ﬁnite regularity for unconﬁned systems wi th screened interactions, Comm. Pure Appl. Math., 71, 537–576, 2018

work page 2018
[5]

Bedrossian, N

J. Bedrossian, N. Masmoudi, Inviscid damping and the asymptotic stability of planar she ar ﬂows in the 2D Euler equations , Publ. Math. Inst. Hautes ´Etudes Sci., 122, 195-300, 2015

work page 2015
[6]

Choi, S.-Y

S.-H. Choi, S.-Y. Ha, H. Lee, Dispersion estimates for the two-dimensional Vlasov-Yuka wa system with small data , J. Diﬀerential Equations 250(1):515–550, 2011

work page 2011
[7]

Duan, Sharp Decay Estimates for the Vlasov-Poisson and Vlasov-Yu kawa Systems with Small Data , Kinet

X. Duan, Sharp Decay Estimates for the Vlasov-Poisson and Vlasov-Yu kawa Systems with Small Data , Kinet. Relat. Models 15(1):119–146, 2022

work page 2022
[8]

Glassey, W

R. Glassey, W. Strauss. Singularity formation in a collisionless plasma could occu r only at high velocities . Arch. Rational Mech. Anal., 92(1):59–90, 1986

work page 1986

Show all 20 references

[9]

Glassey, W

R. Glassey, W. Strauss. High velocity particles in a collisionless plasma . Math. Methods Appl. Sci., 9(1):46–52, 1987

1987
[10]

Glassey, W

R. Glassey, W. Strauss. Large velocities in the relativistic Vlasov-Maxwell equat ions. J. Fac. Sci. Univ. Tokyo Sect. IA Math., 36(3):615–627, 1989

1989
[11]

Hwang, A

H. Hwang, A. Rendall, J. Vel´ azquez, Optimal gradient estimates and asymptotic behaviour for th e Vlasov- Poisson system with small initial data , Arch. Ration. Mech. Anal. 200(1):313–360, 2011

2011
[12]

Iacobelli, S

M. Iacobelli, S. Rossi, K. Widmayer, On the stability of vacuum in the screened Vlasov-Poisson eq uation, arXiv:2410.17978
[13]

Ionescu, H

A. Ionescu, H. Jia, Inviscid damping near the Couette ﬂow in a channel , Comm. Math. Phys., 374:2015- 2096, 2020

2015
[14]

Ionescu, H

A. Ionescu, H. Jia, Nonlinear inviscid damping near monotonic shear ﬂows , Acta Math. 230:321–399, 2023

2023
[15]

A. D. Ionescu, B. Pausader, X. Wang, K. Widmayer, On the Asymptotic Behavior of Solutions to the Vlasov–Poisson System , International Mathematics Research Notices. IMRN, 8865– 8889, 2022

2022
[16]

A. D. Ionescu, B. Pausader, X. Wang, K. Widmayer, Nonlinear Landau damping and wave operators in sharp Gevrey spaces , arXiv:2405.04473

arXiv
[17]

Masmoudi, Z

N. Masmoudi, Z. Zhao, Nonlinear inviscid damping for a class of monotone shear ﬂow s in ﬁnite channel , Ann. of Math. (2) 199:1093–1175, 2024

2024
[18]

Mouhot, C

C. Mouhot, C. Villani, On Landau damping , Acta Math., 207:29–201, 2011

2011
[19]

D. Wei, S. Yang, On the 3D relativistic Vlasov-Maxwell system with large Max well ﬁeld , Comm. Math. Phys. 383(3):2275–2307, 2021

2021
[20]

Yukawa, On the interaction of elementary particles , Proc

H. Yukawa, On the interaction of elementary particles , Proc. Phys. Math. Soc. Japan 17:48–57, 1935. School of Mathematical Sciences, Peking University, Beiji ng 100871, China Email address : jnwdyi@pku.edu.cn

1935

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.

[1] [1]

Bardos, P

C. Bardos, P. Degond, Global existence for the Vlasov–Poisson equation in 3 space variables with small initial data , Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire 2(2):101–118, 1985

work page 1985

[2] [2]

Bedrossian, A brief introduction to the mathematics of Landau damping , arXiv:2211.13707

J. Bedrossian, A brief introduction to the mathematics of Landau damping , arXiv:2211.13707

arXiv

[3] [3]

Bedrossian, N

J. Bedrossian, N. Masmoudi, C. Mouhot, Landau damping: paraproducts and Gevrey regularity , Ann. PDE, 2, pp. Art. 4, 71, 2016

work page 2016

[4] [4]

Bedrossian, N

J. Bedrossian, N. Masmoudi, C. Mouhot, Landau damping in ﬁnite regularity for unconﬁned systems wi th screened interactions, Comm. Pure Appl. Math., 71, 537–576, 2018

work page 2018

[5] [5]

Bedrossian, N

J. Bedrossian, N. Masmoudi, Inviscid damping and the asymptotic stability of planar she ar ﬂows in the 2D Euler equations , Publ. Math. Inst. Hautes ´Etudes Sci., 122, 195-300, 2015

work page 2015

[6] [6]

Choi, S.-Y

S.-H. Choi, S.-Y. Ha, H. Lee, Dispersion estimates for the two-dimensional Vlasov-Yuka wa system with small data , J. Diﬀerential Equations 250(1):515–550, 2011

work page 2011

[7] [7]

Duan, Sharp Decay Estimates for the Vlasov-Poisson and Vlasov-Yu kawa Systems with Small Data , Kinet

X. Duan, Sharp Decay Estimates for the Vlasov-Poisson and Vlasov-Yu kawa Systems with Small Data , Kinet. Relat. Models 15(1):119–146, 2022

work page 2022

[8] [8]

Glassey, W

R. Glassey, W. Strauss. Singularity formation in a collisionless plasma could occu r only at high velocities . Arch. Rational Mech. Anal., 92(1):59–90, 1986

work page 1986