{"id":"9fcd7be7-e26f-42e7-83b5-729984d95f2b","arxiv_id":"2507.08257","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The first IFU observations of RCW86's southwestern Balmer-dominated shock map shock velocities of 500-1500 km/s and show electron-ion equilibration decreasing with shock speed, plus a new intermediate H-alpha component.","lead":"The SDSS-V Local Volume Mapper observed the southwestern shock of the supernova remnant RCW86 and measured how shock speed and electron-to-proton temperature ratio vary across the remnant at roughly 0.3 pc scales. The data support and extend prior conclusions that faster collisionless shocks have lower electron-ion temperature equilibration, and they reveal a new intermediate-width H-alpha component likely produced by a neutral precursor.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Te/Tp trend is inferred with precursor-free vA08 grids, yet §5.5 detects a neutral precursor in the southern shocks anchoring the slow-shock end; precursor-driven Ib/In suppression is not excluded.","rationale":"The paper is an observational study with genuinely new IFU data, and the agreement of BDS-S1-smaller with the earlier Ghavamian et al. (2001) slit measurements is good independent support for the data quality and the broad-component fits. The central claim, however, depends on comparing those measurements to van Adelsberg et al. (2008) models that omit neutral and cosmic-ray precursors. The reader's weakest-assumption analysis already identified model accuracy as the main risk; I agree with that, and would sharpen it: the paper's own detection of an intermediate Hα component in the southern shocks means a precursor is demonstrably present in the very regions that define the slow-shock, high-Te/Tp side of the trend. If precursor effects suppress Ib/In at low velocities, the apparent move toward full equilibration could be an artifact of using a precursor-free model grid. This does not overturn the paper—the result is consistent with previous studies and the authors are transparent about the missing physics—but it does mean the headline inference should be presented as conditional on precursor-free shock models. The existing CONDITIONAL verdict is therefore appropriate; my concern adds a concrete acceptance condition rather than changing the verdict. I set agreement_with_reader to partial because the reader pointed to model accuracy in general, while the most specific load-bearing issue is the internal tension with the precursor detection in §5.5 and the velocity-dependent contamination correction in Eq. (10).","tokens_in":33339,"tokens_out":6434,"duration_ms":77891,"concrete_test":"Recompute the Figure 7 inference using precursor-inclusive shock models, e.g., the Morlino et al. (2012) grids, for vs = 500–1500 km/s and Te/Tp = 0.01–1, applied to the same FWHM and Ib/I'n values. If the inferred Te/Tp versus vs trend persists with similar slope, the concern is resolved; if the slow-shock points shift to lower Te/Tp or the trend flattens, the claim requires qualification. As a cross-check, refit the BDS-S1 and BDS-S1-smaller stacked spectra with a narrow+intermediate+broad model and recompute Ib/I'n while attributing the intermediate component to the precursor, to test whether the slow-shock Ib/In anchor survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central inference in §5.1.3 maps observed Ib/In and FWHM to (vs, Te/Tp) using van Adelsberg et al. (2008) grids, which the paper explicitly states exclude cosmic-ray and neutral precursors. The slow-shock endpoint of the claim is not purely hypothetical: §5.5 reports an intermediate Hα component (FWHM 192–207 km/s) in BDS-S1 and BDS-S1-smaller, attributed to a neutral return-flux precursor. Thus the stacked regions that anchor the near-full-equilibration interpretation at low shock velocities are, by the paper's own analysis, in a regime where the comparison model omits a physical ingredient that can alter Ib/In. Morlino et al. (2012) show precursor heating can lower Ib/In relative to precursor-free predictions; if that happens here, the apparent rise of Te/Tp toward low vs could be a precursor effect rather than a real electron-ion equilibration trend. The radiative-contamination correction in Eq. (10), with a fixed (NII/Hα)_rad = 2/3, is a second velocity-dependent bias because southern fibers lie nearest the bright radiative cloud. The by-eye definition of the BDS groups makes a quantitative sensitivity test especially important. This is a load-bearing model-dependence, not merely a disagreement with the current consensus.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents SDSS-V LVM integral-field spectroscopy of the southwestern part of the supernova remnant RCW86, covering the Balmer-dominated non-radiative shock and adjacent radiative clouds at ~0.3 pc scales. After median stacking and local background subtraction, the authors fit Hα+[NII] with one- and two-component Gaussian models and select 148 spaxels with a statistically significant broad Hα component, using AIC-based criteria validated in Appendix A. From the broad-to-narrow intensity ratio and broad FWHM, they derive shock velocities of ~500–1500 km/s and, by comparison with the van Adelsberg et al. (2008) model grids for Te/Tp = 0.1, 0.5, and 1.0, conclude that faster shocks have lower electron-ion temperature equilibration while slower shocks tend toward full equilibration. Additional results include broad-narrow centroid offsets, enhanced narrow Hα/Hβ ratios attributed to Lyβ trapping, a marginal HeII 4686 detection implying a high preshock neutral fraction, and the first intermediate Hα component in RCW86, interpreted as a neutral precursor. A brief comparison of radiative-cloud line ratios with MAPPINGS-III grids is also presented.","tokens_in":33539,"tokens_out":7576,"duration_ms":84062,"significance":"If the equilibration trend survives the model-dependence tests requested below, the paper delivers a significant observational step: a single SNR covering a wide velocity range, filling the Ib/In–FWHM plane previously populated by a handful of slit observations, and providing the first evidence for a neutral precursor in RCW86. The analysis is careful in several respects: the broad-component selection uses explicit AIC and FWHM-ratio criteria with rejected examples shown in Appendix A, the background-subtraction robustness is tested with three independent fiber sets in Appendix B, and the measurements are consistent with previous slit spectroscopy of the same region (Section 5.3). The main caveat is that the central Te/Tp inference relies on precursor-free models even though the data indicate a neutral precursor in the slow-shock region; this is acknowledged in Section 5.1.3 but not in the abstract or the conclusions.","major_comments":[{"comment":"The central inference that Te/Tp rises toward unity at low shock velocity is made by comparing measured Ib/In and FWHM with the van Adelsberg et al. (2008) model grid, which the authors explicitly state excludes cosmic-ray and neutral precursors. The same dataset, however, shows an intermediate Hα component with FWHM 192–207 km/s in BDS-S1 and BDS-S1-smaller (Section 5.5, Figure 12), attributed to a neutral return-flux precursor, and those southern stacks are precisely the slow-shock regions that anchor the near-full-equilibration endpoint. Because Morlino et al. (2012), cited in Section 5.1.3, show that precursor heating can lower Ib/In relative to precursor-free predictions, the apparent trend toward Te/Tp ~ 0.5–1 at low vs could be a precursor effect rather than genuine electron-ion equilibration. I ask the authors to compare the slow-shock stacks with precursor-inclusive models, at least in the 685–941 km/s range and extending Morlino et al.'s 10^3 km/s grid downward, or to restate the abstract and Section 7 conclusions as an upper limit with this caveat made explicit.","section":"§5.1.3, Figure 7; §5.5"},{"comment":"The radiative-shock contamination correction subtracts a fixed (NII/Hα)_rad = 2/3 of the [NII]6584 intensity from the narrow Hα flux in every fiber. The correction is therefore largest in the southern fibers, which have the highest [NII]/Hα ratios (Figure 5, top panel) because they lie nearest the bright radiative cloud; this differentially increases Ib/In' at low shock velocities and can mimic the rise of Te/Tp toward the slow-shock end. Please report the sensitivity of Figure 7 to (NII/Hα)_rad over a plausible range (e.g., 0.4–1.0) and to the by-eye membership of the BDS-S1 and BDS-S1-smaller stacks, and state explicitly whether the trend persists.","section":"§5.1.3, Eq. (10)"},{"comment":"The paper itself notes in footnote 23 that it cannot verify the Te/Tp-vs relation because it lacks the full model grid, and the comparison in Figure 7 is presented by eye, with the data split into panels according to the Te/Tp value used to convert FWHM to shock velocity. As a result, the statement that electrons and ions are not fully equilibrated in faster shocks but tend toward full equilibration in slower shocks is a qualitative consistency statement rather than a statistical inference. Please add a quantitative agreement measure (e.g., a chi-squared or bootstrap test over the three panels, including the covariance of FWHM and Ib/In and the uncertainties in the model grid) or clearly label the claim as illustrative.","section":"§5.1.3, Figure 7 and footnote 23"}],"minor_comments":[{"comment":"The abstract states shock velocities of 500–900 km/s in the south and 1000–1500 km/s in the north, while Section 4.2 reports broad-component FWHM values of 400–800 km/s and 1000–1500 km/s for the same regions; please clarify whether the quoted numbers are shock velocities or FWHM and make the wording consistent.","section":"Abstract and §4.2"},{"comment":"The BDS-S1-smaller group is defined as four fibers with the most prominent broad Hα components selected visually; please give the fiber count and selection criterion for every BDS group, since the stacked spectra in Figure 11 are the basis for the He and Balmer-decrement results.","section":"§5.3, paragraph defining BDS groups"},{"comment":"The phrase 'the broad component amplitude is forced to be negative' should read 'the fitted broad-component amplitude goes to zero or negative' to avoid implying a physical constraint rather than a fitting behavior.","section":"§5.5, penultimate paragraph"},{"comment":"The reference entry for van Adelsberg et al. (2008) is missing the journal name and article identifier; as printed it appears as a title and DOI only.","section":"References"},{"comment":"Figure 7 would benefit from a legend identifying the three Te/Tp panels and from error bars on a few representative data points; the grayscale overplot makes the comparison hard to evaluate by eye.","section":"Figure 7"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern is well-founded and lands: the slow-shock endpoint of the equilibration claim is anchored in the same southern stacks that show a neutral precursor, while the comparison model grid excludes precursors. This should be addressed by a sensitivity test with precursor-inclusive models or by softening the abstract and conclusions. The paper is otherwise well executed, with transparent selection criteria, robust background tests, and a valuable IFU dataset."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the genuinely new stuff: this gives the first ~0.3 pc IFU view of the southwestern RCW86 shock, 148 fibers with broad H-alpha, a nice velocity gradient from 500-900 km/s in the south to 1000-1500 km/s in the north, and the first intermediate H-alpha component in this remnant. The Balmer-decrement and HeII results are plausible additions. The data themselves are worth having.\n\nThe paper's central physical claim is that Te/Tp falls from near equilibrium at low shock speeds to ~0.1 above 800 km/s. That is consistent with prior work, so the confirmation is fine. But the quantitative inference at the slow end has a specific problem the authors walk right up to and then step around. They use the van Adelsberg et al. (2008) grids, which are precursor-free, and in Section 5.5 they report an intermediate H-alpha component in BDS-S1 and BDS-S1-smaller, the very regions that anchor the near-equilibration points. Morlino et al. (2012) show that neutral precursors can suppress Ib/In relative to precursor-free predictions, and in the vA08 mapping a suppressed ratio moves a point toward higher Te/Tp. So part of the apparent rise in Te/Tp at low vs could be a precursor effect, not genuine equilibration. The authors acknowledge the model limitation in 5.1.3, but the abstract and summary don't carry that caveat.\n\nOther soft spots are milder. The radiative contamination correction uses a fixed (NII/Halpha)_rad = 2/3, which is honest but clearly velocity-dependent; the BDS groups are selected by eye; and there is no mention of data or model-grid release. The background-subtraction test in Appendix B is good, and the agreement with Ghavamian et al. (2001) on BDS-S1-smaller is reassuring.\n\nWho is this for? People working on collisionless shocks, SNR kinematics, and wide-field IFU spectroscopy. It deserves a serious referee, as the dataset alone justifies that. I'd ask for a sensitivity analysis on the radiative correction, the model grids, and either a precursor-inclusive model comparison or a much more cautious statement about the low-velocity equilibration values.\n\nRecommendation: send it to review. The data are real, the analysis is mostly careful, and the main scientific caveat is identifiable and fixable with a softer claim or extra modeling.","headline":"The IFU dataset is a solid new resource, but the slow-shock end of the Te/Tp trend is on shakier ground than the abstract suggests, because the anchoring regions show the very neutral precursor that the comparison models omit.","tokens_in":34313,"tokens_out":3595,"would_cite":true,"duration_ms":38195,"reading_group":"yes","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 supernova remnant maps how collisionless shocks divide heat between electrons and protons.","keywords":["supernova remnants","Balmer-dominated shocks","collisionless shocks","electron-ion equilibration","shock precursors","integral field spectroscopy","RCW86","H-alpha line profiles"],"falsifier":"Recompute the same broad-to-narrow versus FWHM comparison using a model grid that includes neutral-precursor heating and full Lyman-line radiative transfer; if the fast shocks above 800 km/s become consistent with full equilibration once precursors are included, the claimed velocity-dependent equilibration would be an artifact of the simplified models.","tokens_in":1744,"feed_emoji":"🌌","tokens_out":3187,"duration_ms":79980,"temperature":0.7,"pith_summary":"This paper reports the first integral-field spectroscopic map of the southwestern shock front of the supernova remnant RCW86, resolving the Balmer-dominated filaments at roughly 0.3 parsec scales. It argues that the shock's broad and narrow hydrogen-alpha components reveal a clear pattern: faster shocks, above about 800 km/s, keep electrons much cooler than protons, with the temperature ratio approaching 0.1, while slower shocks, below about 500 km/s, approach full electron-proton equilibration. The paper also reports the first intermediate-velocity H-alpha component in RCW86, which it attributes to a neutral precursor ahead of the shock, plus elevated Balmer decrements consistent with Lyman-beta trapping. The significance is that a single remnant now populates the parameter space previously assembled from many remnants, sharpening observational constraints on how collisionless shocks distribute energy between particle species.","feed_headline":"RCW86 shocks show electrons stay cool where speeds top 800 km/s","feed_subtitle":"Integral-field spectra map 148 shock filaments in one remnant, tying electron-ion heating to shock speed.","key_machinery":"The central object is the Balmer-dominated shock, a fast non-radiative collisionless shock in partially ionized gas whose H-alpha line splits into a narrow component from cold pre-shock neutrals and a broad component from fast neutrals produced by charge exchange with hot post-shock protons. The argument is carried by the relation between the broad-to-narrow intensity ratio, the broad-line FWHM, and a model grid that converts these observables into a shock velocity and an electron-to-proton temperature ratio. The integral-field unit provides the new capability: 148 fibers across one shock front now fill the broad-to-narrow versus FWHM plane that earlier single-slit studies had sampled with only a few points from different remnants.","core_discovery":"On the paper's own terms, the discovery is that electron-ion temperature equilibration in collisionless shocks varies systematically with shock speed within one supernova remnant. Comparing the broad-to-narrow H-alpha intensity ratio and the broad-line FWHM against the adopted model grid that maps these observables to shock velocity and electron-to-proton temperature ratio, the authors find that the faster parts of the RCW86 shock, above about 800 km/s, favor a temperature ratio close to 0.1, while the slower parts, below about 500 km/s, favor a ratio close to 0.5 to 1. The broad-component FWHM spans roughly 400 to 800 km/s in the south and 1000 to 1500 km/s in the north, revealing a large velocity gradient along the same shock front. The paper also reports the first intermediate-velocity H-alpha component in RCW86, with FWHM of 192 to 207 km/s, attributed to a neutral-precursor-heated region, and narrow-component Balmer decrements of 3 to 5 that support Lyman-beta trapping in the cold neutrals.","pith_inferences":["The same integral-field method could be applied to other Balmer-dominated remnants, such as Tycho, SN1006, and Kepler, to test whether the velocity-dependent equilibration trend seen in RCW86 is universal or depends on local conditions.","The inferred high neutral fraction in the southern shock, combined with the neutral-precursor evidence, suggests RCW86's forward shock is advancing into partly neutral material, which may sharpen limits on the ionizing flux of its Type Ia progenitor.","If neutral-precursor heating is modeled more completely, part of the scatter in the broad-to-narrow ratio may be absorbed by the precursor, possibly tightening the inferred temperature-ratio relation and exposing residual variations due to magnetic field or density structure.","The ongoing full-remnant survey that this observation previews could map the entire 360-degree shock front of RCW86 and test whether the velocity-equilibration relation is azimuthally uniform or correlated with the surrounding density distribution."],"forward_implications":["If correct, electron-ion equilibration is not a constant but a function of shock velocity, with faster shocks leaving electrons far cooler than protons in a way that supports the idea that electron heating is nearly independent of velocity above roughly 400 km/s.","The single-remnant mapping means Balmer-dominated supernova remnants can calibrate collisionless shock physics over a wide velocity range without the systematic differences that plague cross-remnant comparisons.","The first detection of an intermediate H-alpha component in RCW86 strengthens neutral-precursor models and gives observers a new line-width diagnostic for precursor heating and cosmic-ray effects.","The narrow-component Balmer decrements of 3 to 5 imply that Lyman-line trapping must be included when interpreting Balmer-dominated shock spectra, and that the cold neutral population is the main site of that trapping.","The broad-to-narrow centroid offsets of up to about 100 km/s indicate bulk post-shock motion or non-Maxwellian ion distributions, which complicates but also enriches simple velocity diagnostics."],"supporting_citations":[{"why":"Supplies the model grid that converts broad-line FWHM and broad-to-narrow ratio into shock velocity and electron-to-proton temperature ratio, the central comparison of the paper.","marker":"M. van Adelsberg et al. (2008)"},{"why":"Introduces the theory of two-component Balmer emission from shocks in partially ionized gas that underlies the entire analysis.","marker":"R. A. Chevalier & J. C. Raymond (1978)"},{"why":"Establishes how Lyman-line trapping enhances narrow H-alpha relative to H-beta, the basis for the Balmer-decrement interpretation.","marker":"R. A. Chevalier et al. (1980)"},{"why":"Provides earlier spectroscopy of the same RCW86 southwest shock filaments, including Monte Carlo Lyman-line trapping, and is the direct comparison baseline for the stacked regions.","marker":"P. Ghavamian et al. (2001)"},{"why":"Previously inferred a Te/Tp proportional to inverse shock-velocity-squared relation from Balmer-dominated shocks, the trend this paper reproduces within one remnant.","marker":"P. Ghavamian et al. (2007)"},{"why":"Quantifies the excitation, ionization, and charge-transfer rates that set the broad-to-narrow ratio as a function of velocity and equilibration.","marker":"K. Heng & R. McCray (2007)"},{"why":"Provides the neutral return-flux precursor models used to interpret the newly found intermediate H-alpha component.","marker":"G. Morlino et al. (2012)"},{"why":"First reported the RCW86 Balmer filaments and estimated shock velocities of 500 to 900 km/s, the starting point this work expands.","marker":"K. S. Long & W. P. Blair (1990)"}],"fun_headline_variants":["RCW86 shock-speed map links speed to electron-ion heating","Fast RCW86 shocks keep electrons cooler than ions","Electron-ion equilibration varies along RCW86 shock front","Shock speed sets electron-ion heating in RCW86"],"cache_read_input_tokens":36224,"weakest_assumption_plain":"The load-bearing premise is that the adopted model grid accurately predicts the broad-to-narrow ratio and line width for the plasma conditions in RCW86; if its assumptions about Lyman-photon transport, precursor heating, and pre-shock ionization are biased, the inferred electron-to-proton temperature ratios would shift.","fun_headline_variants_meta":{"raw":{"variants":["RCW86 shock-speed map links speed to electron-ion heating","Fast RCW86 shocks keep electrons cooler than ions","Electron-ion equilibration varies along RCW86 shock front","Shock speed sets electron-ion heating in RCW86"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000681,"raw_usage":{"total_tokens":3207,"prompt_tokens":1170,"completion_tokens":2037,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":786,"completion_tokens_details":{"reasoning_tokens":1969}},"tokens_in":786,"tokens_out":2037,"duration_ms":18117,"temperature":1.0,"reasoning_tokens":1969,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:21:41.380040+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the same broad-to-narrow versus FWHM comparison using a model grid that includes neutral-precursor heating and full Lyman-line radiative transfer; if the fast shocks above 800 km/s become consistent with full equilibration once precursors are included, the claimed velocity-dependent equilibration would be an artifact of the simplified models.","supporting_citations":[],"review_version":1}