Mass Transfer in Tidally Heated Stars Orbiting Massive Black Holes and Implications for Repeating Nuclear Transients

Resumen Técnico: Transferencia de Masa en Estrellas Calentadas por Mareas que Orbitan Agujeros Negros Masivos

Planteamiento del Problema
Las estrellas que orbitan cerca de agujeros negros supermasivos (SMBH, por sus siglas en inglés) están sujetas a interacciones complejas que involucran dinámica gravitatoria e hidrodinámica. Si bien los eventos de disrupción de marea (TDE) ocurren cuando las estrellas pasan dentro del radio de marea ($r_t$), observaciones recientes de transitorios nucleares repetitivos (RNT) —incluyendo erupciones cuasi-periódicas (QPE) y candidatos a pTDE parciales— sugieren que las estrellas pueden sobrevivir a múltiples pasos por el pericentro. En estos escenarios, incluso si la distancia del pericentro ($r_p$) permanece fuera de $r_t$, la estrella experimenta un calentamiento por marea significativo. Trabajos previos indican que este calentamiento puede inflar el radio estelar, iniciando potencialmente una transferencia de masa hacia el SMBH antes de que la órbita decaiga vía ondas gravitacionales (GW) hasta el punto de disrupción. Sin embargo, la dinámica acoplada del calentamiento por marea, la evolución estructural estelar, el decaimiento orbital y la estabilidad de la transferencia de masa resultante sigue siendo poco comprendida, particularmente en lo que respecta a si la transferencia de masa es estable o inestable y qué luminosidad produce.

Metodología
Los autores emplean un enfoque dual que combina modelos analíticos y simulaciones numéricas utilizando Modules for Experiments in Stellar Astrophysics (MESA).

1. Modelos Analíticos: Los autores desarrollan un marco semi-analítico para una estrella similar al Sol (en la secuencia principal con un núcleo radiativo y una envoltura convectiva) que orbita un SMBH ($M_\bullet \gg M_\star$). Consideran dos regímenes orbitales: altamente excéntrico ($1-e \ll 1$) y casi circular ($e \ll 1$).

   • Calentamiento por Marea: Utilizan la teoría de perturbación lineal (Press & Teukolsky 1977; Lai 1997) para órbitas altamente excéntricas y disipación de marea dinámica mediante ondas de gravedad internas para órbitas casi circulares (Linial & Quataert 2024b). La tasa de calentamiento se parametriza mediante factores adimensionales ($f_p$ y $f_c$) para dar cuenta de las incertidumbres en el acoplamiento de modos no lineales y la disipación.
   • Ecuaciones de Evolución: El modelo rastrea la evolución del radio estelar ($R_\star$) y la masa ($M_\star$) impulsada por el equilibrio entre la deposición de energía de marea y la energía de ligadura estelar. Crucialmente, asumen que el momento angular y la energía perdidos mediante la transferencia de masa son eliminados permanentemente del sistema (no devueltos a la órbita), una condición apropiada para binarias con una alta relación de masas.
   • Criterios de Estabilidad: La estabilidad se evalúa mediante la respuesta estelar a la pérdida de masa. Si $R_\star \propto M_\star^\zeta$ con $\zeta \lesssim 1/3$, la transferencia de masa es inestable; de lo contrario, es estable.

2. Simulaciones MESA: Para abordar las limitaciones de los modelos analíticos respecto a la estructura interna estelar, los autores realizan simulaciones 1D con MESA.

   • Configuración: Los modelos consisten en una estrella de la secuencia principal de edad cero (ZAMS) de $1 M_\odot$ orbitando un SMBH.
   • Implementación del Calentamiento por Marea: Dado que la distribución espacial de la disipación de marea es incierta, prueban dos perfiles de calentamiento fenomenológicos: (1) Concentrado Centralmente ($0.1 < r/R_\star < 0.2$), motivado por la ruptura de ondas de gravedad internas, y (2) Uniforme por unidad de masa en el 75% exterior de la estrella.
   • Evolución Orbital: Incorporan el decaimiento orbital por GW (Peters 1964) y términos de evolución orbital por marea en el código de evolución binaria de MESA.
   • Determinación de la Estabilidad: Las simulaciones rastrean explícitamente si el radio estelar se expande o se contrae en respuesta a la pérdida de masa, determinando así la estabilidad de la fase de transferencia de masa.

Contribuciones Clave y Resultados

