TDCOSMO XXVI: Uniform lens modeling of eight doubly imaged quasars

1. Declaração do Problema

A cosmografia de atraso temporal (TDC) é um método primário para medir a constante de Hubble ($H_0$) a fim de resolver a "tensão de Hubble" entre medições do universo primordial (CMB) e do universo tardio (escada de distâncias). Embora quasares com quatro imagens ("quads") tenham sido o padrão para TDC devido às suas ricas restrições geométricas, quasares com duas imagens ("doubles") são aproximadamente quatro vezes mais abundantes no universo. No entanto, os "doubles" historicamente foram subutilizados porque oferecem menos observáveis (apenas duas posições de imagem), tornando-os mais suscetíveis à Degenerescência do Folha de Massa (MSD) e à transformação de folha de massa (MST).

O desafio central abordado neste artigo é a falta de uma estrutura de modelagem uniforme e escalável especificamente projetada para sistemas com duas imagens. Análises anteriores frequentemente dependiam de conjuntos de dados heterogêneos e técnicas de modelagem, impedindo a construção de grandes amostras estatisticamente robustas necessárias para inferência cosmológica bayesiana hierárquica. Além disso, a contribuição específica dos arcos estendidos da galáxia hospedeira (brilho superficial) para quebrar degenerescências em "doubles" não havia sido quantificada sistematicamente.

2. Metodologia

Os autores utilizaram a estrutura de software de código aberto Lenstronomy para desenvolver um pipeline de modelagem uniforme para oito quasares duplamente imagemados observados pelo Telescópio Espacial Hubble (HST) como parte do programa GO-17199.

• Dados: Observações multibanda usando filtros WFC3: F160W (IR), F814W e F475X (UVIS).
• Modelagem de Massa:
  • As galáxias defletoras principais foram modeladas como Elipsoides Isotérmicos Singulares (SIE).
  • Cisalhamento externo ($\gamma_{ext}$) foi incluído para contabilizar perturbações ao longo da linha de visada, com um prior conservador ($\gamma_{ext} \leq 0.15$).
  • Para sistemas com perturbadores próximos (por exemplo, J1001+5027), componentes SIE adicionais foram adicionados.
• Modelagem de Luz:
  • Lente de Luz: Modelada usando perfis Sérsic elípticos. Sistemas complexos (por exemplo, galáxias de perfil lateral como J0602-4335) utilizaram componentes duplos Sérsic.
  • Luz da Fonte: A galáxia hospedeira do quasar foi modelada como um perfil Sérsic elíptico. Para J2325-5229, Shapelets (polinômios Gauss-Hermite) foram adicionados à banda IR para capturar estruturas complexas de arcos.
  • Tratamento da PSF: Funções de Espalhamento de Ponto (PSF) foram reconstruídas usando STARRED a partir de estrelas isoladas. Um processo iterativo foi empregado onde a luz da lente foi subtraída para refinar a PSF com base nos resíduos do quasar. Um mapa de variância da PSF foi gerado para reduzir o peso de pixels onde a incerteza da PSF era alta.
• Sequência de Ajuste:
  1. Otimização Inicial de Enxame de Partículas (PSO) em F160W.
  2. Reconstrução e refinamento iterativos da PSF.
  3. Otimização conjunta em todas as três bandas (F160W, F814W, F475X).
  4. Amostragem de Cadeia de Markov Monte Carlo (MCMC) usando emcee para marginalizar totalmente as incertezas.
• Análise de Pontos Conjugados: Para isolar o valor dos dados de brilho superficial, os autores realizaram uma análise de "ponto conjugado". Este método usou apenas as posições das imagens do quasar (tratando-as como verdades observacionais) para restringir o modelo de massa, excluindo todas as informações de arcos estendidos. Isso foi comparado contra os resultados completos de modelagem de imagem.

3. Principais Contribuições

