Bridging the Gap between Collisional and Collisionless Plasma Shocks: A Simulation Study using OSIRIS

Antoine Bret; Asaf Pe'er; Kevin M. Schoeffler; Lu\'is O. Silva; Yossef Nissim Kindi

REVIEW 4 major objections 4 minor 4 references

Bridging the Gap between Collisional and Collisionless Plasma Shocks: A Simulation Study using OSIRIS

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

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

Pith's one-line read A single plasma number sets the width of shocks in both collisional and collisionless regimes

desk verdict First simulation scan across the collisional-collisionless shock transition, with a plausible ND-controlled crossover, but the quantitative result rests on an unverified mean-free-path normalization and a fitted prefactor. read the letter →

arxiv 2508.11097 v1 pith:CHX6I5VF submitted 2025-08-14 physics.plasm-ph astro-ph.HE

Yossef Nissim Kindi , Asaf Pe'er , Antoine Bret , Luís O. Silva , Kevin M. Schoeffler This is my paper · ORCID

classification physics.plasm-phastro-ph.HE

keywords plasmashockscollisionlesscollisionalionparametershockwidthparticle-in-cellsimulationCoulombcollisionsbreakout

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 tries to show that a single dimensionless number—the ion plasma parameter $N_D$, roughly the number of ions inside a Debye sphere—sets the width of a plasma shock across the entire range from collision-dominated to collisionless behavior. Using particle-in-cell simulations with a Coulomb-collision module, it measures the 10%–90% density-jump width for shocks formed by two interpenetrating plasma slabs and finds a flat branch near one mean free path for $N_D < 1$ and a falling branch $\ell/(M \ell_{\rm mfp}) \sim 0.7 \ln(N_D)/N_D$ for $N_D > 1$, with a smooth bridge at $N_D \approx 1$. The result matters because it turns the old collisional-versus-collisionless classification into one practical criterion, and it gives astrophysical shock-breakout models a single number to watch for when a shock thins and particle acceleration switches on.

What carries the argument

The central object is the ion plasma parameter $N_D = (4\pi/3) n_0 \lambda_{Di}^3$, the mean number of ions in a Debye sphere, which enters through the collision-rate ratio $\nu_{\rm coll}/\omega_p \sim \ln(1+N_D)/N_D$ once the mean-free-path formula of Eq. (4) is used. The measurement is the shock width $\ell$, defined as the distance over which the ion density rises from 10% to 90% of its downstream value, normalized by the Mach number $M$ and the upstream mean free path $\ell_{\rm mfp}$ and plotted against $N_D$. Plotting the data this way maps the two asymptotic theories onto one horizontal line and one logarithmic curve whose intersection sits at $N_D \approx 1$, making the crossover visible as a single crossing of the two branches.

What would settle it

Run the same two-slab shock at $N_D \approx 3$ in three dimensions with a magnetic field, measure the 10%–90% ion-density-jump width, and normalize it by an independently computed mean free path; if the normalized width does not land near $0.7\ln(3)/3 \approx 0.26$ but instead remains near 1, then the crossover scale at $N_D \approx 1$ is an artifact of the quasi-one-dimensional electrostatic setup or of the mean-free-path normalization rather than a general property of plasma shocks.

Watch

Extended reading notes

Core claim

The central claim is that the ion plasma parameter $N_D = (4\pi/3) n_0 \lambda_{Di}^3$ governs the transition between collisional and collisionless shock dissipation in unmagnetized, nonrelativistic plasmas. In the collisional limit $N_D < 1$, the measured shock width obeys $\ell/(M \ell_{\rm mfp}) \approx 1$, matching the Mott-Smith ansatz with a BGK collision operator; in the collisionless limit $N_D > 1$, it follows $\ell/(M \ell_{\rm mfp}) \approx 0.7 \ln(N_D)/N_D$, matching Tidman's electrostatic shock formalism. The intermediate runs around $N_D \approx 1$ show a smooth, monotone crossover, which the paper presents as numerical validation of the proposed unifying scenario and as evidence that $N_D$, not density or temperature alone, is the right metric for classifying shocks.

Load-bearing premise

The load-bearing premise is that the mean-free-path formula of Eq. (4) is the true ion-ion mean free path for every simulated plasma, including the low-$N_D$ runs where fewer than one ion on average occupies a Debye sphere; both axes of the central plot are normalized by this quantity, so if it is wrong there, the flat collisional branch and the crossover at $N_D \approx 1$ could be normalization artifacts rather than real shock physics.

Editorial extensions

If this is right

