Cosmic Ray Magnetohydrodynamics: A New Two-Moment Framework with Numerical Implementation

Resumen Técnico: Magnetohidrodinámica de Rayos Cósmicos: Un Nuevo Marco de Dos Momentos con Implementación Numérica

Planteamiento del Problema
Los rayos cósmicos (RC) son componentes críticos en sistemas astrofísicos, proporcionando retroalimentación a través de una amplia gama de escalas, desde la regulación de la formación estelar hasta el impulso de vientos galácticos. Sin embargo, modelar el transporte de los RC sigue siendo un desafío debido a la naturaleza multiescala del problema y a la complejidad de su microfísica. Si bien los enfoques recientes de hidrodinámica de RC de dos momentos han mejorado respecto a los métodos anteriores de un solo momento (que sufren de inestabilidades numéricas en los extremos de presión), los marcos existentes poseen limitaciones. Específicamente, las formulaciones actuales a menudo asumen que las ondas de Alfvén se propagan en una sola dirección a lo largo del campo magnético, fallando al capturar configuraciones de ondas más generales, como las impulsadas por la Inestabilidad de la Anisotropía de Presión de los Rayos Cósmicos (CRPAI). Además, persisten desafíos numéricos relacionados con las inestabilidades derivadas de los términos de fuente de flujo (streaming) y las ambigüedades en el diseño de los solvers para sistemas acoplados de CR-MHD.

Metodología
Los autores desarrollan una nueva descripción de dos momentos para la hidrodinámica de los RC derivada de los primeros principios, partiendo de la ecuación de Boltzmann no colisional para partículas relativistas. Los pasos metodológicos clave incluyen:

1. Derivación de la Teoría Cinética: El marco comienza con la función de distribución de los RC en un marco que se mueve con la velocidad perpendicular del gas ($\vec{u}_\perp$). Al tomar los momentos de la ecuación cinética, los autores derivan las ecuaciones de evolución para la densidad de energía de los RC ($E_{cr}$) y el flujo de energía ($\vec{F}_{cr}$).
2. Incorporación de Anisotropía y Ondas Bidireccionales: A diferencia de modelos anteriores que asumen implícitamente la propagación unidireccional de ondas, esta formulación contabiliza explícitamente la anisotropía de presión de los RC ($\Delta P_{cr}$) y permite ondas de Alfvén que se propagan en ambas direcciones (hacia adelante y hacia atrás) a lo largo del campo magnético. Esto se logra expandiendo la función de distribución utilizando polinomios de Legendre hasta el segundo orden.
3. Prescripciones de Física de Subrejilla: El marco incorpora coeficientes de dispersión ($\sigma$) derivados de procesos microfísicos. Los autores proporcionan prescripciones para:
   • Inestabilidad de Flujo de Rayos Cósmicos (CRSI): Donde los gradientes de presión impulsan ondas en una dirección específica.
   • Inestabilidad de la Anisotropía de Presión de los Rayos Cósmicos (CRPAI): Donde la invariancia adiabática en campos magnéticos variables impulsa ondas en ambas direcciones.
   • Turbulencia Extrínseca: Donde las ondas no tienen una dirección preferida.
4. Implementación Numérica en Athena++: Las ecuaciones se implementan como una especie fluida acoplada al gas MHD en el código Athena++.
   • Solver de Riemann: Se emplea un nuevo solver de Riemann HLLE. Para mitigar la difusión numérica excesiva, los autores derivan la relación de dispersión completa del sistema CR-MHD para identificar velocidades de señal realistas (incluyendo ondas acústicas de CR) en lugar de depender únicamente del límite de la velocidad reducida de la luz.
   • Integración Temporal: Se utilizan integradores implícitos (Backward Euler e Implicit van Leer 2) para manejar términos de fuente rígidos (stiff) provenientes de la dispersión.
   • Criterio de Estabilidad: Se realiza un análisis de estabilidad de von Neumann para derivar una guía de paso de tiempo diseñada específicamente para suprimir oscilaciones a escala de rejilla causadas por términos de flujo con cambio de signo.

Resultados Clave
Los autores validan el nuevo módulo a través de una serie de pruebas de referencia:

• Pruebas de Flujo (Streaming) 1D: El código reproduce con éxito las soluciones analíticas para el flujo de los RC bajo regímenes de amortiguamiento tanto lineales (ión-neutro) como no lineales (Landau). Las pruebas confirman que el criterio de estabilidad derivado mitiga eficazmente las inestabilidades numéricas cuando los coeficientes de dispersión son grandes.
• Pruebas de Difusión 1D: El módulo modela con precisión el transporte difusivo impulsado por la turbulencia extrínseca, coincidiendo con las soluciones de la función de Green para perfiles gaussianos tanto en fondos de gas estáticos como en movimiento.
• Transporte Anisotrópico: En pruebas 2D con campos magnéticos circulares, el código maneja correctamente la alineación del flujo de los RC con las líneas de campo, demostrando que el término de la fuerza de Lorentz contrarresta con éxito la inercia perpendicular de los RC.
• Coordenadas Curvilíneas: El marco funciona correctamente en coordenadas esféricas y cilíndricas, manejando discontinuidades agudas y efectos geométricos (por ejemplo, campos monopolares y anisotropía de presión cilíndrica) consistentes con las expectativas teóricas.
• Relaciones de Dispersión: La derivación de la relación de dispersión completa identifica la inestabilidad acústica de los CR y otras ramas de ondas, proporcionando referencias rigurosas para el código y permitiendo el uso de velocidades de señal físicamente realistas en el solver de Riemann.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco general y extensible para el modelado de la retroalimentación de los RC que supera las limitaciones direccionales de los esquemas previos de dos momentos. Al acomodar configuraciones arbitrarias de ondas de Alfvén e incluir explícitamente la anisotropía de presión de los RC, el marco puede incorporar diversos mecanismos de transporte, incluyendo CRSI, CRPAI y sus combinaciones. Los autores enfatizan que la inclusión de velocidades de señal realistas en el solver de Riemann y la derivación de un criterio de paso de tiempo específico mejoran significativamente la robustez numérica. Este trabajo tiene como objetivo facilitar estudios más precisos de fenómenos impulsados por los RC, tales como vientos galácticos, interacciones de nubes y la estabilidad térmica de los medios circumgalácticos, al proporcionar una herramienta numéricamente estable que une los procesos de dispersión microfísica con la dinámica de fluidos macroscópica.