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

Technische Samenvatting: Radiatie-hydrodynamica van ster–schijf-botsingen: Van systeemparameters tot outflows en lichtcurves

Probleem en Context
Quasi-periodieke uitbarstingen (QPE's) zijn nucleaire transiënten geassocieerd met supermassieve zwarte gaten (SMBH's), gekenmerkt door herhalende zachte röntgenflares. Hoewel er verschillende modellen bestaan, is het ster–schijf-botsingsscenario—waarbij een ster op een nauwe, excentrische en geïnclineerde baan herhaaldelijk een accretieschijf snijdt—een leidende kandidaat. Dit model verklaart op natuurlijke wijze de herhalingstijden, de afwisselende "lang–kort" intervallen en de "sterk–zwak" flare-patronen waargenomen in bronnen zoals GSN 069. Echter, het verbinden van de lokale fysieke parameters van deze botsingen aan de observeerbare flare-eigenschappen (luminositeit, duur, asymmetrie) vereist een begrip van de gekoppelde evolutie van schokken, outflow-geometrie en radiatieve transport in een expanderend, optisch dik medium. Eerdere studies hebben delen van dit probleem onderzocht met behulp van 1D, 2D of niet-radiatieve benaderingen, maar een toegewijde systematische 3D radiatie-hydrodynamische studie die systeemparameters verbindt met lichtcurves ontbrak.

Methodologie
De auteurs voerden een reeks driedimensionale lokale radiatie-hydrodynamische simulaties uit met behulp van de smoothed-particle hydrodynamics (SPH) code Phantom. De simulaties werden uitgevoerd in een lokaal Cartesiaans kader gecentreerd op het midvlak van de schijf, waarbij de globale orbitale dynamica, algemene relativiteitstheorie en de zelfgravitatie van het gas werden genegeerd, aangezien de botsingstijdsschalen veel korter zijn dan de orbitale tijdsschalen.

De studie maakte gebruik van een fiducieel model gebaseerd op een hoofdreeksster ($R_\star \approx R_\odot$) die botst met een optisch dikke, door stralingsdruk gedomineerde $\alpha$-schijf rond een $10^6 M_\odot$ SMBH. De ster werd gemodelleerd als een rigide sfeer die beweegt met $v_\star = 0.1c$. De auteurs varieerden systematisch vijf sleutelparameters ten opzichte van het fiduciaal model, waarbij ze telkens één parameter tegelijk varieerden:

1. Stellaire snelheid ($v_\star$): Variërend van $0.5$ tot $1.5$ keer de fiduciaal waarde.
2. Schijf-oppervlakte-dichtheid ($\Sigma_d$): Variërend van $1$ tot $100$ keer de fiduciaal waarde, reikend van standaard dunne schijven tot dichte post-tidal disruption event (TDE) schijven.
3. Stellaire radius ($R_\star$): Variërend van $0.5$ tot $4$ keer de fiduciaal waarde, wat effectief de verhouding tussen de schijf-schaalhoogte en de stellaire radius ($H/R_\star$) varieert.
4. Verticale schijf-dichtheidsprofiel: Gemodelleerd als een Gaussische distributie $\rho(z) \propto e^{-z^2/2\sigma^2}$, waarbij de concentratie van massa naar het midvlak wordt gevarieerd ($\sigma/R_\odot$).
5. Lokale botsingshoek ($i$): De hoek tussen de relatieve snelheidsvector en het schijf-midvlak, variërend van $90^\circ$ (loodrecht) tot $30^\circ$ (obliqu).

Radiatieve transport werd behandeld met de flux-beperkte diffusie-benadering met elektron-verstrooiingsopaciteit. De resulterende bolometrische luminositeit werd berekend door integratie van diffusie- en advectieve fluxen door een dynamisch bepaalde fotosfeer.

Belangrijkste Bijdragen en Resultaten
De studie biedt empirische schalingen voor de piek-luminositeit ($L_{\text{peak}}$) en de flare-duur ($\Delta t$) als functies van de systeemparameters. De belangrijkste bevindingen zijn:

• Botsingsdynamica: De ster drijft een bow shock aan, wat twee outflows (forward en backward) creëert aan weerszijden van de schijf. De forward outflow is over het algemeen massiever en energieker. De kwalitatieve morfologie van de bow shock en de outflows blijft grotendeels ongewijzigd bij variaties in $v_\star$ en $\Sigma_d$, hoewel de thermodynamische normalisatie verandert.
• Stellaire snelheid ($v_\star$): Snellere sterren produceren helderdere flares ($L_{\text{peak}} \propto v_\star^2$) vanwege hogere schokverwarmingssnelheden. Ze produceren ook langere flares ($\Delta t \propto v_\star^{1.4-2.2}$) omdat er meer energie wordt geïnjecteerd, hoewel de luminositeit sneller afneemt door sterkere adiabatische verliezen.
• Schijf-oppervlakte-dichtheid ($\Sigma_d$): D seeminglyere schijven verhogen de massa en impuls van de outflows, maar hebben een zwak effect op de piek-luminositeit ($L_{\text{peak}} \propto \Sigma_d^{-0.1}$ voor forward, $\propto \Sigma_d^{0.2}$ voor backward). Dit komt door een competitie tussen verhoogde schokverwarming en sterkere adiabatische verliezen in optisch dikkere ejecta. Echter, dichtere schijven verlengen de flare-duur aanzienlijk ($\Delta t \propto \Sigma_d^{0.5-0.6}$) door langere foton-diffusietijden.
• Stellaire radius ($R_\star$): Grotere sterren interageren met meer gas, wat leidt tot helderdere ($L_{\text{peak}} \propto R_\star^{1.5}$ forward, $\propto R_\star^{0.7}$ backward) en langere ($\Delta t \propto R_\star^{1.5}$) flares. Het vergroten van $R_\star$ versterkt de momentum-asymmetrie tussen de forward en backward outflows.
• Vertical dichtheidsprofiel ($\sigma$): Meer centraal geconcentreerde schijven (kleinere $\sigma$) leveren helderdere en kortere flares op. De dichte centrale lagen absorberen meer schokenergie, terwijl de laag-densiteit buitenste lagen de straling efficiënter laten ontsnappen met verminderde adiabatische verliezen.
• Botsingshoek ($i$): Meer obliqu crossing (kleinere $i$) vermindert de momentum- en luminositeitsasymmetrie tussen de twee outflows. De backward outflow wordt aanzienlijk helderder ($L_{\text{peak}} \propto (\sin i)^{-0.8}$) en de flares worden langer ($\Delta t \propto (\sin i)^{-0.9 \text{ tot } -1}$) door de langere interactietijd en een meer symmetrische herverdeling van het geschokte gas.

Toepassing op GSN 069
De auteurs pasten hun empirische schalingen toe op de goed gekarakteriseerde bron GSN 069. Door de parameterruimte te bemonsteren om de geobserveerde piek-luminositeiten ($\sim 3 \times 10^{42}$ erg s$^{-1}$), duur ($\sim 1$ uur) en de sterk-zwak amplitude ratio te matchen, identificeerden zij twee levensvatbare regio's:

1. Voorkeursoplossing: Een retrograde baan die botst met een dichte, verticaal geconcentreerde post-TDE schijf ($\Sigma_d \sim 45 \times$ fiduciaal, $\sigma/R_\odot \sim 0.3$). Dit model reproduceert de geobserveerde flare-eigenschappen goed.
2. Alternatieve oplossing: Een retrograde baan met een standaard dunne schijf-dichtheid maar een zeer grazende botsingsgeometrie ($i \sim 6^\circ$). Deze oplossing is minder accuraat en vereist extreme parameters.

De auteurs merken op dat de voorkeur voor een dichte, verticaal geconcentreerde schijf consistent is met observationele bewijzen van TDE-achtige activiteit in GSN 069.

Betekenis en Beperkingen
Het artikel claimt de eerste systematische 3D radiatie-hydrodynamische studie te zijn die lokale botsingsparameters verbindt met QPE-lichtcurves. De afgeleide empirische schalingen dienen als een brug tussen lokale simulaties en globale QPE-modellen, waardoor de voorspelling van flare-eigenschappen op basis van systeemparameters mogelijk wordt. De auteurs benadrukken dat hun bevindingen suggereren dat specifieke combinaties van parameters de karakteristieke amplitudes, duur en duty cycles van QPE's kunnen reproduceren.

De studie erkent bescheiden beperkingen: de modellen gaan uit van scheidbare één-parameter schalingen, negeren schijfrotatie en afschuiving (shear), behandelen de ster als een rigide sfeer, en houden geen rekening met kijkhoek-effecten of spectrale evolutie (bijv. Comptonisatie). Bij gevolg is de resultaten gepresenteerd als fysiek gemotiveerde regio's van de parameterruimte in plaats van unieke, definitieve best-fit modellen. De auteurs suggereren dat toekomstig werk realistische stellaire structuren, schijfdynamica en kijkhoek-afhankelijke luminositeiten moet bevatten om het vermogen van het model om de diversiteit van geobserveerde QPE-bronnen te reproduceren volledig te kwantificeren.