{"id":"7e1c275a-2fbc-4d36-84f5-81a6b56b3ffc","arxiv_id":"2502.02678","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Neutral Vlasov-Poisson plasmas can be constructed to exhibit arbitrarily fast polynomial decay of charge density and electric field, at rates t^{-m-3} and t^{-m-2} for any integer m.","lead":"A mathematical study of plasmas shows that, if the electric field decays quickly enough, the plasma's charge density and field can be made to decay at any chosen polynomial speed. This reveals a whole ladder of possible long-time behaviors for neutral collisionless plasmas, from the standard slow fade to much faster ones.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Existence step in Theorem 1.3 invokes scattering map [14] without verifying its smallness hypotheses or the compact-support hypothesis used by Theorem 1.4; central claim rests on this unverified application.","rationale":"I worked through the conditional analysis: the Taylor expansion in (8), the derivative lemmas in Section 2, and the induction in Theorem 1.4 close with the stated rates, and I found no internal inconsistency in that part of the argument. The genuinely load-bearing weakness is the existence half of Theorem 1.3. The proof constructs final states and cites the scattering map of [14], but every result that converts these final states into the claimed decay rates is conditional on compactly supported initial data and on the limiting moment conditions holding for the actual solution. Neither the smallness hypotheses of [14] nor the compatibility of the produced initial data with Theorem 1.4 is checked. This is not a defect in the conditional estimates but an unverified external applicability step, so a conditional verdict is appropriate rather than acceptance or rejection.","tokens_in":39253,"tokens_out":13826,"duration_ms":144168,"concrete_test":"Read [14, Theorem 1.1(ii) and Remark 1.2(5)] and write out the precise hypotheses, including the norm and smallness threshold for the asymptotic state f_+. Then compute that norm for the constructed states f^α_{m,∞} and for their ε-scalings. If some ε > 0 places all f^α_{m,∞} in the allowed class and the theorem guarantees an initial datum compatible with Theorem 1.4 (compact support, or at least the finite-moment substitutes used in the paper), the existence step is sound. Otherwise the existence claim requires a separate argument before the stated decay rates are established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"To prove the existence part of Theorem 1.3, the paper (Section 3, final paragraph) constructs smooth compactly supported asymptotic states f^α_{m,∞} = φ^α_m(x)ψ_m(v) and asserts that [14, Theorem 1.1(ii) and Remark 1.2(5)] guarantees a global VP solution scattering to them. The cited scattering map is not stated with its full hypotheses, and the paper does not verify that these constructed states lie in the required smallness class, nor that the initial data produced by the map are compactly supported, which is the standing hypothesis of Theorem 1.4. This is load-bearing because the sharp rates claimed in Theorem 1.3 are obtained by applying Theorem 1.4 to the scattering-constructed solutions, and Theorem 1.4 requires both compactly supported initial data and the vanishing moment conditions ρℓ,∞ ≡ 0 for ℓ ≤ n to hold for the actual solution, not merely for the prescribed f∞. Compact support of f∞ does not automatically imply compact support of the initial data obtained from the scattering map, and the moment conditions must be verified for the evolved solution. The gap is probably repairable, for instance by small-amplitude rescaling of the asymptotic states if the [14] norm is homogeneous under such a scaling, but as written this application is unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":39478,"tokens_out":5881,"duration_ms":63017,"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":[{"comment":"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":"Section 3, final paragraph (proof of Theorem 1.3)"},{"comment":"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":"Section 3, final paragraph and Theorem 1.4"},{"comment":"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.","section":"Section 3, proof of Theorem 1.3, m=0 case"}],"minor_comments":[{"comment":"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.","section":"Section 3, proof of Theorem 1.3"},{"comment":"Reference [22] contains the typo 'Arsigmave for rational mechanics and analysis'; this should read 'Archive for Rational Mechanics and Analysis'.","section":"Reference list"},{"comment":"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.","section":"Remark 1.4"},{"comment":"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.","section":"Theorem 1.3, m=0 case"}],"recommendation":"major_revision","confidential_remarks":"The analytic core of the paper (Sections 2–4) appears sound and the main theorem is significant if the existence step can be completed. The load-bearing gap is the unverified use of the scattering map theorem from [14]; this is precisely the kind of issue that should be fixed before publication. I would encourage the editor to ask the authors to state the hypotheses of [14] explicitly and either verify them for their constructed asymptotic states or supply a different existence proof. The paper would also benefit from a short proof of the positivity/charge-decomposition constraint in the m=0 construction. The citation practice is generally appropriate, and the self-citation of [31] for the base case is legitimate rather than circular."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's core achievement is a clean inductive framework showing that neutral Vlasov-Poisson solutions can have charge density and field decay at any prescribed polynomial rate, with linear scattering for the distribution function when the rate exceeds the dispersive one. The mechanism is charge cancellation through vanishing spatial moments of the asymptotic state, not dispersion. That is a genuine step beyond the known dispersive rates and beyond the base case in [31]. The Gegenbauer-polynomial construction of states with prescribed vanishing moments is elegant, and the estimates in Sections 2–4 are carefully and coherently derived. The rates are forced by the first nonvanishing moment, so there is no fitting to a target decay rate.\n\nThe soft spot is where the rubber meets the road: the existence part of Theorem 1.3. The authors build smooth compactly supported asymptotic states and then invoke the scattering map from [14] to produce global solutions that scatter to them. But they never check the hypotheses of that theorem. In particular, the scattering map in [14] is stated for small asymptotic states in a suitable norm, and the constructed states have arbitrary amplitude; nothing is rescaled. Also, the scattering map yields initial data that are not obviously compactly supported, while Theorem 1.4—which supplies the rates—explicitly requires compactly supported initial data. The moment conditions ρℓ,∞ = 0 are verified for the asymptotic state by construction, but Theorem 1.4 needs those conditions to hold for the actual solution's limits F^α,ℓ_∞. That is likely true if the scattering construction is valid, but it is not shown. This is load-bearing: without it, Theorem 1.3 is an assertion about asymptotic states, not about genuine solutions.\n\nThe gap is specific and probably fixable—for example, by rescaling the constructed states to a sufficiently small amplitude if the [14] norm is homogeneous, and by checking or working around the compact-support issue. But as written, the existence claim is unproven. The rest of the paper—the inductive estimates, the derivative bounds, the scattering statements conditional on field decay—looks solid and written by people who know the material.\n\nI would send this to a serious referee. The machinery deserves review and likely publication after a major revision that closes the scattering-map gap. For a reader in kinetic theory or asymptotic analysis, this is a valuable contribution even in its current state. I would cite the inductive framework, but only after the existence issue is sorted out.","headline":"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.","tokens_in":40033,"tokens_out":2725,"would_cite":true,"duration_ms":28486,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q83","35B40","82D10"],"pacs":[],"model":"deepseek-v4-flash","headline":"A neutral plasma can be made to decay at any polynomial rate.","keywords":["Vlasov–Poisson system","collisionless plasma","neutral plasma","polynomial decay rates","charge cancellation","linear scattering","polyhomogeneous expansion","asymptotic profiles"],"falsifier":"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.","tokens_in":39026,"feed_emoji":"⚡","tokens_out":9809,"duration_ms":89904,"temperature":0.7,"pith_summary":"This paper proves that a neutral, collisionless plasma in three dimensions can settle to zero at any prescribed polynomial rate, not just the generic dispersive rates $t^{-3}$ for charge density and $t^{-2}$ for electric field. The rate is set by the first nonvanishing moment of the limiting charge density: if the first $n+1$ such moments vanish, the density decays like $t^{-n-4}$ and the field like $t^{-n-3}$. For every integer $m$, the authors construct nontrivial solutions realizing exactly the rates $t^{-m-3}$ and $t^{-m-2}$; for $m\\ge 1$ the distribution functions scatter linearly to free transport at rate $t^{-m}$, while for $m=0$ scattering is modified by a logarithmic correction. This fills the gap between slow dispersive decay on whole space and exponential phase mixing on the torus, and it exhibits infinitely many distinct asymptotic profiles.","feed_headline":"Plasma decay rates can be tuned to any polynomial order","feed_subtitle":"For any m, the charge density can fall like t^{-(m+3)} and the electric field like t^{-(m+2)}.","key_machinery":"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.","core_discovery":"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$.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"It supplies the prior asymptotic and support-growth results that the new induction starts from, including the field-decay assumption and Theorem 1.1.","marker":"[31]"},{"why":"It provides the scattering-map theorem used to convert each constructed asymptotic state into an actual global solution.","marker":"[14]"},{"why":"It supplies inverse modified scattering and polyhomogeneous expansions that complement the existence construction.","marker":"[38]"},{"why":"It gives the small-data linear scattering and asymptotic behavior that the faster decay rates extend.","marker":"[23]"},{"why":"It supplies the Gegenbauer polynomial orthogonality used to build compactly supported functions with prescribed vanishing moments.","marker":"[1]"},{"why":"It establishes exact large-time behavior for spherically symmetric plasmas, a setting in which the field-decay assumption is known to hold.","marker":"[30]"}],"fun_headline_variants":["Arbitrary polynomial decay rates proven for neutral plasmas","First construction of tuned plasma decay beyond dispersive rates","Plasma decay tuned to any polynomial order: proof","Neutral plasmas can decay at any chosen polynomial rate","Polynomial decay of any order achievable in collisionless plasmas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Arbitrary polynomial decay rates proven for neutral plasmas","First construction of tuned plasma decay beyond dispersive rates","Plasma decay tuned to any polynomial order: proof","Neutral plasmas can decay at any chosen polynomial rate","Polynomial decay of any order achievable in collisionless plasmas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1533,"prompt_tokens":1002,"completion_tokens":531,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":618,"completion_tokens_details":{"reasoning_tokens":451}},"tokens_in":618,"tokens_out":531,"duration_ms":5277,"temperature":1.0,"reasoning_tokens":451,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T11:29:44.290208+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the prior asymptotic and support-growth results that the new induction starts from, including the field-decay assumption and Theorem 1.1."},{"cited_title":"et al., Scattering M ap for the Vlasov-Poisson System, Peking Math J","cited_arxiv_id":null,"evidence_quote":"It provides the scattering-map theorem used to convert each constructed asymptotic state into an actual global solution."},{"cited_title":"International Mathematics Research Notices 2022, 2022: 8865-8889","cited_arxiv_id":null,"evidence_quote":"It gives the small-data linear scattering and asymptotic behavior that the faster decay rates extend."},{"cited_title":"and Stegun, I., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables","cited_arxiv_id":null,"evidence_quote":"It supplies the Gegenbauer polynomial orthogonality used to build compactly supported functions with prescribed vanishing moments."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It establishes exact large-time behavior for spherically symmetric plasmas, a setting in which the field-decay assumption is known to hold."}],"review_version":1}