pith. sign in

arxiv: 2605.23741 · v1 · pith:T3KNWLH2new · submitted 2026-05-22 · 🌌 astro-ph.HE · astro-ph.GA

Hydrodynamic model of nonthermal emission from the Fermi bubbles

Pith reviewed 2026-05-25 03:27 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.GA
keywords Fermi bubblesRayleigh-Taylor instabilitiescosmic ray accelerationnonthermal emissiongalactic halostochastic accelerationgamma-ray emissionmicrowave emission
0
0 comments X

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.

The paper proposes that nonthermal microwave and gamma-ray emission from the Fermi Bubbles arises from relativistic electrons accelerated in situ within the bubble shell. It identifies a problem with the standard shock-acceleration picture: the Mach number of the forward shock is too low to produce the required particle energies. Instead, the authors model late-stage Rayleigh-Taylor instabilities that develop as the shell expands into the halo, generating MHD fluctuations capable of stochastic acceleration. They derive the instability spectrum, solve the kinetic equations for the resulting turbulence, and calculate the time required for electrons to reach TeV energies, showing consistency with observed emission.

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

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

  • 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

Figures reproduced from arXiv: 2605.23741 by A.M. Kiselev, C.M. Ko, D.O. Chernyshov, T.S. Fatekhov, V.A. Dogiel.

