---
title: 'CoREAS: Monte Carlo Radio Emission Simulator'
url: https://www.emergentmind.com/topics/coreas
type: topic
---

# CoREAS: Monte Carlo Radio Emission Simulator

to=arxiv_search  天天中彩票这个json code  lower?
to=arxiv_search ाठमाडौं code 《凤凰大参考json
{"query":"CoREAS cosmic ray air shower radio emission", "max_results": 10}
CoREAS is a Monte Carlo code for the simulation of radio emission from extensive air showers, implemented directly in CORSIKA and based on the endpoint formalism for radiation from moving charges [1301.2132]. In this formulation, the radio signal is computed from the full time-dependent radiation field of the electrons and positrons tracked in the air shower, without assumptions on the emission mechanism and without intermediate histogramming of the shower particles [1301.2132]. Across the literature, CoREAS is described as parameter-free, first-principles, and microscopic; it has been used to model MHz-to-GHz emission, to compare simulations with measurements from arrays such as LOPES and LOFAR, and to support radio-based reconstruction of shower energy and shower maximum \(X_{\max}\) [1307.7566].

## 1. Origins, scope, and integration into CORSIKA

CoREAS was introduced as the successor to REAS 3.11 and as the end-point-formalism implementation of the radio-emission module integrated into CORSIKA [1309.5920]. Its defining architectural feature is that the radio calculation is performed during the air-shower simulation itself rather than as an afterburner using histogrammed particle distributions [1301.2132]. CORSIKA propagates charged particles in small steps, and CoREAS computes the radio contribution from the endpoints of each track segment and accumulates the electric field at user-specified observer positions [1301.2132].

This integration has two immediate consequences. First, CoREAS preserves correlations in the simulated particle distributions, including energy, angle, and lateral position, because no intermediate histogramming is performed [1301.2132]. Second, it incorporates the full electromagnetic shower development alongside the radio calculation, including the effects of geomagnetic deflection, charge excess, and atmospheric refractivity [1507.07389]. In the 2013 description, benchmarks indicated a speed-up of up to a factor of ten over the predecessor REAS3.1 [1301.2132].

Later work extended the CoREAS ecosystem beyond the original FORTRAN-based CORSIKA 7 implementation. In CORSIKA 8, the radio code was re-implemented as a modular observer in a C++17 framework designed to support both the CoREAS and ZHS formalisms [2112.11761]. That re-write preserved the endpoint-based physics while introducing a plugin architecture, HDF5/ROOT output, automated validation, and a path toward GPU acceleration [2112.11761]. This suggests that CoREAS has evolved from a specialized CORSIKA extension into a more general radio-emission framework embedded in modern air-shower simulation infrastructure.

## 2. Electrodynamic formulation and emission content

The central theoretical ingredient of CoREAS is the endpoint formalism, derived from the Liénard–Wiechert fields [1301.2132]. In the time-domain presentation, the vector potential for a point charge \(q\) with trajectory \(\mathbf r(t)\) at observer position \(\mathbf x\) is written as
\[
\mathbf A(\mathbf x,t)
=
\frac{\mu_0 q}{4\pi}
\frac{\boldsymbol \beta(t_{\mathrm{ret}})}
{R(t_{\mathrm{ret}})\left(1-\mathbf n\cdot \boldsymbol\beta(t_{\mathrm{ret}})\right)},
\]
with the electric field obtained by
\[
\mathbf E(\mathbf x,t)=-\frac{\partial \mathbf A}{\partial t}.
\]
After Fourier transformation and discretization into track endpoints, CoREAS evaluates a sum over segment endpoints rather than a continuous acceleration integral [1301.2132].

In practical terms, each charged-particle step is represented as a straight segment with a start and an end. The radiation is then computed from the velocity change at those endpoints, and the total field is the coherent superposition of all endpoint contributions from all shower particles [1301.2132]. Because the method is based on the exact particle tracks provided by CORSIKA, CoREAS does not need to impose separate sub-models for different emission channels [1301.2132].

Within this framework, the two dominant radio-emission mechanisms emerge automatically. The geomagnetic contribution arises from the transverse acceleration of shower electrons and positrons in the geomagnetic field and is dominant [1309.5920]. The Askaryan, or charge-excess, contribution arises from the net negative charge buildup in the shower front and is subdominant but non-negligible [1309.5920]. For LOPES energies and geomagnetic field strength, the Askaryan contribution was summarized as approximately \(5\)–\(10\%\) of the geomagnetic contribution [1309.5920]. A South Pole simulation study expressed these components in shower-plane coordinates as
\[
G \equiv |E_x| \approx |E_{\rm geo}|,\qquad
A \equiv |E_y| \approx |E_{\rm Ask}|,
\]
along the \(\mathbf v\times(\mathbf v\times\mathbf B)\) axis, and defined both the amplitude ratio
\[
R \equiv \frac{A}{G}
\]
and the Askaryan fraction
\[
F \equiv \frac{A}{A+G}
\]
for characterizing their relative importance [2108.06336].

