arxiv: 2604.04068 · v2 · submitted 2026-04-05 · 🌌 astro-ph.GA

Recognition: no theorem link

The Distribution of Cosmic Ray Electrons in Star-Forming Galaxies

Anvar Shukurov , Charles Jose

Authors on Pith no claims yet

Pith reviewed 2026-05-13 17:11 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords cosmic ray electronsdiffusion equationGreen's functiongalaxy radio emissionstar-forming galaxiesenergy lossesaxisymmetric distributions

0 comments

The pith

Explicit algebraic expressions give the steady-state density of cosmic ray electrons in galaxies as a function of position and energy.

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

The paper derives explicit algebraic expressions for the number density of cosmic ray electrons at any position and energy in a galaxy. It solves the diffusion equation that includes energy losses using the Green's function method for sources distributed axisymmetrically in galactic radius and height above the plane. The expressions are first derived assuming Gaussian source distributions but are shown to work for general source distributions to within about 10 percent accuracy over broad ranges. This provides a practical way to model cosmic ray electrons without relying on the common assumption that their energy density matches that of the magnetic field.

Core claim

We derive explicit, algebraic expressions for the steady-state number density of cosmic ray electrons as a function of position and energy using Green's function of the diffusion equation with energy losses for axisymmetric distributions of the particle sources in the galactocentric radius r and distance to the mid-plane z. The solution is obtained for a Gaussian distribution of the particle sources in r and z but we show that it can be used for an arbitrary spatial distribution of the sources. The accuracy of our results is about 10% or better in wide ranges of r, z and particle energies. These solutions can be used in the interpretation of radio astronomical observations of galaxies, and a

What carries the argument

Green's function of the diffusion equation with energy losses for axisymmetric source distributions in galactocentric radius and height.

If this is right

The expressions enable direct calculation of cosmic ray electron densities for interpreting radio observations of galaxies.
They serve as a physically motivated alternative to assuming energy equipartition between cosmic rays and magnetic fields.
The formulas support efficient modeling of radio luminosities across large samples of galaxies.
The same Green's function approach applies to general source distributions while retaining useful accuracy.

Where Pith is reading between the lines

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

The expressions could improve magnetic field estimates in galaxies where radio data are combined with independent field measurements.
This framework might be extended to include time-dependent effects or additional loss processes for other cosmic ray species.
Predicted radio emission maps from these densities could be tested against high-resolution observations from next-generation telescopes.

Load-bearing premise

That the explicit solution derived for Gaussian source distributions can be applied to arbitrary spatial distributions with 10% or better accuracy across wide ranges of r, z, and energies.

What would settle it

A direct numerical solution of the diffusion equation for a non-Gaussian source distribution, such as a uniform disk, compared against the algebraic expression to check if deviations exceed 10 percent at some radii, heights or energies.

Figures

Figures reproduced from arXiv: 2604.04068 by Anvar Shukurov, Charles Jose.

