Evaporation of Primordial Black Holes with Multimodal Mass and Extended Spin Distributions: Cosmological Imprints on the Effective Number of Relativistic Species

Sintesi Tecnica: Evaporazione di Buchi Neri Primordiali con Distribuzioni di Massa Multimodali ed Estese

Enunciato del Problema
Si ipotizza che i Buchi Neri Primordiali (PBH) si formino attraverso vari meccanismi nell'Universo primordiale, ognuno dei quali imprime una funzione di massa distinta. Gli studi precedenti sull'evaporazione dei PBH e sulle loro conseguenze cosmologiche si sono ampiamente basati sull'approssimazione "monocromatica" (assumendo una singola massa caratteristica) o, al massimo, su singole distribuzioni estese. Questo approccio rischia di oscurare le firme caratteristiche di specifici meccanismi di formazione. Inoltre, sebbene l'impatto delle distribuzioni di massa estese sul numero effettivo di specie relativistiche ($N_{\text{eff}}$) sia stato esplorato, il caso di una popolazione genuinamente multimodale — dove operano simultaneamente molteplici canali di formazione con pesi derivati dalle probabilità di collasso fisico piuttosto che da parametri di adattamento (fit) — rimane non affrontato. Inoltre, l'interazione tra distribuzioni di massa estese, lo spin iniziale dei PBH e lo spin della radiazione oscura emessa ($s_{\text{DR}}$) richiede un trattamento sistematico per determinare la robustezza delle previsioni osservative.

Metodologia
Gli autori introducono FRISHBEE, un'estensione del codice pubblico FRISBHEE, per modellare l'evaporazione Hawking di popolazioni di PBH con distribuzioni di massa e spin estese e multimodali. La metodologia prevede:

1. Distribuzioni di Massa: Il framework implementa quattro distinte funzioni di massa corrispondenti a specifici canali di formazione fisica:

   • Log-Normale (LN): Derivante da esaltazioni localizzate nello spettro di potenza primordiale (es. ultra-slow-roll inflation).
   • Legge di Potenza (PL): Risultante da spettri invarianti di scala che collassano in un'era a equazione di stato costante.
   • Collasso Critico (CC): Caratterizzato da scaling critico vicino alla soglia di collasso con un taglio esponenziale.
   • Preheating Metrico (MP): Generato dall'amplificazione risonante delle perturbazioni durante le oscillazioni dell'inflatone.
   • Multimodale (MM): Una sovrapposizione lineare dei quattro canali sopra citati.

2. Pesatura Fisica: A differenza dei precedenti approcci fenomenologici, i contributi frazionari di densità energetica ($\alpha_i$) di ciascun canale nella miscela multimodale non sono parametri liberi. Essi sono derivati auto-coerentemente dalla probabilità di collasso primordiale, $\beta'_i$, che dipende dall'ampiezza dello spettro di potenza della curvatura primordiale alla scala caratteristica di ciascun canale. Tre scenari di pesatura benchmark vengono testati: Abbondanza Uguale (EA), Pesatura della Massa Media (MW) e Pesatura della Densità Numerica (NW).

