Counter-streaming heat-flux closure for electron-only collisionless magnetic reconnection

Resumen Técnico: Cierre de Flujo de Calor por Contracorriente para la Reconexión Magnética de Solo Electrones sin Colisiones

Planteamiento del Problema
La reconexión magnética (RM) en plasmas espaciales turbulentos exhibe cada vez más regímenes de "solo electrones", donde las hojas de corriente que se reconectan son más pequeñas que las escalas características de los iones, desacoplando a estos últimos de la dinámica. Aunque las simulaciones cinéticas de Particle-in-Cell (PIC) han revelado que la RM de solo electrones presenta fenómenos distintivos —como chorros super-iónicos-Alfvénicos y tasas de reconexión más rápidas—, los modelos de fluidos convencionales no logran capturar estas dinámicas. Los cierres de fluidos existentes, como el cierre de Hammett-Perkins (que depende de distribuciones casi Maxwellianas) o el cierre de presión de Le et al. (derivado para regímenes de alto $\beta$), son fundamentalmente incapaces de reproducir las distribuciones de electrones fuertemente no-Maxwellianas y multimodales observadas en la RM de solo electrones. Específicamente, los modelos actuales no pueden dar cuenta del "contracorriente de electrones" (counter-streaming), donde subpoblaciones de electrones coherentes con velocidades de fluido distintas ocupan la misma ubicación espacial. Esta limitación impide el modelado eficiente de la RM de solo electrones, lo cual es crítico para comprender entornos como la cola magnética de la Tierra.

Metodología
Los autores utilizan simulaciones cinéticas PIC completas (usando el código OSIRIS) para investigar el origen físico del flujo de calor de los electrones en la RM de solo electrones sin colisiones. Inicializan una hoja de corriente Harris doble de fuerza constante en 2D-3V con parámetros representativos de la cola magnética ($\beta_{Reconn.} = 0.1$, $B_{Guide}/B_{Reconn.} = 0.1$). Al etiquetar los electrones basándose en sus posiciones iniciales por encima o por debajo del plano medio de la hoja de corriente, los autores descomponen la función de distribución total de los electrones $f$ en dos subpoblaciones ($f_1$ y $f_2$).

A través de esta descomposición, los autores derivan expresiones analíticas para el tensor de presión térmica total y el tensor de flujo de calor. Demuestran que el flujo de calor total está dominado por un término de "contracorriente" ($q_{CS}$) que surge de la deriva relativa y las diferencias de temperatura entre las dos subpoblaciones, en lugar de provenir de un comportamiento no adiabático dentro de las propias subpoblaciones.

Basándose en este hallazgo, los autores desarrollan un nuevo cierre de dos fluidos de electrones (2-e-Fluid). En este modelo:

1. La población de electrones se divide en dos fluidos distintos inicializados en lados opuestos de la hoja de corriente.
2. Cada subpoblación se trata como aproximadamente adiabática ($q_{1,ijk} \approx q_{2,ijk} \approx 0$).
3. El flujo de calor total se cierra mediante el movimiento relativo de los dos fluidos ($q_{ijk} \approx q_{CS,ijk}$), eliminando la necesidad de una ecuación de tercer momento.

Los autores implementan este cierre en el código de fluidos Gkeyll, evolucionando dos fluidos de electrones con 10 momentos cada uno (densidad, velocidad y componentes del tensor de presión). Para comparar, también ejecutan una simulación de un solo fluido de electrones (1-e-Fluid) bajo las mismas condiciones.

Resultados Clave

• Origen del Flujo de Calor: El análisis PIC confirma que el flujo de calor de los electrones en la RM de solo electrones surge casi enteramente de la contracorriente entre subpoblaciones provenientes de lados opuestos de la hoja de corriente. Aunque la distribución total es altamente no-Maxwelliana, cada subpoblación mantiene un carácter aproximadamente adiabático a lo largo de la región de reconexión.
• Desempeño del Modelo de Fluidos: El modelo 2-e-Fluid reproduce con éxito características cinéticas clave observadas en las simulaciones PIC, incluyendo:
  • Densidad de Corriente de Reconexión ($J_z$): A diferencia del modelo 1-e-Fluid, que sufre de un agotamiento de densidad no físico en el centro de la hoja de corriente (lo que conduce a una $J_z$ suprimida), el modelo 2-e-Fluid mantiene una hoja de corriente llena al permitir que las subpoblaciones crucen el plano medio.
  • Presión Térmica ($p_{xx}$): El modelo 1-e-Fluid desarrolla acumulaciones de presión espurias en los flujos de salida debido a la falta de canales de flujo de calor. El modelo 2-e-Fluid captura el flujo de calor por contracorriente, el cual transporta la energía térmica lejos de los flujos de salida, evitando la acumulación de presión.
  • Estructura del Flujo de Calor: El modelo 2-e-Fluid reproduce la estructura espacial y el signo del flujo de calor por contracorriente ($q_{CS,xxx}$), el cual constituye la gran mayoría del flujo de calor total.
• Regimen de Validez: Un escaneo de parámetros sobre $\beta_{Reconn.}$ y $B_{Guide}/B_{Reconn.}$ revela que el cierre es más preciso en valores bajos de estos parámetros (relevantes para la cola magnética de la Tierra). La precisión está controlada por la relación entre la velocidad térmica inicial de los electrones y la velocidad de flujo inicial ($v_{th}/v_{fl}$). El modelo se degrada a medida que esta relación aumenta o cuando las subpoblaciones comienzan a circular dentro de las islas magnéticas, rompiendo la simple suposición de contracorriente.

Significancia y Reivindicaciones
El artículo afirma resolver una limitación fundamental en el modelado de la RM de solo electrones sin colisiones al identificar la contracorriente como el origen físico del flujo de calor de los electrones. Al demostrar que las distribuciones complejas y no-Maxwellianas pueden descomponerse en dos subpoblaciones adiabáticas, los autores establecen un cierre de fluidos computacionalmente eficiente que captura la física que antes solo era accesible mediante simulaciones cinéticas.

Los autores sostienen que este enfoque 2-e-Fluid permite el modelado de la RM de solo electrones con una descripción reducida que los modelos de un solo fluido convencionales "no pueden reproducir fundamentalmente". Destacan que este método captura la importancia dinámica del flujo de calor para modificar la presión térmica en escalas de tiempo electrónicas. Aunque el artículo sugiere posibles extensiones a la RM de solo electrones en 3D y a la reconexión turbulenta, señala modestamente que el cierre actual es más robusto en regímenes de bajo $\beta$ y bajo campo guía, y que podría requerir fluidos adicionales para resolver la dinámica compleja dentro de las islas. El trabajo proporciona una vía para el modelado eficiente de la reconexión de solo electrones sin sacrificar la física esencial de la contracorriente.