pith. machine review for the scientific record. sign in

arxiv: 2604.06302 · v1 · submitted 2026-04-07 · 🌌 astro-ph.HE · astro-ph.CO

Recognition: 2 theorem links

· Lean Theorem

Zooming in on radio relics -- II. How relic morphology probes density fluctuations at the edge of galaxy clusters

Joseph Whittingham , Christoph Pfrommer , Maria Werhahn , L\'ena Jlassi , Philipp Girichidis
Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:57 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.CO
keywords radio relicsgalaxy clustersdensity fluctuationspower spectrumshock morphologyintracluster mediumsynchrotron emissioncosmic ray electrons
0
0 comments X

The pith

Radio relic morphology encodes the power spectrum of density fluctuations in the outer intracluster medium.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Radio relics form at merger shocks where electrons are accelerated to produce synchrotron emission. The authors run shock-tube simulations that introduce controlled upstream density fluctuations with varying coherence length, amplitude, and power-law slope, then post-process the results with the Crest solver to generate mock radio relics. They show that distinct morphological features respond separately to each parameter: double-strand structures trace the coherence length through projected shock curvature, relic width and patchiness grow with fluctuation amplitude, and filament number, spacing, and curvature increase with steeper spectral slopes. This link turns observed relic shapes into a probe of gas properties at large cluster radii where X-ray and other methods lose sensitivity. Application to the Toothbrush and Sausage relics yields a typical coherence length near 500 kpc, non-uniform amplitudes, and spectra steeper than Kolmogorov.

Core claim

By systematically varying the upstream density power spectrum in cosmologically motivated shock-tube simulations and post-processing with the Crest cosmic ray electron solver, the authors demonstrate that the downstream radio morphology independently encodes the coherence length, amplitude, and spectral slope of density fluctuations. Curved shock fronts in projection produce double strands whose scale maps directly to the fluctuation coherence length. Higher amplitudes lengthen the downstream extent of the relic, broaden the Mach number distribution, raise integrated flux, and create patchier emission. Steeper slopes make double strands more likely, increase the number and spacing of shock-t

What carries the argument

The mapping from upstream density power spectrum parameters (coherence length, amplitude, power-law slope) to distinct features in downstream radio relic morphology after cosmic-ray acceleration and synchrotron modeling.

If this is right

  • Double-strand features in relics arise from curved shocks in projection and their spatial scale directly indicates the coherence length of upstream density fluctuations.
  • Higher fluctuation amplitudes produce relics that extend farther downstream than idealised thin-shock models predict and create patchier emission through a broader Mach number distribution.
  • Steeper power-law slopes increase the number and spacing of radio filaments oriented parallel to the shock front.
  • At higher Mach numbers, steeper slopes can generate curved radio filaments that trace large eddies.
  • Analysis of the Toothbrush and Sausage relics indicates a typical fluctuation coherence length of approximately 500 kpc, non-uniform amplitudes, and power spectra steeper than Kolmogorov scaling.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Radio relic morphology could provide constraints on outer ICM turbulence without needing direct density measurements at large radii.
  • Combining morphological analysis with X-ray or Sunyaev-Zel'dovich data might help separate the effects of coherence length from amplitude.
  • The same simulation framework could be applied to full cosmological runs to predict how relic properties evolve with cluster mass or redshift.
  • Future high-resolution radio observations could test whether real relics exhibit the predicted trends in filament spacing and curvature.

Load-bearing premise

The assumption that the Crest post-processing accurately models cosmic ray electron acceleration and synchrotron emission without significant contributions from unmodeled effects like magnetic field structure or additional microphysics.

What would settle it

High-resolution radio maps of a relic whose filament scales or widths fail to match the coherence length and amplitude inferred from independent density or turbulence measurements at the same shock location.

Figures

Figures reproduced from arXiv: 2604.06302 by Christoph Pfrommer, Joseph Whittingham, L\'ena Jlassi, Maria Werhahn, Philipp Girichidis.

