Pith. sign in

REVIEW 4 major objections 5 minor 15 references

Euler-Poisson equations with velocity-dependent damping

T0 review · 4 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

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

Pith's one-line read Power-law damping saves Euler-Poisson smoothness only in dimensions 1 and 4.

desk verdict The 1D threshold analysis is solid and new, but the multidimensional theorem's proof is not established: the Floquet claim is applied to the wrong truncated equation, and the conclusion is overbroad. read the letter →

arxiv 2607.18035 v2 pith:26M3FB2M submitted 2026-07-20 math.AP math-phmath.MP

Olga S. Rozanova This is my paper · ORCID
classification math.APmath-phmath.MP MSC 35Q6035L6035L6734M10
keywords Euler-Poissonequationsvelocity-dependentdampingcoldplasmablow-upcriticalthresholdradialsymmetryFloquettheorysmoothnessdomain
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

Repulsive Euler-Poisson equations describe cold plasma, and a key question is which initial data stay smooth for all time. This paper shows that power-law velocity-dependent damping, ν0|V|^k, does something surprising: in one spatial dimension it enlarges the smoothness domain, yet in all other dimensions except four it changes nothing, so arbitrarily small generic perturbations of the zero state still blow up in finite time even though the oscillation amplitude decays. The paper derives a sharp threshold curve in 1D, gives the long-time asymptotics of the solution and its derivatives, and explains the multidimensional failure as a race between polynomial decay from damping and exponential Floquet growth of derivative oscillations.

What carries the argument

The working tool is Radon's lemma, which linearizes the Riccati-type ODEs for the derivatives along characteristics: derivatives are written as ratios p/Q, so blow-up is exactly the first zero of an auxiliary function Q(t). In d dimensions Q is built from a damped Hill-type oscillator. The decisive step is the substitution p1=H exp(-1/2∫a(V)) exp(-(d+2)/2∫F), which converts the linearized equation into H''+JH=0 with J=1-(d+2)/2 ε cos t+o(ε). Floquet theory then decides: for d≠1,4 the characteristic exponent is real, giving exponential growth of H; the damping factor only contributes t^{-κ ε^k}, too weak to stop Q from vanishing. For d=1,4 the exponent is non-real and the polynomial factor wi

What would settle it

Take d=2, k=2, ν0=1 and a radial perturbation with amplitude ε=10^-6; integrate the six characteristic ODEs (23) and (26) numerically for long time. If Q(t) never hits zero, Theorem 4 is false. Conversely, for d=4, any small perturbation that blows up would falsify the claimed exceptional behavior.

Watch

Extended reading notes

Core claim

The central claim is Theorem 4: for d≠1,4, with radial symmetry and damping ν0|V|^k with even k, every sufficiently small generic perturbation of the nontrivial steady state loses smoothness in finite time; damping does not create a neighborhood of smoothness. For d=1 the opposite holds: the threshold curve near the critical point has the form δ=ν0 C_k ε^k, where ε is the initial electric-field perturbation and δ is the excess of the initial density derivative above 1/2, so the safe region genuinely widens. For d=4 the paper argues the behavior matches d=1, though no analytical threshold is given. Theorems 1 and 2 give the 1D asymptotics: velocity oscillations decay like t^{-1/k}, derivative

Load-bearing premise

The proof imports the fact that for d≠1,4 the Hill equation H''+(1-(d+2)/2 ε cos t)H=0 has a real characteristic exponent, and then assumes that the resulting exponential growth of p1 forces Q(t)=1+∫p1 to reach zero; the actual crossing of Q is asserted but not proved.

Editorial extensions

