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

Résumé technique : Modèle d'extinction cohérent pour les supernovae de type Ia dans les galaxies de calibration basées sur les Céphéides et son impact sur $H_0$

Énoncé du problème
La tension actuelle en cosmologie provient d'une divergence significative ($5,2\sigma$) entre la constante de Hubble ($H_0$) mesurée via l'échelle de distance locale (utilisant les Céphéides et les supernovae de type Ia dans le programme SH0ES) et la valeur déduite du fond diffus cosmologique (Planck) en supposant un modèle $\Lambda$CDM standard. Une composante critique de la mesure SH0ES implique la correction des magnitudes des supernovae de type Ia pour l'extinction des galaxies hôtes. L'approche standard (Popovic et al. 2023, ci-après P23) emploie un modèle de poussière probabiliste entraîné sur des supernovae en flux de Hubble ($z > 0,03$) et l'extrapole aux galaxies de calibration (celles présentant des Céphéides observées).

Les auteurs identifient une incohérence systématique dans cette extrapolation. L'échantillon de calibration est biaisé en faveur de galaxies de type tardif, riches en poussière, alors que le modèle P23 assigne des propriétés d'extinction différentes en fonction de la masse stellaire de l'hôte ($M_\star$). Plus précisément, le modèle P23 suppose un coefficient d'extinction totale par rapport à l'extinction sélective faible ($R_B \approx 3,1$) pour les hôtes de haute masse ($M_\star > 10^{10} M_\odot$) et une valeur standard ($R_B \approx 4,0$) pour les hôtes de faible masse. Les auteurs soutiennent que cette hypothèse conduit à une sous-estimation de l'extinction pour les supernovae rougies dans les galaxies de calibration de haute masse, résultant en des magnitudes absolues inférées systématiquement plus faibles et, par conséquent, en une surestimation de $H_0$. De plus, le $R_B \approx 3,1$ du modèle P23 pour les hôtes de haute masse entre en conflit avec la courbe d'extinction de type Voie Lactée ($R_B \approx 4,3$) utilisée pour corriger les couleurs des Céphéides dans ces mêmes galaxies.

Méthodologie
Les auteurs réanalysent la compilation de supernovae Pantheon+ et les modules de distance des Céphéides SH0ES (Riess et al. 2022). Leur méthodologie comprend :

1. Test de cohérence : Ils effectuent des analyses de vraisemblance sur l'échantillon de calibration, en séparant les supernovae en classes de haute masse ($M_\star > 10^{10} M_\odot$) et de faible masse. Ils comparent la magnitude absolue ($M_B$) la mieux ajustée et la $H_0$ dérivée pour ces sous-échantillons par rapport à des échantillons de contrôle aléatoires. Ils constatent que les hôtes de haute masse produisent des $M_B$ systématiquement plus faibles ( $H_0$ plus élevée) que les hôtes de faible masse et les échantillons de contrôle, particulièrement pour les supernovae rouges ($c > 0$).
2. Modification du modèle : Ils proposent une modification « minimaliste » du modèle d'extinction appliquée uniquement aux galaxies de calibration, tout en conservant le modèle P23 pour le flux de Hubble (où le modèle a été entraîné). Le nouveau modèle implique deux changements clés :
   • Distribution uniforme de $R_B$ : Au lieu d'un $R_B$ dépendant de la masse, ils supposent une distribution gaussienne unique, de type Voie Lactée, pour le coefficient d'extinction totale par rapport à l'extinction sélective dans les deux classes de masse, avec une moyenne $\langle R_B \rangle = 4,3$ et un écart-type $\sigma_{R_B} = 0,4$. Cela aligne la correction d'extinction des supernovae avec la courbe utilisée pour les Céphéides.
   • Modification de la forme de la distribution de l'affadissement : Pour préserver la pente effective ($\beta \approx 3,0$) de la relation magnitude de pic de la supernova–couleur et l'affadissement moyen $\langle E(B-V) \rangle$ mesuré dans le flux de Hubble, ils remplacent la distribution exponentielle de l'affadissement de la poussière $E(B-V)$ (utilisée dans P23) par une distribution de Gamma. Le paramètre de forme est ajusté à $\gamma = 3,44$, ce qui déplace le pic de la distribution vers un affadissement non nul tout en maintenant la moyenne correcte.
