Zooming-in on cluster radio relics -- I. How density fluctuations explain the Mach number discrepancy, microgauss magnetic fields, and spectral index variations

It is generally accepted that radio relics are the result of synchrotron emission from shock-accelerated electrons. Current models, however, are still unable to explain several aspects of their formation. In this paper, we focus on three outstanding problems: i) Mach number estimates derived from radio data do not agree with those derived from X-ray data, ii) cooling length arguments imply a magnetic field that is at least an order of magnitude larger than the surrounding intracluster medium (ICM), and iii) spectral index variations do not agree with standard cooling models. To solve these problems, we first identify typical shock conditions in cosmological simulations, using the results to inform significantly higher resolution shock-tube simulations. We apply the cosmic ray electron spectra code CREST and the emission code CRAYON+ to these, thereby generating mock observables ab-initio. We identify that upon running into an accretion shock, merger shocks generate a shock-compressed sheet, which, in turn, runs into upstream density fluctuations in pressure equilibrium. This mechanism directly gives rise to solutions to the three problems: it creates a distribution of Mach numbers at the shock-front, which flattens cosmic ray electron spectra, thereby biasing radio-derived Mach number estimates to higher values. We show that this effect is particularly strong in weaker shocks. Secondly, the density sheet becomes Rayleigh-Taylor unstable at the contact discontinuity, causing turbulence and additional compression downstream. This amplifies the magnetic field from ICM-like conditions up to microgauss levels. We show that synchrotron-based measurements are strongly biased by the tail of the distribution here too. Finally, the same instability also breaks the common assumption that matter is advected at the post-shock velocity downstream, thus invalidating laminar-flow based cooling models.

💡 Research Summary

This paper tackles three long‑standing puzzles of galaxy‑cluster radio relics—(i) the systematic discrepancy between Mach numbers derived from radio spectral slopes and those obtained from X‑ray temperature/density jumps, (ii) the need for micro‑gauss magnetic fields in relics despite the surrounding intracluster medium (ICM) typically hosting fields an order of magnitude weaker, and (iii) the failure of standard synchrotron cooling models (KP, JP, CI, KGJP) to reproduce the observed spectral‑index variations and colour‑colour tracks.
The authors adopt a hybrid approach. First, they analyse a cosmological MHD simulation (using the moving‑mesh code AREPO) to extract typical merger‑shock conditions (Mach numbers ≲3, upstream densities, temperatures, magnetic fields). These conditions inform a suite of high‑resolution shock‑tube experiments that include realistic upstream density fluctuations (δρ/ρ≈0.1–0.3) inferred from X‑ray observations of cluster outskirts.
In the shock‑tube runs, the merger shock compresses the upstream gas into a thin, high‑density sheet. When this sheet encounters the pre‑existing density fluctuations, the shock front becomes corrugated, producing a distribution of local Mach numbers rather than a single value. The authors show that this distribution is especially broad for weak shocks (M≲2). Because diffusive shock acceleration (DSA) predicts a flatter electron spectrum for higher Mach numbers, the radio emission—being weighted toward the brightest, highest‑Mach patches—systematically overestimates the Mach number relative to the X‑ray average. This naturally explains the radio‑X‑ray Mach‑number discrepancy without invoking projection effects.
The dense sheet also becomes Rayleigh‑Taylor (RT) unstable at the contact discontinuity. The RT instability drives turbulence and additional compression downstream, which stretches and winds magnetic field lines. The MHD calculations demonstrate that the magnetic field can be amplified from typical ICM values (~0.1 µG) to the observed micro‑gauss levels (∼1–5 µG) within a few hundred kiloparsecs of the shock. The amplified field further boosts synchrotron emissivity, helping to reconcile the high radio luminosities of relics.
Crucially, the RT‑driven turbulence breaks the assumption of a laminar post‑shock flow. In standard cooling models, the distance from the shock front is taken as a proxy for electron age (t_cool). The simulations reveal that downstream advection speeds become highly variable and that electrons are mixed and re‑accelerated by turbulent eddies. Consequently, the simple distance‑time mapping fails, and the classic KP/JP/CI/KGJP spectral shapes cannot reproduce the observed colour‑colour diagrams. By post‑processing the simulations with the cosmic‑ray electron spectral solver CREST and the synchrotron emission code CRAYON+, the authors generate mock radio maps that successfully match the observed curvature in colour‑colour space.
The paper’s results are summarized as follows:

1. Mach‑number bias: A Mach‑number distribution at the shock front leads to radio‑derived Mach numbers that are 20–30 % higher than the true average, especially for weak shocks.
2. Magnetic‑field amplification: RT instability amplifies magnetic fields by an order of magnitude, reaching µG levels consistent with Faraday‑rotation and relic‑width estimates.
3. Cooling model breakdown: Turbulent downstream flow invalidates the simple advection‑cooling relation, explaining why standard cooling models cannot fit observed spectral‑index gradients.
4. Improved synthetic observables: The combined CREST+CRAYON+ pipeline reproduces both the radio brightness and the detailed spectral‑index behaviour of observed relics.
   The authors acknowledge limitations: the initial spectrum of upstream density fluctuations is not directly constrained, particle‑in‑cell (PIC) physics of electron pre‑acceleration is not included, and the simulations are limited to 2‑D/3‑D with finite resolution, possibly under‑resolving the smallest turbulent scales and reconnection events. Nonetheless, the study provides a unified physical framework—density fluctuations plus RT instability—that simultaneously addresses the three major relic puzzles, offering a compelling direction for future high‑resolution, fully kinetic cluster simulations.