---
title: 'Fermi Bubbles: Galactic Center Outflows'
url: https://www.emergentmind.com/topics/fermi-bubbles-fbs
type: topic
---

# Fermi Bubbles: Galactic Center Outflows

The Fermi Bubbles (FBs) are two giant, bilobate, nonthermal structures centered on the Galactic Center (GC), extending to approximately $\pm 10$ kpc above and below the Galactic plane, or to angular heights of about $50^\circ$–$55^\circ$. They are defined observationally by hard gamma-ray emission with sharp edges and nearly uniform surface brightness, and they are closely associated with a microwave “haze,” polarized radio structures, and diffuse X-ray shells. Since their identification in Fermi-LAT data, the FBs have become a focal problem in Galactic feedback, cosmic-ray (CR) transport, and GC activity, with interpretations ranging from long-duration hadronic calorimetry to collimated AGN-like outbursts, star-formation-driven winds, and turbulence-driven in-situ leptonic acceleration [1008.2658] [1802.03890].

## 1. Discovery, geometry, and global phenomenology

The defining observational picture is stable across the literature. The FBs are bilateral lobes symmetric about the GC, with sharp outer edges, little internal substructure, and a remarkably uniform gamma-ray surface brightness. Their total solid angle is roughly $1\ \mathrm{sr}$, and their angular footprint is concentrated toward the GC, with the most intense gamma-ray emission at low Galactic latitudes. In physical terms, the lobes extend to heights of order $10$ kpc; in angular terms, they subtend about $50^\circ$ in latitude and several tens of degrees in width [1008.2658] [1802.03890].

Their integrated nonthermal energetics are likewise well constrained. The gamma-ray luminosity is $L_\gamma \approx 4\times 10^{37}\ \mathrm{erg\ s^{-1}}$ in the GeV band, while the microwave haze luminosity over $20$–$60$ GHz is $L_{\rm haze}\approx (1$–$5)\times 10^{36}\ \mathrm{erg\ s^{-1}}$. The diffuse X-ray luminosity associated with the same general region is at least an order of magnitude larger than the gamma-ray luminosity. This luminosity hierarchy is one of the basic empirical constraints on formation models, because any successful picture must reproduce the gamma-ray, microwave, and X-ray outputs simultaneously rather than in isolation [1008.2658].

The low-latitude base is not simply a scaled version of the high-latitude lobes. It is brighter and spectrally harder than the higher-latitude FB emission, and it is displaced to negative Galactic longitudes relative to the GC. If that base emission is physically connected to the high-latitude FBs, the westward shift disfavors strictly symmetric SMBH-driven models; alternatively, the base may include a separate component aligned along the same lines of sight [1904.01454].

## 2. Spectral structure and multiwavelength counterparts

The gamma-ray spectrum is hard, with $dN/dE \propto E^{-2}$ as a useful first description, but the detailed phenomenology is latitude-dependent. Integrated analyses find a high-energy exponential cutoff near $E_\gamma \simeq 110$–$113$ GeV, whereas the base at $|b|<10^\circ$ remains hard and extends to $\sim 1$ TeV in Fermi-LAT data without a significant cutoff; the corresponding $95\%$ confidence lower limit on the cutoff energy is about $500$ GeV. At higher latitudes, the intensity decreases and the spectrum softens or cuts off around $E\sim 100$ GeV [2411.14916] [2512.15316] [1904.01454].

Spatially resolved analyses add further structure. A 60-month reanalysis found a robust relative deficit of low-energy gamma-ray flux toward the top of the South Bubble, together with an energy-dependent morphology in which the bubbles are more extended at higher energies. That study also found no spectral softening with longitude in either lobe. Interpreting the high-latitude spectra in either hadronic or leptonic terms requires a low-energy break in the parent CR population [1402.0403].

