Simulation based parameter space for shock in transonic, sub-Keplerian accretion flow onto non-rotating black holes

Sintesi Tecnica: Spazio dei Parametri Basato su Simulazioni per gli Shock in Flussi di Accrezione Transonici, Sub-Kepleriani

Definizione del Problema
I flussi di accrezione transonici, sub-kepleriani e non dissipativi che precipitano su buchi neri sono governati da due parametri conservati: l'energia specifica ($\epsilon$) e il momento angolare specifico ($l$). Mentre i modelli analitici basati sull'analisi di Rankine-Hugoniot (RH) e sull'equilibrio verticale hanno identificato uno spazio di parametri limitato in cui possono formarsi shock stazionari, precedenti simulazioni numeriche multidimensionali hanno suggerito che gli shock possano verificarsi anche al di fuori di questi limiti teorici. Tali simulazioni hanno inoltre indicato che gli shock possono essere dinamici (oscillatori) piuttosto che stazionari. Tuttavia, mancava una mappatura sistematica dello spazio dei parametri degli shock basata su simulazioni. Questo lavoro affronta la discrepanza tra le previsioni analitiche e la realtà numerica, scansionando sistematicamente lo spazio dei parametri $(l, \epsilon)$ per identificare dove avviene la formazione di shock, caratterizzare la natura di tali shock (stazionari vs oscillatori) e analizzare la dinamica del flusso risultante.

Metodologia
Gli autori hanno eseguito oltre 2 ever 280 simulazioni idrodinamiche ideali bidimensionali in un dominio cilindrico $R-Z$, assumendo l'assialsimmetria. Le simulazioni hanno utilizzato un potenziale pseudo-newtoniano ($\Phi(r) = -1/2(r-1)$) per modellare un buco nero non rotante.

• Schema Numerico: Il codice ha impiegato un metodo a volumi finiti con accuratezza del secondo ordine, utilizzando un limitatore di van Leer per la ricostruzione spaziale, un solutore di Riemann HLL per i flussi interfacciali e uno schema Runge-Kutta a due stadi con stabilità preservata (SSP) per l'integrazione temporale.
• Condizioni al Contorno: L'afflusso è stato imposto al limite radiale esterno ($R_{out}$) con variabili primitive costanti derivate da soluzioni analitiche per coppie specifiche di $(l, \epsilon)$, garantendo un'iniezione supersonica per prevenire il feedback acustico. Sono state applicate condizioni di deflusso al limite verticale superiore ($Z_{out}$). Sono state utilizzate condizioni di riflessione sull'asse di rotazione ($R=0$) e sul piano equatoriale ($Z=0$). La materia all'interno di un raggio $r_{BH} = 1.8$ è stata assorbita.
• Scansione dei Parametri: Lo studio ha scansionato un dominio rettangolare di $l \in [1.4, 1.9]$ e $\epsilon \in [0.0001, 0.0061]$ utilizzando 286 punti. I domini di simulazione sono stati adattati per ogni punto in base alla posizione dell'estremità più esterna del punto sonico di tipo X. Le simulazioni sono state eseguite fino a $t=50000$ per garantire il comportamento a regime stazionario o periodico.

Contributi Principali e Risultati
Il contributo primario è l'identificazione di uno spazio di parametri per la formazione di shock significativamente più ampio della regione prevista analiticamente. Gli autori hanno categorizzato lo spazio dei parametri in cinque regioni distinte in base alle caratteristiche delle soluzioni:

1. Regione A (Shock Stazionari): Questa regione contiene parametri in cui si formano shock stazionari e stabili. Lo shock crea uno strato limite supportato dalla barriera centrifuga (CENBOL). Gli autori notano che il confine identificato numericamente per questa regione si estende a valori di $l$ inferiori rispetto alle previsioni analitiche. Questa estensione è attribuita ai moti verticali nella simulazione che accumulano materia vicino all'equatore, aumentando la pressione e facilitando la formazione dello shock anche quando le assunzioni di equilibrio verticale nei modelli analitici lo proibirebbero.
2. Regione E (Nessuno Shock): Situata a valori di $l$ più bassi, queste soluzioni presentano una barriera centrifuga che causa un accumulo di densità e un salto supersonico-supersonico, ma non si forma un vero shock (transizione da supersonico a subsonico). Il flusso rimane supersonico durante tutto il percorso e accresce sul buco nero.
3. Regioni B e C (Shock Oscillatori): Queste regioni contengono parametri in cui si formano shock che però oscillano.
   • Regione B: Caratterizzata da oscillazioni ad alta frequenza e piccola ampiezza, confinate principalmente alla parte inferiore della superficie dello shock (vicino all'equatore). La parte superiore del CENBOL rimane statica.
   • Regione C: Caratterizzata da oscillazioni a grande ampiezza e bassa frequenza che coinvolgono l'intera superficie dello shock. Queste soluzioni presentano vortici turbolenti significativi nella regione post-shock.
4. Regione D (Dominata dall'Outflow): Ad alti valori di $l$, la barriera centrifuga è così forte da impedire l'accrezione. Il flusso sviluppa uno shock che si propaga verso l'esterno raggiungendo il bordo del dominio, lasciando dietro di sé un flusso subsonico. Più del 86% della materia iniettata viene espulsa come outflow.

Dipendenza dai Parametri

• Posizione dello Shock: Per un $\epsilon$ fissato, la posizione dello shock generalmente si sposta più lontano dal buco nero all'aumentare di $l$. La dipendenza da $\epsilon$ per un $l$ fissato è non monotona; in alcuni casi, la posizione dello shock aumenta con $\epsilon$, mentre in altri (particolarmente a $l$ più bassi), rimane costante o diminuisce leggermente.
• Outflow: L'outflow è generato in quasi tutte le simulazioni a causa della combinazione tra la barriera centrifuga e i gradienti di pressione termica. Il rapporto tra outflow e inflow aumenta significamente al crescere di $l$, passando da un regime dominato dall'assorbimento (Regioni A ed E) a un regime dominato dall'outflow (Regione D).

Significatività e Rivendicazioni
L'articolo sostiene che le simulazioni numeriche multidimensionali rivelano uno spazio di parametri per la formazione di shock più ampio di quanto precedentemente compreso attraverso i metodi analitici. Nello specifico, l'inclusione della dinamica verticale permette la formazione di shock a momenti angolari inferiori rispetto a quelli previsti dai modelli di equilibrio verticale.

Gli autori evidenziano la rilevanza astrofisica di queste scoperte:

• Formazione della Corona: Il flusso post-shock compresso e geometricamente spesso (CENBOL) agisce come una regione calda e otticamente sottile capace di compiere l'effetto Compton inverso sui fotoni a bassa energia provenienti da un disco Kepleriano (o auto-Comptonizzazione se il disco è assente). Ciò fornisce un meccanismo naturale per la formazione della corona osservata nei binari X e nei nuclei galattici attivi negli stati spettrali "hard".
• Oscillazioni Quasi-Periodiche (QPO): L'identificazione delle regioni di shock oscillatorio (B e C) supporta il legame teorico tra le oscillazioni dello shock e le QPO osservate nella radiazione ad alta energia.
• Outflow: Le simulazioni producono in modo autoconsistente degli outflow, che sono frequentemente osservati negli stati spettrali "hard", suggerendo che il meccanismo dello shock sia un driver unificato sia per la corona che per l'outflow.

Gli autori mantengono la modestia riguardo all'ambito del loro lavoro, notando che i risultati sono specifici per buchi neri non rotanti utilizzando un potenziale pseudo-newtoniano e l'idrodinamica ideale. Riconoscono che l'inclusione di viscosità, perdite radiative, campi magnetici o effetti relativistici generali per buchi neri rotanti potrebbe alterare i confini dello spazio dei parametri. Tuttavia, affermano che l'attuale lavoro stabilisce una base fondamentale, basata su simulazioni, per la mappatura della formazione di shock che si estende oltre gli attuali limiti analitici.