If this is right

  • In one dimension the safe region is strictly larger than without damping, and its boundary for small data is determined by the power k: δ≈ν0 C_k ε^k, with C_k increasing in k.
  • In dimensions 2, 3, 5, and higher the theorem rules out any neighborhood of zero that stays smooth, even under power-law damping; the blow-up time shrinks as d and k grow.
  • Amplitude decay of velocity does not imply stability of derivatives: oscillations of the solution can look damped while their spatial derivatives are already singular.
  • For constant damping, all dimensions gain a smoothness neighborhood; switching to velocity-dependent damping removes that gain except in dimension 4, so the form of the friction term is mathematically decisive.
  • Dimension 4 is exceptional and qualitatively one-dimensional: damping should enlarge its smoothness domain, but the exact threshold is likely very complicated.

Reading between the lines

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

  • If the even-exponent restriction k=2m is merely technical, as the paper's numerics suggest, the same dichotomy—dimensions 1 and 4 safe, other dimensions unsafe—should hold for non-integer k; this is testable by direct simulation for k=1, the aerodynamic-friction case.
  • The theorem's 'general perturbation' should be read as excluding the measure-zero simple-wave class; the paper itself notes that such initial data may stay smooth, so a sharper statement is that the generic small perturbation blows up.
  • An unexamined consequence is that in physical cold-plasma experiments with collisions, observing damped oscillation amplitudes in 2D or 3D radial setups would not be evidence of stability; derivative blow-up could still occur.
  • The large-amplitude behavior in 1D suggests a different route to stability: a strong initial velocity acts as temporary large damping and delays blow-up, so a threshold curve at large E0 may cease to follow the ε^k law and instead show a return of smoothness; the numerics indicate this direction.
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.

Referee Report

4 major / 5 minor

Summary. The paper studies the repulsive Euler-Poisson system with power-law velocity-dependent damping ν(|V|)=ν0|V|^k. In one space dimension, it derives long-time asymptotics for V and its derivatives along characteristics and, for the case V(0)=0, constructs a threshold curve separating initial data that lead to global smoothness from those that blow up; the leading-order form δ=ν0 C_k ε^k is computed analytically. In the radially symmetric multidimensional case, the paper claims that for d≠1,4 every sufficiently small general perturbation of the zero equilibrium blows up in finite time, despite the oscillation amplitude decaying algebraically, while dimensions 1 and 4 are exceptional. The multidimensional proof uses the Radon lemma linearization and a Floquet argument, citing the author's previous work [11] for the decisive characteristic exponent.

Significance. If the multidimensional claim is correct, it is a striking and nontrivial result: unlike constant damping, power-law velocity damping does not create a smoothness neighborhood in most dimensions. The one-dimensional threshold formula is explicit and parameter-free, and the paper gives a falsifiable prediction about the special role of d=1,4. However, the multidimensional proof is only a sketch, and the key Floquet step is not established in the text. The paper also benefits from explicit numerical illustrations of the threshold curve.

major comments (4)
  1. [§3.1, proof of Theorem 4, after Eq. (32)] The Floquet step is not justified. The proof keeps only \bar J=1-((d+2)/2)ε cos t and cites [11] for the assertion that the characteristic exponent μ is real for d≠1,4. But for this leading-order Hill equation, the coefficient A=(d+2)/2 is nonzero for every d≥1; a Mathieu-type analysis gives a nonzero real μ (of order ε^2 at this resonance) for all d, including d=1 and 4. The dimension-dependence that makes d=1,4 exceptional must come from the O(ε^2) terms discarded as o(ε), in particular -((d-2)(d-4))/4 F^2 in J. The proof as written therefore does not establish the claimed dichotomy; it must include the O(ε^2) terms in the Floquet computation or quote a precise statement from [11] for the full J.
  2. [§3.1, Theorem 4, transition from p1 growth to Q(t) vanishing] The assertion that exponential growth of p1 forces Q(t)=1+∫_0^t p1(τ)dτ to vanish is unproved. At best the argument shows p1(t) behaves like e^{μt} times a periodic factor times a power-decaying factor t^{-κ ε^k}. An oscillatory integrand can remain above -1 if the phase is unfavorable, and the periodic factor has zeros. No control of the phase or of the amplitude relative to the initial value 1 is provided. A rigorous proof requires an asymptotic for Q(t), not merely a lower bound for p1(t).
  3. [Theorem 4 vs. Remark 4] The theorem states that blow-up occurs 'for all arbitrary small general perturbation of nontrivial initial data,' but Remark 4 acknowledges that simple-wave solutions with F=F(G) may be globally smooth, and such data form a codimension-one set in the radially symmetric data space. Since these data can be arbitrarily small, the universal statement is internally inconsistent as written. The theorem must be restricted to a generic set excluding the invariant simple-wave branch, or the wording must be changed accordingly.
  4. [§2.2, proof of Theorem 3] The defined function Q1 does not satisfy the linearized equation (20). Direct substitution gives Q1''+Q1 = (k+1)ν0 s0 ((k+3)/(k+2)) sin^{k+3} t, whereas the forcing obtained from Q0 and |V|^k≈ε^k sin^k t is (k+1)ν0 s0 sin^{k+1} t. The displayed expression is therefore incorrect. Although the value at t=π used for the threshold constant happens to be correct, the derivation must be redone and the claim that \bar Q attains its minimum at t=π must be checked for the correct Q1.
