Radio Morphing: Hybrid Sources & Air Showers

Updated 17 January 2026

Radio morphing is a phenomenon where changes in radio morphology arise from intrinsic and environmental variations, observed in both hybrid radio sources and cosmic air showers.
It employs analytical scaling laws, spatial transformations, and Fourier-domain interpolation to rapidly and accurately simulate electric field footprints.
The technique underpins the optimization of radio observatory arrays and bridges diverse applications from astrophysical synchrotron studies to particle-induced air shower modeling.

Radio morphing refers to a class of physical and computational phenomena in which the observed radio morphology of an astrophysical or particle-induced system changes as a result of variations in intrinsic or environmental parameters. The term is used in two distinct research contexts: (1) in extragalactic astrophysics, to describe hybrid morphology radio sources (HyMoRS) that exhibit a mixture of Fanaroff-Riley I and II radio morphologies, and (2) in astroparticle physics, to denote semi-analytical techniques for surrogate modeling of radio signals from extensive air showers (EAS), enabling rapid computation of simulated electric fields at arbitrary antenna positions. Both usages fundamentally exploit the concept of morphological transformation under systematic physical influences, but differ in scope, methodology, and application domain.

1. Astrophysical Radio Morphing: Hybrid Morphology Radio Sources (HyMoRS)

Radio morphing in extragalactic radio astronomy is exemplified by the discovery and analysis of hybrid morphology radio sources. HyMoRS present a morphological transition where one lobe displays Fanaroff-Riley II (FR II) morphology—well-collimated jets ending in hotspots with backflow—and the opposite lobe exhibits Fanaroff-Riley I (FR I) morphology—diffuse, turbulent plumes lacking terminal hotspots. RG1 is a prototypical case, with multi-band radio imaging and X-ray data confirming this asymmetry (Gasperin, 2017).

Quantitative morphological metrics for RG1 include a projected size of ≃840 kpc, distinct surface brightness profiles for each lobe, region-specific spectral indices, and polarization fractions. For example, the northern (FR II) lobe shows surface brightness jumps by a factor of three and compact peaks at the jet termination, whereas the southern (FR I) lobe displays a smooth brightness decrease and lacks clear terminal features. Polarization and spectral index measurements (e.g., α ≃ –1.20 near the FR II-side core and α ≃ –1.30 at the FR I-side terminus) further differentiate the lobes.

Multi-frequency and polarization analyses demonstrate that such radio morphing is not explainable by central engine properties alone; rather, it results from the interplay of intrinsic parameters (e.g., jet power, accretion regime) and strongly asymmetric environmental factors (e.g., local density gradients, medium-induced jet instabilities). Thus, HyMoRS serve as laboratories for understanding the FR I/II dichotomy as a product of radio morphing processes (Gasperin, 2017).

2. Fundamental Physical Mechanisms in Particle-Induced Air-Shower Radio Morphing

Radio morphing in astroparticle physics is formalized as a semi-analytical procedure for the rapid computation of radio emission footprints from cosmic-ray-induced extensive air showers. The dominant emission mechanisms for EAS radio signals are the geomagnetic effect (transverse current due to v × B drift of charged particles) and Askaryan (charge-excess) emission. The resulting electric field at an antenna is governed by:

Primary energy (linear scaling with particle number)
Geomagnetic angle (sin α dependence)
Air density at shower maximum, $ρ(X_{\rm max})$
Cherenkov-cone geometry (set by the local index of refraction, $n$ )
Observer position (lateral and longitudinal relative to shower axis)
Geomagnetic field vector ( $B$ )

These dependencies are encapsulated in scaling relations that factorize the electric field amplitude into multiplicative correction factors for energy, geomagnetic angle, density, and Cherenkov stretch (Zilles et al., 2018, Zilles et al., 2018, Zilles et al., 2018).

3. Radio Morphing Methodologies in Air-Shower Simulations