The microwave and radio counterparts are not incidental. The WMAP/Planck haze is morphologically coincident with the FBs, especially at low and intermediate latitudes, and polarized synchrotron ridges and filaments outline the lobes and their edges. Joint IC–synchrotron fits in a leptonic framework infer a characteristic magnetic field $B \simeq 8.4\ \mu{\rm G}$ with a plausible range of $\approx 5$–$20\ \mu{\rm G}$. In a hadronic framework, secondary $e^\pm$ can account for a microwave luminosity of order $2\times 10^{36}\ \mathrm{erg\ s^{-1}}$ for $B\gtrsim 10^{-5}\ {\rm G}$, but the spectral hardness of the haze at high latitudes remains a nontrivial constraint on purely secondary models [2411.14916] [1008.2658].

The X-ray counterpart is comparably important. ROSAT revealed broad arcs and shells, and eROSITA later traced larger X-ray bubbles extending to $\approx 14$ kpc in height. Some analyses interpret the FB region as filled by hot, low-density plasma with $n_{\rm H}\approx 0.004$–$0.006\ \mathrm{cm^{-3}}$ and $T\gtrsim 3\times 10^6\ {\rm K}$, while others emphasize shell-like structures with weak or strong shock diagnostics depending on the specific tracer and location. The X-ray morphology therefore constrains both the ambient halo and the dynamical state of the FB edges [1008.2658] [2411.14916].

## 3. Multiphase gas, magnetic structure, and the nuclear wind

Radio and UV spectroscopy show that the FBs are not only nonthermal gamma-ray lobes but also a multiphase nuclear-wind system. High-velocity absorption components inside the FBs span typical LSR ranges of roughly $-300 \lesssim v_{\rm LSR} \lesssim +300\ \mathrm{km\ s^{-1}}$, while compact H I clouds in the inner $\pm 10^\circ$ reach velocities of a few hundred $\mathrm{km\ s^{-1}}$. Spectroscopic analyses infer kinematic ages of $\approx 6$–$9$ Myr for the nuclear wind, warm-ionized mass outflow rates of $\approx 0.2$–$0.4\ M_\odot\ {\rm yr^{-1}}$, neutral mass outflow rates of $\approx 0.1\ M_\odot\ {\rm yr^{-1}}$, and a conservative combined lower limit $\dot M_{\rm tot}\gtrsim 0.3$–$0.5\ M_\odot\ {\rm yr^{-1}}$ for the cool and warm phases alone. A commonly used shell estimator is
$$
\dot{M} \approx \Omega\, R\, \mu m_p\, N_{\rm H}\, v,
$$
with the hot X-ray phase likely carrying a substantial fraction of the total energy budget [1903.04833].

The cool clouds embedded in FB directions do not form a chemically homogeneous population. The first metallicity survey of FB high-velocity clouds found values ranging from less than $20\%$ solar to about $320\%$ solar. Most sampled clouds are sub-solar, but a minority are near-solar or super-solar. This directly challenges the earlier assumption that all cool FB clouds are launched from the GC with solar or super-solar metallicity. Instead, the cloud population appears to have dual origins in both the Milky Way disk and halo, implying that the expanding FB system both entrains disk material and sweeps up circumgalactic gas [2207.08838].

Magnetic structure is likewise part of the FB phenomenology rather than a secondary detail. Polarization measurements reveal a partly ordered magnetic-field geometry extending beyond the gamma-ray lobes. One MHD explosive model reproduces this by stretching and twisting a magnetic loop anchored in a rotating molecular torus around the GC, generating a large-scale spiral field consistent with the observed polarized structures. In that scenario, either confinement of gamma-ray-emitting electrons to deeper ejecta regions or consecutive explosive episodes is required to reconcile the larger radio/X-ray extent with the smaller GeV lobes [1311.6722].

## 4. Physical origin scenarios

The origin debate is not a binary choice between “AGN” and “starburst.” The literature contains several physically distinct classes of models, each reproducing some subset of the phenomenology and confronting specific tensions.

