The Vlasov equation cannot fully account for collisionless shocks

Antoine Bret

REVIEW 3 major objections 3 minor 19 references

The Vlasov equation cannot fully account for collisionless shocks

T0 review · 3 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read The Vlasov equation cannot fully account for collisionless shocks because it conserves entropy, while a shock must generate entropy.

desk verdict Clear restatement of a known entropy objection to Vlasov shocks, but the categorical conclusion is not supported by the argument given. read the letter →

arxiv 2506.01548 v1 pith:GCAUSP46 submitted 2025-06-02 physics.plasm-ph astro-ph.HEastro-ph.SR

Antoine Bret This is my paper

classification physics.plasm-phastro-ph.HEastro-ph.SR PACS 52.35.Tc52.25.Dg

keywords VlasovequationcollisionlessshocksentropyconservationKlimontovichdistributioncontinuumhypothesisplasmakineticsparticle-in-cellsimulationsshock

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

This paper argues that the Vlasov equation, the standard kinetic description of a collisionless plasma, cannot be the right foundation for a theory of collisionless shocks. The reason is an entropy mismatch: the Vlasov equation conserves entropy, while every shock, collisionless or not, must increase entropy because the conservation of matter, momentum, and energy across the front forces a higher-entropy downstream state. The paper therefore proposes that a rigorous mathematical theory of collisionless shocks would have to work at the Klimontovich level, where particles are discrete and no continuum approximation has been made. If correct, this redirects theoretical work on shocks away from smoothed distribution functions and toward discrete-particle or fluctuation-based formalisms, matching what particle-in-cell simulations already do.

What carries the argument

The argument turns on two levels of description: the Klimontovich distribution $F(x,v)=\sum_i \delta(x-x_i)\,\delta(v-v_i)$, which records every particle exactly, and the smoothed kinetic distribution $f(x,v)$ used by the Vlasov equation. The two fixed points are the entropy conservation of Vlasov dynamics and the entropy production required by shock conservation; the proposed escape is the claim that at the shock front the continuum approximation fails, making $f$ inadequate. The constructive machinery is the decomposition $F=f+\delta F$, whose insertion into the Vlasov dynamics produces a fluctuation term like $\langle\delta F\,\delta E\rangle$ that acts as an effective particle-particle collision.

What would settle it

A direct check would be to integrate the Vlasov equation with negligible numerical dissipation for a smooth collisionless shock initial condition and compute $\int f\ln f\,dx\,dv$ across the front. If a steady shock with the predicted downstream entropy emerges from the resolved dynamics, the paper's central claim is wrong; if the solution only produces ever-finer filamentation with no entropy jump in the resolved distribution, the claim is supported.

Watch

Extended reading notes

Core claim

The central claim is a limitation result: because the Vlasov equation conserves the entropy of the smooth distribution function, and because a shock transition necessarily generates entropy through the conservation laws, no solution of the Vlasov equation alone can represent a full collisionless shock. The breakdown is located in the step from the Klimontovich distribution, a sum of delta functions with one term per particle, to the smooth kinetic distribution function. The author argues that at some point along the shock, probably near the front, the real particle distribution becomes too sparse for the continuum hypothesis to hold, and only then can entropy grow in a way that Vlasov dynamics would forbid. The constructive suggestion is to build a theory at the Klimontovich level, for example by splitting the distribution into a kinetic mean plus a fluctuation term, so the resulting correlations play the role of collisions.

Load-bearing premise

The load-bearing premise is that the entropy conserved by the Vlasov equation is the same thermodynamic entropy that must rise across a shock front; if phase-space filamentation and coarse graining can supply the entropy increase, the argument collapses.

Editorial extensions

If this is right

Efforts to prove that the Vlasov equation alone can produce steady collisionless shock solutions would be aiming at an impossible object, because the equation's entropy conservation prohibits the required irreversible transition.
A rigorous theory of collisionless shocks should be formulated for discrete particle distributions, for instance through the Klimontovich equation or a fluctuation expansion around it.
The decomposition $F=f+\delta F$ offers a starting point in which correlated field fluctuations generate an effective collision operator, potentially recovering entropy production without invoking binary collisions.
Particle-in-cell simulations are consistent with the picture because they advance discrete particles with Maxwell's and Newton's equations rather than solving the smoothed Vlasov equation.
Observed entropy generation across the Earth's bow shock and in particle simulations becomes the expected signature of a process that lives at the discrete-particle level, not a puzzle for the continuum model.