The computational radio morphing technique (hereafter RM: Editor's term) is predicated on the universality of air-shower radio footprints and reduces required Monte Carlo simulations by orders of magnitude. The workflow consists of three main steps:

Amplitude scaling: The electric field traces from a reference simulation are scaled according to analytic laws—a combination of $E/E_{\rm ref}$ (energy), $(\sin\alpha)/(\sin\alpha_{\rm ref})$ (geomagnetic), $[\rho_{\rm ref}/\rho]^{½}$ (air density), and Cherenkov-cone kinematics ( $\Theta_C$ ). For example:

$\vec E_B(\vec x_B, t) = j_E\,j_\rho\,j_C\,J_\alpha\,\vec E_A(j_C^{-1}\vec x_A, t)$

where each $j$ factor is defined by the respective physical scaling (Zilles et al., 2018, Zilles et al., 2018).

Spatial transformation (isometry): Reference antenna positions are rotated, scaled (Cherenkov-cone stretch), and translated to coincide with the geometry of the target shower. The transformation is accomplished via rotation matrices constructed from the reference and target shower axes.
Fourier-domain interpolation: The scaled and transformed reference traces are interpolated to the desired antenna grid. This is conducted in the frequency domain for both amplitude and phase, using barycentric or linear interpolation between the nearest reference points (Zilles et al., 2018, Chiche et al., 2022). For 3D topographies, bilinear in-plane and linear longitudinal interpolation are employed.

Quantitative performance, as validated against full ZHAireS/CoREAS simulations, includes amplitude reproduction with sub-10--25% deviations for the vast majority of antennas and signal timing accuracy to better than 1 ns. Table 1 summarizes key benchmark results:

Metric	RM (2026, $\theta>60^\circ$ ) (Chiche et al., 10 Jan 2026)	RM (2022, all zeniths) (Chiche et al., 2022)
Peak amplitude bias (raw)	≤ 0%	<10% (91% of antennas)
RMS deviation (full band)	≤17%	<25% (99% of antennas)
RMS deviation ([30–80] MHz)	≤10%	<10% (91% of antennas)
Timing accuracy	>99% antennas <5 ns	RMS ≲1 ns
CPU speed-up vs full MC	$>10^4$ (RM: 10–30 s/shower; ZHS: days)	$>10^2$ (RM: 1 s vs 10,000 s)

4. Upgrades in Scaling, Interpolation, and Fluctuation Modeling

Recent developments incorporate refined scaling functions for air density, improved decomposition of geomagnetic and charge-excess components, and robust stochastic modeling of shower-to-shower fluctuations in $X_{\rm max}$ and total particle number. The empirical functions $f_{\rm geo}(\rho)$ and $f_{\rm ce}(\rho)$ , derived from ZHAireS fits to 11,000 showers, correct the amplitude for atmospheric density variations (Chiche et al., 2022, Chiche et al., 10 Jan 2026). Timing fidelity is enhanced by explicit correction for geometrical propagation delays.

Fluctuations in $X_{\rm max}$ and shower particle number are included via parameterizations of their mean and variance as a function of primary energy and mass, with corresponding adjustments in electric field via Gaussian sampling. This approach allows RM to realistically reproduce statistical distributions over large simulation campaigns.

5. Applications, Performance Limits, and Future Extensions

Radio morphing is deployed in the simulation chains of next-generation radio arrays (GRAND, AERA, LOFAR, SKA-Radio, RNO-G, IceCube-Gen2-Radio) for sensitivity mapping, detector optimization, and algorithm testing (Chiche et al., 10 Jan 2026, Chiche et al., 2022, Zilles et al., 2018). The method's computational cost scales linearly with the number of target antennas and reference showers, enabling full-coverage studies (≥ $10^6$ antennas, ≥ $10^4$ showers) that are otherwise intractable.

Known limitations include decreased accuracy near the Cherenkov ring (where field gradients are steep), for parameters outside the covered space of reference showers, and a current restriction to $\theta>60^\circ$ for inclined-shower applications (Chiche et al., 10 Jan 2026). Second-order corrections—such as improved scaling for geomagnetic field moduli and inclusion of full Askaryan-geomagnetic interference structure—are under development. Planned extensions include generalization to dense media (ice, sand) and local B-field adjustments, further enhancing RM's universality (Zilles et al., 2018, Chiche et al., 2022).

6. Radio Morphing in Extragalactic Context: Connection to AGN Jet Studies

The terminology "radio morphing" in the context of extragalactic radio sources provides an interpretive bridge between particle-induced and astrophysical synchrotron phenomena. In HyMoRS, the transition in lobe morphology is directly associated with environmental or intrinsic triggers—such as jet deceleration due to ambient medium density discontinuities or the growth of helical instabilities—which lead to a transformation of jet collimation and energetics. This process is accompanied by observable shifts in surface brightness, spectral index, polarization fraction, and orientation of the magnetic field vectors (Gasperin, 2017).

Such hybrid sources empirically demonstrate that large-scale radio morphologies are not solely determined by central power or accretion regime but are morphed by environmental asymmetries, aligning the conceptually distinct usages of the term across astrophysics and astroparticle physics.

7. Concluding Remarks

Radio morphing encapsulates both a physical process (the transformation of radio morphology in extragalactic and galactic jet systems under environmental gradients) and a computational methodology (the high-fidelity, high-speed prediction of air-shower radio emission using reference simulations and analytic scalings) (Gasperin, 2017, Chiche et al., 2022, Chiche et al., 10 Jan 2026). The success of RM in both scientific and computational domains is predicated on the approximate universality of the underlying radio emission mechanisms and the validity of scaling laws governing their dependence on primary, environmental, and observational parameters. Ongoing research seeks to refine these techniques, expand their domain of applicability, and exploit them for both fundamental physics discovery and technological innovation in radio observatory design.