Modelling Galactic neutrino emission: contributions from massive star clusters and interstellar cosmic rays

Resumen Técnico: Modelado de la Emisión de Neutrinos Galácticos de Cúmulos Estelares Masivos y Rayos Cósmicos Interestelares

Planteamiento del Problema
El origen de los rayos cósmicos (RC) más energéticos del Universo, particularmente aquellos que alcanzan energías de escala de PeV (Pevatrones), sigue siendo incierto. Aunque las remanentes de supernovas (SNR) se consideran tradicionalmente como los principales aceleradores, las observaciones de rayos $\gamma$ sugieren que tienen dificultades para alcanzar energías de PeV, lo que requiere clases de fuentes alternativas. Los Cúmulos Estelares Masivos Jóvenes (YMSC, por sus siglas en inglés) han sido propuestos como candidatos viables a Pevatrones debido al aumento de la turbulencia magnética generada por el choque de vientos estelares y múltiples eventos de SNR. Concurrentemente, el Observatorio IceCube ha detectado un flujo de neutrinos de alta energía desde el Plano Galáctico. Sin embargo, los análisis actuales dependen de plantillas que solo dan cuenta de la emisión de neutrinos puramente difusa generada por la interacción del mar de RC con el gas interestelar. Estos modelos no resuelven explícitamente la contribución de los aceleradores hadrónicos no resueltos, tales como la población de YMSC, lo que potencialmente enmascara las firmas específicas de estas fuentes.

Metodología
Los autores presentan un marco teórico novedoso para calcular la emisión de neutrinos esperada de la Vía Láctea, integrando dos componentes distintos:

1. Población Sintética de YMSC: Para superar la limitación de observar los YMSC solo en el vecindario solar local ($\lesssim 2$ kpc), los autores generaron 100 realizaciones independientes de la población de YMSC galáctica. Esta población sintética muestrea masas de cúmulo ($10^3 - 6.3 \times 10^4 M_\odot$), edades ($<30$ Myr) y distribuciones espaciales (siguiendo los brazos espirales, la barra Galáctica y la Espuela Local) basadas en restricciones observacionales.
2. Aceleración de Partículas y Transporte: Para cada cúmulo sintético, el modelo simula la población estelar (Secuencia Principal, Wolf-Rayet y progenitores de Supernova) para calcular la luminosidad mecánica colectiva de los vientos y las supernovas. Los autores modelan las burbujas de viento soplado resultantes, calculando los perfiles de densidad de gas y campos magnéticos.
   • Mecanismos de Aceleración: La aceleración de partículas se modela en dos sitios: el choque de terminación del viento colectivo (WTS) y los choques de remanentes de supernova (SNR).
   • Escenarios de Difusión: La propagación de los RC dentro de las burbujas de los cúmulos se computa bajo tres supuestos de turbulencia: difusión Kolmogorov, Kraichnan y Bohm.
   • Producción de Neutrinos: Los flujos de neutrinos se calculan a partir de las colisiones hadrónicas de protones acelerados con el gas objetivo (primordialmente en la densa cáscara barrida de las burbujas).
3. Componente Galáctico Difuso: El modelo también computa la emisión de neutrinos difusa que surge de la interacción del mar de RC galáctico con el gas interestelar. Esto utiliza modelos de propagación de RC (incorporando el endurecimiento espectral en el interior de la Galaxia) y mapas de gas (GALPROP y un modelo 3D alternativo).
4. Generación de Plantillas: El flujo total predicho es la suma de la contribución de los YMSC y el componente difuso inducido por los RC. Estos se sintetizan en plantillas espaciales y espectrales para su comparación con los datos de IceCube.

Contribuciones Clave

• Nuevas Plantillas de Neutrinos: El artículo proporciona las primeras plantillas de neutrinos integrales que incorporan explícitamente la emisividad de la población de YMSC junto con la componente difusa estándar inducida por los RC.
• Modelado de Población Sintética: Ofrece una estimación estadísticamente robusta de la producción acumulada de neutrinos de toda la población de YMSC galáctica, en lugar de depender únicamente de los pocos cúmulos cercanos y resueltos.
• Análisis de Múltiples Componentes: El marco distingue entre la producción "in-situ" de neutrinos dentro de los cúmulos y la emisión difusa de fondo, permitiendo una separación de las contribuciones de las fuentes.
• Sensibilidad de Parámetros: El estudio explora la dependencia de los flujos de partículas en los regímenes de difusión (Kolmogorov, Kraichnan, Bohm) y el endurecimiento espectral de los RC en el interior de la Galaxia.

Resultados

• Detectabilidad de Fuentes Individuales: Se encuentra que el flujo predicho de los YMSC individuales es marginalmente detectable incluso por detectores de escala de kilómetro cúbico (por ejemplo, IceCube, KM3NeT). En consecuencia, la emisión acumulada de la población de YMSC se espera que aparezca como una componente difusa no resuelta.
• Magnitud del Flujo: La normalización de los modelos es consistente con el mejor ajuste de IceCube de las plantillas galácticas existentes. Los autores encuentran que la contribución no resuelta de la emisión de los cúmulos puede ser no despreciable.
• Dependencia de la Difusión:
  • En el escenario de difusión de Bohm (el confinamiento más eficiente), los YMSC dominan el flujo de neutrinos por debajo de $\sim 50$ TeV, potencialmente sobreprediciendo el flujo en comparación con los datos de IceCube en los escenarios de Kraichnan y Kolmogorov si no se asume un endurecimiento espectral.
  • En los escenarios Kolmogorov y Kraichnan, la componente difusa inducida por los RC generalmente domina o proporciona contribuciones comparables a la componente de los YMSC, particularmente a energías más altas.
• Distribución Espacial: La contribución de los YMSC está concentrada en el rango de longitud galáctica interior ($|l| \le 30^\circ$, la Cresta Galáctica), mientras que la emisión inducida por los RC domina en longitudes mayores.
• Endurecimiento Espectral: Incluir el endurecimiento espectral de los RC en el interior de la Galaxia reduce la contribución relativa de los YMSC en los escenarios de Kolmogorov y Kraichnan, llevando el espectro total predicho a un mejor acuerdo con las observaciones de IceCube a $\sim 100$ TeV.
• Consistencia: Las predicciones del modelo para la contraparte de rayos $\gamma$ de las colisiones hadrónicas siguen siendo consistentes con las observaciones actuales de rayos $\gamma$ difusos en el rango de TeV-PeV.

Significancia y Reivindicaciones
Los autores afirman que su trabajo proporciona una evolución necesaria en el modelado teórico de la emisión de neutrinos galácticos. Al ir más allá de las plantillas puramente difusas para incluir la emisividad específica de la población de YMSC, ofrecen una herramienta para "desentrañar" las dos componentes de la emisión de neutrinos galácticos. El artículo postula que la aplicación de estas nuevas plantillas a futuros conjuntos de datos de neutrinos permitirá a la comunidad constreñir el papel de los cúmulos estelares en la aceleración extrema de RC y la producción de neutrinos. Los autores enfatizan que sus resultados sugieren que la contribución no resuelta de los cúmulos es probablemente no despreciable y que un ajuste conjunto de las contribuciones de las fuentes y la componente difusa es crucial para identificar los aceleradores galácticos y resolver el origen de los rayos cósmicos. El trabajo no pretende haber identificado definitivamente a los YMSC como la fuente de la señal de IceCube, sino que proporciona el marco teórico requerido para probar esta hipótesis contra los datos observacionales.