Reading between the lines

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

One consequence the author leaves implicit is that the obstruction is not specific to Vlasov: any phase-space continuum description that conserves entropy, including reduced kinetic or fluid closures without a dissipation term, would face the same difficulty at a shock transition.
An extension suggested by the logic is to reconstruct $\int f\ln f\,dx\,dv$ from a PIC shock simulation as the particle number per cell is varied; if the fine-grained entropy production vanishes as the initial distribution becomes smoother, the discreteness is indeed the mechanism, while a finite limit would point back to the continuum.
The proposed decomposition also invites a quantitative test: compute $\langle\delta F\,\delta E\rangle$ from simulation data and compare its integral across the front with the measured entropy jump, turning the qualitative argument into a closure relation.

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 restates a well-known tension rather than resolving it. The observation that Vlasov conserves fine-grained entropy while shocks are entropy-generating is in the textbooks, and the suggested way out—working at the Klimontovich level or with fluctuations—appears in standard kinetic fluctuation literature. What the paper does well is lay out the three levels of description (Klimontovich, Vlasov, fluid) with a clean figure and make the argument easy to grasp quickly. It also honestly cites the Parks et al. result that entropy density across the bow shock matches Vlasov predictions, which is the right piece of evidence to raise.

The soft spots are load-bearing. The central claim relies on identifying the entropy conserved by Vlasov with the thermodynamic entropy that jumps at a shock. That identification is exactly what needs scrutiny. Vlasov dynamics conserves fine-grained entropy, but coarse-grained entropy can grow through phase-space filamentation, and this is the standard way people get irreversible behavior while keeping Vlasov. The paper notes this distinction in passing but does not engage with it, and the categorical conclusion that the Klimontovich level is required would need a proof that no coarse-graining of a Vlasov solution can supply the entropy jump. No such proof appears.

The second soft spot is the assertion in Section 3 that the continuum hypothesis fails at the shock front. That is plausible in some regime, but it is asserted, not derived. There is no scaling argument or estimate of when the phase-space density becomes too sparse. Smooth kinetic models used in shock simulations routinely represent very fine structure, so one cannot just assume a breakdown without more specifics.

All that said, the paper is a clear and honest statement of a concern. It is not a dishonest paper, and its conclusion might be correct in a qualitative sense. But as it stands, it is a perspective piece, not a proof. It does not give the reader a new equation, a numerical experiment, or a quantitative estimate that would move the discussion forward. I would not send it to a serious referee; it is too thin and the central argument has a known gap. It might work as a comment or a blog-style piece, but for a journal it needs either a rigorous treatment of the coarse-graining issue or a concrete example where Vlasov demonstrably cannot reproduce the shock entropy jump.

Referee Report

3 major / 3 minor

Summary. The manuscript argues that the Vlasov equation cannot fully account for collisionless shocks because, while the Vlasov equation conserves entropy, a shock necessarily produces entropy. The author consequently proposes that a rigorous mathematical theory of collisionless shocks may require working at the Klimontovich (discrete-particle) level. Section 2 presents the entropy-conservation argument; Section 3 asserts that the continuum hypothesis fails near the shock front, making the Klimontovich level necessary; the conclusion cites PIC simulations as evidence that particle-level descriptions can produce entropy growth.

Significance. If the central claim were correct, it would overturn the standard kinetic description of collisionless shocks and redirect theoretical work toward discrete-particle formalisms. The paper is concise and references relevant literature, including the Earth's bow-shock entropy measurements and recent PIC studies. However, the claim is not supported by the argument presented. The entropy-conservation point is a textbook statement about fine-grained Vlasov entropy, and standard phase-mixing theory shows that coarse-grained entropy can increase even when fine-grained entropy is conserved. The asserted breakdown of the continuum hypothesis is given without any scaling estimate. The paper therefore does not establish its headline conclusion, although the proposed avenue of studying fluctuations about the Vlasov mean (Eq. 3) is a legitimate research direction.

major comments (3)