| Scenario | Core ingredients | Representative tensions or signatures |
|---|---|---|
| Saturation-based hadronic wind | Long-term GC star formation, nuclear wind, pp losses, secondaries | Requires multi-Gyr confinement; predicts neutrinos |
| Mixed hadronic–leptonic model | pp gamma rays plus primary-electron IC/synchrotron | Secondary radio alone is too soft without additional primaries |
| Routine TDE-driven turbulent shell | Repeated TDEs, RT turbulence, stochastic electron acceleration | Weak outer shock; hadronic-only microwave fit is disfavored |
| Collimated outburst / jet | Narrow, near-vertical injection, forward shocks | Isotropic injection becomes too spherical |
| Star-formation Galactic wind | Sustained central SFR, weak forward shock, in-situ electron acceleration | Requires extended episode and CGM conditions consistent with X-rays |

In the long-duration hadronic picture, CR protons and heavier ions injected by GC star formation are advected into the bubble volumes and confined for multi-Gyr timescales. The relevant transport equation is
$$
\frac{\partial N_p}{\partial t}
=
\nabla\cdot(D\nabla N_p)
-
\frac{\partial}{\partial E}[b_p(E)N_p]
+
Q_p(r,E,t),
$$
with hadronic loss time
$$
t_{pp}=\frac{1}{n\kappa \sigma_{pp} c}.
$$
For $n\approx 0.005$–$0.01\ \mathrm{cm^{-3}}$, one obtains $t_{pp}\approx (5$–$10)\times 10^9\ {\rm yr}$. In the saturation limit, $N_p \simeq \dot Q_p\, t_{pp}$ and $L_\gamma \simeq \dot Q_p/3$, so the observed luminosity follows naturally if $t_{\rm inj}\gtrsim t_{pp}$. In that framework, the FBs are multi-billion-year-old reservoirs of GC CRs, with $t_{\rm inj}\approx 8\times 10^9\ {\rm yr}$ and $n_{\rm H}\approx 0.005\ \mathrm{cm^{-3}}$ as representative values [1008.2658].

A major criticism of purely hadronic models is the radio/microwave slope. If secondary electrons cool radiatively through synchrotron and IC, their steady-state spectrum is too soft to reproduce the observed radio index; a mixed model with a primary electron component is then required. One such analysis found that, for a primary spectrum $\propto E_e^{-2}$, the permitted magnetic-field range is $2$–$7\ \mu{\rm G}$, and the pp contribution can reach about $80\%$ of the total gamma-ray flux. Outside that $B$ range, the model fails to reproduce the radio and gamma-ray data simultaneously. The same work noted that a purely hadronic solution could be recovered only if adiabatic losses dominate, but then the required magnetic field strength and CR source power become much higher than those followed from observations [1411.6395].

A different class of models treats the FBs as the product of repeated TDEs at Sgr A*. In one recent implementation, TDEs at a rate $R\simeq 10^{-5}$–$10^{-4}\ {\rm yr^{-1}}$ and energies $E_{\rm TDE}\simeq 10^{52}$–$10^{53}\ {\rm erg}$ supply an average power $P\simeq 3\times 10^{41}\ \mathrm{erg\ s^{-1}}$. The outer shock is weak, with $v_{\rm sh}\simeq 340\ \mathrm{km\ s^{-1}}$ and Mach number $\simeq 1.5$, so the observed nonthermal emission is attributed to RT-driven turbulence, MHD fluctuations via the Lighthill mechanism, and stochastic electron re-acceleration in a thin shell. That model fits the gamma-ray spectrum with $dF_\gamma/dE_\gamma \propto E_\gamma^{-1.87}\exp(-E_\gamma/113\ {\rm GeV})$ and a leptonic electron cutoff $E_{{\rm cut},e}\simeq 1.25\ {\rm TeV}$, while finding that a purely hadronic origin is disfavored because the secondary synchrotron is too soft and too faint by a factor $\gtrsim 3$–$4$ to account for the microwave haze [2411.14916].

