🔭 astrophysics

Revisiting the Radio Lateral Distribution Function: An amplitude dependence on $X_{\rm max}$ and primary composition

This paper demonstrates that the amplitude of the radio lateral distribution function at ground level exhibits a strong dependence on the atmospheric depth of the shower maximum ( $X_{\rm max}$ ) due to competing density and geometric scaling effects, thereby encoding primary cosmic ray composition information directly in the signal amplitude and offering a potential method for event-by-event composition inference without explicit $X_{\rm max}$ reconstruction.

Original authors: Washington R. Carvalho Jr., Lech Wiktor Piotrowski

Published 2026-03-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Washington R. Carvalho Jr., Lech Wiktor Piotrowski

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is constantly raining down invisible particles called cosmic rays. When one of these high-speed particles hits our atmosphere, it doesn't just stop; it explodes into a massive, cascading shower of billions of smaller particles. Scientists call this an Extensive Air Shower.

Think of this shower like a giant, invisible snowball rolling down a mountain. As it rolls, it picks up more snow (particles) and grows bigger until it reaches its largest size, then it starts to shrink. The point where the snowball is biggest is called $X_{max}$ (Shower Maximum).

For decades, scientists have tried to figure out what the original "snowball" was made of. Was it a light, fluffy piece of ice (a proton)? Or was it a heavy, dense rock (an iron nucleus)?

The Old Way: Looking at the Shape

Traditionally, scientists used radio antennas on the ground to listen to the "crackling" sound these showers make. They realized that the shape of the sound pattern on the ground changes depending on whether the shower is made of light or heavy stuff. It's like looking at the shadow of a tree: a tall, thin tree casts a different shadow than a short, bushy one. By measuring the shape of the radio "shadow," they could guess the composition.

The New Discovery: Listening to the Volume

This paper reveals a secret that was hiding in plain sight. The authors discovered that it's not just the shape of the radio signal that matters; it's also the loudness (amplitude).

Here is the twist: The loudness of the radio signal depends heavily on where the shower reached its maximum size ( $X_{max}$ ), even if we account for how much energy the shower had.

To understand why, the authors propose a battle between two opposing forces, like two people pulling on a rope:

1. The "Closer is Louder" Effect (Distance Scaling)

Imagine you are standing in a field watching a fireworks display. If the fireworks explode right above your head, they look and sound huge. If they explode high up in the sky, they look smaller and quieter.

The Physics: Lighter particles (protons) tend to penetrate deeper into the atmosphere before exploding. This means the "explosion" happens closer to the ground. Because the radio signal has less distance to travel, it arrives louder.
The Metaphor: This is the $1/R$ rule. The closer the source ( $R$ ), the louder the signal.

2. The "Thin Air is Quiet" Effect (Density Scaling)

Now, imagine the fireworks are set off in a thick fog versus a thin, high-altitude mist. In the thick fog, the sound waves bounce around and stay strong. In the thin mist, the sound dissipates quickly.

The Physics: As the shower goes deeper, the air gets thicker (denser). In thick air, the charged particles in the shower get "dragged" by collisions with air molecules. This drag actually helps them stay in sync, creating a stronger radio signal. In thin air (high up), they fly wilder and get out of sync, weakening the signal.
The Metaphor: This is the $1/\rho$ rule. The denser the air ( $\rho$ ), the stronger the signal.

The Great Tug-of-War

The paper explains that the final volume of the radio signal is the result of these two forces fighting each other:

Force A (Distance): Wants the signal to be loud because the shower is closer to the ground.
Force B (Density): Wants the signal to be quiet because the shower is in thinner air (if it's high up) or loud because it's in thick air (if it's low down).

The Winner Depends on the Angle:

At low angles (looking straight down): The "Closer is Louder" force wins. Since protons go deeper (closer), they sound louder than iron showers.
At high angles (looking sideways): The "Thin Air" effect gets complicated. The air is so thin that the particles get so wild they lose their "coherence" (they stop marching in step). This "loss of coherence" acts like a volume knob turning down the signal. In some cases, this effect is so strong that it flips the script, making the "heavy" iron showers sound louder than the "light" proton ones.

Why This Matters

The authors realized that this "volume" difference is a goldmine of information.

Direct Composition Detection: You don't even need to calculate exactly where the shower peaked ( $X_{max}$ ) to guess what it's made of. You can just listen to the volume of the radio signal. If it's loud, it's likely a light proton; if it's quiet, it's likely a heavy iron nucleus.
Better Energy Estimates: Scientists use these radio signals to calculate the energy of the cosmic ray. If they ignore this "volume vs. depth" relationship, they might get the energy wrong. It's like trying to guess how loud a speaker is without knowing if it's close to you or far away.
The "Magic Angle": There is a specific angle where the shape of the signal looks identical for all showers, making it impossible to tell them apart by shape. But at this exact angle, the volume still changes! This means the "loudness" is the only clue left to solve the mystery.

The Bottom Line

This paper is like finding a new instrument in an orchestra. For years, scientists were only listening to the melody (the shape of the signal) to identify the instrument. This research shows that the volume (the amplitude) is also a distinct fingerprint that tells us exactly what the instrument is. By understanding the tug-of-war between distance and air density, we can finally hear the true composition of the universe's most energetic particles.

1. Problem Statement

The detection of Ultra-High Energy Cosmic Rays (UHECRs) via radio emission from Extensive Air Showers (EAS) is a mature technique. While it is well-established that the shape of the radio Lateral Distribution Function (LDF) depends on the depth of the shower maximum ( $X_{max}$ ) and primary mass composition (due to geometric projection effects of the Cherenkov cone), the dependence of the electric field amplitude on $X_{max}$ has been historically overlooked or attributed solely to differences in the electromagnetic (EM) energy of the showers.

Previous reconstruction methods (e.g., LOFAR-like $\chi^2$ comparisons) often treat amplitude variations as systematic errors to be corrected or normalized away. However, the authors observed that even after normalizing for EM energy, a strong, systematic dependence of the peak radio amplitude on $X_{max}$ (and thus primary composition) persists. The paper aims to:

Quantify this residual amplitude dependence.
Provide a complete physical explanation based on competing scaling laws.
Demonstrate how this dependence can be exploited for event-by-event composition discrimination without explicit $X_{max}$ reconstruction.

2. Methodology

The authors utilized a combination of machine learning analysis, full-scale simulations, and semi-analytical modeling:

Simulations: They employed the ZHAireS simulation package (combining AIRES for shower development and ZHS for radio emission calculation).
- Sites: Simulations were performed for the Giant Radio Array for Neutrino Detection (GRAND) site (high geomagnetic field: 56.4 $\mu$ T) and the Pierre Auger Observatory (lower field: 24.0 $\mu$ T).
- Parameters: Proton and Iron primaries were simulated at various zenith angles ( $42^\circ$ – $85^\circ$ ) and energies (1.25 EeV for GRAND, 5 EeV for Auger).
- Normalization: To isolate the $X_{max}$ effect from energy effects, all simulated electric fields were normalized by the Electromagnetic (EM) energy of the shower, not the primary energy. This ensures that any remaining amplitude differences are due to shower development geometry and density effects, not intrinsic energy variations.
Machine Learning (Motivation): The study was motivated by a Random Forest classifier that achieved ~82% accuracy in distinguishing protons from iron using only antenna distance and peak amplitude, bypassing $X_{max}$ reconstruction. Feature importance analysis revealed that the amplitude of the closest antenna was the most critical discriminator.
Analytical Modeling: The authors derived scaling laws for the electric field based on the ZHS formalism, considering:
- Distance scaling ( $1/R$ ).
- Air density scaling ( $1/\rho$ ).
- Coherence loss effects at low densities/high magnetic fields.

3. Key Contributions and Physical Mechanism

The core contribution is the identification and explanation of two competing scalings that govern the radio amplitude at ground level:

Distance Scaling ( $1/R$ ):
- $R$ is the distance from the shower maximum ( $X_{max}$ ) to the ground observer.
- Deeper showers (higher $X_{max}$ , typical of protons) have a smaller $R$ .
- Effect: As $X_{max}$ increases, $R$ decreases, leading to a stronger electric field amplitude. This favors lighter primaries (protons).
Density Scaling ( $1/\rho$ ):
- $\rho$ is the air density at $X_{max}$ .
- Deeper showers occur in denser air.
- Effect: Higher density reduces the drift velocity of charged particles, decreasing the transverse current and thus the geomagnetic emission. As $X_{max}$ increases, $\rho$ increases, leading to a weaker electric field amplitude. This favors heavier primaries (iron, which have shallower $X_{max}$ ).
Coherence Loss Factor ( $J(\theta)$ ):
- At high zenith angles (low densities) and high geomagnetic fields, the Lorentz force causes large particle deflections.
- This introduces time delays between particle tracks, causing a loss of coherence in the radio emission.
- The authors introduce a non-linearity factor $J(\theta)$ to modify the density scaling to $(1/\rho)J(\theta)$ .
- Crucial Finding: At very high zenith angles and high magnetic fields (e.g., GRAND), coherence loss becomes so severe that it can reverse the density scaling trend, effectively suppressing the amplitude of deeper showers even further.

Net Result:

Low Zenith Angles: The $1/R$ scaling dominates. Protons (deeper showers) have higher amplitudes than Iron.
High Zenith Angles (Auger): The $1/\rho$ scaling dominates. Iron (shallower showers) has higher amplitudes than Protons.
High Zenith Angles (GRAND): Due to the strong magnetic field, coherence loss ( $J$ ) suppresses the density scaling so heavily that the $1/R$ scaling remains dominant even at high zeniths, meaning Protons still have higher amplitudes.

4. Results

Amplitude vs. $X_{max}$ : The simulations confirm a strong, monotonic dependence of the peak LDF amplitude on $X_{max}$ even after EM energy normalization.
Site Dependence:
- At Auger (lower $B$ ), there is a transition angle (~ $72^\circ$ ) where the dominant scaling flips from $1/R$ (proton-favored) to $1/\rho$ (iron-favored).
- At GRAND (high $B$ ), the coherence loss is so significant that the $1/R$ scaling dominates across all zenith angles; protons consistently show higher amplitudes than iron.
Magic Angle: Near the "magic angle" (~ $82^\circ$ ), the shape of the LDF becomes nearly independent of $X_{max}$ , but the amplitude retains a strong dependence. This highlights that amplitude carries unique information not present in the footprint shape.
Machine Learning Validation: The physical explanation aligns with the Random Forest results: the classifier successfully used the amplitude difference (driven by $X_{max}$ ) to discriminate composition.

5. Significance and Implications

New Discrimination Method: The amplitude dependence allows for event-by-event primary composition discrimination without reconstructing $X_{max}$ first. This is particularly valuable for experiments like GRAND where $X_{max}$ reconstruction via footprint shape is difficult at high zeniths.
Refining $X_{max}$ Reconstruction: Current "black-box" $\chi^2$ methods (like LOFAR) often use a free energy scaling factor that normalizes away amplitude differences. The authors argue that fixing or constraining this scaling factor could improve $X_{max}$ resolution by utilizing both shape and amplitude information, especially near the magic angle.
Energy Estimation Bias: If radio fluence (used to estimate EM energy) depends on $X_{max}$ and composition, failing to account for this could introduce biases in energy reconstruction. The paper suggests that existing empirical corrections for distance and density in energy reconstruction methods are physically rooted in these competing scalings.
Site Selection: For experiments targeting nearly horizontal showers, the geomagnetic field strength is less critical than previously thought because coherence loss suppresses the field dependence at high zeniths. This opens up potential sites with lower magnetic fields for such experiments.

In summary, the paper provides a unified physical framework explaining why radio amplitude varies with shower depth, revealing that this variation is a robust observable for determining cosmic ray composition and refining reconstruction techniques.