minor comments (5)
  1. [§2, Eq. (7) and Lemma 1] The damping coefficient sign is inconsistent: the paper defines ν(|V|)=ν0|V|^k>0, but then writes ν(|V|)=-ν0|V|^k in (7) and in Lemma 1. Presumably the minus sign belongs in the momentum equation, not in the definition of ν.
  2. [§2, Eq. (16)] The transformed equation appears to have a typo: the standard change p2=uv gives \ddot u+(1-φ)u=0, not \ddot u+(1-φ)\dot u=0.
  3. [§2, Eq. (20)] Equation (20) contains a sign error: from (15), Q''+a(V)Q'+Q-1+s0=0, so the term -1-s0 should be -1+s0.
  4. [§3, Lemma 2] The claim that S>0 for d=2 is false for G close to 1/2, since then 1-2G is small positive and ln(1-2G)+1 is negative. The lemma should be restricted to small perturbations, which is the only case used later.
  5. [§3, equation for p1 in Theorem 4] The term (da(V)-2ν(V))ν a'(V)\dot V is dimensionally unclear and appears typographically corrupted. The subsequent definition of J also mixes a(V), ν(V), and a'(V)\dot V in a way that should be written consistently.

Circularity Check

1 steps flagged · score 3.0 of 10

Multidimensional blow-up theorem leans on the author's prior [11] for the decisive Floquet exponent; no fitted/definitional circularity found in the 1D threshold derivation.

  1. self citation load bearing [Section 3.1, Proof of Theorem 4 (after Eq. (32)); see also last paragraph of Introduction]
    "This paper makes extensive use of the methodology developed in [11] ... We do not provide all the details to avoid repeating the steps taken in previous works. As is shown in [11], the characteristic exponent μ for J=Jbar in the cases d,1,d,4 is real. Thus, the amplitude of oscillations of H rises exponentially."

    Theorem 4's d≠1,4 blow-up classification is not proved in the present paper. The decisive step is the reality of the Floquet exponent μ for the truncated Hill equation Jbar=1−((d+2)/2)ε cos t, which is asserted by citation to the author's own [11]. The exponential growth of H, and hence the conclusion that p1 grows and Q(t) must vanish, is inherited from that citation; the paper adds only the polynomial decay factor from damping. Thus the multidimensional result's central premise is load-bearing self-citation rather than an independent derivation, although no fitted parameter is renamed as a prediction.

full rationale

