The Sensitivity of Substructure Lensing to SIDM Core-collapse Model Variation

Resumo Técnico: A Sensibilidade do Lenteamento de Subestrutura à Variação do Modelo de Colapso do Núcleo de Matéria Escura Auto-Interagente (SIDM)

Declaração do Problema
O lenteamento gravitacional forte emergiu como uma sonda poderosa para a subestrutura de matéria escura, oferecendo um teste específico para modelos de Matéria Escura Auto-Interagente (SIDM). Em cenários SIDM com colisões elásticas, a transferência de calor leva à formação de núcleos isotérmicos, seguidos por uma "catástrofe gravotérmica" (colapso do núcleo) que produz núcleos centrais densos. Esses halos compactos são previstos para deixar marcas distintas nos observáveis de lenteamento. No entanto, a previsão precisa dessas assinaturas é dificultada pela dificuldade em modelar a linha do tempo do colapso do núcleo e o estado final dos subhalos. As abordagens atuais frequentemente dependem de suposições simplificadas, como tratar o colapso do núcleo como uma transição instantânea ou usar atribuições probabilísticas para estados de colapso. Os autores identificam uma desconexão entre as pipelines de inferência de lenteamento existentes e os modelos físicos subjacentes, observando que pequenas variações na forma como os subhalos em colapso do núcleo são modelados podem levar a diferenças significativas no poder de lenteamento. Sem quantificar essas incertezas de modelagem, a capacidade das observações de lenteamento de restringir a natureza da matéria escura é limitada.

Metodologia
Os autores quantificam esse desafio de modelagem testando várias modificações aos perfis de densidade de halos em colapso do núcleo SIDM no contexto do lenteamento de subestrutura. Sua metodologia envolve:

1. Simulação e Calibração: Eles utilizam dados de simulação $N$-corpos (gerados com o código Arepo e o módulo SIDM da Ref. [20]) para calibrar um tratamento aproximado dos efeitos de maré no colapso gravotérmico do núcleo. Eles derivam um "fator de aceleração de colapso" para contabilizar a diferença nas escalas de tempo de colapso entre halos isolados e subhalos evoluindo em um potencial hospedeiro.
2. Realizações de Lenteamento: Usando o pacote de código aberto pyHalo, eles geram 200 realizações aleatórias de sistemas de lentes em escala galáctica (em $z_{lens}=0.5$) contendo subhalos do defletor principal e halos da linha de visada (LOS). Eles substituem perfis padrão de Matéria Escura Fria (CDM) por perfis SIDM baseados em uma solução de fluido gravotérmico calibrada para simulações $N$-corpos.
3. Variações de Modelo: O estudo testa várias escolhas de modelagem:
   • Contínuo vs. Instantâneo: Comparando perfis de densidade evoluindo suavemente contra modelos binários onde os halos estão ou totalmente não-colapsados ou totalmente colapsados.
   • Limiares de Colapso: Variando o momento em que o perfil de densidade é "congelado" (parado) ou "impulsionado" para o estado final, testando diferentes estados finais do colapso.
   • Probabilístico vs. Determinístico: Comparando o rastreamento da evolução individual de halos contra atribuição aleatória de estados de colapso baseados em faixas de massa.
   • Perfis Iniciais: Testando halos não-colapsados modelados tanto como NFW Truncado padrão (TNFW) quanto como tendo atingido a fase de núcleo maximamente desenvolvido.
4. Análise Estatística: Eles calculam o potencial de lenteamento gravitacional efetivo de múltiplos planos, calculando a divergência ($\kappa_{div}$) e o rotacional ($\kappa_{curl}$) do campo de deflexão efetivo. A estatística primária é a função de correlação de dois pontos desses campos, decomposta em momentos multipolares (monopolo e quadrupolo) usando bases ortonormais para desvendar os efeitos da LOS do lente principal.

Principais Contribuições

