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

Probleemstelling
Scalar Field Dark Matter (SFDM), vaak aangeduid als Fuzzy Dark Matter (FDM) wanneer zelfinteracties worden verwaarloosd, stelt dat donkere materie bestaat uit ultralichte bosonen ($m \sim 10^{-22}$ eV/$c^2$). Hoewel standaard FDM-modellen met succes de vorming van solitonische kernen in galactische halos verklaren, blijft de impact van zelfinteracties (SI) op de structuur en evolutie van deze halos een open vraag. Eerdere literatuur heeft schaalrelaties vastgesteld voor niet-interagerende FDM (bijvoorbeeld kern-halo-massa-relaties), maar het is onduidelijk hoe repulsieve of attractieve zelfinteracties deze relaties wijzigen, met name wat betreft de stabiliteit van kernen, de massa-radiussrelatie en het potentieel voor gravitationele instorting tot supermassieve zwarte gaten. Bovendien, hoewel het Thomas-Fermi (TF)-regime van sterk repulsieve SFDM is bestudeerd met behulp van vloeistofbenaderingen, zijn directe numerieke simulaties van soliton-samensmeltingen in dit regime zonder vloeistofbenaderingen schaars.

Methodologie
De auteurs voeren driedimensionale numerieke simulaties uit van meerdere soliton-samensmeltingen om SFDM-halos te onderzoeken over een breed scala aan zelfinteractiestraken en -tekens.

• Sturende Vergelijkingen: De studie lost de niet-relativistische Gross-Pitaevskii-Poisson (GPP)-vergelijkingen op, die de evolutie van de scalarveldgolfunctie $\psi$ gekoppeld aan zwaartekracht beschrijven. De vergelijkingen omvatten een kubische zelfinteractieterm (evenredig met de verstrooiingslengte $a_s$) en een hogere-orde kwintische term.
• Numeriek Schema: De auteurs maken gebruik van een pseudo-spectrale tijdsplitsingsmethode die is geoptimaliseerd voor GPU's. Deze benadering behandelt de zelfinteractietermen als een effectief potentieel, waardoor het gebruik van een kick-drift-kick (KDK)-schema mogelijk wordt. Voor attractieve interacties wordt een hogere-orde (O(4)) KDK-methode geïmplementeerd om stabiliteit te waarborgen.
• Simulatieopstelling: Simulaties omvatten de samensmelting van $N_{sol}$ solitons (variërend van 10 tot 50) binnen een statische, periodieke simulatiebox van $L=100$ kpc. De deeltjesmassa is vastgezet op $m = 10^{-22}$ eV/$c^2$.
• Parameter Ruimte: De studie verkent drie distincte regimes:
  1. Vrij (FDM): $a_s = 0$.
  2. Repulsieve SI: $a_s > 0$, inclusief het sterk repulsieve Thomas-Fermi-regime ($\Lambda \gg 1$).
  3. Attractieve SI: $a_s < 0$, tot waarden waarbij het systeem gravitationele instorting nadert.
• Analyse: De auteurs analyseren dichtheidsprofielen, kernmassa ($M_c$), kernstraal ($r_c$) en totale energie ($E$) om schaalrelaties af te leiden. Ze vergelijken simulatieresultaten met analytische voorspellingen die zijn afgeleid uit energieminimalisatie en zelfgelijkende oplossingen.

Belangrijkste Bijdragen en Resultaten

1. Kernstructuur en Dichtheid:

   • Repulsieve SI: Leidt tot de vorming van massievere en uitgebreidere kernen met lagere centrale dichtheden in vergelijking met het vrije FDM-geval. De repulsieve interactie biedt extra drukondersteuning, waardoor de kern niet zo compact wordt als in het niet-interagerende geval.
   • Attractieve SI: Verhoogt de centrale dichtheden. Voor voldoende sterke attractieve interacties divergeert de kern dichtheid, wat wijst op gravitationele instorting. De auteurs identificeren een kritische massaschaal ($M_{max}$) voor attractieve modellen; als de kernmassa deze drempel overschrijdt, gaat het evenwicht verloren, wat potentieel kan leiden tot de vorming van supermassieve zwarte gaten.
   • Evolutie: Repulsieve kernen blijven stabiel en uitgebreid gedurende de simulatie ($t = 10\tau_{dyn}$), terwijl attractieve kernen een monotoon toenemende centrale dichtheid vertonen tot instorting.

