The effect of inverse Compton losses on particle acceleration in three-dimensional relativistic reconnection

Resumen Técnico: El Efecto de las Pérdidas por Compton Inverso en la Aceleración de Partículas en la Reconexión Relativista Tridimensional

Planteamiento del Problema
La reconexión magnética relativista es un mecanismo primario para disipar la energía magnética y acelerar partículas en entornos astrofísicos, tales como magnetosferas de objetos compactos y chorros relativistas. Simulaciones recientes de primeros principios de tipo partícula-en-celda (PIC) han establecido un modelo para la aceleración de partículas en la reconexión 3D (Zhang et al. 2023), distinguiendo entre una "fase libre" donde las partículas se aceleran activamente en la región aguas arriba (upstream) y una "fase atrapada" donde las partículas están confinadas en cuerdas de flujo aguas abajo (downstream) sin una energización significativa. Este modelo predice un espectro de ley de potencia universal ($dN/d\gamma \propto \gamma^{-1}$) para partículas libres y un espectro más pronunciado ($dN/d\gamma \propto \gamma^{-2}$) para partículas atrapadas en ausencia de enfriamiento radiativo.

Sin embargo, el comportamiento de este mecanismo de aceleración en presencia de enfriamiento radiativo —específicamente las pérdidas por Compton Inverso (IC)— permanece poco explorado. No está claro si las pérdidas por enfriamiento alteran la tasa de reconexión, la estructura de la capa o la física fundamental de las fases libre y atrapada, particularmente en el régimen de "enfriamiento débil" donde el factor de Lorentz limitado por la reacción de radiación ($\gamma_{\text{rad}}$) excede la magnetización ($\sigma$).

Metodología
Los autores realizaron una serie de simulaciones PIC 3D utilizando el código TRISTAN-MP para investigar la reconexión relativista con pérdidas por enfriamiento IC. Las simulaciones utilizaron una configuración de lámina de corriente de Harris con un campo guía ($B_g = 0.1 B_0$) y condiciones de contorno periódicas a lo largo de la dirección de la corriente.

Los parámetros clave de la simulación incluyeron:

• Composición del Plasma: Pares electrón-positrón, con enfriamiento radiativo aplicado únicamente a los electrones. Los positrones se dejaron sin enfriar para imitar el comportamiento de los iones pesados en plasmas electrón-ión y para asegurar una alteración mínima del equilibrio de fuerzas de la lámina, permitiendo una comparación directa con el modelo sin enfriamiento de Zhang et al. (2023).
• Magnetización: Una magnetización fija de $\sigma = 10$.
• Régimen de Enfriamiento: El estudio se centró en el régimen de enfriamiento débil ($\gamma_{\text{rad}} > \sigma$). Se probaron tres valores específicos para el factor de Lorentz limitado por la reacción de radiación: $\gamma_{\text{rad}} = 30, 43$ y $60$. Estos corresponden a rupturas de enfriamiento ($\gamma_{\text{cool}}$) de aproximadamente $3.5, 7.2$y$14.2$, respectivamente.
• Análisis: Los autores rastrearon grandes muestras de partículas (0.5 millones de electrones por ejecución) para distinguir entre poblaciones "libres" (aguas arriba) y "atrapadas" (aguas abajo) utilizando un criterio de factor de mezcla. Analizaron las tasas de aceleración, los tiempos de escape y los espectros de energía, comparándolos con las predicciones analíticas derivadas del modelo de Zhang et al. (2023) extendido para incluir el enfriamiento.

Contribuciones Clave y Resultados

1. Robustez de la Dinámica de la Reconexión: Las simulaciones demuestran que las pérdidas por enfriamiento de los electrones no alteran apreciablemente la tasa de reconexión ni la estructura macroscópica de la lámina de corriente. La tasa de reconexión en estado estacionario se establece en $\eta_{\text{rec}} \approx 0.05$, consistente con los casos sin enfriamiento. Aunque un enfriamiento más fuerte conduce a láminas ligeramente más delgadas y plasmoides más pequeños, la física fundamental de la lámina permanece inalterada a través de las intensidades de enfriamiento probadas.

