Using Strong Lensing to Detect Subhalos with Steep Inner Density Profiles

Probleemstelling

Het koude donkere materie (CDM)-model voorspelt dat subhalo's van donkere materie een universeel Navarro-Frenk-White (NFW) dichtheidsprofiel hebben, gekenmerkt door een binnenste helling ($\beta$) van 1 ($\rho \propto r^{-1}$). Echter, alternatieve donkere materie-modellen, zoals Zelf-interagerende Donkere Materie (SIDM), voorspellen een veel grotere diversiteit in binnenste dichtheidsprofielen. SIDM kan leiden tot zowel "gecorede" profielen ($\beta \sim 0$, door kernuitbreiding) als zeer steile profielen ($\beta > 2$, door gravitationele kerninstorting).

De huidige uitdaging is het testen van deze modellen op sub-galactische schalen. Sterk gravitationeel lenzen is een krachtige tool om lage-massa subhalo's te detecteren die te donker zijn om direct te zien. Echter, de detecteerbaarheid van deze subhalo's hangt sterk af van hun binnenste dichtheidsprofiel. Bestaande studies nemen vaak alleen NFW-profielen aan, wat de mogelijkheid om alternatieve donkere materie-modellen te onderscheiden beperkt. Het is onduidelijk hoe de detecteerbaarheid verandert als subhalo's steilere binnenste hellingen hebben, en hoe dit wordt beïnvloed door de complexiteit van het lensmodel (bijv. multipool-storingen) en de kwaliteit van de data.

Methodologie

De auteurs hebben een uitgebreide reeks van "mock" (gesimuleerde) waarnemingen gegenereerd met de open-source software PyAutoLens om de detecteerbaarheid van subhalo's te analyseren onder verschillende omstandigheden.

1. Simulatieset:

   • Ze simuleerden een "macro"-systeem bestaande uit een lensgalaxie (roodverschil $z=0.2$) en een brongalaxie ($z=1.0$) in een kruis-vierling (quad) configuratie.
   • Ze introduceerden één subhalo-stoorder binnen de lensgalaxie met een generalized NFW (gNFW) dichtheidsprofiel.
   • Drie verschillende binnenste hellingen ($\beta$) werden getest als proxies voor SIDM-fysica:
     • Cored: $\beta = 0.2$ (kernuitbreiding).
     • NFW: $\beta = 1.0$ (standaard CDM).
     • Steep: $\beta = 2.2$ (gravitationele kerninstorting).
   • De massa van de subhalo's varieerde logaritmisch tussen $10^{7.5}$ en $10^{10} M_\odot$.

2. Observatie Suites:
   De simulaties werden uitgevoerd voor vijf verschillende suites die variëren in instrumentele eigenschappen en subhalo-positie:

   • HST On-Ring: Gesimuleerd als Hubble Space Telescope-data, met de subhalo op de Einstein-ring (maximale interactie met de bron).
   • HST Offset: Dezelfde data, maar de subhalo staat ver van de Einstein-ring.
   • HST On-Ring (Higher SNR): HST-data met een verdubbelde signaal-ruisverhouding (SNR).
   • Euclid On-Ring: Gesimuleerd als Euclid-ruimtevaartuig-data (lagere resolutie).
   • JWST On-Ring: Gesimuleerd als James Webb Space Telescope-data (hoge resolutie en SNR).

