Hydrodynamic model of nonthermal emission from the Fermi bubbles
Pith reviewed 2026-05-25 03:27 UTC · model grok-4.3
The pith
Rayleigh-Taylor instabilities in the Fermi Bubbles shell accelerate cosmic ray electrons to TeV energies via stochastic processes without strong shocks.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
At late evolutionary stages the Fermi Bubbles shell develops Rayleigh-Taylor instabilities whose associated MHD fluctuations stochastically accelerate cosmic-ray electrons to the TeV energies needed to explain the observed gamma-ray and microwave emission; this process operates without requiring strong shock fronts.
What carries the argument
Spectrum of Rayleigh-Taylor instabilities that drives MHD fluctuations for stochastic cosmic-ray acceleration, together with the derived kinetic equations and acceleration timescale.
If this is right
- The model removes the need for high-Mach-number shocks inside the Fermi Bubbles.
- Electron acceleration occurs in situ within the shell rather than at a distant shock front.
- Derived spectra of MHD fluctuations and kinetic equations can be compared with future radio or gamma-ray observations of the bubble edges.
- The required energy release from the Galactic Center can be reassessed under the lower shock-strength constraint.
Where Pith is reading between the lines
- Similar late-stage instability acceleration could operate in other expanding galactic structures such as supernova remnants or superbubbles.
- If confirmed, the model predicts a characteristic spatial gradient in particle spectra across the bubble shell that could be searched for with high-resolution gamma-ray telescopes.
- The approach may reduce the total energy budget attributed to past activity at the Galactic Center.
Load-bearing premise
Rayleigh-Taylor instabilities at late stages of the bubble shell generate enough MHD fluctuations to accelerate electrons to TeV energies on the available timescale.
What would settle it
A direct measurement of the turbulence spectrum or acceleration timescale in the bubble shell that falls short of the TeV energies required to match the observed gamma-ray and microwave fluxes.
Figures
read the original abstract
We suggest a model of Fermi Bubbles (FBs) in the Galactic halo of the altitude about 7-8 kpc, which is seen in non-thermal microwave and gamma-ray ranges. It was assumed that this emission is generated by relativistic electrons of cosmic rays whose origin is still under debate. It has been assumed that the FB shell is generated in the halo by the release of energy, generated by the routine capture of stars at the central black hole of the Galactic Centre (GC). In this case cosmic ray electrons (CR) in the shells of the FBs of sufficiently high energies are generated by the standard shock acceleration. However, one of the problems of this model is that the Mach number of the FB shock is not high enough to generate the observed non-thermal radiation from the halo. We propose an alternative model of stochastic CR acceleration by Rayleigh-Taylor (RT) instabilities in the shell of the FB at the late stages of the evolution of the shell in the halo. Unlike the shock model of CR acceleration, the RT model of in-situ acceleration in the FBs does not require strong shock fronts. In our model, we derived the spectrum of RT instabilities and estimate the spectra of kinetic equations for MHD-fluctuations needed for acceleration of CRs. We assessed the time of CR electron acceleration up to TeV energies that needed to interpret the observed data of gamma-ray and microwave emission from the envelope of FBs.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a hydrodynamic model for nonthermal emission from the Fermi Bubbles at ~7-8 kpc altitude, attributing the observed microwave and gamma-ray signals to relativistic electrons accelerated stochastically by Rayleigh-Taylor instabilities in the expanding shell at late evolutionary stages. This is presented as an alternative to shock acceleration, which the authors note suffers from insufficient Mach numbers; they derive the RT instability spectrum, estimate the resulting MHD fluctuation spectra for the kinetic equations, and assess the time required for electrons to reach TeV energies.
Significance. If the time assessment and spectral mapping hold, the work supplies a concrete in-situ acceleration channel that avoids the need for strong shocks, addressing a documented limitation of prior FB models. The explicit derivation of the RT wavenumber spectrum and its connection to Fokker-Planck coefficients for CR electrons constitutes a falsifiable framework that could be tested against multi-wavelength data.
major comments (2)
- [Time-assessment paragraph (abstract and corresponding derivation section)] The central claim that RT-driven acceleration reaches TeV energies within the FB shell lifetime at 7-8 kpc rests on the quantitative mapping from the derived RT instability spectrum to the stochastic acceleration rate; the abstract states an assessment was performed, but the growth-rate saturation level and resulting diffusion coefficient in the kinetic equation are not shown to produce t_acc shorter than the dispersal timescale when the shell has expanded into the stratified halo.
- [Section deriving RT spectrum and MHD fluctuation spectra] The spectrum of MHD fluctuations is obtained from the RT instability spectrum, yet the transition from hydrodynamic growth rates to the amplitude of the wave spectrum that enters the Fokker-Planck coefficients for electron acceleration is stated without an explicit efficiency factor or saturation amplitude; this step is load-bearing for the claim that the RT model operates without strong shocks.
minor comments (2)
- [Abstract] The abstract contains repetitive phrasing ('It was assumed' appears twice in the first paragraph) and minor grammatical awkwardness that reduces readability.
- [Derivation sections] Notation for the wavenumber spectrum of RT instabilities and its conversion to the MHD fluctuation power spectrum is introduced without a clear equation label or reference to prior work on RT-MHD coupling, which would aid readers.
Simulated Author's Rebuttal
We thank the referee for the thoughtful review and for highlighting the importance of making the saturation assumptions and efficiency factors explicit. We address each major comment below and agree that revisions are needed to strengthen the quantitative claims.
read point-by-point responses
-
Referee: [Time-assessment paragraph (abstract and corresponding derivation section)] The central claim that RT-driven acceleration reaches TeV energies within the FB shell lifetime at 7-8 kpc rests on the quantitative mapping from the derived RT instability spectrum to the stochastic acceleration rate; the abstract states an assessment was performed, but the growth-rate saturation level and resulting diffusion coefficient in the kinetic equation are not shown to produce t_acc shorter than the dispersal timescale when the shell has expanded into the stratified halo.
Authors: The manuscript derives the RT wavenumber spectrum and maps it to a stochastic acceleration timescale via the resulting MHD fluctuations, with the abstract assessment based on that mapping. However, the referee is correct that the saturation amplitude of the growth rates and the explicit diffusion coefficient D(p) entering the Fokker-Planck equation are not displayed numerically. We will revise the derivation section to include the saturation level (taken from nonlinear RT theory) and the computed t_acc value, demonstrating it remains below the dispersal timescale at 7-8 kpc. revision: yes
-
Referee: [Section deriving RT spectrum and MHD fluctuation spectra] The spectrum of MHD fluctuations is obtained from the RT instability spectrum, yet the transition from hydrodynamic growth rates to the amplitude of the wave spectrum that enters the Fokker-Planck coefficients for electron acceleration is stated without an explicit efficiency factor or saturation amplitude; this step is load-bearing for the claim that the RT model operates without strong shocks.
Authors: The transition assumes that a fraction of the RT mode energy converts to MHD wave energy. The referee correctly notes the absence of an explicit efficiency factor or saturation amplitude. In revision we will introduce a dimensionless efficiency parameter η (with a plausible range justified by RT simulations in stratified media) and show how it scales the wave amplitude that enters the Fokker-Planck coefficients, thereby making the mapping fully transparent. revision: yes
Circularity Check
No circularity: RT spectrum and acceleration timescale derived independently from first principles
full rationale
The paper claims to derive the RT instability spectrum from the hydrodynamic evolution of the FB shell and then estimate the resulting MHD fluctuation spectrum for the kinetic equations, followed by an explicit assessment of the electron acceleration timescale to TeV energies. No equations are shown that define the output spectrum in terms of a fitted input parameter, rename a known result, or reduce the central prediction to a self-citation chain. The abstract presents the RT model as an alternative whose viability is tested against the observed emission, without load-bearing self-citations or ansatzes smuggled from prior work by the same authors. This satisfies the default expectation that the derivation chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Rayleigh-Taylor instabilities develop in the FB shell at late evolutionary stages and generate the MHD fluctuations required for stochastic acceleration
Reference graph
Works this paper leans on
-
[1]
Ackermann, M., Albert, A., Atwood, W. B. et al. 2014, ApJ, 793, 64
work page 2014
-
[2]
K., Theory of Homogeneous Turbulence , Cambridge University Press, New York, 1970
Batchelor, G. K., Theory of Homogeneous Turbulence , Cambridge University Press, New York, 1970
work page 1970
- [3]
-
[4]
Berezinskii, V. S., Bulanov, S. V., Dogiel, V. A., Ginzburg, V. L., & Ptuskin, V. S. 1990, Astrophysics of Cosmic Rays , ed. V.L.Ginzburg, (Norht-Holland, Amsterdam)
work page 1990
- [5]
- [6]
- [7]
-
[8]
Cheng, K.-S., Chernyshov, D. O. & Dogiel, V. A. 2007, A&A, 473, 351
work page 2007
-
[9]
Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., Ko, C. M., & Ip, W.-H. 2011, ApJL, 731, L17
work page 2011
-
[10]
Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., & Ko, C. M. 2014, ApJ, 790, 23
work page 2014
-
[11]
Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., & Ko, C. M. 2015a, ApJ, 799, 112
-
[12]
Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., & Ko, C. M. 2015b, ApJ, 804, 135
- [13]
-
[14]
Chernyshov, D. O., Dogiel, V. A., Ivlev, A. V. et al. 2022, ApJ, 937, 107
work page 2022
- [15]
-
[16]
Cho, J., Lazarian, A., & Vishniac, E. T. 2002, ApJ, 564, 291
work page 2002
-
[17]
Cordes, J. M., Weisberg, J. M., Frail, D. A. et al. 1991, Nature, 354, 121
work page 1991
- [18]
-
[19]
Crocker, R. M., Johnes, & Aharonian, F. 2011, MNRAS, 413, 763
work page 2011
-
[20]
Dogiel, V. A. 2000, A&A, 357, 66
work page 2000
-
[21]
Dogiel, V. A., Ivlev, A. V., Chernyshov, D. O., & Ko, C. M. 2020, ApJ, 903, 135
work page 2020
- [22]
- [23]
- [24]
-
[25]
Fujita, Y., Takizawa, M., & Sarazin C. L. 2003, ApJ, 584, 190
work page 2003
-
[26]
Gaensler, B. M., Madsen, G. J., Chatterjee, S., & Mao, S. A. 2008, PASA, 25, 184
work page 2008
-
[27]
Ginzburg, V. L., & Syrovatskii, S. I. 1964, The Origin of Cosmic Rays , New York: Macmillan (New York: Macmillan)
work page 1964
-
[28]
Ginzburg, V. L. Theoretical physics and astrophysics , Oxford, Pergamon Press, (International Series in Natural Philosophy. Volume 99), 1979
work page 1979
- [29]
- [30]
-
[31]
Goldwurm, A., Clavel, M., Gabici, S., & Terrier, R. 2026, arXiv:2602.14148
-
[32]
Grotova, I., Rau, A., Baldini, P. et al. 2025, A&A, 697, A159
work page 2025
-
[33]
Gurevich, A. V. 1960, Sov. Phys. JETP, 38, 1150
work page 1960
- [34]
- [35]
- [36]
- [37]
- [38]
- [39]
-
[40]
M., Breitschwerdt, D., Chernyshov, D
Ko, C. M., Breitschwerdt, D., Chernyshov, D. O., Cheng, H., Dai, L., & Dogiel, V. A. 2020, ApJ , 904, 46
work page 2020
-
[41]
N., 1941, C R Acad Sci URSS, 30, 201
Kolmogoroff, A. N., 1941, C R Acad Sci URSS, 30, 201
work page 1941
-
[42]
Kompaneets, A. S. 1960, Akademiia Nauk SSSR, Doklady (DoSSR, in Russian), 130, 5
work page 1960
- [43]
- [44]
-
[45]
Ligthill, M. J. 1952, Proc, R. Soc. Lond. A, 211, 564
work page 1952
- [46]
- [47]
-
[48]
Miller, J.A., & Roberts, D.A. ApJ. 1995, 452, 912
work page 1995
- [49]
- [50]
- [51]
-
[52]
Nakashima, S., Inoue, Y., Yamasaki, N. et al. 2018, ApJ, 862, 34
work page 2018
-
[53]
Nakashima, S., Koyama, K., Wang, Q. D. et al. 2019, ApJ, 875, 32
work page 2019
- [54]
- [55]
-
[56]
Obukhov, A. M. 1941, Dokl. Akad. Nauk SSSR, 32, 22
work page 1941
-
[57]
Parker, E. N. 1964, ApJ, 140, 1170
work page 1964
- [58]
-
[59]
Peretti, E., Morlino, G., Blasi, P. & Cristofari, P. 2022, MNRAS, 511, 1336
work page 2022
-
[60]
Planck Collaboration IX, 2013, A&A, 554, A139
work page 2013
- [61]
-
[62]
Sarkar, K. C. (2024), A&A Rev., 32, 1
work page 2024
-
[63]
, Schulreich, M.M., Schurer, D
Scheffler, T. , Schulreich, M.M., Schurer, D. P. P. R., & Breitschwerdt, D. 2025, A&A, 695A, 34
work page 2025
- [64]
- [65]
- [66]
-
[67]
Snowden, S. L., Egger, R., Freyberg, M. J. et al. 1997, ApJ, 485, 125
work page 1997
- [68]
-
[69]
Sreenivasan, K. R. 1995, Phys. Fluids, 7, ll
work page 1995
- [70]
- [71]
- [72]
-
[73]
Vladimirov, A. E., Digel, S. W., J \'o hannesson, G., et al. 2011, CoPhC, 182, 1156
work page 2011
-
[74]
Tseng, P.-H., Yang, H. -Y. K., Chen, Ch.-Y. et al 2024, ApJ, 970, 146
work page 2024
- [75]
-
[76]
Yang, C. & Li, S. 2013, ApJ, 773, 138
work page 2013
-
[77]
Yang, H. Y. K., Ruszkowski, M., Ricker, P. M., Zweibel, E., & Lee, D. 2012, ApJ, 761, 185
work page 2012
-
[78]
Yang, H.-Y., Ruszkowski, M., & Zweibel, E. 2018, Galaxies, 6, 29
work page 2018
-
[79]
K., Ruszkowski, M., & Zweibel, E
Yang, H.-Y. K., Ruszkowski, M., & Zweibel, E. G. 2022, Nature Astronomy, 6, 584
work page 2022
-
[80]
Zhou, Ye, & Matthaeus W. H. 1990, J. Geophys. Res., 95, 14
work page 1990
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.