Different effects of the Lorentz and Gaussian bump functions on the formation of primordial black holes and secondary gravitational waves

Sintesi Tecnica: Diversi Effetti delle Funzioni di Bump di Lorentz e Gaussiana su Buchi Neri Primordiali e Onde Gravitazionali Secondarie

Definizione del Problema
I Buchi Neri Primordiali (PBH) sono candidati validi per la materia oscura fredda (DM) e potenziali sorgenti di eventi di onde gravitazionali (GW). La loro formazione è spesso attribuita all'amplificazione delle perturbazioni scalari durante l'inflazione, che può essere ottenuta introducendo un "bump" locale al potenziale di inflazione base, $V(\phi) = V_b(\phi)(1 + f(\phi))$. Sebbene varie funzioni di bump (Gaussiana, iperbolica, Lorentz) siano state utilizzate in letteratura, non esiste un principio fisico sottostante che detti la scelta specifica di $f(\phi)$. Di conseguenza, rimane poco chiaro come la specifica forma funzionale del bump influenzi lo spettro di potenza della curvatura risultante, l'abbondanza di PBH e le onde gravitazionali secondarie indotte (SIGW). Questo articolo affronta la mancanza di un'analisi comparativa tra le funzioni di bump ampiamente utilizzate, specificamente le forme Lorentz e Gaussiana, sotto identiche condizioni di base.

Metodologia
Gli autori utilizzano il modello di inflazione di Starobinsky come potenziale base $V_b(\phi)$, scelto per la sua compatibilità con le osservazioni del Fondo Cosmico a Microonde (CMB). Introducono correzioni locali di tipo bump $f(\phi)$ in due forme:

1. Bump di Lorentz: $f_L(\phi) = \frac{b}{1 + (\frac{\phi-\phi_0}{c})^2}$
2. Bump Gaussiano: $f_G(\phi) = b e^{-\frac{(\phi-\phi_0)^2}{2c^2}}$

Per isolare gli effetti della forma funzionale, gli autori mantengono identici il parametro di ampiezza ($b$), la posizione del picco ($\phi_0$) e il parametio di larghezza ($c$) per entrambi i casi. Risolvono numericamente l'equazione di Mukhanov-Sasaki per calcolare lo spettro di potenza scalare $\mathcal{P}_\zeta$. Utilizzando il formalismo di Press-Schechter, calcolano l'abbondanza di PBH ($Y_{PBH}$) e lo spettro di densità di energia delle onde gravitazionali indotte da scalari ($\Omega_{GW}$). Vengono testati tre distinti set di parametri (I, II, III) per verificare la generalità dei risultati.

Contributi Chiave e Risultati

• Potenziamento dello Spettro di Potenza: Lo studio rivela una differenza nuova e significativa nella capacità delle due funzioni di potenziare lo spettro di potenza. Il bump di Lorentz esibisce una "coda più grassa" rispetto al bump Gaussiano, che decade esponenzialmente. Questa differenza strutturale crea un plateau più ampio e piatto nel potenziale di inflazione, portando a un periodo più esteso di ultra-slow roll. Di conseguenza, il caso Lorentz produce uno spettro di potenza di picco ($\mathcal{P}_\zeta \sim \mathcal{O}(10^{-2})$) che è di ordini di grandezza superiore rispetto al caso Gaussiano ($\mathcal{P}_\zeta \sim \mathcal{O}(10^{-7})$) a parità di parametri.
• Abbondanza di Buchi Neri Primordiali: A causa dello spettro di potenza significativamente più elevato, il bump di Lorentz genera un'abbondanza sostanziale di PBH, mentre il bump Gaussiano produce un'abbondanza trascurabile (effettivamente zero negli scenari testati).
  • Set di Parametri I: Genera PBH di massa stellare ($\sim 37 M_\odot$), potenzialmente rilevanti per gli eventi di fusione di binari di buchi neri rilevati da LIGO.
  • Set di Parametri II: Produce PBH con masse $\sim 10^{-5} M_\odot$, che potrebbero spiegare gli eventi di microlensing ultrabrevi di tipo OGLE.
  • Set di Parametri III: Produce PBH con masse $\sim 10^{-12} M_\odot$ con un'abbondanza $Y_{PBH} \approx 1$, suggerendo che potrebbero costituire tutta la materia oscura.
• Onde Gravitazionali Indotte da Scalari (SIGWs): Le perturbazioni scalari potenziate nel caso Lorentz inducono una densità di energia delle SIGW significativamente più alta rispetto al caso Gaussiano. Le frequenze di picco di questi segnali rientrano negli intervalli di sensibilità di vari rilevatori attuali e futuri:
  • Set I: Rilevabile da array di timing di pulsar (EPTA, PPTA, SKA) e NANOGrav.
  • Set II: Rilevabile da SKA.
  • Set III: Rilevabile da interferometri spaziali (LISA, Taiji, TianQin, DECIGO).

Significatività e Rivendicazioni
L'articolo sostiene che la scelta della funzione di bump non è una mera comodità matematica ma ha profonde conseguenze fisiche. Il decadimento di legge di potenza della funzione di Lorentz permette un potenziamento delle perturbazioni più efficace rispetto al decadimento esponenziale della funzione Gaussiana. Questa distinzione offre un potenziale elemento osservativo: futuri esperimenti che rilevano specifici intervalli di massa dei PBH o segnali SIGW potrebbero distinguere tra questi due scenari teorici.

Gli autori notano con modestia che l'uso del formalismo di Press-Schechter è una semplificazione e che la loro analisi assume statistiche gaussiane, trascurando potenziali effetti di non-gaussianità nelle code della funzione di densità di probabilità. Suggeriscono che lavori futuri che incorporino la teoria dei picchi e la non-gaussianità, nonché l'investigazione delle correzioni a un loop sulle scale CMB, fornirebbero una comprensione più accurata e completa dei meccanismi di formazione dei PBH.