Self-Limited Accretion onto Embedded Binaries in a Uniform Medium

1. Problématique et Contexte

L'article s'intéresse au problème de l'accrétion de gaz uniforme au repos sur une binaire de masses égales, souvent appelé « problème de Bondi binaire ». La question centrale est de savoir si une binaire enfouie dans un milieu gazeux accrète de manière coopérative (comme un objet unique de masse totale $M$) ou de manière indépendante (comme deux objets séparés).

La réponse dépend de deux paramètres adimensionnels :

• La compacité $\xi \equiv R_B / a$, où $R_B$ est le rayon de Bondi de la binaire et $a$ la séparation des composantes.
• L'indice adiabatique $\gamma$ du gaz.

L'hypothèse traditionnelle suggère que pour $\xi \gg 1$ (binaire profondément enfouie), l'accrétion devrait être coopérative avec une efficacité $\eta \equiv \dot{M}_{\text{binaire}} / \dot{M}_{\text{Bondi}} \to 1$. Cependant, les auteurs montrent que cette attente échoue pour tout gaz avec $\gamma > 1$.

2. Méthodologie

Les auteurs ont combiné une analyse théorique rigoureuse et des simulations numériques tridimensionnelles.

• Simulations Numériques :

  • Utilisation du code hydrodynamique Godunov SIMBI.
  • Résolution des équations d'Euler et de continuité avec une équation d'état polytropique ($P \propto \rho^\gamma$).
  • Étude de trois régimes thermodynamiques : isotherme ($\gamma = 1$), relativiste ($\gamma = 4/3$) et monoatomique ($\gamma = 5/3$).
  • Exploration d'une gamme de compacités $\xi \in \{0.1, 1, 10\}$.
  • Utilisation d'une grille télescopique (Fixed Mesh Refinement) pour résoudre les échelles allant de la séparation binaire ($r \sim a$) au rayon de Bondi ($r \sim R_B$).
  • Conditions aux limites relaxantes pour maintenir un flux de masse infini sans réflexion artificielle.

• Analyse Théorique :

  • Dérivation d'un critère de stabilité basé sur la génération d'entropie par les ondes de choc.
  • Identification du point de selle de Bondi comme séparatrice du flux transsonique.
  • Établissement d'une relation entre la génération d'entropie, la convection et la destruction de la structure d'accrétion stationnaire.

3. Contributions Clés et Résultats

A. Le Mécanisme d'Auto-Limitation (Instabilité Entropique)

Pour les gaz adiabatiques ($\gamma > 1$), l'orbite binaire agit comme un « propulseur » (propeller) supersonique par rapport au gaz local, générant des ondes de choc.

• Génération d'entropie : Ces chocs créent de l'entropie de manière irréversible. Contrairement au cas isotherme où l'entropie est radiée instantanément, le gaz adiabatique conserve cette entropie.
• Instabilité convective : La distribution non-axisymétrique de l'entropie crée des gradients de densité instables (violation du critère de Schwarzschild), déclenchant une turbulence convective à grande échelle.
• Destruction du point de selle : Cette turbulence accumule de l'entropie jusqu'à ce qu'elle « annihile » le point de selle (saddle point) de Bondi, qui est la condition nécessaire pour un flux d'accrétion transsonique stable. Une fois ce point détruit, le flux d'accrétion global est fortement supprimé.

B. Résultats selon l'Indice Adiabatique ($\gamma$)

1. Cas Isotherme ($\gamma = 1$) :

   • Les chocs ne génèrent pas d'entropie (ou celle-ci est immédiatement dissipée).
   • Le flux reste laminaire et stationnaire.
   • Pour $\xi \gg 1$, une « passerelle » de haute densité se forme entre les deux composantes, permettant une accrétion coopérative efficace ($\eta \to 1$).

