Accretion of Primordial Black Holes in Stellar Interiors

Sintesi Tecnica: Accrezione di Buchi Neri Primordiali nei Cuori Stellari

Problematica
I Buchi Neri Primordiali (PBH) rimangono un candidato valido per la materia oscura, particolarmente nell'intervallo di massa sub-lunare scarsamente vincolato ($10^{-16} M_\odot \lesssim M_{BH} \lesssim 10^{-10} M_\odot$). Le stelle forniscono un laboratorio unico per sondare i PBH in questo regime, poiché un PBH catturato e immerso nel nucleo stellare accreerebbe il gas circostante, alterando potenzialmente l'evoluzione della stella. Studi precedenti (ad esempio, Bellinger et al. 2023) hanno modellato questo scenario utilizzando codici di evoluzione stellare 1D, ma si sono basati su assunzioni derivate dalla teoria dei dischi di accrezione sottili. Nello specifico, assumevano un'efficienza radiativa fissa $\eta \approx 0,08$ e applicavano il limite di luminosità di Eddington a tutte le masse. Tuttavia, per i PBH in questo intervallo di massa, il raggio di Bondi ($r_B$) è di ordini di grandezza inferiore al raggio di circolarizzazione stellare, impedendo la formazione del disco. Invece, l'accrezione è sferica, e il gas viene compresso e riscaldato a temperature relativistiche ($T \sim 10^{10}–10^{12}$ K) vicino al raggio di Schwarzschild ($r_S$). In questo regime, la microfisica (bremsstrahlung relativistica, processi di coppia, scattering Klein–Nishina) e l'opacità del flusso differiscono fondamentalmente dal paradigma del disco sottile, rendendo necessaria una determinazione auto-coerente del tasso di accrezione ($\dot{M}$), dell'efficienza radiativa ($\eta$) e della struttura del flusso.

Metodologia
Gli autori risolvono le equazioni di Euler sferiche dipendenti dal tempo con un termine di sorgente di raffreddamento implicito per modellare il flusso di accrezione su un PBH immerso in un nucleo di tipo solare.

• Microfisica: Il termine di raffreddamento $\epsilon(\rho, T)$ include la bremsstrahlung relativistica (elettrone-ione ed elettrone-elettrone) e l'annichilazione elettrone-positrone. L'opacità tiene conto dello scattering Thomson, dell'assorbimento free-free e della soppressione Klein–Nishina alle alte energie fotoniche.
• Approccio Numerico: Gli autori utilizzano radbondi, un pacchetto Python open-source, per evolvere il flusso da un profilo di Bondi adiabatico iniziale verso uno stato stazionario auto-coerente. Questo approccio dipendente dal tempo è necessario perché, nel regime dominato dal raffreddamento, il punto sonico diventa un focus (autovalori complessi) piuttosto che un punto di sella, rendendo invalidi i classici metodi di shooting per ODE stazionarie.
• Gestione dei Regimi:
  • Otticamente Sottile ($M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): Il solutore di Euler calcola direttamente l'efficienza radiativa fisica e il tasso di accrezione.
  • Otticamente Spesso ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): La sfera di Bondi diventa otticamente spessa, entrando in un regime di intrappolamento dei fotoni. Gli autori adottano la prescrizione analitica di Begelman (1979), in cui i fotoni prodotti localmente all'interno del raggio di intrappolamento ($r_{trap}$) vengono trascinati verso il buco nero, mentre la luminosità in uscita è autoregolata a $L_{esc} \lesssim 0,6 L_{Edd}$.
• Feedback: Lo studio modella in modo auto-coerente il feedback radiativo, dove la luminosità depositata al di fuori della sfera di Bondi riscalda il plasma ambiente. Questo riscaldamento è moderato dal trasporto convettivo (Teoria della Lunghezza di Misura), prevenendo gli estremi aumenti di temperatura predetti dalla pura diffusione radiativa.

Risultati Chiave
Lo studio identifica tre distinti regimi di accrezione sferica e deriva i seguenti risultati:

1. Tre Regimi di Accrezza:

   • Hot Bondi ($M_{BH} \lesssim 10^{-14} M_\odot$): Il raffreddamento è dinamicamente trascurabile. Il flusso è quasi adiabatico, raggiungendo temperature $T \sim 10^{11}$ K vicino a $r_S$. Il flusso interno diventa collisionless prima di raggiungere le temperature di produzione di coppie, sopprimendo l'equilibrio di coppie. La bremsstrahlung domina il raffreddamento.
   • Raffreddamento da Bremsstrahlung ($10^{-14} \lesssim M_{BH} \lesssim 5 \times 10^{-13} M_\odot$): Il raffreddamento diventa efficiente, drenando l'energia termica prima che il gas raggiunga le temperature di produzione di coppie. Il flusso si avvicina a uno stato isotermico ($T \approx T_\infty$) e il tasso di accrezione è potenziato di un fattore $\sim 2–7$ rispetto al tasso di Bondi adiabatico.
   • Intrappolamento dei Fotoni ($M_{BH} \gtrsim 5 \times 10^{-13} M_\odot$): La sfera di Bondi è otticamente spessa. Sebbene la luminosità emessa localmente aumenti rapidamente con la massa, la luminosità in uscita è limitata dal limite di Begelman ($L_{esc} \lesssim 0,6 L_{Edd}$). Fondamentalmente, il tasso di accrezione rimane vicino al valore di Bondi ($\dot{M} \propto M_{BH}^2$) perché non esiste un limite di Eddington sferico per l'afflusso; la pressione di radiazione modula solo la gravità effettiva.

2. Efficienza Radiativa: L'efficienza radiativa calcolata auto-coerentemente è $\eta \approx 10^{-2}$ nell'intervallo otticamente sottile, un ordine di grandezza inferiore al valore canonico del disco sottile ($\eta \approx 0,08$) assunto nella letteratura precedente. L'efficienza raggiunge il picco vicino a $10^{-15} M_\odot$ e poi diminuisce nei limiti di raffreddamento e isotermico.

3. Traiettorie di Crescita: Poiché il limite di Eddington sferico non si applica, la crescita dei PBH rimane super-esponenziale ($\dot{M} \propto M_{BH}^2$) durante tutto il regime di intrappolamento dei fotoni. Il collo di bottiglia della crescita è determinato interamente dalla lenta fase Hot Bondi alle masse più basse.

4. Massa Iniziale Critica: La massa iniziale critica richiesta affinché un PBH consumi una stella di massa solare entro un tempo di Hubble è $M_{0,crit} \sim 10^{-16} M_\odot$. Questo è cinque ordini di grandezza inferiore alle stime precedenti ($\sim 10^{-11} M_\odot$) che assumevano un'efficienza maggiore e un collo di bottiglia di Eddington.

5. Campi Magnetici: Gli autori riscontrano che l'accumulo di flusso magnetico è trascurabile. Il campo magnetico ambiente nel nucleo solare è troppo debole per arrestare dinamicamente il flusso, e la riconnessione magnetica limita l'accumulo di flusso a livelli molto al di sotto dell'importanza dinamica ($P_{mag}/P_{ram} \ll 1$).

Significato e Rivendicazioni
Il lavoro sostiene di fornire il primo calcolo auto-coerente dell'efficienza radiativa e dei tassi di accrezione per l'accrezione sferica di PBH attraverso l'intero intervallo di massa sub-lunare. Sostituendo le assunzioni del disco sottile con un trattamento microfisico del flusso sferico, gli autori dimostrano che:

• L'efficienza radiativa è significativamente inferiore a quanto assunto in precedenza.
• Non esiste una fase limitata dall'effetto Eddington per l'accrezione sferica; il flusso rimane super-Eddington in termini di flusso di massa senza essere bloccato dalla pressione di radiazione.
• Di conseguenza, i vincoli sui PBH come materia oscura dalle osservazioni stellari sono molto più stretti di quanto precedentemente pensato, poiché i PBH con masse iniziali anche di soli $10^{-16} M_\odot$ possono crescere fino a consumare le loro stelle ospiti entro un tempo di Hubble.

Gli autori notano che il loro lavoro si concentra sulla microfisica dell'accrezione e non include la reazione inversa sulla struttura stellare globale (che richiederebbe l'accoppiamento con un codice di evoluzione stellare come MESA), ma forniscono le relazioni $\eta(M_{BH})$ e $\dot{M}(M_{BH})$ necessarie per consentire tali studi futuri. Inoltre, rimandano la discussione sulla formazione del disco (che avviene a masse più elevate o in stelle con rotazione più rapida) a un articolo complementare.