Time delay measurements with Broken Power Law model

Título: Medições de Atraso Temporal com o Modelo de Lei de Potência Quebrada (Broken Power Law - BPL)

1. Problema e Contexto

A medição precisa da constante de Hubble ($H_0$) através da cosmografia de lentes gravitacionais fortes enfrenta um desafio fundamental: a degenerescência da folha de massa (Mass-Sheet Degeneracy - MSD). Esta degenerescência permite que diferentes perfis de massa do lente (a galáxia que distorce a luz) produzam as mesmas posições e fluxos das imagens, mas resultem em atrasos temporais e distâncias de atraso temporal ($D_{\Delta t}$) diferentes.

A maioria dos estudos anteriores assume um perfil de massa de lei de potência simples (como o modelo EPL - Elliptical Power Law), que é aproximadamente isoterma. No entanto, modelos físicos mais complexos, como aqueles envolvendo Dark Matter auto-interagente (SIDM) ou feedback bariônico, podem sugerir perfis com "núcleos" (cores) ou inclinações internas diferentes. A questão central é: quão sensíveis são as inferências cosmológicas ($H_0$) às suposições sobre o perfil de massa radial do lente, especialmente na região interna?

2. Metodologia

Os autores desenvolveram e implementaram uma abordagem robusta para investigar essa sensibilidade:

• Implementação do Modelo BPL: Eles integraram o modelo de Lei de Potência Quebrada (Broken Power Law - BPL) no framework de modelagem de lentes Lenstronomy. O modelo BPL permite que a inclinação da densidade de massa mude de um valor interno ($\alpha_c$) para um valor externo ($\alpha$) em um raio de quebra ($r_c$). Isso oferece maior flexibilidade para capturar desvios de um perfil de lei de potência único, simulando núcleos de densidade ou transições suaves.
• Solução Numérica para Potencial de Fermat: Como não existe uma forma analítica fechada para o potencial de deflexão do modelo BPL, os autores desenvolveram um solucionador numérico para integrar o campo de deflexão e calcular o potencial de Fermat, necessário para prever os atrasos temporais.
• Validação: O código foi validado comparando resultados de dados simulados com o modelo EPL (caso limite do BPL) e verificando a precisão da integração numérica contra soluções analíticas conhecidas.
• Aplicação ao Sistema WGD 2038–4008: O estudo foi aplicado ao quasar quadriplicado WGD 2038–4008. Foram utilizados:
  • Imagens do HST (Hubble Space Telescope) em múltiplas bandas.
  • Dados de cinemática estelar (dispersão de velocidade) da galáxia lente.
  • Estimativas de convergência externa ($\kappa_{ext}$) baseadas no ambiente da linha de visão (LoS).
• Análise Bayesiana: Foi realizada uma inferência Bayesiana comparando as famílias de modelos EPL e BPL, utilizando o Critério de Informação Bayesiano (BIC) para ponderar diferentes configurações e marginalizar sobre incertezas de modelagem.

3. Contribuições Principais

• Implementação do BPL no Lenstronomy: Disponibilização de um modelo de massa mais flexível e numericamente estável para a comunidade de lentes gravitacionais.
• Quantificação da Incerteza Sistemática: Demonstração clara de que a escolha do perfil de massa (EPL vs. BPL) introduz um viés sistemático significativo na inferência de $H_0$, mesmo quando ambos os modelos ajustam bem os dados de imagem.
• Análise de Degenerescência Interna: Ilustração de como a flexibilidade radial do modelo BPL interage com a degenerescência da folha de massa interna ($\lambda_{int}$) e a convergência externa, alterando a escala de tempo efetiva.

4. Resultados Chave

• Ajuste de Imagem: Tanto o modelo EPL quanto o BPL ajustam os dados de imagem (posições, morfologia e fluxos) com qualidade estatística comparável. Os dados de imagem sozinhos não conseguem distinguir entre os dois perfis de massa, pois são mais sensíveis à massa próxima ao raio de Einstein do que ao perfil interno profundo.
• Diferenças nos Atrasos Temporais: Ao incorporar a cinemática estelar e a convergência externa para quebrar a degenerescência da folha de massa, os modelos produzem previsões de atraso temporal sistematicamente diferentes:
  • O modelo BPL tende a prever atrasos temporais maiores (em magnitude) do que o EPL.
  • A diferença nos atrasos temporais relativos aumenta de ~9-16% (apenas imagem) para ~20-24% após a correção cinemática/LoS.
• Inferência de $H_0$: A diferença nos atrasos temporais se traduz diretamente em uma diferença na constante de Hubble inferida (assumindo $\Omega_m = 0.3$ e cosmologia $\Lambda$CDM plana):
  • Modelo EPL: $H_0 = 61.1^{+19.2}_{-13.2}$ km s$^{-1}$ Mpc$^{-1}$.
  • Modelo BPL: $H_0 = 75.3^{+23.1}_{-16.3}$ km s$^{-1}$ Mpc$^{-1}$.
  • Nota: Mesmo em um caso diagnóstico onde o fator de folha de massa interna é fixado em 1, a diferença persiste ($66.1$ vs $74.2$), embora seja menor.
• Fatores de Folha de Massa: O fator de folha de massa total inferido ($\Lambda_{tot} = (1-\kappa_{ext})\lambda_{int}$) difere entre os modelos, mostrando que a suposição do perfil de massa interno altera como a degenerescência é resolvida pelos dados de cinemática.

5. Significado e Conclusões

• Sensibilidade a Modelos: O estudo demonstra que a cosmografia de atraso temporal é altamente sensível às suposições sobre o perfil de massa do lente. A diferença de ~14 km s$^{-1}$ Mpc$^{-1}$ na inferência de $H_0$ entre os dois modelos é significativa e comparável à tensão atual de Hubble.
• Não é uma Detecção Física Direta: A diferença nos resultados não deve ser interpretada como uma detecção definitiva de um núcleo de matéria escura (core) na galáxia lente, pois os dados atuais não restringem fortemente o perfil interno (sem imagem central detectada). Em vez disso, reflete como diferentes famílias de modelos exploram o espaço de degenerescência permitido.
• Importância da Modelagem Flexível: Os resultados enfatizam que, para atingir a precisão necessária para resolver a tensão de Hubble, é crucial adotar modelos de massa mais flexíveis (como o BPL) e marginalizar explicitamente sobre as incertezas do perfil de massa, em vez de confiar em modelos de lei de potência única rígidos.
• Futuro: A precisão atual de um único sistema de lente ainda é limitada por essas incertezas sistemáticas. A combinação hierárquica de múltiplos lentes com modelagem flexível e tratamento consistente de cinemática e efeitos da linha de visão será essencial para futuras medições robustas de $H_0$.

Em resumo, o artigo alerta que a "escolha do modelo" não é apenas um detalhe técnico, mas uma fonte dominante de incerteza sistemática na cosmografia de lentes gravitacionais, exigindo uma abordagem mais sofisticada para a modelagem de massa.