• Dominancia de la Escala de Tiempo de Calentamiento por Marea: El hallazgo principal es que, para estrellas que orbitan SMBHs, la tasa de transferencia de masa estable ($\dot{M}_\star$) está típicamente determinada por la escala de tiempo de calentamiento por marea ($\tau_{\text{heat}} \sim E_\star / \dot{E}_{\text{tides}}$) más que por la escala de tiempo de decaimiento orbital por GW ($\tau_{\text{GW}}$).

  • Para órbitas altamente excéntricas, $\dot{M}_\star \sim 10^{-5} - 10^{-3} M_\odot \text{yr}^{-1}$.
  • Para órbitas casi circulares, $\dot{M}_\star \sim 10^{-6} - 10^{-4} M_\odot \text{yr}^{-1}$ (potencialmente mayor si la excitación de ondas inerciales es eficiente).
  • Esta tasa es significativamente mayor que la predicha únicamente por el decaimiento orbital impulsado por GW, lo que implica una fase de transferencia de masa estable mucho más corta ($\sim 10^3 - 10^4$ años) comparada con los escenarios impulsados por GW ($\sim 10^6 - 10^7$ años).

• Dependencia de la Estabilidad según la Ubicación del Calentamiento: La estabilidad de la transferencia de masa es altamente sensible a dónde se disipa la energía de marea dentro de la estrella.

  • Calentamiento Concentrado Centralmente: Provoca que la estrella se vuelva totalmente convectiva. En los modelos de MESA, esto conduce a una transferencia de masa inestable, caracterizada por una expansión descontrolada y un aumento brusco de $\dot{{M}}_\star$, que el código no puede resolver numéricamente.
  • Calentamiento Uniforme/Exterior: Aumenta la entropía en las capas exteriores, inhibiendo la convección y manteniendo un núcleo radiativo. Esto conduce a una transferencia de masa estable, donde la estrella se contrae o se expande lo suficientemente lento como para mantener el equilibrio, permitiendo una fase sostenida de pérdida de masa.

• Evolución Orbital: Aunque el calentamiento por marea infla la estrella, la evolución orbital impulsada por las mareas es generalmente subdominante al decaimiento por GW antes del inicio de la transferencia de masa para los parámetros fiduciarios. Sin embargo, una vez que comienza la transferencia de masa, la evolución orbital es lenta en comparación con la escala de tiempo de pérdida de masa estelar.

Significancia e Implicaciones
El artículo argumenta que la transferencia de masa calentada por marea ofrece un mecanismo viable para alimentar una población de núcleos galácticos activos (AGN) de baja luminosidad y ciertos transitorios nucleares repetitivos.

• Transitorios Nucleares Repetitivos (RNTs):

  • Candidatos a pTDE: El modelo proporciona una versión más "suave" de los TDE parciales para estrellas en órbitas donde $r_p$ se aproxima a $r_t$ lentamente (por ejemplo, ASASSN-14ko). Los autores señalan que sus modelos por defecto predicen tasas ligeramente menores a las requeridas para ASASSN-14ko, lo que sugiere que muchos de estos sistemas podrían ser más tenues de lo detectado actualmente.
  • QPEs: La fase de transferencia de masa estable, particularmente para órbitas casi circulares con tasas de calentamiento mejoradas (por ejemplo, mediante ondas inerciales), podría alimentar QPEs de período corto (por ejemplo, GSN 069, eRO-QPE2). El modelo explica cualitativamente la emisión de disco quiescente de larga duración y la falta de líneas anchas de AGN.
  • AGN de Baja Luminosidad: Las tasas de transferencia de masa predichas son suficientes para alimentar AGN de baja luminosidad en muchas galaxias, representando potencialmente una población significativa pero actualmente no detectada.

• Avance Teórico: Este trabajo destaca un bucle de retroalimentación crítico: el calentamiento por marea altera la estructura estelar, la cual a su vez modifica la tasa de calentamiento por marea y la estabilidad de la transferencia de masa. Los autores enfatizan que las conclusiones robustas requieren cálculos consistentes de este acoplamiento, que sus modelos 1D actuales aproximan pero no resuelven completamente.

Limitaciones
Los autores declaran explícitamente las siguientes limitaciones:

• Física del Calentamiento: La tasa de calentamiento total y la distribución espacial de la disipación son inciertas y tratadas de forma fenomenológica.
• Aproximaciones Orbitales: Los modelos analíticos asumen una evolución orbital simplificada, y los modelos MESA aproximan la deposición de energía de marea como ráfagas promediadas en el tiempo en lugar de ráfagas impulsivas en el pericentro.
• Hidrodinámica: Los modelos son hidrostáticos y no capturan las pulsaciones dependientes del tiempo o la eyección dinámica de masa que podrían ocurrir durante fases no lineales.
• Interacciones con el Disco: Se ignora la interacción entre la estrella y los discos de acreción preexistentes (de otros TDE), aunque se reconoce como un factor potencial para la supervivencia estelar.

En conclusión, el artículo establece que el calentamiento por marea puede iniciar la transferencia de masa en estrellas que orbitan SMBHs a tasas y escalas de tiempo distintas de la evolución impulsada por GW, y que la estabilidad de esta transferencia depende críticamente de la disipación interna de la energía de marea. Este mecanismo ofrece una explicación plausible para un subconjunto de transitorios repetitivos observados y sugiere una población mayor de sistemas más tenues y no detectados.