• Quantificação da Incerteza de Modelagem: O artigo fornece uma quantificação sistemática de como suposições simplificadas específicas na modelagem do colapso do núcleo SIDM afetam as estatísticas de lenteamento.
• Calibração de Maré: Introduz um método probabilístico para contabilizar a aceleração/desaceleração por maré do colapso do núcleo em subhalos, calibrado contra simulações $N$-corpos.
• Comparação Binária vs. Contínua: Compara rigorosamente a aproximação de "colapso binário" (comum na literatura) contra um modelo gravotérmico contínuo e evolutivo no tempo, destacando onde a aproximação falha.
• Sensibilidade do Perfil: O estudo demonstra que a escolha do perfil de densidade inicial (TNFW vs. com núcleo) e o método de atribuição de estados de colapso alteram significativamente o sinal de lenteamento previsto.

Resultados

• Sensibilidade à Densidade vs. Tempo: A função de correlação de dois pontos (especificamente os momentos multipolares) é altamente sensível à densidade central dos halos colapsados, mas relativamente insensível à escala de tempo exata do colapso. Variações no tempo de colapso (por exemplo, parar a evolução em 20% do tempo total) não produzem resultados significativamente diferentes em comparação com a evolução completa, desde que o perfil de densidade final seja preciso.
• Supressão por Modelos Binários: Modelos de colapso do núcleo binário suprimem consistentemente o poder em pequena escala nas estatísticas de lenteamento em comparação com o modelo de evolução temporal completa. Essa supressão é mais severa quando:
  • Halos não-colapsados são modelados com perfis iniciais com núcleo em vez de perfis TNFW.
  • Estados de colapso são atribuídos aleatoriamente em vez de baseados na evolução individual do halo.
• Dependência de Escala: As diferenças entre os modelos são mais distinguíveis em pequenas escalas de comprimento (abaixo de 0,1 segundos de arco). O momento quadrupolo do campo de divergência ($\xi_{div,2}$) permanece sensível a variações do modelo em escalas maiores do que os momentos monopolo ou rotacional.
• Correlação entre Colapso e Poder de Lenteamento: O estudo encontra que halos com um efeito de lenteamento mais pronunciado (maior massa/concentração) tendem a colapsar mais rapidamente. A aleatorização dos estados de colapso suprime o aumento SIDM porque desacopla o estado de colapso das propriedades do halo que impulsionam o sinal de lenteamento.
• Validação $N$-corpos: O modelo aproximado de aceleração por maré dos autores produz estatísticas de lenteamento consistentes com simulações $N$-corpos dentro de um desvio padrão, validando seu uso como uma aproximação de primeira ordem para este estudo.

Significado e Alegações
O artigo alega que a função de correlação de dois pontos é uma sonda sensível das simplificações específicas feitas na modelagem do colapso do núcleo SIDM. Os autores afirmam que:

1. Robustez do Estado Final: Incertezas sistemáticas modestas no tempo do colapso do núcleo podem não afetar significativamente as previsões de lenteamento, desde que o perfil de densidade final do halo colapsado seja modelado com precisão.
2. Limitações das Aproximações Binárias: Embora os modelos de colapso binário sejam computacionalmente eficientes, eles podem introduzir vieses significativos ao suprimir o poder em pequena escala, particularmente se assumirem condições iniciais com núcleo ou aleatorizarem estados de colapso.
3. Necessidade de Motivação Física: Futuras análises de lenteamento visando testar SIDM devem fazer escolhas intencionais e motivadas fisicamente em relação aos perfis de subestrutura e halos da LOS. Os autores recomendam que tais trabalhos incluam testes análogos de suas escolhas de perfil e comparem seus perfis de ângulo de deflexão contra modelos de evolução completa motivados fisicamente.
4. Dependência Observacional: A detectabilidade dessas diferenças de modelagem depende fortemente da sonda observacional e das escalas sendo analisadas. O artigo não alega que dados atuais podem definitivamente descartar modelos SIDM específicos, mas sim que a interpretação desses dados depende fortemente da robustez dos modelos de subestrutura subjacentes.

Os autores concluem que o refinamento de modelos de colapso do núcleo SIDM em ambientes realistas é essencial, pois previsões testáveis podem depender fortemente do estado final dos halos com colapso do núcleo. Eles sugerem que sua estrutura pode ser usada para implementar esses modelos de evolução completa em previsões para anomalias de razão de fluxo e outros observáveis de lenteamento em trabalhos futuros.