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

Technische Zusammenfassung: Strahlungshydrodynamik von Stern–Disk-Kollisionen: Von Systemparametern zu Ausflüssen und Lichtkurven

Problem und Kontext
Quasi-periodische Eruptionen (QPEs) sind nukleare Transienten, die mit supermassereichen Schwarzen Löchern (SMBHs) assoziiert werden und durch sich wiederholende weiche Röntgenausbrüche charakterisiert sind. Während mehrere Modelle existieren, ist das Szenario der Stern–Disk-Kollision – bei dem ein Stern auf einer engen, exzentrischen und geneigten Umlaufbahn wiederholt eine Akkretionsscheibe durchquert – ein führender Kandidat. Dieses Modell erklärt auf natürliche Weise die Wiederholungszeiten, die abwechselnden „lang–kurz“-Intervalle und die „stark–schwach“-Ausbruchsmuster, wie sie in Quellen wie GSN 069 beobachtet werden. Die Verknüpfung der lokalen physikalischen Parameter dieser Kollisionen mit den beobachtbaren Ausbruchseigenschaften (Leuchtkraft, Dauer, Asymmetrie) erfordert jedoch ein Verständnis der gekoppelten Entwicklung von Schocks, Ausflussgeometrie und Strahlungstransport in einem expandierenden, optisch dichten Medium. Vorherige Studien haben Teile dieses Problems mittels 1D-, 2D- oder nicht-strahlungsbasierten Ansätzen untersucht, aber eine dedizierte systematische 3D-Strahlungshydrodynamik-Studie, welche die Systemparameter mit Lichtkurven verknüpft, fehlte bislang.

Methodik
Die Autoren führten eine Reihe von dreidimensionalen lokalen Strahlungs-Hydrodynamik-Simulationen unter Verwendung des Smoothed-Particle-Hydrodynamics-Codes (SPH) Phantom durch. Die Simulationen wurden in einem lokalen kartesischen Bezugssystem zentriert auf der Mittelscheibe der Disk durchgeführt, wobei globale Orbitaldynamik, allgemeine Relativitätstheorie und die Selbstgravitation des Gases vernachlässigt wurden, da die Kollisions-Zeitskalen wesentlich kürzer sind als die orbitalen Zeitskalen.

Die Studie nutzte ein Fiduzialmodell basierend auf einem Hauptreihenstern ($R_\star \approx R_\odot$), der mit einer optisch dichten, strahlungsdruckdominierten $\alpha$-Disk um ein $10^6 M_\odot$ SMBH kollidiert. Der Stern wurde als starrer Körper modelliert, der sich mit $v_\star = 0.1c$ bewegt. Die Autoren variierten systematisch fünf Schlüsselparameter relativ zum Fiduzialmodell, indem sie jeweils einen Parameter gleichzeitig veränderten:

1. Stellare Geschwindigkeit ($v_\star$): Bereich von $0.5$ bis $1.5$ mal dem Fiduzialwert.
2. Disk-Oberflächendichte ($\Sigma_d$): Bereich von $1$ bis $100$ mal dem Fiduzialwert, was von Standard-Dünnscheiben bis hin zu dichten Post-TDE-Disk-Modellen reicht.
3. Stellarer Radius ($R_\star$): Bereich von $0.5$ bis $4$ mal dem Fiduzialwert, was effektiv das Verhältnis der Skalenhöhe der Disk zum stellaren Radius ($H/R_\star$) variiert.
4. Vertikales Disk-Dichteprofil: Modelliert als Gauß-Verteilung $\rho(z) \propto e^{-z^2/2\sigma^2}$, wobei die Massenkonzentration zur Mittelscheibe hin variiert wird ($\sigma/R_\odot$).
5. Lokaler Kollisionswinkel ($i$): Der Winkel zwischen dem relativen Geschwindigkeitsvektor und der Disk-Mittelscheibe, reichend von $90^\circ$ (senkrecht) bis $30^\circ$ (schräg).

Der Strahlungstransport wurde mittels der Flux-Limited-Diffusion-Approximation mit Elektronenstreu-Opazität behandelt. Die resultierende bolometrische Leuchtkraft wurde durch Integration der Diffusions- und Advektionsflüsse durch eine dynamisch bestimmte Photosphäre berechnet.

Wesentliche Beiträge und Ergebnisse
Die Studie liefert empirische Skalierungen für die Spitzenleuchtkraft ($L_{\text{peak}}$) und die Ausbruchsdauer ($\Delta t$) als Funktionen der Systemparameter. Die Hauptergebnisse sind:

• Kollisionsdynamik: Der Stern treibt einen Bugschock voran, der zwei Ausflüsse (vorwärts und rückwärts) auf gegenüberliegenden Seiten der Disk erzeugt. Der Vorwärtsausfluss ist im Allgemeinen massereicher und energetischer. Die qualitative Morphologie des Bugschocks und der Ausflüsse bleibt über Variationen in $v_\star$ und $\Sigma_d$ weitgehend unverändert, wenngleich sich die thermodynamische Normierung ändert.
• Stellare Geschwindigkeit ($v_\star$): Schnellere Sterne erzeugen hellere Ausbrüche ($L_{\text{peak}} \propto v_\star^2$) aufgrund höherer Schockheizungsraten. Sie erzeugen auch längere Ausbrüche ($\Delta t \propto v_\star^{1.4-2.2}$), da mehr Energie injiziert wird, obwohl die Leuchtkraft aufgrund stärkerer adiabatischer Verluste schneller abfällt.
• Disk-Oberflächendichte ($\Sigma_d$): Dichtere Disks erhöhen die Masse und den Impuls der Ausflüsse, haben aber einen schwachen Effekt auf die Spitzenleuchtkraft ($L_{\text{peak}} \propto \Sigma_d^{-0.1}$ für Vorwärts, $\propto \Sigma_d^{0.2}$ für Rückwärts). Dies liegt am Wettbewerb zwischen erhöhter Schockheizung und stärkeren adiabatischen Verlusten in optisch dickeren Ejekta. Dichtere Disks verlängern jedoch die Ausbruchsdauer signifikant ($\Delta t \propto \Sigma_d^{0.5-0.6}$) aufgrund längerer Photondiffusionszeiten.
• Stellarer Radius ($R_\star$): Größere Sterne interagieren mit mehr Gas, was zu helleren ($L_{\text{peak}} \propto R_\star^{1.5}$ vorwärts, $\propto R_\star^{0.7}$ rückwärts) und längeren ($\Delta t \propto R_\star^{1.5}$) Ausbrüchen führt. Eine Erhöhung von $R_\star$ verstärkt die Impulsasymmetrie zwischen Vorwärts- und Rückwärtsausfluss.
• Vertikales Dichteprofil ($\sigma$): Zentral konzentriertere Disks (kleineres $\sigma$) liefern hellere und kürzere Ausbrüche. Die dichten zentralen Schichten absorbieren mehr Schockenergie, während die ungeordneten äußeren Schichten die Strahlung effizienter mit reduzierten adiabatischen Verlusten entweichen lassen.
• Kollisionswinkel ($i$): Obliquere Durchgänge (kleineres $i$) reduzieren die Impuls- und Leuchtkraftasymmetrie zwischen den beiden Ausflüssen. Der Rückwärtsausfluss wird signifikant heller ($L_{\text{peak}} \propto (\sin i)^{-0.8}$) und die Ausbrüche werden länger ($\Delta t \propto (\sin i)^{-0.9 \text{ bis } -1}$) aufgrund der längeren Interaktionszeit und der symmetrischeren Umverteilung des geschockten Gases.

Anwendung auf GSN 069
Die Autoren wendeten ihre empirischen Skalierungen auf die gut charakterisierte Quelle GSN 069 an. Durch Abtastung des Parameterraums, um die beobachteten Spitzenleuchtkräfte ($\sim 3 \times 10^{42}$ erg s$^{-1}$), Dauern ($\sim 1$ h) und das Stark-Schwach-Amplitudenverhältnis abzugleichen, identifizierten sie zwei lebensfähige Regionen:

1. Bevorzugte Lösung: Eine retrograde Umlaufbahn, die mit einer dichten, vertikal konzentrierten Post-TDE-Disk kollidiert ($\Sigma_d \sim 45 \times$ Fiduzial, $\sigma/R_\odot \sim 0.3$). Dieses Modell reproduziert die beobachteten Ausbruchseigenschaften gut.
2. Alternative Lösung: Eine retrograde Umlaufbahn mit einer Standard-Dünnscheiben-Dichte, aber einer sehr flachen Kollisionsgeometrie ($i \sim 6^\circ$). Diese Lösung ist weniger genau und erfordert extreme Parameter.

Die Autoren stellen fest, dass die Bevorzugung einer dichten, vertikal konzentrierten Disk konsistent mit den Beobachtungsdaten zu TDE-ähnlicher Aktivität in GSN 069 ist.

Bedeutung und Einschränkungen
Das Paper beansprucht, die erste systematische 3D-Strahlungs-Hydrodynamik-Studie zu sein, die lokale Kollisionsparameter mit QPE-Lichtkurven verknüpft. Die abgeleiteten empirischen Skalierungen dienen als Brücke zwischen lokalen Simulationen und globalen QPE-Modellen, was die Vorhersage von Ausbruchseigenschaften basierend auf Systemparametern ermöglicht. Die Autoren betonen, dass ihre Ergebnisse darauf hindeuten, dass spezifische Kombinationen von Parametern die charakteristischen Amplituden, Dauern und Duty Cycles von QPEs reproduzieren können.

Die Studie räumt bescheiden Einschränkungen ein: Die Modelle nehmen separierbare Ein-Parameter-Skalierungen an, vernachlässigen die Rotations- und Scherungseffekte der Disk, behandeln den Stern als starren Körper und berücksichtigen keine Sichtwinkel-Effekte oder die spektrale Evolution (z. B. Comptonisierung). Folglich werden die Ergebnisse als physikalisch motivierte Parameterbereiche präsentiert und nicht als eindeutige, definitiv beste Anpassungsmodelle. Die Autoren schlagen vor, dass zukünftige Arbeiten realistische Sternstrukturen, Diskdynamiken und Sichtwinkel-abhängige Leuchtkräfte einbeziehen sollten, um die Fähigkeit des Modells, die Vielfalt beobachteter QPE-Quellen zu reproduzieren, vollständig zu quantifizieren.