A key consequence of the endpoint treatment is that no explicit choice between “geomagnetic,” “geosynchrotron,” and “charge-excess” descriptions is required at the algorithmic level [1301.2132]. The literature repeatedly emphasizes that both the dominant geomagnetic signal and the charge-excess signal are recovered from the same microscopic electrodynamics [1301.2132, 1307.7566]. At GHz frequencies, however, a 2013 study noted that the observed “clover-leaf” polarization structure might indicate that a direct geosynchrotron acceleration component becomes non-negligible at the shortest wavelengths, while also stressing that this point remained under study [1301.2132]. That constitutes a clear example where the simulations reproduce a pattern but the physical interpretation is treated cautiously.

## 3. Atmospheric refractivity, coherence, and observable structure

A defining physical feature of CoREAS is its treatment of the atmospheric refractive index. In early descriptions, CoREAS used a simple exponential parameterization,
\[
n(h)=1+(n_0-1)\exp(-h/h_0),
\]
with \(n_0-1 \simeq 2.92\times 10^{-4}\) at sea level and scale height \(h_0 \simeq 8\ {\rm km}\) [1309.5920]. In another presentation, this appears as
\[
n(h)=1+\delta n_0 \exp(-h/H),
\]
with \(H\approx 8\,\mathrm{km}\) and \(\delta n_0 \simeq 3\times 10^{-4}\) [1301.2132]. The refractive index enters both the retarded-time calculation and the coherence condition.

The Cherenkov-like enhancement follows from the denominator \(1-\mathbf n\cdot \boldsymbol\beta\) in the endpoint expression approaching zero near the Cherenkov angle [1301.2132]. In simplified form, the local Cherenkov condition is written as
\[
\cos\theta_C(h)=\frac{1}{\beta\,n(h)},
\]
or, for \(\beta \simeq 1\),
\[
\cos\theta_C \simeq \frac{1}{n}
\]
[1301.2132, 1309.5920]. This time compression is responsible for the well-known “bump” around an axis distance of approximately \(120\ {\rm m}\) in the lateral distribution and for rising lateral profiles in some LOPES events [1309.5920]. The 2013 LOPES comparison stated explicitly that only simulations including the refractive index can reproduce rising lateral distributions [1309.5920].

At MHz frequencies, CoREAS predicts the familiar asymmetric footprint produced by the vector superposition of geomagnetic and Askaryan emission [1301.2132]. In the 43–74 MHz band, individual showers were described by a two-parameter exponential lateral distribution,
\[
\epsilon(r)=\epsilon_{100}\exp[-\eta(r-100\,{\rm m})],
\]
where \(\epsilon_{100}\) is the electric-field amplitude at \(100\ {\rm m}\) and \(\eta\) the slope parameter [1309.5920]. In that same study, \(\epsilon_{100}\) was reported to scale nearly linearly with primary energy, \(\epsilon_{100}\propto E^p\) with \(p=0.95\ldots1.05\), and the slope \(\eta\) was found to correlate with zenith angle and \(X_{\max}\) [1309.5920].

At higher frequencies, the refractive index has even stronger consequences. In the 300–1200 MHz range, CoREAS predicts Cherenkov rings with field strengths \(2\)–\(3\) times larger than off-ring values when the shower geometry places the Cherenkov cone on the ground [1301.2132]. In the 3.4–4.2 GHz band, the Cherenkov ring persists but the polarization pattern changes qualitatively: the north–south component exhibits a characteristic “clover-leaf” quadrupole structure [1301.2132].

The same atmospheric refractivity that shapes coherence also contributes materially to reconstruction systematics. A dedicated LOFAR study found that a realistic \(4\%\) variation in refractivity induces a systematic error in inferred \(X_{\max}\) between \(3.5\) and \(11\ \mathrm{g/cm^2}\) in the 30–80 MHz band for proton showers with zenith angles from \(15^\circ\) to \(50^\circ\), and between \(10\) and \(22\ \mathrm{g/cm^2}\) in the 120–250 MHz band [1701.07338]. A later implementation of time-dependent event-specific atmospheres using GDAS in CORSIKA and CoREAS showed that under normal conditions the systematic shift in reconstructed \(X_{\max}\) was about \(2\ \mathrm{g/cm^2}\), but under extreme weather conditions it could reach \(15\ \mathrm{g/cm^2}\) [2006.02228]. For highly inclined showers, refractivity also causes a displacement of the radio-emission footprint with respect to the Monte Carlo shower impact point, reaching approximately \(1500\ {\rm m}\) in the ground plane for \(85^\circ\) zenith [2005.06775].

