{"id":"04ccb3f4-244b-4db6-8fe4-cc94b84798cf","arxiv_id":"2412.03530","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using a reduced ion-to-electron mass ratio of 100 instead of the realistic 1836 in 1D3V PIC simulations of strong parallel shocks reverses electron acceleration behavior and changes electron and ion heating at low and high Alfvén Mach numbers.","lead":"A particle-in-cell simulation study shows that using a reduced ion-to-electron mass ratio instead of the real value (1836) changes how energy is split between electrons and ions in strong collisionless shocks, affecting both heating and particle acceleration. The study compares mass ratios of 100 and 1836 at two Alfvén Mach numbers and concludes that the realistic ratio is needed to capture electron acceleration and the electron-to-ion temperature ratio.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The high-MA comparison does not isolate mass-ratio scale separation from the factor-4.3 difference in upstream electron thermal velocity; the causal attribution is therefore not yet controlled.","rationale":"I agree with the reader that the upstream electron thermal-spread difference is the most load-bearing unresolved issue. The setup fixes energy temperatures, so the mr=100 and mr=1836 runs differ not only in ion-to-electron mass ratio but also in v_te by about a factor of 4.3 (Sec. 2 and footnote 3). The low-MA result has an independent mechanistic anchor: the intermediate-scale instability is stable for mr=100 and unstable for mr=1836 (Sec. 3.3), so that finding is robust to the confound. But the high-MA result, where no instability threshold explains the difference, relies entirely on the mr=100 versus mr=1836 pair; without holding v_te fixed, the 'over-efficient electron acceleration' attributed to scale separation could be a seed-temperature effect. I considered the two-point 'independence' of Te/Ti as a candidate concern, but it is a secondary claim and would not overturn the main mass-ratio conclusion. I also weighed the finite simulation time; the paper's caveat about simulating only the first ~1.5 minutes is explicit, and it weakens quantitative convergence claims, not the qualitative early-time comparison. Therefore the conditional verdict is appropriate, pending the proposed control run.","tokens_in":13224,"tokens_out":12507,"duration_ms":129318,"concrete_test":"Run the MA=21.3 case with mr=100 and upstream electron temperature multiplied by 1836/100 so that v_te matches the mr=1836 run, keeping B0, vu, and the numerical setup fixed (Ms remains strongly supersonic, ~117). Compare the downstream electron spectrum, Te/Ti, and ion high-momentum flux with the existing mr=100 and mr=1836 results. If the mr=100 curves move toward the realistic run, the high-MA discrepancy is partly due to the initial electron thermal spread; if they are unchanged, the scale-separation attribution is supported. A complementary scan varying upstream electron temperature in both mass-ratio runs would bound the effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central comparison (Sec. 2) fixes kBTe = kBTi = 4e-8 mi c^2 while varying mr, so the mr=100 upstream electron thermal velocity is v_te = sqrt(4e-8 mr) c = 0.002c, whereas the mr=1836 runs have v_te = 0.0086c, a factor 4.3 difference. Footnote 3 acknowledges this but does not control for it. The high-MA conclusion that mr=100 produces over-efficient electron acceleration and a nonphysical high-momentum ion flux rests on exactly this pair of runs, and electron injection and heating can plausibly depend on v_te/v_sh. The low-MA suppression has a separate mechanistic anchor in the intermediate-scale instability threshold (Sec. 3.3), so it is less threatened. But no mechanism is given for the high-MA sensitivity, and without a run that changes mr while holding v_te fixed, the paper cannot distinguish 'scale separation' from 'initial electron thermal spread' as the cause of the high-MA differences. This is load-bearing because the title and abstract attribute the observed differences to ion-to-electron scale separation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses 1D3V particle-in-cell simulations with the SHARP code to compare reduced (mr=100) and realistic (mr=1836) ion-to-electron mass ratios in strong, parallel, non-relativistic electron-ion shocks at Alfvén Mach numbers MA=5.3 and 21.3. It reports that the reduced mass ratio leads to qualitatively different particle acceleration and thermal energy dissipation: at low MA it suppresses electron acceleration, while at high MA it produces over-efficient electron acceleration and an unrealistically high ion flux at high momentum. The paper further claims that with the realistic mass ratio the downstream electron-to-ion temperature ratio is independent of the upstream magnetic field, a result not found in reduced-mass-ratio runs. The simulations are shown for long durations (~4e4 omega_i^-1) and made publicly available.","tokens_in":13440,"tokens_out":3401,"duration_ms":35188,"significance":"If the conclusions hold, the paper is an important methodological caution for the PIC simulation community, and it provides a concrete physical result (the approximate constancy of the electron-to-ion temperature ratio with upstream magnetic field) that could be tested by other codes and extended to higher dimensions. The study's strengths include the use of a code specifically designed to minimize numerical heating over long runtimes, the long simulation durations relative to typical kinetic shock simulations, and the public release of visualization movies. The low-MA suppression of electron acceleration is anchored in a specific instability-threshold argument (Section 3.3) and is therefore a credible, mechanism-based claim. The quantitative acceleration efficiencies (11-12% into ions) are consistent with hybrid-simulation results, which lends plausibility. However, the high-MA comparison and the field-independence claim are currently under-controlled, as detailed below.","major_comments":[{"comment":"The high-Mach-number comparison does not isolate the ion-to-electron scale separation from the change in upstream electron thermal velocity. The upstream electron temperature is fixed in energy units (kBTe = 4e-8 mi c^2) in all runs, so when mr is increased from 100 to 1836 the upstream electron thermal velocity increases by a factor sqrt(1836/100) ~ 4.3. Footnote 3 acknowledges this difference, but no simulation controls for it. Since electron injection, heating, and wave-particle interactions can depend on v_te/v_sh, the differences attributed to scale separation in Section 4.2 (electron acceleration and ion flux at MA=21.3) and Section 4.3 (excessive electron heating) could in part be caused by the different initial electron thermal spread. A run with mr=100 and kBTe scaled as 1/mr (so that the electron thermal velocity matches the mr=1836 case), or a scan over v_te at fixed mr, is needed before the title-level attribution to scale separation is justified. The low-MA conclusion is less affected because it is supported by the intermediate-scale instability threshold in Section 3.3, but the high-MA claim is load-bearing and currently under-controlled.","section":"Section 2, footnote 3; Sections 4.2-4.3"},{"comment":"The claim that the electron-to-ion temperature ratio is independent of the upstream magnetic field rests on only two Mach numbers (MA=5.3 and 21.3) with the realistic mass ratio. No error bars or run-to-run variability are provided, and the time series in Figure 5 show secular evolution, so the difference between the two curves could reflect a monotonic trend rather than independence. This claim is highlighted as a central result in the abstract and conclusion. The authors should either soften the wording to 'weakly dependent' or add at least one intermediate-MA simulation and a quantitative estimate of the temporal and sampling uncertainty to support the independence statement.","section":"Section 4.3, Figure 5"},{"comment":"The normalization of momenta, Mach numbers, and temperatures relies on assuming a compression ratio R=4, but the actual downstream compression ratio is not reported, and the text and figures use inconsistent Alfvén Mach numbers (MA=5.3 and 21.3 in the text, MA=5.1 and 21 in Figures 1 and 3). If the actual R differs from 4, then the normalized momenta in Figure 2 and the thermal-energy normalization in Figure 4 shift correspondingly, which could affect the quantitative comparison. The paper should state the measured downstream compression ratio for each run and reconcile the Mach-number labels in the text and figures.","section":"Equations (1) and (10); Figures 1-3"}],"minor_comments":[{"comment":"The phrase 'necessitate' in the abstract should be 'necessitates' to agree with the singular subject 'complexity'.","section":"Abstract and Introduction"},{"comment":"The non-thermal energy threshold ps > 5 ps,max is a free parameter; a brief sensitivity check (e.g., using 3 and 10 ps,max) would strengthen the quantitative acceleration-efficiency claims, even though the qualitative conclusions appear robust.","section":"Section 4.4, Figure 6"},{"comment":"The concluding statement that the time evolution of electron acceleration efficiency is 'contrary to what is incorrectly inferred from simulations run for much shorter physical times in the literature' is a strong claim that should be supported by a specific citation or a direct comparison to the earlier work.","section":"Section 5"},{"comment":"The issue of differing upstream electron thermal velocity between mass-ratio runs is currently confined to a footnote; given its potential impact on the high-MA results, it should be addressed in the main text with a clear statement of the confound.","section":"Section 2, footnote 3"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful methodological contribution and the low-MA instability-suppression argument is well grounded. The main reservation is the uncontrolled change in upstream electron thermal velocity between the mass-ratio runs, which directly affects the high-MA conclusions, and the over-strong 'independence' claim based on two data points. These are fixable with additional simulations or careful reformulation, so major revision rather than rejection seems appropriate. I also note the discrepancy between the Mach numbers quoted in the text and those in the figures, which the authors should correct during revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does a genuinely useful thing: it runs 1D3V PIC simulations of strong parallel shocks with mr=100 and mr=1836 at two Alfvén Mach numbers and shows that the mass ratio changes the sign of electron acceleration efficiency—suppressing it at low MA, over-producing it at high MA. That high-MA comparison is new to me, and the code (SHARP) is built to minimize numerical heating, which matters for realistic mass ratio runs. The low-MA electron suppression also has a concrete mechanistic anchor in the intermediate-scale instability threshold from their earlier work, so that part is on solid ground. I also credit the honest caveat that they simulate only the first ~1.5 minutes of shock evolution and can't speak to final maximum energies.\n\nThe soft spot, and it is load-bearing for the abstract's claim, is the confound the stress-test note flags. They fix kBTe in energy units while changing mr, so the upstream electron thermal velocity is sqrt(4e-8*mr) c—a factor 4.3 larger in the mr=1836 runs. The high-MA difference in electron acceleration and ion flux is attributed to ion-to-electron scale separation, but those two runs differ in both scale separation and initial electron thermal spread. Footnote 3 acknowledges the confound but doesn't control it. A run with mr=100 and a lower upstream electron temperature to match v_te of the mr=1836 run would settle this. Without that, the high-MA conclusion is suggestive, not decisive.\n\nTwo smaller issues. The 'temperature ratio independent of upstream magnetic field' claim rests on exactly two Mach numbers with no run-to-run variability, so it's a trend, not a demonstrated independence. And the non-thermal threshold (p > 5 pmax) is standard but arbitrary; it's fine as a diagnostic, just not something to over-interpret.\n\nOverall: the paper is a legitimate caution to the PIC shock community, and the central qualitative result—mass ratio changes electron heating and acceleration in ways that can flip sign—will likely hold up. But the mechanism attribution needs the controlled run and more Mach numbers or error bars.\n\nA serious referee should see this if it goes through revision; I'd send it to review but with the confound and the two-point independence claim flagged as required changes.","headline":"A useful warning that reduced mass ratio can flip electron acceleration efficiency in PIC shock runs, but the headline causal attribution to scale separation is not yet controlled because the upstream electron thermal speed changes with mr.","tokens_in":13953,"tokens_out":1087,"would_cite":true,"duration_ms":13728,"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":"Reduced mass ratio in PIC shock simulations produces qualitatively wrong particle acceleration and heating.","keywords":["particle-in-cell simulations","collisionless shocks","ion-to-electron mass ratio","electron acceleration","thermal energy dissipation","Alfvén Mach number","intermediate-scale instability","cosmic ray acceleration"],"falsifier":"Run a series of 1D3V particle-in-cell simulations at Alfvén Mach numbers 5.3 and 21.3 with a fixed upstream electron energy but varying the ion-to-electron mass ratio across 100, 400, 900, and 1836, and measure the downstream electron-to-ion temperature ratio and the non-thermal electron fraction; if these quantities change smoothly with mass ratio or match the realistic-ratio result already at intermediate ratios, the claim that only the full realistic ratio recovers the correct physics would be weakened. Additionally, perform a control run at mass ratio 100 with the electron thermal velocity artificially matched to that of the mass-ratio-1836 run; if the reduced-mass-ratio anomalies persist, the paper's attribution to ion-to-electron scale separation would be supported rather than the confound of different initial electron thermal spread.","tokens_in":12996,"feed_emoji":"⚡","tokens_out":3009,"duration_ms":31416,"temperature":0.7,"pith_summary":"This paper argues that the common cost-cutting trick of lowering the ion-to-electron mass ratio in particle-in-cell simulations of strong collisionless shocks does not merely introduce quantitative error: it changes the physics. Comparing runs with mass ratio 100 against runs with the realistic 1836 in one-dimensional shock simulations, the author finds that the reduced ratio completely suppresses electron acceleration at low Alfvén Mach number and, at high Alfvén Mach number, produces an unrealistically strong ion flux and over-efficient electron acceleration. These differences appear in both thermal energy partitioning and non-thermal acceleration efficiencies. The paper's central claim is that correctly resolving the ion-to-electron scale separation is necessary to capture how shock kinetic energy is dissipated into heating and particle acceleration.","feed_headline":"Shock simulations need the real proton-to-electron mass ratio","feed_subtitle":"Reducing the mass ratio to 100 flips electron heating and acceleration results at both low and high Alfvén Mach numbers.","key_machinery":"The key machinery is the comparison between two mass-ratio settings across two Alfvén Mach numbers, with the ion-to-electron scale separation expressed through the ratio of ion to electron skin depths and gyro-radii. The argument hinges on the intermediate-scale instability, a wave mode whose growth and resonant interaction with electrons is only captured when the mass ratio is large enough that the electron gyro-radius can match the unstable wavelength. The realistic-ratio runs resolve this instability, whereas the reduced-ratio runs either suppress it (low Mach number) or produce spurious resonant behavior (high Mach number), explaining the qualitative differences in electron acceleration and in the electron-to-ion temperature ratio.","core_discovery":"The paper demonstrates that in 1D3V particle-in-cell simulations of strong, parallel, non-relativistic electron-ion shocks, the choice of ion-to-electron mass ratio qualitatively determines the simulated dissipation channels. At Alfvén Mach number 5.3, the reduced mass ratio (100) suppresses electron acceleration entirely, while the realistic ratio (1836) permits efficient electron acceleration through the destabilization of intermediate-scale unstable wave modes. At Alfvén Mach number 21.3, the reduced mass ratio yields more efficient electron acceleration and a nonphysical high-momentum ion flux relative to the realistic ratio. In the thermal sector, the reduced ratio over-heats electrons and under-heats ions (at low Mach number), so the downstream electron-to-ion temperature ratio comes out high in both reduced-mass-ratio runs, whereas the realistic-ratio runs show this temperature ratio to be independent of the upstream magnetic field. The paper concludes that reduced-mass-ratio simulations misrepresent both heating and acceleration, and that realistic mass ratios are required to infer the physics of astrophysical collisionless shocks.","pith_inferences":["If the electron-to-ion temperature ratio is independent of upstream magnetic field in realistic mass ratio simulations, the dominant heating balance may be set by the sonic Mach number and the compression ratio rather than by the Alfvénic Mach number; this could be tested by varying the sonic Mach number while keeping the mass ratio fixed.","The mass-ratio sensitivity at high Alfvén Mach number suggests that a resonant scale (likely the electron gyroradius relative to the wavelength of the unstable gyroscale wave) controls the electron acceleration efficiency; running simulations at intermediate mass ratios (e.g., 400 or 900) would reveal whether the transition is continuous or abrupt.","The upstream electron thermal velocity is larger in the reduced-mass-ratio runs because the temperature is fixed in energy units; an explicit control run matching the electron thermal velocity between mass ratios would isolate whether the observed differences come from scale separation or from the initial electron thermal spread.","The conclusion that the ion-to-electron temperature ratio is independent of the upstream magnetic field, if confirmed in higher-dimensional simulations, could be used as a calibration target for sub-grid models of collisionless shock heating in astrophysical fluid simulations."],"forward_implications":["Simulations with reduced mass ratios cannot be used to infer electron acceleration efficiencies or electron-to-ion temperature ratios in strong non-relativistic parallel shocks.","The downstream electron-to-ion temperature ratio is claimed to be independent of the upstream magnetic field strength when a realistic mass ratio is used, providing a new observational constraint for shock models.","Published results relying on reduced mass ratios for non-relativistic shock microphysics should be revisited, especially for predicting the electron injection into diffusive shock acceleration.","The study shows that roughly 78% of upstream kinetic energy is converted to downstream thermal energy (mostly ions) in realistic-ratio runs, with about 11-12% going to non-thermal ions, regardless of Alfvén Mach number.","Since the simulated physical time is only about 1.5 minutes for supernova-remnant conditions, the maximum particle energies reached are not yet in the regime of long-lived shocks, so conclusions about final maximum energies must wait for longer runs."],"supporting_citations":[{"why":"Provides the previous demonstration that realistic mass ratios destabilize intermediate-scale instabilities and greatly increase electron acceleration efficiency at low Alfvén Mach number, the baseline this paper extends.","marker":"Shalaby et al. 2022"},{"why":"Introduces the intermediate-scale instability and its growth rates, which the paper uses to predict the wavelength and resonant particle momentum for electron acceleration.","marker":"Shalaby et al. 2021"},{"why":"Supplies the linear dispersion and polynomial equation for the intermediate-scale instability, used here to estimate the fastest-growing wave number and the expected resonant momentum.","marker":"Shalaby et al. 2023"},{"why":"Hybrid simulation result that about 10-12% of upstream kinetic energy goes into non-thermal ions, which the paper compares against its measured ion acceleration efficiency.","marker":"Caprioli & Spitkovsky 2014"},{"why":"Describes the SHARP particle-in-cell code and its higher-order interpolation that maintains exact momentum conservation and low numerical heating, enabling the realistic-mass-ratio runs.","marker":"Shalaby et al. 2017b"},{"why":"Provides the method for defining the non-thermal energy fraction as energy in particles with momentum above five times the peak of the distribution, which the paper adopts for its acceleration-efficiency calculation.","marker":"Xu et al. 2020"}],"fun_headline_variants":["Mass ratio choice flips shock electron acceleration results","Reduced ion-to-electron ratio skews shock heating and acceleration","Realistic mass ratio essential for accurate shock simulations","Shock sims: reduced mass ratio inverts electron acceleration"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The upstream electron temperature is set to the same energy value in the reduced and realistic mass-ratio runs, so the electrons in the reduced-mass-ratio simulation start with a larger thermal velocity; the paper attributes all downstream differences to the mass-ratio scale separation, but this uncontrolled difference in initial electron thermal spread could contribute to the observed effects.","fun_headline_variants_meta":{"raw":{"variants":["Mass ratio choice flips shock electron acceleration results","Reduced ion-to-electron ratio skews shock heating and acceleration","Realistic mass ratio essential for accurate shock simulations","Shock sims: reduced mass ratio inverts electron acceleration"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000552,"raw_usage":{"total_tokens":2713,"prompt_tokens":1104,"completion_tokens":1609,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":720,"completion_tokens_details":{"reasoning_tokens":1544}},"tokens_in":720,"tokens_out":1609,"duration_ms":13305,"temperature":1.0,"reasoning_tokens":1544,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:16:58.782820+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a series of 1D3V particle-in-cell simulations at Alfvén Mach numbers 5.3 and 21.3 with a fixed upstream electron energy but varying the ion-to-electron mass ratio across 100, 400, 900, and 1836, and measure the downstream electron-to-ion temperature ratio and the non-thermal electron fraction; if these quantities change smoothly with mass ratio or match the realistic-ratio result already at intermediate ratios, the claim that only the full realistic ratio recovers the correct physics would be weakened. Additionally, perform a control run at mass ratio 100 with the electron thermal velocity artificially matched to that of the mass-ratio-1836 run; if the reduced-mass-ratio anomalies persist, the paper's attribution to ion-to-electron scale separation would be supported rather than the confound of different initial electron thermal spread.","supporting_citations":[{"cited_title":"2022, ApJ, 932, 86","cited_arxiv_id":null,"evidence_quote":"Provides the previous demonstration that realistic mass ratios destabilize intermediate-scale instabilities and greatly increase electron acceleration efficiency at low Alfvén Mach number, the baseline this paper extends."},{"cited_title":"2021, ApJ, 908, 206","cited_arxiv_id":null,"evidence_quote":"Introduces the intermediate-scale instability and its growth rates, which the paper uses to predict the wavelength and resonant particle momentum for electron acceleration."},{"cited_title":"2023, Journal of Plasma Physics, 89, 175890603","cited_arxiv_id":null,"evidence_quote":"Supplies the linear dispersion and polynomial equation for the intermediate-scale instability, used here to estimate the fastest-growing wave number and the expected resonant momentum."}],"review_version":1}