Bridging Roche Lobe Overflow and micro-TDEs: The Runaway Evolution of Eccentric Mass Transfer in Star-Black Hole Binaries

Resumen Técnico: Cerrando la brecha entre el desbordamiento del lóbulo de Roche y los micro-TDEs

Planteamiento del problema
Los sistemas binarios que contienen agujeros negros de masa estelar (sBH) y estrellas en órbitas excéntricas presentan una brecha teórica entre dos regímenes establecidos: el desbordamiento de Roche (RLOF) estable en binarias circulares y los eventos de disrupción de marea (TDE) violentos de un solo paso. Mientras que los modelos semianalíticos han intentado describir la evolución secular de las binarias excéntricas, estos dependen de ecuaciones promediadas en la órbita y tasas de transferencia de masa parametrizadas que pueden fallar al capturar la respuesta estructural no lineal de la estrella. Por el contrario, las simulaciones hidrodinámicas previas se han centrado principalmente en la fase terminal de encuentros de marea profundos o pasajes únicos de pericentro, omitiendo las trayectorias evolutivas a largo plazo que determinan si un sistema se establece en una transferencia de masa estable o si sufre una disrupción descontrolada. La región de transición crítica, gobernada por la competencia entre la expansión estelar adiabática y el ensanchamiento orbital, sigue siendo insuficientemente comprendida.

Metodología
Los autores emplean la Hidrodinámica de Partículas Suavizadas (SPH) utilizando el código Phantom para simular la interacción entre una estrella de la secuencia principal de $1 M_\odot$ (similar al Sol) y un agujero negro de masa estelar de $10 M_\odot$. Las simulaciones utilizan una ecuación de estado adiabática ($\gamma = 5/3$) y omiten el enfriamiento radiativo, un supuesto justificado dado que la escala de tiempo de la interacción (periodos orbitales) es órdenes de magnitud más corta que la escala de tiempo térmica estelar. Esto asegura que la disipación de marea y el calentamiento por choque queden atrapados dentro de la envoltura estelar, impulsando la expansión adiabática.

El estudio rastrea sistemas durante decenas o más de 100 periodos orbitales ($P_0$), iniciando la simulación notablemente fuera del radio de marea para capturar el inicio de la interacción. El espacio de parámetros se explora variando la excentricidad inicial ($e_0 = 0.30, 0.55, 0.70$) y la distancia de pericentro adimensional ($b_0 = r_{p,0}/r_{tide,0} = 3.33, 3.45, 3.57$). Los autores distinguen entre material ligado y no ligado mediante un criterio del parámetro de Bernoulli que incluye la entalpía específica, y rastrean la evolución orbital utilizando métodos tanto dinámicos (energía/momento angular) como geométricos.

Resultos clave
Las simulaciones revelan dos trayectorias evolutivas distintas dictadas por la distancia inicial del pericentro ($b_0$):

1. Disrupción descontrolada (ej., Run A, $b_0 = 3.33$):

   • Mecanismo: La pérdida de masa en el pericentro desencadena la expansión adiabática de la envoltura estelar. Debido a que el exponente de masa-radio adiabático de la estrella ($\zeta_{ad}$) es más negativo que la respuesta del lóbulo de Roche ($\zeta_L$), el radio estelar se expande más rápido que el lóbulo de Roche.
   • Bucle de retroalimentación: Esto crea un bucle de retroalimentación positiva donde el aumento del factor de llenado del lóbulo de Roche impulsa mayores tasas de pérdida de masa, lo que a su vez impulsa una mayor expansión.
   • Resultado: El sistema experimenta una disrupción catastrófica en aproximadamente $\sim 40 P_0$. Los escombros forman un toro de acreción térmicamente calentado por choques y geométricamente grueso.
   • Acreción: El agujero negro experimenta una tasa de acreción super-Eddington que alcanza un pico de $\sim 5 \times 10^6 \dot{M}_{Edd}$, seguido de un decaimiento que escala como $t^{-9/4}$.

2. Transferencia de masa estable (ej., Run B, $b_0 = 3.57$):

   • Mecanismo: Aunque la pérdida de masa impulsa inicialmente la expansión estelar, la pérdida asociada de momento angular orbital provoca un ensanchamiento significativo de la órbita binaria.
   • Regulación: El ensanchamiento orbital amplía el lóbulo de Roche lo suficiente como para contrarrestar la expansión estelar. El sistema alcanza un equilibrio cuasi-estable donde la tasa de pérdida de masa se estabiliza y el factor de llenado permanece por debajo de la unidad.
   • Resultado: El sistema persiste en una fase de transferencia de masa de larga duración (simulada hasta $150 P_0$) sin disrupción total.
   • Acreción: La tasa de acreción se mantiene, pero es menor, alcanzando un pico alrededor de $10^5 \dot{M}_{Edd}$, exhibiendo una modulación cuasi-periódica.

Sensibilidad del umbral
La transición entre estos resultados es sensible a $b_0$. La Carrera C ($b_0 = 3.45$) demuestra una "disrupción descontrolada retardada", donde el sistema imita una transferencia estable durante $\sim 70 P_0$ antes de que el bucle de retroalimentación positiva eventualmente domine. Las variaciones en la excentricidad ($e_0$) para un $b_0$ fijo muestran que las excentricidades más bajas en $b_0=3.33$ conducen a una disrupción más temprana, mientras que $b_0=3.57$ permanece estable incluso con excentricidades más bajas ($e_0=0.30$).

Significado y afirmaciones
El artículo afirma proporcionar una imagen física unificada que conecta el desbordamiento de Roche excéntrico con el inicio de los "micro-TDEs" (eventos de disrupción de marea que involucran agujeros negros de masa estelar). Los autores argumentan que la evolución a largo plazo de las binarias excéntricas no está predeterminada, sino que depende de la competencia específica entre la respuesta estructural estelar y la dinámica orbital.

• Implicaciones observacionales: Se propone que el canal de disrupción descontrolada (Run A) alimenta transitorios rápidos de rayos X/UV o azules/ópticos mediante la acreción super-Eddington sobre el sBH, manifestándose potencialmente como fuentes de rayos X ultra luminosas (ULX) aparentes. El canal de transferencia de masa estable (Run B) se sugiere que se manifestará como llamaradas repetitivas y cuasi-periódicas.
• Contribución teórica: El trabajo cierra la brecha entre los modelos seculares semianalíticos y las simulaciones de encuentros únicos violentos, demostrando que el tratamiento hidrodinámico completo es necesario para capturar la ruptura de las aproximaciones durante la transición de la transferencia de masa estable a la inestable.

Los autores mantienen la modestia respecto a sus conclusiones, señalando que el umbral de transición está moldeado tanto por $b_0$ como por $e_0$, en lugar de un único valor universal, y que se requiere trabajo futuro que incorpore transferencia radiativa y campos magnéticos para conectar plenamente estos resultados hidrodinámicos con las observaciones multi-longitud de onda.