Radiation-hydrodynamics of star-disc collisions: From system parameters to outflows and lightcurves

Résumé technique : Hydrodynamique de la radiation des collisions étoile-disque : des paramètres du système aux flux sortants et aux courbes de lumière

Problème et contexte
Les éruptions quasi-périodiques (QPE) sont des transitoires nucléaires associés à des trous noirs supermassifs (SMBH), caractérisés par des éruptions de rayons X mous répétitives. Bien que plusieurs modèles existent, le scénario de collision étoile-disque — où une étoile sur une orbite étroite, excentrique et inclinée intersecte de manière répétée un disque d'accrétion — est un candidat de premier plan. Ce modèle explique naturellement les temps de récurrence, les intervalles alternés « long–court » et les motifs d'éruption « fort–faible » observés dans des sources comme GSN 069. Cependant, relier les paramètres physiques locaux de ces collisions aux propriétés observables des éruptions (luminosité, durée, asymétrie) nécessite de comprendre l'évolution couplée des chocs, de la géométrie des flux sortants et du transport radiatif dans un milieu en expansion et optiquement épais. Des études précédentes ont exploré certaines parties de ce problème en utilisant des approches 1D, 2D ou non radiatives, mais une étude systématique en 3D de la radiation-hydrodynamique reliant les paramètres du système aux courbes de lumière faisait défaut.

Méthodologie
Les auteurs ont réalisé une série de simulations de radiation-hydrodynamique tridimensionnelles locales en utilisant le code de smoothed-particle hydrodynamics (SPH) Phantom. Les simulations ont été effectuées dans un référentiel cartésien local centré sur le plan médian du disque, négligeant la dynamique orbitale globale, la relativité générale et l'auto-gravité du gaz, car les échelles de temps des collisions sont bien plus courtes que les échelles de temps orbitales.

L'étude a utilisé un modèle fiduciaire basé sur une étoile de la séquence principale ($R_\star \approx R_\odot$) entrant en collision avec un disque $\alpha$ optiquement épais et dominé par la pression de radiation autour d'un SMBH de $10^6 M_\odot$. L'étoile a été modélisée comme une sphère rigide se déplaçant à $v_\star = 0,1c$. Les auteurs ont fait varier systématiquement cinq paramètres clés par rapport au modèle fiduciaire, en ne faisant varier qu'un paramètre à la fois :

1. Vitesse stellaire ($v_\star$) : allant de $0,5$à$1,5$ fois la valeur fiduciaire.
2. Densité de surface du disque ($\Sigma_d$) : allant de $1$à$100$ fois la valeur fiduciaire, couvrant des disques fins standards jusqu'à des disques denses post-événement de disruption tidale (TDE).
3. Rayon stellaire ($R_\star$) : allant de $0,5$à$4$ fois la valeur fiduciaire, faisant varier efficacement le rapport entre la hauteur d'échelle du disque et le rayon stellaire ($H/R_\star$).
4. Profil de densité verticale du disque : modélisé comme une distribution gaussienne $\rho(z) \propto e^{-z^2/2\sigma^2}$, faisant varier la concentration de masse vers le plan médian ($\sigma/R_\odot$).
5. Angle de collision local ($i$) : l'angle entre le vecteur de vitesse relative et le plan médian du disque, allant de $9^\circ$ (perpendiculaire) à $30^\circ$ (oblique).

Le transport de la radiation a été traité par l'approximation de la diffusion à limite de flux avec une opacité de diffusion électronique. La luminosité bolométrique émergente a été calculée en intégrant les flux de diffusion et d'advection à travers une photosphère déterminée dynamiquement.

Contributions clés et résultats
L'étude fournit des mises à l'échelle empiriques pour la luminosité de crête ($L_{\text{peak}}$) et la durée de l'éruption ($\Delta t$) en fonction des paramètres du système. Les principales conclusions sont les suivantes :

