Constraints on Dark Matter Annihilation from Starlight-Induced Inverse Compton Emission in Omega Centauri

Technical Summary: Constraints on Dark Matter Annihilation from Starlight-Induced Inverse Compton Emission in Omega Centauri

Problem Statement
The identification of the particle nature of Dark Matter (DM) remains a primary challenge in modern physics. While Weakly Interacting Massive Particles (WIMPs) are a leading candidate, indirect detection efforts have traditionally focused on high-energy gamma rays from DM annihilation or secondary radiation from interactions with the Cosmic Microwave Background (CMB). However, CMB-induced limits are often conservative. Furthermore, while attention has shifted toward the sub-GeV mass regime, conclusive signals remain elusive. This work addresses the need for more sensitive probes in dense stellar environments where the local radiation field significantly exceeds the CMB energy density, potentially enhancing secondary signals in the X-ray and gamma-ray bands.

Methodology
The authors investigate the inverse Compton (IC) scattering of electron-positron pairs ($e^+e^-$) injected by DM annihilation off the intense starlight radiation field within Omega Centauri ($\omega$ Cen), the most massive globular cluster in the Milky Way. The study proceeds through the following steps:

1. Astrophysical Modeling of $\omega$ Cen:

   • Dark Matter Distribution: The authors adopt the Navarro–Frenk–White (NFW) profile as a benchmark. To account for uncertainties in the cluster's DM content, three representative parameter sets (Profiles 1–3) from previous literature are utilized, yielding J-factors ranging from $\sim 10^{23}$ to $5 \times 10^{23} \, \text{GeV}^2 \text{cm}^{-5}$.
   • Starlight Radiation Field: A realistic spectral and spatial distribution of starlight is constructed using the PARSEC stellar population synthesis model (age 12 Gyr, metallicity [Fe/H] = -1.7) normalized to the observed V-band luminosity ($L_V \approx 1.086 \times 10^6 L_\odot$). The spatial distribution of luminosity is assumed to follow the stellar mass density. The photon number density is calculated by integrating the production rate over the cluster volume.
   • Propagation Physics: The transport of $e^+e^-$ pairs is modeled using a diffusion equation solved via the Green's function method (RX-DMFIT package). The diffusion coefficient is parameterized as $D(E) = D_0 E^{1/3}$. Energy losses include IC scattering (off both CMB and starlight), synchrotron radiation, Coulomb interactions, and bremsstrahlung. A magnetic field profile ($B_0 = 5 \, \mu\text{G}$) is included for synchrotron calculations.

2. Signal Calculation:

   • The injection rate of $e^+e^-$ is derived from the DM annihilation cross-section $\langle \sigma v \rangle$ and the DM density profile.
   • The resulting IC emissivity is calculated by convolving the steady-state electron distribution with the starlight photon field.
   • The observable flux is obtained by integrating the emissivity along the line of sight.

3. Observational Constraints:

   • Sub-GeV Regime (Chandra): For DM masses below $\sim 1$ GeV, where IC emission peaks in X-rays, the authors utilize Chandra observations. They apply a $2\sigma$ upper limit on the integrated luminosity of unresolved sources in the 0.5–4.5 keV band ($L_X < 3 \times 10^{30} \, \text{erg s}^{-1}$) to constrain the monochromatic IC luminosity at $\sim 2$ keV.
   • GeV–TeV Regime (Fermi-LAT): For masses above $\sim 1$ GeV, the authors analyze Fermi-LAT gamma-ray data. They construct a background model comprising magnetospheric emission from millisecond pulsars (MSPs) and IC emission from MSP-injected electrons. Using a $\chi^2$ profile likelihood method, they fit the total background and set 95% confidence level upper limits on the DM-induced IC component.

Key Contributions

• Starlight-Induced IC Mechanism: The paper establishes a theoretical framework specifically for IC scattering off the local starlight in a globular cluster, rather than relying solely on the CMB. It demonstrates that the starlight energy density in $\omega$ Cen exceeds the CMB by orders of magnitude, shifting the secondary emission to higher frequencies and significantly boosting the flux.
• Multi-Wavelength Approach: The study uniquely combines X-ray constraints (Chandra) for the MeV range and gamma-ray constraints (Fermi-LAT) for the GeV–TeV range to cover a continuous DM mass spectrum from MeV to TeV.
• Systematic Uncertainty Analysis: The authors quantify the impact of astrophysical uncertainties, including three distinct DM density profiles, a wide range of diffusion coefficients ($D_0 = 10^{26} \text{--} 10^{30} \, \text{cm}^2 \text{s}^{-1}$), and extreme stellar temperature scenarios.

Results

• Flux Enhancement: The predicted IC flux from starlight is approximately three orders of magnitude higher than that induced by the CMB for the same DM parameters.
• New Constraints:
  • MeV Range: In the tens-of-MeV mass range, the constraints derived from Chandra data improve upon existing limits (including CMB and X-ray limits based on CMB IC) by approximately two orders of magnitude, reaching $\langle \sigma v \rangle \sim 10^{-31} \, \text{cm}^3 \text{s}^{-1}$.
  • GeV–TeV Range: Constraints derived from Fermi-LAT data for DM masses from several GeV to the TeV scale are stronger than existing limits from CMB observations, cosmic ray measurements (Voyager), and analyses of dwarf galaxies.
• Parameter Sensitivity: The derived limits scale inversely with the J-factor (DM density). Variations in the diffusion coefficient $D_0$ cause the constraints to fluctuate by roughly two orders of magnitude; however, even under conservative assumptions (high $D_0$), the limits remain competitive with CMB-based constraints.

Significance
The paper claims that starlight-induced IC emission is a critical, previously underutilized component for indirect DM searches in dense stellar systems. By leveraging the unique environment of Omega Centauri, the study demonstrates that globular clusters can serve as powerful probes for DM, particularly in the sub-GeV to MeV regime where traditional gamma-ray searches are less sensitive. The results suggest that $\omega$ Cen provides stringent limits that surpass current benchmarks, establishing it as a prime target for future DM investigations. The work underscores the necessity of incorporating realistic stellar radiation fields into DM propagation models to fully exploit the potential of dense stellar environments.