Self-Limited Accretion onto Embedded Binaries in a Uniform Medium

Resumen Técnico: Acreción Auto-Limitada en Binarias Incrustadas

1. Planteamiento del Problema

El estudio aborda la pregunta fundamental sobre cómo acreta gas una binaria de masas iguales inmersa en un medio gaseoso uniforme en reposo. El problema se define por la relación entre el radio de Bondi del sistema binario ($R_B$) y la separación entre los componentes ($a$), cuantificada mediante el parámetro de compactidad $\xi \equiv R_B/a$.

• Expectativa tradicional: Se asumía que para $\xi \gg 1$ (binaria profundamente incrustada en una esfera de Bondi compartida), el sistema acreta de manera cooperativa con una eficiencia $\eta \equiv \dot{M}_{binaria}/\dot{M}_{Bondi} \to 1$.
• Hipótesis del trabajo: Los autores proponen que esta expectativa falla para cualquier gas con índice adiabático $\gamma > 1$. La órbita binaria genera choques que producen entropía, creando turbulencia convectiva que suprime la acreción y enmascara la señal orbital.

2. Metodología

Los autores realizaron simulaciones hidrodinámicas tridimensionales utilizando el código de dinámica de gases tipo Godunov SIMBI.

• Parámetros de la simulación:
  • Compactidad: $\xi \in \{0.1, 1, 10\}$.
  • Índice adiabático: $\gamma \in \{1 \text{ (isotérmico)}, 4/3, 5/3\}$.
  • Relación de masas: $q = 1$ (masas iguales).
• Configuración numérica:
  • Se utilizó una malla telescópica (Fixed Mesh Refinement) para resolver desde la escala global del radio de Bondi ($r \sim \xi$) hasta los embudos de acreción individuales ($r \sim 1$).
  • Se resolvieron las ecuaciones de Euler con un solucionador Riemann HLLE y reconstrucción lineal/parabólica.
  • Se implementaron condiciones de frontera de relajación para mantener un suministro de masa infinito y evitar reflexiones artificiales.
  • Se definió una eficiencia de acreción normalizada para cancelar sesgos sistemáticos de los sumideros (sink particles).

3. Contribuciones Clave y Marco Teórico

El artículo deriva un criterio de estabilidad analítico desde primeros principios:

1. Escala Geométrica ($\xi_m$): Se define un umbral de compactidad $\xi_m = 4/(5-3\gamma)$ donde las superficies sónicas individuales de cada componente se tocan geométricamente.
   • Para $\gamma = 5/3$, $\xi_m \to \infty$, lo que implica que nunca existe una superficie sónica desconectada; el flujo es siempre transónico o sónico en la superficie de los acretores.
2. Mecanismo de Inestabilidad (Entropía):
   • En gas adiabático ($\gamma > 1$), el movimiento orbital actúa como un "propulsor" supersónico que genera choques. Estos choques producen un salto de entropía $\Delta s \propto (\gamma - 1)$.
   • La generación de entropía no es axisimétrica, creando gradientes que violan el criterio de estabilidad de Schwarzschild ($ds/dr < 0$), sembrando inestabilidad convectiva.
3. Aniquilación del Punto de Silla:
   • La superficie sónica actúa como la separatrix del punto de silla de Bondi. La acumulación de entropía en la región confinada (dentro de la superficie sónica) destruye este punto de silla tras un número de órbitas $N$:
     $$N \propto \frac{1}{\gamma - 1} \left( \sqrt{\frac{\xi}{\xi_m}} - 1 \right)$$
   • Condición de ruptura (Breakout): Si el tiempo de enfriamiento $t_{cool} \gg N T$ (donde $T$ es el periodo orbital), la entropía se acumula, la turbulencia destruye la estructura transónica y la acreción se suprime. Si $t_{cool} \ll NT$, el gas se comporta como isotérmico y la acreción es cooperativa.

4. Resultados Principales

A. Morfología del Flujo y Regímenes:

• Gas Isotérmico ($\gamma = 1$): La acreción es cooperativa y estable. Para $\xi \gg 1$, se forma un "puente" de alta densidad laminar que conecta los componentes, logrando $\eta \to 1$. La superficie sónica se mantiene estable y filtra la modulación orbital.
• Gas Adiabático ($\gamma = 4/3$): Por debajo de $\xi_m$, el flujo es organizado. Por encima de $\xi_m$, la entropía generada por los choques orbitales desencadena turbulencia convectiva. La estructura espiral coherente se destruye, reemplazada por fluctuaciones caóticas. La eficiencia cae a $\eta \approx 0.3$.
• Gas Monatómico ($\gamma = 5/3$): No existe superficie sónica desconectada. El flujo es inestable y turbulento en todos los $\xi$, incluso para separaciones grandes. La eficiencia es muy baja ($\eta \approx 0.1$).

B. Estructura Radial y Turbulencia:

• En los casos adiabáticos, la velocidad radial de infallo se suprime en 1-2 órdenes de magnitud respecto a la caída libre.
• La dispersión de velocidad turbulenta ($\delta v$) domina sobre el flujo coherente en el interior.
• El perfil de densidad en el interior se aplana a $\rho \propto r^{-1}$ (en lugar de $r^{-3/2}$ de Bondi), mientras que el perfil de temperatura (sonido) mantiene la estratificación de Bondi.
• El espectro de potencia de la entropía sigue una ley de potencia $k^{-3.5}$ (o $k^{-7/2}$), más pronunciada que la cascada de Kolmogorov ($k^{-5/3}$), indicando estructuras coherentes (frentes de choque) dominantes.

C. Variabilidad Temporal y Señales Observables:

• $\gamma = 5/3$: La señal orbital está completamente enterrada en ruido de banda ancha.
• $\gamma = 4/3$: La señal orbital (a $2\Omega$) sobrevive parcialmente, pero está amortiguada por el ruido turbulento.
• $\gamma = 1$: Se observa una periodicidad coherente, pero no a la frecuencia orbital fundamental ($2\Omega$) ni a la mitad ($\Omega$). Debido a la no linealidad del flujo, la potencia se transfiere a armónicos superiores ($4\Omega$ y $8\Omega$). Un observador inferiría incorrectamente que el periodo orbital es la mitad del real.

D. Eficiencia de Acreción ($\eta$):

• Isotérmico: $\eta$ aumenta monótonamente de 0.5 (acreción independiente) a 1.0 (cooperativa).
• Adiabático: $\eta$ cae drásticamente una vez que $\xi > \xi_m$. La turbulencia suprime la acreción por debajo de la tasa de acreción independiente.

5. Significado e Implicaciones Astrofísicas

1. Solución a la "Falta" de Binarias Periódicas: Explica por qué las búsquedas de binarias de agujeros negros supermasivos (SMBH) basadas en variabilidad periódica en curvas de luz de AGN han tenido resultados negativos (falsos positivos o ausencia de señales). En entornos adiabáticos (ineficientes radiativamente), la señal se borra o se desplaza a frecuencias inesperadas ($4\Omega$).
2. Problema del Último Parsec: En entornos calientes (post-fusión de galaxias), la acreción de gas para endurecer la binaria es mucho menos eficiente de lo que predice la teoría de Bondi laminar. Las estimaciones actuales de endurecimiento por gas podrían estar sobreestimadas por un factor de 3 a 10 ($1/\eta$), lo que agrava el problema del último parsec.
3. Nueva Ruta a la Acreción Dominada por Convección: Identifica un mecanismo puramente hidrodinámico (sin necesidad de viscosidad o campos magnéticos) para generar convección en flujos de acreción, impulsado por la energía orbital de la binaria.
4. Criterio de Diagnóstico: La relación entre el tiempo de enfriamiento y el tiempo de ruptura ($t_{cool}$ vs $NT$) determina si un sistema binario incrustado crecerá eficientemente (regímen isotérmico) o se ahogará en su propia turbulencia (regímen adiabático).

En conclusión, el trabajo demuestra que las binarias incrustadas en gas adiabático son auto-limitantes: la misma órbita que impulsa la acreción genera la entropía que la suprime, transformando el flujo ordenado en un reservorio turbulento con baja eficiencia de acreción y señales periódicas enmascaradas.