3. Modelleerprocedure:

   • Ze gebruikten een Bayesiaanse niet-lineaire zoektocht (met Nautilus en Dynesty) om lensmodellen te fitten.
   • Stap 1: Fitten van een "smooth" model (zonder subhalo) om de lens- en bronparameters te bepalen.
   • Stap 2: Voegen van een subhalo-component toe aan het model en berekenen van de Bayesiaanse bewijskracht ($\mathcal{Z}$).
   • Detectiecriteria: Een subhalo wordt als gedetecteerd beschouwd als de toename in logaritmische bewijskracht $\Delta \ln \mathcal{Z} > 50$ (Jeffreys' schaal voor een beslissende detectie).
   • Complexiteit: Ze testten zowel een "Fiducial" model (elliptisch power-law) als een complexer model met multipool-storingen (orde $m=1, 3, 4$) om te zien of macro-model complexiteit de detectie van subhalo's verstoort.
   • Validatie: Een subset werd opnieuw geanalyseerd met een pixelized source reconstruction (Voronoi-mesh) in plaats van parametrische profielen, om de robuustheid te testen tegen complexe bronmorfologieën.

Belangrijkste Bijdragen

1. Kwantificering van de invloed van $\beta$: Het paper demonstreert dat de binnenste helling $\beta$ de meest kritische parameter is voor de detecteerbaarheid van subhalo's, vaak belangrijker dan de massa of concentratie.
2. Robuustheid tegen Macro-model Complexiteit: Een unieke bevinding is dat subhalo's met steile profielen ($\beta \sim 2.2$) veel minder gevoelig zijn voor degeneratie met multipool-storingen in het lensmodel dan NFW- of cored-subhalo's.
3. Instrumentele Vergelijking: Een systematische vergelijking van de potentie van HST, Euclid en JWST voor het detecteren van substructuren met verschillende dichtheidsprofielen.

Resultaten

• Detecteerbaarheid en Massa-drempel:

  • Steep subhalo's ($\beta=2.2$): Zijn aanzienlijk makkelijker te detecteren. In de HST On-Ring suite zijn ze detecteerbaar tot een massa van $\sim 5.4 \times 10^8 M_\odot$.
  • NFW subhalo's ($\beta=1.0$): Vereisen een massa van $\sim 6.0 \times 10^9 M_\odot$ voor detectie (meer dan een orde van grootte zwaarder dan steile profielen).
  • Cored subhalo's ($\beta=0.2$): Blijven ondetecteerbaar tot ten minste $10^{10} M_\odot$ in HST-achtige data.
  • JWST: Drukt de detectiedrempel voor steile subhalo's verder naar beneden tot $\sim 8.3 \times 10^7 M_\odot$.

• Invloed van Positie:

  • De superioriteit van steile profielen is het grootst wanneer de subhalo zich op de Einstein-ring bevindt, waar de bronemissie direct wordt verstoord.
  • Als de subhalo ver van de ring staat ("Offset"), neemt het verschil in detecteerbaarheid tussen steile en NFW-profielen af, omdat de afbuigingssterkte op grotere afstanden minder afhankelijk is van de binnenste helling.

• Invloed van Multipool-storingen:

  • Het toevoegen van multipool-storingen aan het lensmodel (om de complexiteit van de lensgalaxie te modelleren) maakt NFW- en cored-subhalo's vaak ondetecteerbaar (ze worden verward met macro-model variaties).
  • Cruciaal: Steep subhalo's op de Einstein-ring blijven sterk detecteerbaar ($\Delta \ln \mathcal{Z} > 50$) zelfs met multipool-storingen. Hun signatuur is uniek en kan niet worden nagebootst door standaard multipool-perturbaties.

• Parameterherwinning:

  • Voor detecteerbare steile subhalo's kan de binnenste helling $\beta$ robuust worden gereconstrueerd, zelfs bij lage massa's.
  • Er bestaat een degeneratie tussen concentratie en $\beta$, maar $\beta$ is over het algemeen beter beperkt dan de concentratie.

• Pixelized Source Reconstruction:

  • Het gebruik van een flexibeler, pixelized bronmodel verlaagt de detectiekracht ($\Delta \ln \mathcal{Z}$) voor alle subhalo's, omdat subhalo-effecten kunnen worden verward met morfologische variaties van de bron.
  • Desondanks blijven steep subhalo's significant detecteerbaarder dan NFW-subhalo's, en kunnen hun eigenschappen nog steeds robuust worden hersteld.

Significantie

De resultaten van dit onderzoek hebben belangrijke implicaties voor het testen van donkere materie-modellen:

1. Testen van SIDM: Omdat SIDM voorspelt dat sommige subhalo's een steile binnenste helling hebben (door kerninstorting), suggereert dit dat huidige en toekomstige waarnemingen (HST, Euclid, JWST) gevoelig genoeg zijn om deze specifieke subhalo's te detecteren, zelfs als ze lichter zijn dan de huidige detectielimieten voor CDM-subhalo's.
2. Onderscheid van Systematiek: Een van de grootste obstakels in sterk lensen is het onderscheiden van echte subhalo's van ongemodelleerde complexiteit in de lensgalaxie. Dit paper toont aan dat steile subhalo's een "veilig" signaal hebben dat minder vatbaar is voor deze systematische fouten.
3. Toekomstperspectief: Met de verwachte ontdekking van tienduizenden sterke lenzen door Euclid en de hoge resolutie van JWST, biedt deze studie een nieuwe route om de micro-fysica van donkere materie te beperken. Het suggereert dat het zoeken naar subhalo's met steile profielen een krachtige manier kan zijn om CDM te onderscheiden van alternatieve modellen zoals SIDM.

Kortom, dit werk benadrukt dat de binnenste dichtheidsstructuur van donkere materie subhalo's niet slechts een theoretisch detail is, maar een fundamentele factor die bepaalt of we ze überhaupt kunnen zien, en dat steile profielen een unieke, robuuste signatuur bieden voor toekomstige waarnemingen.