Figure 2
Figure 2. Figure 2: Histograms of the electron number density, calculated using all gas cells in our initial high-resolution upstream region. Each distribution is log-normal and has a mean density set to 3.5×10−5 cm−3 . The standard value employed in our simulations is 𝜎/𝜇 = 0.4. in both M = 2 and M = 3 configurations. To do so, we set the initial downstream pressure to be 𝑃2 = 1 × 10−12 dyne cm−2 and 𝑃2 = 2.32 × 10−12 dyne c… view at source ↗
Figure 3
Figure 3. Figure 3: Slices through our M = 3 simulations at 𝑡 = 244 Myr, where the shock travels from left to right and colours indicate gas density. Dashed regions indicate secondary instabilities (see text). In each column, we vary only one characteristic of the upstream initial conditions (see [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left panel: Shock-dissipated energy in our M = 3, 𝐿coh = 75 kpc simulation at 𝑡 = 190 Myr, projected over a depth of 150 kpc. Green arrows indicate the rough position of ‘knots’, where the morphology narrows. The coherence length of the upstream density fluctuations is indicated by the vertical barred line. Middle panel: As previous, but for our 𝐿coh = 150 kpc simulation. Both simulations otherwise have 𝜎/… view at source ↗
Figure 5
Figure 5. Figure 5: Left panel: Slice showing gas temperature in our M = 2, high fluctuation amplitude (𝜎/𝜇 = 0.8) simulation at 𝑡 = 170 Myr. Middle panel: As previous, but showing shock-frame gas speed divided by sound speed (see text). Right panel: Thin projection of the emission-weighted Mach number with depth 35 kpc. Grey indicates no shock found. In our simulations, higher temperatures have higher sound speeds, which red… view at source ↗
Figure 6
Figure 6. Figure 6: Left-most panels: Synchrotron intensity maps at 150 MHz and spectral index maps taken between 1.5 GHz and 150 MHz for the Mach 2 high fluctuation amplitude (𝜎/𝜇) variations at 𝑡 = 244 Myr. The projection depth is 300 kpc. We have masked spectral index values where no emission takes place along the line of sight. Increasing 𝜎/𝜇 increases the extent of the downstream emission and the average intensity. A var… view at source ↗
Figure 7
Figure 7. Figure 7: As [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Top row: Synchrotron intensity maps at 1.5 GHz for our Mach 3 𝐿coh = 150 kpc simulation (with 𝛿 = −5/3) shown at 50 Myr intervals. Bottom row: as previous, but for our 𝛿 = −10/3 simulation. The dashed white box indicates a radio spur (see text). In contrast to the standard paradigm, in our scenario the radio relic only starts forming when the merger shock collides with an accretion shock. This leads to an … view at source ↗
Figure 9
Figure 9. Figure 9: Top left: Synchrotron intensity map at 1.5 GHz for our Mach 3 𝛿 = −10/3 simulation at 𝑡 = 150 Myr. Top right: As previous, but show￾ing intensity at 150 MHz. Bottom panels: Corresponding maps blurred with a Gaussian kernel to mimic current observational capabilities (see text). At lower frequencies, the downstream appears more extended, which reduces the coherence of ‘double strand’ features. We predict th… view at source ↗
read the original abstract

Gas properties in the outer intracluster medium (ICM) are not well-constrained, as traditional probes lose sensitivity at Mpc distances. We show that the morphology of radio relics effectively encodes the power spectrum of the surrounding density fluctuations, and that they hence represent a new observational window. To demonstrate this, we use cosmologically motivated shock-tube simulations in which we systematically vary the coherence length, amplitude, and power-law slope of the upstream density power spectra. We then post-process our simulations with the cosmic ray electron spectral solver, Crest, thereby producing a suite of mock radio relics. We find that the downstream morphology of our simulated relics is independently sensitive to each of the aforementioned parameters. Specifically, we show that observed 'double strand' features can be formed by curved shock fronts in projection, and that the scale of these features maps directly to the fluctuation coherence length. Increasing the fluctuation amplitude, meanwhile, progressively lengthens the downstream extent of the relic, thus explaining why relics are observed to be broader than the idealised expectation. It also broadens the Mach number distribution across the shock, which simultaneously increases the integrated radio flux density and produces patchier emission. Finally, steepening the power-law slope makes 'double strand' features more likely, and additionally increases both the number of radio filaments oriented parallel to the shock front and their spacing. At higher Mach numbers, steepening the power-law slope can further lead to the production of curved radio filaments, which trace large eddies. We apply our analysis to the Toothbrush and Sausage relics, and find evidence for a typical fluctuation coherence length of ~500 kpc, a non-uniform amplitude, and power spectra that are steeper than a Kolmogorov-like scaling.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper claims that the morphology of radio relics encodes the power spectrum of surrounding density fluctuations in the ICM, demonstrated through cosmologically motivated shock-tube simulations varying coherence length, amplitude, and power-law slope of upstream density fluctuations. Post-processing these with the Crest cosmic ray electron spectral solver produces mock relics showing independent morphological sensitivities: double strands to coherence length, downstream extent and patchiness to amplitude, filament number/spacing/curvature to slope. Application to the Toothbrush and Sausage relics suggests a typical coherence length of ~500 kpc, non-uniform amplitude, and steeper than Kolmogorov power spectra.

Significance. If the results hold, this represents a potentially important new observational window on density fluctuations at the edges of galaxy clusters, where traditional probes are insensitive. The strength lies in the forward-modeling approach with systematic variation of input parameters to map directly to observable features, providing falsifiable predictions for relic morphologies. This could complement other ICM studies if the density-specific sensitivities are confirmed to be robust.

major comments (1)
  1. [Simulation and post-processing methodology] The demonstration that morphological features map independently to density power spectrum parameters (e.g., double strands to coherence length, extended emission to amplitude) is based on varying only those parameters in shock-tube setups and applying Crest. However, Crest models CR electron spectra and synchrotron emission without self-consistent treatment of magnetic field geometry, shock amplification, or post-shock turbulence. These omitted processes can produce similar patchiness, filaments, and broadened emission, raising the possibility of degeneracies. This is central to the claim of a new window, as the abstract's application to Toothbrush/Sausage infers specific parameters without addressing such alternatives.
minor comments (2)
  1. The abstract states that 'double strand' features can be formed by curved shock fronts in projection, but it would be helpful to clarify in the main text how projection effects are quantified separately from the density fluctuation effects.
  2. [Application to observed relics] Details on the quantitative comparison or fitting procedure used to infer the ~500 kpc coherence length and other parameters from the Toothbrush and Sausage relics are not evident; including error estimates or sensitivity tests would strengthen the observational application.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive and insightful comments on our manuscript. We have carefully considered the major comment regarding the simulation and post-processing methodology and provide a point-by-point response below. Where appropriate, we have revised the manuscript to strengthen the discussion of limitations.

read point-by-point responses

  1. Referee: [Simulation and post-processing methodology] The demonstration that morphological features map independently to density power spectrum parameters (e.g., double strands to coherence length, extended emission to amplitude) is based on varying only those parameters in shock-tube setups and applying Crest. However, Crest models CR electron spectra and synchrotron emission without self-consistent treatment of magnetic field geometry, shock amplification, or post-shock turbulence. These omitted processes can produce similar patchiness, filaments, and broadened emission, raising the possibility of degeneracies. This is central to the claim of a new window, as the abstract's application to Toothbrush/Sausage infers specific parameters without addressing such alternatives.

    Authors: We agree that our shock-tube simulations with Crest post-processing represent an idealized setup that does not incorporate self-consistent magnetic field evolution, shock amplification, or post-shock turbulence. These processes could indeed contribute to morphological features such as patchiness, filaments, and broadened emission, potentially introducing degeneracies. Our methodology is deliberately designed to isolate the effects of upstream density fluctuations on shock morphology and Mach number distribution, which in turn drive variations in cosmic ray acceleration and synchrotron emission. The Crest solver assumes a uniform magnetic field strength for emission calculations, but the key morphological sensitivities (e.g., double strands from projected curved shocks, downstream extent from amplitude) arise primarily from the density-induced variations in shock properties. We have revised the manuscript by adding a dedicated paragraph in the Discussion section (new Section 5.3) that explicitly addresses these limitations, discusses possible degeneracies with magnetic field and turbulence effects, and notes that full MHD simulations would be required for a complete treatment. We have also updated the abstract and conclusions to temper the language around the applications to the Toothbrush and Sausage relics, emphasizing that the inferred parameters (coherence length ~500 kpc, non-uniform amplitude, steeper-than-Kolmogorov spectra) are illustrative and subject to the assumptions of the simplified model. This revision clarifies the scope of our claims without altering the core simulation results. revision: partial

Circularity Check

0 steps flagged

No circularity: forward modeling of varied density spectra produces independent morphological diagnostics

full rationale

The derivation proceeds by running shock-tube simulations in which coherence length, amplitude, and power-law index of upstream density fluctuations are independently varied as inputs, followed by post-processing with Crest to generate mock relics whose downstream features (strand scale, extent, filament spacing, patchiness) are then inspected for sensitivity. No equation defines an output quantity in terms of itself or a fitted parameter; no prediction is obtained by renaming a fitted input; no load-bearing step rests on a self-citation whose content is unverified or whose uniqueness theorem is imported from the same authors. The mapping from input power-spectrum parameters to observed morphology is therefore an empirical result of the simulation suite rather than a tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the fidelity of the simulation framework and post-processing; no free parameters are fitted to data in the claim itself, but the power spectrum parameters are inputs varied to demonstrate sensitivity.

axioms (2)
  • domain assumption Density fluctuations upstream of the shock can be fully characterized by a power spectrum with coherence length, amplitude, and power-law slope.
    This is the input varied systematically in the shock-tube simulations.
  • domain assumption The Crest cosmic ray electron spectral solver accurately produces mock radio relics from the simulated shocks without missing microphysical processes.
    Invoked in the post-processing step to generate the suite of mock relics used for the morphology analysis.

pith-pipeline@v0.9.0 · 5637 in / 1250 out tokens · 132547 ms · 2026-05-10T18:57:42.837414+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Revisiting radio synchrotron diagnostics in star-forming galaxies

    astro-ph.GA 2026-04 conditional novelty 7.0

    Advection-only galactic wind models fail to reproduce observed vertical radio profiles without unrealistic velocities, synchrotron spectra are biased toward young electrons in dense regions, and bremsstrahlung/Coulomb...

Reference graph

Works this paper leans on

6 extracted references · 1 canonical work pages · cited by 1 Pith paper

  1. [1]

    Angelinelli, M., Ettori, S., Vazza, F., & Jones, T. W. 2021, A&A, 653, A171 Axford, W. I., Leer, E., & Skadron, G. 1977, in International Cosmic Ray Con- ference, Vol. 11, International Cosmic Ray Conference, 132 Article number, page 14 of 20 Joseph Whittingham et al.: Zooming in on radio relics II. Bagchi, J., Sirothia, S. K., Werner, N., et al. 2011, Ap...

    work page arXiv 2021

  2. [2]

    At the edge of the main relic, there is a region of low-brightness (≈15μJy beam−1) emission, which bridges some of the gap

    and so is likely to have a common origin. At the edge of the main relic, there is a region of low-brightness (≈15μJy beam−1) emission, which bridges some of the gap. We highlight this with a dashed white ellipse. Moreover, this ‘bridge’ region has a lower average spectralindex(vanWeerenetal.2010).Thisistobeexpectedin a scenario where this feature forms th...

    2010

  3. [3]

    The dashed white ellipse indicates a potentially missing part of the relic (cf. Fig. 5). If produced by shock projection, the ‘knots’ along the relic imply a coherence length of∼500 kpc. This is consistent with our analysis of the Toothbrush relic in Fig. 4, and is difficult to explain through alternative mechanisms. 0.000 0.025 0.050 0.075 0.100 0.125 0....

    2022

  4. [4]

    This is because stronger shocks already have relatively flat spectral slope, and the consequent increaseinnormalisationislimited(seeSect.4.3ofW26).This, inturn,limitstheoverallabilityofdensityfluctuationstoincrease relic luminosity up to observed levels. E.2. Absolute intensity Observedradiorelicsreachfluxdensitiesbetweenoneandafew hundred millijanskys at...

    2011

  5. [5]

    A major additional boost could be realised through Fermi-I re-acceleration

    This will naturally be accounted for in future studies when we will attempt to replicate specific examples. A major additional boost could be realised through Fermi-I re-acceleration. As discussed in Sect. 2, we have assumed an initially thermal state for our CR electrons. This was chosen for twomainreasons:i)thetruedistributionofnon-thermalelectrons in c...

    2019

  6. [6]

    9, in particular)

    normalisation, and hence flux density, by a factor of103–104 in Mach 2 shocks and by a factor of 10 in Mach 3 shocks (see their Fig. 9, in particular). These values are also consistent with re- acceleration in so-called multi-shock scenarios, as modelled by Smolinskietal.(2023).IftheToothbrushrelichasameanMach numberofM=1.23,asquotedpreviously,itmaystandt...

    2023