Simulation of PBH formation in a matter-dominated universe

Sintesi Tecnica: Simulazione della Formazione di PBH in un Universo Dominato dalla Materia

Definizione del Problema
La formazione di Buchi Neri Primordiali (PBH) è altamente sensibile all'equazione di stato della componente cosmica dominante. Mentre la formazione di PBH durante l'era dominata dalla radiazione è relativamente ben compresa, il processo durante un primo' epoca dominata dalla materia rimane meno certo. In un'era dominata dalla materia, i gradienti di pressione sono trascurabili, permettendo alle regioni sovradense di collassare più facilmente. Tuttavia, il criterio di formazione è significativamente influenzato dalla dinamica non sferica, dall'inomogeneità e dall'ellitticità. Le precedenti stime analitiche si basano spesso su approssimazioni perturbative o sull'approssimazione di Zel'dovich, che potrebbero non essere valide nel regime altamente non lineare del collasso gravitazionale. Una sfida principale nell'investigare questo regime è l'implementazione della simmetria non sferica nelle simulazioni numeriche, particolarmente quando si trattano fluidi di polvere (dust) in cui le singolarità di incrocio dei gusci (shell-crossing) causano crash numerici prima che gli orizzonti dei buchi neri possano essere pienamente risolti.

Metodologia
Gli autori utilizzano la relatività numerica non lineare completa tramite il formalismo BSSN (Baumgarte-Shapiro-Shibata-Nakamura) per simulare la formazione di PBH in un background in espansione dominato dalla materia.

• Condizioni Iniziali: Le simulazioni utilizzano una forma funzionale della perturbazione di curvatura $\zeta$ che include l'ellitticità, creando configurazioni triassiche. I dati iniziali sono costruiti utilizzando soluzioni di espansione del gradiente a lunghezze d'onda elevate per un fluido di polvere, caratterizzate da un parametro di ampiezza $\mu$ e parametri di ellitticità ($e, p$).
• Gauge e Coordinate: Viene adottata una nuova condizione di slicing temporale, che transita da un gauge basato sul numero di e-folding a uno slicing di tipo "1+log" per evitare le singolarità man mano che l'ingresso dell'orizzonte si avvicina. Viene utilizzato un sistema di coordinate non cartesiano per ingrandire la regione centrale della scatola di simulazione, migliorando la risoluzione dove avviene il collasso.
• Descrizioni della Materia: Vengono testati due distinti modelli di materia:
  1. Fluido di Polvere (Dust Fluid): Una descrizione di un fluido privo di pressione. Gli autori notano che questa descrizione conduce inevitabilmente a singolarità di incrocio dei gusci dove la densità del fluido diverge, causando il crash del codice numerico.
  2. Particelle Collisionless: Per aggirare il fallimento dell'approssimazione del fluido, gli autori implementano un approccio Particle-in-Cell (PIC) all'interno della relatività numerica completa. Alle particelle vengono assegnate masse per riprodurre la desiderata inomogeneità di densità, e il loro moto è governato dalle equazioni geodetiche. Ciò permette alla simulazione di proseguire oltre il punto in cui gli elementi del fluido incrocerebbero le proprie traiettorie.

Risultati Chiave

• Validazione della Simmetria Sferica: Il codice è stato validato rispetto alla soluzione Lemaître-Tolman-Bondi (LTB) per il collasso della polvere con simmetria sferica. I risultati numerici hanno riprodotto con successo la soluzione analitica fino al punto di formazione dell'orizzonte e hanno identificato correttamente la soglia tra singolarità nude locali (formazione di PBH) e singolarità nude globali.
• Simulazioni di Fluido di Polvere (Non Sferiche):
  • Per il collasso triassico, le simulazioni del fluido di polvere falliscono a causa delle singolarità di incrocio dei gusci prima che possa formarsi un orizzonte apparente, a meno che l'ampiezza iniziale $\mu$ non sia significativamente elevata.
  • L'ampiezza soglia necessaria per osservare la formazione dell'orizzonte prima del crash del codice è risultata essere circa il doppio del valore previsto dalle precedenti stime analitiche basate sull'approssimazione di Zel'dovich e sulla congettura dell'anello (hoop conjecture).
  • Ciò suggerisce che, entro l'intervallo di validità dell'approssimazione del fluido, la formazione di PBH è più difficile da ottenere rispetto a quanto precedentemente ipotizzato.
• Simulazioni di Particelle Collisionless:
  • Modellando lo stesso sistema con particelle collisionless, le simulazioni non subiscono crash a causa dell'incrocio dei gusci.
  • La formazione di buchi neri (indicata dalla formazione dell'orizzonte apparente) è osservata per ampiezze iniziali così basse come $\mu \approx 0.045$ (per ellitticità $e=0.2$).
  • Questa soglia è un ordine di grandezza inferiore rispetto alle precedenti stime analitiche ed è significativamente più bassa della soglia derivata dalle simulazioni del fluido di polvere.
  • Tuttavia, le simulazioni delle particelle mostrano significative violazioni del vincolo hamiltoniano dopo l'intersezione delle particelle, indicando che, sebbene il collasso proceda, l'accuratezza quantitativa dell'evoluzione tardiva è limitata dal disallineamento tra la descrizione delle particelle e la geometria.

Significatività e Rivendicazioni
Il documento sostiene di dimostrare una grande incertezza teorica nel criterio di formazione dei PBH durante un'epoca dominata dalla materia, dipendente dalla descrizione fisica del campo di materia:

1. Limite dell'Approssimazione del Fluido: Se la formazione di singolarità associate all'incrocio dei gusci in una descrizione di fluido viene interpretata come un meccanismo di soppressione fisica (o se l'approssimazione del fluido decade prima del collasso), la soglia per la formazione di PBH è significativamente più alta delle previsioni analitiche, rendendo la formazione di PBH più difficile.
2. Dinamica delle Particelle: Se la descrizione delle particelle collisionless è accettata come un'estensione valida oltre l'approssimazione del fluido a singolo flusso, la formazione di PBH diventa molto più facile (soglia più bassa) rispetto a quanto precedentemente previsto.
3. Incertezza Teorica: La discrepanza tra i risultati del fluido e delle particelle evidenzia come i criteri analitici standard (approssimazione di Zel'dovich) possano non catturare accuratamente la dinamica non lineare del collasso triassico. Gli autori concludono che il valore soglia è altamente sensibile al trattamento del campo di materia ad alte densità.

Il lavoro non pretende di risolvere definitivamente l'abbondanza di PBH, ma piuttosto di esporre la sensibilità delle soglie di formazione alla descrizione della materia e al trattamento numerico delle singolarità. Si identifica come necessario un lavoro futuro per raffinare i trattamenti ad alta densità, migliorare la preservazione dei vincoli nelle simulazioni delle particelle e quantificare la risultante abbondanza di PBH e le firme di onde gravitazionali.