The Effects of Cosmic Ray Protons on Galactic Nonthermal Filaments

Resumen Técnico: Los Efectos de los Protones de Rayos Cósmicos en los Filamentos No Térmicos Galácticos

Planteamiento del Problema
El Centro Galáctico (CG) alberga una población de filamentos no térmicos (NTF, por sus siglas en inglés) que se extienden decenas de pársecs de longitud, iluminados por radiación sincrotrón proveniente de rayos cósmicos (CR). El origen físico de estas estructuras sigue sin resolverse. Se proponen dos mecanismos de inyección primarios: (1) los NTF son alimentados por chorros de nebulosas de viento de púlsar, lo que resulta en una población dominada por leptones; o (2) los NTF son alimentados por partículas aceleradas por choques interestelares, lo que resulta en una población dominada por protones (con una relación protón-leptón de aproximadamente 100:1). Si bien la población de CR galáctica exhibe una distribución de energía de ley de potencia consistente con las observaciones de los NTF, la conexión causal específica entre la fuente de inyección y la morfología del filamento no está clara. Además, la topología del campo magnético poloidal dentro de la Zona Molecular Central (CMZ) es compleja, con muchos filamentos que parecen no lineales, lo que sugiere el papel potencial de la turbulencia interestelar. Este estudio tiene como objetivo distinguir entre los escenarios dominados por leptones y dominados por protones mediante la simulación de la propagación de CR, el calentamiento y el flujo de plasma bajo diversas condiciones del Centro Galáctico.

Metodología
Los autores emplearon el código de magnetohidrodinámica (MHD) Athena++, modificado para incluir términos de pérdida radiativa y colisional, para simular la propagación de CR en un tubo de flujo 1D que representa un filamento. La simulación modela un medio de fondo uniforme perturbado por la inyección de energía de CR en el punto medio del filamento.

• Procesos Físicos: El modelo incorpora difusión (dispersión por ondas de Alfvén) y streaming (convección a la velocidad de Alfvén). Los mecanismos de pérdida para los leptones incluyen radiación sincrotrón, dispersión de Compton inverso, pérdidas de Coulomb y bremsstrahlung. La producción secundaria se modela donde los protones primarios interactúan con protones ambientales para generar electrones y positrones secundarios.
• Parámetros: El estudio varía parámetros clave incluyendo la fuerza del campo magnético ($B$), la densidad del plasma ($n_p$), el coeficiente de difusión de CR ($\kappa$) y la energía total inyectada ($E_{\pm}$).
• Conjuntos de Simulación: Cada conjunto de parámetros fue simulado dos veces: una vez solo con leptones y otra vez con una mezcla de protón-leptón (escalando la tasa de inyección total por 100 para contabilizar la relación protón-leptón). Los parámetros fiduciarios adoptados fueron $B = 200\,\mu\text{G}$, $n_p = 100\,\text{cm}^{-3}$ y $\kappa = 1 \times 10^{26}\,\text{cm}^2\,\text{s}^{-1}$.
• Análisis: Los autores analizaron los perfiles espaciales de densidad de energía, los flujos de gas inducidos, las tasas de calentamiento y el índice espectral y la luminosidad sincrotrón resultantes a lo largo del tiempo (hasta 40 kyr).

Resultos Clave

1. Indistinguibilidad de los Mecanismos de Inyección: Las simulaciones revelaron pocas diferencias observables entre los casos dominados por protones y los dominados por leptones. Aunque la inyección de protones impulsa flujos de gas más fuertes y proporciona más calentamiento debido a la mayor densidad de energía total, estos efectos permanecen por debajo de los umbrales detectables. Incluso en las densidades más bajas consideradas, las velocidades inducidas ($\sim 2.5\,\text{km s}^{-1}$) y las tasas de calentamiento son insignificantes en comparación con el enfriamiento radiativo y los anchos de línea térmicas típicas.
2. Restricciones de Propagación: El brillo superficial observado, que presenta una meseta uniforme, sugiere un coeficiente de difusión significativamente menor que el valor galáctico típico ($\kappa \sim 10^{28}\,\text{cm}^2\,\text{s}^{-1}$). Los autores encuentran que es necesario un coeficiente de difusión más bajo ($\kappa \sim 10^{26}\,\text{cm}^2\,\text{s}^{-1}$) para reproducir los perfiles de luminosidad observados.
3. Longitud y Edad del Filamento: Las simulaciones indican que, sin un mecanismo para truncar la propagación, los filamentos continuarían expandiéndose indefinidamente. La longitud de un NTF es, por lo tanto, dependiente de su edad. Las pérdidas radiativas y colisionales, incluso combinadas con la producción secundaria, no crean una longitud de estado estacionario porque los leptones secundarios nacen más lejos del sitio de inyección y se difunden de manera similar a las partículas primarias.
4. Variabilidad del Índice Espectral: El estudio exploró mecanismos para la variación observada en los índices espectrales de sincrotrón (que oscilan entre -2 y +0.5).
   • Difusión dependiente de la energía: Si la difusividad escala con la energía ($\kappa \propto \gamma^q$), el índice espectral varía a lo largo del filamento, endureciéndose generalmente al alejarse del sitio de inyección.
   • Envejecimiento de sincrotrón: Las partículas de alta energía pierden energía más rápido, lo que conduce a un endurecimiento espectral con el tiempo.
   • Producción secundaria: Endurece el espectro en tiempos tardíos.
   • A pesar de estos mecanismos, los autores señalan que ninguna combinación única de estos efectos explica completamente las variaciones de índice espectral irregulares y asimétricas observadas en los NTF reales.

Significado y Conclusiones
La contribución principal de este trabajo es la demostración de que la presencia de una gran población de protones de rayos cósmicos no produce firmas distintas y observables en la morfología, el flujo o el calentamiento de los filamentos que permitan a los astrónomos distinguir entre los orígenes impulsados por púlsares (leptones) o por choques (protones) basándose en las simulaciones actuales de MHD 1D.

Consecuentemente, los autores argumentan que la incapacidad de distinguir entre estos dos modelos de inyección estándar, combinada con las imperfecciones estructurales observadas de los NTF (quiebres, divisiones y alineaciones no lineales), motiva la consideración de un tercer mecanismo de formación. El artículo propone que los NTF pueden surgir de estructuras intermitentes en la turbulencia del Centro Galáctico. En este escenario, la turbulencia MHD inhomogénea podría amplificar los campos magnéticos y estirar los remolinos en filamentos delgados, creando las estructuras observadas sin requerir una fuente de inyección de CR específica y localizada. Los autores concluyen que el trabajo futuro debería centrarse en simular estos efectos turbulentos para comprender mejor el origen de los NTF.