Core-halo scaling relations in self-interacting scalar field dark matter

Enunciato del Problema
La Materia Oscura a Campo Scalare (SFDM), spesso indicata come Materia Oscura "Fuzzy" (FDM) quando le auto-interazioni sono trascurate, postula che la materia oscura sia composta da bosoni ultraleggeri ($m \sim 10^{-22}$ eV/$c^2$). Sebbene i modelli FDM standard spieghino con successo la formazione di nuclei solitonici negli aloni galattici, l'impatto delle auto-interazioni (SI) sulla struttura e sull'evoluzione di tali aloni rimane una questione aperta. La letteratura precedente ha stabilito relazioni di scala per la FDM non interagenti (ad esempio, relazioni massa-nucleo), ma non è chiaro come le auto-interazioni repulsive o attrattive modifichino queste relazioni, in particolare per quanto riguarda la stabilità dei nuclei, la relazione massa-raggio e il potenziale di collasso gravitazionale in buchi neri supermassicci. Inoltre, sebbene il regime di Thomas-Fermi (TF) della SFDM fortemente repulsiva sia stato studiato utilizzando approssimazioni fluide, le simulazioni numeriche dirette di fusioni di solitoni in questo regime senza approssimazioni fluide sono state scarse.

Metodologia
Gli autori eseguono simulazioni numeriche tridimensionali di fusioni multiple di solitoni per indagare gli aloni SFDM su un'ampia gamma di intensità e segni delle auto-interazioni.

• Equazioni Governative: Lo studio risolve le equazioni non relativistiche di Gross-Pitaevskii-Poisson (GPP), che descrivono l'evoluzione della funzione d'onda del campo scalare $\psi$ accoppiata alla gravità. Le equazioni includono un termine di auto-interazione cubico (proporzionale alla lunghezza di scattering $a_s$) e un termine di ordine superiore quintico.
• Schema Numerico: Gli autori utilizzano un metodo pseudo-spettrale con splitting temporale ottimizzato per GPU. Questo approccio tratta i termini di auto-interazione come un potenziale efficace, permettendo l'uso di uno schema kick-drift-kick (KDK). Per le interazioni attrattive, viene implementato un metodo KDK di ordine superiore (O(4)) per garantire la stabilità.
• Configurazione della Simulazione: Le simulazioni coinvolgono la fusione di $N_{sol}$ solitoni (da 10 a 50) all'interno di una scatola di simulazione statica e periodica di $L=100$ kpc. La massa della particella è fissata a $m = 10^{-22}$ eV/$c^2$.
• Spazio dei Parametri: Lo studio esplora tre regimi distinti:
  1. Libero (FDM): $a_s = 0$.
  2. SI Repulsiva: $a_s > 0$, incluso il regime di Thomas-Fermi fortemente repulsivo ($\Lambda \gg 1$).
  3. SI Attrattiva: $a_s < 0$, fino a valori in cui il sistema si avvicina al collasso gravitazionale.
• Analisi: Gli autori analizzano i profili di densità, la massa del nucleo ($M_c$), il raggio del nucleo ($r_c$) e l'energia totale ($E$) per derivare relazioni di scala. Confrontano i risultati delle simulazioni con le previsioni analitiche derivate dalla minimizzazione dell'energia e dalle soluzioni autosimili.

Contributi e Risultati Chiave

1. Struttura e Densità del Nucleo:

   • SI Repulsiva: Porta alla formazione di nuclei più massicci ed estesi con densità centrali più basse rispetto al caso FDM libero. L'interazione repulsiva fornisce un supporto di pressione aggiuntivo, impedendo al nucleo di diventare compatto come nel caso non interagenti.
   • SI Attrattiva: Aumenta le densità centrali. Per interazioni attrattive sufficientemente forti, la densità del nucleo diverge, indicando un collasso gravitazionale. Gli autori identificano una scala di massa critica ($M_{max}$) per i modelli attrattivi; se la massa del nucleo supera questa soglia, l'equilibrio viene perso, portando potenzialmente alla formazione di buchi neri supermassicci.
   • Evoluzione: I nuclei repulsivi rimangono stabili ed estesi per tutta la durata della simulazione ($t = 10\tau_{dyn}$), mentre i nuclei attrattivi mostrano un aumento monotono della densità centrale fino al collasso.