## 4. Frequency coverage, wavefront structure, and geometry dependence

CoREAS has been used across a wide radio-frequency range, from tens of MHz to several GHz [1301.2132, 1307.7566]. In the tens-of-MHz regime relevant for LOPES, LOFAR, and AERA, the footprint is dominated by the interference of geomagnetic and Askaryan fields, with a pronounced east–west asymmetry [1301.2132]. In this band, proton and iron primaries produce different lateral-distribution slopes because iron showers develop higher and therefore yield flatter footprints [1301.2132].

The wavefront structure of the radio pulse is another observable directly studied with CoREAS. For LOPES conditions, the wavefront was found to be approximately hyperbolic rather than spherical or conical [1507.07753]. In the notation of that study, the pulse-arrival time relative to the shower-core arrival time can be fitted by
\[
c\,\tau(d,z_s)=\sqrt{(d\sin\rho)^2+(cb)^2}+z_s\cos\rho+cb,
\]
where \(d\) is the orthogonal distance to the shower axis, \(z_s\) is the orthogonal distance above or below the shower plane, \(\rho\) is the cone-angle parameter, and \(b\) is a time offset empirically found to be approximately \(-3\ {\rm ns}\) [1507.07753]. The steepness of this hyperbolic wavefront was shown to be sensitive to the shower maximum.

After correcting for zenith-angle dependence, the same study obtained the empirical relation
\[
X_{\max}\,[{\rm g/cm^2}] \simeq 2.5\times 10^4 \, (\rho/{\rm rad})\,\cos^{-1.5}\theta,
\]
with an \(X_{\max}\) resolution of approximately \(25\ \mathrm{g/cm^2}\) under ideal simulated conditions [1507.07753]. CoREAS also predicted a slight east–west asymmetry in the wavefront due to geomagnetic–Askaryan interference, but the arrival-time asymmetry was only at the \(\sim 1\ {\rm ns}\) level at \(r\approx 200\ {\rm m}\), much weaker than the amplitude asymmetry [1507.07753]. This suggests that symmetric hyperbolic fits are usually sufficient for practical reconstruction even though the fully microscopic calculation produces a weak asymmetry.

Shower geometry also strongly affects the size of the radio footprint. CoREAS simulations showed that the illuminated ground area increases strongly from near-vertical to very inclined showers [1307.7566]. Reported characteristic scales include footprint radii of about \(100\)–\(200\ {\rm m}\) at \(10^{18}\ {\rm eV}\) for \(\theta=30^\circ\), a major-axis radius of about \(600\ {\rm m}\) for \(\theta=50^\circ\), and extensions to approximately \(3\ {\rm km}\) at \(10^{18}\ {\rm eV}\) and \(5\ {\rm km}\) at \(10^{19.5}\ {\rm eV}\) along the major axis for \(\theta=75^\circ\) [1307.7566]. These results underlie the frequent observation that very inclined air showers are especially promising for sparse radio arrays [1307.7566].

## 5. Comparison with measurements and resolution of the amplitude-scale tension

The most prominent experimental controversy associated with CoREAS concerned the absolute amplitude scale in comparisons with LOPES. In the 2013 LOPES comparison, CoREAS amplitudes were found to be lower than the measured \(\epsilon_{100}\) by roughly a factor of two, with a mean deviation of \(-50\%\pm10\%\) even after allowing for the quoted systematic uncertainties of \(\pm 35\%\) in the LOPES calibration and \(\pm 20\%\) in the KASCADE-Grande energy scale [1309.5920]. By contrast, REAS 3.11 agreed with the then-current LOPES amplitude scale to within about \(10\%\) on average [1309.5920]. At that stage, the origin of the CoREAS-versus-REAS amplitude difference was explicitly described as not understood [1309.5920].

This tension was resolved in 2015 through a revised absolute amplitude calibration of LOPES. The original external reference source had been characterized under free-field conditions with a horizontal reflective ground plane, whereas cosmic-ray radio pulses required free-space calibration [1507.07389, 1508.03471]. After the manufacturer re-measured the same reference source under free-space conditions, all previously published LOPES electric-field amplitudes had to be divided by a factor
\[
C_{\rm scale}=2.6\pm0.2,
\qquad
E_{\rm new}=E_{\rm old}/2.6
\]
[1507.07389]. The remaining one-sigma uncertainty on the absolute scale became about \(\pm16\%\) [1507.07389].

