{"id":"96875f46-ceaa-443a-a985-67bbd3840e20","arxiv_id":"2608.00253","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Adaptive boundary feedback holds collisionless plasma shocks nearly stationary in hybrid simulations, shrinking simulation domains by about 1000x while retaining quasi-steady ion spectra and shock profiles.","lead":"A new algorithm lets hybrid plasma simulations pin a collisionless shock in place instead of chasing it across a huge grid, cutting the required domain by about three orders of magnitude. The method is demonstrated in 2D for five magnetic-field angles, yielding shock profiles the authors compare to interplanetary shock observations and angle-dependent ion acceleration spectra.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The adaptive outflow boundary reflects rather than absorbs outgoing ions at low local density, so it may artificially return energetic ions and inflate acceleration spectra; without a same-parameter moving-shock comparison, the reported efficiencies are not established.","rationale":"The reader's weakest assumption is that the AFORA feedback boundary is physically equivalent to an open downstream reservoir; I agree this is the load-bearing premise, and I sharpen it by identifying the specular-reflection branch as a concrete mechanism that could artificially return energetic ions and inflate the acceleration spectra. The paper provides convincing internal evidence that the method is numerically stationary and efficient: the shock front and particle count flatten after a transient (Fig. 1), spectra converge after the cutoff (Fig. 2), and the reported profiles are qualitatively consistent with IP shock observations. These are real supporting facts. However, quasi-stationarity of x_sh and N_p is enforced by the controller—not independently validated—so it cannot certify physical fidelity. The absence of any same-parameter moving-shock comparison or sensitivity study over alpha, Δt_sh, and K(θ_B) leaves a clear path by which the boundary condition, rather than shock physics, sets the accelerated-ion distributions. This is not an internal inconsistency, but it is an unresolved correctness risk that directly affects the central claim of converged spectra and efficiency rankings. The proposed test directly compares AFORA against a standard moving-shock simulation and against a modified boundary without specular reflection; it would settle whether the artifact is real. Until such a check is done, a conditional verdict is appropriate.","tokens_in":17309,"tokens_out":4110,"duration_ms":48582,"concrete_test":"Run a benchmark case (θB=40°, M_A=3, β_i=0.383, same resolution) in three configurations: (i) AFORA as described; (ii) a standard moving-shock hybrid simulation with a large enough domain to reach spectral convergence; (iii) AFORA with the specular-reflection branch disabled, instead using a purely absorbing outflow plus a Rankine-Hugoniot particle injection rate at the boundary to maintain downstream density. Compare time-averaged density and magnetic profiles, r_n, r_B, and downstream EDFs over E/E0 ∈ [1,20]. If (i) matches (ii)/(iii) within the few-percent scatter quoted in Table 1, the boundary artifact is ruled out; if the high-energy tail or efficiency ranking shifts, AFORA spectra are controller-dependent and not yet physically validated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim requires that the adaptive semi-transparent boundary in Sec. 3 (Eqs. 3–4) behaves as an open downstream reservoir. The implementation reflects every outgoing particle specularly when the local cell density is below the controller threshold n_d. In a physical downstream reservoir, particles that leave do not return; specular reflection at the outflow plane is a nonphysical source of returning ions. Since DSA and SDA efficiencies depend on ions returning upstream for repeated shock crossings, this boundary can inflate the high-energy tail by artificially recirculating ions that should have escaped. The feedback keeps N_p and x_sh stationary by construction, so the stationarity diagnostics in Fig. 1 and the spectral convergence in Fig. 2 cannot distinguish this artifact from true shock physics. No same-parameter comparison with a standard moving-shock simulation is presented, and the claimed agreement with IP shock observations (Sec. 5.2) is qualitative. Thus the paper's claim that AFORA yields converged shock dynamics and spectra is not yet established. The controller also introduces unstated parameters (alpha, Δt_sh, K(θ_B)) whose sensitivity is untested, and the release of data does not substitute for a physical-equivalence validation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces AFORA, an adaptive-frame algorithm for hybrid (kinetic-ion/fluid-electron) simulations of collisionless shocks. A shock is first generated in the downstream frame, then the simulation is Lorentz/ Galilean transformed to the shock frame, and a semi-transparent downstream boundary with a density threshold n_d is continually adjusted via Eq. (4) to keep the shock front stationary. The authors demonstrate 2D runs at five magnetic inclination angles using the event-driven HYPERS code, report quasi-steady shock structures and downstream ion spectra, argue that the results agree with IP shock observations, and conclude that the method reduces the domain size by roughly three orders of magnitude, making multi-dimensional and eventually 3D shock simulations practical.","tokens_in":17635,"tokens_out":7359,"duration_ms":80235,"significance":"If the central claim is valid, this is a potentially important methodological contribution: holding the shock in the simulation frame would remove a serious bottleneck in kinetic shock modeling and could enable 3D hybrid simulations and self-consistent seed-spectra inputs for SEP/Fokker-Planck codes. The paper also has concrete positive features: it tests five shock geometries, reports run diagnostics, makes reduced data and processing scripts available on Zenodo, and relies on an asynchronous, event-driven solver with published predecessors. However, the central claim is not yet established because the adaptive boundary condition may itself control the physics it is supposed to simulate, and the supporting comparisons are partly qualitative or post hoc.","major_comments":[{"comment":"The central load-bearing concern is the physical equivalence of the adaptive outflow boundary to an open downstream reservoir. When the local cell density is below n_d, outgoing particles are reflected specularly; when above, they are absorbed. In a physical downstream region, a high-energy particle crossing the outflow plane downstream of the shock does not return. Here, a right-going downstream ion can be reflected, re-enter the downstream region, and possibly recross the shock, artificially contributing to the DSA/SDA cycle. Because Eq. (4) deliberately adjusts n_d to cancel shock-front displacements, the stationarity shown in Fig. 1 and the spectral convergence shown in Fig. 2 cannot by themselves distinguish this artifact from true shock physics. The claim that AFORA yields converged spectra and acceleration efficiencies requires a same-parameter comparison with a moving-shock simul","section":"Sec. 3, Eqs. (3)-(4)"},{"comment":"The paper states \"excellent agreement\" with IP shock observations, but no observed profile is plotted, overlaid, or scored. The comparison is purely qualitative and based on visual similarity of wave packets, ramps, and downstream oscillations. Given that reproducing observed IP shock profiles is one of the paper's validation pillars, this comparison needs to be quantitative: overlay representative observed and simulated profiles, define an error metric (e.g., amplitude, wavelength, and decay of downstream oscillations), and state the physical parameters of the events used.","section":"Sec. 5.2 and Figs. 4-6"},{"comment":"The cutoff time Omega_ci t = 750 is chosen post hoc, after the shock front has supposedly settled. Figure 2(c) shows that the energy-averaged scatter depends on the start time of the averaging window, but the paper does not test whether the acceleration efficiencies in Table 1 or the spectral ordering in Fig. 3 are robust to moving the cutoff earlier or later. Since the discarded transient contains the largest front excursions (Fig. 1), and since the controller is active throughout, the claim of quasi-stationarity should be supported by a sensitivity scan over the cutoff time rather than a single fixed value.","section":"Sec. 5, Fig. 2 and Table 1"},{"comment":"The empirical factor K(theta_B) is used to modify the shock-frame transformation, with values 1.3, 1.2, 1.1, 1.0, 1.0 for theta_B = 0-80 degrees, and the controller gain alpha in Eq. (4) is stated only as \"selected to guarantee alpha|Delta x_sh| << 1\". No sensitivity study is presented for either parameter. More importantly, since V_sh = K V_sh0 enters Eq. (3), the effective Mach number in the continuity estimate changes with K, while Section 4 still defines M_A = V_u/V_A = 3. The manuscript should clarify how K affects the physical shock parameters (compression ratio, upstream flow speed) and demonstrate that the reported spectra are not sensitive to reasonable variations of K and alpha.","section":"Sec. 3 and Sec. 4"},{"comment":"The convergence claim is based on runs with 25 macro-particles per cell and a 400x400 grid. No particle-number or resolution convergence test is reported. The high-energy tail in Fig. 3 is sampled by at most 44 out of 111 shots at the curve endpoints, which is marginal for a claim of converged spectra. A short convergence study (e.g., 50 or 100 particles per cell, or a finer mesh) would substantially strengthen the paper.","section":"Sec. 4, simulation setup"}],"minor_comments":[{"comment":"Typo: \"shock fame\" should be \"shock frame\".","section":"Sec. 3, paragraph after Eq. (2)"},{"comment":"The symbol M_A is used both for the nominal Mach number M_A = V_u/V_A and for M_A = V_u/|V_sh| in Eq. (3). These definitions differ when V_sh is modified by K(theta_B). Please standardize notation.","section":"Eq. (3) and Sec. 4"},{"comment":"The fitted power-law slopes (-4.8, -5.9, -6.9) are quoted in the text but the fits are not shown in Fig. 3. State the fitting interval and the uncertainty of each slope.","section":"Sec. 5.3 and Fig. 3"},{"comment":"The sentence \"These spectra were found not to be sensitive to the length of their recording box in the x-direction, up to 0.3 times the length of the simulation domain\" is not accompanied by any figure or table. Either show the test or move this claim to a supplemental section.","section":"Sec. 5.3"},{"comment":"The comparison with Caprioli & Spitkovsky (2014) uses L = 10^5 d_i, but the physical parameters of that reference differ from those used here. The factor-of-1000 domain reduction should be stated with the caveat that the comparison is across different runs/parameters.","section":"Sec. 2.2 and Sec. 6"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern is real and is the main reason I cannot recommend acceptance at this stage. The AFORA idea is promising and the paper is well organized, but the key claim that the adaptive boundary is equivalent to an open downstream reservoir needs a direct test. A same-parameter moving-shock benchmark, or at minimum a boundary-reflection tracer analysis, would resolve the issue. If the authors can provide that, the paper would be a strong contribution to the shock-simulation methodology literature. There are no concerns about citation patterns or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThis one is worth a look. The authors propose a frame-co-moving algorithm (AFORA) that holds a hybrid shock stationary by adaptively adjusting a downstream particle reflection/absorption threshold. If it works as advertised, it could cut simulation domain sizes by orders of magnitude and make 2D/3D parameter scans feasible. That's a real contribution.\n\nWhat's new: the density-threshold feedback (Eq. 4) is a genuine departure from Bennett & Ellison's velocity-based semi-transparent boundary and from Simon et al.'s faux shock. They demonstrate it in 2D for five magnetic angles, show shock-front and particle-count stationarity, plus spectral convergence diagnostics. Data and scripts are on Zenodo. The paper is clearly written and does not hide the empirical K(theta_B) factor.\n\nThe soft spots are real. The outflow boundary reflects particles specularly when the boundary-cell density drops below n_d(t). That is not an open downstream reservoir; it is a controller that injects returning particles to hold density up. Because DSA relies on ions returning from downstream, this can inflate the high-energy tail. The stationarity diagnostics (x_sh, N_p) are exactly what the controller forces, so they do not independently validate the physics. The claimed \"excellent agreement\" with IP shock profiles is qualitative: no observed profiles are plotted or scored. K(theta_B) and alpha are free parameters whose sensitivity is untested. The spectra cutoff at Omega_ci t = 750 is reasonable but post hoc, though the scatter-error analysis helps.\n\nNone of this kills the idea. The controller might approximate a physical downstream reservoir well enough, especially since downstream turbulence isotropizes reflected particles. But that equivalence has to be shown, not assumed. A referee should ask for a matched moving-shock run (same M_A, theta_B, beta) and a quantitative profile comparison. If that passes, this could become a widely used tool.\n\nMy advice: send to peer review, conditional on those checks. It is a serious paper from a serious group; the burden is on the authors to prove the boundary does not manufacture spectra.","headline":"AFORA is a clever and potentially important shock-frame method, but until its reflective outflow boundary is benchmarked against a moving-shock run, the claimed spectra should be treated as provisional.","tokens_in":18094,"tokens_out":4383,"would_cite":true,"duration_ms":46395,"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 feedback-controlled outflow boundary holds simulated collisionless plasma shocks stationary in the simulation frame, allowing hybrid runs that match observations at roughly one-thousandth of the usual domain length.","keywords":["collisionless shocks","hybrid simulation","shock frame","adaptive boundary condition","ion acceleration","solar energetic particles","quasi-parallel shock","event-driven simulation"],"falsifier":"Run a controlled benchmark for the same physical parameters (M_A = 3, beta_i = 0.383, beta_e = 0.5, one magnetic inclination angle) in a conventional moving-shock domain long enough to reach a quasi-steady spectrum, then compare the downstream magnetic compression, wave spectra, and time-averaged ion energy distribution against the AFORA stationary-frame run. A systematic deviation beyond statistical scatter would mean the boundary controller is producing the reported spectra rather than the shock physics.","tokens_in":17227,"feed_emoji":"⚡","tokens_out":8205,"duration_ms":71598,"temperature":0.7,"pith_summary":"This paper aims to remove a major computational bottleneck in kinetic plasma physics: simulating collisionless plasma shocks, the kind that accelerate particles at coronal mass ejections and planetary bow shocks, currently requires large moving-shock domains and long run times. The authors introduce AFORA, an adaptive frame-of-reference algorithm that holds the shock front nearly stationary in the simulation box by continuously adjusting a semi-transparent downstream boundary. Using an event-driven hybrid code, they demonstrate in two dimensions that a small domain (a few hundred ion inertial lengths) yields quasi-steady shock structure and converged energetic-ion spectra over thousands of ion cyclotron periods, with magnetic profiles matching interplanetary shock observations. If the approach extends to three dimensions, it would make high-resolution 3D hybrid shock simulations and systematic parameter scans practical, and provide seed-ion distributions for solar energetic particle transport models.","feed_headline":"Adaptive boundary holds shocks still, shrinking run domains ~1,000x","feed_subtitle":"Stationary-frame runs match observed shock profiles and ion spectra at a fraction of the usual cost.","key_machinery":"The central mechanism is the AFORA feedback loop, built around a semi-transparent downstream boundary. Outgoing ions are reflected specularly if the local plasma density at their cell is below a threshold n_d, and absorbed otherwise. The threshold is initialized from the mass continuity equation, n_d = n_u(M_A + 1), and then every control interval Delta_t_sh it is updated by n_new = n_old exp(alpha Delta_x_sh), where Delta_x_sh is the measured displacement of the shock front from its target position. This converts the boundary into an adaptive open reservoir that cancels front drift by automatically changing the proportion of reflected vs. absorbed ions. The front position itself is tracked","core_discovery":"The paper's central claim is that a kinetic simulation of a collisionless shock can be run in the shock frame, with the shock front held artificially stationary, without losing the physical fidelity of a much larger moving-shock simulation. The AFORA algorithm first creates a back-propagating shock by reflecting the incoming plasma off a wall; after the shock speed is measured, the system is Galilean-transformed into the shock frame and the wall is switched to a controlled outflow. The controller reflects or absorbs outgoing ions depending on whether the local density is below or above an adaptive threshold n_d, seeded from the mass continuity equation and updated as n_new = n_old exp(alpha","pith_inferences":["A direct benchmark against a long, large-domain moving-shock run for identical physical parameters would isolate whether the boundary controller, rather than the shock physics, is responsible for the reported spectra.","The same density-feedback controller could in principle be adapted to full electromagnetic (full-PIC) models, potentially extending the domain reduction to electron-scale shock physics.","If the efficiency ranking holds in 3D, it makes a testable SEP prediction: quasi-parallel CME shocks should dominate the near-injection-energy seed population, while oblique shocks should contribute disproportionately to the highest-energy seed ions.","The exponential update rule behaves like a proportional controller on the downstream density; replacing it with a more physics-informed boundary that also controls momentum flux could improve robustness for very high Mach numbers or strong reflected-ion pressures."],"forward_implications":["Simulation domains and run times for converged shock dynamics can be reduced by nearly three orders of magnitude, making multi-dimensional parameter scans in Mach number, magnetic inclination, and beta computationally feasible.","Time-averaged, quasi-stationary ion spectra can be computed with a few-percent scatter, so acceleration efficiencies can be quantified reliably rather than inferred from transients.","Downstream spectra in the solar wind frame provide ready seed-ion inputs for Fokker-Planck transport models of solar energetic particles.","Resolved upstream wave packets, sharp ramps, overshoots, and downstream oscillations in the simulated magnetic profiles match observed interplanetary shock profiles for oblique and quasi-perpendicular geometries."],"fun_headline_variants":["Shock-frame fix cuts plasma sim cost ~1,000x","Pinning shocks still shrinks sim size, matches data","Adaptive frame holds shock, reduces run domain ~1,000x","Stationary shock frames slash simulation cost, keep accuracy","New algorithm keeps shock fixed, slashes sim size and time"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The method's physical faithfulness rests on the assumption that the semi-transparent outflow boundary — reflecting or absorbing ions based on an adaptively adjusted density threshold — behaves like a genuine open downstream reservoir for a supercritical shock, so that forcing the front to stand still does not distort the downstream turbulence, compression ratios, or the return of energetic ions to the shock front.","fun_headline_variants_meta":{"raw":{"variants":["Shock-frame fix cuts plasma sim cost ~1,000x","Pinning shocks still shrinks sim size, matches data","Adaptive frame holds shock, reduces run domain ~1,000x","Stationary shock frames slash simulation cost, keep accuracy","New algorithm keeps shock fixed, slashes sim size and time"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1254,"prompt_tokens":813,"completion_tokens":441,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":355}},"tokens_in":557,"tokens_out":441,"duration_ms":5863,"temperature":1.0,"reasoning_tokens":355,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T00:53:21.857063+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a controlled benchmark for the same physical parameters (M_A = 3, beta_i = 0.383, beta_e = 0.5, one magnetic inclination angle) in a conventional moving-shock domain long enough to reach a quasi-steady spectrum, then compare the downstream magnetic compression, wave spectra, and time-averaged ion energy distribution against the AFORA stationary-frame run. A systematic deviation beyond statistical scatter would mean the boundary controller is producing the reported spectra rather than the shock physics.","supporting_citations":[],"review_version":1}