2. Relazioni di Scala Massa-Raggio e Nucleo-Alone:

   • Relazione Massa-Raggio: I risultati delle simulazioni per la relazione massa-raggio del nucleo ($M_c$ vs. $r_c$) si allineano bene con le previsioni analitiche (Chavanis 2011; Padilla et al. 2021) anche in presenza di forti auto-interazioni. Le interazioni repulsive spostano la relazione verso l'alto (permettendo più massa per un dato raggio), mentre le interazioni attrattive la spostano verso il basso.
   • Relazione Massa-Nucleo/Alone: Gli autori indagano la relazione tra il rapporto massa-nucleo/massa-alone ($M_c/M$) e la quantità invariante $\Xi = |E| / [M^3(Gm/\hbar^2)]$. Scoprono che, sebbene la pendenza di questa relazione rimanga relativamente universale (simile al caso FDM $M_c/M \propto \Xi^{1/3}$), la normalizzazione dipende dall'intensità dell'auto-interazione. Ciò suggerisce che il meccanismo che alimenta il nucleo di massa è principalmente gravitazionale e universale una volta che gli aloni sono virializzati, mentre le SI determinano localmente il supporto di pressione e il raggio del nucleo.
   • Relazione Dimensione-Energia: La relazione tra la dimensione del nucleo ($r_c$) e l'energia totale ($E$) è sensibile alle SI. Le interazioni repulsive risultano in pendenze più dolci (nuclei più grandi per una data energia), mentre le interazioni attrattive risultano in pendenze più ripide.

3. Decomposizione dell'Energia:

   • Nei casi repulsivi, l'energia di auto-interazione domina la regione del nucleo, spiegando la minore densità e la maggiore dimensione.
   • Nei casi attrattivi e liberi, l'energia potenziale gravitazionale domina il nucleo.
   • Gli autori confermano che le relazioni di scala non sono strettamente universali, ma dipendono dallo stadio evolutivo dell'alone e dai specifici parametri delle SI.

Significato e Affermazioni
Il lavoro afferma che le auto-interazioni forniscono un meccanismo naturale per regolare le proprietà dei nuclei degli aloni SFDM. Lo studio dimostra che:

• Le relazioni di scala nucleo-alone derivate per la FDM libera sono casi limite all'interno di una famiglia più ampia di modelli SFDM che includono auto-interazioni.
• Le interazioni repulsive possono stabilizzare i nuclei contro il collasso e produrre strutture più grandi e diffuse, il che potrebbe avere implicazioni per il problema "nucleo-cuspide" nella formazione galattica.
• Le interazioni attrattive possono spingere i nuclei verso l'instabilità e il collasso, offrendo un potenziale canale di formazione per i semi di buchi neri supermassicci nei centri galattici.
• Il metodo pseudo-spettrale è validato come strumento robusto per simulare la SFDM nel regime di Thomas-Fermi senza ricorrere ad approssimazioni fluide, catturando con successo la transizione dal comportamento solitonico a quello fluido.

Gli autori concludono che, sebbene la dinamica gravitazionale fondamentale della crescita del nucleo sia universale, le proprietà strutturali specifiche (densità, dimensione, stabilità) sono altamente sensibili al segno e all'entità delle auto-interazioni e alla storia evolutiva dell'alone. Questi risultati suggeriscono che i vincoli osservativi sui nuclei galattici potrebbero essere utilizzati per vincolare la lunghezza di scattering e la massa della particella della SFDM.