Radiation-hydrodynamics of star-disc collisions for quasi-periodic eruptions

Sintesi Tecnica: Idrodinamica della Radiazione delle Collisioni Stella-Disco per le Eruzioni Quasi-Periodiche

Definizione del Problema
Le eruzioni quasi-periodiche (QPE) sono transienti di raggi X luminosi che si verificano in prossimità di buchi neri supermassicci (SMBH), caratterizzati da brillanti lampi separati da ore o giorni. Un modello teorico principale ipotizza che queste eruzioni siano il risultato di una stella su un'orbita inclinata che collide ripetutamente con un disco di accrezione. Sebbene studi precedenti abbiano esplorato le interazioni stella-disco, molti mancavano di un trasporto di radiazione dettagliato o si concentravano su regimi in cui la stella è significativamente perturbata o il tempo di attraversamento è troppo breve affinché il gas possa fluire attorno alla stella. Inoltre, l'origine fisica dei pattern osservati di lampi alternati "forte-debole" in alcune sorgenti QPE rimane non ancora pienamente spiegata dalla dinamica della collisione. Questo studio affronta la necessità di un trattamento idrodinamico della radiazione tridimensionale (3D) della collisione nel regime in cui la stella rimane ampiamente non perturbata e il gas scosso ha tempo sufficiente per fluire attorno ad essa.

Metodologia
Gli autori hanno eseguito una simulazione idrodinamica della radiazione 3D utilizzando il codice di smoothed particle hydrodynamics (SPH) Phantom. La simulazione è stata condotta in un sistema di riferimento che si muove solidalmente con il disco di accrezione per isolare la dinamica locale della collisione.

• Configurazione: Una stella sferica rigida ($R_\star = R_\odot$) si muove perpendicolarmente attraverso un disco di accrezione uniforme, otticamente spesso e dominato dalla pressione di radiazione ($H = 3R_\star$, $\rho_d = 3.7 \times 10^{-8}$ g cm$^{-3}$) a una velocità di $v_\star = 0.1c$. I parametri del sistema corrispondono a un SMBH di massa $10^6 M_\odot$ e a un periodo QPE di 4 ore.
• Fisica: Il gas segue un'equazione di stato adiabatica ($\Gamma = 5/3$). Il trasporto della radiazione è trattato tramite l'approssimazione della diffusione limitata dal flusso, assumendo l'equilibrio termodinamico locale (LTE). La stella è modellata come un corpo rigido che interagisce elasticamente con le particelle di gas; le forze gravitazionali dell'SMBH e l'autogravità sono trascurate poiché il tempo scala della collisione è molto più breve dei tempi scala orbitali.
• Ambito: Lo studio si concentra su una sezione localizzata del disco ($12R_\star \times 2H$) per risolvere la struttura dello shock e la formazione dell'outflow senza il costo computazionale di una simulazione globale.

Contributi Chiave
Questo lavoro presenta la prima simulazione idrodinamica della radiazione 3D di una collisione stella-disco mirata specificamente allo spazio dei parametri rilevante per la maggior parte delle QPE (dove il disco è abbastanza spesso da permettere al gas di fluire attorno alla stella). Lo studio fornisce un'analisi dettagliata di:

1. La formazione e l'evoluzione dell'onda d'urto (bow shock) e degli shock secondari.
2. La ridistribuzione asimmetrica del momento e dell'energia iniettati dalla stella.
3. Le proprietà distinte degli outflow verso l'avanti (forward) e verso il retro (backward).
4. Le curve di luce bolometriche e le temperature efficaci risultanti.

Risultati

• Dinamica dello Shock: La collisione genera un bow shock quasi paraboloidale che viaggia con la stella ed espande lateralmente. Uno shock secondario si forma nella scia della stella mentre i flussi di gas convergono. L'espansione laterale del fronte d'urto segue una scalatura a legge di potenza ($r_\perp \propto t^{0.4}$), leggermente più debole della previsione di stato stazionario $t^{0.5}$, indicando che il sistema è ancora in fase di evoluzione verso uno stato stazionario.
• Outflow Asimmetrici: La collisione produce due outflow distinti su lati opposti del disco:
  • Outflow Forward (direzione del moto): Questo flusso è guidato da un breakout dello shock quasi sferico mentre la stella esce dal disco. Si espande più velocemente, contiene molta più massa (circa 10 volte l'outflow backward), trasporta più energia e momento, ed è approssimativamente due volte più luminoso dell'outflow backward.
  • Outflow Backward: Questo flusso manca di una componente di breakout quasi sferico. È formato da gas che fluisce attorno alla stella ed è reindirizzato dallo shock secondario o scappa dal bordo superiore del disco. È meno massiccio e meno luminoso.
• Luminosità e Curve di Luce: L'outflow forward raggiunge una luminosità bolometrica di picco di $\sim 10^{41}$ erg s$^{-1}$, circa il doppio dell'outflow backward. La curva di luce forward mostra una salita rapida seguita da una rottura verso un declino più lento, interpretata come un passaggio dalla radiazione che emerge dall'ejecta iniziale del breakout dello shock alla radiazione dal gas che scappa dalla cavità dopo un ritardo maggiore. Le temperature efficaci raggiungono il picco a $\sim 29$ eV per il flusso forward e valori leggermente inferiori per il flusso backward.
• Iniezione di Momento: La stella inietta momento preferenzialmente lungo la sua direzione di moto. Il momento verticale totale del gas scosso asintota a $P_z \approx -0.45 M_d v_\star$, riflettendo la natura anisotropa dell'iniezione di energia e gli angoli obliqui di riflessione del gas.

Significato e Rivendicazioni
Gli autori sostengono che l'asimmetria intrinseca nella dinamica della collisione spiega naturalmente i pattern di lampi "forte-debole" osservati in diverse sorgenti QPE (es. GSN 069). Poiché l'outflow forward è significativamente più luminoso dell'outflow backward, un osservatore che guarda il sistema da una specifica angolazione vedrebbe lampi alternati brillanti e deboli corrispondenti ai successivi passaggi della stella attraverso il disco.

Lo studio conclude che il modello di collisione stella-disco, nel regime in cui la stella rimane non perturbata, può riprodurre le luminosità di picco osservate ($10^{41}-10^{43}$ erg s$^{-1}$), le scale temporali (ore) e i pattern di luminosità alternati delle QPE. Gli autori notano limitazioni, tra cui l'uso di una stella rigida (trascurando lo stripping di massa e la deformazione), l'assunzione di LTE (che può venire meno in regimi poveri di fotoni) e l'omissione della rotazione globale del disco e della gravità dell'SMBH. Suggeriscono che lavori futuri che incorporino la fisica non-LTE e strutture stellari realistiche siano necessari per vincolare completamente il modello e affrontare le discrepanze nella durezza spettrale (simulata a $\sim 30$ eV rispetto agli osservati $\sim 100-200$ eV).