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

Resumen Técnico: Restricciones a la aniquilación de materia oscura mediante la emisión de Compton inverso inducida por la luz estelar en Omega Centauri

Planteamiento del problema
La identificación de la naturaleza de partícula de la Materia Oscura (DM, por sus siglas en inglés) sigue siendo un desafío primordial en la física moderna. Si bien las Partículas Masivas que Interactúan Débilmente (WIMPs) son un candidato principal, los esfuerzos de detección indirecta se han centrado tradicionalmente en los rayos gamma de alta energía provenientes de la aniquilación de la DM o en la radiación secundaria resultante de interacciones con el Fondo Cósmico de Microondas (CMB). Sin embargo, los límites inducidos por el CMB suelen ser conservadores. Además, aunque la atención se ha desplazado hacia el régimen de masa sub-GeV, las señales concluyentes siguen siendo esquivas. Este trabajo aborda la necesidad de sondas más sensibles en entornos estelares densos donde el campo de radiación local excede significativamente la densidad de energía del CMB, potenciando potencialmente las señales secundarias en las bandas de rayos X y gamma.

Metodología
Los autores investigan la dispersión Compton inverso (IC) de los pares electrón-positrón ($e^+e^-$) inyectados por la aniquilación de la DM contra el intenso campo de radiación de la luz estelar dentro de Omega Centauri ($\omega$ Cen), el cúmulo globular más masivo de la Vía Láctea. El estudio procede a través de los siguientes pasos:

1. Modelado astrofísico de $\omega$ Cen:

   • Distribución de Materia Oscura: Los autores adoptan el perfil Navarro–Frenk–White (NFW) como referencia. Para dar cuenta de las incertidumbres en el contenido de DM del cúmulo, se utilizan tres conjuntos de parámetros representativos (Perfiles 1–3) de la literatura previa, lo que arroja factores J que oscilan entre $\sim 10^{23}$ y $5 \times 10^{23} \, \text{GeV}^2 \text{cm}^{-5}$.
   • Campo de radiación de la luz estelar: Se construye una distribución espectral y espacial realista de la luz estelar utilizando el modelo de síntesis de población estelar PARSEC (edad 12 Gyr, metalicidad [Fe/H] = -1.7) normalizado a la luminosidad en la banda V observada ($L_V \approx 1.086 \times 10^6 L_\odot$). La distribución espacial de la luminosidad se asume que sigue la densidad de masa estelar. La densidad numérica de fotones se calcula integrando la tasa de producción sobre el volumen del cúmulo.
   • Física de propagación: El transporte de los pares $e^+e^-$ se modela mediante una ecuación de difusión resuelta mediante el método de la función de Green (paquete RX-DMFIT). Las pérdidas de energía incluyen la dispersión IC (contra el CMB y la luz estelar), radiación de sincrotrón, interacciones de Coulomb y bremsstrahlung. Se incluye un perfil de campo magnético ($B_0 = 5 \, \mu\text{G}$) para los cálculos de sincrotrón.

2. Cálculo de la señal:

   • La tasa de inyección de $e^+e^-$ se deriva de la sección eficaz de aniquilación de la DM $\langle \sigma v \rangle$ y del perfil de densidad de la DM.
   • La emisividad de IC resultante se calcula convolucionando la distribución de electrones en estado estacionario con el campo de fotones de la luz estelar.
   • El flujo observable se obtiene integrando la emisividad a lo largo de la línea de visión.

3. Restricciones observacionales:

   • Régimen sub-GeV (Chandra): Para masas de DM por debajo de $\sim 1$ GeV, donde la emisión de IC alcanza su máximo en rayos X, los autores utilizan observaciones de Chandra. Aplican un límite superior de $2\sigma$ a la luminosidad integrada de fuentes no resueltas en la banda de 0.5–4.5 keV ($L_X < 3 \times 10^{30} \, \text{erg s}^{-1}$) para restringir la luminosidad de IC monocromática a $\sim 2$ keV.
   • Régimen GeV–TeV (Fermi-LAT): Para masas por encima de $\sim 1$ GeV, los autores analizan datos de rayos gamma de Fermi-LAT. Construyen un modelo de fondo que comprende la emisión magnetosférica de púlsares de milisegundo (MSPs) y la emisión IC de electrones inyectados por los MSP. Utilizando un método de verosimilitud de perfil $\chi^2$, ajustan el fondo total y establecen límites superiores al nivel de confianza del 95% sobre el componente de IC inducido por la DM.

Contribuciones clave

• Mecanismo de IC inducido por la luz estelar: El artículo establece un marco teórico específicamente para la dispersión IC contra la luz estelar local en un cúmulo globular, en lugar de depender únicamente del CMB. Demuestra que la densidad de energía de la luz estelar en $\omega$ Cen supera al CMB por órdenes de magnitud, desplazando la emisión secundaria hacia frecuencias más altas y aumentando significativamente el flujo.
• Enfoque multi-longitud de onda: El estudio combina de manera única las restricciones de rayos X (Chandra) para el rango de MeV y las de rayos gamma (Fermi-LAT) para el rango de GeV–TeV para cubrir un espectro de masa de DM continuo desde el MeV hasta el TeV.
• Análisis de incertidumbre sistemática: Los autores cuantifican el impacto de las incertidumbres astrofísicas, incluyendo tres perfiles distintos de densidad de DM, un amplio rango de coeficientes de difusión ($D_0 = 10^{26} \text{--} 10^{30} \, \text{cm}^2 \text{s}^{-1}$) y escenarios extremos de temperatura estelar.

Resultados

• Mejora del flujo: El flujo de IC predicho por la luz estelar es aproximadamente tres órdenes de magnitud mayor que el inducido por el CMB para los mismos parámetros de DM.
• Nuevas restricciones:
  • Rango de MeV: En el rango de masas de decenas de MeV, las restricciones derivadas de los datos de Chandra mejoran los límites existentes (incluyendo los límites de CMB y de rayos X basados en IC del CMB) en aproximadamente dos órdenes de magnitud, alcanzando $\langle \sigma v \rangle \sim 10^{-31} \, \text{cm}^3 \text{s}^{-1}$.
  • Rango GeV–TeV: Las restricciones derivadas de los datos de Fermi-LAT para masas de DM desde varios GeV hasta la escala de TeV son más fuertes que los límites actuales de observaciones del CMB, mediciones de rayos cósmicos (Voyager) y análisis de galaxias enanas.
• Sensibilidad de parámetros: Los límites derivados escalan inversamente con el factor J (densidad de DM). Las variaciones en el coeficiente de difusión $D_0$ causan que las restricciones fluctúen aproximadamente dos órdenes de magnitud; sin embargo, incluso bajo supuestos conservadores (alto $D_0$), los límites siguen siendo competitivos con las restricciones basadas en el CMB.

Significancia
El artículo sostiene que la emisión de IC inducida por la luz estelar es un componente crítico y previamente subutilizado para las búsquedas indirectas de DM en sistemas estelares densos. Al aprovechar el entorno único de Omega Centauri, el estudio demuestra que los cúmulos globulares pueden servir como potentes sondas para la DM, particularmente en el régimen sub-GeV a MeV donde las búsquedas tradicionales de rayos gamma son menos sensibles. Los resultados sugieren que $\omega$ Cen proporciona límites estrictos que superan los puntos de referencia actuales, estableciéndolo como un objetivo primordial para futuras investigaciones de DM. El trabajo subraya la necesidad de incorporar campos de radiación estelar realistas en los modelos de propagación de la DM para explotar plenamente el potencial de los entornos estelares densos.