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

Sintesi Tecnica: Idrodinamica della radiazione delle collisioni tra stella e disco: dai parametri del sistema agli outflow e alle curve di luce

Problema e Contesto
Le esplosioni quasi periodiche (QPE) sono transienti nucleari associati a buchi neri supermassicci (SMBH), caratterizzati da ripetuti flare di raggi X "soft". Sebbene esistano diversi modelli, lo scenario della collisione tra stella e disco — in cui una stella su un'orbita stretta, eccentrica e inclinata interseca ripetutamente un disco di accrezione — è un candidato principale. Questo modello spiega naturalmente i tempi di ricorrenza, gli intervalli alternati "lungo-corto" e i pattern di flare "forte-debole" osservati in sorgenti come GSN 069. Tuttavia, collegare i parametri fisici locali di queste collisioni alle proprietà osservabili dei flare (luminosità, durata, asimmetria) richiede la comprensione dell'evoluzione accoppiata di shock, geometria degli outflow e trasporto radiativo in un mezzo in espansione e otticamente spesso. Studi precedenti hanno esplorato parti di questo problema utilizzando approcci 1D, 2D o non radiativi, ma mancava uno studio sistematico di idrodinamica della radiazione in 3D dedicato a collegare i parametri del sistema alle curve di luce.

Metodologia
Gli autori hanno eseguito una serie di simulazioni di idrodinamica della radiazione 3D locali utilizzando il codice di smoothed-particle hydrodynamics (SPH) Phantom. Le simulazioni sono state condotte in un sistema di riferimento cartesiano locale centrato sul piano medio del disco, trascurando la dinamica orbitale globale, la relatività generale e l'autogravità del gas, poiché le scale temporali delle collisioni sono molto più brevi delle scale temporali orbitali.

Lo studio ha utilizzato un modello fiduciale basato su una stella di sequenza principale ($R_\star \approx R_\odot$) che collide con un disco $\alpha$ otticamente spesso e dominato dalla pressione di radiazione attorno a un SMBH di $10^6 M_\odot$. La stella è stata modellata come una sfera rigida che si muove a $v_\star = 0.1c$. Gli autori hanno variato sistematicamente cinque parametri chiave rispetto al modello fiduciale, variandone uno alla volta:

1. Velocità stellare ($v_\star$): Da $0.5$a$1.5$ volte il valore fiduciale.
2. Densità superficiale del disco ($\Sigma_d$): Da $1$a$100$ volte il valore fiduciale, spaziando dai normali dischi sottili a densi dischi post-evento di distruzione tidale (TDE).
3. Raggio stellare ($R_\star$): Da $0.5$a$4$ volte il valore fiduciale, variando efficacemente il rapporto tra l'altezza di scala del disco e il raggio stellare ($H/R_\star$).
4. Profilo di densità verticale del disco: Modellato come una distribuzione gaussiana $\rho(z) \propto e^{-z^2/2\sigma^2}$, variando la concentrazione di massa verso il piano medio ($\sigma/R_\odot$).
5. Angolo di collisione locale ($i$): L'angolo tra il vettore di velocità relativa e il piano medio del disco, da $90^\circ$ (perpendicolare) a $30^\circ$ (obliquo).

Il trasporto della radiazione è stato trattato utilizzando l'approssimazione della diffusione a limite di flusso con opacità di scattering elettronico. La luminosità bolometrica emergente è stata calcolata integrando i flussi di diffusione e avvezione attraverso una fotosfera determinata dinamicamente.

Contributi Chiave e Risultati
Lo studio fornisce scalature empiriche per la luminosità di picco ($L_{\text{peak}}$) e la durata del flare ($\Delta t$) in funzione dei parametri del sistema. Le principali scoperte sono:

• Dinamica della collisione: La stella guida uno shock bow (a prua), creando due outflow (avanti e dietro) su lati opposti del disco. L'outflow anteriore è generalmente più massiccio ed energetico. La morfologia qualitativa dello shock bow e degli outflow rimane ampiamente invariata attraverso le variazioni di $v_\star$ e $\Sigma_d$, sebbene la normalizzazione termodinamica cambi.
• Velocità stellare ($v_\star$): Stelle più veloci producono flare più luminosi ($L_{\text{peak}} \propto v_\star^2$) a causa di tassi di riscaldamento per shock più elevati. Producono anche flare più lunghi ($\Delta t \propto v_\star^{1.4-2.2}$) perché viene iniettata più energia, sebbene la luminosità decada più rapidamente a causa di perdite adiabatiche più forti.
• Densità superficiale del disco ($\Sigma_d$): Dischi più densi aumentano la massa e la quantità di moto degli outflow ma hanno un effetto debole sulla luminosità di picco ($L_{\text{peak}} \propto \Sigma_d^{-0.1}$ per l'anteriore, $\propto \Sigma_d^{0.2}$ per il posteriore). Ciò è dovuto alla competizione tra l'aumento del riscaldamento da shock e le più forti perdite adiabatiche in ejecta otticamente più spessi. Tuttavia, i dischi più densi prolungano significativamente la durata del flare ($\Delta t \propto \Sigma_d^{0.5-0.6}$) a causa dei tempi di diffusione dei fotoni più lunghi.
• Raggio stellare ($R_\star$): Stelle più grandi interagiscono con più gas, producendo flare più luminosi ($L_{\text{peak}} \propto R_\star^{1.5}$ anteriore, $\propto R_\star^{0.7}$ posteriore) e più lunghi ($\Delta t \propto R_\star^{1.5}$). L'aumento di $R_\star$ accentua l'asimmetria di quantità di moto tra gli outflow anteriore e posteriore.
• Profilo di densità verticale ($\sigma$): Dischi più concentrati centralmente ( $\sigma$ minore) producono flare più luminosi e brevi. Gli strati centrali densi assorbono più energia da shock, mentre gli strati esterni a bassa densità permettono alla radiazione di uscire più efficientemente con ridotte perdite adiabatiche.
• Angolo di collisione ($i$): Attraversamenti più obliqui ( $i$ minore) riducono l'asimmetria di quantità di moto e luminosità tra i due outflow. L'outflow posteriore diventa significativamente più luminoso ($L_{\text{peak}} \propto (\sin i)^{-0.8}$) e i flare diventano più lunghi ($\Delta t \propto (\sin i)^{-0.9 \text{ a } -1}$) a causa del tempo di interazione più lungo e di una redistribuzione più simmetrica del gas riscaldato.

Applicazione a GSN 069
Gli autori hanno applicato le loro scalature empiriche alla sorgente ben caratterizzata GSN 069. Campionando lo spazio dei parametri per far corrispondere le luminosità di picco osservate ($\sim 3 \times 10^{42}$ erg s$^{-1}$), le durate ($\sim 1$ h) e il rapporto di ampiezza forte-debole, hanno identificato due regioni percorribili:

1. Soluzione Preferita: Un'orbita retrograda che collide con un disco post-TDE denso e verticalmente concentrato ($\Sigma_d \sim 45 \times$ il fiduciale, $\sigma/R_\odot \sim 0.3$). Questo modello riproduce bene le proprietà del flare osservate.
2. Soluzione Alternativa: Un'orbita retrograda con una densità di disco sottile standard ma una geometria di collisione estremamente radente ($i \sim 6^\circ$). Questa soluzione è meno accurata e richiede parametri estremi.

Gli autori notano che la preferenza per un disco denso e verticalmente concentrato è coerente con l'evidenza osservativa di attività simile a un TDE in GSN 069.

Significato e Limitazioni
L'articolo sostiene di fornire il primo studio sistematico di idrodinamica della radiazione in 3D che colleghi i parametri locali di collisione alle curve di luce delle QPE. Le scalature empiriche derivate servono come ponte tra le simulazioni locali e i modelli globali di QPE, permettendo di predire le proprietà dei flare basandosi sui parametri del sistema. Gli autori sottolineano che i loro risultati suggeriscono che combinazioni specifiche di parametri possono riprodurre le caratteristiche ampiezze, durate e duty cycle delle QPE.

Lo studio riconosce modestamente alcune limitazioni: i modelli assumono scalature a un parametro separabili, trascurano la rotazione e il taglio (shear) del disco, trattano la stella come una sfera rigida e non tengono conto degli effetti dell'angolo di visuale o dell'evoluzione spettrale (es. Comptonizzazione). Di conseguenza, i risultati sono presentati come regioni fisicamente motivate dello spazio dei parametri piuttosto che come modelli unici e definitivi. Gli autori suggeriscono che i lavori futuri dovrebbero incorporare strutture stellari realistiche, dinamiche del disco e luminosità dipendenti dall'angolo di visuale per quantificare pienamente la capacità del modello di riprodurre la diversità delle sorgenti QPE osservate.