2. Cas Adiabatique ($\gamma = 4/3$ et $5/3$) :

   • Seuil de compacité ($\xi_m$) : L'instabilité se déclenche lorsque la surface sonique englobe la binaire. Le seuil théorique est $\xi_m = 4/(5-3\gamma)$.
     • Pour $\gamma = 4/3$, $\xi_m = 4$.
     • Pour $\gamma = 5/3$, $\xi_m \to \infty$ (le flux est toujours instable car aucune surface sonique détachée n'existe).
   • Suppression de l'accrétion : Au-delà du seuil, la turbulence détruit la structure cohérente.
     • Pour $\gamma = 4/3$, l'efficacité chute à $\eta \approx 0.3$.
     • Pour $\gamma = 5/3$, l'efficacité chute à $\eta \approx 0.1$.
   • Profil de densité : La densité interne ne suit plus la loi de chute libre ($\rho \propto r^{-3/2}$) mais s'aplatit vers $\rho \propto r^{-1}$, caractéristique des écoulements dominés par la convection.

C. Variabilité Temporelle et Signaux Observables

L'étude de la variabilité temporelle du taux d'accrétion ($\dot{M}$) révèle des signatures distinctes :

• $\gamma = 5/3$ : Le signal orbital est complètement noyé dans un bruit large bande (turbulence chaotique). Aucune périodicité n'est détectable.
• $\gamma = 4/3$ : Une résonance résiduelle à $2\Omega$ (deux fois la fréquence orbitale) survit mais est partiellement enterrée dans le bruit.
• $\gamma = 1$ : Le signal est cohérent mais non linéaire. Au lieu de la fréquence fondamentale $2\Omega$ attendue pour une binaire égale, l'accrétion montre des pics dominants à $4\Omega$ et $8\Omega$. Un observateur déduirait donc une période orbitale deux fois plus courte que la réalité.

D. Critère de Stabilité et Échelle de Temps

Les auteurs dérivent un critère de stabilité fondamental basé sur le temps de refroidissement $t_{\text{cool}}$ par rapport au temps orbital $T$ et au nombre d'orbites nécessaires pour détruire le point de selle $N$ :

• Si $t_{\text{cool}} \ll N T$ : Le gaz se refroidit avant que l'entropie ne s'accumule. Le système reste coopératif (régime isotherme effectif).
• Si $t_{\text{cool}} \gg N T$ : L'entropie s'accumule, déclenchant la turbulence et l'auto-limitation (régime adiabatique).

4. Signification et Implications Astrophysiques

1. Le Problème du Dernier Parsec :
   Les modèles de synthèse de population pour les trous noirs supermassifs (SMBH) binaires utilisent souvent des taux d'accrétion de type Bondi laminaire pour estimer le durcissement orbital par le gaz. Ce papier montre que dans les environnements chauds et radiativement inefficaces (typiques des noyaux galactiques), l'accrétion est réduite d'un facteur 3 à 10 ($\eta \ll 1$). Cela pourrait aggraver le « problème du dernier parsec » en ralentissant la fusion des binaires.

2. Recherche de Binaires Périodiques :
   L'absence de détection de binaires SMBH périodiques dans les courbes de lumière des AGN (Active Galactic Nuclei) peut s'expliquer par ce mécanisme :

   • Soit le signal est totalement effacé par la turbulence (cas adiabatique).
   • Soit le signal est décalé vers des harmoniques supérieurs ($4\Omega$), rendant la détection difficile si l'on cherche uniquement la fréquence fondamentale.

3. Nouvelle Voie vers la Convection Dominée :
   L'article identifie une troisième voie menant à des écoulements dominés par la convection (CDAF - Convective-Dominated Accretion Flows), distincte de la dissipation visqueuse ou magnétique. Ici, l'énergie provient du mouvement orbital de la binaire elle-même, générant de l'entropie via des chocs hydrodynamiques purs.

Conclusion

L'article démontre que l'accrétion sur des binaires enfouies dans un gaz adiabatique est intrinsèquement auto-limitée. L'orbite qui devrait alimenter l'accrétion génère en réalité l'entropie qui la supprime. La distinction entre un système qui accrète efficacement et un système qui étouffe son propre carburant dépend d'une seule comparaison : le temps de refroidissement du gaz par rapport au temps nécessaire pour que la turbulence détruise la structure d'accrétion. Cela remet en question les estimations actuelles de la croissance des trous noirs binaires et de la détection de leurs signatures périodiques.