Dynamical Systematics for Time Delay Lenses and the Impact on the Hubble Constant

Technische Samenvatting: Dynamische Systematiek voor Tijdvertraging-lenzen en de Impact op de Hubble-constante

Probleemstelling
Tijdvertraging-gravitatielenzen bieden een methode om de Hubble-constante ($H_0$) onafhankelijk van de kosmische afstandsladder te meten. Echter, deze metingen lijden onder een fundamentele degeneratie waarbij $H_0 \propto 1 - \kappa_E$, waarbij $\kappa_E$ de convergentie nabij de Einstein-straal is. Hoewel lensgegevens de Einstein-straal ($R_E$) en een dimensieloze grootheid $\xi$ gerelateerd aan de massaprofielhelling beperken, kunnen ze $\kappa_E$ niet bepalen zonder aanvullende aannames. Stellaire dynamische metingen, specifiek centrale snelheidsdispersies ($\sigma$), worden gebruikt om deze degeneratie te doorbreken. De nauwkeurigheid van de uit $H_0$ afgeleide waarde wordt beperkt door de precisie van de kinematische metingen; het bereiken van een precisie van 2% in $H_0$ vereist dat de systematische onzekerheden in de gekwadrateerde snelheidsdispersie ($\Delta\sigma^2/\sigma^2$) kleiner zijn dan 2%. Dit artikel onderzoekt het bereik van potentiële systematische onzekerheden die acht tijdvertraging-lenzen beïnvloeden die in kosmologische analyses worden gebruikt, en betoogt dat de huidige schattingen van $H_0$ vertekend kunnen zijn door onverantwoorde systematiek in de dynamische modellering.

Methodologie
De auteurs analyseren acht tijdvertraging-lenzen (DES J0408−5354, HE 0435−1223, PG 1115+080, RX J1131−1231, SDSS J1206+4332, B1608+656, WFI 2033−4723, en WGD 2038−4008) met behulp van een sferisch symmetrisch dynamisch model. De studie maakt gebruik van een Hernquist stellaire distributie ingebed in een Singular Isothermal Sphere (SIS) massamodel om systematische effecten te isoleren. De analyse richt zich op vier primaire categorieën van systematiek:

1. Point Spread Function (PSF) Modellering: De auteurs testen de gevoeligheid van snelheidsdispersie-metingen voor fouten in de PSF Full Width at Half Maximum (FWHM), apertuur-miscentrering, en de aanname van een Gaussische PSF versus meer realistische Moffat-profielen.
2. Gemeten versus Ware Dispersie: De studie maakt onderscheid tussen de gemeten snelheidsdispersie ($\sigma_*$, afgeleid van het fitten van een Gaussische curve aan spectrale lijnen) en de root-mean-square snelheid ($v_{rms}$) die vereist is voor de Jeans-vergelijkingen. De auteurs kwantificeren de bias geïntroduceerd door niet-Gaussische snelheidsdistributies, gekarakteriseerd door het Gauss-Hermite moment $h_4$.
3. Anisotropie-modellen: De auteurs vergelijken constante anisotropie-modellen, het Osipkov-Merritt (O-M) model, en alternatieve profielen (Mamon & Lokas; Cuddeford) om te bepalen hoe de keuze van orbitale anisotropie de afgeleide snelheidsdispersie en $H_0$ beïnvloedt.
4. Fotometrische en Stellaire Populatiemodellen: De analyse onderzoekt de impact van het gebruik van verschillende stellaire dichtheidsprofielen (Sérsic versus Hernquist versus Jaffe), de keuze van de effectieve straal-as (major versus intermediair), en de complicaties die voortvloeien uit kleurgradiënten en meerdere stellaire populaties in vroege-type sterrenstelsels.

Daarnaast behandelt het artikel de statistische behandeling van de lens-steekproef, waarbij wordt betoogd dat vroege-type sterrenstelsels een homogene populatie vormen wat betreft dynamica en stellaire populaties, wat impliceert dat bepaalde systematische fouten over de gehele steekproef gemarginaliseerd moeten worden in plaats van behandeld te worden als onafhankelijke willekeurige fouten voor elke lens afzonderlijk.

Belangrijkste Bijdragen en Resultaten
De auteurs stellen vast dat systematische onzekerheden in de dynamische modellering gemakkelijk de drempel van 2% kunnen overschrijden die vereist is voor competitieve $H_0$-metingen. De specifieke bevindingen zijn:

