Electric-Field Reconstruction for Radio Detection of Inclined Air Showers in Three Polarizations

Chao Zhang; Kewen Zhang; Lukas G\"ulzow; Mat\'ias Tueros; Pengfei Zhang; Pengxiong Ma; Ramesh Koirala; Tim Huege; Xin Xu; Yi Zhang (for the GRAND collaboration)

REVIEW 4 major objections 5 minor 13 references

Electric-Field Reconstruction for Radio Detection of Inclined Air Showers in Three Polarizations

T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2507.06874 v2 pith:LGNU47EK submitted 2025-07-09 astro-ph.IM astro-ph.HE

Kewen Zhang , Lukas Gülzow , Tim Huege , Ramesh Koirala , Pengxiong Ma , Matías Tueros , Xin Xu , Chao Zhang

show 2 more authors

Pengfei Zhang Yi Zhang (for the GRAND Collaboration)

This is my paper · ORCID

classification astro-ph.IMastro-ph.HE

keywords electric-fieldreconstructionradiodetectionofairshowersinclinedthree-polarizationantennaschi-squareminimizationenergyfluenceZHAireSsimulationsantennaresponse

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

Referee Report

4 major / 5 minor

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.

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 (4)

[Section 3.2, Eq. (1)] 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.
[Abstract; Sections 4.1 and 4.2, Figs. 3 and 4] 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 4, first paragraph] 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 4.1 and 4.2, Figs. 3 and 4] 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.

minor comments (5)

[Throughout] 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 3.2, Eq. (1)] 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.
[Figures 3 and 4] 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 2.3] 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.
[Reference [7]] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the reconstruction is a closed-form weighted least-squares inversion validated against independent ZHAireS simulations.

full rationale

The derivation chain is self-contained. The analytic chi-squared minimization solves E = (H^T sigma_V^{-1} H)^{-1} H^T sigma_V^{-1} V directly from the voltage, antenna-response, and noise covariance inputs (Eq. 1), with no parameters fitted to the validation outcomes. The reported 4% and 6% standard deviations are empirical scatter of reconstructed versus simulated quantities, where the simulation ground truth is generated independently by ZHAireS (Section 2.1). No equation in the paper reduces the predicted PEA or fluence to a fitted value, and no load-bearing claim relies solely on a self-citation; reference [3] is cited as inspiration for the approach rather than as proof of this method's performance. The E_r=0 assumption is an explicit physical modeling choice, and while it is a limitation that could bias the reconstruction if violated, an untested assumption is not circularity under the criteria used here. The validation metrics compare against the full simulated field, so the assumption is externally checkable; the paper's central claim remains an independently supported empirical result.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The reconstruction rests on weighted least squares and four explicit domain assumptions. No free parameters are fit to data and no new entities are introduced. The E_r = 0 assumption is the most fragile: the paper does not quantify the radial component in its own ZHAireS data, so a non-negligible E_r would introduce a systematic bias in the two-component solution.

assumptions (5)

domain assumption The air-shower radio signal is a transverse wave, so the radial electric-field component vanishes: E_r = 0.
Invoked in Section 3.2 before Eq. (1): 'E_r is neglected here as it corresponds to the propagation direction of the transverse wave.' The entire two-component solution assumes this.
domain assumption The measurement error is dominated by Gaussian background noise that is uncorrelated between polarization channels, with a known power spectrum.
Stated near Eq. (1): 'the covariance matrix is constructed as a diagonal matrix diag(sigma_V1, sigma_V2, sigma_V3)'; noise spectrum is taken from LFmap at LST=18h. Correlated or mis-modeled noise would bias the weights.
domain assumption The antenna response matrix H is known exactly from HFSS simulations, including ground effects.
Section 2.2 uses simulated responses for the dipole and HORIZON antennas. Any response error propagates directly into the reconstructed field.
domain assumption The input signal is the open-circuit voltage V_oc, and electronics noise between the antenna and digitizer is negligible.
End of Section 2.3: 'We assume that the input to the reconstruction algorithm is the open-circuit voltage V_oc... Electronic noise ... is not modeled here.' Real-data performance will include this noise.
standard math Standard linear algebra: the normal equations (H^T Sigma^-1 H) E = H^T Sigma^-1 V give the minimizing solution.
Used in Section 3.2 to obtain E=(H^T Sigma^-1 H)^-1 H^T Sigma^-1 V.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Electric-Field Reconstruction for Radio Detection of Inclined Air Showers in Three Polarizations." pith.science (2026). https://pith.science/paper/LGNU47EK