Below $N_D \approx 1$, shock widths sit near one mean free path, so the dissipative structure is set by binary ion-ion collisions.
Above $N_D \approx 1$, widths fall as $0.7\ln(N_D)/N_D$ times $M \ell_{\rm mfp}$, so collective electrostatic fields take over the dissipation.
The crossover itself is smooth, meaning there is no critical point separating the two regimes, only a gradual takeover that is already underway at $N_D \approx 1$.
Since density and temperature variations that leave $N_D$ fixed do not change the width, $N_D$ can serve as the organizing control parameter for shock experiments and simulations.
In astrophysical shock breakouts, the upstream density drop pushes $N_D$ across unity, so the shock width collapses and this transition marks the onset of efficient non-thermal particle acceleration.

Reading between the lines

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

A direct test of universality would be a three-dimensional, magnetized run across the same $N_D$ range: if the two-branch collapse survives with a possibly shifted crossover, the $N_D$ criterion is more general than the quasi-1D electrostatic case studied here.
The same criterion suggests a practical rule for laboratory plasma-shock experiments: choose densities and temperatures that put $N_D$ on either side of unity to select the dissipation mechanism on demand.
One useful byproduct would be an interpolation formula for the shock width valid for all $N_D$, which numerical codes could use instead of switching between collisional and collisionless prescriptions.
Because only two electron-to-ion temperature ratios were tested, a continuous scan of this ratio could show whether the crossover sharpens or shifts; a stable crossover at $N_D \approx 1$ would strengthen the paper's central conclusion.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

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

Punchline: this is the first PIC study to scan the shock-width crossover near ND≈1 with a Coulomb collision module, and the qualitative picture—a smooth transition controlled by the ion plasma parameter—looks credible. But the quantitative confirmation is softer than the abstract implies: the collisional-branch normalization uses a Daligault mean free path that is never checked against the collision operator OSIRIS actually implements, and the collisionless-branch curve is a least-squares fit with A=0.7 to the same data it is said to confirm.

What's new: the simulations cover ND from about 1e-2 to 1e5, including the previously unexplored ND=1 regime, with a realistic mass ratio and Debye-resolved grids. The authors show that neither density nor temperature alone orders the shock widths; ND does. That is a genuinely useful consolidation for the shock community.

The good: the setup is clean, the table gives full run parameters, and the authors are upfront about the 2D, quasi-1D, unmagnetized, low-Mach limitations. The claim that the crossover sits near ND~1 is consistent with the data's trend, even if the scatter is large.

Soft spots: first, Eq. (4) is load-bearing. The measured widths are normalized by the Daligault mfp to make the low-ND points sit at y≈1. But the OSIRIS collision module sets its scattering rate via Eq. (2) with an unspecified Coulomb logarithm; nothing in the paper shows that the mfp realized in the simulation matches Eq. (4) across the full ND range. At ND<1, Daligault's formula is an extrapolation into strong coupling, and a mismatch that varies with ND could manufacture the apparent bridge. This is the main issue. Second, A=0.7 is fitted, so the collisionless branch is not an independent test; the functional form is theoretical, but the prefactor is not. Third, there are no error bars or convergence studies (particle number, cell size, run duration); the claim that 4000 1/ωp is sufficient is asserted, not shown.

None of this kills the paper's basic thesis—the smooth transition is probably real—but the quantitative scaling should be treated as provisional. Who is it for? Plasma shock theorists and observers modeling supernova breakout or CME shocks will want to know it exists. It deserves a serious referee, but the authors should be asked to verify the mfp consistency, add error estimates, and run at least one convergence test before the ND criterion is presented as universal.

Recommendation: send to peer review with an expectation of major revision.

Referee Report

4 major / 4 minor

Summary. This manuscript reports 2D PIC simulations with the OSIRIS code including a Coulomb collision module, for two interpenetrating hydrogen plasma slabs. The authors vary density, temperature, and relative drift to cover ion plasma parameters from N_D ≈ 0.009 to 8×10^4, and measure the shock width ℓ as the 10%-90% density rise. They normalize ℓ by M ℓ_mfp and plot versus N_D. The main empirical claim is that the normalized width is ≈1 for N_D < 1 and follows A ln(N_D)/N_D for N_D > 1 with A ≈ 0.7, with a smooth transition near N_D ≈ 1, consistent with the Mott–Smith/BGK prediction in the collisional regime and Tidman's electrostatic-shock scaling in the collisionless regime. The paper proposes N_D as a practical metric for identifying when shock dissipation transitions from collisional to collisionless.

