Electron-positron Pair Production in Global GRMHD Simulations of Black Hole Accretion Flows

Resumen Técnico: Producción de Pares en Simulaciones GRMHD de Agujeros Negros

1. Planteamiento del Problema

Existe un interés creciente en determinar si los pares electrón-positrón juegan un papel crucial en la dinámica y termodinámica de los agujeros negros acrecientes. En el límite ópticamente delgado, se cree que estos pares son esenciales para apantallar el campo eléctrico y conducir corrientes en los chorros (jets) de Blandford-Znajek. Sin embargo, los mecanismos para suministrar este plasma de pares a los chorros han sido un tema de debate (modelos de aceleración en brechas vs. modelos de "drizzle" o goteo).

A pesar de que la física de pares ha sido estudiada en modelos analíticos y simulaciones de partículas (PIC), su impacto en la dinámica global de los flujos de acreción en simulaciones de Magnetohidrodinámica Relativista General (GRMHD) ha permanecido poco explorado. Las preguntas clave son:

• ¿Cuál es la escala de tiempo para alcanzar el equilibrio de pares?
• ¿Cómo se distribuyen espacial y temporalmente los pares?
• ¿Respetan los pares el modelo de equilibrio de una sola zona (one-zone)?
• Si no lo hacen, qué mecanismos (como la advección) impiden el equilibrio local?

2. Metodología

Los autores realizaron simulaciones globales tridimensionales de GRMHD utilizando el código de código abierto iharm3d.

• Enfoque de Física de Pares:

  • Se modeló la densidad de masa de positrones ($\rho_+$) como un escalar pasivo que se advecta con el campo de velocidad común del fluido (protones, electrones y positrones).
  • Se asumió una temperatura constante para los electrones y positrones ($T_e$), variando entre $10^9$ y $10^{10.5}$ K.
  • La producción y aniquilación de pares se resolvieron mediante un esquema de división de operadores (operator split), calculando las tasas de creación ($\dot{n}_C$) y aniquilación ($\dot{n}_A$) basadas en procesos de colisión partícula-partícula, fotón-partícula y fotón-fotón.
  • Se impuso un límite crítico de profundidad óptica de dispersión de protones ($\tau_{crit}$). Cuando $\tau_p > \tau_{crit}$, la producción de pares se desactiva manualmente (permitiendo solo aniquilación) para evitar soluciones inestables y no físicas (la rama dominada por pares), justificando esto mediante el enfriamiento eficiente por colisiones de Coulomb en discos ópticamente gruesos.

• Configuración de las Simulaciones:

  • Estado: Acreción en estado SANE (Stable and Normal Evolution).
  • Geometría: Se probaron diferentes relaciones de aspecto del disco ($H/r = 0.1, 0.2, 0.3$) para modelar diferentes tasas de acreción.
  • Parámetros: Se varió la masa del agujero negro ($10^8 M_\odot$ y $10 M_\odot$) y la tasa de acreción ($\dot{M}$) hasta cerca del límite de Eddington.
  • Enfriamiento: Se incluyó una función de enfriamiento radiativo para mantener la altura del disco y simular diferentes regímenes de acreción.

3. Contribuciones Clave

• Implementación de Escalar Pasivo: A diferencia de modelos anteriores de fluido único, este estudio rastrea explícitamente la densidad de positrones, evitando la inyección artificial de masa debida a los "floors" (suelos) de densidad en regiones magnetizadas.
• Análisis de Escalas de Tiempo: Se cuantificó la competencia entre la escala de tiempo de equilibrio de pares y las escalas de tiempo dinámicas (advección, orbital, de infalo).
• Mapeo de Fuentes de Pares: Se identificaron las regiones dominantes de producción de pares (fotón-fotón vs. partícula-partícula) y cómo la advección transporta estos pares desde el disco hacia la corona y los chorros.

4. Resultados Principales