[Section 2] The central syllogism conflates the fine-grained Vlasov entropy S_fg = -∫ f ln f dx dv, which is conserved by Vlasov dynamics, with the thermodynamic entropy that increases across a shock. In Vlasov dynamics, phase-space filamentation moves entropy to unresolved scales, so a coarse-grained entropy computed on any finite phase-space cell can increase even though the fine-grained entropy is constant. The paper itself notes this distinction for the Earth's bow shock when it says that the measured quantity is the entropy density, not the total entropy, but then drops the distinction and uses plain 'entropy' to make the categorical claim. To support the conclusion, the author would need to prove that no physically appropriate coarse-graining of a Vlasov solution can produce the required entropy jump across the shock; no such proof is given, and standard kinetic theory of collisionless shocks indicates that coarse-grained entropy does increase through phase mixing and wave-particle interactions.
[Section 3] The claim that 'the Klimontovich distribution ... has to become so sparse that the continuum hypothesis fails' at the shock front is asserted without any quantitative support. No definition of 'sparse' is given, and no estimate is provided for the number of particles per relevant phase-space volume at the scales of the shock transition. A smooth continuum distribution can represent arbitrarily fine structure, and Vlasov solutions are routinely used to model collisionless shocks without invoking a breakdown of the continuum hypothesis. Without a concrete criterion for the failure of the continuum approximation, the conclusion that a rigorous theory 'requires working at the Klimontovich level' does not follow from the paper's reasoning.
[Section 3] The statement that 'PIC simulations do work at the Klimontovich, particle, level' is inaccurate as a description of how PIC codes operate: standard PIC uses finite-size macroparticles and grid interpolation, not delta-function point particles, and introduces numerical approaches to discretization and smoothing. While PIC does not strictly solve the Vlasov equation, the reason entropy can increase in PIC simulations is more subtle than the paper suggests, involving numerical coarse-graining and finite particle statistics. This is a supporting point, but it further weakens the illustrative argument.

minor comments (3)

[Introduction] There are grammatical slips: 'observe particles acceleration' should be 'observe particle acceleration' or 'observe particles' acceleration', and 'via de Vlasov equation' in Section 2 should be 'via the Vlasov equation'.
[Section 2] The phrase 'entropy generation has been measured across the Earth's bow shock, and found consistent with a Vlasov model of the entropy density' deserves a more explicit explanation of how a Vlasov model yields a nonzero entropy density increase if the total Vlasov entropy is conserved; this is a point where the paper could have clarified the coarse-graining issue rather than leaving it implicit.
[Figure 1 and Section 2] The distinction between the 'Klimontovich' and 'kinetic' levels is explained in the text but the figure caption could state that the 'kinetic' level assumes a continuum distribution, to avoid confusion when the text notes that the literature often uses 'kinetic' for both levels.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the central entropy argument rests on external standard references, and self-citations are illustrative, not load-bearing.

full rationale

The paper's chain is: Vlasov conserves entropy (Landau-Lifshitz), shocks increase entropy (Zel'dovich-Raizer and Landau-Lifshitz), so Vlasov alone cannot supply the required rise; the proposed resolution is that the continuum approximation fails at the shock front, requiring Klimontovich. No step is defined in terms of the conclusion. The author's prior Refs. [7] and [11] are used only as illustrative support, not to forbid alternatives or import a uniqueness theorem. No fitted parameter is renamed as a prediction. The empirical bow-shock measurement is explicitly said to concern entropy density, not total entropy. The main weakness is an unsupported assertion about the continuum hypothesis failing, which is a correctness risk rather than circularity.

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

No free parameters or invented entities. The paper rests on the entropy-conservation property of Vlasov, the entropy-increasing property of shocks, and an unproven assertion that the continuum picture fails at the shock front.

assumptions (4)

standard math The Vlasov equation conserves entropy
Cited from Landau-Lifshitz, Physical Kinetics §27; used as premise in Sect. 2.
domain assumption A shock must raise entropy
Standard thermodynamics, cited from Landau-Lifshitz Fluid Mechanics §85 and Zel'dovich and Raizer; used as premise in Sect. 2.
ad hoc to paper The Klimontovich distribution becomes so sparse at the shock front that the continuum hypothesis fails
Stated in Sect. 3 as 'At one point of the motion along the shock, probably about the shock front, the Klimontovich distribution ... has to become so sparse that the continuum hypothesis fails.' No quantitative estimate or proof is given.
domain assumption PIC simulations operate at the Klimontovich level and can therefore capture entropy production
Argued in Sect. 3 based on the discrete nature of PIC; supported by Ref. [7] which includes the author.