Significance. Should the central claim hold, this would be the first simulation scan to bridge collisional and collisionless plasma shocks with a single dimensionless parameter, with direct implications for shock breakout and particle-acceleration diagnostics. The paper has genuine strengths: a broad parameter scan, a realistic electron-ion mass ratio, use of a well-established collision module, and explicit comparison with two theoretical scalings. The asymptotic forms are plausible and the paper states its own limitations clearly in Section V. However, the quantitative confirmation is currently incomplete: the collisionless curve uses a prefactor fitted to the same measurements, the mean-free-path normalization is not verified against the simulation's own collision operator, no error bars or convergence tests are reported, and the transition region is sparsely sampled. These are addressable in revision.

major comments (4)

[Section III.C, Figure 2] The collisionless-branch comparison curve ℓ/(M ℓ_mfp) = A ln(N_D)/N_D with A = 0.7 is obtained by a least-squares fit to the same measurements it is then used to confirm. Consequently, the reported 'agreement' in that regime tests only the functional form, not the prefactor. Please derive or otherwise fix A from the theory of Bret and Pe'er, state explicitly which points entered the fit, report the fitted value with its uncertainty, and show residuals so the reader can judge whether any N_D-dependence beyond the assumed functional form is present.
[Section III.B, Eq. (4)] The y-axis normalization uses the Daligault mean-free-path formula, Eq. (4), but the paper does not demonstrate that this formula describes the ion-ion mean free path actually realized by the OSIRIS collision module. The collision frequency in Eq. (2) depends on a Coulomb logarithm whose value is not reported, and Eq. (4) originates from diffusion in a strongly coupled one-component plasma. At N_D < 1, where strong-coupling corrections are largest, an extrapolation error in Eq. (4) would shift the low-N_D points off the collisional branch and could create or destroy the apparent crossover near N_D ≈ 1. Please measure the effective mean free path directly in OSIRIS (for example, by test-particle velocity relaxation or momentum decorrelation) across the simulated N_D range and overlay it with Eq. (4).
[Section II.B and Table I] The shock widths are quoted without error bars, and no convergence tests are shown. The 10%-90% definition and the finite simulation time imply both statistical and systematic uncertainties; without these, the scatter of points in Figure 2 cannot be interpreted. Please provide uncertainty estimates for each width, and show at least one convergence test (particles per cell, cell size, number of transverse cells, run duration) for a representative run in each regime, for example one N_D < 1, one N_D ≈ 1, and one N_D > 1 run.
[Section III.C and Table I] The transition region around N_D ≈ 1 is currently sampled by only three runs (N_D = 1.90, 2.78, 6.01) and with no error bars. The word 'smooth' in the abstract and in Section V is therefore not yet supported by the data. Additional runs concentrated at N_D ≈ 0.5-10 are needed to establish the continuity and monotonicity of the bridge between the two asymptotic branches.

minor comments (4)

[Figure 1 caption and Section III.A] The caption calls the N_D = 2.78 case a 'collisional regime' shock while the text describes it as the 'intermediary regime'; please reconcile the terminology.
[Section II.B] The sentence 'All simulations ran for 40001/ωp' appears to contain a typo; presumably 4000 ω_p^{-1} is intended. Please correct.
[Footnote 25] The numerical check stated in footnote 25, that shock width is insensitive to density and temperature separately when N_D is held fixed, is only described in words; please include a supporting figure or table, since this check underlies the claim that N_D is the controlling parameter.
[Section III.B] The simplified scaling ν_coll/ω_p ∼ ln(1 + N_D)/N_D drops the constants from Eq. (4); consider writing the full expression with the constants to avoid confusion with the actual Ln(1 + 6.477 N_D) form.

Circularity Check

1 steps flagged · score 6.0 of 10

The collisionless-branch curve in Fig. 2 is a least-squares fit to the same simulation widths it is then used to validate; the functional form and the collisional branch are independent.

fitted input called prediction [Section III.C, 'Shock width dependence on the plasma parameter', and Fig. 2]
"A least–squares fit to our measurements yields A≃ 0.7; the red dashed line in the figure adopts this value. The measured points track the two analytic branches in their respective limits and exhibit a smooth bridge around ND∼ 1, thereby validating the crossover scenario proposed by Bret and Pe'er9."

The collisionless reference curve ℓ/(M·lmfp)=A ln(ND)/ND is drawn with A obtained by least-squares fitting the very measurements plotted in Fig. 2. Because the Bret & Pe'er theory leaves A undetermined, the red dashed curve is not a parameter-free prediction: its vertical placement is constructed from the data, so the agreement of the high-ND points with this curve is statistically enforced rather than independently confirmed. The functional form ln(ND)/ND remains an external prediction, and the collisional branch y=1 is parameter-free, so the circularity is limited to the amplitude of the collisionless branch and to the statement that the fitted curve 'validates' the scenario.

