The Sensitivity of Substructure Lensing to SIDM Core-collapse Model Variation

Technische Samenvatting: De Gevoeligheid van Substructuurlensings voor Variatie in SIDM-kerninstortingsmodellen

Probleemstelling
Sterke gravitationele lensing is een krachtige sonde voor donkere-materiesubstructuur geworden, die een specifieke test biedt voor modellen van Zelf-interagerende Donkere Materie (SIDM). In SIDM-scenario's met elastische botsingen leidt warmtetransport tot de vorming van isotherme kernen, gevolgd door een "gravothermale catastrofe" (kerninstorting) die dichte centrale kernen produceert. Deze compacte halos voorspellen dat ze duidelijke afdrukken achterlaten op lensingsobservabelen. Echter, een accurate voorspelling van deze signatuur wordt gehinderd door de moeilijkheid om de tijdslijn van kerninstorting en de eindtoestand van subhalos te modelleren. Huidige benaderingen vertrouwen vaak op vereenvoudigende aannames, zoals het behandelen van kerninstorting als een instantane overgang of het gebruik van probabilistische toewijzingen voor instortingstoestanden. De auteurs identificeren een disconnectie tussen bestaande lensingsinferentiepijplijnen en de onderliggende fysische modellen, en merken op dat kleine variaties in de modellering van kerninstortende subhalos kunnen leiden tot aanzienlijke verschillen in lensingskracht. Zonder kwantificering van deze modelleringseonzekerheden is het vermogen van lensingsobservaties om de aard van donkere materie te beperken, beperkt.

Methodologie
De auteurs kwantificeren deze modelleringstevredenheid door verschillende modificaties van SIDM-dichtheidsprofielen van kerninstortende halos te testen binnen de context van substructuurlensing. Hun methodologie omvat:

1. Simulatie en Kalibratie: Zij gebruiken $N$-lichaamssimulatiegegevens (gegenereerd met de code Arepo en de SIDM-module uit Ref. [20]) om een benaderde behandeling van getij-effecten op gravothermale kerninstorting te kalibreren. Zij leiden een "instortingsversnellingsfactor" af om het verschil in instortingstijdschalen tussen geïsoleerde halos en subhalos die evolueren in een gastheerpotentiaal, in aanmerking te nemen.
2. Lensing-realiteiten: Met behulp van het open-source pakket pyHalo genereren zij 200 willekeurige realisaties van lenssystemen op galaxieschaal (bij $z_{lens}=0.5$) die subhalos van de hoofddeflector en halos langs de zichtlijn (LOS) bevatten. Zij vervangen standaard Koude Donkere Materie (CDM)-profielen door SIDM-profielen gebaseerd op een gravothermale vloeistofoplossing die is gekalibreerd aan $N$-lichaamssimulaties.
3. Modelvariaties: De studie test verschillende modelleringkeuzes:
   • Continue versus Instantaan: Het vergelijken van glad evoluerende dichtheidsprofielen met binaire modellen waarbij halos óf volledig niet-geinstort óf volledig geinstort zijn.
   • Instortingsdrempels: Variëren van het tijdstip waarop het dichtheidsprofiel wordt "bevroren" (gestopt) of "geboost" naar de eindtoestand, en het testen van verschillende eindtoestanden van de instorting.
   • Probabilistisch versus Deterministisch: Het vergelijken van individuele halo-evolutietracking tegen willekeurige toewijzing van instortingstoestanden op basis van massabins.
   • Initiële Profielen: Het testen van niet-geinstorte halos die gemodelleerd zijn als óf standaard afgeknipte NFW (TNFW) óf als halos die de maximaal gekerde fase hebben bereikt.
4. Statistische Analyse: Zij berekenen het effectieve meerlagige gravitationele lensingspotentiaal, en berekenen de divergentie ($\kappa_{div}$) en de rotatie ($\kappa_{curl}$) van het effectieve afbuigingsveld. De primaire statistiek is de tweepuntscorrelatiefunctie van deze velden, ontbonden in multipoolmomenten (monopool en quadrupool) met behulp van orthonormale bases om LOS-effecten van de hoofdlenst te ontwarren.

Belangrijkste Bijdragen

• Kwantificering van Modelleringseonzekerheid: Het artikel biedt een systematische kwantificering van hoe specifieke vereenvoudigende aannames in SIDM-kerninstortingsmodellering lensingsstatistieken beïnvloeden.
• Getijkalibratie: Het introduceert een probabilistische methode om versnelling/vertraging van kerninstorting in subhalos door getijdenkrachten in aanmerking te nemen, gekalibreerd tegen $N$-lichaamssimulaties.
• Vergelijking Binaire versus Continue Modellen: Het vergelijkt rigoureus de "binaire instorting"-benadering (veelvoorkomend in de literatuur) met een continue, tijd-evoluerend gravothermaal model, en wijst aan waar de benadering faalt.
• Gevoeligheid van Profielen: De studie toont aan dat de keuze van het initiële dichtheidsprofiel (TNFW versus gekerd) en de methode van toewijzing van instortingstoestanden het voorspelde lensingsignaal aanzienlijk verandert.