3. Évaluation statistique : Ils recomputent les corrections de biais ($\delta$) et les incertitudes de plancher ($\sigma_{floor}$) pour l'échantillon de calibration en utilisant le nouveau modèle. Ils mettent à jour la matrice de covariance pour refléter la réduction de la dispersion intrinsèque résultant du $\sigma_{R_B}$ plus faible. Enfin, ils effectuent des ajustements conjoints des données de calibration et de flux de Hubble pour dériver une nouvelle $H_0$.

Contributions clés et résultats

• Identification d'un biais systématique : L'étude démontre que l'hypothèse du $R_B \approx 3,1$ du modèle P23 pour les galaxies de calibration de haute masse crée une tension de $2,0\sigma$ à $2,3\sigma$ dans les magnitudes absolues dérivées entre les hôtes de haute et de faible masse. Ce biais est plus prononcé pour les supernovae rouges.
• Amélioration de l'ajustement du modèle : Le nouveau modèle d'extinction proposé produit un ajustement nettement meilleur aux données de calibration, avec une amélioration du critère d'information bayésien de $\Delta \text{BIC} = -11,0$ par rapport au modèle P23. Cette amélioration est portée à la fois par une réduction de $\chi^2_{min}$ (due aux biais corrigés) et par un changement dans la normalisation de la vraisemblance (due à la réduction des incertitudes de covariance).
• Révision de la constante de Hubble : L'application du nouveau modèle à l'ensemble de données conjoint (calibration et flux de Hubble) aboutit à une constante de Hubble plus faible :
  $$H_0 = 70,5 \pm 1,0 \text{ km s}^{-1} \text{ Mpc}^{-1}$$
  Cela représente une réduction d'environ $2,9 \text{ km s}^{-1} \text{ Mpc}^{-1}$ par rapport à la valeur SH0ES standard ($73,4 \pm 1,0$).
• Réduction de la tension : La $H_0$ révisée réduit la tension statistique avec la mesure Planck $\Lambda$CDM de $5,2\sigma$ à $2,8\sigma$. Le résultat est également cohérent avec les mesures indépendantes basées sur la pointe de la branche des géantes rouges (TRGB).

Signification et affirmations
L'article affirme que la « tension de Hubble » peut être partiellement ou largement pilotée par des erreurs systématiques non prises en compte dans la modélisation de l'extinction des galaxies de calibration. Les auteurs soutiennent que l'approche standard échoue car elle extrapole un modèle entraîné sur une population mixte de galaxies à un échantillon de calibration composé exclusivement de systèmes de type tardif et formateurs d'étoiles.

La signification de leur travail réside dans la démonstration que :

1. Le biais de sélection importe : La sélection spécifique des galaxies de calibration (nécessitant des Céphéides observables) introduit un biais qui n'est pas capturé par une simple sélection de galaxies de « type tardif » dans le flux de Hubble, car cette dernière sélection peut ne pas correspondre précisément aux environnements de poussière locaux des supernovae de calibration.
2. Physique de l'extinction : L'hypothèse d'un $R_B$ faible ($\approx 3,1$) pour les hôtes de haute masse est physiquement incohérente avec l'extinction de type Voie Lactée requise pour les Céphéides dans ces mêmes galaxies et avec les contraintes observationnelles sur les galaxies formant des étoiles.
3. Effets de second ordre : La forme de la distribution de l'affadissement (spécifiquement le second moment) est critique. Les auteurs montrent qu'une distribution de Gamma ($\gamma \approx 3,4$) décrit mieux les données qu'un modèle exponentiel standard, impliquant une corrélation plus forte entre la poussière et la position des supernovae dans ces galaxies spécifiques.

Les auteurs concluent que leur modèle d'extinction révisé fournit un cadre physique plus cohérent pour l'échelle de distance, atténuant la divergence apparente entre les hôtes de haute et de faible masse, et abaissant significativement la tension $H_0$ sans invoquer de nouvelle physique. Ils notent que bien que leur modèle soit fortement favorisé par les données actuelles, un raffinement supplémentaire utilisant des observations indépendantes dans l'infrarouge proche et une modélisation directe complète des propriétés de la poussière est nécessaire pour résoudre pleinement les résidus restants.