2. Massa-straal- en Kern-halo Schaalrelaties:

   • Massa-straal Relatie: De simulatieresultaten voor de kern massa-straal relatie ($M_c$ vs. $r_c$) sluiten goed aan bij analytische voorspellingen (Chavanis 2011; Padilla et al. 2021), zelfs in aanwezigheid van sterke zelfinteracties. Repulsieve interacties verschuiven de relatie omhoog (meer massa toestaan voor een gegeven straal), terwijl attractieve interacties deze naar beneden verschuiven.
   • Kern-halo Massa Relatie: De auteurs onderzoeken de relatie tussen de kern-tot-halo-massa verhouding ($M_c/M$) en de invariant grootheid $\Xi = |E| / [M^3(Gm/\hbar^2)]$. Ze vinden dat de helling van deze relatie relatief universeel blijft (gelijk aan het FDM-geval $M_c/M \propto \Xi^{1/3}$), maar dat de normalisatie afhangt van de zelfinteractiestrakte. Dit suggereert dat het mechanisme dat massa aan de kern toevoert primair gravitationeel en universeel is zodra halos virialiseren, terwijl SI lokaal de drukondersteuning en straal van de kern bepaalt.
   • Grootte-Energie Relatie: De relatie tussen kerngrootte ($r_c$) en totale energie ($E$) is gevoelig voor SI. Repulsieve interacties resulteren in vlakkere hellingen (grotere kernen voor een gegeven energie), terwijl attractieve interacties resulteren in steilere hellingen.

3. Energie Decompositie:

   • In repulsieve gevallen domineert de zelfinteractie-energie het kerngebied, wat de lagere dichtheid en grotere grootte verklaart.
   • In attractieve en vrije gevallen domineert de gravitationele potentiële energie de kern.
   • De auteurs bevestigen dat schaalrelaties niet strikt universeel zijn, maar afhankelijk van het evolutiestadium van de halo en de specifieke SI-parameters.

Betekenis en Claims
Het artikel claimt dat zelfinteracties een natuurlijk mechanisme bieden om de eigenschappen van SFDM-halo-kernen te reguleren. De studie toont aan dat:

• De voor vrij FDM afgeleide kern-halo schaalrelaties limietgevallen zijn binnen een bredere familie van SFDM-modellen die zelfinteracties omvatten.
• Repulsieve interacties kernen kunnen stabiliseren tegen instorting en grotere, meer diffuse structuren kunnen produceren, wat implicaties kan hebben voor het "kern-cusp"-probleem in galactische vorming.
• Attractieve interacties kernen naar instabiliteit en instorting kunnen drijven, wat een potentieel vormingskanaal biedt voor zaden van supermassieve zwarte gaten in galactische centra.
• De pseudo-spectrale methode is gevalideerd als een robuust hulpmiddel voor het simuleren van SFDM in het Thomas-Fermi-regime zonder terug te grijpen op vloeistofbenaderingen, en slaagt erin de overgang van solitonisch naar vloeistofachtig gedrag succesvol vast te leggen.

De auteurs concluderen dat hoewel de fundamentele gravitationele dynamiek van kerngroei universeel is, de specifieke structurele eigenschappen (dichtheid, grootte, stabiliteit) zeer gevoelig zijn voor het teken en de grootte van zelfinteracties en de evolutiegeschiedenis van de halo. Deze bevindingen suggereren dat observationele beperkingen op galactische kernen kunnen worden gebruikt om de verstrooiingslengte en de deeltjesmassa van SFDM te beperken.