{"id":"d1218ba0-492d-400c-ab16-2acea0c3f3fb","arxiv_id":"2504.16873","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"Applying LOFAR's reconstruction method to full SKA-Low simulations yields 5 to 8 g/cm2 precision on the air-shower maximum between 10^16.6 and 10^18 eV, with beamforming extending the range down to 10^16 eV.","lead":"This paper simulates how the future SKA-Low radio telescope would detect cosmic-ray air showers and estimates how precisely it could measure the shower maximum, a key proxy for cosmic-ray composition. It reports that SKA-Low could beat current LOFAR precision by roughly a factor of three and reach lower energies, which matters for mapping where galactic cosmic rays come from.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Beamforming emulation assumes sqrt(N) SNR gain for a sky-noise-dominated array; correlated Galactic noise at SKA-Low baselines likely reduces this, leaving the 10^16 eV threshold claim unproven.","rationale":"The reader identified the unmodeled full system response as the weakest assumption. That is a real and honestly disclosed limitation, and the paper's 'lower limit' framing partly insulates it. My stress-test pass found a more specific and less acknowledged assumption: the beamforming SNR gain. The paper repeatedly emphasizes that SKA-Low is sky-noise dominated, yet the beamforming emulation in Sections II D and III B boosts SNR by sqrt(N), which presupposes uncorrelated antenna noise. For an array where the dominant noise is diffuse Galactic emission, the spatial coherence of that noise over the 1-3 m station baselines means that beamforming will not deliver the full sqrt(N) improvement, especially at the lower end of the 50-350 MHz band. This directly undermines the claim that beamforming lowers the energy threshold to 10^16 eV without losing resolution, which is one of the two headline results (the other being the 5-8 g/cm2 Xmax precision). The concern is concrete and testable: one can generate noise traces with the correct spatial covariance from the same sky model and perform actual patch beamforming on a few simulated showers. Because the paper presents a planning baseline for SKA-Low, an overestimated beamforming gain would mislead system-design decisions and event-rate estimates. The high-energy Xmax precision result, including the dynamic-range caveat, is more robust to this concern, so the paper still merits conditional acceptance as a baseline study, provided the beamforming emulation is validated against correlated sky noise. Hence the reader's CONDITIONAL verdict is unchanged, but the required condition should include this beamforming check.","tokens_in":16655,"tokens_out":10604,"duration_ms":110878,"concrete_test":"Generate spatially correlated sky-noise voltage traces for 4- and 16-antenna patches using the same GSM2016 sky model, SKALA4 antenna response, and the true station geometry, so that the noise covariance between antennas reflects the diffuse Galactic background. Add these noise realizations to simulated CoREAS traces for a 10^16.5 eV shower, time-shift and sum each patch toward the known shower direction, and measure the realized SNR gain against a single antenna. Compare with sqrt(4) and sqrt(16); if the realized gain is below, say, 80% of sqrt(N), re-derive the energy-threshold and Xmax precision curves of Fig. 6 using the measured gains for N = 4, 16, and 64. If the low-energy reconstruction threshold shifts above 10^16 eV, the Summary's beamforming claim must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II D and Section III B emulate patch beamforming by using 1 in N antennas and increasing the SNR by sqrt(N). This scaling is correct only if the noise in the combined antennas is uncorrelated. The paper states that the SKA system is sky-noise dominated (Section II A, Fig. 2). At SKA-Low baselines of 1-3 m and frequencies 50-350 MHz, the Galactic background noise is partially spatially correlated within a station: at 50 MHz (lambda = 6 m) closely spaced antennas see nearly the same noise waveform, and at 150 MHz correlation over ~2 m baselines is still non-negligible. In a real beam sum, the signal amplitude adds as N, but the correlated fraction of the sky noise also adds coherently, so the SNR gain lies between 1 and sqrt(N), not at sqrt(N). The emulation therefore overestimates the beamforming benefit, and the claimed low-energy threshold of 10^16 eV (Fig. 6, Summary) is optimistic. A secondary effect is that time-shifting and summing finite-bandwidth traces with slightly different wavefront arrival times broadens the pulse, which the simple SNR-boost of a single antenna does not reproduce. The dynamic-range boundary and the statement that beamforming introduces no significant resolution loss near threshold inherit the same issue. The high-energy Xmax precision claim without beamforming is less affected, but the low-energy half of the central claim is not supported by the current simulation.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents end-to-end simulations of radio signals from cosmic-ray air showers as they would be recorded by the SKA-Low array. Showers are simulated with CORSIKA/CoREAS at a single primary energy of 10^17 eV and scaled to other energies; antenna voltages are produced with the SKALA4 model, and Galactic noise from the GSM2016 model is added, with a unit-gain de-dispersion filter standing in for the not-yet-available full system response. The LOFAR fluence-footprint reconstruction is applied in a leave-one-out mode to an ensemble of 140 simulated showers. The authors report an Xmax precision of 5 to 8 g/cm2 with negligible bias in the range 10^16.6 to 10^18 eV, a lowering of the detection threshold to 10^16 eV via patch beamforming emulated by downsampling and a sqrt(N) SNR boost, an analysis of dynamic-range restrictions showing degradation above about 10^17.5 eV for a 32-sigma clip limit, and a composition-statistics study indicating that 1000 to 3000 showers per energy bin would be sufficient to become systematics-limited.","tokens_in":16924,"tokens_out":8815,"duration_ms":90067,"significance":"If the quoted numbers hold, SKA-Low would offer an Xmax resolution roughly 2.5 to 4 times better than LOFAR across the knee-ankle range and would reach lower primary energies than current radio arrays, making it a potentially powerful instrument for composition measurements in an energy region relevant to the Galactic-to-extragalactic transition. The paper is valuable as a baseline study: it adapts a validated reconstruction pipeline to SKA, uses the official SKALA4 antenna model and a high-precision interpolation method, and is transparent about several limitations, including the absence of a full system response model, the ad hoc treatment of near-zero fluence uncertainties, and the arbitrary choice of beamforming threshold energies. The leave-one-out validation is an appropriate internal check, and the event-rate and data-volume estimates are useful for instrument planning. The main caveat is that the headline statements are conditional on the simulation model, and the beamforming emulation in particular rests on an SNR scaling that may be optimistic.","major_comments":[{"comment":"The beamforming emulation assumes that the noise in different antennas is uncorrelated, so that combining n antennas raises the SNR by a factor sqrt(n). This assumption is not justified for SKA-Low station-scale baselines in the 50-350 MHz band. At 50 MHz the wavelength is about 6 m, so baselines of a few meters are sub-wavelength, and the diffuse Galactic emission that dominates the system temperature has substantial large-scale structure; the sky-noise waveforms in closely spaced antennas are therefore partially correlated. In a real sum the SNR gain lies between 1 and sqrt(n), and the emulation by downsampling plus SNR boost overestimates the benefit. Since the 10^16 eV threshold claim in Fig. 6 and the Summary depends directly on this emulation, the low-energy half of the central claim is not established by the present simulations. I recommend simulating actual beam sums over realistic station layouts with a sky-noise model, or at minimum presenting the results as a band that brackets the correlated-noise case.","section":"Sections II D and III B"},{"comment":"The paper states that 'there is currently no model for the full system response' and therefore uses a unit-gain de-dispersion filter after the SKALA4 antenna model. This means that all results, including the fluence calibration, the dynamic-range boundary, and the beamforming threshold, are computed with a partial instrument model. A different phase or amplitude response in the real signal chain would change the pulse shapes, the fluence estimates, and the Xmax fits, and would also affect the SNR-based energy threshold. The precision and bias numbers should therefore be stated as explicitly conditional on this assumed response, and the sensitivity to plausible response variations should be quantified or at least discussed as a dominant systematic uncertainty.","section":"Section II A"}],"minor_comments":[{"comment":"The sentence 'The SNR would rise by a factor sqrt(n) when combining n antennas' should explicitly state the assumption of uncorrelated noise; as written it is presented as a general result rather than as an ideal-case approximation.","section":"Section II D"},{"comment":"The factor 5 applied to the zero-fluence uncertainty is an ad hoc choice. I recommend showing that the Xmax results are robust to this factor, for example by repeating the reconstruction with a different value, since the text notes that the weakest measurements can otherwise bias the footprint fit.","section":"Section II B"},{"comment":"The precision curves are presented without uncertainty bars. With about 130 showers after the 3.5% tail cuts, the statistical uncertainty on a quoted precision value is of order 10%, so small differences between the three zenith-angle curves should not be overinterpreted.","section":"Section III A and Fig. 5"},{"comment":"The data-volume estimate of 4.5 GB per shower should state explicitly that it includes two polarizations; the formula as written evaluates to 2.4 GB per polarization, and this distinction would help reproducibility.","section":"Section IV B and Eq. (6)"},{"comment":"The vertical threshold lines are described as 'chosen somewhat arbitrarily'; the exact threshold energies used for the different beamforming group sizes should be listed in the text or caption so that the figure is self-contained.","section":"Fig. 6 caption and Section III B"}],"recommendation":"major_revision","confidential_remarks":"The beamforming concern is the main technical obstacle to acceptance. If the authors can replace the downsampling-plus-SNR-boost emulation with a realistic beamformer simulation or provide a convincing bound on the correlated-noise effect, I would support publication. The paper fits the journal's scope and is a useful baseline study for a future instrument."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first complete end-to-end simulation study of SKA-Low as a cosmic-ray air shower array, applying the established LOFAR Xmax reconstruction to CoREAS simulations with the SKALA4 antenna model, Galactic noise, interpolation to the full core, beamforming emulation, dynamic range effects, and composition bootstrap. As a planning baseline it is genuinely useful, and the authors are transparent about most of their simplifications.\n\nThe high-energy half of the central claim survives scrutiny. The leave-one-out reconstruction within a 140-shower ensemble is internally consistent, and the 5–8 g/cm2 precision between roughly 10^16.6 and 10^18 eV is credible: it matches what you would expect from a very dense array with LOFAR-proven methods. The bias check, the SKALA2/SKALA4 comparison, and the dynamic range study are all honest and informative. I believe the paper delivers a reliable lower-limit estimate for the non-beamformed case.\n\nThe low-energy half is weaker. The beamforming emulation takes 1 in N antennas and boosts SNR by sqrt(N), which assumes the noise is uncorrelated between antennas in the summed group. But the text states the system is sky-noise dominated, and at 50–350 MHz with baselines of a few meters the Galactic noise is partially correlated across a station. So the real SNR gain from beamforming is between 1 and sqrt(N), and the 10^16 eV threshold in Fig. 6 is probably optimistic. The stress-test note is right about this. Also, time-shifting and summing finite-bandwidth traces broadens the pulse; a simple SNR boost does not capture that. This does not kill the paper, but it does mean the summary's beamforming claim should be read as a best-case scenario, not a baseline.\n\nOther soft spots are minor and mostly acknowledged: the energy is scaled from a single 10^17 eV simulation, the factor 5 on near-zero fluence uncertainty is ad hoc, and the removal of the 5 highest and lowest Xmax showers could flatter the precision. These are acceptable for a lower-limit study, though they mean the exact numbers should not be quoted without error bars.\n\nWho is this for? People planning SKA-Low's cosmic-ray mode, and anyone comparing radio Xmax techniques. It deserves a serious referee: I would send it to review, with the request that the beamforming simulation be improved or the threshold claim softened. I'd cite it for the high-energy precision baseline.","headline":"A solid, honest simulation baseline for SKA-Low cosmic-ray Xmax, where the high-energy precision claim holds up but the beamforming-enabled low-energy threshold rests on an optimistic SNR scaling.","tokens_in":17591,"tokens_out":2190,"would_cite":true,"duration_ms":23870,"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":"Simulations show SKA-Low measures cosmic-ray shower depth to 5-8 g/cm².","keywords":["cosmic rays","air showers","SKA-Low","radio detection","shower maximum","mass composition","beamforming","LOFAR"],"falsifier":"Once the real SKA-Low signal chain is characterized, measure the end-to-end amplitude and phase response of a SKALA4 element and receiver across 50 to 350 MHz and feed those measurements into the same simulation and reconstruction code; if the Xmax precision degrades beyond about 10 g/cm², the bias exceeds a few g/cm², or the beamformed energy threshold rises above $10^{16}$ eV, the paper's central performance claims would be contradicted.","tokens_in":16410,"feed_emoji":"📡","tokens_out":5446,"duration_ms":46655,"temperature":0.7,"pith_summary":"The paper argues that SKA-Low, the low-frequency Square Kilometre Array core with roughly 57,000 antennas packed into about one square kilometer, can reconstruct the depth of shower maximum (Xmax), the main mass-composition observable for cosmic-ray air showers, with a precision of 5 to 8 g/cm² between $10^{16}$.6 and $10^{18}$ eV and essentially no bias. This is roughly 2.5 to 4 times sharper than LOFAR, the current most precise radio air-shower instrument. Using the existing LOFAR reconstruction pipeline on full simulations, the authors show that patch-wise beamforming of groups of antennas extends the energy range down to $10^{16}$ eV, and that even with an 8-bit dynamic range restriction the resolution stays better than LOFAR's typical 20 g/cm² up to about $10^{17}$.5 eV. They further show that with 1000 to 3000 measured showers per energy bin, an Xmax-only mass composition analysis becomes limited by systematics rather than statistics, which motivates reconstructing the full longitudinal shower profile instead of just its maximum.","feed_headline":"Simulations show SKA-Low measures shower depth to 5-8 g/cm²","feed_subtitle":"Dense radio array beats LOFAR's precision by 2.5-4x and reaches 10^16 eV with beamforming.","key_machinery":"The load-bearing machinery is the fluence-footprint fit: simulated CoREAS radio pulses, convolved with the SKALA4 antenna model and a unit-gain de-dispersion filter, are interpolated across the full core using the high-precision interpolation method of [36]; after noise whitening and a 24 ns fluence window, each mock shower is fit against an ensemble of 140 simulated showers, and the Xmax of the best-fitting model is read off from a parabola fit to the lower envelope of chi-squared versus Xmax. Beamforming enters as a patch-wise coherent sum of 4, 16, or 64 neighboring antennas, raising signal-to-noise by the square root of the group size, which effectively lowers the detectable energy threshold.","core_discovery":"The central claim is that a first complete simulation of SKA-Low air-shower signals, reconstructed with the LOFAR technique, already yields a baseline performance far beyond current instruments: Xmax precision of 5 to 8 g/cm² from $10^{16}$.6 to $10^{18}$ eV, bias below 1.5 g/cm², an energy threshold of about $10^{16}$ eV when 4-, 16-, or 64-antenna beamforming groups are used, and only a moderate degradation under ADC clipping above roughly $10^{17}$.5 eV. The precision floor is not set by signal-to-noise or antenna count but by shower-to-shower variations in longitudinal development beyond Xmax, implying that finer details of the shower profile become measurable with SKA-Low. The paper frames these numbers as a conservative lower limit, since newer methods exploiting the full 50 to 350 MHz bandwidth and additional longitudinal parameters are expected to improve on them.","pith_inferences":["Because the precision floor is dominated by neglected shape parameters, the path to pushing below 5 g/cm² is likely wider bandwidth and profile-aware reconstruction rather than more antennas; this is a testable prediction that profile-based methods should beat the baseline reported here.","The paper's robustness check with the older SKALA2 antenna model suggests that small antenna-model differences are not critical, which weakens the otherwise main concern about the missing full system response; this could be confirmed by re-running the pipeline with measured responses once the real signal chain is characterized.","The clipping study implies an observable signature: if SKA-Low operates with a lower noise level than assumed, the measured Xmax bias at high energies would appear as an apparent lightening of composition above 10^17.5 eV, a cross-check that could be performed once real data accumulate.","The interpolation-based footprint method used here could in principle be applied to data from any dense radio array, so the same pipeline might serve as a community baseline for other high-density instruments, though the paper itself does not claim this."],"forward_implications":["If the simulated performance holds, SKA-Low will measure cosmic-ray mass composition across the knee-to-ankle energy range with roughly 2.5 to 4 times the Xmax precision of LOFAR, directly testing models of a Galactic-to-extragalactic transition.","The demonstration that 1000 to 3000 showers per energy bin make an Xmax-only composition analysis systematics-limited means the main return from SKA-Low's event rate will come from methods that extract the full longitudinal shower profile, not from more statistics alone.","Beamforming extends radio air-shower measurements down to 10^16 eV, a range previously hard to reach with other radio arrays, potentially enabling composition studies around the knee with a single instrument.","The dynamic-range analysis gives concrete input to SKA system design: with an 8-bit ADC and noise set at the third bit, clipping begins to bias Xmax above about 10^17.5 eV, so mitigation strategies or a different noise-level setting are needed for the highest energies.","The estimated event rate of roughly 15,000 showers per observing year above 10^16 eV means SKA-Low could, in one year, accumulate enough Xmax measurements to reach the systematics-dominated regime in several energy bins simultaneously."],"supporting_citations":[{"why":"Supplies the two-dimensional radio-intensity profile reconstruction method that the paper adapts to SKA-Low.","marker":"[4]"},{"why":"Establishes the LOFAR Xmax precision and composition results that serve as the performance baseline and systematic budget.","marker":"[16]"},{"why":"CORSIKA Monte Carlo code underlying the air-shower simulations.","marker":"[33]"},{"why":"CoREAS code that generates the simulated air-shower radio pulses.","marker":"[34]"},{"why":"High-precision interpolation method used to produce fluence footprints across the full SKA-Low core.","marker":"[36]"},{"why":"NuRadioMC framework used to add realistic Galactic and electronic noise to the simulated voltage traces.","marker":"[38]"},{"why":"Galactic radio emission model (GSM2016) providing the sky-noise level.","marker":"[39]"},{"why":"Sibyll 2.3d hadronic interaction model assumed in the mass-composition bootstrap analysis.","marker":"[52]"}],"fun_headline_variants":["SKA-Low simulations hit 5-8 g/cm² cosmic-ray depth precision","Beamforming lets SKA-Low probe cosmic rays down to 10^16 eV","SKA-Low beats LOFAR precision in cosmic-ray shower simulations","Conservative estimate: SKA-Low measures shower depth to 5-8 g/cm²","Simulated SKA-Low reaches 5-8 g/cm² Xmax precision with beamforming"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central assumption is that the simulated SKALA4 antenna response plus the simplified de-dispersion filter faithfully represents the real SKA-Low signal chain, since no complete model of the full system response exists yet.","fun_headline_variants_meta":{"raw":{"variants":["SKA-Low simulations hit 5-8 g/cm² cosmic-ray depth precision","Beamforming lets SKA-Low probe cosmic rays down to 10^16 eV","SKA-Low beats LOFAR precision in cosmic-ray shower simulations","Conservative estimate: SKA-Low measures shower depth to 5-8 g/cm²","Simulated SKA-Low reaches 5-8 g/cm² Xmax precision with beamforming"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000947,"raw_usage":{"total_tokens":4039,"prompt_tokens":940,"completion_tokens":3099,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":2986}},"tokens_in":556,"tokens_out":3099,"duration_ms":20323,"temperature":1.0,"reasoning_tokens":2986,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:54:10.642513+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Once the real SKA-Low signal chain is characterized, measure the end-to-end amplitude and phase response of a SKALA4 element and receiver across 50 to 350 MHz and feed those measurements into the same simulation and reconstruction code; if the Xmax precision degrades beyond about 10 g/cm², the bias exceeds a few g/cm², or the beamformed energy threshold rises above $10^{16}$ eV, the paper's central performance claims would be contradicted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"CORSIKA Monte Carlo code underlying the air-shower simulations."},{"cited_title":"A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers","cited_arxiv_id":"2306.13514","evidence_quote":"High-precision interpolation method used to produce fluence footprints across the full SKA-Low core."}],"review_version":1}