An 18-year Fermi-LAT stacking limit on GeV $γ$-ray emission from particle-accelerating colliding-wind binaries

Resumen Técnico: Un límite de apilamiento de 18 años de Fermi-LAT sobre la emisión de $\gamma$ en el rango de GeV de binarias de vientos colisionantes que aceleran partículas

Planteamiento del problema
Las binarias de vientos colisionantes (CWB, por sus siglas en inglés) masivas poseen regiones de colisión de vientos (WCR) donde los vientos estelares supersónicos colisionan, generando choques fuertes capaces de acelerar partículas a energías relativistas. Aunque la emisión de sincrotrón de radio no térmica confirma la aceleración de partículas en estos sistemas, y $\eta$ Carinae es detectada inequívocamente como una fuente de $\gamma$ en el rango de GeV, la emisión de GeV a nivel de población de las CWB permanece desconocida. Las investigaciones dedicadas previas, con exposición limitada, no encontraron emisión significativa de sistemas individuales (con la excepción de $\gamma^2$ Velorum, que es menos segura que $\eta$ Car debido a la zona congestionada de Vela). La cuestión central abordada es si la población más amplia de CWB que aceleran partículas (PACWB) produce una señal colectiva de GeV, individualmente por debajo del umbral de detección.

Metodología
Los autores analizaron 17.8 años de datos de Fermi-LAT (agosto de 2008 a mayo de 2026) por encima de 1 GeV. Se eligió este umbral de energía porque la función de dispersión de punto (PSF) se estrecha significativamente ($\lesssim 0.8^\circ$ en contención del 68%), lo que permite la separación de las binarias de sus vecinos en el plano galáctico, a diferencia de la respuesta más amplia por debajo de 1 GeV.

1. Selección de la muestra: El estudio utilizó las 61 PACWB confirmadas ("lista A") del catálogo de De Becker & Raucq.

   • Exclusiones: Los dos sistemas detectados individualmente ($\eta$ Car y $\gamma^2$ Vel) se apartaron para buscar la señal colectiva de la población no detectada.
   • Limpieza: Se aplicó una máscara estricta de 1.0$^\circ$ alrededor de fuentes brillantes en el catálogo de 16 años de Fermi-LAT (FL16Y, Signif$_{Avg} \ge 5$). Los sistemas dentro de esta máscara fueron excluidos porque sus PSF se mezclan con fuentes brillantes no relacionadas (por ejemplo, regiones de formación estelar, fuentes del Centro Galáctico), lo que hace que las mediciones de flujo por fuente sean poco fiables.
   • Fusión: Los sistemas separados por menos de 1.5$^\circ$ se fusionaron para evitar el apilamiento correlacionado, conservando el miembro con el mayor Test Statistic (TS).
   • Muestra final: Este proceso produjo una muestra "limpia" de $N=6$ sistemas aislados mutuamente en latitudes más altas.

2. Análisis de apilamiento (Stacking):

   • Se realizó un escaneo de verosimilitud bidimensional sobre el índice de fotones ($\Gamma \in [1.0, 3.5]$) y el flujo integrado ($\log_{10} F \in [-13, -6]$).
   • El Test Statistic conjunto ($TS_{stack}$) se calculó sumando las razones de verosimilitud individuales de las $N$ fuentes.
   • Campos de control: Para establecer una distribución nula, se generaron 200 campos de control. Estos se emparejaron con la muestra objetivo en latitud galáctica ($|b|$), densidad de fuentes local (suma de Signif$_{Avg}$ de las fuentes de FL16Y dentro de 2$^\circ$) e intensidad difusa galáctica ($> 1$ GeV). Crucialmente, se aplicó la misma máscara de catálogo de 1.0$^\circ$ a los campos de control para asegurar que la distribución nula tuviera en cuenta las coincidencias fortuitas con fuentes brillantes.

Resultos clave

1. No detección: El apilamiento limpio de los 6 sistemas produjo un TS acumulado de 3.91. Esto es consistente con el nulo de los campos de control emparejados por el entorno (mediana de TS = 10.0), lo que resulta en un valor p de 0.83 ($-0.9\sigma$). No hay evidencia de emisión colectiva de GeV de la población de PACWB.
2. Límites superiores: Asumiendo un índice de fotones $\Gamma = 2$, el límite superior del 95% para el flujo medio por fuente es:
   $$F(> 1 \text{ GeV}) \lesssim 1.1 \times 10^{-11} \text{ ph cm}^{-2} \text{ s}^{-1}$$
   Esto corresponde a un límite de luminosidad de $L_\gamma \lesssim 1.3 \times 10^{31} (d/\text{kpc})^2 \text{ erg s}^{-1}$.
3. Restricciones de eficiencia: Para una potencia cinética de viento representativa ($L_{wind} \sim 3.3 \times 10^{36} \text{ erg s}^{-1}$), la eficiencia de producción de $\gamma$ se restringe a:
   $$\eta = L_\gamma / L_{wind} \lesssim 4 \times 10^{-6} (d/\text{kpc})^2$$
   A la distancia mediana de la muestra ($\sim 1.7$ kpc), esto es aproximadamente dos órdenes de magnitud inferior a la eficiencia de $\eta$ Car ($\eta \sim 10^{-3}$).
4. Identificación de artefactos: Cuando se eliminó la máscara del catálogo y se apiló la población completa (59 sistemas), apareció un exceso espurio de $\gtrsim 3.5\sigma$. Los autores demuestran que esto es un artefacto de las coincidencias fortuitas entre las PACWB y las fuentes brillantes del catálogo en el plano galáctico, resaltando una advertencia metodológica crítica para el apilamiento en regiones congestionadas.

Restricciones sobre la aceleración de partículas
Utilizando un modelo de una sola zona, los autores traducen el límite de flujo en restricciones sobre la aceleración de partículas:

• Electrones: El límite limita la eficiencia de aceleración de electrones ($\eta_e$) a $\sim 10^{-3}$ para campos magnéticos ($B$) inferidos de los datos de radio (pocos mG a $\sim 1$ G). Para lograr eficiencias canónicas ($\eta_e \sim 0.01–0.1$), el campo magnético tendría que ser comparable o superior al valor de equipartición magnetismo-fotón ($B \gtrsim 3–7$ G), que es mayor que los valores típicamente inferidos para las CWB.
• Protones: El límite restringe la aceleración de protones de forma mucho más débil ($\eta_p \lesssim 0.11$) porque los protones relativistas probablemente escapan de la WCR antes de radiar mediante la desintegración $pp \to \pi^0$.

Significado y afirmaciones
El artículo concluye que $\eta$ Car parece ser un objeto singular con un poder de viento y una excentricidad extremos, en lugar de ser el miembro más brillante de una población común de CWB que emiten en el rango de GeV. La no detección implica que la mayoría de las PACWB son significativamente menos eficientes en la conversión de la potencia del viento a rayos $\gamma$ de GeV o poseen propiedades espectrales (por ejemplo, cortes) que suprimen la emisión en la banda de LAT.

El estudio enfatiza la necesidad de utilizar campos de control emparejados no solo en latitud, sino también en densidad de fuentes local e intensidad difusa al realizar apilamientos en el plano galáctico. Establece que, sin un enmascaramiento idéntico para objetivos y controles, los análisis de apilamiento en campos congestionados pueden producir falsos positivos impulsados por coincidencias con el catálogo. Las restricciones resultantes son robustas para el subconjunto de la población en latitudes altas y no contaminado, aunque los candidatos más prometedores (aquellos en regiones densas de formación estelar) permanecen sin explorar debido a la máscara necesaria.