The 1D analysis (Theorems 1–3) is largely self-contained: the threshold constant C_k is computed from an explicit asymptotic expansion Q=Q0+ε^k Q1 and the condition Qbar=0, with no fitted input and no quantity defined in terms of the claimed conclusion. No statistical forcing or renaming of known results appears there. The main circularity-type risk is Theorem 4: its d≠1,4 dichotomy is imported from the author's prior work [11], and the paper explicitly says it omits details already contained in [11],[12]. This is a load-bearing self-citation, but it is not an equivalence by construction: [11] is a prior published result, and the damping analysis adds an independent comparison of polynomial decay versus exponential growth. Hence no definitional circularity is established; the score reflects the heavy self-reliance in the multidimensional proof rather than a circular reduction.

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

No fitted constants enter the central claims; ν0 and k are model inputs and C_k is a computed integral. The main implicit inputs are imported mathematical facts from the author's prior work and formal asymptotics without convergence proof. No new physical entities are postulated.

assumptions (5)
  • domain assumption Floquet characteristic exponent μ for J̄=1−(d+2)ε cos t/2 is real for d≠1,4, yielding exponential growth of H
    Imported from [11]; decisive for Theorem 4 and not derived in this paper.
  • domain assumption Formal asymptotic series (11), (12) may be differentiated termwise and describe the true asymptotics
    Theorems 1–2 give an expansion with secular terms compensated by t^{-n}, but no convergence or error estimates are provided.
  • domain assumption Small deviations satisfy F=ε sin t+o(ε), G=ε cos t+o(ε) uniformly, so Floquet analysis of the linearized periodic equation controls the nonlinear problem
    Used after (31)–(32); uniform-in-time justification is absent.
  • standard math Radon lemma linearization: derivative blow-up is equivalent to vanishing of Q(t)
    Standard Riccati linearization, stated in [7],[8] and applied in (13),(26).
  • standard math Bihari–LaSalle inequality / comparison principle for the amplitude estimates
    Used in Lemmas 1–2 to derive decay bounds.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Euler-Poisson equations with velocity-dependent damping." pith.science (2026). https://pith.science/paper/26M3FB2M

@misc{pith2026260718035,
  author       = {Pith},
  title        = {Pith review of: Euler-Poisson equations with velocity-dependent damping},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/26M3FB2M}},
  note         = {Machine review of arXiv:2607.18035}
}
read the original abstract

We consider repulsive Euler-Poisson equations in both one-dimensional space and in the multidimensional case with radial symmetry, assuming a power-law velocity-dependent damping coefficient. We show that under this assumption, in the one-dimensional case the set of initial data corresponding to a globally time-smooth solution expands, whereas in other dimensions (except for dimension 4) the damping does not influence on improving of smooth properties of the solution. Namely, any arbitrary small perturbation of the steady state blow up, despite of its amplitude of oscillations decays.

Figures

Figures reproduced from arXiv: 2607.18035 by the authors.

Figure 1
Figure 1. k = 2, V0 = 0, ν0 = 1. The threshold curve close to s0 = 0.5 (left) and close to s0 = 1 (right) on the plane (s0,E0)) for E0 > 0. The blow-up domain is under the curve. The picture is symmetric with respect to E0 = 0. Fig.2 shows the decay of the amplitude of E with time for large E0 (the behavior of V is similar, but the graph starts from zero), and the respective behavior of Q(t). Here E0 = 85, s0 = 0.9. Note that… view at source ↗
Figure 2
Figure 2. k = 2, V0 = 0, E0 = 85, s0 = 0.9, ν0 = 1. The behavior of E (left) and Q (right). The initial data that correspond to these solutions are (V,E)| t=0 = (V0(r),E0(r)) = (F0(r)r,G0(r)r), (F0(r),G0(r)) ∈ C 2 (R¯ +). (21) A natural question arises about the class of initial data that lead to a blow-up of the solution in the case of a velocity-dependent damping coefficient for radially symmetric oscillations in a space of… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