3. Distribuzioni di Spin: La distribuzione congiunta massa-spin è fattorizzata. Vengono considerati due profili di spin: un profilo Gaussiano (centrato su $a_\star = 0.7$ o vicino all'estremo $a_\star = 0.99$) e la distribuzione universale guidata dalla fusione di Fishbach et al. (con picco vicino a $a_\star \approx 0.7$).

4. Equazioni di Evoluzione: Il codice risolve le equazioni accoppiate di Friedmann-Boltzmann in una formulazione di spazio delle fasi comobile. Traccia l'evoluzione della massa ($M$) e dello spin ($a_\star$) del PBH utilizzando i tassi di emissione Hawking di Kerr (inclusa la superradianza per $s_{\text{DR}} \geq 1$) e il contenuto di particelle del Modello Standard. Il trattamento tiene esplicitamente conto del disaccoppiamento dei neutrini, concentrandosi sulla finestra pre-BBN ($M_{\text{PBH}}^{\text{in}} \in \{10^5, 10^7, 10^8\}$ g) dove l'evaporazione avviene molto prima del disaccoppiamento dei neutrini ($T_{\text{evap}} \gg T_{\text{dec}}$), garantendo $\Delta N_{\text{eff}} > 0$.

Contributi Chiave

• Framework Multimodale: La prima implementazione sistematica di una popolazione di PBH multimodale in cui le frazioni dei sottopopolazioni sono fisicamente derivate dai parametri del modello inflazionario (ampiezze dello spettro di potenza) anziché adattate tramite fit.
• Sviluppo del Codice: Il rilascio di FRISHBEE, che incorpora le quattro funzioni di massa e il modulo di miscelazione multimodale, validato rispetto ai risultati monocromatici di FRISBHEE.
• Interazione Spin-Massa: Una quantificazione dettagliata di come la larghezza della distribuzione di massa moduli l'incremento indotto dallo spin di $\Delta N_{\text{eff}}$, in particolare in presenza di superradianza.

Risultati
Lo studio calcola $\Delta N_{\text{eff}}$ per varie distribuzioni, scale di massa e configurazioni di spin:

• Incremento della Distribuzione di Massa: Le funzioni di massa estese potenziano sistematicamente $\Delta N_{\text{eff}}$ rispetto al benchmark monocromatico nella finestra pre-BBN. La gerarchia di incremento è preservata attraverso le scale di massa: Log-Normale ($\sim 1.84\times$) > Multimodale ($\sim 1.77\times$) > Legge di Potenza ($\sim 1.69\times$) > Preheating Metrico ($\sim 1.15\times$) > Collasso Critico ($\sim 1.03\times$). Questa gerarchia si mantiene per tutti e tre gli scenari di pesatura (EA, MW, NW), confermando la robustezza del segnale multimodale.
• Effetti di Spin:
  • Per la radiazione oscura scalare ($s_{\text{DR}}=0$) o particelle a basso spin, l'aumento dello spin iniziale ($a_\star$) produce incrementi modesti ($R_i \lesssim 1.7$).
  • Per la radiazione oscura di spin-2 ($s_{\text{DR}}=2$) e spin quasi-estremi ($a_\star = 0.99$), la superradianza aumenta significativamente l'emissione. Tuttavia, l'incremento è diluito nelle distribuzioni di massa ampie perché i PBH più leggeri nella coda perdono rapidamente lo spin prima di evaporare.
  • Inversione: Per $s_{\text{DR}}=2$ e $a_\star=0.99$, il rapporto di incremento $\Delta N_{\text{eff}} / \Delta N_{\text{eff}}^{\text{mono}}$ scende sotto l'unità per le distribuzioni ampie (Log-Normale, Legge di Potenza). In questo specifico regime, l'approssimazione monocromatica sovrastima $\Delta N_{\text{eff}}$ perché assume che tutti i PBH mantengano uno spin quasi-estremo, mentre le popolazioni estese perdono lo spin più velocemente in media.
• Significatività Osservativa:
  • Per la radiazione oscura scalare ($s_{\text{DR}}=0$) con $M_{\text{PBH}}^{\text{in}} = 10^7$ g, la previsione monocromatica ($\Delta N_{\text{eff}} \approx 0.057$) è al di sotto della sensibilità progettata per CMB-S4 ($\sigma \approx 0.06$).
  • Al contrario, le distribuzioni estese Log-Normale, Legge di Potenza e Multimodale yield $\Delta N_{\text{eff}} \approx 0.096\text{--}0.105$, raggiungendo il livello di rilevamento di $\sim 1.6\text{--}1.7\sigma$ per CMB-S4.
  • Per la radiazione oscura di spin-2 ($s_{\text{DR}}=2$), il segnale assoluto è soppresso di un ordine di grandezza, rendendo tutte le distribuzioni non rilevabili alle attuali scale di massa.

Significatività
L'articolo stabilisce che la forma della funzione di massa dei PBH è un driver primario dell'impronta cosmologica dell'evaporazione di Hawking, non un semplice dettaglio secondario. Collegando le frazioni della popolazione multimodale direttamente all'ampiezza dello spettro di potenza primordiale, il framework trasforma $\Delta N_{\text{eff}}$ da un semplice limite di esistenza a un potenziale discriminatore tra scenari di formazione dei PBH a canale singolo e a canali multipli. I risultati indicano che i prossimi esperimenti CMB (CMB-S4, Simons Observatory) potrebbero potenzialmente distinguere tra popolazioni di PBH monocromatiche ed estese/multimodali, a condizione che la radiazione oscura sia scalare. Il lavoro evidenzia che l'osservabilità del segnale di Hawking è contingente sia dal canale di formazione che dallo spin della radiazione oscura emessa, rendendo necessaria una transizione oltre l'approssimazione monocromatica per ottenere previsioni fenomenologiche accurate.