Resultaten

• Gevoeligheid voor Dichtheid versus Tijd: De tweepuntscorrelatiefunctie (specifiek de multipoolmomenten) is zeer gevoelig voor de centrale dichtheid van de geinstorte halos, maar relatief ongevoelig voor de exacte instortingstijdschaal. Variaties in de instortingstijd (bijvoorbeeld het stoppen van evolutie op 20% van de volledige tijd) produceren niet significant verschillende resultaten vergeleken met de volledige evolutie, mits het uiteindelijke dichtheidsprofiel accuraat is.
• Onderdrukking door Binaire Modellen: Binaire kerninstortingsmodellen onderdrukken consistent de kracht op kleine schaal in lensingsstatistieken vergeleken met het volledige tijd-evoluerende model. Deze onderdrukking is het ernstigst wanneer:
  • Niet-geinstorte halos worden gemodelleerd met gekerde initiële profielen in plaats van TNFW-profielen.
  • Instortingstoestanden willekeurig worden toegewezen in plaats van gebaseerd op individuele halo-evolutie.
• Schaalafhankelijkheid: De verschillen tussen modellen zijn het meest onderscheidbaar op kleine lengteschalen (onder de 0,1 boogseconde). Het quadrupoolmoment van het divergentieveld ($\xi_{div,2}$) blijft gevoelig voor modelvariaties op grotere schalen dan het monopool- of rotatiemoment.
• Correlatie van Instorting en Lensingskracht: De studie vindt dat halos met een meer uitgesproken lensingseffect (hogere massa/concentratie) sneller instorten. Het randomiseren van instortingstoestanden onderdrukt de SIDM-versterking omdat het de instortingstoestand ontkoppelt van de halo-eigenschappen die het lensingsignaal aandrijven.
• N-lichaam Validatie: Het benaderde getijversnellingsmodel van de auteurs produceert lensingsstatistieken die consistent zijn met $N$-lichaamssimulaties binnen één standaardafwijking, wat de bruikbaarheid ervan als eerste-orde benadering voor deze studie valideert.

Betekenis en Claims
Het artikel claimt dat de tweepuntscorrelatiefunctie een gevoelige sonde is van de specifieke vereenvoudigingen die worden gemaakt in SIDM-kerninstortingsmodellering. De auteurs stellen het volgende:

1. Robuustheid van de Eindtoestand: Matige systematische onzekerheden in de timing van kerninstorting kunnen lensingsvoorspellingen niet significant beïnvloeden, mits de einddichtheidsprofiel van de geinstorte halo accuraat wordt gemodelleerd.
2. Beperkingen van Binaire Benaderingen: Hoewel binaire instortingsmodellen computatie-efficiënt zijn, kunnen ze aanzienlijke bias introduceren door onderdrukking van kracht op kleine schaal, vooral als ze gekerde initiële condities aannemen of instortingstoestanden randomiseren.
3. Behoefte aan Fysische Motivatie: Toekomstige lensingsanalyses die SIDM willen testen, moeten intentionele, fysisch gemotiveerde keuzes maken met betrekking tot substructuur- en LOS-haloprofielen. De auteurs raden aan dat dergelijke werken vergelijkbare tests van hun profielkeuzes opnemen en hun afbuigingshoekprofielen vergelijken met fysisch gemotiveerde, volledig evoluerende modellen.
4. Observatie-afhankelijkheid: De detecteerbaarheid van deze modelleringverschillen hangt sterk af van de observationele sonde en de schalen die worden geanalyseerd. Het artikel claimt niet dat huidige data specifieke SIDM-modellen definitief kan uitsluiten, maar wel dat de interpretatie van dergelijke data sterk afhankelijk is van de robuustheid van de onderliggende substructuurmodellen.

De auteurs concluderen dat het verfijnen van modellen van SIDM-kerninstorting in realistische omgevingen essentieel is, aangezien testbare voorspellingen sterk kunnen afhangen van de eindtoestand van kerngeinstorte halos. Zij suggereren dat hun kader kan worden gebruikt om deze volledig evoluerende modellen te implementeren in voorspellingen voor flux-ratio anomalieën en andere lensingsobservabelen in toekomstig werk.