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

Sintesi Tecnica: Modello di Estinzione Consistente per le Supernovae di Tipo Ia nelle Galassie di Calibrazione basate su Cefeidi e il suo Impatto su $H_0$

Enunciato del Problema
L'attuale tensione in cosmologia deriva da una significativa discrepanza ($5,2\sigma$) tra la costante di Hubble ($H_0$) misurata tramite la scala delle distanze locale (utilizzando Cefeidi e supernovae di tipo Ia nel programma SH0ES) e il valore inferito dalla Radiazione Cosmica di Fondo (Planck) assumendo un modello standard $\Lambda$CDM. Una componente critica della misurazione SH0ES riguarda la correzione delle magnitudini delle supernovae di tipo Ia per l'estinzione della galassia ospite. L'approccio standard (Popovic et al. 2023, qui in seguito P23) impiega un modello di polvere probabilistico addestrato su supernovae in flusso di Hubble ($z > 0,03$) ed estrapola tale modello alle galassie di calibrazione (quelle con Cefeidi osservate).

Gli autori identificano un'inconsistenza sistematica in questa estrapolazione. Il campione di calibrazione è distorto verso galassie di tipo tardivo e ricche di polvere, mentre il modello P23 assegna diverse proprietà di estinzione in base alla massa stellare dell'ospite ($M_\star$). Nello specifico, il modello P23 assume un coefficiente di estinzione totale-selettiva più basso ($R_B \approx 3,1$) per gli ospiti ad alta massa ($M_\star > 10^{10} M_\odot$) e un valore standard ($R_B \approx 4,0$) per gli ospiti a bassa massa. Gli autori sostengono che tale assunzione porti a una sottostima dell'estinzione per le supernovae arrossate in galassie di calibrazione ad alta massa, risultando in magnitudini assolute inferite sistematicamente più deboli e, di conseguenza, in una $H_0$ sovrastimata. Inoltre, l'assunzione $R_B \approx 3,1$ del modello P23 per gli ospiti ad alta massa entra in conflitto con la curva di estinzione simile alla Via Lattea ($R_B \approx 4,3$) utilizzata per correggere i colori delle Cefeidi nelle stesse galassie.

Metodologia
Gli autori rianalizzano la compilazione di supernovae Pantheon+ e i moduli di distanza delle Cefeidi SH0ES (Riess et al. 2022). La loro metodologia prevede:

1. Test di Consistenza: Effettuano analisi di verosimiglianza sul campione di calibrazione, separando le supernovae in bin di alta massa ($M_\star > 10^{10} M_\odot$) e bassa massa. Confrontano la magnitudine assoluta ($M_B$) migliore e la derivata $H_0$ per questi sottocampioni rispetto a campioni di controllo casuali. Riscontrano che gli ospiti ad alta massa producono una $M_B$ sistematicamente più debole (maggiore $H_0$) rispetto agli ospiti a bassa massa e ai campioni di controllo, particolarmente per le supernovae rosse ($c > 0$).
2. Modifica del Modello: Propongono una modifica "minimalista" al modello di estinzione applicato esclusivamente alle galassie di calibrazione, mantenendo il modello P23 per il flusso di Hubble (dove il modello è stato addestrato). Il nuovo modello prevede due cambiamenti chiave:
   • Distribuzione Uniforme di $R_B$: Invece di una $R_B$ dipendente dalla massa, assumono una singola distribuzione gaussiana simile a quella della Via Lattea per il coefficiente di estinzione totale-selettiva sia nei bin di massa, con una media $\langle R_B \rangle = 4,3$ e dispersione $\sigma_{R_B} = 0,4$. Ciò allinea la correzione dell'estinzione delle supernovae con la curva utilizzata per le Cefeidi.
   • Modifica della Forma della Distribuzione di Arrossamento: Per preservare la pendenza effettiva ($\beta \approx 3,0$) della relazione colore-magnitudine di picco della supernova e il reddening medio $\langle E(B-V) \rangle$ misurato nel flusso di Hubble, sostituiscono la distribuzione esponenziale del reddening della polvere $E(B-V)$ (usata in P23) con una distribuzione Gamma. Il parametro di forma è tarato su $\gamma = 3,44$, il che sposta il picco della distribuzione verso un reddening non nullo pur mantenendo la corretta media.
