{"id":"6ea6f6ee-d18f-4581-88a2-a20baad562e2","arxiv_id":"2507.06874","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"An analytic chi-squared minimization reconstructs air-shower electric fields from three-polarization radio antennas with 4-6% precision in simulations.","lead":"A new analytic chi-squared method reconstructs the full electric field of inclined cosmic-ray air showers from radio antennas with three polarizations. It achieves 4% precision on the field amplitude and 6% on energy fluence in simulations, improving on standard matrix inversion.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The E_r=0 assumption is load-bearing: the paper never quantifies the radial field in its ZHAireS library, and Section 3.1 itself shows a small r-component corrupting the 3-polarization matrix inversion.","rationale":"The derivation of the weighted least-squares solution is standard and internally consistent: for a known antenna response and diagonal noise covariance, Eq. (1) gives the maximum-likelihood estimate of the fitted components. The validation on two antenna models, dipole and HORIZON, is a genuine stress test, and the reported histograms show a real improvement over the 2x2 matrix-inversion baseline. The concern is not with the algebra but with the physical truncation to two components. The paper's own Section 3.1 states that a small r-component degrades the three-polarization matrix inversion, so E_r is not exactly zero in the very simulations used for validation; yet the new method's reliance on E_r=0 is never checked against those simulations. This absence is directly relevant to the headline precision on total field and fluence, and it can be settled with an analysis of the already-stored ZHAireS electric-field traces, without new simulations. Because the requested check is feasible and the current paper lacks the necessary quantification, the CONDITIONAL verdict is appropriate and unchanged by this stress test.","tokens_in":9704,"tokens_out":11811,"duration_ms":127469,"concrete_test":"From the existing 4,160 ZHAireS events, take a random subset of about 200 events spanning the full zenith range (63 to 87 degrees) and core distances. For each selected antenna, compute the simulated radial component E_r = E·r_hat along the unit vector from the emission region to the antenna, and form the Hilbert-peak ratio |E_r|_peak / sqrt(|E_theta|_peak^2 + |E_phi|_peak^2). Rerun the authors' two-component chi2 reconstruction and regress the per-trace total-PEA and fluence relative errors against this ratio and against zenith/core distance. If the errors grow systematically with the radial fraction, or if the median ratio exceeds about 5-10%, the E_r=0 assumption is the dominant limitation of the closed-form solution and the 4%/6% claim is not supported as stated. If the ratio is uniformly below about 1% and the errors show no trend with it, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central 4%/6% precision claim rests on the closed-form solution of Eq. (1), which fits only E_theta and E_phi and sets E_r=0 on transverse-wave grounds. The validation metric, however, compares this two-component reconstruction with the full simulated electric field, which contains E_r. If E_r is not negligible, the closed-form estimate is biased in exactly the total-PEA and fluence metrics used to support the abstract claim; if E_r is negligible, that should be demonstrated, but the paper reports no |E_r|/|E| distribution anywhere. The manuscript itself gives evidence that E_r is present: in Section 3.1, including the third polarization in the conventional matrix inversion 'slightly degrades the stronger phi component due to artifacts introduced by the small r-component'. The chi2 method removes those artifacts by projecting E_r out, but any true radial component is then lost from the total field or leaks into the fitted components. The reported dipole bias (~5% PEA, ~3% fluence) is left unexplained and could be a projection loss rather than an antenna-response effect. This matters most for inclined, near-core antennas, where the analysis intentionally keeps the innermost 16 antennas per arm and where curved wavefronts and near-field contributions make E_r largest. Thus the E_r=0 assumption is the least secure condition for the central claim, and the claim that the method is free of assumptions about signal properties is overstated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytic chi-squared minimization method for reconstructing the electric field of inclined air showers from antenna voltage traces in three polarizations. The method solves a frequency-domain weighted least-squares problem for the two transverse field components E_theta and E_phi, setting the radial component E_r to zero, and is validated against ZHAireS simulations with both a simple dipole and the realistic HORIZON antenna, including simulated galactic noise. The authors report standard deviations better than 4% for the Hilbert peak envelope amplitude (PEA) and better than 6% for the energy fluence, with an antenna-response-dependent bias, and claim improved performance when a vertical polarization channel is included.","tokens_in":9917,"tokens_out":3263,"duration_ms":34899,"significance":"If the central claims hold, the method provides a simple, closed-form, and assumption-light reconstruction of the electric field for inclined air showers, which is directly relevant to GRAND and other three-polarization radio-detection experiments. The derivation is transparent and parameter-free, and the validation uses a large, independent simulation library (4,160 ZHAireS events) with realistic antenna responses and noise. The explicit comparison against the conventional matrix-inversion method clarifies the advantage of the chi-squared approach. The main caveats—the unquantified radial field component, the omitted theta-component resolution numbers in the abstract, and the unexplained dipole bias—currently prevent the claims from being fully supported.","major_comments":[{"comment":"The model reduces the electric field to two components by setting E_r = 0, but the validation compares the reconstructed two-component field with the full simulated electric field that contains E_r. The manuscript never reports the distribution of |E_r|/|E| in the ZHAireS library, so the magnitude of the resulting projection bias is unknown. Section 3.1 itself states that the small r-component introduces artifacts in the three-polarization matrix inversion, demonstrating that a nonzero radial component is present in the simulated data. The total-PEA and fluence metrics used in the abstract are exactly those most sensitive to such a bias. The authors should quantify E_r (for example, histograms of radial-to-total amplitude versus zenith angle and core distance) and either show that it is negligible for the selected events or extend the fit to include E_r.","section":"Section 3.2, Eq. (1)"},{"comment":"The abstract claims standard deviations better than 4% for PEA and better than 6% for energy fluence, but Figures 3 and 4 show that the theta-component reconstruction for the HORIZON antenna has a standard deviation of 12% for PEA and 19% for fluence. These exceptions are acknowledged in Sections 4.1 and 4.2 but are not reflected in the abstract or the conclusions, where the blanket statement 'reducing the typical standard deviation to 4% for PEA and less than 6% for energy fluence' appears without qualification. The abstract and conclusions should explicitly state that the headline precision applies to the total and phi components, and report the theta-component numbers.","section":"Abstract; Sections 4.1 and 4.2, Figs. 3 and 4"},{"comment":"The text states that the dependence of the method on arrival direction was studied and that 'it has a good performance in the zenith range from 63 up to 80 degrees,' but no figure, table, or quantitative result for this zenith dependence appears anywhere in the manuscript. The simulation library extends to 87.1 degrees, so the behavior beyond 80 degrees is also unaddressed. The claim is therefore unsupported as written; the authors should add a zenith-binned performance plot or table, or explicitly label this as a stated intent rather than a demonstrated result.","section":"Section 4, first paragraph"},{"comment":"The dipole antenna shows a consistent bias of about 5% in PEA and about 3% in fluence (an underestimation), while the HORIZON antenna shows near-zero bias. The manuscript attributes this to antenna response, frequency coverage, and directional sensitivity, but states that it is 'still under investigation.' Because the abstract advertises an antenna-response-dependent bias as a feature of the method, this unexplained systematic is load-bearing for the claimed accuracy. The authors should either provide a mechanism and correction for the bias or demonstrate that it does not arise from the E_r projection loss discussed in the first major comment.","section":"Section 4.1 and 4.2, Figs. 3 and 4"}],"minor_comments":[{"comment":"There are several typographical errors: 'minization' in the Section 3.2 heading, 'thesimulated' in Section 4.1, 'refinehe' in Section 5, and 'PAE' instead of 'PEA' in the first paragraph of Section 4.1. These should be corrected in a final pass.","section":"Throughout"},{"comment":"The notation for the noise covariance is inconsistent: sigma_V is first defined as the background noise level and then as the diagonal covariance matrix with squared noise spectra. Please use distinct symbols for the noise standard deviation and the covariance matrix to avoid confusion.","section":"Section 3.2, Eq. (1)"},{"comment":"The captions and legends are difficult to read because the color-matching between the histograms and the statistical summaries is not explicit. Please add direct labels to each histogram (e.g., 'HORIZON' and 'Dipole') and clarify which distribution corresponds to which antenna.","section":"Figures 3 and 4"},{"comment":"The manuscript notes that electronic noise between V_oc and V_ADC is not modeled. This is a limitation for real-data application, and it would be helpful to state explicitly in the conclusions that the covariance matrix in Eq. (1) would need to incorporate electronic noise for experimental data.","section":"Section 2.3"},{"comment":"Reference [7] is cited as 'ARENA. 8, 2022,' which is not a standard citation format. Please provide the full proceedings information or arXiv identifier, as is done for other references.","section":"Reference [7]"}],"recommendation":"major_revision","confidential_remarks":"This is an ICRC proceedings contribution, so the level of detail is appropriate for a conference paper, but the central numerical claims must be represented accurately. The main technical concern—the unquantified E_r component—is not a flaw in the least-squares derivation itself, but it directly affects the headline precision numbers, and the authors have the simulation data already in hand to test it. I would encourage the editor to request the additional E_r quantification and the zenith-dependence plot before final acceptance, as these are feasible within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the analytic chi-squared minimization for three-polarization electric-field reconstruction. The normal-equations solution is textbook linear algebra, but no one in the cited radio-detection literature (AERA matrix inversion, forward folding, information field theory, LOPES-3D) has formulated it this way. The paper validates it on a large ZHAireS library, 4,160 events, with two antenna models and simulated galactic noise. That is real, reproducible work, and the cross-correlation gains over matrix inversion are convincing, especially the theta-component improvement for the dipole. The total-field PEA and fluence histograms do support the 4%/6% standard-deviation claims. The soft spots are real but manageable. The abstract says 'better than 4% and 6%' and then hides the HORIZON theta-component exceptions, which are 12% for PEA and 19% for fluence. That is an overstatement by any standard, and it needs to be fixed in the abstract. The claimed zenith-range performance from 63 to 80 degrees is also not shown anywhere in the paper; the dataset spans that range, but there is no figure or table demonstrating the dependence. The dipole shows a consistent ~5% PEA bias and ~3% fluence bias that the authors say is under investigation; I would want that understood before relying on the method for absolute energy calibration. On the E_r=0 assumption: the stress-test is right that the paper never quantifies the radial component in its own simulation library. But the validation compares the reconstruction to the full 3D simulated electric field, which contains whatever E_r actually exists. If E_r were large, the 4%/6% numbers would not hold. So the empirical agreement suggests E_r is small for these inclined showers, but the paper should still show the |E_r|/|E| distribution and discuss projection loss. The unexplained dipole bias could in principle be a projection artifact, so quantifying E_r would also help there. This is a caveat, not a fatal flaw. This paper is a solid methods contribution for GRAND and other inclined-radio arrays. It ships no code, but the derivation is self-contained and the validation is against independent simulated ground truth. The lack of electronic noise modeling is acknowledged and appropriate for a simulation study. For peer review: yes, it deserves refereeing. A serious editor should send it out, with the expectation of an abstract revision and a request for the missing zenith-dependence figure. I would not desk-reject it.","headline":"A genuinely useful, straightforward three-polarization E-field reconstruction method with solid simulation validation, but the abstract overstates the theta-component and zenith-range results.","tokens_in":678,"tokens_out":1000,"would_cite":true,"duration_ms":33110,"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 closed-form chi^2 minimization of three-polarization voltage traces reconstructs inclined air-shower electric fields with better-than-4% peak envelope accuracy and better-than-6% energy-fluence accuracy.","keywords":["electric-field reconstruction","radio detection of air showers","inclined air showers","three-polarization antennas","chi-square minimization","energy fluence","ZHAireS simulations","antenna response"],"falsifier":"Compute the radial field component $E_r$ in the same ZHAireS simulation traces and compare its Hilbert peak envelope with the transverse components; if $E_r$ reaches a few percent of the total field, the $E_r=0$ assumption will bias the peak envelope beyond the claimed 4%. Then re-run the reconstruction with $E_r$ included as a third fitted unknown using the three polarization channels: if the three-component closed-form solution gives materially different peak-envelope and fluence statistics, the paper's two-component solution is the one that is incomplete.","tokens_in":9486,"feed_emoji":"📡","tokens_out":8763,"duration_ms":85070,"temperature":0.7,"pith_summary":"The radio technique for detecting cosmic-ray air showers needs the electric field at each antenna, but antennas record voltage traces, so the field must be recovered by inverting the antenna response. This paper proposes an analytic $chi^{2}$ minimization that performs that inversion for antennas with two or, more usefully, three polarizations, in one closed-form noise-weighted least-squares step. Tested on a ZHAireS simulation library of inclined showers with realistic antenna responses and galactic noise, the method reproduces the Hilbert peak envelope amplitude with standard deviations below 4% and the square-root energy fluence below 6%, with an antenna-dependent bias. The authors argue this is a model-independent alternative to forward folding and matrix inversion, and that adding a vertical polarization channel meaningfully improves the reconstruction for inclined events.","feed_headline":"Three-polarization fit recovers air-shower radio fields to 4%","feed_subtitle":"One noise-weighted inversion of all three antenna polarizations also pins the energy fluence to under 6 percent.","key_machinery":"The machinery is the closed-form $chi^{2}$ minimizer $E = (H^T \\sigma_V^{-1} H)^{-1} H^T \\sigma_V^{-1} V$, a noise-weighted least-squares inversion performed per frequency bin. The matrix $H$ is the antenna response relating the two retained field components ($E_\\theta$, $E_\\varphi$) to the three voltage channels, and $\\sigma_V$ is the per-channel noise power. The noise weighting suppresses the low-gain or low-SNR frequency bins that degrade the standard matrix inversion, and dropping $E_r$ via the transverse-wave assumption reduces the inversion from three unknowns to two.","core_discovery":"The central claim is that the electric field of an inclined air shower can be recovered from three-polarization voltage traces by minimizing a single frequency-domain $chi^{2}$ whose analytic solution is $E = (H^T \\sigma_V^{-1} H)^{-1} H^T \\sigma_V^{-1} V$. Here $H$ maps the two transverse field components $E_\\theta$ and $E_\\varphi$ to the three antenna channels, $\\sigma_V$ is the diagonal noise spectrum per channel, and the radial component $E_r$ is set to zero because the wave is transverse. In the ZHAireS validation over thousands of traces, the method achieves standard deviations better than 4% for the Hilbert peak envelope amplitude and better than 6% for the square-root energy fluence, and it cross-correlates with the true field much better than the conventional matrix-inversion method, particularly for the weaker $\\theta$ component. The paper also shows that including the vertical polarization channel is what makes the precise three-polarization reconstruction possible.","pith_inferences":["The paper does not quantify $E_r$ in its own ZHAireS library; a direct measurement of that component, which would settle the transverse-wave assumption, is a natural next step.","If the transverse-wave assumption holds for other geometries, the same closed-form inversion should transfer to any three-polarization radio detector with a known antenna response.","The 12-19% scatter on the weak $\\theta$ component suggests that direction-refined iteration, which the authors list as future work, may be required before polarization-resolved physics uses that channel.","The consistent few-percent dipole bias implies a per-antenna calibration factor could absorb most of the residual, and testing that correction on experimental data would be a direct confirmation."],"forward_implications":["Three-polarization antennas become directly useful for inclined-shower arrays, since the added vertical channel measurably tightens the reconstructed field.","A single closed-form matrix product per frequency bin recovers the field, so the method can run in real time without iterative fitting.","Energy fluence estimates with better-than-6% typical standard deviation bring primary-energy measurements closer to the statistical limit set by the detector.","Frequency-domain noise weighting makes the reconstruction stable at low SNR and low antenna gain without assuming a signal template.","The verified zenith range of 63 to 80 degrees covers the geometrically most important window for very inclined air-shower radio detection."],"supporting_citations":[{"why":"Supplies the conventional two-polarization matrix-inversion method that serves as the baseline.","marker":"[1]"},{"why":"Provides the forward-folding chi^2 framework whose analytic form the new method is inspired by.","marker":"[3]"},{"why":"Reports a prior three-polarization study that found no gain with simple weighting, motivating a better method.","marker":"[6]"},{"why":"Generates the ZHAireS simulation library of inclined air showers used to validate the method.","marker":"[8]"},{"why":"Defines the realistic three-polarization antenna model used alongside the simple dipole.","marker":"[10]"},{"why":"Simulates the antenna responses, including ground effects, that form the measurement matrix.","marker":"[12]"},{"why":"Provides the galactic-noise model used to create realistic voltage traces.","marker":"[15]"}],"fun_headline_variants":["Three-polarization chi^2 cuts air-shower field error to 4%","Inclined air-shower radio fields pinned to 4% with 3-polarization fit","Analytic chi^2 gives 4% field, 6% fluence for inclined showers","Three-axis radio reconstruction beats 4% for inclined air showers","Polarization trio sharpens air-shower electric-field reconstruction"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes the radio wave is purely transverse, so the electric field along the propagation direction (the radial component) is zero; if a measurable radial component exists, the closed-form solution is systematically biased, and the paper does not quantify that component in its simulations.","fun_headline_variants_meta":{"raw":{"variants":["Three-polarization chi^2 cuts air-shower field error to 4%","Inclined air-shower radio fields pinned to 4% with 3-polarization fit","Analytic chi^2 gives 4% field, 6% fluence for inclined showers","Three-axis radio reconstruction beats 4% for inclined air showers","Polarization trio sharpens air-shower electric-field reconstruction"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001117,"raw_usage":{"total_tokens":4658,"prompt_tokens":963,"completion_tokens":3695,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":579,"completion_tokens_details":{"reasoning_tokens":3591}},"tokens_in":579,"tokens_out":3695,"duration_ms":30500,"temperature":1.0,"reasoning_tokens":3591,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:52:29.343112+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the radial field component $E_r$ in the same ZHAireS simulation traces and compare its Hilbert peak envelope with the transverse components; if $E_r$ reaches a few percent of the total field, the $E_r=0$ assumption will bias the peak envelope beyond the claimed 4%. Then re-run the reconstruction with $E_r$ included as a third fitted unknown using the three polarization channels: if the three-component closed-form solution gives materially different peak-envelope and fluence statistics, the paper's two-component solution is the one that is incomplete.","supporting_citations":[{"cited_title":"Welling, C","cited_arxiv_id":null,"evidence_quote":"Provides the forward-folding chi^2 framework whose analytic form the new method is inspired by."},{"cited_title":"Huber,Analysing the electric field vector of air shower radio emission","cited_arxiv_id":null,"evidence_quote":"Reports a prior three-polarization study that found no gain with simple weighting, motivating a better method."},{"cited_title":"Alvarez-Muñiz, W","cited_arxiv_id":null,"evidence_quote":"Generates the ZHAireS simulation library of inclined air showers used to validate the method."},{"cited_title":"Ansys HFSS, High Frequency Structure Simulator","cited_arxiv_id":null,"evidence_quote":"Simulates the antenna responses, including ground effects, that form the measurement matrix."},{"cited_title":"PolisenskyLong Wavelength Array Memo Series111(2007) 515","cited_arxiv_id":null,"evidence_quote":"Provides the galactic-noise model used to create realistic voltage traces."}],"review_version":1}