how reviews work

0 comments

Cite this review

Pith. "Pith review of The Vlasov equation cannot fully account for collisionless shocks." pith.science (2026). https://pith.science/paper/GCAUSP46

@misc{pith2026250601548,
  author       = {Pith},
  title        = {Pith review of: The Vlasov equation cannot fully account for collisionless shocks},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GCAUSP46}},
  note         = {Machine review of arXiv:2506.01548}
}

read the original abstract

It is argued that the Vlasov equation cannot fully account for collisionless shocks since it conserves entropy, while a shock does not. A rigorous mathematical theory of collisionless shocks could require working at the Klimontovich level.

Figures

Figures reproduced from arXiv: 2506.01548 by the authors.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages

[1]

Birdsall and A

C. Birdsall and A. Langdon. Plasma Physics via Computer Simulation . Series in Plasma Physics and Fluid Dynamics. Taylor & Franci s, 2004. 4

work page 2004
[2]

R. D. Blandford and J. P. Ostriker. Particle acceleratio n by astrophys- ical shocks. Astrophysical Journal Letters , 221:L29–L32, Apr. 1978

work page 1978
[3]

J. Cronin. Cosmic rays: the most energetic particles in t he universe. Reviews of Modern Physics , 71(2, SI):S165–S172, MAR 1999

work page 1999
[4]

R. J. Ewart, M. L. Nastac, P. J. Bilbao, T. Silva, L. O. Silv a, and A. A. Schekochihin. Relaxation to universal non-maxwellian equ ilibria in a collisionless plasma. Proceedings of the National Academy of Sciences , 122(17):e2417813122, 2025

work page 2025
[5]

Fiuza, G

F. Fiuza, G. F. Swadling, A. Grassi, H. G. Rinderknecht, D . P. Hig- ginson, D. D. Ryutov, C. Bruulsema, R. P. Drake, S. Funk, S. Gl enzer, G. Gregori, C. K. Li, B. B. Pollock, B. A. Remington, J. S. Ross , W. Rozmus, Y. Sakawa, A. Spitkovsky, S. Wilks, and H. S. Park. Elec- tron acceleration in laboratory-produced turbulent colli sionless shocks. Nature ...

work page 2020
[6]

Grošelj, L

D. Grošelj, L. Sironi, and A. Spitkovsky. Long-term Evol ution of Rela- tivistic Unmagnetized Collisionless Shocks. The Astrophysical Journal Letters, 963(2):L44, Mar. 2024

work page 2024
[7]

C. C. Haggerty, A. Bret, and D. Caprioli. Kinetic simulat ions of strongly magnetized parallel shocks: deviations from MHD jump condi tions. Monthly Notices of the Royal Astronomical Society , 509(2):2084–2090, Jan. 2022

work page 2022
[8]

Keshet, B

U. Keshet, B. Katz, A. Spitkovsky, and E. Waxman. Magneti c Field Evolution in Relativistic Unmagnetized Collisionless Sho cks. The As- trophysical Journal Letters , 693:L127–L130, Mar. 2009

work page 2009

Show all 19 references

[9]

Landau and E

L. Landau and E. Lifshitz. Course of Theoretical Physics, Physical Kinetics, volume 10. Elsevier, Oxford, 1981

1981
[10]

Landau and E

L. Landau and E. Lifshitz. Fluid Mechanics: Volume 6 . Butterworth- Heinemann, 1987

1987
[11]

Marcowith, A

A. Marcowith, A. Bret, A. Bykov, M. E. Dieckman, L. O’C Dr ury, B. Lembège, M. Lemoine, G. Morlino, G. Murphy, G. Pelletier, I. Plot- nikov, B. Reville, M. Riquelme, L. Sironi, and A. Stockem Nov o. The microphysics of collisionless shock waves. Reports on Progress in Physics,...

