Hybrid Simulations of Supersonic Shear Flows: II) Cosmic Ray Viscosity

Resumen Técnico: Simulaciones Híbridas de Flujos de Cizalla Supersónicos: II) Viscosidad de Rayos Cósmicos

1. Planteamiento del Problema
Los flujos de cizalla (shear flows) son omnipresentes en entornos astrofísicos, como en los límites de chorros (jets), discos de acreción y el medio interestelar (ISM). En plasmas colisionales, la viscosidad microscópica disipa la energía cinética libre. Sin embargo, en plasmas astrofísicos colisionales, la viscosidad atómica es despreciable, y el transporte de momento suele estar dominado por inestabilidades y turbulencia.

El problema central de este trabajo es comprender el papel dinámico de las partículas no térmicas (rayos cósmicos, CRs), ya sean preexistentes o generadas espontáneamente, en la inducción de una viscosidad cinética efectiva. A diferencia de los enfoques fluidos que solo pueden disipar energía cinética directamente en calor, los CRs, debido a sus grandes radios de giro, actúan como mensajeros a larga distancia que facilitan el intercambio de momento entre capas de flujo distantes. El objetivo es cuantificar cómo esta "viscosidad de CRs" afecta la disipación de la cizalla, la partición de la energía y la aceleración de partículas en flujos que van desde subsónicos hasta supersónicos.

2. Metodología
Los autores utilizan simulaciones híbridas de partículas en celda (PIC), donde los iones se tratan como partículas cinéticas y los electrones como un fluido neutro de masa nula.

• Configuración: Se simula un flujo de cizalla sinusoidal tipo Kolmogorov ($U_x(y) \propto \sin(2\pi y/L)$) en una caja 2D periódica.
• Parámetros: Se exploran diferentes densidades de CRs ($n_{CR}$), momentos isotrópicos ($p_{iso}$) y números de Mach de Alfvén ($M_A$) que varían de 0.5 a 40.
• Escalas: Las simulaciones cubren regímenes donde el radio de giro de los CRs ($r_g$) es menor o comparable a la escala de la cizalla ($L$), así como casos donde $r_g \gg L$.
• Comparativa: Se comparan casos con CRs preexistentes ("semillas") contra casos sin ellos (Run N5) y contra el caso de referencia (Run B, $M_A=20$).

3. Contribuciones Clave

• Viscosidad Cinética Inducida por CRs: Se demuestra cuantitativamente que los CRs actúan como un mecanismo de viscosidad cinética que acelera la disipación de la estructura de velocidad, incluso cuando no son energéticamente dominantes.
• Condición de Acoplamiento: Se establece que la eficiencia de esta viscosidad depende críticamente de que el radio de giro de los CRs sea menor o comparable a la escala de longitud de la cizalla ($r_g \lesssim L$). Si $r_g$ es demasiado grande, los CRs no se acoplan eficazmente con el plasma térmico ni con las fluctuaciones magnéticas generadas.
• Partición de Energía No Termodinámica: Se revela que la transferencia de energía entre el gas térmico y los CRs no es un simple proceso de equipartición termodinámica, sino un proceso cinético controlado por la capacidad de los CRs con radios de giro adecuados para extraer energía de la cizalla en descomposición.

4. Resultados Principales

• Disipación de la Cizalla:

  • La presencia de CRs reduce significativamente el tiempo de disipación de la cizalla (tiempo de reducción a la mitad, $\tau_{50}$).
  • Existe un umbral crítico (aprox. $n_{CR} \sim 1\%$ y densidad de energía $\eta_{CR} \lesssim 0.1$) por encima del cual la disipación se acelera drásticamente.
  • Aumentar la densidad de CRs siempre acelera la disipación, pero aumentar su energía (momento) solo es efectivo si $r_g < L$. Si $r_g$ excede la escala de cizalla, el efecto de viscosidad satura o disminuye.

• Aceleración de Partículas:

  • Efecto de Semilla: Los CRs preexistentes favorecen una aceleración mucho más rápida de los iones térmicos, generando colas no térmicas más pronunciadas y alcanzando energías máximas hasta 4 veces mayores que en ausencia de CRs.
  • Compensación (Trade-off): Existe una competencia entre la aceleración y la vida útil de la cizalla. Una viscosidad de CRs muy alta disipa la cizalla tan rápidamente que reduce el tiempo disponible para la aceleración continua, lo que puede suprimir las colas no térmicas en casos de densidad de CRs muy alta.
  • Límite de Hillas: La energía máxima alcanzable ($\bar{E}_{max}$) tiende a acercarse al límite de Hillas ($E_H \approx 40 E_0$), pero a menudo no lo alcanza debido a la rápida disipación de la fuente de energía (la cizalla).

• Partición de Energía:

  • La energía cinética libre inicial se redistribuye entre el gas térmico, los CRs y los campos magnéticos.
  • En flujos subsónicos o de bajo Mach, la energía puede fluir desde los CRs hacia el gas térmico (calentamiento).
  • En flujos supersónicos ($M_A \gtrsim 20$), la energía fluye predominantemente desde la cizalla hacia la aceleración de partículas y la amplificación del campo magnético.
  • La energía en fluctuaciones magnéticas generadas auto-consistentemente alcanza aproximadamente el 1% de la energía libre inicial, suficiente para acoplar dinámicamente a todas las especies.

• Dependencia del Número de Mach:

  • A medida que aumenta el número de Mach ($M_A$), la eficiencia de aceleración aumenta, pasando de un calentamiento térmico (bajo $M_A$) a la formación de espectros de potencia duros no térmicos (alto $M_A$).
  • La fracción de energía en partículas no térmicas puede alcanzar hasta el 40% del total para $M_A \gtrsim 20$.

5. Significado e Implicaciones
Este trabajo demuestra que la viscosidad inducida por partículas no térmicas es una propiedad intrínseca de los plasmas astrofísicos colisionales.

• Más allá de la MHD: Los enfoques de fluidos o MHD estándar no pueden capturar estos efectos cinéticos, ya que no distinguen entre calentamiento y aceleración, ni pueden modelar el intercambio de momento mediado por partículas con grandes radios de giro.
• Aplicaciones Astrofísicas: Los resultados sugieren que en entornos como choques de jets, discos de acreción y turbulencia intra-cúmulo, la presencia de CRs (ya sean acelerados in situ o preexistentes) es fundamental para determinar la tasa de disipación de la cizalla, la amplificación del campo magnético y la eficiencia de la aceleración de partículas a altas energías.
• Futuro: El estudio subraya la necesidad de considerar estos efectos cinéticos al extrapolar simulaciones a escalas astrofísicas reales, donde la inyección continua de energía podría permitir que las colas de potencia alcancen el límite de Hillas y que la turbulencia magnética crezca hacia la equipartición.

En resumen, el papel de los rayos cósmicos no es solo ser un producto de la aceleración, sino un agente activo que modifica la dinámica macroscópica del plasma a través de una viscosidad cinética efectiva.