---
title: 'Fermi Bubbles: Galactic Gamma-Ray Lobes'
url: https://www.emergentmind.com/topics/fermi-bubbles
type: topic
---

# Fermi Bubbles: Galactic Gamma-Ray Lobes

The Fermi Bubbles are two giant gamma-ray lobes centered on the Galactic Center, extending up to \(55^\circ\) above and below the Galactic plane, or roughly \(10\ \mathrm{kpc}\) into each Galactic hemisphere. They are characterized by spectrally hard gamma-ray emission, sharp edges, and approximate bilateral symmetry, and they are spatially associated with the microwave haze, polarized radio lobes, and large-scale X-ray structures. Their existence is firmly established, but both the dominant radiation mechanism and the driving engine—past activity of Sgr A\(^*\), nuclear star formation, or repeated smaller outbursts—remain debated [1508.06592] [1903.04833] [1802.03890].

## 1. Discovery and defining observational properties

The modern picture of the bubbles is explicitly multiwavelength. The central review of spectroscopic work notes radio detection as early as Sofue & Handa (1984), X-ray and mid-IR detection by Bland-Hawthorn & Cohen (2003), the microwave “WMAP haze” identified by Finkbeiner (2004), and gamma-ray identification by Su et al. (2010) and Dobler et al. (2010) [1903.04833]. In gamma rays, the structures appear as two large, homogeneous regions of spectrally hard emission, centered on the Galactic Center and extending to about \(\pm 55^\circ\) in latitude [1508.06592].

The gross morphology is unusually simple for such a large Galactic structure. The bubbles are approximately symmetric about the Galactic plane, have smooth large-scale surfaces, and show sharp boundaries rather than diffuse fading into the halo. Their total solid angle is about \(1\ \mathrm{sr}\), and their gamma-ray surface brightness is approximately uniform up to \(|b| \approx 55^\circ\) [1802.03890] [1704.02629]. These properties immediately suggest a nuclear origin and make the bubbles a nearby laboratory for Galactic feedback, cosmic-ray transport, and magnetized outflows [1903.04833].

## 2. Multiwavelength counterparts and large-scale environment

A major empirical anchor of bubble studies is the close association between the gamma-ray lobes and the microwave haze. In a Fermi-LAT analysis of 50 months of data, the inverse-Compton interpretation was shown to be compatible with the WMAP and Planck haze if the same electron population produces synchrotron emission in a magnetic field \(B = 8.4 \pm 0.2\,{\rm [stat]}^{+11.2}_{-3.5}\,{\rm [syst]}\ \mu{\rm G}\), with an allowed approximate range of \(5\)–\(20\,\mu{\rm G}\) [1407.7905]. This connection is central because the gamma rays mainly probe electrons above \(\sim 100\) GeV, whereas the microwave haze probes mainly \(10\)–\(30\) GeV electrons [1407.7905].

The radio environment is broader than the GeV morphology alone. Polarized radio lobes seen in S-PASS extend to approximately \(\pm 60^\circ\), and the associated structures include the North and South ridges, the Galactic Center spur, and limb-brightened spurs [1903.04833]. Several MHD studies interpret this as evidence for a partly ordered field threading a larger bubble complex than the gamma-ray-defined lobes [1311.6722] [1311.6159].

At X-ray energies, the bubble context is again larger. The 2021 MeV study notes that the eROSITA X-ray bubbles are larger, more spherical, and edge-brightened relative to the gamma-ray bubbles, with average \(0.6\)–\(1\) keV surface brightness \((2-4)\times 10^{-15}\,\mathrm{erg\,cm^{-2}\,s^{-1}\,arcmin^{-2}}\), or roughly \((1.5-3)\times 10^{4}\,\mathrm{eV\,cm^{-2}\,s^{-1}\,sr^{-1}}\), nearly two orders of magnitude above the Fermi Bubble gamma-ray intensity [2111.10362]. This disparity suggests that the X-ray, radio, microwave, and gamma-ray structures need not trace exactly the same plasma phase or acceleration epoch.

## 3. Gamma-ray spectrum, substructure, and high-energy constraints

