Nuclei in high-energy neutrino sources: A multimessenger study of in-source propagation

Resumen Técnico: Núcleos en Fuentes de Neutrinos de Alta Energía: Un Estudio Multimensajero de la Propagación Intra-fuente

Planteamiento del Problema
El descubrimiento de neutrinos astrofísicos de alta energía por parte de IceCube ha planteado interrogantes críticos sobre la naturaleza de sus fuentes. Mientras que los modelos estándar asumen que estos neutrinos se originan en interacciones protón-protón ($pp$) o protón-fotón ($p\gamma$), la falta de contrapartidas de rayos gamma en el rango TeV–GeV sugiere que estas fuentes son "ocultas" u opacas a los fotones de alta energía. Estudios previos han asumido mayoritariamente una composición de protones puros para las partículas aceleradas. Sin embargo, la composición química de las partículas inyectadas (por ejemplo, núcleos ligeros frente a pesados) sigue siendo una cuestión abierta. Si la fuente contiene una fracción significativa de núcleos, la propagación intra-fuente implica cascadas nucleares complejas (fotodesintegración) y procesos electromagnéticos que alteran el rendimiento resultante de neutrinos y rayos gamma. El problema central abordado es cómo, al relajar la hipótesis de protones puros y considerar la composición nuclear, se afecta la fenomenología multimensajera, específicamente la relación entre los flujos de neutrinos y rayos gamma y la luminosidad requerida de la fuente.

Metodología
Los autores realizan simulaciones detalladas de Monte Carlo de cascadas nucleares y electromagnéticas dentro de una fuente astrofísica de alta energía, utilizando la galaxia Seyfert NGC 1068 como un caso de estudio representativo. La simulación se divide en dos etapas:

1. Interacciones Fotonucleares y Cascada Nuclear: El código rastrea la evolución de núcleos primarios (desde $Z=1$ hasta $Z=90$) a través de un campo de fotones de fondo. Incluye:
   • Fotodesintegración (PD): Modelada utilizando secciones eficaces TALYS/TENDL-2023 para núcleos pesados y parametrizaciones dedicadas para núcleos ligeros ($A \le 9$). Esto incluye canales de ruptura exclusiva (p. ej., $(\gamma, n)$, $(\gamma, p)$, $(\gamma, \alpha)$) y fluorescencia de resonancia nuclear.
   • Producción de Fotomesones (PM): Modelada utilizando el generador de eventos SOPHIA con una aproximación de superposición para núcleos.
   • Producción de Pares Bethe-Heitler (BH): Tratada como un proceso de pérdida de energía continua para núcleos cargados.
   • Producción de Pares con Captura (PPC): Una inclusión novedosa de la producción de pares donde el electrón es capturado en la capa L (PPCL), seguido de la deexcitación atómica inmediata que emite fotones de línea.
   • Decaimiento Nuclear: Incluido para completitud, pero hallado como insignificante para las escalas de tiempo consideradas.
2. Cascada Electromagnética: Los fotones secundarios, electrones y positrones de la etapa (i) se propagan a través del campo de radiación. Esto incluye la Dispersión Inversa de Compton (ICS) y la Producción de Pares Electrónicos (EPP) estándar, así como procesos poco comunes como la producción de pares doble (DPP) y la producción de pares de muones (MPP), aunque este último se señala como insignificante en los rangos de energía estudiados.

La fuente se modela como una esfera homogénea de radio $R$ embebida en un campo de fotones de disco-corona (disco UV + corona de rayos X). El estudio varía el tamaño de la fuente ($R$) y la composición nuclear inyectada (protones, neutrones y diversos núcleos como Carbono, Hierro y Circonio) para mapear las profundidades ópticas y la partición de energía.

Contribuciones Clave

• Marco Numérico: El artículo presenta un nuevo tratamiento numérico de las cascadas inducidas por núcleos que extiende los procesos fotonucleares hasta isótopos de $Z=90$ y contabiliza explícitamente la producción de pares con captura en la capa L, un proceso previamente omitido en tales contextos.
• Análisis Sistemático de la Composición: Proporciona un análisis exhaustivo de cómo la composición nuclear inyectada altera la relación de rendimiento entre neutrinos y rayos gamma, demostrando que esta relación es no monotónica y altamente dependiente de la profundidad óptica de la fuente (densidad de fotones y tamaño).
• Reanálisis de COMPTEL: Los autores realizan un reanálisis dedicado de los datos de archivo de COMPTEL para NGC 1068 utilizando marcos modernos (GammaLib/ctools), derivando límites actualizados de confianza del 95% en la banda de MeV que son aproximadamente un factor de 2 más estrictos que los límites previos en los bins de menor energía.