2016
[12]

Nicholson

D. Nicholson. Introduction to Plasma Theory . Plasma Physics Series. Wiley, 1983

1983
[13]

G. K. Parks, E. Lee, S. Y. Fu, N. Lin, Y. Liu, and Z. W. Yang. Shocks in collisionless plasmas. Reviews of Modern Plasma Physics , 1(1):1, Dec. 2017. 5

2017
[14]

G. K. Parks, E. Lee, M. McCarthy, M. Goldstein, S. Y. Fu, J . B. Cao, P. Canu, N. Lin, M. Wilber, I. Dandouras, H. Réme, and A. Fazak erley. Entropy generation across earth’s collisionless bow shock . Phys. Rev. Lett., 108:061102, Feb 2012

2012
[15]

M. Pohl, M. Hoshino, and J. Niemiec. PIC simulation meth ods for cos- mic radiation and plasma instabilities. Progress in Particle and Nuclear Physics, 111:103751, Mar. 2020

2020
[16]

R. Z. Sagdeev. Cooperative Phenomena and Shock Waves in Collision- less Plasmas. Reviews of Plasma Physics , 4:23, 1966

1966
[17]

Schroedter and M

E. Schroedter and M. Bonitz. Classical and quantum theo ry of ﬂuctu- ations for many-particle systems out of equilibrium. Contributions to Plasma Physics , 64(5):e202400015, 2024

2024
[18]

Spitkovsky

A. Spitkovsky. Simulations of relativistic collision less shocks: shock structure and particle acceleration. In AIP Conference Proceedings , volume 801, page 345–350. AIP, 2005

2005
[19]

Y. B. Zel’dovich and Y. P. Raizer. Physics of shock waves and high- temperature hydrodynamic phenomena. Dover Publications, 2002. 6

2002

Pith tools

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

[1] [1]

Birdsall and A

C. Birdsall and A. Langdon. Plasma Physics via Computer Simulation . Series in Plasma Physics and Fluid Dynamics. Taylor & Franci s, 2004. 4

work page 2004

[2] [2]

R. D. Blandford and J. P. Ostriker. Particle acceleratio n by astrophys- ical shocks. Astrophysical Journal Letters , 221:L29–L32, Apr. 1978

work page 1978

[3] [3]

J. Cronin. Cosmic rays: the most energetic particles in t he universe. Reviews of Modern Physics , 71(2, SI):S165–S172, MAR 1999

work page 1999

[4] [4]

R. J. Ewart, M. L. Nastac, P. J. Bilbao, T. Silva, L. O. Silv a, and A. A. Schekochihin. Relaxation to universal non-maxwellian equ ilibria in a collisionless plasma. Proceedings of the National Academy of Sciences , 122(17):e2417813122, 2025

work page 2025

[5] [5]

Fiuza, G

F. Fiuza, G. F. Swadling, A. Grassi, H. G. Rinderknecht, D . P. Hig- ginson, D. D. Ryutov, C. Bruulsema, R. P. Drake, S. Funk, S. Gl enzer, G. Gregori, C. K. Li, B. B. Pollock, B. A. Remington, J. S. Ross , W. Rozmus, Y. Sakawa, A. Spitkovsky, S. Wilks, and H. S. Park. Elec- tron acceleration in laboratory-produced turbulent colli sionless shocks. Nature ...

work page 2020

[6] [6]

Grošelj, L

D. Grošelj, L. Sironi, and A. Spitkovsky. Long-term Evol ution of Rela- tivistic Unmagnetized Collisionless Shocks. The Astrophysical Journal Letters, 963(2):L44, Mar. 2024

work page 2024

[7] [7]

C. C. Haggerty, A. Bret, and D. Caprioli. Kinetic simulat ions of strongly magnetized parallel shocks: deviations from MHD jump condi tions. Monthly Notices of the Royal Astronomical Society , 509(2):2084–2090, Jan. 2022

work page 2022

[8] [8]

Keshet, B

U. Keshet, B. Katz, A. Spitkovsky, and E. Waxman. Magneti c Field Evolution in Relativistic Unmagnetized Collisionless Sho cks. The As- trophysical Journal Letters , 693:L127–L130, Mar. 2009

work page 2009