**Figure 2.** Figure 2: The compensated energy spectra of relativistic electrons obtained from equation ( [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

**Figure 3.** Figure 3: As Fig [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗

**Figure 4.** Figure 4: The relative deviations of the approximate CRE spatial distributions [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗

read the original abstract

We derive explicit, algebraic expressions for the steady-state number density of cosmic ray electrons as a function of position and energy using Green's function of the diffusion equation with energy losses for an axisymmetric distributions of the particle sources in the galactocentric radius $r$ and distance to the mid-plane $z$. The solution is obtained for a Gaussian distribution of the particle sources in $r$ and $z$ but we show that it can be used for an arbitrary spatial distribution of the sources. The accuracy of our results is about 10\% or better in wide ranges of $r$, $z$ and particle energies. These solutions can be used in the interpretation of radio astronomical observations of galaxies, particularly in the studies of the radio luminosities for large galaxy samples, and represent a physically justifiable and efficient alternative to the assumption of the energy equipartition between cosmic rays and interstellar magnetic fields.

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.

Desk Editor's Note private letter to a colleague

Closed-form Green's function expressions for cosmic ray electron density in galaxies give a usable analytical alternative to equipartition, but the 10% accuracy claim for arbitrary source distributions needs explicit checks.

read the letter

The main point is that Shukurov and Jose have produced explicit algebraic expressions for the steady-state number density of cosmic ray electrons versus position and energy. They solve the diffusion-loss equation with Green's functions in an axisymmetric geometry, starting from Gaussian source profiles in galactocentric radius and height above the plane, then argue the same form applies to arbitrary axisymmetric sources at roughly 10% accuracy or better across wide ranges of r, z, and particle energy. The goal is to give radio observers a direct, efficient way to model synchrotron emission without defaulting to equipartition between cosmic rays and magnetic fields, especially for large galaxy samples. That is the practical advance. The derivation rests on the standard transport equation, which keeps the foundation familiar and the algebra tractable. If the full paper shows the Green's function construction and the superposition step clearly, readers can judge how much work is saved versus running full numerical codes. The soft spot is the accuracy statement for non-Gaussian sources. The abstract asserts 10% or better without showing the truncation or approximation procedure, without error bounds plotted against radius or energy, and without side-by-side comparisons to numerical solutions for exponential or other realistic disks. If those tests are absent or only sketched, the claim stays unverified and will need tightening in review. This paper is for people who model cosmic-ray transport or fit radio data to galaxy populations. A reader who needs quick analytic profiles for parameter studies or survey work will find the expressions worth testing. It is coherent on its own terms and formally grounded enough to merit referee time rather than a desk reject, even if the validation section requires expansion. I would send it to peer review.

Referee Report

1 major / 0 minor

Summary. The paper derives explicit algebraic expressions for the steady-state number density of cosmic ray electrons as a function of position and energy using the Green's function of the diffusion-loss equation for axisymmetric source distributions in galactocentric radius r and height z. The solution is obtained explicitly for Gaussian source profiles but is asserted to apply to arbitrary axisymmetric source distributions with an accuracy of ~10% or better over wide ranges of r, z, and particle energy; the expressions are positioned as a physically motivated alternative to energy equipartition for interpreting radio observations of galaxies.

Significance. If the accuracy claim for arbitrary distributions holds after validation, the work supplies a useful closed-form analytical tool that could streamline modeling of cosmic-ray electron distributions for large galaxy samples and radio-luminosity studies. The Green's-function approach for the diffusion-loss equation is a clear technical strength, offering reproducibility and efficiency where numerical solutions are otherwise required.

major comments (1)

[Abstract] Abstract: The central claim that the Gaussian-derived Green's function solution applies to arbitrary axisymmetric source distributions with 10% or better accuracy lacks any derivation of the approximation procedure, quantitative error bounds, or direct comparisons to numerical solutions for non-Gaussian profiles such as exponential disks. This assertion is load-bearing for the paper's claimed utility and must be supported by explicit validation in the manuscript.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thoughtful review and for highlighting the need for stronger validation of our approximation. We agree that explicit quantitative support is essential for the claimed utility and will incorporate the requested comparisons and error analysis in the revised manuscript.

read point-by-point responses

Referee: [Abstract] Abstract: The central claim that the Gaussian-derived Green's function solution applies to arbitrary axisymmetric source distributions with 10% or better accuracy lacks any derivation of the approximation procedure, quantitative error bounds, or direct comparisons to numerical solutions for non-Gaussian profiles such as exponential disks. This assertion is load-bearing for the paper's claimed utility and must be supported by explicit validation in the manuscript.

Authors: We acknowledge the referee's point that the manuscript would be strengthened by explicit validation. The original text derives the closed-form solution for Gaussian sources and states that the same Green's function form can be applied to arbitrary axisymmetric profiles via direct integration, with the ~10% accuracy asserted on the basis of the smoothness of the diffusion-loss kernel. However, we agree that this requires quantitative demonstration. In the revision we will add a dedicated subsection (new Section 4.3) that (i) outlines the approximation procedure (convolving the Green's function with the source profile and comparing to the exact integral), (ii) presents error maps and maximum relative errors for exponential-disk and other non-Gaussian profiles across r, z, and E ranges relevant to radio observations, and (iii) shows that the error remains below 10% for the majority of the parameter space. These additions will be supported by new figures and a brief appendix with the numerical integration method used for validation. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper starts from the standard diffusion equation with energy losses and constructs explicit algebraic solutions via the Green's function method for Gaussian source distributions in (r,z), then asserts (without shown derivation of error bounds) that the same form applies to arbitrary axisymmetric sources at ~10% accuracy. No step reduces a claimed prediction to a fitted parameter by construction, no self-citation is invoked as a load-bearing uniqueness theorem, and no ansatz is smuggled via prior work. The central algebraic result is obtained directly from the diffusion-loss operator and remains independent of the target observables or fitted quantities.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard domain assumptions about cosmic ray transport; no free parameters or invented entities are mentioned in the abstract.

axioms (2)

domain assumption The diffusion equation with energy losses governs steady-state cosmic ray electron transport in galaxies
Invoked to obtain the Green's function solution for number density.
domain assumption Steady-state condition applies to the electron number density
Used to derive the equilibrium distribution from the time-dependent diffusion equation.

pith-pipeline@v0.9.0 · 5449 in / 1378 out tokens · 61886 ms · 2026-05-13T17:11:32.663188+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages · 2 internal anchors

[1]

and Blasi, P

Amato, E. and Blasi, P. (2018). Cosmic ray transport in the Galaxy: A review. Adv. Space Res. 62, 2731–2749. doi:10.1016/j.asr.2017.04.019 Atoyan, A. M., Aharonian, F. A., and Völk, H. J. (1995). Electrons and positrons in the galactic cosmic rays. Phys. Rev. D 52, 3265–3275. doi:10.1103/PhysRevD.52.3265 Frontiers 17 Shukurov and Jose Cosmic Ray Electrons...

work page doi:10.1016/j.asr.2017.04.019 2018
[2]

and Amato, E

doi:10.1007/ s00159-013-0070-7 Blasi, P. and Amato, E. (2012). Diffusive propagation of cosmic rays from supernova remnants in the Galaxy. I: spectrum and chemical composition. J. Cosmology Astropart. Phys. 2012,

work page 2012
[3]

A., Kovalenko, V

doi:10.1088/1475-7516/2012/01/010 Dogiel, V. A., Kovalenko, V. M., and Prishchep, V. L. (1980). Cosmic-ray electrons in the diffusion- convection model of particle propagation. Sov. Astron. Lett. 6, 366–368 Draine, B. T. (2011). Physics of the Interstellar and Intergalactic Medium (Princeton: Princeton Univ. Press) Evoli, C., Gaggero, D., Grasso, D., and ...

work page doi:10.1088/1475-7516/2012/01/010 2012
[4]

doi:10.1088/1475-7516/2008/10/018 Ginzburg, V. L. and Syrovatskii, S. I. (1964). The Origin of Cosmic Rays (Oxford: Pergamon) Hanasz, M., Strong, A. W., and Girichidis, P. (2021). Simulations of cosmic ray propagation. Living Reviews in Computational Astrophysics 7,

work page doi:10.1088/1475-7516/2008/10/018 2008
[5]

P., Lagos, C

doi:10.1007/s41115-021-00011-1 Hansen, S. P., Lagos, C. D. P., Bonato, M., Cook, R. H. W., Davies, L. J. M., Delvecchio, I., et al. (2024). Modelling the galaxy radio continuum from star formation and active galactic nuclei in the SHARK semi-analytic model. MNRAS 531, 1971–1987. doi:10.1093/mnras/stae1235 Heesen, V. (2021). The radio continuum perspective...

work page doi:10.1007/s41115-021-00011-1 2024
[6]

Cosmic Rays on Galaxy Scales: Progress and Pitfalls for CR-MHD Dynamical Models

doi:10.1007/s10509-021-04026-1 Heesen, V., Krause, M., Beck, R., Adebahr, B., Bomans, D. J., Carretti, E., et al. (2018). Radio haloes in nearby galaxies modelled with 1D cosmic ray transport using SPINNAKER. MNRAS 476, 158–183. doi:10.1093/mnras/sty105 Hopkins, P. F. (2025). Cosmic rays on galaxy scales: progress and pitfalls for CR–MHD dynamical models....

work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/s10509-021-04026-1 2018
[7]

doi:10.3390/galaxies12030022 Jose, C., Chamandy, L., Shukurov, A., Subramanian, K., Rodrigues, L. F. S., and Baugh, C. M. (2024). Understanding the radio luminosity function of star-forming galaxies and its cosmological evolution. MNRAS 532, 1504–1521. doi:10.1093/mnras/stae1426 Kroupa, P. (2001). On the variation of the initial mass function. MNRAS 322, ...

work page doi:10.3390/galaxies12030022 2024
[8]

(Cambridge: Cambridge Univ

Stars, the Galaxy and the Interstellar Medium. (Cambridge: Cambridge Univ. Press) Lorimer, D. R. (2004). The Galactic population and birth rate of radio pulsars. In Young Neutron Stars and Their Environments, eds. F. Camilo and B. M. Gaensler. IAU Symposium 218, p. 105–

work page 2004
[9]

D., Fletcher, A., Beck, R., Mitra, D., and Scaife, A

doi:10.48550/arXiv.astro-ph/0308501 Mulcahy, D. D., Fletcher, A., Beck, R., Mitra, D., and Scaife, A. M. M. (2016). Modelling the cosmic ray electron propagation in M51. Astron. Astrophys. 592, A123. doi:10.1051/0004-6361/ 201628446 Popescu, C. C., Yang, R., Tuffs, R. J., Natale, G., Rushton, M., and Aharonian, F. (2017). A radiation transfer model for th...

work page internal anchor Pith review doi:10.48550/arxiv.astro-ph/0308501 2016
[10]

doi:10.1093/mnras/stae138 Schleicher, D

MNRAS 528, 1596–1603. doi:10.1093/mnras/stae138 Schleicher, D. R. G. and Beck, R. (2013). A new interpretation of the far-infrared–radio correlation and the expected breakdown at high redshift. Astron. Astrophys. 556, A142. doi:10.1051/ 0004-6361/201321707 Schlickeiser, R. (2002). Cosmic Ray Astrophysics (Berlin: Springer) Schober, J., Sargent, M. T., Kle...

work page doi:10.1093/mnras/stae138 2013
[11]

doi:10.3390/galaxies7020045 Stein, M., Heesen, V., Dettmar, R.-J., Stein, Y., Brüggen, M., Beck, R., et al. (2023). CHANG-ES. XXVI. Insights into cosmic-ray transport from radio halos in edge-on galaxies. Astron. Astrophys. 670, A158. doi:10.1051/0004-6361/202243906 Strong, A. W., Moskalenko, I. V., and Ptuskin, V. S. (2007). Cosmic-ray propagation and in...

work page doi:10.3390/galaxies7020045 2023
[12]

doi:10.3847/1538-4357/ad7385 Frontiers 20

work page doi:10.3847/1538-4357/ad7385