Resultos
El estudio identifica dos regímenes distintos basados en la compacidad de la fuente (densidad de fotones):

1. Régimen Compacto/Denso ($R \lesssim 2 R_S$):

   • La fuente es ópticamente gruesa para la fotodesintegración. Los núcleos pesados se fragmentan rápidamente en nucleones libres (protones y neutrones) antes de escapar.
   • Una fracción significativa de la energía se transfiere a neutrones libres. Dado que los neutrones no sufren pérdidas de energía por Bethe-Heitler (BH), la energía disponible para la producción de neutrinos se preserva.
   • Resultado: La inyección de núcleos pesados es más eficiente en la producción de neutrinos por unidad de luminosidad de rayos cósmicos que los protones. En consecuencia, igualar el flujo de neutrinos observado por IceCube requiere una luminosidad de rayos cósmicos total menor para núcleos pesados. La cascada electromagnética resultante (rayos gamma de MeV–GeV) es más débil (en un ~28–35%) comparada con el caso de los protones, relajando así las restricciones de las observaciones de rayos gamma.

2. Régimen Extendido/Rarefacto ($R \gtrsim 5 R_S$):

   • La fotodesintegración es incompleta; una fracción sustancial de la energía inyectada permanece bloqueada en fragmentos nucleares cargados durante el cruce de la fuente.
   • Estos fragmentos pesados experimentan pérdidas de energía por BH aumentadas (escalando como $Z^2/A$) en comparación con los protones a la misma energía por nucleón.
   • Resultado: La inyección de núcleos pesados es menos eficiente en la producción de neutrinos. Para igualar el flujo de IceCube, se requiere una luminosidad de rayos cósmicos mayor (un aumento del 20–30%). Además, las pérdidas de BH aumentadas y la mayor luminosidad conducen a una cascada electromagnética significativamente más fuerte (un factor de ~2.0–2.5 superior en el flujo de rayos gamma), causando que los escenarios de núcleos pesados sobrepasen los límites existentes de rayos gamma de forma más severa que los escenarios de protones.

Significancia y Reivindicaciones
El artículo sostiene que la composición nuclear de las partículas aceleradas es un parámetro crítico, aunque a menudo pasado por alto, en el modelado de fuentes de neutrinos de alta energía. La significancia del trabajo radica en demostrar que la interpretación de los datos multimensajeros (neutrinos vs. rayos gamma) está intrínsecamente ligada a la compacidad de la fuente y a la composición inyectada.

Específicamente, los autores argumentan que:

• La banda de rayos gamma de MeV es una herramienta de diagnóstico poderosa para fuentes de neutrinos "ocultas", ya que la cascada electromagnética reprocesa los fotones de alta energía hacia esta banda.
• Los datos actuales de rayos gamma (específicamente los límites de COMPTEL reanalizados) ya proporcionan restricciones leves pero significativas sobre la composición y el tamaño de la fuente.
• En fuentes compactas, los núcleos pesados pueden explicar el flujo de neutrinos con menos potencia total y menos emisión de rayos gamma, resolviendo potencialmente la tensión con las no detecciones de rayos gamma. Por el contrario, en fuentes menos compactas, los núcleos pesados son desfavorecidos por los límites de rayos gamma.
• Las futuras instalaciones de rayos gamma de MeV (p. ej., COSI, AMEGO-X, e-ASTROGAM) serán esenciales para agudizar estas sondas y distinguir entre escenarios de inyección de protones y de núcleos en la era emergente de la astrofísicaia multimensajera.

Los autores mantienen la modestia, señalando que su modelo es un "modelo de juguete" simplificado del sistema disco-corona y que el trabajo futuro deberá incorporar pérdidas por sincrotrón y geometrías más complejas (2D/3D) para capturar plenamente la física de la fuente.