Resolving the (Debate About) Nozzle Shocks in Tidal Disruption Events

Resumen Técnico: Resolución de la (discusión sobre los) choques de boquilla en eventos de disrupción de marea

Planteamiento del problema
Los eventos de disrupción de marea (TDE, por sus siglas en inglés) ocurren cuando una estrella pasa dentro del límite de Roche de un agujero negro supermasivo (SMBH), lo que resulta en la disrupción de la estrella y la formación de un flujo de escombos ligados. Mientras que la posterior circularización y acreción de estos escombos impulsan las llamaradas observables, el mecanismo de circularización sigue siendo objeto de debate. Una fase crítica, aunque mal resuelta, ocurre en el pericentro durante el retorno del flujo de escombos, donde el enfoque balístico crea una compresión vertical extrema, formando un "choque de boquilla" (nozzle shock). Simulaciones 3D previas han arrojado resultados contradictorios respecto al impacto de este choque: algunas sugieren que impulsa una rápida expansión y circularización del flujo, mientras que otras encuentran una expansión mínima. Estas discrepancias se atribuyen a menudo a la disipación numérica en las simulaciones 3D, que no logran resolver la extrema compresión vertical del choque (escalas de $\sim 10^{-4} R_\odot$). Además, estudios previos han omitido en gran medida los efectos termodinámicos de la recombinación de Hidrógeno (H) y la formación de hidrógeno molecular (H$_2$), que pueden alterar significativamente el presupuesto energético y la ecuación de estado (EOS) de los escombos.

Metodología
Para resolver de manera inequívoca el choque de boquilla e incorporar una termodinámica realista, los autores desarrollaron un marco de modelado híbrido computacionalmente eficiente que combina tres enfoques distintos:

1. Simulaciones 3D SPH (phantom): Utilizadas para simular la disrupción inicial de una estrella de masa solar (politropo $n=3$) por un SMBH de $10^6 M_\odot$ con un factor de penetración $\beta=2$. Esta fase captura la formación del flujo de escombos y su dispersión de energía inicial.
2. Modelo Afín Semi-analítico: Utilizado para evolucionar el flujo de escombos en la primera órbita. Este modelo asume una evolución homóloga y rastrea las dimensiones del flujo (anchura vertical $H$, anchura en el plano $\Delta$) y la entropía. Se inicializa con instantáneas de la simulación SPH.
3. Hidrodinámica de Volumen Finito 1D (Athena++): Utilizada para resolver el choque de boquilla. A medida que el flujo se acerca al pericentro, el modelo afín falla debido a la formación del choque. Los autores transicionan a simulaciones 1D que representan la columna central de la sección transversal del flujo, acopladas a la evolución en el plano del modelo afín.

Contribuciones e Innovaciones Clave

• Ecuación de Estado (EOS) Realista: A diferencia de la mayoría de los estudios de TDE que asumen una ley $\gamma$ simple ($\gamma=5/3$), este trabajo implementa una EOS compleja (basada en Tomida et al. 2013) que tiene en cuenta la recombinación de H, la formación/disociación de H$_2$, la ionización de He y la presión de radiación.
• Resolución del Choque de Boquilla: Al utilizar simulaciones 1D con alta resolución ($\Delta z \sim 10^{-7} R_\odot$), los autores evitan la disipación numérica inherente a las simulaciones globales 3D, permitiendo un cálculo preciso de la disipación de energía en el choque.
• Acoplamiento Híbrido: La transición fluida entre el modelo afín, la física del choque 1D y el regreso al modelo afín permite el seguimiento de las propiedades del flujo desde la disrupción a través de la segunda órbita y hasta el punto de autointersección.

Resultos Clave

1. Ensanchamiento del Flujo Pre-Boquilla: Los procesos químicos alteran significativamente la evolución del flujo antes de alcanzar la boquilla. La recombinación de H ensancha el flujo por un factor de $\sim 5$, y la formación de H$_2$ añade un factor adicional de orden unitario. Este ensanchamiento es impulsado por la liberación de calor latente, lo que mantiene temperaturas altas y reduce el índice adiabático efectivo, potenciando los gradientes de presión.
2. Dinámica del Choque de Boquilla: El flujo experimenta una compresión vertical extrema hasta $\sim 10^{-4} R_\odot$ antes de que los gradientes de presión reviertan el colapso. El choque ioniza el gas y disipa la energía cinética, aumentando la entropía específica. La disipación es significetivamente mayor para la EOS realista en comparación con el caso de la ley $\gamma$, debido al ensanchamiento previo en la boquilla impulsado por la recombinación.
3. Evolución Post-Boquilla y Puntos Focales: La disipación en la boquilla altera la trayectoria y el perfil de presión del flujo. Crucialmente, los gradientes de presión incrementados (debido a la disipación y los procesos químicos) resisten el colapso en el segundo punto focal (ubicado a radios mayores que la boquilla). Esto conduce a una expansión significativa del flujo en la segunda órbita, haciendo que el ancho del flujo en la autointersección sea altamente sensible a la cantidad de disipación en la boquilla.
4. Geometría de Autointersección:
   • Desplazamiento Vertical: El ensanchamiento del flujo aumenta la probabilidad de autointersección a pesar de la precesión nodal relativista. Sin embargo, para giros altos del BH ($\chi \gtrsim 0.6$), los desplazamientos verticales aún podrían prevenir la intersección.
   • Desajuste de Anchura: Se encuentra que el flujo saliente es más delgado que el flujo entrante por un factor de pocos. Esto sugiere que solo una fracción de los escombros entrantes participa en el choque de autointersección, reduciendo la energía disponible para la disipación.
   • Precesión Apsidal Diferencial (DAP): Los autores señalan que la DAP probablemente domina el ancho en el plano en la autointersección, un factor que no se incluyó completamente en su modelo actual pero que se identifica como una consideración crítica para el futuro.
5. Implicaciones de Emisión: La luminosidad generada en la autointersección, calculada para el caso optimista de flujos idénticos con desplazamiento vertical cero, está ligeramente por debajo de las luminosidades pico observadas en TDE ($\sim 10^{42}$ erg/s vs. $10^{44}$ erg/s). El ensanchamiento derivado de los procesos químicos solo aumenta ligeramente esta luminosidad. Los autores sugieren que la circularización secundaria y los choques de acreción en el flujo de salida pueden proporcionar la contribución dominante a la emisión observada.

Significancia y Reivindicaciones
El artículo afirma resolver el debate sobre los choques de boquilla demostrando que los resultados conflictivos previos fueron probablemente artefactos de la resolución numérica y la omisión de la termodinámica química. Los autores aseveran que:

• El choque de boquilla no circulariza directamente los escombos mediante gradientes de presión en el plano; más bien, modifica las propiedades del flujo (anchura y entropía).
• Estas propiedades modificadas aumentan sustancialmente la probabilidad de autointersección del flujo en la segunda órbita, incluso en presencia de precesión relativista.
• El papel de la recombinación de H y la formación de H$_2$ es esencial para modelar con precisión la dinámica del flujo de escombos y la circularización temprana.

Este trabajo proporciona un modelo idealizado fundacional para la evolución del flujo de escombos que clarifica el papel del choque de boquilla, ofreciendo una base más realista para futuros modelos de circularización y emisión. Los autores declaran explícitamente que sus resultados están limitados por la exclusión de la precesión apsidal diferencial y los campos magnéticos, los cuales deben ser analizados en trabajos futuros.