full rationale

The only reduction found is the fitted prefactor A in the collisionless branch. The paper explicitly discloses the fit, but then labels the resulting curve as 'predicted by Bret and Pe'er' and uses it as validation, which is fitted input called prediction for that branch. The functional form ln(ND)/ND is an independent theoretical scaling (Tidman plus the Daligault mean-free-path formula), and the collisional branch ℓ/(M lmfp)=1 is a parameter-free comparison that the low-ND data genuinely test. The crossover location near ND≈1 is not forced by the fit: the reference curves for A=0.7 do not intersect, so the observed departure of the data from y=1 is an empirical feature. The Daligault mean-free-path formula and the OSIRIS collision module are external inputs; whether Eq. (4) matches the actual ion-ion mfp realized in OSIRIS is a validity/correctness concern, not circularity. The Bret & Pe'er (2021) predictions are a same-author prior theory, but the simulations provide an independent, externally falsifiable test of it, so the self-citation is not load-bearing in a circular sense. The paper's own caveat that temperature-ratio choices could shift the precise crossover ND is a limitation, not a circular step. Overall: partial circularity in the collisionless amplitude, with the central crossover claim retaining independent content; score 6.

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

The central claim rests on the OSIRIS collision model, the Daligault mean-free-path formula used to normalize the comparison, and the quasi-1D electrostatic setup. The only free parameter is the prefactor A fitted to the collisionless branch. No new physical entities are introduced.

free parameters (1)

A = 0.7
Prefactor in the collisionless scaling ℓ/(M lmfp) = A ln N_D/N_D; obtained by least-squares fit to the 14 measured shock widths. The theory (Bret & Pe'er 2021) left A unspecified, so the agreement of the data with the resulting curve is partly by construction.

assumptions (5)

domain assumption OSIRIS's Coulomb collision module correctly reproduces the Landau collision operator in the small-angle limit (Eq. 2).
The collision operator is implemented via random small-angle deflections per Takizuka-Abe; the paper relies on this to model collisional dissipation.
domain assumption The Daligault mean-free-path formula (Eq. 4) applies to ion-ion scattering in this two-species plasma across the full N_D range, including N_D < 1.
The formula was derived for a one-component plasma; the paper extends it to the proton-electron plasma and to the strongly coupled regime where it may not be valid. This formula normalizes both axes of Figure 2.
domain assumption A quasi-1D, 2D, unmagnetized, nonrelativistic electrostatic treatment captures the shock physics relevant to the claimed transition.
Stated in the caveats (Section V); electromagnetic instabilities and 3D effects are excluded, which may modify the crossover.
domain assumption The 10%-90% density-ramp width is a faithful measure of the shock dissipation width and is converged at simulation time 4000/ω_p.
No convergence study or sensitivity analysis of the width definition is provided; density oscillations in the downstream could bias the 90% threshold.
domain assumption 36 particles per cell and cell size 0.3 λ_De are sufficient for numerical convergence.
No convergence tests are shown; particle noise may affect the measured widths, especially for the narrow collisionless shocks.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Bridging the Gap between Collisional and Collisionless Plasma Shocks: A Simulation Study using OSIRIS." pith.science (2026). https://pith.science/paper/CHX6I5VF

@misc{pith2026250811097,
  author       = {Pith},
  title        = {Pith review of: Bridging the Gap between Collisional and Collisionless Plasma Shocks: A Simulation Study using OSIRIS},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/CHX6I5VF}},
  note         = {Machine review of arXiv:2508.11097}
}

abstract

Shock waves in plasmas can be characterized by the mechanisms behind their formation. When binary collisions are frequent, dissipation is collision-driven and the shock width is a few mean free paths. In contrast, collisionless shocks rely on collective plasma processes to establish dissipation on scales well below the mean free path. We bridge these regimes with particle-in-cell simulations using OSIRIS with a Coulomb-collision module, varying parameters that control collisionality. We find a smooth transition of the shock width in the intermediate region where the ion plasma parameter $N_D \approx 1$. Our results recover the asymptotic predictions: a collisional-regime width consistent with the Mott-Smith ansatz with a BGK operator, and the collisionless limit consistent with Tidman's classical formalism. The ion plasma parameter thus serves as a practical metric for identifying when shocks shift from fluid-like, mean-free-path scales to collisionless, sub-mean-free-path scales. We discuss implications for astrophysical environments, where shock breakout changes the shock width and marks the onset of efficient particle acceleration.

Figures

Figures reproduced from arXiv: 2508.11097 by the authors.

