{"id":"d2ca4598-a244-4948-83e3-4868ab6bfbbb","arxiv_id":"2508.11097","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Simulations show that the ion plasma parameter N_D, the number of particles in an ion Debye sphere, controls a smooth transition of shock width from many mean free paths to a collisionless scale around N_D = 1.","lead":"This paper uses computer simulations of colliding plasma flows to show that the width of a plasma shock changes smoothly from a collision-dominated scale to a collisionless scale as a single plasma parameter crosses 1. The result gives astrophysicists a simple way to tell whether a shock in space, such as a supernova blast wave, should be treated as a fluid or as a collisionless particle accelerator.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Figure 2's y-axis is normalized by the Daligault mfp (Eq. 4), but the paper never verifies that this formula equals the ion-ion mfp actually realized by the OSIRIS collision module; a mismatch could make the reported collapse and the N_D≈1 crossover an artifact of the normalization.","rationale":"The paper presents a useful simulation campaign across a genuinely difficult intermediate regime, and the qualitative picture—shock widths of order the mfp at low N_D and far below the mfp at high N_D—is likely robust to the precise mfp definition. The use of the established OSIRIS code and a dedicated parameter scan are strengths. However, the headline claim is quantitative: it asserts recovery of the Mott-Smith and Tidman asymptotic width scalings. The only unfitted comparison is the collisional branch y=1, and both branches' functional forms in Fig. 2 are tied to Eq. (4) through the y-axis and through the logarithmic form of the collisionless curve. Since the paper does not specify the Coulomb logarithm used in the collision module or verify that Eq. (4) matches the simulation's actual mean free path, the collapse could in principle be a consequence of the chosen normalization rather than an emergent shock property. The concrete test I propose directly measures the mfp within the same collision model, which would settle whether the normalization is correct. Because the reader's verdict is already CONDITIONAL and identifies the Daligault mfp as the weakest assumption, my read does not change the verdict; it sharpens the reason for caution and sets a precise acceptance check.","tokens_in":164,"tokens_out":8233,"duration_ms":100493,"concrete_test":"Use the same OSIRIS setup (cell size, time step, particles per cell, and the exact Coulomb-log implementation used in the runs) in a uniform, unmagnetized plasma at N_D values representative of the table, e.g., 0.06, 0.28, 2.8, 100, and 8×10^4. Measure the ion-ion mean free path directly, for example by fitting the exponential relaxation of a small ion temperature anisotropy or by tracking the spatial diffusion of a tagged ion population, and compare the inferred lmfp with Eq. (4). If the ratio deviates by more than about 20% at any N_D, recompute the y-axis in Fig. 2 using the simulation-derived mfp and check whether the data still collapse onto the y=1 and y=0.7 ln N_D/N_D curves.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the collapse of the measured shock widths in Fig. 2 onto the collisional branch y=1 and the collisionless branch y=0.7 ln N_D/N_D. The y-axis normalizes the measured width by the Daligault mfp, Eq. (4). This is the only parameter-free comparison: the collisionless branch uses a fitted prefactor A, so the collisional branch y=1 is the only independent theoretical test. The paper never demonstrates that Eq. (4) is the ion-ion mfp actually realized by the OSIRIS collision module, whose collision frequency is set by Eq. (2) with an unspecified Coulomb logarithm. Daligault's formula comes from diffusion in a strongly coupled one-component plasma, not automatically from the two-species, electron-screened PIC plasma in OSIRIS. At N_D<1, strong coupling makes both the Coulomb logarithm and the diffusive mfp concept coupling-dependent, so an extrapolation error is plausible. If the true simulation mfp differs from Eq. (4) by a factor that varies with N_D, the low-N_D points will no longer sit at y=1, and the apparent smooth bridge near N_D≈1 could be an artifact of the normalization rather than an emergent property of the shocks. The reader flagged the extrapolation; the sharper issue is the missing consistency check between the external formula and the simulation's actual collision operator.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8021,"tokens_out":5945,"duration_ms":65486,"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":[{"comment":"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":"Section III.C, Figure 2"},{"comment":"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":"Section III.B, Eq. (4)"},{"comment":"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":"Section II.B and Table I"},{"comment":"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.","section":"Section III.C and Table I"}],"minor_comments":[{"comment":"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":"Figure 1 caption and Section III.A"},{"comment":"The sentence 'All simulations ran for 40001/ωp' appears to contain a typo; presumably 4000 ω_p^{-1} is intended. Please correct.","section":"Section II.B"},{"comment":"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":"Footnote 25"},{"comment":"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.","section":"Section III.B"}],"recommendation":"major_revision","confidential_remarks":"This is a well-motivated study within the scope of Physics of Plasmas, and the authors' caveats in Section V show appropriate awareness of the model's limitations. The main concern is that the quantitative claims are currently under-supported: the fitted prefactor A, the unverified mean-free-path normalization, missing error bars, and the sparse sampling of the transition region all need to be addressed. I believe these points can be fixed with additional analysis and a modest number of new runs, so I recommend major revision rather than rejection. The self-citations to Bret and Pe'er are appropriate given that the paper directly tests their theoretical predictions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nWhat'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.\n\nThe 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.\n\nSoft 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.\n\nNone 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.\n\nRecommendation: send to peer review with an expectation of major revision.","headline":"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.","tokens_in":8645,"tokens_out":4098,"would_cite":true,"duration_ms":41955,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single plasma number sets the width of shocks in both collisional and collisionless regimes","keywords":["plasma shocks","collisionless shocks","collisional shocks","ion plasma parameter","shock width","particle-in-cell simulation","Coulomb collisions","shock breakout"],"falsifier":"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.","tokens_in":7528,"feed_emoji":"⚡","tokens_out":13734,"duration_ms":120217,"temperature":0.7,"pith_summary":"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.","feed_headline":"Shock width obeys one plasma number across all regimes","feed_subtitle":"Particle-in-cell runs show width stays near a mean free path below N_D=1 and drops as ln(N_D)/N_D above it.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the proposed crossover scenario and the asymptotic shock-width scalings that the simulations are compared with.","marker":"[9]"},{"why":"Provides the Mott-Smith collisional shock model behind the few-mean-free-path width expected at low $N_D$.","marker":"[10]"},{"why":"Provides the classical electrostatic collisionless shock formalism behind the high-$N_D$ width scaling.","marker":"[12]"},{"why":"Supplies the mean-free-path formula (Eq. 4) used to normalize measured widths and to convert $N_D$ into a collision-rate ratio.","marker":"[24]"},{"why":"Original binary Coulomb collision model by random small-angle scattering that the collision module's deflection statistics follow.","marker":"[17]"},{"why":"Describes the Coulomb-collision module added to the particle-in-cell code, making the collisional runs possible.","marker":"[16]"},{"why":"Describes the particle-in-cell code in which all shock simulations were run.","marker":"[14]"},{"why":"Provides the two-slab geometry and hot-electron, cold-ion initial conditions used to form the shocks.","marker":"[21]"},{"why":"Guides the choice of electron-dominated sound speed and upstream conditions for the collisionless runs.","marker":"[26]"}],"fun_headline_variants":["Ion plasma parameter sets shock width transition","N_D predicts shock width from collisional to collisionless","One plasma number unifies shock dissipation regimes","Shock width scaling follows ion plasma parameter","Plasma parameter controls shock width across collisionality"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Ion plasma parameter sets shock width transition","N_D predicts shock width from collisional to collisionless","One plasma number unifies shock dissipation regimes","Shock width scaling follows ion plasma parameter","Plasma parameter controls shock width across collisionality"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000179,"raw_usage":{"total_tokens":1297,"prompt_tokens":939,"completion_tokens":358,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":287}},"tokens_in":555,"tokens_out":358,"duration_ms":3907,"temperature":1.0,"reasoning_tokens":287,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:28:18.805930+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":2}