Chasing the neutrino blazar candidates II: SED modeling with hadronic model

Resumen Técnico: Persiguiendo candidatos a blázares de neutrinos II: Modelado SED con modelo hadrónico

Enunciado del Problema
Aunque los blázares se consideran candidatos prometedores para fuentes de neutrinos de alta energía, el origen físico de su emisión de neutrinos sigue siendo incierto. Aunque el escenario leptónico (dispersión Compton inversa) es el marco estándar para modelar las distribuciones de energía espectral (SED) de los blázares, la detección de neutrinos de alta energía sugiere un papel significativo para los hadrones relativistas. Un trabajo anterior de los autores (Documento I) identificó 199 candidatos a blázares de neutrinos (NBC) utilizando un enfoque de aprendizaje con pocos ejemplos. Este estudio tiene como objetivo extender ese trabajo investigando las propiedades físicas de estos candidatos mediante modelado SED de banda ancha dentro de un marco hadrónico, específicamente para estimar la producción máxima posible de neutrinos y restringir la composición de partículas de los chorros de los blázares.

Metodología
Los autores construyeron SEDs de banda ancha para 103 NBC (seleccionados de los 199 originales basándose en una cobertura multibanda suficiente) y la fuente confirmada TXS 0506+056. El modelado empleó un código hadrónico dependiente del tiempo de una sola zona (AM3) con las siguientes suposiciones y procedimientos clave:

• Dominancia Hadrónica: Para estimar el límite superior de la emisión de neutrinos, el modelo asume que los fotones de alta energía (HE) y muy alta energía (VHE) se producen principalmente a través de interacciones $p\gamma$ y las cascadas electromagnéticas asociadas. La contribución de la dispersión Compton inversa (IC) leptónica se suprime fuertemente.
• Configuración Física: La región de emisión se modela como una mancha esférica que se mueve con un factor de Lorentz de masa $\Gamma$ en un campo magnético $B$. Electrones y protones no térmicos se inyectan con espectros de ley de potencia.
• Procesos Modelados: El modelo considera la radiación sincrotrón, la producción de fotopiones ($p\gamma \to n\pi^+$ o $p\pi^0$), la producción de pares Bethe-Heitler ($p\gamma \to pe^+e^-$), la producción de pares fotón-fotón y la desintegración subsiguiente de piones y muones para producir neutrinos. Se aplica la atenuación de la luz de fondo extragaláctica (EBL) a los rayos gamma de alta energía.
• Estrategia de Ajuste: Los autores utilizaron el algoritmo de Optimización por Enjambre de Partículas (PSO) a través del paquete PySwarms para ajustar nueve parámetros libres (por ejemplo, intensidad del campo magnético, tamaño de la región, luminosidades de inyección de partículas e índices espectrales) contra datos multibanda que abarcan desde radio hasta bandas de rayos $\gamma$. Los datos de radio por debajo de $10^{11}$ Hz se excluyeron debido a limitaciones de autoabsorción sincrotrón en el modelo de una sola zona. Se aplicó un refinamiento de "ajuste a simple vista" para asegurar representaciones SED aceptables.
• Análisis de Correlación: Se realizó un análisis de correlación parcial para investigar las relaciones entre la luminosidad de neutrinos ($L_\nu$) y las luminosidades electromagnéticas (radio, óptico, rayos X, rayos $\gamma$), así como las propiedades del motor central (masa del agujero negro $M_{BH}$ y luminosidad del disco de acreción $L_{disk}$), controlando la dependencia del corrimiento al rojo.

Contribuciones y Resultados Clave

