Consistent extinction model for type Ia supernovae in Cepheid-based calibration galaxies and its impact on $H_{0}$

Technische Samenvatting: Consistent Extinctiemodel voor Type Ia Supernovae in Cepheïde-gekalibreerde Sterrenstelsels en de Impact op $H_0$

Probleemstelling
De huidige spanning in de kosmologie komt voort uit een significante discrepantie ($5,2\sigma$) tussen de Hubble-constante ($H_0$) gemeten via de lokale afstandsladder (met behulp van Cepheïden en Type Ia supernovae in het SH0ES-programma) en de waarde die is afgeleid van de kosmische achtergrondstraling (Planck) onder aanname van een standaard $\Lambda$CDM-model. Een cruciaal onderdeel van de SH0ES-meting betreft de correctie van de magnitudes van Type Ia supernovae voor de extinctie door stof in het gaststelsel. De standaardbenadering (Popovic et al. 2023, hierna P23) hanteert een probabilistisch stofmodel dat is getraind op Hubble-flow supernovae ($z > 0,03$) en extrapoleert dit naar kalibratiestelsels (die met geobserveerde Cepheïden).

De auteurs identificeren een systematische inconsistentie in deze extrapolatie. De kalibratieset is vertekend richting late-type, stofrijke stelsels, terwijl het P23-model verschillende extinctie-eigenschappen toekent op basis van de stellaire massa ($M_\star$). Specifiek gaat het P23-model uit van een lage totale-tot-selectieve extinctiecoëfficiënt ($R_B \approx 3,1$) voor gasrijke gasten met een hoge massa ($M_\star > 10^{10} M_\odot$) en een standaardwaarde ($R_B \approx 4,0$) voor gasten met een lage massa. De auteurs stellen dat deze aanname leidt tot een onderschatting van de extinctie voor rode supernovae in hoogmassieve kalibratiestelsels, wat resulteert in systematisch zwakkere afgeleide absolute magnitudes en, bijgevolg, een overschatte $H_0$. Bovendien conflicteert de $R_B \approx 3,1$ van het P23-model voor hoogmassieve gasten met de Melkweg-achtige extinctiecurve ($R_B \approx 4,3$) die wordt gebruikt om de kleuren van Cepheïden in diezelfde stelsels te corrigeren.

Methodologie
De auteurs heranalyseren de Pantheon+-supernovacompilatie en de SH0ES Cepheïde-afstandmoduli (Riess et al. 2022). Hun methodologie omvat:

1. Consistentietesten: Ze voeren likelihood-analyses uit op de kalibratieset, waarbij supernovae worden gescheiden in hoogmassieve ($M_\star > 10^{10} M_\odot$) en laagmassieve categorieën. Ze vergelijken de best-fit absolute magnitude ($M_B$) en de daaruit afgeleide $H_0$ voor deze subgroepen tegenover willekeurige controlegroepen. Ze stellen vast dat hoogmassieve gasten systematisch zwakkere $M_B$ (hogere $H_0$) opleveren vergeleken met laagmassieve gasten en controlegroepen, met name voor rode supernovae ($c > 0$).
2. Modelmodificatie: Ze stellen een "minimalistische" modificatie voor van het extinctiemodel die uitsluitend wordt toegepast op de kalibratiestelsels, terwijl het P23-model behouden blijft voor de Hubble-flow (waar het model voor is getraind). De nieuwe methode behelst twee belangrijke wijzigingen:
   • Uniforme $R_B$-distributie: In plaats van een massafunctie-afhankelijke $R_B$, nemen ze een enkele, Melkweg-achtige Gaussische distributie aan voor de totale-tot-selectieve extinctiecoëfficiënt in beide massacategorieën, met een gemiddelde $\langle R_B \rangle = 4,3$ en een spreiding $\sigma_{R_B} = 0,4$. Dit stemt de extinctiecorrectie van de supernovae af op de curve die voor Cepheïden wordt gebruikt.
   • Gewijzigde vorm van de roodingsdistributie: Om de effectieve helling ($\beta \approx 3,0$) van de supernova peak magnitude–kleurrelatie en de gemiddelde roding $\langle E(B-V) \rangle$ gemeten in de Hubble-flow te behouden, vervangen ze de exponentiële distributie van stofroding $E(B-V)$ (gebruikt in P23) door een Gamma-distributie. De vormparameter is afgestemd op $\gamma = 3,44$, wat het piekpunt van de distributie naar een niet-nul roding verschuift terwijl het juiste gemiddelde behouden blijft.