**Figure 2.** FIG. 2. Measured shock width (in units of the upstream mean-free-path times the Mach num [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

4 extracted references · 3 canonical work pages

[1]

1Y. B. Zel’dovich and Y. P. Raizer,Physics of Shock Waves and High-Temperature Hydro- dynamic Phenomena, dover ed. ed. (Dover Publications, Mineola, New York, 2002). 2L. D. Landau and E. M. Lifshitz,Fluid Mechanics, 2nd ed., Course of Theoretical Physics, Vol. 6 (Pergamon Press, Oxford, 1987). 12 3O. Buneman, Physics of Fluids7, S3 (1964). 4R. Z. Sagdeev,...

work page 1964
[3]

14R. A. Fonseca, L. O. Silva, F. S. Tsung, V. K. Decyk, W. Lu, C. Ren, W. B. Mori, S. Deng, S. Lee, T. Katsouleas, and J. C. Adam, inComputational Science – ICCS 2002, Lecture Notes in Computer Science, Vol. 2331, edited by P. M. A. Sloot, A. G. Hoekstra, C. J. K. Tan, and J. J. Dongarra (Springer, Berlin, Heidelberg,

work page 2002
[1964]

5S. D. Bale, F. S. Mozer, and T. S. Horbury, Phys. Rev. Lett.91, 265004 (2003). 6S. J. Schwartz, E. Henley, J. Mitchell, and V. Krasnoselskikh, Phys. Rev. Lett.107, 215002 (2011). 7A. Bamba, R. Yamazaki, M. Ueno, and K. Koyama, Astrophys. J. 589, 827 (2003), arXiv:astro-ph/0302174 [astro-ph]. 8A. Bret, A. Stockem, F. Fiuza, C. Ruyer, L. Gremillet, R. Nara...

work page Pith review arXiv 2003
[2002]

pp. 342–351. 15R. A. Fonseca, F. S. Tsung, S. D. Martins,et al., Lect. Notes Comput. Sci. 2331, 342 (2005). 16M. Fiore, F. Fiúza, M. Marti, R. A. Fonseca, and L. O. Silva, Journal of Plasma Physics 76, 813 (2010). 17T. Takizuka and H. Abe, Journal of Computational Physics25, 205 (1977). 18K. Nanbu, Phys. Rev. E55, 4642 (1997). 19F. Peano, M. Marti, L. O. ...

arXiv 2005

Pith tools

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

[1] [1]

1Y. B. Zel’dovich and Y. P. Raizer,Physics of Shock Waves and High-Temperature Hydro- dynamic Phenomena, dover ed. ed. (Dover Publications, Mineola, New York, 2002). 2L. D. Landau and E. M. Lifshitz,Fluid Mechanics, 2nd ed., Course of Theoretical Physics, Vol. 6 (Pergamon Press, Oxford, 1987). 12 3O. Buneman, Physics of Fluids7, S3 (1964). 4R. Z. Sagdeev,...

work page 1964

[2] [3]

14R. A. Fonseca, L. O. Silva, F. S. Tsung, V. K. Decyk, W. Lu, C. Ren, W. B. Mori, S. Deng, S. Lee, T. Katsouleas, and J. C. Adam, inComputational Science – ICCS 2002, Lecture Notes in Computer Science, Vol. 2331, edited by P. M. A. Sloot, A. G. Hoekstra, C. J. K. Tan, and J. J. Dongarra (Springer, Berlin, Heidelberg,

work page 2002

[3] [1964]

5S. D. Bale, F. S. Mozer, and T. S. Horbury, Phys. Rev. Lett.91, 265004 (2003). 6S. J. Schwartz, E. Henley, J. Mitchell, and V. Krasnoselskikh, Phys. Rev. Lett.107, 215002 (2011). 7A. Bamba, R. Yamazaki, M. Ueno, and K. Koyama, Astrophys. J. 589, 827 (2003), arXiv:astro-ph/0302174 [astro-ph]. 8A. Bret, A. Stockem, F. Fiuza, C. Ruyer, L. Gremillet, R. Nara...

work page Pith review arXiv 2003

[4] [2002]

pp. 342–351. 15R. A. Fonseca, F. S. Tsung, S. D. Martins,et al., Lect. Notes Comput. Sci. 2331, 342 (2005). 16M. Fiore, F. Fiúza, M. Marti, R. A. Fonseca, and L. O. Silva, Journal of Plasma Physics 76, 813 (2010). 17T. Takizuka and H. Abe, Journal of Computational Physics25, 205 (1977). 18K. Nanbu, Phys. Rev. E55, 4642 (1997). 19F. Peano, M. Marti, L. O. ...

arXiv 2005