A Consistent Implementation of Cluster Strong Lensing in Cosmological Simulation Light Cones

Resumo Técnico: Uma Implementação Consistente de Lente Gravitacional Forte de Aglomerados em Cone de Luz de Simulações Cosmológicas

Definição do Problema
Embora o efeito de lente gravitacional forte em aglomerados de galáxias seja um pilar da cosmologia de precisão, as previsões teóricas robustas historicamente ficaram atrás das observações de alta qualidade. Essa lacuna decorre de dois problemas primários: o desajuste geométrico entre os padrões de caixas de simulação cúbicas e a geometria alongada necessária para cones de luz, e a tensão fundamental entre o grande volume necessário para abrigar aglomerados massivos e a alta resolução exigida para resolver galáxias fontes de fundo e subestruturas.

Os métodos de previsão atuais frequentemente dependem de abordagens híbridas, extraindo lentes individuais de simulações e combinando-as com populações de fontes analíticas ou observadas. Esses métodos frequentemente falham em preservar correlações através do redshift ou omitem estruturas correlacionadas e não correlacionadas na linha de visada (LoS - line-of-sight). Consequentemente, há uma carência de imagens de lente forte simuladas geradas consistentemente a partir de saídas de simulação, onde a lente primária, as fontes e toda a matéria interveniente são extraídas do mesmo volume cosmológico de maneira fiel ao paradigma padrão da matéria escura fria.

Metodologia
Os autores apresentam um procedimento totalmente baseado em simulação para gerar imagens de lente forte diretamente de dados de partículas, utilizando as simulações hidrodinâmicas cosmológicas IllustrisTNG. A metodologia consiste em três componentes principais:

1. Remapeamento de Volume: Para abordar as restrições geométricas, os autores empregam uma técnica de remapeamento de caixa (usando o código BoxRemap) para transformar o volume de simulação cúbico em um prisma retangular de base quadrada. Esta geometria apresenta dois eixos curtos (representando o plano do céu) e um eixo longo (representando a linha de visada). Esta transformação preserva o volume total e a estrutura local sem introduzir estruturas repetidas dentro de um único cone de luz, desde que o campo de visão seja escolhido apropriadamente em relação ao tamanho da caixa. O remapeamento é aplicado consistentemente a cada snapshot da simulação.
2. Extração de Planos e Suavização: Dos boxes remapeados, os autores extraem informações de massa e luz para criar uma série de planos de lente e de fonte correspondentes aos snapshots da simulação.
   • Planos de Massa: As partículas são divididas em cascas radiais baseadas nos redshifts dos snapshots. Mapas de massa são gerados suavizando partículas de matéria escura, estrelas e gás usando um kernel de vizinho mais próximo (32 vizinhos, kernel Wendland-C2). Para contabilizar a massa difusa em regiões esparsas e evitar o ruído de discretização (shot noise), a massa total da casca é adicionada como uma folha uniforme, e a densidade média de matéria do universo é subtraída para criar um mapa de contraste de massa.
   • Planos de Fonte: Partículas estelares são suavizadas em mapas de fluxo para filtros observacionais específicos (ex: JWST f200w). O fluxo é estimado com base na idade estelar, massa inicial e metalicidade, deslocado para o vermelho (redshifted) para o plano da fonte, e convertido em magnitude aparente. Subhalos com contagem insuficiente de partículas são substituídos por perfis de Sersic analíticos para garantir a completude.
3. Traçado de Raios Multi-Plano: Os planos extraídos são utilizados em um algoritmo de traçado de raios multi-plano (usando o código GLAMER). A luz de planos de fonte em $z_s$ é rastreada retroativamente através de todos os planos de lente em redshifts $z_l < z_s$. As imagens resultantes são empilhadas para produzir uma previsão teórica final. O pipeline não requer injeção manual de fontes ou perfis analíticos; a cosmologia e o modelo de formação de galáxias da simulação servem como o prior.

Contribuições Principais e Resultados
Utilizando a simulação TNG300-1, os autores geraram cones de luz centrados em grupos de galáxias massivos ($M_{200m} \ge 5 \times 10^{14} M_\odot$) em um redshift de lente de $z_G \approx 0.38$, estendendo-se até um redshift de fonte de $z_{max} \approx 4.18$.

• Geração de Imagens: O pipeline produz com sucesso características de lente forte, incluindo configurações de múltiplas imagens, arcos e curvas críticas, diretamente dos dados da simulação sem injeção de fontes ad-hoc. As imagens resultantes mostram características típicas, como arcos tangenciais e morfologias complexas de curvas críticas impulsionadas por fusões de aglomerados e subestruturas.
• Impacto da Estrutura na Linha de Visada (LoS): O estudo quantifica o impacto da estrutura de LoS não correlacionada comparando o traçado de raios de cone de luz completo contra o traçado de raios através apenas do plano da lente primária.
  • Posições de Imagem: A estrutura de LoS não correlacionada pode deslocar as posições relativas das imagens lensed em vários segundos de arco (até $\sim 10''$ em casos extremos), excedendo significativamente os resíduos típicos de modelagem de massa de lentes fortes ($\lesssim 1''$).
  • Morfologia da Curva Crítica: A inclusão de material de LoS introduz um espalhamento de $\sim 6\%$ na área da curva crítica primária de um aglomerado. Em alguns casos, a área da curva crítica primária foi aumentada em 30% ou diminuída em 50%.
  • Área Crítica Total: Dentro de $100''$ do mínimo do potencial do aglomerado, a área crítica total em $z_s = 4$ muda em $16^{+20}_{-14}\%$ quando a estrutura de LoS é incluída. Este aumento surge tanto da criação de novas curvas críticas por objetos de LoS quanto da modificação das curvas existentes.

Significância e Limitações
O artigo afirma fornecer a primeira demonstração de imagens de lente forte de aglomerados geradas de forma autossuficiente a partir de um cone de luz de uma simulação hidrodinâmica cosmológica. Esta abordagem oferece um método mais automatizado e "manual livre" para produzir grandes catálogos de observações simuladas, que podem ser usados para estudar sistemáticas em modelagem de lentes, vieses em inferência de massa e a tensão da seção de choque de lente forte galáxia-galáxia (GGSL).

Os autores observam várias limitações:

• Potência de Grande Escala: O método não introduz modos de potência de grande escala ausentes na caixa de simulação inicial, o que pode enviesar a amostragem do espectro de potência de grande escala.
• Resolução: O método é limitado pela resolução de massa da simulação, que tem dificuldade em resolver simultaneamente aglomerados massivos e subestruturas/fontes de baixa massa.
• Massa Difusa: A massa difusa é aproximada como uma folha uniforme em vez de preservar sua estrutura específica, embora isso seja necessário para evitar o ruído de discretização.
• Custo Computacional: Gerar um único cone de luz requer o processamento de cada partícula na caixa de simulação em cada snapshot, o que é computacionalmente caro e gera requisitos significativos de armazenamento de dados.
• Cisalhamento Externo (External Shear): A implementação atual não inclui cisalhamento externo da estrutura de grande escala fora do cone de luz, embora o efeito esperado seja estimado como pequeno ($\gamma_{ext} \sim 0.02$).

Os autores concluem que incorporar a informação completa do cone de luz é crucial para previsões teóricas e que considerar cuidadosamente a estrutura de LoS é necessário ao interpretar observações, pois a estrutura não correlacionada pode alterar significativamente as configurações de imagem e as propriedades das curvas críticas.