3. Statistische Evaluatie: Ze berekenen de bias-correcties ($\delta$) en de vloer-onzekerheden ($\sigma_{floor}$) voor de kalibratieset opnieuw met behulp van het nieuwe model. Ze actualiseren de covariantie-matrix om de verminderde intrinsieke spreiding te reflecteren als gevolg van de lagere $\sigma_{R_B}$. Ten slotte voeren ze gezamenlijke fits uit van de kalibratie- en Hubble-flowgegevens om een nieuwe $H_0$ af te leiden.

Belangrijkste Bijdragen en Resultaten

• Identificatie van Selectiebias: De studie toont aan dat de aanname van het P23-model van $R_B \approx 3,1$ voor hoogmassieve kalibratiestelsels een $2,0\sigma$ tot $2,3\sigma$ spanning creëert in de afgeleide absolute magnitudes tussen hoogmassieve en laagmassieve gasten. Deze bias is het sterkst bij rode supernovae.
• Verbeterde Model Fit: Het voorgestelde nieuwe extinctiemodel levert een aanzienlijk betere fit aan de kalibratiedata op, met een verbetering in de Bayesian Information Criterion van $\Delta \text{BIC} = -11,0$ vergeleken met het P23-model. Deze verbetering wordt gedreven door zowel een reductie in $\chi^2_{min}$ (door gecorrigeerde biases) als een verandering in de likelihood-normalisatie (door verminderde covariantie-onzekerheden).
• Herziene Hubble-constante: Het toepassen van het nieuwe model op de gezamenlijke kalibratie- en Hubble-flowdataset resulteert in een lagere best-fit Hubble-constante:
  $$H_0 = 70,5 \pm 1,0 \text{ km s}^{-1} \text{ Mpc}^{-1}$$
  Dit vertegenwoordigt een reductie van ongeveer $2,9 \text{ km s}^{-1} \text{ Mpc}^{-1}$ vergeleken met de standaard SH0ES-waarde ($73,4 \pm 1,0$).
• Reductie van de Spanning: De herziene $H_0$ vermindert de statistische spanning met de Planck $\Lambda$CDM-meting van $5,2\sigma$ naar $2,8\sigma$. Het resultaat is consistent met onafhankelijke metingen gebaseerd op de Tip of the Red Giant Branch (TRGB).

Betekenis en Claims
Het artikel beweert dat de "Hubble-spanning" deels of grotendeels wordt gedreven door niet-verantwoorde systematische fouten in de extinctiemodellering van kalibratiestelsels. De auteurs stellen dat de standaardaanpak faalt omdat deze een model extrapoleert dat is getraind op een gemengde populatie van stelsels naar een kalibratieset die uitsluitend bestaat uit late-type, stervormende stelsels.

De betekenis van hun werk ligt in het aantonen dat:

1. Selectiebias ertoe doet: De specifieke selectie van kalibratiestelsels (die observeerbare Cepheïden vereisen) introduce wen een bias die niet simpelweg wordt gevangen door louter "late-type" stelsels te selecteren in de Hubble-flow, aangezien die laatste selectie de lokale stofomgevingen van de kalibratiesupernovae mogelijk niet nauwkeurig evenaart.
2. Extinctiefysica: De aanname van een lage $R_B$ ($\approx 3,1$) voor hoogmassieve gasten is fysiek inconsistent met de Melkweg-achtige extinctie die vereist is voor Cepheïden in diezelfde stelsels en met observationele beperkingen op stervormende stelsels.
3. Tweede-orde Effecten: De vorm van de rodingsdistributie (specifiek het tweede moment) is cruciaal. De auteurs tonen aan dat een Gamma-distributie ($\gamma \approx 3,4$) de data beter beschrijft dan het standaard exponentiële model, wat impliceert dat er een sterkere correlatie is tussen stof en de positie van supernovae in deze specifieke stelsels.

De auteurs concluderen dat hun herziene extinctiemodel een meer consistente fysische basis biedt voor de afstandsladder, waardoor de schijnbare discrepantie tussen hoog- en laagmassieve gasten wordt verkleind en de $H_0$-spanning aanzienlijk wordt verlaagd zonder nieuwe fysica aan te roepen. Ze merken op dat hoewel hun model sterk wordt ondersteund door de huidige data, verdere verfijning met behulp van onafhankelijke nabij-infraroodobservaties en volledige forward-modellering van stofeigenschappen noodzakelijk is om de resterende residuen volledig op te lossen.