Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment

Título: Implementación del algoritmo CRESP en el código MHD PIERNIK. II. Propagación de núcleos primarios y secundarios en un entorno magnetohidrodinámico.

1. El Problema

La comprensión de la propagación de los Rayos Cósmicos (RC) y su interacción con el Medio Interestelar (MIE) es fundamental para la astrofísica de partículas. Los modelos tradicionales (como GALPROP) suelen tratar el MIE como un entorno estacionario, utilizando ecuaciones de tipo "caja con fugas" o ecuaciones de Fokker-Planck en un campo magnético fijo. Sin embargo, estos enfoques a menudo ignoran la retroalimentación dinámica entre los RC, el gas térmico y el campo magnético galáctico.

Existe una necesidad crítica de modelos autoconsistentes que acoplen la propagación de múltiples especies de RC (núcleos primarios y secundarios) con la evolución dinámica del MIE. Específicamente, se requiere validar algoritmos numéricos capaces de simular la producción de núcleos secundarios (como Litio, Berilio y Boro) a través de reacciones de espalación, y cómo estos procesos afectan las relaciones observables (como B/C o $^{10}\text{Be}/^{9}\text{Be}$) en un entorno donde los RC impulsan vientos galácticos y modifican la estructura del campo magnético.

2. Metodología

Los autores extienden el algoritmo CRESP (Cosmic Ray Energy SPectrum) dentro del código de Magnetohidrodinámica (MHD) PIERNIK.

• Enfoque Numérico: Se utiliza un método de dos momentos (densidad numérica y densidad de energía cinética) con una aproximación de ley de potencia por tramos (piece-wise power-law) para resolver la ecuación de Fokker-Planck en un espacio de momentos discretizado (16 bins).
• Especies Simuladas:
  • Primarios: Se aceleran en remanentes de supernova (SN) y se inyectan como núcleos de Carbono ($^{12}\text{C}$), Nitrógeno ($^{14}\text{N}$) y Oxígeno ($^{16}\text{O}$). Los protones se tratan en una aproximación "gris" (no resuelta espectralmente) para impulsar la dinámica del gas.
  • Secundarios: Se generan mediante colisiones de espalación de los primarios contra protones térmicos del MIE, produciendo isótopos de Litio ($^{7}\text{Li}$), Berilio ($^{9}\text{Be}$, $^{10}\text{Be}$) y Boro ($^{10}\text{B}$, $^{11}\text{B}$).
• Física Incluida:
  • Ecuaciones MHD resistivas acopladas dinámicamente a los RC.
  • Difusión anisotrópica dependiente de la rigidez ($D \propto R^\delta$).
  • Advección con el gas térmico.
  • Pérdidas de energía por expansión adiabática y enfriamiento Coulombiano.
  • Decaimiento radiactivo (para $^{10}\text{Be}$).
• Configuración de la Simulación:
  • Se utiliza una caja 3D estratificada por gravedad que representa un parche local del MIE cerca del Sistema Solar (0.5 x 1 x 8 kpc).
  • Se inyectan supernovas aleatoriamente que proporcionan el 10% de su energía en protones RC.
  • Se realiza un estudio de parámetros variando la tasa de supernovas, el coeficiente de difusión paralelo ($D_{\parallel}$) y el índice de potencia de la difusión ($\delta$).

3. Contribuciones Clave

• Extensión del Algoritmo CRESP: Se ha implementado exitosamente la capacidad de manejar múltiples especies nucleares hadrónicas (primarias y secundarias) con resolución espectral dentro de un entorno MHD vivo, superando las limitaciones de modelos anteriores que solo trataban electrones o protones en aproximación gris.
• Acoplamiento Dinámico Completo: El modelo demuestra cómo la presión de los RC y sus efectos de flotabilidad (buoyancy) modifican la estructura del campo magnético y generan vientos galácticos, los cuales a su vez retroalimentan la propagación de los RC.
• Validación de Canales de Espalación: Se ha verificado la producción correcta de núcleos secundarios y la evolución de sus espectros en comparación con los primarios, incluyendo el tratamiento de isótopos inestables como $^{10}\text{Be}$.

4. Resultados Principales

El estudio de parámetros reveló comportamientos complejos y no intuitivos:

• Impacto de la Tasa de Supernovas (SN):

  • Una mayor tasa de SN genera vientos más fuertes y una mayor advección vertical.
  • Esto reduce el tiempo de residencia de los RC primarios en el disco denso, disminuyendo la producción de secundarios.
  • Resultado: Una mayor tasa de SN conduce a una menor relación B/C (Boro/Carbono).

• Impacto del Coeficiente de Difusión ($D_{\parallel}$) - Hallazgo Novel:

  • Contrario a los modelos de "caja con fugas" o estacionarios (donde un mayor $D$ implica una menor relación B/C debido a una fuga más rápida), este modelo autoconsistente predice lo contrario.
  • Un mayor coeficiente de difusión distribuye a los RC de manera más uniforme a lo largo de las líneas de campo horizontales, reduciendo la inyección local de presión que causa la flotabilidad.
  • Esto resulta en un campo magnético vertical más débil (menor inclinación).
  • Un campo vertical más débil atrapa a los primarios por más tiempo cerca del plano del disco, aumentando la producción de secundarios.
  • Resultado: En este modelo, B/C aumenta al aumentar el coeficiente de difusión.

• Impacto del Índice de Potencia ($\delta$):

  • La relación B/C es sensible a $\delta$. Un valor de $\delta = 0.5$ (modelo C1) mostró un mejor ajuste cualitativo con los datos de AMS-02 en energías relativistas (>10 GeV/n) en comparación con $\delta = 0.3$, aunque la amplitud absoluta requiere un ajuste fino de otros parámetros.

• Espectros y Ratios:

  • Se observó que las pendientes espectrales de los núcleos primarios y secundarios difieren debido a la mezcla de partículas viejas (enfriadas) y nuevas (inyectadas), así como a las pérdidas por Coulomb en energías no relativistas.
  • La relación $^{10}\text{Be}/^{9}\text{Be}$ se ajusta bien a los datos en el disco, pero muestra desviaciones significativas a grandes altitudes debido a la desintegración radiactiva y la menor densidad de primarios.

5. Significado e Implicaciones

• Cambio de Paradigma en la Interpretación de Datos: Este trabajo demuestra que los observables de RC (como B/C) no dependen únicamente de los parámetros de difusión de entrada, sino que son el resultado de una retroalimentación dinámica entre los RC, el gas y el campo magnético. Ignorar este acoplamiento puede llevar a interpretaciones erróneas de los parámetros de transporte.
• Validación de Modelos Autoconsistentes: El éxito de PIERNIK+CRESP en reproducir tendencias observacionales (aunque con simplificaciones) valida el uso de simulaciones MHD completas para estudiar la evolución galáctica y la física de los RC.
• Limitaciones y Futuro: Los autores reconocen que la resolución actual (62.5 pc) y la falta de un MIE multifase o de efectos de "streaming" (flujo de RC) son limitaciones. Sin embargo, el estudio establece una base sólida para futuros modelos que incluyan turbulencia a pequeña escala y un MIE más realista, mejorando la precisión en la determinación de los tiempos de escape de los RC de la Galaxia.

En resumen, el paper establece que la estructura del campo magnético inducida por los propios RC es un factor crítico que modifica la relación entre la difusión y la abundancia de núcleos secundarios, desafiando las intuiciones derivadas de modelos estacionarios simplificados.