Radiation-hydrodynamics of star-disc collisions: From system parameters to outflows and lightcurves

Resumen Técnico: Hidrodinámica de radiación de colisiones entre estrella y disco: desde los parámetros del sistema hasta los flujos y curvas de luz

Problema y Contexto
Las erupciones cuasi-periódicas (QPEs) son transitorios nucleares asociados con agujeros negros supermasivos (SMBHs), caracterizados por destellos repetitivos de rayos X blandos. Aunque existen varios modelos, el escenario de colisión estrella-disco —donde una estrella en una órbita estrecha, excéntrica e inclinada intersecta repetidamente un disco de acreción— es un candidato principal. Este modelo explica naturalmente los tiempos de recurrencia, los intervalos alternos "largos-cortos" y los patrones de destello "fuertes-débiles" observados en fuentes como GSN 069. Sin embargo, conectar los parámetros físicos locales de estas colisiones con las propiedades observables de los destellos (luminosidad, duración, asimetría) requiere comprender la evolución acoplada de los choques, la geometría de los flujos y el transporte de radiación en un medio en expansión y ópticamente grueso. Estudios previos han explorado partes de este problema utilizando enfoques 1D, 2D o no radiativos, pero se carecía de un estudio sistemático de hidrodinámica de radiación en 3D que vinculara los parámetros del sistema con las curvas de luz.

Metodología
Los autores realizaron una serie de simulaciones de hidrodinámica de radiación tridimensional local utilizando el código de hidrodinámica de partículas suavizadas (SPH) Phantom. Las simulaciones se llevaron a cabo en un marco cartesiano local centrado en el plano medio del disco, omitiendo la dinámica orbital global, la relatividad general y la autogravedad del gas, ya que las escalas de tiempo de la colisión son mucho más cortas que las escalas de tiempo orbitales.

El estudio utilizó un modelo fiduciario basado en una estrella de la secuencia principal ($R_\star \approx R_\odot$) colisionando con un disco-$\alpha$ ópticamente grueso y dominado por la presión de radiación alrededor de un SMBH de $10^6 M_\odot$. La estrella se modeló como una esfera rígida moviéndose a $v_\star = 0.1c$. Los autores variaron sistemáticamente cinco parámetros clave respecto al modelo fiduciario, variando uno a la vez:

1. Velocidad estelar ($v_\star$): Desde $0.5$ hasta $1.5$ veces el valor fiduciario.
2. Densidad superficial del disco ($\Sigma_d$): Desde $1$ hasta $100$ veces el valor fiduciario, abarcando desde discos delgados estándar hasta discos densos post-evento de disrupción de marea (TDE).
3. Radio estelar ($R_\star$): Desde $0.5$ hasta $4$ veces el valor fiduciario, variando efectivamente la relación entre la escala de altura del disco y el radio estelar ($H/R_\star$).
4. Perfil de densidad vertical del disco: Modelado como una distribución gaussiana $\rho(z) \propto e^{-z^2/2\sigma^2}$, variando la concentración de masa hacia el plano medio ($\sigma/R_\odot$).
5. Ángulo de colisión local ($i$): El ángulo entre el vector de velocidad relativa y el plano medio del disco, que va desde $90^\circ$ (perpendicular) hasta $30^\circ$ (oblicuo).

El transporte de radiación se trató mediante la aproximación de difusión de flujo limitado con opacidad de dispersión de electrones. La luminosidad bolométrica emergente se calculó integrando los flujos de difusión y advección a través de una fotosfera determinada dinámicamente.

Contribuciones Clave y Resultados
El estudio proporciona escalamientos empíricos para la luminosidad máxima ($L_{\text{peak}}$) y la duración del destello ($\Delta t$) como funciones de los parámetros del sistema. Los principales hallazgos son:

• Dinámica de la colisión: La estrella impulsa un choque de proa (bow shock), creando dos flujos (hacia adelante y hacia atrás) en lados opuestos del disco. El flujo hacia adelante es generalmente más masivo y energético. La morfología cualitativa del choque de proa y los flujos permanece mayormente sin cambios a través de las variaciones en $v_\star$ y $\Sigma_d$, aunque la normalización termodinámica cambia.
• Velocidad estelar ($v_\star$): Las estrellas más rápidas producen destellos más brillantes ($L_{\text{peak}} \propto v_\star^2$) debido a tasas de calentamiento por choque más altas. También producen destellos más largos ($\Delta t \propto v_\star^{1.4-2.2}$) porque se inyecta más energía, aunque la luminosidad decae más rápido debido a pérdidas adiabáticas más fuertes.
• Densidad superficial del disco ($\Sigma_d$): Los discos más densos aumentan la masa y el momento de los flujos pero tienen un efecto débil en la luminosidad máxima ($L_{\text{peak}} \propto \Sigma_d^{-0.1}$ para el flujo hacia adelante, $\propto \Sigma_d^{0.2}$ para el hacia atrás). Esto se debe a una competencia entre el aumento del calentamiento por choque y las pérdidas adiabáticas más fuertes en el eyecta ópticamente más grueso. Sin embargo, los discos más densos prolongan significativamente la duración del destello ($\Delta t \propto \Sigma_d^{0.5-0.6}$) debido a los tiempos de difusión de fotones más largos.
• Radio estelar ($R_\star$): Estrellas más grandes interactúan con más gas, produciendo destellos más brillantes ($L_{\text{peak}} \propto R_\star^{1.5}$ hacia adelante, $\propto R_\star^{0.7}$ hacia atrás) y más largos ($\Delta t \propto R_\star^{1.5}$). El aumento de $R_\star$ potencia la asimetría de momento entre los flujos hacia adelante y hacia atrás.
• Perfil de densidad vertical ($\sigma$): Discos más concentrados centralmente (menor $\sigma$) producen destellos más brillantes y cortos. Las capas centrales densas absorben más energía del choque, mientras que las capas exteriores de baja densidad permiten que la radiación escape más eficientemente con menores pérdidas adiabáticas.
• Ángulo de colisión ($i$): Los cruces más oblicuos (menor $i$) reducen la asimetría de momento y luminosidad entre los dos flujos. El flujo hacia atrás se vuelve significativamente más brillante ($L_{\text{peak}} \propto (\sin i)^{-0.8}$) y los destellos se vuelven más largos ($\Delta t \propto (\sin i)^{-0.9 \text{ a } -1}$) debido al tiempo de interacción más largo y una redistribución más simétrica del gas sacudido.

Aplicación a GSN 069
Los autores aplicaron sus escalamientos empíricos a la fuente bien caracterizada GSN 069. Al muestrear el espacio de parámetros para coincidir con las luminosidades máximas observadas ($\sim 3 \times 10^{42}$ erg s$^{-1}$), duraciones ($\sim 1$ h) y la relación de amplitud fuerte-débil, identificaron dos regiones viables:

1. Solución Preferida: Una órbita retrógrada colisionando con un disco post-TDE denso y verticalmente concentrado ($\Sigma_d \sim 45 \times$ el valor fiduciario, $\sigma/R_\odot \sim 0.3$). Este modelo reproduce bien las propiedades observadas de los destellos.
2. Solución Alternativa: Una órbita retrógrada con una densidad de disco delgado estándar pero con una geometría de colisión muy rasante ($i \sim 6^\circ$). Esta solución es menos precisa y requiere parámetros extremos.

Los autores señalan que la preferencia por un disco denso y verticalmente concentrado es consistente con la evidencia observacional de actividad tipo TDE en GSN 069.

Significancia y Limitaciones
El artículo afirma proporcionar el primer estudio sistemático de hidrodinámica de radiación en 3D que vincula los parámetros locales de colisión con las curvas de luz de las QPEs. Los escalamientos empíricos derivados sirven como un puente entre las simulaciones locales y los modelos globales de QPE, permitiendo la predicción de las propiedades de los destellos basadas en los parámetros del sistema. Los autores enfatizan que sus hallazgos sugieren que combinaciones específicas de parámetros pueden reproducir las amplitudes, duraciones y ciclos de trabajo características de las Q QPEs.

El estudio reconoce modestamente ciertas limitaciones: los modelos asumen escalamientos de un solo parámetro separables, omiten la rotación y cizalladura del disco, tratan a la estrella como una esfera rígida y no consideran los efectos del ángulo de visión o la evolución espectral (por ejemplo, Comptonización). Por consiguiente, los resultados se presentan como regiones físicamente motivadas del espacio de parámetros en lugar de modelos únicos y definitivos. Los autores sugieren que trabajos futuros deberían incorporar estructuras estelares realistas, dinámica de discos y luminosidades dependientes del ángulo de visión para cuantificar plenamente la capacidad del modelo para reproducir la diversidad de las fuentes de QPE observadas.