A benchmark spectral analysis used 50 months of Fermi-LAT data between 100 MeV and 500 GeV above \(|b|>10^\circ\) and found that the bubble spectrum is not well described by a simple power law. It is well fit either by a log parabola or by a power law with exponential cutoff,
$$
\frac{dN}{dE} = I \left(\frac{E}{1\,{\rm GeV}}\right)^{-\gamma} e^{-E/E_{\rm cut}},
$$
with
$$
\gamma = 1.87 \pm 0.02 \,{\rm [stat]}^{+0.14}_{-0.17}\,{\rm [syst]},\qquad
E_{\rm cut} = 113 \pm 19 \,{\rm [stat]}^{+45}_{-53}\,{\rm [syst]}\ {\rm GeV}.
$$
The abstract summarizes this as index \(1.9\pm 0.2\) and cutoff \(110\pm 50\) GeV, and a simple power law is excluded at more than \(7\sigma\) significance [1407.7905].

The same study measured a total gamma-ray luminosity for \(|b|>10^\circ\) and 100 MeV–500 GeV of
$$
L_\gamma = (4.4 \pm 0.1\,{\rm [stat]}^{+2.4}_{-0.9}\,{\rm [syst]}) \times 10^{37}\ {\rm erg\,s^{-1}},
$$
and estimated a median boundary width
$$
\Delta = 3.4 \pm 2.0 \,{\rm [stat]}^{+3.1}_{-1.7}\,{\rm [syst]} \ {\rm deg},
$$
reported in the abstract as \(3.4^{+3.7}_{-2.6}\ {\rm deg}\) [1407.7905]. It also confirmed a south-eastern “cocoon,” with \(\mathrm{TS}\) values between 95 and 975 depending on foreground model, but found no significant evidence for a jet [1407.7905].

Spatial-spectral uniformity remains partly unsettled. The 50-month analysis reported no significant spectral variation across the bubbles [1407.7905]. By contrast, a later 60-month analysis found, robustly, a relative deficit of the flux at low energies toward the top of the South Bubble, and also found that at high energies the bubbles are more extended [1402.0403]. The same work stressed that the detailed spectral behavior of each slice is somewhat dependent on the assumed background model, and that the North Bubble is much less secure because of Loop I confusion [1402.0403].

At TeV energies, non-detections have become an important constraint. The 2021 MeV study summarizes that HAWC found no significant northern-bubble emission above 1.2 TeV, and H.E.S.S. found no significant low-latitude TeV emission, constraining the cutoff to \(700\) GeV–\(1.7\) TeV [2111.10362]. An earlier HAWC conference proceeding had already framed the logic explicitly: TeV observations are crucial because “a steeper cutoff will favor a leptonic model” [1508.06592].

## 4. Emission mechanisms and spectral diagnostics

In the leptonic inverse-Compton picture, the gamma rays are produced by relativistic electrons scattering the interstellar radiation field. A gamma-ray fit to the Fermi-LAT data implies an electron spectrum
$$
\frac{dn_e}{dE} \propto E^{-n} e^{-E/E_{\rm cut}},
$$
with
$$
n = 2.17 \pm 0.05\,{\rm [stat]}^{+0.33}_{-0.89}\,{\rm [syst]},\qquad
E_{\rm cut} = 1.25 \pm 0.13\,{\rm [stat]}^{+1.73}_{-0.68}\,{\rm [syst]}\ {\rm TeV},
$$
and a total electron energy above 1 GeV
$$
W_e = (1.0 \pm 0.2\,{\rm [stat]}^{+6.0}_{-1.0}\,{\rm [syst]}) \times 10^{52}\ {\rm erg}.
$$
This scenario fits the gamma-ray data very well and also naturally connects to the microwave haze [1407.7905]. Its main physical difficulty is cooling: for 1 TeV electrons in a \(5\,\mu{\rm G}\) field and the ISRF at 5 kpc, the cooling time is only about \(5\times 10^5\ {\rm yr}\), or \(\sim 1\) Myr if only IC losses are considered, implying either very fast transport or in-situ reacceleration [1407.7905].

In the hadronic picture, gamma rays arise from inelastic proton-proton collisions and \(\pi^0\)-decay. The corresponding proton spectrum is
$$
\frac{dn(p)}{dp} \propto p^{-n} e^{-pc/E_{\rm cut}},
$$
with
$$
n = 2.13 \pm 0.01\,{\rm [stat]}^{+0.15}_{-0.52}\,{\rm [syst]},\qquad
E_{\rm cut} = 14 \pm 7\,{\rm [stat]}^{+6}_{-13}\,{\rm [syst]}\ {\rm TeV}.
$$
Assuming \(n_{\rm H}=0.01\ {\rm cm^{-3}}\), the required CR proton energy above 1 GeV is
$$
W_p = (3.5 \pm 0.1\,{\rm [stat]}^{+4.7}_{-3.0}\,{\rm [syst]}) \times 10^{55} \left(\frac{0.01\ {\rm cm^{-3}}}{n_{\rm H}}\right)\ {\rm erg},
$$
with proton-proton collision timescale
$$
t_p \approx 3.5\times 10^9 \left(\frac{0.01\ {\rm cm^{-3}}}{n_{\rm H}}\right)\ {\rm yr}.
$$
Primary hadronic gamma rays alone fit the data less well than IC, with IC preferred over primary-only hadronic models at at least \(4.9\sigma\), but hadronic models including IC emission from secondary leptons can also fit the gamma-ray data well [1407.7905].