3. Valutazione Statistica: Ricalcolano le correzioni di bias ($\delta$) e le incertezze di soglia ($\sigma_{floor}$) per il campione di calibrazione utilizzando il nuovo modello. Aggiornano la matrice di covarianza per riflettere la riduzione della dispersione intrinseca derivante dalla minore $\sigma_{R_B}$. Infine, eseguono fit congiunti dei dati di calibrazione e di flusso di Hubble per derivare una nuova $H_0$.

Contributi Chiave e Risultati

• Identificazione del Bias Sistematico: Lo studio dimostra che l'assunzione del modello P23 di $R_B \approx 3,1$ per le galassie di calibrazione ad alta massa crea una tensione di $2,0\sigma$ a $2,3\sigma$ nelle magnitudini assolute derivate tra ospiti ad alta e bassa massa. Questo bias è più pronunciato per le supernovae rosse.
• Miglioramento del Fit del Modello: Il nuovo modello di estinzione proposto produce un fit significativamente migliore ai dati di calibrazione, con un miglioramento del Criterio di Informazione Bayesiana di $\Delta \text{BIC} = -11,0$ rispetto al modello P23. Questo miglioramento è guidato sia da una riduzione di $\chi^2_{min}$ (dovuta ai bias corretti) sia da un cambiamento nella normalizzazione della verosimiglianza (dovuto alle ridotte incertezze di covarianza).
• Hubble Costante Rivista: L'applicazione del nuovo modello al dataset congiunto di calibrazione e flusso di Hubble risulta in una costante di Hubble con un valore migliore di:
  $$H_0 = 70,5 \pm 1,0 \text{ km s}^{-1} \text{ Mpc}^{-1}$$
  Ciò rappresenta una riduzione di circa $2,9 \text{ km s}^{-1} \text{ Mpc}^{-1}$ rispetto al valore standard SH0ES ($73,4 \pm 1,0$).
• Riduzione della Tensione: La $H_0$ rivista riduce la tensione statistica con la misurazione Planck $\Lambda$CDM da $5,2\sigma$ a $2,8\sigma$. Il risultato è inoltre coerente con le misurazioni indipendenti basate sul Tip of the Red Giant Branch (TRGB).

Significatività e Rivendicazioni
L'articolo sostiene che la "tensione di Hubble" possa essere parzialmente o ampiamente guidata da errori sistematici non contabilizzati nella modellazione dell'estinzione delle galassie di calibrazione. Gli autori argomentano che l'approccio standard fallisce perché estrapola un modello addestrato su una popolazione mista di galassie a un campione di calibrazione composto esclusivamente da sistemi a formazione stellare e di tipo tardivo.

La significatività del loro lavoro risiede nel dimostrare che:

1. Il Bias di Selezione Conta: La specifica selezione delle galassie di calibrazione (che richiedono Cefeidi osservabili) introduce un bias che non viene catturato semplicemente selezionando galassie di "tipo tardivo" nel flusso di Hubble, poiché tale selezione potrebbe non corrispondere precisamente agli ambienti di polvere locali delle supernovae di calibrazione.
2. Fisica dell'Estinzione: L'assunzione di una bassa $R_B$ ($\approx 3,1$) per gli ospiti ad alta massa è fisicamente incoerente con l'estinzione simile alla Via Lattea richiesta per le Cefeidi nelle stesse galassie e con i vincoli osservativi sulle galassie a formazione stellare.
3. Effetti di Secondo Ordine: La forma della distribuzione di arrossamento (specificamente il secondo momento) è critica. Gli autori mostrano che una distribuzione Gamma ($\gamma \approx 3,4$) descrive meglio i dati rispetto al modello esponenziale standard, implicando una correlazione più forte tra polvere e posizioni delle supernovae in queste specifiche galassie.

Gli autori concludono che il loro modello di estinzione rivisto fornisce un quadro fisico più coerente per la scala delle distanze, alleviando la discrepanza apparente tra ospiti ad alta e bassa massa e abbassando significativamente la tensione su $H_0$ senza invocare nuova fisica. Notano che, sebbene il loro modello sia fortemente favorito dai dati attuali, è necessaria un'ulteriore raffinazione utilizzando osservazioni indipendenti nel vicino infrarosso e una modellazione diretta completa delle proprietà della polvere per risolvere completamente i residui rimanenti.