Radiation-hydrodynamics of star-disc collisions for quasi-periodic eruptions

Resumen Técnico: Hidrodinámica de Radiación de Colisiones Estrella-Disco para Erupciones Cuasi-Periódicas

Planteamiento del Problema
Las erupciones cuasi-periódicas (QPE, por sus siglas en inglés) son transitorios de rayos X luminosos que ocurren cerca de agujeros negros supermasivos (SMBH), caracterizados por estallidos brillantes separados por horas o días. Un modelo teórico principal postula que estas erupciones resultan de una estrella en una órbita inclinada que colisiona repetidamente con un disco de acreción. Aunque estudios previos han explorado las interacciones estrella-disco, muchos carecieron de un transporte de radiación detallado o se centraron en regímenes donde la estrella es significativamente perturbada o el tiempo de cruce es demasiado corto para que el gas fluya alrededor de la estrella. Este estudio aborda la necesidad de un tratamiento tridimensional (3D) de la hidrodinámica de radiación de la colisión en el régimen donde la estrella permanece mayormente inalterada y el gas impactado tiene tiempo suficiente para fluir alrededor de ella.

Metodología
Los autores realizaron una simulación de hidrodinámica de radiación 3D utilizando el código de partículas suavizadas (smoothed particle hydrodynamics, SPH) Phantom. La simulación se llevó a cabo en un marco de referencia que se mueve conjuntamente con el disco de acreción para aislar la dinámica local de la colisión.

• Configuración: Una estrella esférica rígida ($R_\star = R_\odot$) se mueve perpendicularmente a través de un disco de acreción uniforme, ópticamente grueso y dominado por la presión de radiación ($H = 3R_\star$, $\rho_d = 3.7 \times 10^{-8}$ g cm$^{-3}$) a una velocidad de $v_\star = 0.1c$. Los parámetros del sistema corresponden a un SMBH de $10^6 M_\odot$ y un periodo de QPE de 4 horas.
• Física: El gas sigue una ecuación de estado adiabática ($\Gamma = 5/3$). El transporte de radiación se trata mediante la aproximación de difusión de flujo limitado, asumiendo un equilibrio termodinámico local (LTE). La estrella se modela como un cuerpo rígido que interactúa elásticamente con las partículas de gas; se descartan las fuerzas gravitatorias del SMBH y la autogravedad, ya que la escala de tiempo de la colisión es mucho más corta que las escalas de tiempo orbitales.
• Alcance: La simulación se centra en una sección localizada del disco ($12R_\star \times 2H$) para resolver la estructura del choque y la formación de la salida (outflow) sin el costo computacional de una simulación global.

Contribuciones Clave
Este trabajo presenta la primera simulación de hidrodinámica de radiación 3D de una colisión estrella-disco dirigida específicamente al espacio de parámetros relevante para la mayoría de las QPE (donde el disco es lo suficientemente grueso para que el gas fluya alrededor de la estrella). El estudio proporciona un análisis detallado de:

1. La formación y evolución del choque de proa (bow shock) y los choques secundarios.
2. La redistribución asimétrica del momento y la energía inyectados por la estrella.
3. Las propiedades distintivas de los flujos de salida (outflows) hacia adelante y hacia atrás.
4. Las curvas de luz bolométricas y temperaturas efectivas resultantes.

Resultados

• Dinámica de Choque: La colisión genera un choque de proa casi paraboloide que viaja con la estrella y se expande lateralmente. Un choque secundario se forma en la estela de la estrella a medida que las corrientes de gas convergen. La expansión lateral del frente de choque sigue una escala de ley de potencia ($r_\perp \propto t^{0.4}$), ligeramente más superficial que la predicción de estado estacionario de $t^{0.5}$, lo que indica que el sistema aún está evolucionando hacia un estado estacionario.
• Flujos de Salida Asimétricos: La colisión produce dos flujos de salida distintos en lados opuestos del disco:
  • Flujo hacia adelante (dirección del movimiento): Este flujo es impulsado por una ruptura de choque (shock breakout) casi esférica a medida que la estrella sale del disco. Se expande más rápido, contiene significativamente más masa (aproximadamente 10 veces el flujo hacia atrás), transporta más energía y momento, y es aproximadamente dos veces más luminoso que el flujo hacia atrás.
  • Flujo hacia atrás: Este flujo carece de un componente de ruptura de choque casi esférica. Es formado por gas que fluye alrededor de la estrella y es redirigido por el choque secundario o escapa desde el borde superior del disco. Es menos masivo y menos luminoso.
• Luminosidad y Curvas de Luz: El flujo hacia adelante alcanza una luminosidad bolométrica máxima de $\sim 10^{41}$ erg s$^{-1}$, aproximadamente el doble que la del flujo hacia atrás. La curva de luz hacia adelante exhibe un ascenso rápido seguido de una ruptura hacia un declive más lento, interpretado como un cambio de la radiación que emerge de los ejectas iniciales de la ruptura del choque a la radiación del gas que escapa de la cavidad tras un retraso mayor. Las temperaturas efectivas alcanzan su máximo en $\sim 29$ eV para el flujo hacia adelante y valores ligeramente menores para el flujo hacia atrás.
• Inyección de Momento: La estrella inyecta momento preferencialmente a lo largo de su dirección de movimiento. El momento vertical total del gas impactado asciende asintóticamente a $P_z \approx -0.45 M_d v_\star$, lo que refleja la naturaleza anisotrópica de la inyección de energía y los ángulos oblicuos de reflexión del gas.

Significado y Reivindicaciones
Los autores afirman que la asimetría intrínseca en la dinámica de la colisión explica naturalmente los patrones de destellos "fuertes-débiles" observados en varias fuentes de QPE (por ejemplo, GSN 069). Debido a que el flujo hacia adelante es significativamente más luminoso que el flujo hacia atrás, un observador que vea el sistema desde un ángulo específico vería destellos alternos brillantes y tenues correspondientes a los pases sucesivos de la estrella a través del disco.

El estudio concluye que el modelo de colisión estrella-disco, en el régimen donde la estrella permanece inalterada, puede reproducir las luminosidades pico observadas ($10^{41}-10^{43}$ erg s$^{-1}$), las escalas de tiempo (horas) y los patrones de luminosidad alternos de las QPE. Los autores señalan limitaciones, incluyendo el uso de una estrella rígida (omitiendo el desprendimiento de masa y la deformación), la suposición de LTE (que puede fallar en regímenes con escasez de fotones) y la omisión de la rotación global del disco y la gravedad del SMBH. Sugieren que trabajos futuros que incorporen física de no-LTE y estructuras estelares realistas son necesarios para constreñir completamente el modelo y abordar las discrepancias en la dureza espectral (simulado $\sim 30$ eV frente a observado $\sim 100-200$ eV).