The principal discriminant is no longer the GeV spectrum alone. The same 2014 analysis emphasized that a pure hadronic explanation of both gamma rays and the microwave haze does not work well, because synchrotron from secondary \(e^\pm\) is too faint by a factor of roughly 3–4 and too soft spectrally [1407.7905]. The 2018 origin review made the same point more generally: purely hadronic models fail to reproduce the microwave haze without an additional population of primary electrons [1802.03890].

The MeV band has therefore become strategically important. In the 2021 study, the benchmark leptonic model remains bright across roughly \(100\,\mathrm{keV}\) to \(100\,\mathrm{MeV}\), whereas the hadronic benchmark is weaker and strongly magnetic-field dependent. COSI is projected to detect a hadronic component in the range \(0.1\)–\(0.6\) MeV during the prime mission if the bubble magnetic field is \(\lesssim 1\,\mu\mathrm{G}\), while an AMEGO-X-like mission would detect a pure leptonic component through almost the entire MeV gap, roughly \(100\,\mathrm{keV}\)–\(100\,\mathrm{MeV}\) [2111.10362].

## 5. Gas kinematics, multiphase structure, and magnetic organization

Imaging establishes the bubbles’ shape, but spectroscopy has provided the strongest direct evidence that they are also a multiphase nuclear wind. Radio spectroscopy shows a deficiency of H I within 2 kpc of the Galactic Center and, inside this cavity, about 200 compact H I clouds. These clouds lie within roughly \(\pm 10^\circ\) of the Galactic plane, have temperatures of about \(8000\ \mathrm{K}\), typical neutral hydrogen column densities of about \(10^{19}\ \mathrm{cm^{-2}}\), and are interpreted as neutral clumps entrained in a biconical nuclear wind. Existing radio analyses give a neutral-gas mass outflow rate of about \(0.1\ M_\odot\,\mathrm{yr^{-1}}\) [1903.04833].

Ultraviolet spectroscopy reveals the ionized phases of the same flow. Warm-ionized tracers at \(T\sim 10^4\ \mathrm{K}\) include O I, N I, C II, Si II, Si III, S II, and Fe II, while C IV, Si IV, and N V trace gas at \(T\sim 10^5\ \mathrm{K}\). High-velocity absorption is more common along sightlines through the bubbles than outside them, and a foreground-background stellar pair demonstrates absorbing gas beyond \(7\ \mathrm{kpc}\), with additional absorption between \(-260\) and \(-60\ \mathrm{km\,s^{-1}}\) in the background star alone. Current UV studies imply a kinematic age of roughly \(6\)–\(9\ \mathrm{Myr}\) and warm-ionized mass outflow rates of \(\sim 0.2\)–\(0.4\ M_\odot\,\mathrm{yr^{-1}}\), but whether the wind accelerates or decelerates with distance remains unresolved [1903.04833].

Magnetic structure is a second major theme. A 3D MHD study showed that an explosion from a few million years ago can push and shear a surrounding magnetic loop anchored in the molecular torus around the Galactic Center, generating a spiral, partly ordered magnetic structure around the bubble region. In that framework, the bubble lifetime is estimated as \(2\times 10^6\) yr, with total injected energy \(>10^{55}\) erg, and consecutive explosive events may match the morphology better than a single one [1311.6722]. A related AGN-jet MHD model found that a single population of cosmic-ray leptons can simultaneously explain the gamma-ray and microwave spectra, while a second, more recent pair of jets embedded in the bubbles helps reproduce the centrally peaked microwave emission; that model also predicts strongly polarized radio emission and a layer of enhanced rotation measure in the shock-compressed region [1311.6159].

## 6. Formation scenarios and unresolved questions

