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

Technische Zusammenfassung: Konsistentes Extinktionsmodell für Typ Ia Supernovae in Cepheiden-basierten Kalibrierungsgalaxien und dessen Auswirkung auf $H_0$

Problemstellung
Die aktuelle Spannung in der Kosmologie ergibt sich aus einer signifikanten Diskrepanz ($5,2\sigma$) zwischen der Hubble-Konstante ($H_0$), die über die lokale Entfernungsleiter (unter Verwendung von Cepheiden und Typ Ia Supernovae im SH0ES-Programm) gemessen wurde, und dem Wert, der aus dem kosmischen Mikrowellenhintergrund (Planck) unter Annahme eines Standard-$\Lambda$CDM-Modells abgeleitet wurde. Eine kritische Komponente der SH0ES-Messung ist die Korrektur der Magnituden von Typ Ia Supernovae hinsichtlich der Wirtsgalaxien-Extinktion. Der Standardansatz (Popovic et al. 2023, im Folgenden P23) verwendet ein probabilistisches Staubmodell, das an Hubble-Flow-Supernovae ($z > 0,03$) trainiert wurde und auf Kalibrierungsgalaxien (die beobachtbare Cepheiden aufweisen) extrapoliert wird.

Die Autoren identifizieren eine systematische Inkonsistenz in dieser Extrapolation. Die Kalibrierungsprobe ist zu spät-Typ-Galaxien mit hoher Staubkonzentration verzerrt, während das P23-Modell unterschiedliche Extinktionseigenschaften basierend auf der stellaren Masse ($M_\star$) des Wirts zuweist. Speziell nimmt das P23-Modell einen niedrigen totalen-zu-selektiven Extinktionskoeffizienten ($R_B \approx 3,1$) für massereiche Wirte ($M_\star > 10^{10} M_\odot$) und einen Standardwert ($R_B \approx 4,0$) für massearme Wirte an. Die Autoren argumentieren, dass diese Annahme zu einer Unterschätzung der Extinktion bei gerötten Supernovae in massereichen Kalibrierungsgalaxien führt, was zu systematisch schwächeren abgeleiteten absoluten Magnituden und folglich zu einer überschätzten $H_0$ führt. Darüber hinaus steht die Annahme des P23-Modells von $R_B \approx 3,1$ für massereiche Wirte im Widerspruch zur Milchstraßen-ähnlichen Extinktionskurve ($R_B \approx 4,3$), die zur Korrektur der Cepheiden-Farben in denselben Galaxien verwendet wird.

Methodik
Die Autoren analysieren die Pantheon+-Supernova-Kompilation und die SH0ES-Cepheiden-Entfernungsmodule (Riess et al. 2022) neu. Ihre Methodik umfasst:

1. Konsistenzprüfung: Sie führen Likelihood-Analysen an der Kalibrierungsprobe durch, wobei sie die Supernovae in hochmassige ($M_\star > 10^{10} M_\odot$) und massearme Bins unterteilen. Sie vergleichen die am besten passende absolute Magnitude ($M_B$) und die daraus abgeleitete $H_0$ dieser Teilproben gegen zufällige Kontrollproben. Sie stellen fest, dass massereiche Wirte systematisch schwächere $M_B$ (höheres $H_0$) liefern als massearme Wirte und Kontrollproben, insbesondere bei roten Supernovae ($c > 0$).
2. Modifikationsmodell: Sie schlagen eine „minimalistische“ Modifikation des Extinktionsmodells vor, die ausschließlich auf die Kalibrierungsgalaxien angewendet wird, während das P23-Modell für den Hubble-Flow (wo das Modell trainiert wurde) beibehalten wird. Das neue Modell beinhaltet zwei wesentliche Änderungen:
   • Uniforme $R_B$-Verteilung: Anstatt einer massenabhängigen $R_B$ nehmen sie eine einzige, Milchstraßen-ähnliche Gauß-Verteilung für den totalen-zu-selektiven Extinktionskoeffizienten in beiden Massen-Bins an, mit einem Mittelwert $\langle R_B \rangle = 4,3$ und einer Streuung $\sigma_{R_B} = 0,4$. Dies stimmt die Supernova-Extinktionskorrektur mit der Kurve ab, die für Cepheiden verwendet wird.
   • Modifizierte Form der Rötungsverteilung: Um die effektive Steigung ($\beta \approx 3,0$) der Supernova-Peak-Magnituden-Farbe-Relation und die im Hubble-Flow gemessene mittlere Rötung $\langle E(B-V) \rangle$ beizubehalten, ersetzen sie die exponentielle Verteilung der Staub-Rötung $E(B-V)$ (verwendet in P23) durch eine Gamma-Verteilung. Der Formparameter wird auf $\gamma = 3,44$ abgestimmt, was das Maximum der Verteilung zu einer Nicht-Null-Rötung verschiebt, während der korrekte Mittelwert erhalten bleibt.