15 extracted references

  1. [11]

    Rozanova,On the behavior of multidimensional radially symmetric solutions of the repulsive Euler-Poisson equations, Physica D: Nonlinear Phenomena Vol

    O.S. Rozanova,On the behavior of multidimensional radially symmetric solutions of the repulsive Euler-Poisson equations, Physica D: Nonlinear Phenomena Vol. 443 (2023), 133578

    2023

  2. [1]

    Astarita, G

    G. Astarita, G. Marucci. Principles of non-Newtonian Fluid Mechanics. McGraw Hill, New York, 1974

    1974

  3. [2]

    Stability Theory of Differential Equations,

    R. Bellman, “Stability Theory of Differential Equations,” Dover Books on Mathematics, Mineola, 1953

    1953

  4. [3]

    Chizhonkov, Mathematical aspects of modelling oscillations and wake waves in plasma, CRC Press, 2019

    E.V. Chizhonkov, Mathematical aspects of modelling oscillations and wake waves in plasma, CRC Press, 2019

    2019

  5. [4]

    Choi, D.-H

    Y.-P . Choi, D.-H. Kim, D. Koo, E. Tadmor, Critical thresholds in pressureless Euler-Poisson equa- tions with background states. Ann. Inst. H. Poincar´e Anal. Non Lin´eaire 43 (2026), no. 1, 203–237

    2026

  6. [5]

    C.M.Dafermos, Hyperbolic conservation laws in continuum physics, The 4th Edition, Berlin- Heidelberg: Springer, 2016

    2016

  7. [6]

    Engelberg, H

    S. Engelberg, H. Liu, and E. Tadmor,Critical thresholds in Euler-Poisson equations, Indiana Univer- sity Mathematics Journal, Vol.50 (2001), 109–157

    2001

  8. [7]

    Freiling,A survey of nonsymmetric Riccati equations, Linear Algebra and its Applications, Vol.351-352 (2002), 243–270

    G. Freiling,A survey of nonsymmetric Riccati equations, Linear Algebra and its Applications, Vol.351-352 (2002), 243–270

    2002

Show all 15 references
  1. [8]

    Riccati differential equations,

    W. T. Reid, “Riccati differential equations,” Academic Press, New York, 1972

    1972

  2. [9]

    Rozanova, E

    O. Rozanova, E. Chizhonkov, and M. Delova,Exact thresholds in the dynamics of cold plasma with electron-ion collisions,AIP Conference Proceedings, Vol.2302 (1) (2020), 060012. EULER-POISSON EQUATIONS WITH VELOCITY-DEPENDENT DAMPING 15

    2020

  3. [10]

    Rozanova,Suppression of singularities of solutions of the Euler-Poisson system with density- dependent damping, Physica D: Nonlinear Phenomena Vol.429 (2022), 133077

    O.S. Rozanova,Suppression of singularities of solutions of the Euler-Poisson system with density- dependent damping, Physica D: Nonlinear Phenomena Vol.429 (2022), 133077

    2022

  4. [12]

    Rozanova, and M.I

    O.S. Rozanova, and M.I. Delova,On radially symmetric oscillations of a collisional cold plasma, Math- ematical Methods in the Applied Sciences, Vol.47 (2024), 8385-8399

    2024

  5. [13]

    O.S. Rozanova,Criterion of singularity formation for radial solutions of the pressureless Euler-Poisson equations in exceptional dimension, Journal of Mathematical Analysis and Applications, Vol. 548 (2025) 129394

    2025

  6. [14]

    Rozanova,The repulsive Euler-Poisson equations with variable doping profile, Physica D: Non- linear Phenomena Vol

    O.S. Rozanova,The repulsive Euler-Poisson equations with variable doping profile, Physica D: Non- linear Phenomena Vol. 472 (2025), 134454

    2025

  7. [15]

    Tan,Eulerian dynamics in multidimensions with radial symmetry, SIAM Journal on Mathematical Analysis, Vol

    C. Tan,Eulerian dynamics in multidimensions with radial symmetry, SIAM Journal on Mathematical Analysis, Vol. 53 (3) (2021), 3040–3071

    2021

Pith tools

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