• PSF-Systematiek: Voor observaties vanaf de grond kunnen biases in de breedte en vorm van de PSF variaties in $\sigma^2$ van 1–5% induceren. Het onderschatten van de PSF-breedte of het aannemen van een Gaussisch profiel wanneer de ware PSF uitgebreide vleugels heeft (Moffat), leidt tot onderschatte snelheidsdispersies, wat modellen richting minder centraal geconcentreerde massadistributies en lagere $H_0$-waarden drijft. Apertuur-miscentrering heeft een kleiner effect (doorgaans <0,3%), maar vertekent de resultaten nog steeds naar beneden.
• Definitie van Dispersie: Het verschil tussen de gemeten $\sigma_*$ en de vereiste $v_{rms}$ is significant. Voor isotrope modellen varieert dit verschil tussen de ~1,8% en 5,6% over de steekproef. Indien vroege-type sterrenstelsels matig radiaal anisotroop zijn (zoals verwacht), is de ware $v_{rms}$ groter dan de gemeten $\sigma_*$, wat leidt tot onderschattingen van $H_0$ indien niet gecorrigeerd.
• Anisotropie-modellen: De keuze van het anisotropie-model heeft een diepgaande impact. Terwijl constante anisotropie-modellen en O-M modellen met grote breukradii binnen een smalle enveloppe vallen, kunnen modellen die matige radiale anisotropie in de centrale regio's toestaan (bijv. Cuddeford-modellen met $\beta_0 > 0$) de afgeleide snelheidsdispersie met 2–18% verhogen. Dit verschuift de massadistributie naar een minder geconcentreerde staat, wat $H_0$ potentieel naar boven kan vertekend als dit niet goed wordt geconstreind. Bovendien vereisen modellen die hogere $h_4$-waarden produceren (consistent met observaties van vroege-type sterrenstelsels) grotere correcties tussen $\sigma_*$ en $v_{rms}$.
• Fotometrische Profielen: Mismatches tussen het aangenomen stellaire massaprofiel en de ware distributie (bijv. het gebruik van een Sérsic-profiel met een incorrecte index) kunnen $\Delta\sigma^2/\sigma^2$-variaties van 1–10% veroorzaken. Daarnaast impliceert de aanwezigheid van kleurgradiënten dat de stellaire getalsdichtheid $n$ in de Jeans-vergelijking gewogen moet worden door de equivalentbreedte van spectrale lijnen in plaats van het breedbandige lichtprofiel, een factor die momenteel niet wordt meegenomen in standaardanalyses.
• Homogeniteit van de Populatie: Het artikel betoogt dat het behandelen van parameters zoals de anisotropieradius of de stellaire template-keuze als onafhankelijke willekeurige variabelen voor elke lens de totale onzekerheid onderschat. Omdat vroege-type sterrenstelsels een homogene populatie zijn, zijn systematische fouten in deze parameters gecorreleerd over de steekproef. Het marginaliseren over deze parameters op populatieniveau (in plaats van per lens) voorkomt de kunstmatige reductie van onzekerheid met een factor $N^{-1/2}$.

Betekenis en Claims
Het artikel claimt dat de huidige precisie van $H_0$-schattingen uit tijdvertraging-lenzen waarschijnlijk wordt beperkt door systematische onzekerheden in de dynamische modellering in plaats van door statistische ruis. De auteurs demonstreren dat biases in $\sigma^2$ van 1–18% mogelijk zijn, afhankelijk van de modelkeuzes, wat vergelijkbaar is met of groter dan de vereiste precisie van 2% voor $H_0$.

De auteurs benadrukken dat:

1. Absolute Kalibratie Cruciaal Is: In tegenstelling tot zuivere stellaire dynamische studies die afhankelijk zijn van ratio's, hangt lensing af van de absolute kalibratie van snelheidsdispersies, wat het zeer gevoelig maakt voor systematische fouten in PSF-modellering en de definitie van dispersie.
2. Modelafhankelijkheid: De keuze van het anisotropie-model en de behandeling van de snelheidverdelingsfunctie (specifiek de $h_4$-correctie) zijn kritiek. Huidige modellen slagen er vaak niet in de geobserveerde $h_4$-waarden van vroege-type sterrenstelsels te reproduceren, wat suggereert dat huidige correcties onderschattingen zijn.
3. Statistische Behandeling: De aanname van een heterogene lens-populatie leidt tot een onderschatting van de uiteindelijke onzekerheden. Een homogeen populatiemodel vereist het marginaliseren van systematische parameters over de gehele steekproef, wat de bijdrage van de systematische fout behoudt, ongeacht de omvang van de steekproef.

Het artikel concludeert dat het bereiken van een 2%-meting van $H_0$ vereist dat deze systematische onzekerheden, met name de absolute kalibratie van snelheidsdispersies en de correcte modellering van orbitale anisotropie en stellaire populaties, onder controle worden gebracht. Hoewel ruimelijk opgeloste kinematica (2D-data) een weg bieden om degeneraties te doorbreken, moeten de fundamentele systematische biases betreffende de definitie van snelheidsdispersie en de homogeniteit van de populatie nog steeds worden geadresseerd.