3. Statistische Evaluierung: Sie berechnen die Bias-Korrekturen ($\delta$) und die Floor-Unsicherheiten ($\sigma_{floor}$) für die Kalibrierungsprobe unter Verwendung des neuen Modells neu. Sie aktualisieren die Kovarianzmatrix, um die reduzierte intrinsische Streuung infolge der geringeren $\sigma_{R_B}$ widerzuspiegeln. Schließlich führen sie gemeinsame Fits der Kalibrierungs- und Hubble-Flow-Daten durch, um eine neue $H_0$ abzuleiten.

Wesentliche Beiträge und Ergebnisse

• Identifizierung eines systematischen Bias: Die Studie zeigt, dass die Annahme des P23-Modells von $R_B \approx 3,1$ für massereiche Kalibrierungsgalaxien eine $2,0\sigma$ bis $2,3\sigma$ Spannung in den abgeleiteten absoluten Magnituden zwischen massereichen und massearmen Wirten erzeugt. Dieser Bias ist bei roten Supernovae am ausgeprägtesten.
• Verbessertes Modell-Fit: Das vorgeschlagene neue Extinktionsmodell liefert einen signifikant besseren Fit an die Kalibrierungsdaten, mit einer Verbesserung des Bayesianischen Informationskriteriums von $\Delta \text{BIC} = -11,0$ gegenüber dem P23-Modell. Diese Verbesserung wird sowohl durch eine Reduktion von $\chi^2_{min}$ (aufgrund korrigierter Biases) als auch durch eine Änderung der Likelihood-Normalisierung (aufgrund reduzierter Kovarianz-Unsicherheiten) getrieben.
• Revidierte Hubble-Konstante: Die Anwendung des neuen Modells auf den kombinierten Kalibrierungs- und Hubble-Flow-Datensatz führt zu einer niedrigeren Best-Fit-Hubble-Konstante:
  $$H_0 = 70,5 \pm 1,0 \text{ km s}^{-1} \text{ Mpc}^{-1}$$
  Dies stellt eine Reduktion von etwa $2,9 \text{ km s}^{-1} \text{ Mpc}^{-1}$ gegenüber dem Standard-SH0ES-Wert ($73,4 \pm 1,0$) dar.
• Reduktion der Spannung: Die revidierte $H_0$ reduziert die statistische Spannung mit der Planck-$\Lambda$CDM-Messung von $5,2\sigma$ auf $2,8\sigma$. Das Ergebnis ist zudem konsistent mit unabhängigen Messungen basierend auf dem Tip of the Red Giant Branch (TRGB).

Bedeutung und Behauptungen
Das Paper behauptet, dass die „Hubble-Spannung“ teilweise oder weitgehend durch nicht berücksichtigte systematische Fehler in der Extinktionsmodellierung von Kalibrierungsgalaxien getrieben werden könnte. Die Autoren argumentieren, dass der Standardansatz fehlschlägt, weil er ein Modell, das an einer gemischten Population von Galaxien trainiert wurde, auf eine Kalibrierungsprobe extrapoliert, die ausschließlich aus spät-Typ, sternbildenden Systemen besteht.

Die Bedeutung ihrer Arbeit liegt in der Demonstration, dass:

1. Selektionsbias entscheidend ist: Die spezifische Auswahl von Kalibrierungsgalaxien (die beobachtbare Cepheiden erfordern) führt einen Bias ein, der nicht einfach durch die Auswahl von „spät-Typ“ Galaxien im Hubble-Flow erfasst wird, da die letztgenannte Auswahl die lokalen Staubumgebungen der Kalibrierungs-Supernovae nicht präzise widerspiegeln muss.
2. Extinktionsphysik: Die Annahme einer niedrigen $R_B$ ($\approx 3,1$) für massereiche Wirte ist physikalisch inkonsistent mit der Milchstraßen-ähnlichen Extinktion, die für Cepheiden in denselben Galaxien erforderlich ist, sowie mit den Beobachtungsdaten für Sternentstehungsgalaxien.
3. Sekundäreffekte: Die Form der Rötungsverteilung (speziell das zweite Moment) ist entscheidend. Die Autoren zeigen, dass eine Gamma-Verteilung ($\gamma \approx 3,4$) die Daten besser beschreibt als das Standard-Exponentialmodell, was auf eine stärkere Korrelation zwischen Staub und Supernova-Positionen in diesen spezifischen Galaxien hindeutet.

Die Autoren kommen zu dem Schluss, dass ihr revidiertes Extinktionsmodell einen konsistenteren physikalischen Rahmen für die Entfernungsleiter bietet, die scheinbare Diskrepanz zwischen massereichen und massearmen Wirten mildert und die scheinbare $H_0$-Spannung signifikant senkt, ohne neue Physik einzuführen. Sie merken an, dass, obwohl ihr Modell durch die aktuellen Daten stark bevorzugt wird, weitere Verfeinerungen unter Verwendung unabhängiger Nahinfrarot-Beobachtungen und einer vollständigen Vorwärtsmodellierung der Staub-Eigenschaften notwendig sind, um die verbleibenden Residuen vollständig aufzulösen.