Time delay measurements with Broken Power Law model

Titolo: Misurazioni dei ritardi temporali con il modello Broken Power Law (BPL)

1. Il Problema

La determinazione precisa della costante di Hubble ($H_0$) rimane una delle sfide più critiche della cosmologia moderna, caratterizzata dalla "tensione di Hubble" tra le misurazioni dell'universo primordiale (CMB, Planck) e quelle dell'universo locale (supernove di Tipo Ia, SH0ES). La cosmografia tramite lenti gravitazionali forti offre un metodo geometrico indipendente per misurare $H_0$ sfruttando i ritardi temporali tra le immagini multiple di una sorgente variabile (es. quasar).

Tuttavia, un ostacolo fondamentale è la degenerazione del foglio di massa (Mass-Sheet Degeneracy, MSD). Le osservazioni di imaging (posizioni e rapporti di flusso delle immagini) sono invarianti sotto trasformazioni di scala della convergenza di massa, rendendo difficile distinguere tra diversi profili di densità radiale. In particolare, i modelli standard assumono spesso un profilo di potenza ellittico (EPL, Elliptical Power Law), che potrebbe non catturare la complessità fisica reale, come la presenza di nuclei di materia oscura "piatti" (core) o feedback barionico. La rigidità di questi modelli può introdurre errori sistematici significativi nelle inferenze cosmologiche.

2. Metodologia

Gli autori hanno sviluppato e implementato un nuovo approccio all'interno del framework Lenstronomy:

• Modello Broken Power Law (BPL): È stato introdotto un profilo di massa che permette una transizione nella pendenza radiale della densità tra una regione interna (con pendenza $\alpha_c$) e una esterna (con pendenza $\alpha$). Questo modello è più flessibile dell'EPL e può phenomenologicamente descrivere strutture con "core" (tipiche di simulazioni di Materia Oscura Interagente, SIDM) o gradienti ripidi (CDM).
• Soluzione Numerica: Poiché il potenziale di Fermat per il modello BPL non ha una forma analitica chiusa, gli autori hanno sviluppato un solutore numerico per integrare il campo di deflessione e calcolare il potenziale, validando l'accuratezza confrontandolo con la soluzione analitica nota per l'EPL.
• Caso di Studio: Il modello è stato applicato al quasar quadruplamente lensato WGD 2038–4008.
• Dati e Vincoli: L'analisi ha combinato:
  • Immagini multi-banda HST (WFC3).
  • Cinematica stellare (dispersione di velocità $\sigma_{los}$) della lente.
  • Stima della convergenza esterna ($\kappa_{ext}$) basata sull'ambiente lungo la linea di vista (LoS).
• Inferenza Bayesiana: Sono stati eseguiti confronti tra le famiglie di modelli EPL e BPL utilizzando il Criterio di Informazione Bayesiano (BIC) pesato, combinando i dati di imaging, cinematica e LoS per vincolare i fattori di scala del foglio di massa ($\lambda_{int}$ e $\kappa_{ext}$).

3. Contributi Chiave

• Implementazione del BPL in Lenstronomy: Integrazione di un profilo di massa a doppia pendenza con un solutore numerico robusto per i ritardi temporali.
• Validazione su Dati Simulati: Dimostrazione che l'uso di un modello EPL su dati simulati con un vero profilo BPL introduce bias significativi, mentre il modello BPL recupera correttamente i parametri di input.
• Analisi Comparativa su WGD 2038–4008: Prima applicazione sistematica del modello BPL a un sistema reale con dati cinematici completi, confrontandolo direttamente con l'EPL standard.
• Quantificazione dell'Incertezza Sistematica: Evidenziazione di come le assunzioni sul profilo di massa interno influenzino l'inferenza di $H_0$ anche quando i dati di imaging sono ugualmente ben adattati da entrambi i modelli.

4. Risultati

• Adattamento ai Dati: Sia il modello EPL che il BPL adattano i dati di imaging e cinematica con qualità statisticamente comparabile (valori di $\chi^2$ ridotti simili). Non ci sono prove osservative decisive (come un'immagine centrale rilevata) che favoriscano un modello sull'altro basandosi solo sull'imaging.
• Differenze Strutturali:
  • L'EPL richiede una pendenza radiale singola ($\gamma \approx 2.30$).
  • Il BPL preferisce un profilo spezzato con una regione interna più piatta ($\alpha_c \approx 1.26$) e una regione esterna più ripida ($\alpha \approx 2.80$), con un raggio di rottura di circa $0.4''$.
• Impatto sui Ritardi Temporali: Sebbene le posizioni delle immagini siano simili, i profili di massa diversi portano a differenze sistematiche nei ritardi temporali predetti (fino al 20-24% dopo correzioni cinematiche).
• Inferenza di $H_0$:
  • Con il fattore di foglio di massa interno ($\lambda_{int}$) libero di variare:
    • Modello BPL: $H_0 = 75.3^{+23.1}_{-16.3}$ km s$^{-1}$ Mpc$^{-1}$.
    • Modello EPL: $H_0 = 61.1^{+19.2}_{-13.2}$ km s$^{-1}$ Mpc$^{-1}$.
  • In un caso diagnostico con $\lambda_{int}$ fissato a 1, la differenza persiste ma è leggermente ridotta ($74.2$ vs $66.1$).
• Interpretazione: La differenza nell'inferenza di $H_0$ non deriva da una misurazione diretta del nucleo, ma da come la flessibilità radiale del modello BPL ridistribuisce la massa, interagendo con la degenerazione del foglio di massa e i vincoli cinematici.

5. Significato e Conclusioni

Questo studio sottolinea che le assunzioni sul profilo di massa della lente costituiscono una fonte di incertezza sistematica non trascurabile nella cosmografia a ritardi temporali.

• La differenza di circa $14$ km s$^{-1}$ Mpc$^{-1}$ tra i risultati EPL e BPL per un singolo sistema dimostra che la scelta del modello fisico può spostare significativamente le stime cosmologiche, anche con dati di alta qualità.
• Attualmente, con le incertezze statistiche, questo scarto non favorisce uno scenario cosmologico specifico rispetto all'altro, ma evidenzia la necessità di modelli di massa flessibili e di una marginalizzazione esplicita sulle incertezze del profilo radiale.
• Per risolvere la tensione di Hubble con la precisione richiesta, le future analisi dovranno integrare modelli più complessi (come il BPL) con dati cinematici risolti spazialmente e osservazioni di immagini centrali per vincolare direttamente la densità interna delle lenti.

In sintesi, il lavoro di Rui et al. fornisce uno strumento metodologico cruciale (il modello BPL in Lenstronomy) e una prova empirica che la modellazione della massa è un fattore limitante critico per la precisione della cosmologia con lenti gravitazionali.