With this recalibration, the event-by-event comparison to CoREAS moved from a factor-of-two discrepancy to consistency at the absolute amplitude level. One report quoted
\[
\langle \epsilon_{100,\rm meas}/\epsilon_{100,\rm CoREAS}\rangle = 1.02
\]
for proton primaries and \(0.91\) for iron primaries [1507.07389], while another quoted
\[
\langle R'\rangle =0.98 \ \text{for proton simulations},\qquad
\langle R'\rangle =1.09 \ \text{for iron simulations}
\]
with \(R'=\epsilon_{100}'/\epsilon_{100}^{\rm CoREAS}\) [1508.03471]. In the latter analysis, a histogram of \((\epsilon_{100}'-\epsilon_{100}^{\rm CoREAS})/\sigma_i\) had a Gaussian width of about unity and the \(\eta\) distributions matched nearly perfectly, leaving only a small zenith-angle dependence at large \(\theta\), plausibly due to imperfect antenna directivity modeling rather than a flaw in CoREAS [1508.03471].

The recalibration therefore changed the interpretation of CoREAS fundamentally. What had appeared to be a significant amplitude deficit became, after calibration correction, evidence that the absolute field strengths predicted by the current simulations agree with experimental data [1507.07389]. The same recalibrated comparison placed REAS 3.11 in tension with the measurements, requiring an ad hoc scaling factor \(f\approx 0.43\) to match the data [1507.07389]. This episode is central to the historiography of CoREAS because it transformed a perceived simulation deficiency into a calibration correction and reinforced the status of CoREAS as the more accurate of the two codes.

## 6. Parametrizations, computational surrogates, and calibration applications

Because full CoREAS simulations are computationally demanding, substantial work has focused on compressed parameterizations and surrogate models derived from CoREAS outputs. For LOFAR, a two-dimensional pulse-power footprint was parameterized with a double Gaussian requiring five fit parameters [1402.2872]. Those parameters were found to correlate strongly with shower energy, arrival direction, and \(X_{\max}\), and the parameterization reproduced CoREAS power footprints at the \(\sim 10\)–\(20\%\) level while also fitting LOFAR measurements with reduced \(\chi^2 \approx 1.3\)–\(1.6\) [1402.2872].

A more explicitly simulation-accelerating approach is the template-synthesis method based on the universality of emission from small shower slices [1908.09543]. In this framework, a high-statistics CoREAS shower is decomposed into slant-depth slices, and the synthesized field for a target shower is formed as
\[
E^{\rm Synth}_{\rm simple}(r,t)=\sum_X
\left[\frac{N^{\rm Real}(X)}{N^{\rm Temp}(X)}\right]
E^{\rm Temp}_{\rm slice}(X,r,t),
\]
with an extended frequency-domain formulation that also rescales the slice spectrum as a function of \(X_{\max}\) [1908.09543]. The advanced synthesis reduced residuals to \(\lesssim 3\%\) at all distances and both polarizations for the ensemble studied, while reducing run times from \(O(10^4\,{\rm CPU\,s})\) per shower to \(O(10\,{\rm CPU\,s})\) [1908.09543]. This suggests that CoREAS has not only served as a forward model, but also as the training ground for fast emulators usable in dense-array reconstruction.

CoREAS has also been central to radio calorimetry. A dedicated study of the radiation energy in the 30–80 MHz band derived an efficient one-axis integration method and found that, after corrections for geomagnetic geometry and atmospheric density, the corrected radiation energy scales nearly quadratically with the electromagnetic shower energy [1606.01641]. The fitted relation
\[
S_{\rm RD}^{\rho}=A\left(\frac{E_{\rm em}}{10^{18}\,{\rm eV}}\right)^B\times 10^7\,{\rm eV}
\]
yielded \(A=1.683\pm0.004\) and \(B=2.006\pm0.001\), with an intrinsic uncertainty of \(4\%\) in determining the electromagnetic energy when \(X_{\max}\) is not accessible [1606.01641]. In a later Auger Engineering Radio Array analysis using 844 events and per-event CoREAS simulations, the ratio of radio to fluorescence-detector energy scales was measured as
\[
\frac{E_{\rm radio\,scale}}{E_{\rm FD\,scale}}
=
1.120^{+0.009}_{-0.008}\,(\mathrm{stat})
\pm 0.173\,(\mathrm{sys}),
\]
corresponding to radio energies \(12\%\) higher than the fluorescence-detector scale but consistent within uncertainties [2507.08500].

These developments indicate that CoREAS has become more than a simulation code for waveform prediction. It underpins energy-scale determinations, atmosphere-aware \(X_{\max}\) reconstruction, parameterized footprint models, and fast synthesis schemes [1606.01641, 2006.02228, 1908.09543, 2507.08500]. A plausible implication is that the scientific role of CoREAS has broadened from validating the physics of radio emission to providing a reference model for precision calibration and high-throughput inference in radio detection of cosmic rays.

Source: https://www.emergentmind.com/topics/coreas