@misc{pith2026250706874,
  author       = {Pith},
  title        = {Pith review of: Electric-Field Reconstruction for Radio Detection of Inclined Air Showers in Three Polarizations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LGNU47EK}},
  note         = {Machine review of arXiv:2507.06874}
}

abstract

Accurate reconstruction of the electric field produced by extensive air showers is essential for the radio-detection technique, as the key parameters of interest of the primary particles that generated the showers are the amplitude, polarization, frequency spectrum, and energy fluence carried by the electric field at each receiving radio antenna. Conventional electric-field reconstruction methods primarily focus on antennas with two horizontal polarizations. In this work, we introduce an analytic $\chi^2$ minimization method that is applicable to both two and three polarizations. This solution has been verified for simple and realistic antenna responses, with a particular focus on inclined air showers. Our method achieves standard deviations better than 4\% and 6\% for the estimation of the Hilbert peak envelope amplitude of the electric field and the energy fluence, respectively, with an antenna-response-dependent bias. Additionally, we have studied the dependence of the method with arrival direction showing that it has a good performance in the zenith range from 63$^\circ$ up to 80$^\circ$. This work also demonstrates that incorporating vertically polarized antennas enhances the precision of the reconstruction, leading to a more accurate and reliable electric-field estimation for inclined air showers.

Figures

Figures reproduced from arXiv: 2507.06874 by the authors.

**Figure 1.** Comparison of the matrix-inversion method, for two (2pol) and three (3pol) polarizations, and the analytic 𝜒 2 minization method (lsq), for the simple dipole antenna. Top: Time traces of the electric field. Bottom: Electric-field spectrum in the frequency domain. field using both two and three polarizations The normalized cross-correlation values [16], ranging between -1 (perfectly anticorrelated) and +1 (perfectly … view at source ↗

**Figure 2.** Comparison of the matrix inversion method (inv) and the analytic 𝜒 2 minimization (lsq) reconstruction methods in three polarizations applied to our example electric-field trace (sim), for the HORIZON antenna. Top: time traces of the electric field. Bottom: electric-field spectrum in the frequency domain. 0.2 0.1 0.0 0.1 0.2 0.3 (E tot PEA,sim - E tot PEA,rec)/E tot PEA,rec) 0 10000 20000 30000 40000 Number of trac… view at source ↗

**Figure 3.** Relative error of reconstructed (rec) PEA w.r.t. simulated (sim) PEA for total (left), 𝜑 (middle), and 𝜃 (right) electric-field components. Results are shown for the HORIZON (green) and dipole (black) antennas, with statistical summaries in matching colors. better than 5%, with the exception of the 𝜃-component reconstruction of the HORIZON antenna, which yields a standard deviation of 12%. This is due to the weaker … view at source ↗

Figures from the paper (1 more)

**Figure 4.** Figure 4: Relative error of the fluence (Φ) reconstruction, comparing the reconstructed √ Φrec with the simulated √ Φsim for total (left), 𝜑 (middle), and 𝜃 (right) electric-field components. Results are shown for the HORIZON (green) and dipole (black) antennas, with statistical…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

[2]

Glaseret al

C. Glaseret al. Eur. Phys. J. C79no. 6, (2019) 464

work page 2019
[3]

Welling, C

C. Welling, C. Glaser, and A. NellesJCAP2019no. 10, (Oct., 2019) 075–075

work page 2019
[4]

Wellinget al

C. Wellinget al. JCAP2021no. 04, (Apr., 2021) 071

work page 2021
[5]

Electric field reconstruction with information field theory,

S. Strähnzet al., “Electric field reconstruction with information field theory,” p. 056, ARENA. 11, 2024

work page 2024
[6]

Huber,Analysing the electric field vector of air shower radio emission

D. Huber,Analysing the electric field vector of air shower radio emission. PhD thesis, Karlsruher Institut für Technologie (KIT), 2014