Collimated-outburst models push in a different direction. If the FB edges are treated as strong forward shocks, isotropic or non-directed energy release becomes too round and too wide once the shock expands into the halo. In contrast, a narrow jet launched nearly perpendicular to the Galactic disc can reproduce the thin, elongated edges. One such study derived favored ballistic parameters $\theta \simeq 4^\circ$, $\beta_{-2}\gtrsim 0.4$, $E\gtrsim 2\beta_{-2}^2\times 10^{55}\ {\rm erg}$, and age $\mathbb{T}\simeq 3.3\beta_{-2}^{-1}\ {\rm Myr}$, with slowdown solutions requiring a ballistic stage reaching $z_s\gtrsim 5\ {\rm kpc}$. A related MHD explosive scenario found that a few-Myr-old event, or consecutive episodes, can also explain the large-scale magnetic structure, with total injected energy $\gtrsim 10^{55}\ {\rm erg}$ and a lifetime estimate $\approx 2\times 10^6\ {\rm yr}$ [2109.03834] [1311.6722].

Star-formation-driven wind models remain viable as well. In one hydrodynamical realization, an episode of star formation with $\mathrm{SFR}\sim 0.5\ M_\odot\ {\rm yr^{-1}}$ for $\sim 30$ Myr drives a Galactic wind of mechanical luminosity $L_w\approx 5\times 10^{40}\ \mathrm{erg\ s^{-1}}$ into a hot CGM with $n_{h,0}\approx (2.2$–$3.5)\times 10^{-3}\ \mathrm{cm^{-3}}$. The forward shock is weak, with Mach number $\approx 1.1$, and the gamma rays and microwave haze are attributed to in-situ accelerated electrons in a leptonic IC/synchrotron framework; in that model the hadronic contribution is $\lesssim 1\%$ of the observed gamma-ray signal [1505.03634].

## 5. High-energy and multimessenger diagnostics

Hadronic models generically predict neutrinos, and the predicted fluxes are highly sensitive to the proton index and cutoff. One IceCube analysis showed that, if the FB gamma rays are hadronic and the parent proton spectrum is hard, with $\alpha \approx 2.1$ and $E_{\rm cut}\gtrsim 10$ PeV, then up to about $6$–$7$ of IceCube’s $37$ three-year HESE events could originate from the FBs. In the same framework, a 10-year exposure at $E_{\rm th}\simeq 40$ TeV and $E_{\rm cut}=30$ PeV yields about $23$ FB signal events versus about $3$ atmospheric events, corresponding to a $\gtrsim 4\sigma$ excess by simple counting; for $\alpha \approx 2.3$ or lower cutoffs, the significance drops sharply [1412.6240].

Earlier hadronic estimates were even more optimistic for idealized conditions. For a proton spectrum $\propto E^{-2.0}$ with a cutoff at $\sim 1$ PeV, one prediction gave $\sim 40$ signal events per year above $10$ TeV in a Northern Hemisphere $\mathrm{km}^3$-class detector against $\sim 100$ background events, or about $4\sigma$ in one year. Any steeper spectrum or lower cutoff strongly suppresses that rate. A low neutrino flux, conversely, strengthens the case for predominantly leptonic emission [1008.2658].

The MeV band is a particularly sharp hadronic–leptonic discriminator because the two scenarios diverge strongly there. A hadronic origin predicts the broad $\pi^0$-decay feature in the $10$–$200$ MeV range, whereas a leptonic IC origin predicts a smooth continuum through the MeV gap. One forecast concluded that COSI can detect a hadronic component from the FBs at approximately $0.1$–$0.6$ MeV for magnetic fields $\lesssim 1\ \mu{\rm G}$, or set informative upper limits below about $1$ MeV. The same study found that a survey Compton–pair telescope such as AMEGO-X would detect a pure leptonic IC component across almost the entire MeV gap and could identify spectral breaks in hybrid scenarios [2111.10362].

