{"id":"a1358f9c-975c-409e-85de-fffea873cd2e","arxiv_id":"2607.11426","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.5,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"HI+radio-continuum selection of cluster galaxies intercepts the early ram-pressure stage, yielding systems ~2.7\times more star-forming than field analogs of equal stellar and HI mass, while stacking shows progressive HI deficiency toward the cluster center.","lead":"MeerKAT HI and 1.4 GHz continuum maps of three nearby clusters catch 61 galaxies still gas-rich but already more star-forming than field counterparts. The dual selection isolates the brief early ram-pressure phase before gas is stripped, when compression boosts both star formation and radio luminosity.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The central claim rests on a clean observational selection and direct comparison to an identically cut xGASS control sample. The reader correctly flags the radio-excess interpretation as the softest ancillary point, but that interpretation is not required for the 2.7\times SFR excess or for the stacking-derived HI-deficiency gradient. Both of those results survive the usual observational caveats (projection, resolution, modest thermal contamination). No load-bearing flaw that would move the verdict away from ACCEPT was found.","tokens_in":26960,"tokens_out":414,"duration_ms":4759,"concrete_test":"Recompute the mean SFR ratio after (i) excluding the three GASP jellyfish galaxies that dominate the high-SFR tail and (ii) applying a uniform 6 % thermal subtraction to every continuum detection; if the factor remains >2 within the quoted uncertainties the dual-selection claim is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The dual-selection claim (HI + 1.4 GHz continuum) that the 61-galaxy sample is ~2.7\times more star-forming than mass-matched field galaxies is supported by the cumulative distributions (Fig. 3), KS tests, and the multi-panel comparisons to field scaling relations (Figs. 4–5, 8–9). The reader’s residual concern about the origin of any radio excess (old electrons vs. residual AGN / thermal free-free / SFR-timescale mismatch) is real but secondary: the paper already states thermal free-free is expected <6 % and that AGN templates were omitted from the SED fits, yet the primary 2.7\times SFR excess is measured from the SED-derived SFRs themselves, not from the radio luminosities. Projection effects, resolution differences between HI and continuum, and the modest sample size are acknowledged and do not reverse the ordering of the distributions or the stacking gradient. No internal inconsistency or hidden circularity is present.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript presents MeerKAT L-band observations of three local clusters (IIZW108, A4059, A3558) that simultaneously map HI 21 cm and 1.4 GHz continuum emission out to ~2 R200. From 116 HI detections the authors assemble a clean sample of 61 galaxies detected in both tracers, derive stellar masses and 100-Myr-averaged SFRs via Bagpipes SED fitting of UV-to-IR photometry, and classify objects by phase-space location and optical morphology (jellyfish/unwinding). They report that these dual-selected cluster galaxies are on average ~2.7 times more star-forming than xGASS field galaxies of comparable stellar and HI mass, interpret the excess as the signature of early ram-pressure compression before substantial gas removal, and use the LR/LSFR_R ratio together with offsets from the main sequence and HI–M* relation to sketch successive early-processing stages. Unwinding candidates are found to be systematically more massive and HI-rich than jellyfish systems. Stacking of 677 optically selected members further shows that the average HI mass is already ~0.4–0.5 dex below field values at R200 and declines further inside R500, reaching a factor ~3 deficit near the centre.","tokens_in":27209,"tokens_out":1338,"duration_ms":18337,"significance":"The work exploits MeerKAT’s simultaneous HI+continuum capability to isolate a rare, short-lived phase of environmental processing that optical or single-tracer radio surveys miss. The dual-selection result, the quantitative stacking gradient against mass-matched xGASS controls, and the tentative UW-to-JF evolutionary ordering are concrete, falsifiable contributions that will be directly testable with SKA-Mid. The analysis is transparent (KS tests, S/N≥5 stacks, explicit detection limits) and the caveats (solar RFI loss, ~40-arcsec resolution, projection effects) are stated. If the early-compression interpretation holds, the paper supplies a practical selection method for next-generation multi-phase ISM studies of cluster galaxies.","major_comments":[{"comment":"Section 3 and Figure 3: the headline factor “~2.7 times more star-forming” is central to the abstract and discussion, yet the precise statistical definition is not given. Is it the ratio of median SFRs after matching both M* and MHI, the ratio of means, or a KS-derived shift? Please state the matching procedure, the uncertainty on the factor, and whether the result survives when the xGASS comparison sample is further restricted to the same radio-luminosity floor used for the cluster sample.","section":null},{"comment":"Section 4.1 and Figures 8–9: the five-stage evolutionary sequence rests on the interpretation that LR/LSFR_R > 1 at Δlog10 SFR < 0 is produced by old relativistic electrons. The manuscript already notes that thermal free-free is expected <6 % and that AGN templates were omitted from the SED fits, but does not quantify how residual low-level AGN or a systematic mismatch between the 100-Myr SED SFR and the ~10-Myr radio tracer would move objects across the stage boundaries. A short quantitative estimate (or an explicit statement that the primary 2.7 factor is independent of the radio excess) is needed to keep the stage scheme from over-interpreting the data.","section":null},{"comment":"Section 4.2: the claim that unwinding morphology “may result from their large baryonic mass” is plausible but currently rests on a small subsample (10 UW vs 13 JF) whose stellar-mass distributions barely overlap. Given the acknowledged orientation bias (UW face-on, JF edge-on) and the modest numbers, the evolutionary sequence UW → JF should be presented more cautiously as a working hypothesis rather than a demonstrated time sequence, or supported by a simple survival-analysis or mass-matched comparison.","section":null}],"minor_comments":[{"comment":"Abstract and Section 5: “can features a radio luminosity excess” → “can feature”.","section":null},{"comment":"Table 1: the HI band and ΔV rows are useful; adding the final continuum RMS and HI column-density limits already quoted in the text would make the table self-contained.","section":null},{"comment":"Figure 1 / A.1: the 3–24σ contour levels are clear, but a scale bar in kpc (in addition to the 6-arcmin cutout size) would help readers unfamiliar with the redshift.","section":null},{"comment":"Section 2.1: the ~30–60 % data loss to solar RFI is stated; a brief note on whether the lost data preferentially affect one cluster or one side of the primary beam would reassure readers that the sample is not spatially biased.","section":null},{"comment":"Section 2.3 / Table 2: the xGASS matching is performed by drawing random realisations that match the stellar-mass distribution; stating the number of realisations (already 10) and the resulting scatter on MxGASS_HI in the table caption would improve reproducibility.","section":null},{"comment":"Figure 6 and Table 4: AHI and ARC are computed relative to their own barycentres and at different resolutions; a one-sentence reminder of these caveats in the caption would prevent over-interpretation of the modest ARC difference between UW and JF.","section":null},{"comment":"References: a few recent MeerKAT/ASKAP HI-deficiency papers (e.g., Reynolds et al. 2022 is cited; check completeness of 2023–2025 works) could be added for context, but this is optional.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is solid and well within A&A scope. The dual-selection result is the main novelty; the stacking gradient is confirmatory but cleanly executed. I see no circularity or hidden free parameters. The three major points are clarifications rather than new analyses and should be addressable in a short revision. I would be comfortable accepting after those points are tightened."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The headline result is real: dual HI + 1.4 GHz continuum selection on three local clusters yields 61 galaxies that sit ~2.7\times above the field SFR at fixed M* and MHI. That is the cleanest observational filter I have seen for the short-lived compression stage before significant gas removal. The stacking of 677 members then gives a clean radial HI-deficiency profile that reaches ~3\times poorer than mass-matched xGASS inside R500. Both numbers are new and useful for simulators and for planning follow-up.\n\nWhat they did well is straightforward. Sample cuts are explicit, SED fitting with Bagpipes is standard and transparent, phase-space bins follow Rhee et al., and the xGASS matching is done properly (repeated random draws). KS tests and S/N ≥ 5 stacks are reported without drama. The UW-versus-JF mass and asymmetry differences are modest but honestly presented; they do not over-claim an evolutionary sequence. Self-citations are to prior GASP and MeerKAT work that is actually used, not decorative.\n\nSoft spots exist but stay secondary. Solar RFI cost them one A3558 track and ~30–60 % of the data; resolution is ~40 arcsec so morphology is global only; thermal free-free is assumed <6 % and AGN templates are omitted from the SEDs. Those caveats matter for the radio-excess staging sequence (old electrons versus residual AGN or timescale mismatch), but they do not touch the primary 2.7\times SFR excess, which is measured from the SED SFRs themselves. Projection effects and modest sample size are acknowledged and do not reverse the ordering of the distributions.\n\nThis is for people who work on environmental quenching, multi-phase ISM, or SKA-pathfinder survey design. It does not rewrite the basic RPS picture, but it gives a concrete selection tool and a quantitative gradient that were missing. I would send it to a serious referee without hesitation; the data and comparisons are solid enough to survive normal observational scrutiny. Worth reading and, for anyone calibrating RPS models or planning HI+continuum follow-up, worth citing.","headline":"Solid MeerKAT dual-selection paper that cleanly isolates the early RPS compression phase and delivers a quantitative HI-deficiency gradient; main 2.7\times SFR excess is robust, radio-excess interpretation secondary.","tokens_in":27883,"tokens_out":539,"would_cite":true,"duration_ms":7414,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Dual HI and radio selection catches cluster galaxies while external pressure still boosts star formation before gas is stripped.","keywords":["galaxy clusters","ram-pressure stripping","HI deficiency","radio continuum","star formation","jellyfish galaxies","unwinding galaxies","MeerKAT"],"falsifier":"Deep multi-frequency radio continuum spectra of the same sample showing flat rather than steep spectral indices would favor ongoing shock acceleration or compression over residual old electrons, falsifying the claimed origin of the radio excess.","tokens_in":27893,"feed_emoji":"🌌","tokens_out":925,"duration_ms":11194,"temperature":0.7,"pith_summary":"Galaxies that fall into clusters lose gas and eventually stop forming stars, but the earliest moments of that process have been hard to catch. This paper uses MeerKAT observations of three nearby clusters to select 61 galaxies that still retain detectable neutral hydrogen and radio continuum emission. Those galaxies form stars roughly 2.7 times faster than field galaxies of the same stellar and HI mass, showing that the dual selection has intercepted the brief compression phase when external pressure has not yet removed the gas but has already raised the star-formation rate and radio luminosity. Unwinding spiral morphologies appear mainly among the most massive, gas-rich systems, while stacking of hundreds of optically selected members shows that average HI content falls by a factor of about three toward the cluster center. The result supplies a practical way to isolate the start of environmental quenching and to time the sequence of morphological and radio changes that follow.","feed_headline":"Cluster galaxies form stars 2.7\times faster before gas is stripped","feed_subtitle":"Dual HI and radio selection catches the brief compression phase that boosts star formation","key_machinery":"The joint HI 21 cm plus 1.4 GHz continuum selection, combined with the luminosity ratio L_R / L_R^SFR and the offset from the star-forming main sequence. Together these quantities order galaxies into successive early stages of ram-pressure compression, SFR burst, and subsequent decline.","core_discovery":"Cluster galaxies detected simultaneously in HI and 1.4 GHz continuum are on average ~2.7 times more star-forming than field galaxies of matched stellar and HI mass. The dual selection therefore intercepts the earliest stages of environmental processing, when external pressure compresses the interstellar medium and temporarily boosts both star formation and radio luminosity before significant gas removal occurs. Galaxies at this stage can also show radio excess relative to their current SFR, attributed to old relativistic electrons still permeating the disk.","pith_inferences":["If the dual selection cleanly isolates the compression phase, the same cut applied to SKA-Mid Band 2 surveys should yield statistical samples large enough to measure the duration of the SFR boost as a function of cluster mass and orbital parameters.","The higher radio asymmetry of jellyfish relative to unwinding systems implies that non-thermal plasma is stripped later than neutral gas, offering a timing diagnostic for multi-phase ISM-ICM coupling.","The radio excess attributed to old electrons predicts steep spectral indices at lower frequencies; a non-detection of that steepening would reopen the possibility of low-level AGN or free-free contamination."],"forward_implications":["HI+radio continuum selection becomes a practical filter for assembling samples of galaxies still in the pre-stripping compression phase.","The L_R / L_R^SFR ratio can separate galaxies observed before, during, and immediately after the compression-induced SFR boost.","Unwinding morphologies preferentially mark massive, baryon-dominated systems whose spiral structure resists rapid disruption.","Stacking confirms that average HI deficiency already exists near R_200 and grows by a factor of ~3 inside R_500, quantifying the radial progress of environmental gas loss."],"fun_headline_variants":["HI+radio galaxies form stars 2.7x faster before cluster gas strip","Dual HI continuum selection intercepts early ISM compression SF boost","Cluster HI+radio galaxies show 2.7 times higher star formation early","MeerKAT spots temporary 2.7x SF and radio excess before stripping","Early pressure phase leaves cluster galaxies 2.7x more star-forming"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"That the radio continuum excess relative to the SED-derived star-formation rate is produced by old relativistic electrons (or mild compression) rather than residual AGN, uncorrected thermal free-free emission, or systematic mismatch between the 100-Myr SED average and the ~10-Myr radio tracer.","fun_headline_variants_meta":{"raw":{"variants":["HI+radio galaxies form stars 2.7x faster before cluster gas strip","Dual HI continuum selection intercepts early ISM compression SF boost","Cluster HI+radio galaxies show 2.7 times higher star formation early","MeerKAT spots temporary 2.7x SF and radio excess before stripping","Early pressure phase leaves cluster galaxies 2.7x more star-forming"]},"model":"grok-4.5","effort":"low","cost_usd":0.008408,"raw_usage":{"total_tokens":2092,"prompt_tokens":940,"num_sources_used":0,"completion_tokens":83,"cost_in_usd_ticks":84080000,"prompt_tokens_details":{"text_tokens":940,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1069,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":940,"tokens_out":83,"duration_ms":8543,"temperature":1.0,"reasoning_tokens":1069,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T05:39:45.900741+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Deep multi-frequency radio continuum spectra of the same sample showing flat rather than steep spectral indices would favor ongoing shock acceleration or compression over residual old electrons, falsifying the claimed origin of the radio excess.","supporting_citations":[],"review_version":1}