{"id":"913adaef-c13f-4baf-a642-309ab6c4a89a","arxiv_id":"2507.20555","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"An IFT+JAX Bayesian pipeline infers the space-time current-density field of an air shower from synthetic antenna signals, reproducing the main structure of the simulated shower.","lead":"The paper shows that Bayesian inference on radio antenna data can reconstruct the evolving electric currents in a simulated cosmic-ray air shower. It combines information field theory with JAX and GPU computing, a first step toward daytime imaging of air showers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Benchmark is an inverse crime: synthetic traces are generated with the same forward operator and prior used in reconstruction, so it cannot validate the physical fidelity of the Green's-function model; independent simulation or data is required.","rationale":"After reading the paper through, the strongest and most defensible objection is the one the Reader flagged: the benchmark is circular with respect to the forward model. I agree with the Reader's weakest_assumption. The paper is careful in many ways—it implements a nontrivial IFT/JAX pipeline, acknowledges the benchmark is simplified, and includes posterior uncertainty—but the evidence it presents is exactly the kind of self-consistency check that cannot carry the weight of 'enables a complete characterization.' I did not find an internal mathematical contradiction in the Green's function derivation; the dimensional structure appears consistent, and the method plausibly solves the stated inverse problem in simulation. The concern is thus about external validity and about quantitative calibration, not about fraud or sloppiness. The unexplained underestimation in Fig. 7 strengthens, rather than replaces, the inverse-crime concern. The one test that would settle the matter is an out-of-sample validation against an independent simulator or real data, together with a reporting of coverage and bias. Because the paper is explicitly a 'functional prototype' and claims only an 'initial pathway' in the abstract, the conditional verdict is appropriate; independent validation would upgrade it, and its absence leaves the central claim unproven rather than disproven.","tokens_in":10421,"tokens_out":10694,"duration_ms":125854,"concrete_test":"Run the inference pipeline unchanged on mock voltage traces produced by an independent, physically validated air-shower radio-emission simulator (e.g., CoREAS or CORSIKA 8) for a known shower geometry and unknown currents, then compare the posterior mean and 1σ intervals of j_k and per-voxel current vectors against the simulator's ground truth. If the reconstructed j_k is systematically biased or the intervals do not cover the truth, the inverse-crime concern lands and the 'complete characterization' claim must be withdrawn or substantially qualified. As a cheaper auxiliary check, inject a model-mismatch version of the benchmark by replacing the straight-line n_eff in Eq. (3.1) with a curved-ray travel time and see whether the inferred currents shift by more than the posterior uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that antenna voltage traces can be inverted into the time-dependent current density of a real air shower—requires that the forward operator R used in reconstruction (Eqs. 2.7–2.9 and 3.1–3.5) be a faithful model of air-shower radio emission, and that the Sec. 5 Gaussian stationary prior not dominate the answer. The Sec. 7 benchmark cannot test either condition. The synthetic shower is generated by drawing ξ_synt from the same Gaussian-process prior used in inference and applying the identical R; the observed traces therefore lie exactly in the model's prior support, so a good fit (Fig. 6) demonstrates only that MGVI can invert the generative model, not that R is physically correct. Errors in R—e.g., the straight-line effective-refractive-index time delay, the far-field transverse projection encoded by D, or the stationarity and isotropy assumptions on the currents—would be invisible to this test. Independent support is also missing: no code or data release is referenced, and no comparison to an external simulator or measured air-shower data is made. A further internal warning is that even in this idealized self-consistency test, the only quantitative summary j_k is visibly underestimated relative to the synthetic shower (Fig. 7), yet the paper does not quantify the bias, report coverage of the 1σ band, or explain why it occurred. Since the abstract and summary claim 'complete characterization' of the electromagnetic component, the evidence offered is currently insufficient to distinguish prior-driven model inversion from genuine physical reconstruction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a Bayesian inference pipeline for reconstructing the space-time dependent electric current density of an air shower from antenna voltage traces. The forward model is built on the Lorentz reciprocity theorem and a Green's function with an effective refractive index, and the shower is discretized as a moving disk of current-carrying voxels. The reconstruction uses NIFTy with Metric Gaussian Variational Inference (MGVI), with JAX providing gradients and GPU acceleration. The method is benchmarked on a synthetic shower generated from the same model and prior, and the paper shows reconstructed current maps (Fig. 4), antenna traces (Fig. 6), and a depth-dependent magnitude curve (Fig. 7). The central claim is that this is a functional prototype for imaging air showers from radio data.","tokens_in":10694,"tokens_out":3922,"duration_ms":42862,"significance":"If it works as claimed, the pipeline is a useful technical contribution: it combines IFT/MGVI with JAX for a high-dimensional four-dimensional inverse problem in astroparticle physics, and the forward model is derived from electrodynamics rather than treated as a black box. The authors are appropriately modest in calling this a prototype, and the derivation of the reciprocity-based forward operator is a clear strength. However, the validation is purely a self-consistency test on synthetic data, and the physical significance of the reconstructed currents remains unestablished. The paper would be substantially strengthened by an external forward-model comparison and by quantitative posterior diagnostics; as it stands, the evidence supports the feasibility of the numerical inversion, not the physical fidelity of the reconstruction.","major_comments":[{"comment":"The benchmark in Sec. 7 is an inverse-crime test: the synthetic data are generated by drawing ξ_synt from the same Gaussian-process prior (5.3) and applying the same operator R (Sec. 3, Eq. 4.3) that is later used in the reconstruction. The good agreement in Fig. 6 therefore demonstrates only that MGVI can invert the generative model used in the paper; it cannot validate the straight-line effective-refractive-index Green's function, the transverse-projection matrix D, or the stationarity and isotropy assumptions of the prior. The statement in Sec. 7 that 'this benchmark provides evidence of the functionality of the inference method' should be limited to numerical self-consistency, and the Sec. 8 claim of 'complete characterization of the electromagnetic component of particle showers' is not supported by the presented benchmark.","section":"Sec. 7"},{"comment":"The only quantitative reconstruction metric, j_k in Fig. 7, lies systematically below the synthetic curve over essentially the entire depth range, but the paper neither quantifies this bias, reports the coverage of the 1 sigma band, nor discusses its origin, for example prior shrinkage through the power-spectrum hyperparameters ξ_Φ or the circular mask. Since this plot is the quantitative basis for the claim of functionality, the missing bias analysis is a load-bearing omission.","section":"Fig. 7 and Sec. 7"},{"comment":"No comparison with an independent forward simulation such as CORSIKA 8/CoREAS or with measured air-shower data is provided, and no code or data release is referenced. Without such an external reference, the benchmark cannot distinguish an accurate physical forward model from one that is merely self-consistent, so the paper's central claim about inferring real air-shower currents is not yet established. This is the key missing validation step for a revision.","section":"Sec. 7 and Sec. 8"}],"minor_comments":[{"comment":"The text states that the covariance matrices 'are derived from the data during inference,' but the model then assumes a Gaussian, translation-invariant prior with power spectra; only the hyperparameters, not the functional form, are inferred. This wording should be corrected.","section":"Sec. 5"},{"comment":"Equation (3.2) uses t in the interpolation weight without defining whether t is the emission time, the receiver time, or one of the boundary bins; the notation should be clarified.","section":"Eq. (3.2)"},{"comment":"The lower panel of Fig. 6 is described in the caption as residuals normalized to the estimated uncertainty, but the vertical-axis label reads '(E_reco - E_synt)/reco'; the label or the caption should be corrected to specify the normalization.","section":"Fig. 6"},{"comment":"The benchmark description does not state the voxel size, the shower-disk dimensions, the time-bin widths, or the number of antennas beyond the 15 by 15 grid and the 500 m spacing; these parameters are needed to assess the resolution and information content of the reconstruction.","section":"Sec. 7"},{"comment":"The word 'probablistically' appears in Sec. 4.1 and should be corrected to 'probabilistically.'","section":"Sec. 4.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable proof-of-concept, but the validation is entirely self-referential and the wording overstates what is shown. For a journal like JCAP, I would want either a comparison against an external air-shower simulation or an explicit reframing of the paper as a numerical-methods demonstration with the physical claims removed. Also, no code or data release is mentioned, which limits reproducibility; the authors should consider adding a reproducibility statement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Max,\n\nThe headline is: this paper does something new and does it competently as a numerical exercise, but its benchmark cannot validate the physics, and the summary goes beyond what the evidence supports.\n\nWhat is new: the combination of IFT/NIFTy with JAX to infer a space-time current-density field of an air shower from synthetic antenna traces. The Green's function from Lorentz reciprocity (via Riegler & Windischhofer) is adapted to this problem, and the time-interpolation between shower frame and observer frame is handled carefully. The GPU implementation and MGVI setup look sound. For a proof-of-concept, that is a real contribution.\n\nWhat is soft: the benchmark in Sec. 7 generates the synthetic shower by drawing from the same Gaussian-process prior and applying the same forward operator used in reconstruction. So it demonstrates that MGVI can invert the authors' generative model, not that the model describes real air-shower radio emission. The forward model's simplifications—straight-line effective-refractive-index propagation, Gaussian stationary prior—are acknowledged, but no independent simulator (e.g., CoREAS) or measured data is used to break the circularity, and no code or data release is mentioned. That is the load-bearing gap. The abstract and summary claim 'complete characterization' of the electromagnetic component, which overstates what a self-consistency test with a visible underestimation bias (Fig. 7) can show. The bias is not quantified, and there is no coverage statistic for the posterior band.\n\nI agree with the stress-test note. This should not be read as evidence that the method works on real air showers; it should be read as evidence that the inference machinery converges and produces plausible fields on data from the same model. The authors are close to saying that, but not clearly enough.\n\nProportionate judgment: the methodology is plausible and the paper is worth a serious referee, but only with major revisions: add validation on an independent simulator, release code/data, quantify the reconstruction bias, and temper the language. As a reading-group paper, it is excellent for discussing inverse crimes in Bayesian imaging.","headline":"Novel IFT/JAX proof-of-concept for air-shower current imaging, but the benchmark is a self-consistency test on the model's own prior, so the physical validation is still missing.","tokens_in":11238,"tokens_out":2591,"would_cite":true,"duration_ms":30592,"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":"Antenna radio traces can be inverted into the time-resolved current field of an air shower.","keywords":["air showers","radio emission","Bayesian inference","information field theory","current density reconstruction","Lorentz reciprocity","variational inference","GPU acceleration"],"falsifier":"Take real antenna traces from an air-shower radio array with independent fluorescence measurements of the depth of shower maximum; if the reconstructed longitudinal current profile disagrees with the fluorescence value beyond the posterior uncertainty, or if the forward model fails to predict held-out antennas within noise, the central claim is refuted.","tokens_in":10222,"feed_emoji":"📡","tokens_out":10627,"duration_ms":98020,"temperature":0.7,"pith_summary":"This paper claims that the radio signals measured by a ground-based antenna array contain enough information to reconstruct the full time-dependent electric current field of a cosmic-ray air shower, not just its direction or energy. It sets up the measurement through the Lorentz reciprocity theorem, which turns each antenna into a Green's-function receiver for the shower's moving current density, and then uses Bayesian field inference to invert the measured voltage traces. In a synthetic benchmark the inversion recovers the orientation of the current vectors slice by slice and tracks the shower's longitudinal development, though it underestimates the total current magnitude. If this holds for real showers, radio arrays could image the electromagnetic cascade continuously, at any time of day, where fluorescence telescopes only work on clear nights.","feed_headline":"Antenna radio traces map an air shower's electric currents","feed_subtitle":"Bayesian inversion of voltage traces recovers the shower's current density in space and time in a synthetic benchmark.","key_machinery":"The load-bearing object is the reciprocity Green's function $\\vec K(\\vec x,t)=K'(\\vec x,t)\\,D(\\vec x)$, built from Eqs. (2.7)-(2.9): $K'$ contains the differentiated antenna response $h'(t-n_{\\mathrm{eff}}r/c)$ with retarded-time travel at an effective refractive index, and $D$ is a $3\\times3$ matrix encoding the dipole sensitivities, $1/r^2$ scaling, and directional weights. This function turns the shower current density $\\vec j(\\vec x,T)$ into the electric field at each antenna by a space-time convolution, discretized onto a moving disk of voxels and mapped to receiver time bins by the interpolation $\\tau$. Everything else in the paper serves to invert this forward map: the current field is represented as a Gaussian process with a power-spectrum covariance, standardized to white latent variables, and the inversion is done by variational Bayes.","core_discovery":"The central discovery is a working inversion pipeline from antenna signal traces to the space-time current density of an air shower. The electric field at three orthogonal dipole antennas is written as an integral over the shower volume of a reciprocity Green's function $\\vec K = K' D$ acting on the current density $\\vec j$, where $K'$ encodes the antenna response and the retarded time $t - n_{\\mathrm{eff}} r/c$ at an effective refractive index, and the matrix $D$ encodes the dipole orientation and $1/r^2$ falloff. The shower is modeled as a thin disk of voxels, each carrying a three-component current density drawn from a Gaussian process with a spectral covariance, and the disk is moved at near-light speed through an atmosphere with a density gradient; a dedicated interpolation $\\tau$ maps emission times in the shower frame to observation times in antenna bins. A variational inference scheme approximates the posterior over the standardized latent variables by a Gaussian whose covariance is estimated from the Fisher information. On a synthetic benchmark with $5\\times4$ antennas, the reconstruction reproduces the directions of the current vectors, its total-magnitude-depth profile follows the synthetic shower (with an underestimation), and posterior uncertainties are reported; the authors state that the benchmark \"provides evidence of the functionality of the inference method.\"","pith_inferences":["If the forward model holds for real events, existing recordings from radio arrays could be reprocessed to produce a movie of each shower's cascade development, a data product the paper does not itself construct.","The straight-line effective-refractive-index Green's function is the most likely point of failure; a natural stress test is to compare reconstructions made with curved ray tracing or frequency-dependent refractivity on the same synthetic events.","The stationary Gaussian-process prior will tend to smooth sharp shower-front features; replacing it with a non-stationary or sparser prior could sharpen the reconstructed front and test how much structure the data actually constrain.","The same reciprocity-Green's-function-plus-Bayesian-inversion construction transfers to any transient radio emitter with a compact current distribution, such as lightning or laboratory discharges, once geometry and priors are adapted."],"forward_implications":["Radio-only imaging of air showers becomes feasible in daylight and poor weather, because the method reconstructs the longitudinal current profile without fluorescence light.","The reconstructed current-density field comes with posterior uncertainties, enabling event-by-event statistical statements about shower geometry and time evolution.","The GPU-accelerated inference ran in roughly 600 seconds per shower at benchmark resolution, suggesting near-real-time processing of array events is within reach.","Tracking where the reconstructed current magnitude peaks gives a radio-based route to the depth of shower maximum, a key observable for cosmic-ray composition.","Absolute current magnitude is underestimated in the benchmark, so recovering overall normalization requires further work even as shape and direction are reproduced."],"supporting_citations":[{"why":"Supplies the exact Maxwell-theory derivation of the Lorentz-reciprocity Green's function on which the forward signal model is built.","marker":"[14]"},{"why":"Provides the macroscopic treatment of air-shower radio emission and supports the straight-line effective-refractive-index propagation approximation.","marker":"[8]"},{"why":"Establishes the information-field-theory formalism (Bayesian Hamiltonians, priors on fields) used to formulate the reconstruction.","marker":"[12]"},{"why":"Implements the numerical field-inference machinery that runs the variational reconstruction.","marker":"[11]"},{"why":"Provides the GPU-accelerated automatic differentiation and vectorized array operations that make the high-dimensional inversion tractable.","marker":"[13]"},{"why":"Supplies the variational inference algorithm that approximates the posterior over shower parameters.","marker":"[21]"},{"why":"Supplies the correlated-field and power-spectrum construction used to build the spatio-temporal prior on current densities.","marker":"[16]"},{"why":"Supplies the far-field dipole expressions used in the weighting-field formula for the reciprocity Green's function.","marker":"[15]"}],"fun_headline_variants":["Bayesian inversion reconstructs air-shower currents from antenna traces","JAX and info field theory map shower currents from radio signals","Antenna traces inverted to reveal cosmic-ray shower current flows","Reconstructing air-shower currents from ground antenna data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the simplified forward model—radio waves traveling along straight lines at a single averaged refractive index, and shower currents obeying a smooth Gaussian prior—truly describes real air showers; the benchmark cannot test this, since its synthetic data come from the same model.","fun_headline_variants_meta":{"raw":{"variants":["Bayesian inversion reconstructs air-shower currents from antenna traces","JAX and info field theory map shower currents from radio signals","Antenna traces inverted to reveal cosmic-ray shower current flows","Reconstructing air-shower currents from ground antenna data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00021,"raw_usage":{"total_tokens":1411,"prompt_tokens":947,"completion_tokens":464,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":563,"completion_tokens_details":{"reasoning_tokens":395}},"tokens_in":563,"tokens_out":464,"duration_ms":5465,"temperature":1.0,"reasoning_tokens":395,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:40:52.306503+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take real antenna traces from an air-shower radio array with independent fluorescence measurements of the depth of shower maximum; if the reconstructed longitudinal current profile disagrees with the fluorescence value beyond the posterior uncertainty, or if the forward model fails to predict held-out antennas within noise, the central claim is refuted.","supporting_citations":[{"cited_title":"Signals induced on electrodes by moving charges, a general theorem for Maxwell's equations based on Lorentz-reciprocity","cited_arxiv_id":"2001.10592","evidence_quote":"Supplies the exact Maxwell-theory derivation of the Lorentz-reciprocity Green's function on which the forward signal model is built."},{"cited_title":"Macroscopic Treatment of Radio Emission from Cosmic Ray Air Showers based on Shower Simulations","cited_arxiv_id":"0712.2517","evidence_quote":"Provides the macroscopic treatment of air-shower radio emission and supports the straight-line effective-refractive-index propagation approximation."},{"cited_title":"Bradbury, R","cited_arxiv_id":null,"evidence_quote":"Provides the GPU-accelerated automatic differentiation and vectorized array operations that make the high-dimensional inversion tractable."},{"cited_title":"Arras, P","cited_arxiv_id":null,"evidence_quote":"Supplies the correlated-field and power-spectrum construction used to build the spatio-temporal prior on current densities."},{"cited_title":"Kraus, Antennas, McGraw-Hill, New York, NY, 2nd ed","cited_arxiv_id":null,"evidence_quote":"Supplies the far-field dipole expressions used in the weighting-field formula for the reciprocity Green's function."}],"review_version":2}