• Dynamique de collision : L'étoile génère un choc de proue, créant deux flux sortants (avant et arrière) de part et d'autre du disque. Le flux avant est généralement plus massif et énergétique. La morphologie qualitative du choc de proue et des flux sortants reste largement inchangée à travers les variations de $v_\star$ et $\Sigma_d$, bien que la normalisation thermodynamique change.
• Vitesse stellaire ($v_\star$) : Les étoiles plus rapides produisent des éruptions plus brillantes ($L_{\text{peak}} \propto v_\star^2$) en raison de taux de chauffage par choc plus élevés. Elles produisent également des éruptions plus longues ($\Delta t \propto v_\star^{1,4-2,2}$) car plus d'énergie est injectée, bien que la luminosité décroisse plus rapidement en raison de pertes adiabatiques plus fortes.
• Densité de surface du disque ($\Sigma_d$) : Les disques plus denses augmentent la masse et la quantité de mouvement des flux sortants mais ont un effet faible sur la luminosité de crête ($L_{\text{peak}} \propto \Sigma_d^{-0,1}$ pour l'avant, $\propto \Sigma_d^{0,2}$ pour l'arrière). Ceci est dû à une compétition entre l'augmentation du chauffage par choc et des pertes adiabatiques plus fortes dans les éjectas optiquement plus épais. Cependant, les disques plus denses prolongent considérablement la durée de l'éruption ($\Delta t \propto \Sigma_d^{0,5-0,6}$) en raison de temps de diffusion des photons plus longs.
• Rayon stellaire ($R_\star$) : Les étoiles plus grandes interagissent avec plus de gaz, produisant des éruptions plus brillantes ($L_{\text{peak}} \propto R_\star^{1,5}$ avant, $\propto R_\star^{0,7}$ arrière) et plus longues ($\Delta t \propto R_\star^{1,5}$). L'augmentation de $R_\star$ accentue l'asymétrie de quantité de mouvement entre les flux avant et arrière.
• Profil de densité verticale ($\sigma$) : Les disques plus concentrés au centre (plus petit $\sigma$) produisent des éruptions plus brillantes et plus courtes. Les couches centrales denses absorbent plus d'énergie de choc, tandis que les couches externes de faible densité permettent à la radiation de s'échapper plus efficacement avec des pertes adiabatiques réduites.
• Angle de collision ($i$) : Les passages plus obliques (plus petit $i$) réduisent l'asymétrie de quantité de mouvement et de luminosité entre les deux flux. Le flux arrière devient nettement plus brillant ($L_{\text{peak}} \propto (\sin i)^{-0,8}$) et les éruptions deviennent plus longues ($\Delta t \propto (\sin i)^{-0,9 \text{ à } -1}$) en raison du temps d'interaction plus long et d'une redistribution plus symétrique du gaz choqué.

Application à GSN 069
Les auteurs ont appliqué leurs mises à l'échelle empiriques à la source bien caractérisée GSN 069. En échantillonnant l'espace des paramètres pour faire correspondre les luminosités de crête observées ($\sim 3 \times 10^{42}$ erg s$^{-1}$), les durées ($\sim 1$ h) et le rapport d'amplitude fort-faible, ils ont identifié deux régions viables :

1. Solution préférée : Une orbite rétrograde collisionnant avec un disque post-TDE dense et verticalement concentré ($\Sigma_d \sim 45 \times$ la valeur fiduciaire, $\sigma/R_\odot \sim 0,3$). Ce modèle reproduit bien les propriétés observées des éruptions.
2. Solution alternative : Une orbite rétrograde avec une densité de disque fin standard mais une géométrie de collision très rasante ($i \sim 6^\circ$). Cette solution est moins précise et nécessite des paramètres extrêmes.

Les auteurs notent que la préférence pour un disque dense et verticalement concentré est cohérente avec les preuves observationnelles d'une activité de type TDE dans GSN 069.

Signification et limites
L'article affirme fournir la première étude systématique de radiation-hydrodynamique 3D reliant les paramètres de collision locale aux courbes de lumière des QPE. Les mises à l'échelle empiriques dérivées servent de pont entre les simulations locales et les modèles globaux de QPE, permettant de prédire les propriétés des éruptions basées sur les paramètres du système. Les auteurs soulignent que leurs résultats suggèrent que des combinaisons spécifiques de paramètres peuvent reproduire les amplitudes, les durées et les cycles de service caractéristiques des QPE.

L'étude reconnaît modestement ses limites : les modèles supposent des mises à l'échelle à un paramètre séparables, négligent la rotation et le cisaillement du disque, traitent l'étoile comme une sphère rigide, et ne tiennent pas compte des effets de l'angle de vue ou de l'évolution spectrale (par exemple, la Comptonisation). Par conséquent, les résultats sont présentés comme des régions physiquement motivées de l'espace des paramètres plutôt que comme des modèles uniques et définitifs. Les auteurs suggèrent que les travaux futurs devraient incorporer des structures stellaires réalistes, une dynamique de disque, et des luminosités dépendantes de l'angle de vue pour quantifier pleinement la capacité du modèle à reproduire la diversité des sources QPE observées.