Supermassive Black Hole Growth in Massive Galaxies at Cosmic Dawn

Sintesi Tecnica: Crescita di Buchi Neri Supermassicci in Galassie Massicce all'Alba Cosmica

Enunciato del Problema
Recenti osservazioni del Telescopio Spaziale James Webb (JWST) hanno rivelato una popolazione di galassie massicce e luminose all'Alba Cosmica ($z \gtrsim 9$) con alte masse stellari ($M_\star \gtrsim 10^8 M_\odot$) e tassi di formazione stellare (SFR). Diverse di queste galassie esibiscono caratteristiche spettrali indicative di nuclei galattici attivi (AGN), suggerendo la presenza di buchi neri supermassicci (SMBH) sovradimensionati rispetto alle relazioni di scala locali. I modelli standard di formazione galattica faticano a spiegare la rapida crescita di questi SMBH e gli elevati regimi di efficienza della formazione stellare osservati. Sebbene proposte teoriche invocano semi pesanti o accrezione super-Eddington, l'interazione tra la crescita degli SMBH, il feedback stellare e il feedback degli AGN nel mezzo interstellare (ISM) turbolento e multifase delle prime galassie massicce rimane scarsamente compresa.

Metodologia
Gli autori hanno condotto una serie di 11 simulazioni idrodinamiche cosmologiche "zoom-in" utilizzando il codice di raffinamento adattivo della griglia (AMR) ramses. Le simulazioni si concentrano su una specifica sovradensità che collassa in un alone di $\sim 10^{11} M_\odot$ entro $z \sim 9$, ospitando una galassia massiccia (Massive Dawn Galaxy o MDG) con $M_\star \gtrsim 10^8 M_\odot$.

Componenti metodologiche chiave:

• Modello di Formazione Stellare: A differenza delle simulazioni standard che utilizzano la legge di Kennicutt-Schmidt, questo lavoro impiega un modello subgrid multi-freefall basato sulla turbolenza e fisicamente motivato. Ciò consente efficienze di formazione stellare locale (SFE) variabili guidate dall'energia cinetica turbolenta dell'ISM, capaci di raggiungere $\gtrsim 10\%$.
• Implementazione degli SMBH: Gli SMBH sono modellati come particelle "sink" seminate a $z \sim 15$ quando le masse locali di gas e stelle superano $10^5 M_\odot$ all'interno di un alone di $10^{11} M_\odot$. Sono stati testati semi con masse di $10^4 M_\odot$ e $10^5 M_\odot$.
• Fisica dell'Accrezione: L'accrezione segue una formula di Bondi-Hoyle-Lyttleton modificata con un fattore di boost per tenere conto delle fluttuazioni di densità non risolte. Il tasso di accrezione è limitato da un limite di Eddington, che viene variato nelle simulazioni tramite un fattore $\lambda_{\text{edd}}$ (valori: 1, 3 e $\infty$ per nessun limite).
• Meccanismi di Feedback:
  • Feedback Stellare: Include fotoionizzazione termica (tipo precoce) e supernovae di Tipo II (tipo tardivo). Una variante a "basso feedback stellare" è stata testata rimuovendo la fotoionizzazione per isolarne l'impatto sulla struttura dell'ISM.
  • Feedback degli AGN: Modellato come iniezione di energia termica isotropa nella sfera del sink con un'efficienza di accoppiamento $\epsilon_c$ (testata al 0% e al 15%).
• Suite di Simulazioni: Gli 11 run variano sistematicamente la forza del feedback degli AGN, la massa del seme, i limiti massimi di accrezione e l'efficienza del feedback stellare per isolare i loro effetti sulla co-evoluzione SMBH-galassia.

Risultati Chiave

1. Accrezione Stocastica Guidata dall'ISM (Festa e Carestia):
   Nell'ambiente MDG fiduciale, il feedback stellare genera un ISM turbolento e multifase composto da patch fredde-dense ($n_H \sim 10^3 \text{ cm}^{-3}$) e calde-diffuse ($n_H \sim 10 \text{ cm}^{-3}$). La sfera del sink dell'SMBH campiona stocasticamente queste fasi.

   • Modalità Festa (Feast Mode): Quando il sink entra in una patch fredda-densa, la velocità del suono scende, permettendo all'accrezione di raggiungere il limite (super-)Eddington.
   • Modalità Carestia (Starvation Mode): Quando il sink entra in una patch calda-diffusa, la velocità del suono aumenta e il tasso di accrezione di Bondi scende significativamente.
   • Nel modello fiduciale, l'SMBH trascorre circa il 50% del tempo in modalità carestia ($\alpha \approx 0.5$) dopo l'inizio della formazione stellare ($t \sim 300$ Myr), sopprimendo efficacemente il tasso di crescita medio della metà rispetto a un ISM uniforme.