work page 2014
[7]

LOPES 3D -- studies on the benefits of EAS-radio measurements with vertically aligned antennas

D. Huberet al., “LOPES 3D – studies on the benefits of EAS-radio measurements with vertically aligned antennas,” ARENA. 8, 2022.arXiv:2208.14664 [astro-ph.IM]

work page Pith review arXiv 2022
[8]

Alvarez-Muñiz, W

J. Alvarez-Muñiz, W. R. Carvalho, and E. ZasAstroparticle Physics35no. 6, (Jan., 2012) 325–341

work page 2012
[9]

D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz, and T. Thouw, CORSIKA: a Monte Carlo code to simulate extensive air showers.1998. [10]GRANDCollaboration, J. Álvarez-Muñizet al. Sci. China Phys. Mech. Astron.63no. 1, (2020) 219501. [11]GRANDCollaboration, R. Alves Batistaet al.7, 2023.arXiv:2308.00120 [hep-ex]

work page Pith review arXiv 2020

Show all 13 references

[12]

Ansys HFSS, High Frequency Structure Simulator

Ansys Inc., “Ansys HFSS, High Frequency Structure Simulator.” Version 2021 r1, https://www.ansys.com/products/electronics/ansys-hfss, 2021. [13]GRANDCollaboration, J. Álvarez-Muñizet al. Science China Physics, Mechanics and Astronomy63no. 1, (Aug., 2019)

2021
[14]

Büsken, T

M. Büsken, T. Fodran, and T. HuegeAstronomy and Astrophysics679(Nov., 2023) A50

2023
[15]

PolisenskyLong Wavelength Array Memo Series111(2007) 515

E. PolisenskyLong Wavelength Array Memo Series111(2007) 515

2007
[16]

Corstanjeet al

A. Corstanjeet al. Journal of Instrumentation18no. 09, (Sept., 2023) P09005

2023
[17]

Emergences

C. Glaseret al. JCAP09(2016) 024. 8 Electric-Field Reconstruction in Three PolarizationsKewen Zhang Full Author List: GRAND Collaboration J. Álvarez-Muñiz1, R. Alves Batista2,3, A. Benoit-Lévy4, T. Bister5,6, M. Bohacova7, M. Bustamante8, W. Carvalho9, Y. Chen10,11, L.Cheng 12...

2016

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.

[1] [2]

Glaseret al

C. Glaseret al. Eur. Phys. J. C79no. 6, (2019) 464

work page 2019

[2] [3]

Welling, C

C. Welling, C. Glaser, and A. NellesJCAP2019no. 10, (Oct., 2019) 075–075

work page 2019

[3] [4]

Wellinget al

C. Wellinget al. JCAP2021no. 04, (Apr., 2021) 071

work page 2021

[4] [5]

Electric field reconstruction with information field theory,

S. Strähnzet al., “Electric field reconstruction with information field theory,” p. 056, ARENA. 11, 2024

work page 2024

[5] [6]

Huber,Analysing the electric field vector of air shower radio emission

D. Huber,Analysing the electric field vector of air shower radio emission. PhD thesis, Karlsruher Institut für Technologie (KIT), 2014

work page 2014

[6] [7]

LOPES 3D -- studies on the benefits of EAS-radio measurements with vertically aligned antennas

D. Huberet al., “LOPES 3D – studies on the benefits of EAS-radio measurements with vertically aligned antennas,” ARENA. 8, 2022.arXiv:2208.14664 [astro-ph.IM]

work page Pith review arXiv 2022

[7] [8]

Alvarez-Muñiz, W

J. Alvarez-Muñiz, W. R. Carvalho, and E. ZasAstroparticle Physics35no. 6, (Jan., 2012) 325–341

work page 2012

[8] [9]

D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz, and T. Thouw, CORSIKA: a Monte Carlo code to simulate extensive air showers.1998. [10]GRANDCollaboration, J. Álvarez-Muñizet al. Sci. China Phys. Mech. Astron.63no. 1, (2020) 219501. [11]GRANDCollaboration, R. Alves Batistaet al.7, 2023.arXiv:2308.00120 [hep-ex]

work page Pith review arXiv 2020