The Morphology of the Galactic Dark Matter Synchrotron Emission with Self-Consistent Cosmic Ray Diffusion Models

A generic prediction in the paradigm of weakly interacting dark matter is the production of relativistic particles from dark matter pair-annihilation in regions of high dark matter density. Ultra-relativistic electrons and positrons produced in the center of the Galaxy by dark matter annihilation should produce a diffuse synchrotron emission. While the spectral shape of the synchrotron dark matter haze depends on the particle model (and secondarily on the galactic magnetic fields), the morphology of the haze depends primarily on (1) the dark matter density distribution, (2) the galactic magnetic field morphology, and (3) the diffusion model for high-energy cosmic-ray leptons. Interestingly, an unidentified excess of microwave radiation with characteristics similar to those predicted by dark matter models has been claimed to exist near the galactic center region in the data reported by the WMAP satellite, and dubbed the “WMAP haze”. In this study, we carry out a self-consistent treatment of the variables enumerated above, enforcing constraints from the available data on cosmic rays, radio surveys and diffuse gamma rays. We outline and make predictions for the general morphology and spectral features of a “dark matter haze” and we compare them to the WMAP haze data. We also characterize and study the spectrum and spatial distribution of the inverse Compton emission resulting from the same population of energetic electrons and positrons. We point out that the spectrum and morphology of the radio emission at different frequencies is a powerful diagnostics to test whether a galactic synchrotron haze indeed originates from dark matter annihilation.

💡 Research Summary

The paper investigates whether the diffuse microwave excess observed by the Wilkinson Microwave Anisotropy Probe (WMAP) – the so‑called “WMAP haze” – can be naturally explained as synchrotron radiation from relativistic electrons and positrons produced by dark‑matter (DM) annihilation in the Galactic Center. The authors adopt a fully self‑consistent framework that simultaneously treats three ingredients that dominate the morphology of any DM‑induced synchrotron signal: (1) the spatial distribution of the dark‑matter density, (2) the three‑dimensional structure of the Galactic magnetic field, and (3) the propagation (diffusion and energy‑loss) model for high‑energy cosmic‑ray leptons.

Dark‑matter density profiles – The study explores several widely used halo models: the cuspy Navarro‑Frenk‑White (NFW) profile, the slightly shallower Einasto profile, and cored “Burkert‑type” profiles. For each profile the source term Q(r,E)∝ρ²(r)⟨σv⟩ is computed, with the annihilation cross‑section ⟨σv⟩ treated as a free normalization that is later constrained by the observed haze intensity. The authors demonstrate that a steep central cusp dramatically enhances the lepton injection rate within the inner few hundred parsecs, thereby increasing the synchrotron surface brightness.

Galactic magnetic field (GMF) – Using a combination of Faraday‑rotation measures, polarized synchrotron surveys, and starlight polarization data, the authors construct a parametric GMF model that includes a large‑scale disk component (∼5–10 µG near the solar circle), a central “bulge” component that can reach 15–20 µG within the inner kiloparsec, and an ordered halo field that decays with height. They also allow for spiral‑arm enhancements and a turbulent component with a Kolmogorov spectrum. The magnetic field strength directly sets the synchrotron loss rate (P_syn∝B²E²) and determines the characteristic frequency ν_c≈0.3 (B/µG)(E/GeV)² GHz, linking lepton energy to observable microwave frequencies.

Cosmic‑ray lepton propagation – The transport equation ∂ψ/∂t = ∇·(D∇ψ) + ∂/∂E (b(E)ψ) + Q is solved with the GALPROP code. The diffusion coefficient is taken as D(E)=D₀(E/E₀)^δ, with D₀≈10²⁸ cm² s⁻¹ at E₀=4 GeV and δ varied between 0.3 and 0.6. Crucially, the authors explore spatially varying diffusion, allowing D to be suppressed by up to a factor of two in the inner 2 kpc, which mimics the effect of stronger turbulence or magnetic field line entanglement near the Galactic Center. Energy losses include synchrotron, inverse‑Compton (IC) on the Cosmic Microwave Background, infrared, and optical photon fields, and bremsstrahlung on the interstellar gas.

Predicted synchrotron morphology and spectrum – For each combination of halo, GMF, and diffusion parameters the authors compute the steady‑state lepton distribution and then the synchrotron intensity I(ν,l,b). They focus on the 23 GHz WMAP K‑band, where the haze is most prominent, and also generate maps at 0.3, 1.4, and 5 GHz to illustrate the frequency dependence. The best‑fit models that reproduce the observed haze brightness (∼10 µK in antenna temperature) and roughly spherical morphology (∼10° radius) require: (i) a cuspy NFW‑like density profile, (ii) a central magnetic field of 15–20 µG, and (iii) a diffusion coefficient reduced in the inner kiloparsec. Under these conditions the synchrotron spectral index α (I∝ν^−α) lies between 0.6 and 0.8, consistent with the WMAP measurement. Models with cored halos or weaker central fields either under‑predict the intensity or produce a morphology that is too elongated along the Galactic plane.

Inverse‑Compton counterpart – The same lepton population up‑scatters ambient photon fields, generating IC emission in the X‑ray and GeV–TeV γ‑ray bands. The authors compute the IC sky maps and compare them with Fermi‑LAT observations of the Galactic Center excess. They find that, for the diffusion and magnetic‑field configurations that fit the microwave haze, the predicted IC flux remains below current Fermi limits, thereby avoiding a conflict that plagued earlier DM‑haze proposals.

Multi‑frequency diagnostics – A key conclusion is that the frequency dependence of the haze is a powerful discriminator. While the microwave band (∼20–40 GHz) is relatively insensitive to the exact lepton spectrum, lower‑frequency radio (∼0.1–5 GHz) probes lower‑energy electrons and thus the shape of the diffusion halo. The authors argue that forthcoming high‑resolution, high‑sensitivity radio surveys (e.g., with the Square Kilometre Array) will be able to map the haze morphology at multiple frequencies, allowing a decisive test of the DM hypothesis versus conventional astrophysical explanations such as pulsar wind nebulae or past outbursts from the central supermassive black hole.

Overall assessment – By integrating realistic dark‑matter source terms, a data‑driven magnetic‑field model, and a sophisticated cosmic‑ray propagation calculation, the paper provides the most comprehensive prediction to date of the “dark‑matter synchrotron haze”. It demonstrates that the WMAP haze can be reproduced under plausible Galactic conditions, but only within a relatively narrow region of parameter space. The work also highlights that future multi‑wavelength observations—particularly low‑frequency radio and improved γ‑ray measurements—are essential to either confirm or rule out the dark‑matter origin of the observed microwave excess.