2. Validación de la Física de la Fase Libre: El estudio confirma que la "fase libre" de aceleración es insensible a las pérdidas por enfriamiento. Los electrones libres se aceleran consistentemente a la tasa óptima ($\dot{\gamma}_{\text{acc}} \sim \eta_{\text{rec}}\omega_c$) y exhiben un espectro de ley de potencia de $dN_{\text{free}}/d\gamma \propto \gamma^{-1}$, idéntico al caso sin enfriamiento. El corte de este espectro está determinado por $\gamma_{\text{rad}}$ en lugar del tamaño del sistema.

3. Modificación del Espectro de la Fase Atrapada: Un hallazgo crítico es la modificación del espectro de las partículas atrapadas en presencia de enfriamiento.

   • En el caso sin enfriamiento, el espectro atrapado escala como $dN/d\gamma \propto \gamma^{-2}$.
   • En el régimen de enfriamiento débil, para energías por encima de la ruptura de enfriamiento ($\gamma_{\text{cool}} < \gamma < \gamma_{\text{rad}}$), el espectro de los electrones atrapados se vuelve más pronunciado hacia $dN/d\gamma \propto \gamma^{-3}$.
   • Este aumento de la pendiente confirma el supuesto fundamental del modelo de Zhang et al. (2023): que no ocurre una energización significativa durante la fase atrapada. El aumento de la pendiente surge porque las partículas atrapadas se enfrían en una escala de tiempo comparable o menor que el tiempo de advección, reduciendo su población en energías más altas en relación con el caso sin enfriamiento.

4. Cinemática de Partículas: Las simulaciones verifican que las partículas libres se mueven casi a lo largo del campo eléctrico de reconexión ($\beta_z \approx -1$), facilitando una aceleración rápida, mientras que las partículas atrapadas exhiben una distribución simétrica de $\beta_z$ alrededor de cero, consistente con la falta de una energización sistemática.

Significancia e Implicaciones Astrofísicas
El artículo valida el modelo analítico propuesto por Zhang et al. (2023) para la reconexión relativista 3D en presencia de un enfriamiento radiativo débil. La implicación principal es que los modelos de emisión radiativa para fuentes impulsadas por reconexión (como los Brotes de Rayos Gamma y los chorros de Blazares) deben emplear una estructura de dos zonas:

• Zona 1 (Libre): Partículas acelerándose rápidamente con un espectro $dN/d\gamma \propto \gamma^{-1}$.
• Zona 2 (Atrapada): Partículas enfriándose pasivamente con un espectro que se vuelve más pronunciado de $\gamma^{-2}$ a $\gamma^{-3}$ por encima de la ruptura de enfriamiento.

Los autores señalan que, si bien la mayoría de las partículas emisoras están atrapadas, las partículas libres se vuelven más abundantes cerca del factor de Lorentz máximo ($\gamma_{\text{rad}}$). Debido a que las partículas libres se proyectan preferencialmente a lo largo del campo eléctrico de reconexión, su emisión podría resultar en firmas observacionales únicas, como la variabilidad temporal rápida.

Limitaciones y Advertencias
Los autores declaran explícitamente que sus resultados se aplican al régimen de enfriamiento débil donde $\gamma_{\text{rad}} > \sigma$. El modelo no es aplicable a escenarios de enfriamiento fuerte donde $\gamma_{\text{rad}} < \sigma$. Además, el rango dinámico entre $\gamma_{\text{cool}}$ y $\gamma_{\text{rad}}$ en las simulaciones es menor que en sistemas astrofísicos realistas. Se sugiere trabajo futuro para extender estos hallazgos a plasmas electrón-protón con relaciones de masa realistas y a la reconexión con campos guía significativos.