• Pipeline Uniforme: Estabeleceu o primeiro pipeline padronizado e de código aberto adaptado especificamente para modelagem de quasares duplamente imagemados, permitindo o processamento de grandes amostras de futuras pesquisas (LSST, Roman, Euclid).
• Quantificação de Informações de Arcos: Demonstrou sistematicamente que o brilho superficial dos arcos da galáxia hospedeira lenteada é o principal motor da precisão do modelo de massa em "doubles", e não apenas as posições de fontes pontuais.
• Benchmarking de Pontos Conjugados: Fornecia uma comparação rigorosa entre restrições apenas de fonte pontual e modelagem completa de brilho superficial, quantificando o "ganho de informação" fornecido por imagens de alta resolução.
• Novas Medições: Fornecia as primeiras medições uniformes de raio de Einstein e separação de imagem para toda a amostra, incluindo as primeiras medições para J0602-4335.

4. Principais Resultados

• Precisão do Modelo:
  • Raios de Einstein modelados ($\theta_E$) concordam com a literatura anterior em uma média de 1.5$\sigma$. A maior discrepância (4.5$\sigma$) ocorreu em J0407-5006, atribuída à pequena separação da Imagem B do centróide da lente, o que desafia a separação de imagens baseada em solo.
  • Separações de imagem modeladas diferem das medições do Gaia DR2 por um RMS de apenas 3.6 mas (excluindo o outlier J0407-5006), validando a precisão astrométrica do pipeline.
• Brilho do Arco vs. Precisão:
  • Uma forte correlação positiva (Pearson $r = 0.97$) foi encontrada entre o brilho superficial do arco lenteado e a precisão do potencial de Fermat ($\Delta\Phi$).
  • Arcos mais brilhantes (por exemplo, J1001+5027, $m_{arc} \approx 18.96$) produzem precisões de potencial de Fermat de $\sim0.87\%$, enquanto arcos mais fracos (por exemplo, J1515+1511, $m_{arc} \approx 21.6$) resultam em incerteza de $\sim2.89\%$.
• Ponto Conjugado vs. Modelagem Completa:
  • O hipervolume posterior (espaço de incerteza) dos parâmetros de massa na análise de ponto conjugado foi, em média, 109 vezes maior do que na modelagem completa de imagem.
  • Existe uma forte anticorrelação ($r = -0.87$) entre a magnitude do arco e o hipervolume do parâmetro de massa, confirmando que arcos mais brilhantes quebram significativamente a degenerescência massa-cisalhamento.
• Alinhamento Massa-Luz:
  • Na banda de alto S/N F160W, os ângulos de posição de massa e luz estão geralmente bem alinhados (dentro do prior de 10$^\circ$), embora desvios em alguns sistemas sugiram a influência de perturbadores próximos não incluídos no modelo simples SIE+cisalhamento.

5. Significado e Perspectivas Futuras

Este trabalho estabelece que quasares duplamente imagemados são candidatos viáveis para cosmografia de precisão, desde que imagens de alta resolução possam resolver os arcos estendidos da galáxia hospedeira. O estudo prova que:

1. Escalabilidade: Um pipeline uniforme pode lidar com a heterogeneidade dos "doubles", preparando o campo para as milhares de lentes esperadas de pesquisas futuras.
2. Estratégia Observacional: A forte dependência da precisão do potencial de Fermat no brilho do arco dita que amostras cosmológicas futuras devem priorizar seguimento profundo e de alta resolução para detectar arcos estendidos. Sem isso, a degenerescência da folha de massa não pode ser suficientemente quebrada usando apenas dados de fonte pontual.
3. Caminho Cosmológico: As saídas uniformes (raios de Einstein, potenciais de Fermat e incertezas) geradas aqui servem como entradas fundamentais para uma subsequente análise bayesiana hierárquica (a ser apresentada em um artigo futuro) que combinará atrasos temporais e cinemática estelar para inferir $H_0$ com precisão aprimorada e sistemáticos reduzidos.

Em resumo, este artigo preenche a lacuna entre a abundância de lentes duplas e sua utilidade na cosmologia, demonstrando que, com modelagem moderna e imagens profundas, os "doubles" podem complementar os "quads" na resolução da tensão de Hubble.