Figure 1
Figure 1. Figure 1: Dependence of the transformed time y˜ on the dimensionless time t˜. and define the time t˜ as a function of the variable y˜, t˜(˜y) =   3 2 Z y˜ 0 q V˜ (˜y)dy˜   2/3 . (12) The function y˜(t˜) is presented in [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Left, Position of the top of the bubble as function of the dimensionless time t˜. Right, Acceleration of the FB envelope at the top the shell. The acceleration at the top of the bubble, z¨u(y), may produce the RT instability if the envelope acceleration is high enough. Inside the envelope the RT instabilities are excited between the dense shell ρsh and the hot interior ρin (see Baumgartner & Breitschwerdt,… view at source ↗
Figure 3
Figure 3. Figure 3: Left.The amplitude of the Rayleigh-Taylor instabilities η(t) and the thickness of the bubble d(t). Right. The ratio of the luminosity of RT turbulence to the total luminosity at the GC derived from LRT(t˜)/L. where ρ(zu) is given by Eq. (2), vturb = v0 and lturb = 2π/k0 (see below). We estimate the mass of turbulent gas as the mass of the bubble above the height with maximum radius, mturb = 2π Z zu z0 ρ0e … view at source ↗
Figure 4
Figure 4. Figure 4: The thick dashed line shows the normalized spectru [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: η˜ 2/3/η˜˙2 as function of dimensionless time t˜. the acceleration time can be estimated as τacc = p 2 4D(p) ∼ 2 2/3 c π η 2/3 η˙(t) 2  η˙(t) vA(t) 1/3  pc eB 1/3 . (48) As η = Hη˜ and η˙ = Hη/t ˜˙ 0, τacc = 2 2/3 c π η˜ 2/3 η˜˙(t) 2  η˙(t) vA(t) 1/3 r 1/3 L t 2 0 H4/3 , (49) where rL is the Larmor radius of electron. The acceleration time depends on the combination η 2/3/η˙ 2 at the position of the … view at source ↗
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.

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

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [Abstract] The abstract contains repetitive phrasing ('It was assumed' appears twice in the first paragraph) and minor grammatical awkwardness that reduces readability.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

Review performed on abstract only; full text unavailable so ledger entries are limited to statements explicit in the abstract.

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
    Central modeling choice stated in the abstract as the basis for the alternative mechanism.

pith-pipeline@v0.9.0 · 5813 in / 1182 out tokens · 53147 ms · 2026-05-25T03:27:43.523632+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

80 extracted references · 80 canonical work pages

  1. [1]

    Ackermann, M., Albert, A., Atwood, W. B. et al. 2014, ApJ, 793, 64

  2. [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

  3. [3]

    2013, A&A, 557, 140

    Baumgartner, V., & Breitschwerdt, D. 2013, A&A, 557, 140

  4. [4]

    S., Bulanov, S

    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)

  5. [5]

    2025, arXiv: 20628

    Bian, X., Betti, R., Zhao, D, & Aluie, H. 2025, arXiv: 20628

  6. [6]

    S., & Silich, S

    Bisnovatyi-Kogan, G. S., & Silich, S. A. 1995, RvMP, 67, 661

  7. [7]

    2018, JCAP, 07 063

    Biswas, S., & Gupta, N. 2018, JCAP, 07 063

  8. [8]

    Cheng, K.-S., Chernyshov, D. O. & Dogiel, V. A. 2007, A&A, 473, 351

  9. [9]

    O., Dogiel, V

    Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., Ko, C. M., & Ip, W.-H. 2011, ApJL, 731, L17

  10. [10]

    O., Dogiel, V

    Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., & Ko, C. M. 2014, ApJ, 790, 23

  11. [11]

    O., Dogiel, V

    Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., & Ko, C. M. 2015a, ApJ, 799, 112

  12. [12]

    O., Dogiel, V

    Cheng, K.-S., Chernyshov, D. O., Dogiel, V. A., & Ko, C. M. 2015b, ApJ, 804, 135

  13. [13]

    O., Dogiel, V

    Chernyshov, D. O., Dogiel, V. A. & Ko, C.-M. 2012, ApJ, 759, 113

  14. [14]

    O., Dogiel, V

    Chernyshov, D. O., Dogiel, V. A., Ivlev, A. V. et al. 2022, ApJ, 937, 107

  15. [15]

    2002, PRL, 88, 245001

    Cho, J., & Lazarian, A. 2002, PRL, 88, 245001

  16. [16]

    Cho, J., Lazarian, A., & Vishniac, E. T. 2002, ApJ, 564, 291

  17. [17]

    M., Weisberg, J

    Cordes, J. M., Weisberg, J. M., Frail, D. A. et al. 1991, Nature, 354, 121

  18. [18]

    M., & Aharonian, F

    Crocker, R. M., & Aharonian, F. 2011, PhRvL, 106, 101102

  19. [19]

    M., Johnes, & Aharonian, F

    Crocker, R. M., Johnes, & Aharonian, F. 2011, MNRAS, 413, 763

  20. [20]

    Dogiel, V. A. 2000, A&A, 357, 66

  21. [21]

    A., Ivlev, A

    Dogiel, V. A., Ivlev, A. V., Chernyshov, D. O., & Ko, C. M. 2020, ApJ, 903, 135

  22. [22]

    A., & Ko, C

    Dogiel, V. A., & Ko, C. M., 2024, Universe, 10, 424

  23. [23]

    A., & Breitschwerdt, D

    Dorfi, E. A., & Breitschwerdt, D. 2012, A&A, 540, A77

  24. [24]

    A., & Henriksen, R

    Eilek, J. A., & Henriksen, R. N. 1984, ApJ, 277, 820

  25. [25]

    Fujita, Y., Takizawa, M., & Sarazin C. L. 2003, ApJ, 584, 190

  26. [26]

    M., Madsen, G

    Gaensler, B. M., Madsen, G. J., Chatterjee, S., & Mao, S. A. 2008, PASA, 25, 184

  27. [27]

    L., & Syrovatskii, S

    Ginzburg, V. L., & Syrovatskii, S. I. 1964, The Origin of Cosmic Rays , New York: Macmillan (New York: Macmillan)

  28. [28]

    Ginzburg, V. L. Theoretical physics and astrophysics , Oxford, Pergamon Press, (International Series in Natural Philosophy. Volume 99), 1979

  29. [29]

    1995, ApJ, 438, 763

    Goldreich, P., & Sridhar, S. 1995, ApJ, 438, 763

  30. [30]

    1997, ApJ, 485, 680

    Goldreich, P., & Sridhar, S. 1997, ApJ, 485, 680

  31. [31]

    2026, arXiv:2602.14148

    Goldwurm, A., Clavel, M., Gabici, S., & Terrier, R. 2026, arXiv:2602.14148

  32. [32]

    Grotova, I., Rau, A., Baldini, P. et al. 2025, A&A, 697, A159

  33. [33]

    Gurevich, A. V. 1960, Sov. Phys. JETP, 38, 1150

  34. [34]

    1964, Sov

    Iroshnikov P.S. 1964, Sov. Astron., 7, 566

  35. [35]

    N., Bridle, A

    Hendriksen, R. N., Bridle, A. H., & Chan, K. L. 1982, ApJ, 257, 63

  36. [36]

    2013, ApJ, 779, 57

    Kataoka, J., Tahara, M., Totani, T., et al. 2013, ApJ, 779, 57

  37. [37]

    2015, ApJ, 807, 77

    Kataoka, J., Tahara, M., Totani, T., et al. 2015, ApJ, 807, 77

  38. [38]

    1968, PASJ, 20, 59

    Kato, S. 1968, PASJ, 20, 59

  39. [39]

    2026, arXiv:2601.14369

    Keshet, U., & Ghosh, A. 2026, arXiv:2601.14369

  40. [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

  41. [41]

    N., 1941, C R Acad Sci URSS, 30, 201

    Kolmogoroff, A. N., 1941, C R Acad Sci URSS, 30, 201

  42. [42]

    Kompaneets, A. S. 1960, Akademiia Nauk SSSR, Doklady (DoSSR, in Russian), 130, 5

  43. [43]

    H., & Nagarajan, S

    Kraichnan, R. H., & Nagarajan, S. 1967, Phys. Fluids, 10, 859

  44. [44]

    1955, ApJ, 121, 461

    Kulsrud, R. 1955, ApJ, 121, 461

  45. [45]

    Ligthill, M. J. 1952, Proc, R. Soc. Lond. A, 211, 564

  46. [46]

    2011, PRL, 107, 091101

    Mertsch, P., & Sarkar, S. 2011, PRL, 107, 091101

  47. [47]

    2019, A&A, 622, A203

    Mertsch, P., & Petrosian, V. 2019, A&A, 622, A203

  48. [48]

    Miller, J.A., & Roberts, D.A. ApJ. 1995, 452, 912

  49. [49]

    J., & Bregman, J

    Miller, M. J., & Bregman, J. N. 2016, ApJ, 829, 9

  50. [50]

    1998, ApJ, 493, 694

    Moskalenko, I.V., & Strong A.W. 1998, ApJ, 493, 694

  51. [51]

    2015, ApJ, 811, 37

    Mou, G., Yuan, F., Gan, Zh., Sun, M. 2015, ApJ, 811, 37

  52. [52]

    Nakashima, S., Inoue, Y., Yamasaki, N. et al. 2018, ApJ, 862, 34

  53. [53]

    Nakashima, S., Koyama, K., Wang, Q. D. et al. 2019, ApJ, 875, 32

  54. [54]

    2018, MNRAS, 478, L127

    Nayakshin, S., & Zubovas, K. 2018, MNRAS, 478, L127

  55. [55]

    E., Cordes, J

    Nordgren, T. E., Cordes, J. M., & Terzian, Y. 1992, AJ, 104, 1465

  56. [56]

    Obukhov, A. M. 1941, Dokl. Akad. Nauk SSSR, 32, 22

  57. [57]

    Parker, E. N. 1964, ApJ, 140, 1170

  58. [58]

    & East, W.E

    Petrosian, V. & East, W.E. 2008, ApJ, 682, 175

  59. [59]

    & Cristofari, P

    Peretti, E., Morlino, G., Blasi, P. & Cristofari, P. 2022, MNRAS, 511, 1336

  60. [60]

    Planck Collaboration IX, 2013, A&A, 554, A139

  61. [61]

    A., Becker, W

    Predehl, P., Sunyaev, R. A., Becker, W. et al. 2020, Nature, 588, 227

  62. [62]

    Sarkar, K. C. (2024), A&A Rev., 32, 1

  63. [63]

    , Schulreich, M.M., Schurer, D

    Scheffler, T. , Schulreich, M.M., Schurer, D. P. P. R., & Breitschwerdt, D. 2025, A&A, 695A, 34

  64. [64]

    2022, MNRAS, 509, 716

    Schulreich, M., & Breitschwerdt, D. 2022, MNRAS, 509, 716

  65. [65]

    2025, PhRvL, 134, 5201

    Schroer, B., Caprioli, D., & Blasi, P. 2025, PhRvL, 134, 5201

  66. [66]

    2024, ApJ, 973, 78

    Shimoda, J., & Asano, K. 2024, ApJ, 973, 78

  67. [67]

    L., Egger, R., Freyberg, M

    Snowden, S. L., Egger, R., Freyberg, M. J. et al. 1997, ApJ, 485, 125

  68. [68]

    1977, A&A, 60, 327

    Sofue, Y. 1977, A&A, 60, 327

  69. [69]

    Sreenivasan, K. R. 1995, Phys. Fluids, 7, ll

  70. [70]

    1994, ApJ, 432, 612

    Sridhar, S., & Goldreich, P. 1994, ApJ, 432, 612

  71. [71]

    C., & Metzger, B

    Stone, N. C., & Metzger, B. D. 2016, MNRAS, 455, 859

  72. [72]

    R., & Finkbeiner, D

    Su, M., Slatyer, T. R., & Finkbeiner, D. P. 2010, ApJ, 724, 1044

  73. [73]

    E., Digel, S

    Vladimirov, A. E., Digel, S. W., J \'o hannesson, G., et al. 2011, CoPhC, 182, 1156

  74. [74]

    Tseng, P.-H., Yang, H. -Y. K., Chen, Ch.-Y. et al 2024, ApJ, 970, 146

  75. [75]

    2006, ApJ, 638, 125

    Wolfe, B., & Melia, F. 2006, ApJ, 638, 125

  76. [76]

    Yang, C. & Li, S. 2013, ApJ, 773, 138

  77. [77]

    Yang, H. Y. K., Ruszkowski, M., Ricker, P. M., Zweibel, E., & Lee, D. 2012, ApJ, 761, 185

  78. [78]

    2018, Galaxies, 6, 29

    Yang, H.-Y., Ruszkowski, M., & Zweibel, E. 2018, Galaxies, 6, 29

  79. [79]

    K., Ruszkowski, M., & Zweibel, E

    Yang, H.-Y. K., Ruszkowski, M., & Zweibel, E. G. 2022, Nature Astronomy, 6, 584

  80. [80]

    Zhou, Ye, & Matthaeus W. H. 1990, J. Geophys. Res., 95, 14