• Distribución Espacial y "Vacío de Pares" (Pair Void):

  • En modelos de alta tasa de acreción ($H/r = 0.1$), se observa una fracción de pares máxima de $\sim O(0.01)$, consistente con observaciones de binarias de rayos X.
  • Se identifica un "vacío de pares" que se extiende desde el horizonte hasta unos pocos radios gravitacionales en el plano medio. Esto se debe a que, en esta región densa, la profundidad óptica supera el umbral crítico, suprimiendo la producción neta y permitiendo solo la aniquilación rápida.
  • La fracción de pares alcanza su máximo en una tira delgada en la base de la corona, justo fuera de la región de inestabilidad (plunging region), donde $\tau_p \approx 1$.

• Equilibrio vs. Desequilibrio y el Rol de la Advección:

  • Regiones de Alta Densidad (Plano Medio): Los pares están cerca del equilibrio local debido a tiempos de equilibrio cortos y alta profundidad óptica.
  • Regiones de Baja Densidad (Corona Superior y Jets): La fracción de pares puede exceder el valor de equilibrio local por órdenes de magnitud. Esto no se debe a producción local, sino a la advección: los pares generados en la tira de la base de la corona son transportados hacia arriba por el flujo global antes de tener tiempo de aniquilarse.
  • La advección domina sobre los procesos locales de pares en la corona y los jets, actuando como el mecanismo principal de inyección de pares.

• Dependencia de Parámetros:

  • Tasa de Acreción: A mayor tasa de acreción, la tira de máxima producción de pares se desplaza hacia el agujero negro y aumenta su ángulo de apertura, facilitando la captura de pares por los vientos y chorros.
  • Masa del Agujero Negro: Los resultados son insensibles a la masa del agujero negro (escala de 7 órdenes de magnitud probada), sugiriendo que los hallazgos son aplicables tanto a AGN como a binarias de rayos X.

• Densidad de Carga en los Jets:

  • La densidad de pares a lo largo de los polos magnéticamente dominados supera la densidad de Goldreich-Julian ($n_{GJ}$) en algunos modelos. Esto sugiere que la inyección de masa de pares por advección desde el disco podría ser un mecanismo viable para apantallar el campo eléctrico en los chorros, incluso en agujeros negros de baja luminosidad.

• Mecanismos de Producción:

  • Cerca del plano medio, la producción de pares está dominada por colisiones fotón-fotón.
  • En la corona superior y los jets (baja $\tau_p$), la producción es dominada por colisiones partícula-partícula, aunque la tasa absoluta es baja.

5. Significado e Implicaciones

• Regulación de la Temperatura de la Corona: Los resultados sugieren que la física de pares podría actuar como un "termostato" que regula la temperatura del plasma. En regiones donde la escala de tiempo de equilibrio es comparable a la de enfriamiento por colisiones de Coulomb, la producción de pares podría limitar la temperatura máxima del plasma.
• Origen de Pares en Binarias de Rayos X: El mecanismo de advección propuesto ofrece una explicación plausible para la presencia de positrones en ciertas binarias de rayos X (como 1E 1740.7–2942) y sus firmas de rayos gamma, sin necesidad de mecanismos de aceleración en brechas locales en el jet.
• Limitaciones y Futuro: El estudio reconoce limitaciones, como la temperatura electrónica fija y la falta de un modelo de transferencia radiativa autoconsistente (se usa una función de enfriamiento ad hoc). Sin embargo, establece una base sólida para futuras simulaciones que incluyan termodinámica de dos temperaturas y transferencia radiativa completa.

Conclusión:
El estudio demuestra que la distribución de pares electrón-positrón en flujos de acreción no está determinada únicamente por el equilibrio local, sino por una competencia compleja entre la producción local (controlada por la profundidad óptica) y el transporte global (advección). La advección desde una región de producción eficiente en la base de la corona es un mecanismo crítico para alimentar los chorros y la corona superior con pares, superando en muchos casos los valores de equilibrio local.