2. Condizioni per l'Autoregolazione:
   Il feedback degli AGN può autoregolare la crescita dell'SMBH solo se il buco nero diventa abbastanza massiccio che il suo riscaldamento domini il raffreddamento all'interno della sfera del sink ($M_{\text{sink}} > M_{\text{sink,crit}}^{\text{cool}}$).

   • Questa massa critica dipende dalla densità del gas e dalla metallicità. Nell'ambiente turbolento dell'MDG, l'SMBH deve superare questa soglia per $\gtrsim 50\%$ del tempo trascorso nelle sacche calde-diffuse per ottenere un riscaldamento incontrollato e slegare il gas dall'alone.
   • Le simulazioni con accrezione limitata da Eddington ($\lambda_{\text{edd}}=1$) e masse di seme standard spesso non raggiungono l'autoregolazione perché la soglia della massa critica sale con la metallicità più velocemente di quanto l'SMBH possa crescere, o l'SMBH trascorre tempo insufficiente sopra tale soglia.
   • L'autoregolazione è raggiunta nelle simulazioni con masse di seme più elevate ($10^5 M_\odot$) o accrezione super-Eddington ($\lambda_{\text{edd}} \ge 3$), a condizione che l'efficienza di accoppiamento termico sia sufficientemente bassa ($\epsilon_c \lesssim 0.15\%$) per corrispondere alle masse osservate degli SMBH.

3. Mancanza di Spegnimento su Scala Galattica:
   Contrariamente alle aspettative dei modelli a basso redshift, il feedback degli AGN in queste simulazioni non spegne il tasso di formazione stellare globale (SFR) né arresta la crescita della massa stellare della galassia ospite.

   • L'alone accretesce gas freddo dai filamenti cosmici a un tasso ($\dot{M}_{\text{halo}}$) che supera significativamente il tasso massimo di crescita dell'SMBH ($\dot{M}_{\text{SMBH}}$).
   • Poiché la scala temporale di crescita dell'alone ($\tau_{\text{halo}}$) è comparabile o più breve della scala temporale di crescita dell'SMBH, il gas freddo fresco ripristina continuamente l'ISM, sovrastando qualsiasi riscaldamento dal feedback degli AGN. La galassia continua a crescere esponenzialmente indipendentemente dall'attività dell'SMBH.

4. Confronto con le Osservazioni:
   La suite di simulazioni che meglio corrisponde alla popolazione osservata di SMBH sovradimensionati all'Alba Cosmica (es. UHZ-1, GN-z11) è la run super_edd_no_AGN ($M_{\text{seed}} \sim 10^4 M_\odot$, $\lambda_{\text{edd}} = 3$, assenza di feedback degli AGN).

   • Per far corrispondere le masse osservate con SMBH autoregolati, l'efficienza di accoppiamento termico deve essere estremamente bassa ($\epsilon_c \lesssim 0.15\%$).
   • I risultati suggeriscono che l'accrezione super-Eddington sia probabilmente necessaria affinché gli SMBH crescano abbastanza da entrare nel regime autoregolato senza richiedere masse di seme così elevate da $10^6 M_\odot$.

Significato e Conclusioni
L'articolo sostiene che la crescita degli SMBH nelle galassie massicce all'Alba Cosmica è regolata principalmente dalla natura stocastica dell'ISM turbolento e multifase guidato dal feedback stellare, piuttosto che dal solo feedback degli AGN. Gli autori concludono che:

• Gli SMBH in questi ambienti sono frequentemente "affamati" dalle sacche di gas caldo-diffuso, limitando la loro crescita a meno che non possano sostenere un'accrezione super-Eddington.
• Il feedback degli AGN è inefficace nello spegnimento della formazione stellare su scale galattiche durante questa epoca perché l'rapida accrezione di gas nell'alone ripristina continuamente la scorta di combustibile.
• Raggiungere i massicci SMBH osservati con JWST richiede probabilmente una combinazione di accrezione super-Eddington e un accoppiamento inefficiente del feedback degli AGN, o semi molto massicci, per superare la crescente soglia della massa critica di raffreddamento prima che la galassia evolva in uno stato più stabile, simile a quello a basso redshift.

Gli autori notano che i loro risultati dipendono dalla risoluzione e che lavori futuri che incorporino la chimica termica fuori equilibrio, il trasferimento radiativo e una risoluzione spaziale più elevata ($\sim 1$ pc) sono necessari per validare pienamente queste scoperte e affinare le scale fisiche del deposito del feedback.