The formation literature spans a wide dynamic range in age, power source, and transport physics. One class of models treats the bubbles as long-lived hadronic reservoirs supplied by sustained star formation in the Galactic Center. In that picture, the relevant timescale is \(\gtrsim 8\) Gyr, the characteristic bubble density is \(n_H \simeq 0.005\ {\rm cm}^{-3}\), the star formation rate is \(\simeq 0.08\,M_\odot\,{\rm yr}^{-1}\), and the GC injects \(\sim 1.3\times 10^{39}\ {\rm erg\,s^{-1}}\) in cosmic rays; the bubbles then act as calorimeters for a hard proton population [1008.2658].

A second class invokes repeated stellar capture or tidal disruption by Sgr A\(^*\). Early versions adopted capture rates of \(3\times 10^{-5}\ {\rm yr}^{-1}\) or \(\lesssim 10^{-5}\ {\rm yr}^{-1}\), with \(\sim 3\times 10^{52}\) erg or \(\sim 10^{52}\) erg per event, heating Galactic Center gas to \(\sim 10\) keV, launching a wind with velocity \(\sim 10^8\ {\rm cm\,s^{-1}}\), and generating many shocks that accelerate electrons to \(\sim\)TeV or tens of TeV [1103.1002] [1109.6087]. A 2024 synthesis retained the repeated-TDE engine, estimating ten to hundred events per million years and an average power of \(3\times 10^{41}\ {\rm erg\,s^{-1}}\), then coupled it to Kompaneets-type expansion, Rayleigh-Taylor disruption of the shell, hydrodynamic turbulence in the envelope, Lighthill-generated MHD fluctuations, and stochastic re-acceleration of electrons; in that framework the preferred nonthermal emission is again leptonic [2411.14916].

A third class minimizes the need for a special accelerator by treating the bubbles as passive reservoirs illuminated by ordinary Galactic cosmic rays. In the diffusive-injection model, the gamma rays are produced by Galactic CR protons leaking into slowly expanding bubbles, with an energy-independent diffusion coefficient inside the bubbles
$$
D_{\rm b}\simeq 2.6\times 10^{27}\ \mathrm{cm^2\,s^{-1}},
$$
best-fit \(K=0.26\), and expansion speed constrained by the spectrum to
$$
U<180\ \mathrm{km\,s^{-1}},
$$
implying
$$
t_{\rm age}>2.44\times 10^7\ \mathrm{yr}.
$$
That model reproduces the observed intensity profile and the integrated spectrum between \(\sim 1\) and \(100\) GeV without invoking additional particle production processes [1304.6972].

A fourth class places acceleration inside the bubbles. In the stochastic-acceleration model, electrons undergo volume-filling 2nd-order Fermi acceleration by plasma turbulence, with characteristic parameters \(L=2\,\mathrm{kpc}\), \(U=2.6\times 10^8\,\mathrm{cm\,s^{-1}}\), and total energy in electrons above 100 MeV of \(\sim 10^{51}\,\mathrm{erg}\); the model predicts nearly constant surface brightness, sharp edges, and limb-brightening at \(\sim 500\) GeV, while matching the microwave haze requires \(B\sim 15\,\mu\mathrm{G}\) rather than the fiducial \(4\,\mu\mathrm{G}\) [1104.3585]. A related hadronic “scaled-up SNR” model treats the bubbles as forward shocks in the halo, with a favored bubble age of \(10^7\) yr and the specific prediction that the apparent bubble size should be larger in the TeV band than in the GeV band [1308.5228].

Purely geometrical and shell-dynamical interpretations also exist. In a superbubble treatment, a thin-shell momentum-conserving model propagating in a medium with inverse-square density decline gave a stated reliability of \(\epsilon_{\rm obs}=90.71\%\) when compared with a digitized bubble section, outperforming the thermal and self-gravitating reference models [1806.09092]. This suggests that some large-scale aspects of the morphology can be captured without specifying the nonthermal particle physics in detail.

The field’s unresolved issues are therefore not the existence or basic phenomenology of the bubbles, but the coupling between dynamics, transport, and radiation. A 2018 review concluded that purely hadronic models are incomplete because they do not explain the microwave haze without extra primary electrons; that the simplest in-situ acceleration models struggle with the flat projected intensity profile; and that leptonic jet models currently reproduce the broad gamma-ray, microwave, X-ray, and radio-polarization phenomenology particularly well, but require past accretion rates \(\gtrsim 10\%\) of Eddington [1802.03890]. Future progress is expected from denser UV and radio spectroscopy, X-ray/UV metallicity and kinematic measurements, deeper TeV limits, neutrino searches, and especially MeV observations, where leptonic and hadronic spectra diverge most strongly [1903.04833] [2111.10362].

Source: https://www.emergentmind.com/topics/fermi-bubbles