1. Modelado SED y Restricciones de Parámetros: El estudio restringió con éxito nueve parámetros clave para 103 NBC. El modelado revela que la intensidad del campo magnético tiende a ser alta (el 66% de las fuentes tienen $\log B \ge 1.0$) para disipar eficientemente la energía de los electrones a través de la radiación sincrotrón, suprimiendo así la emisión IC y permitiendo que los protones dominen el pico de alta energía. La potencia transportada por protones relativistas ($P_p$) domina tanto sobre la potencia de electrones ($P_e$) como sobre la potencia del campo magnético ($P_B$), con una relación promedio $\log(P_p/P_B) \simeq 1.67$.
2. Predicción de la Banda MeV: Una característica distintiva del modelo es la predicción de una emisión sincrotrón de protones prominente que alcanza su máximo en la banda MeV. Para 99 de los 103 NBC, el pico sincrotrón de protones cae dentro de 0.1–100 MeV. En 90 de estas fuentes, el flujo sincrotrón de protones en la banda MeV supera el flujo de producción de pares Bethe-Heitler en la banda GeV.
3. Correlaciones con la Emisión Electromagnética:
   • Rayos Gamma y Rayos X: Se encontraron correlaciones positivas fuertes entre la luminosidad de neutrinos y las luminosidades de rayos $\gamma$ ($r=0.78$) y rayos X ($r=0.52$). Esto apoya un origen hadrónico donde los protones de alta energía producen tanto neutrinos como fotones de alta energía.
   • Óptico (Banda R): Existe una correlación débil/moderada ($r=0.30$) entre la luminosidad de neutrinos y la luminosidad óptica en la banda R. Los autores sugieren que esto puede surgir porque los fotones sincrotrón ópticos del chorro sirven como fotones semilla para la producción de fotopiones, o potencialmente debido a una conexión con el disco de acreción.
   • Radio: No se encontró ninguna correlación significativa entre la luminosidad de neutrinos y la luminosidad de radio a 1.4 GHz.
   • Motor Central: No se encontró ninguna correlación estadísticamente significativa entre la luminosidad de neutrinos y la masa del agujero negro ($M_{BH}$) o la luminosidad del disco de acreción ($L_{disk}$), aunque visualmente era aparente una tendencia positiva débil con $L_{disk}$.
4. Perspectivas de Detección: Basado en los flujos máximos de neutrinos predichos por el modelo:
   • IceCube: 3 NBC son potencialmente detectables.
   • Instalaciones Futuras: Hasta 22 fuentes podrían ser detectadas por KM3NeT, 45 por NEON y 62 por TRIDENT.
   • El flujo de neutrinos sumado de la muestra no supera los límites actuales de apilamiento de IceCube, lo que sugiere que el escenario de "máximo neutrino" es físicamente viable.

Significado y Afirmaciones
El artículo afirma que sus resultados proporcionan predicciones comprobables para futuras observaciones de mensajería múltiple y ofrecen nuevas perspectivas sobre la composición de los chorros de los blázares. Específicamente:

• La Ventana MeV: Los autores destacan la banda MeV como una ventana observacional crítica para distinguir entre escenarios leptónicos y hadrónicos, ya que el modelo hadrónico predice un pico sincrotrón de protones distinto en este rango que está ausente en los modelos leptónicos estándar.
• Composición del Chorro: El modelado respalda un escenario donde los protones relativistas dominan el presupuesto energético de los chorros en estos candidatos, consistente con la evidencia de mensajería múltiple de TXS 0506+056.
• Límites Superiores: Los resultados representan un límite superior en la emisión de neutrinos; los flujos reales pueden ser más bajos, lo que explica por qué muchos candidatos aún no han sido confirmados como emisores de neutrinos a pesar de estar por encima de las curvas de sensibilidad actuales en el modelo.
• Observaciones Futuras: El estudio enfatiza que las misiones de rayos $\gamma$ de MeV de próxima generación (por ejemplo, e-ASTROGAM, AMEGO) y los detectores de neutrinos de próxima generación (KM3NeT, NEON, TRIDENT) son esenciales para verificar estos escenarios hadrónicos y restringir la composición de partículas de los chorros de los blázares.

Los autores permanecen modestos respecto al origen de la correlación óptico-neutrino, señalando que podría derivar de efectos de selección de la muestra o múltiples mecanismos (chorro frente a disco), y reconocen que el requisito de potencias de chorro super-Eddington en algunas fuentes podría requerir modelos de múltiples zonas para una reconciliación completa.