At TeV energies, present constraints are already informative. H.E.S.S. and HAWC have not detected the FBs, and the low-latitude H.E.S.S. limits suggest a possible cutoff near $\sim 1$ TeV. CTAO simulations indicate that the low-latitude base is the highest-value target: with the $525$-hour GC survey, CTAO-South should recover the base spectrum from $\sim 60$ GeV up to $\sim 5$ TeV for an injected power law with exponential cutoff at $1$ TeV, and it should distinguish a PLEC from a pure power law at $\ge 3\sigma$ for true cutoff energies in the approximate range $\sim 1$ TeV to $\sim 150$ TeV, depending on spectral index. Even a $50$-hour exposure is expected to detect the base from $\sim 200$ GeV to $\sim 1$ TeV [2512.15316].

## 6. Current synthesis, disputes, and open problems

Several points now appear robust. The FBs are GC-centered, sharp-edged, multiwavelength structures coupled to a multiphase nuclear wind. They are associated with hot X-ray-emitting plasma, warm and neutral high-velocity clouds, synchrotron-emitting electrons, and a hard GeV gamma-ray component. Spectroscopic work has shown that the cool gas is chemically diverse rather than uniformly metal-rich, so the FB system must interact dynamically with both disk material and pre-existing halo gas [1903.04833] [2207.08838].

The central disputes concern age, engine, and dominant radiation channel. Age estimates range from $\approx 2\times 10^6$ yr in some explosive MHD scenarios, to a few Myr in collimated jet models, to $\approx 6$–$9$ Myr from spectroscopy of the nuclear wind, to $\approx 20$–$30$ Myr in some star-formation wind or X-ray-shell interpretations, and up to $\gtrsim 8\times 10^9$ yr in saturation-based hadronic calorimetry [1311.6722] [2109.03834] [1903.04833] [1505.03634] [1008.2658]. These are not minor parameter shifts; they encode fundamentally different transport regimes and different meanings of the gamma-ray luminosity.

There is also no single uncontested view of the edges. Some studies infer weak shocks, with $v_{\rm sh}\simeq 340\ \mathrm{km\ s^{-1}}$ and Mach number $\simeq 1.5$, whereas others model the edges explicitly as strong forward shocks with Mach number $\gtrsim 4$–$5$. A recent nested-bubble interpretation extends the strong-shock picture to both the FBs and the larger eROSITA bubbles, arguing that both structures can be produced by very similar collimated GC outbursts, each with $\sim 10^{55}$ erg, half-opening angle $\sim 4^\circ$, and launch speed near $2000\ \mathrm{km\ s^{-1}}$ at $\sim 100$ pc, separated by $\sim 10$ Myr. In that framework, the east–west asymmetry favors an eastern ambient-density gradient rather than a western crosswind [2411.14916] [2109.03834] [2602.00226].

A further misconception corrected by recent work is that the microwave haze automatically validates a purely leptonic or purely hadronic model. In reality, the microwave spectrum and normalization are among the strongest discriminants. Purely hadronic models can match the gamma rays only if they also solve the radio problem, and purely leptonic models must still reproduce the luminosity ratios, edge sharpness, latitude dependence, and cooling constraints without excessive fine-tuning. This is why MeV spectroscopy, TeV measurements of the base, spatially resolved X-ray spectroscopy, and neutrino searches remain decisive rather than supplementary [1411.6395] [2111.10362] [2512.15316].

The FBs therefore remain a composite problem in CR transport, shock and turbulence physics, GC feeding and feedback, and multiphase halo structure. What is established is the phenomenology: giant bilateral lobes, hard gamma rays, microwave synchrotron, X-ray shells, sharp boundaries, and chemically diverse high-velocity clouds. What remains unsettled is whether the dominant engine is a long-lived nuclear wind, repeated TDEs, a collimated AGN-like outburst, or some hybrid sequence, and whether the gamma-ray power is mainly hadronic, mainly leptonic, or mixed in an energy- and latitude-dependent way.

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