Cooling, conduction, compact objects: Gravothermal evolution of dissipative self-interacting dark matter halos

Resumo Técnico: Evolução Gravotérmica de Halos de Matéria Escura Autointeragente Dissipativa

Enunciado do Problema
Modelos de matéria escura autointeragente (SIDM) frequentemente predizem processos radiativos que dissipam energia, como o bremsstrahlung ou o decaimento de estados ligados colisionalmente excitados. Embora a SIDM elástica seja conhecida por impulsionar a evolução gravotérmica (formação de núcleo seguida de colapso), o impacto específico da dissipação de energia nesse processo permanece subexplorado em simulações numéricas. Métodos N-body existentes para espalhamentos frequentes de pequeno ângulo (fSIDM) tipicamente assumem interações elásticas, enquanto modelos dissipativos frequentemente dependem de aproximações de fluido ou suposições de seção de choque isotrópica (que podem não capturar a física de mediadores leves, que são frequentemente direcionados para a frente). Além disso, observações recentes de um objeto compacto no sistema de lente forte JVAS B1938+666 sugerem um halo de SIDM em colapso, mas modelos elásticos exigem seções de choque irrealisticamente grandes para explicar o colapso pelo redshift observado ($z=0.881$). Este trabalho investiga se as autointerações dissipativas podem alterar a evolução gravotérmica de halos isolados para explicar tais observações com seções de choque mais moderadas.

Metodologia
Os autores apresentam a primeira extensão do formalismo N-body de fSIDM (Fischer et al. 2021) para incluir a dissipação efetiva.

• Implementação Numérica: O método trata os espalhamentos frequentes de pequeno ângulo como uma força de arrasto e difusão de momento transversal. Os autores derivam uma extensão dissipativa onde a força de arrasto é aumentada por um parâmetro de dissipação $r_{\text{diss}}$, e o chute de velocidade estocástico é modificado para garantir que o sistema perca energia consistente com a energia irradiada $k_0$ no processo de física de partículas subjacente. O passo de tempo é reduzido por um fator de $r_{\text{diss}}$ para manter a estabilidade.
• Configuração da Simulação: Eles simulam halos NFW (Navarro-Frenk-White) isolados ($M \sim 10^{10} M_\odot$) usando o código OpenGadget3. Duas séries de simulações são realizadas:
  1. Variando a taxa de dissipação enquanto mantém a taxa de condução de calor constante.
  2. Variando a taxa de condução de calor (via seção de choque de transferência $\sigma_T/m_\chi$) enquanto mantém a taxa de resfriamento constante.
• Comparações: Os resultados são comparados contra:
  • Um modelo de fluido gravotérmico (aumentado com um termo de resfriamento) para avaliar a validade.
  • Simulações de "esparramento raro de autointeração" (rSIDM) isotrópicas para testar a influência da dependência angular na seção de choque.
• Aplicação Observacional: Os autores projetam os resultados das simulações para derivar perfis de densidade superficial e massas contidas, comparando-os com a distribuição de massa inferida para o perturbador de JVAS B1938+666 por Vegetti et al. (2026).

Principais Resultados

1. Mudança Qualitativa na Evolução: A dissipação não apenas acelera o colapso; ela altera qualitativamente a evolução. No caso elástico, a condução de calor eventualmente inverte a direção (para fora) à medida que o núcleo se torna isotérmico, impulsionando o colapso. No caso dissipativo, o resfriamento central suficientemente forte impede a formação de um núcleo isotérmico. O gradiente de dispersão de velocidade permanece positivo no centro, mantendo a condução de calor direcionada para dentro durante toda a evolução.
2. Influxo vs. Outfluxo de Massa: Em halos elásticos, o halo externo é aquecido pela condução, levando ao outflow de massa. Em halos dissipativos, as regiões externas resfriam eficientemente (dominadas pela dissipação sobre a condução externa), levando ao influxo de massa. Isso resulta em uma região maior de influxo de massa e uma indentação menos pronunciada entre o núcleo e o halo externo no perfil de densidade final.
3. Dependência de Parâmetros:
   • Dissipação: Aumentar a taxa de resfriamento (controlada por $r_{\text{diss}}$) acelera significativamente o colapso (por quase duas ordens de magnitude para os parâmetros testados) e produz núcleos mais densos e menos estendidos.
   • Condução: O efeito da seção de choque de transferência é não monotônico. Em seções de choque baixas (longo caminho livre médio), o aumento da condução acelera o colapso. Em seções de choque altas (curto caminho livre médio), o aumento da condução suprime a eficiência do transporte de calor, retardando o colapso.
   • Dependência Angular: Os autores descobrem que a evolução de halos isolados é amplamente insensível à dependência angular da seção de choque (direcionada para a frente vs. isotrópica), desde que as taxas de resfriamento e condução sejam correspondidas.
4. Aplicação JVAS B1938+666: Os autores demonstram que halos de SIDM fracamente dissipativos podem reproduzir o perfil de massa inferido para o perturbador de JVAS B1938+666. Comparado ao caso elástico, modelos dissipativos podem explicar a observação com:
   • Tempos de evolução significativamente mais curtos (permitindo o colapso até $z=0.881$ com seções de choque moderadas).
   • Ou, equivalentemente, seções de choque de autointeração menores para o mesmo tempo de evolução.
   • As seções de choque necessárias ($\sigma_T/m_\chi \sim 10^2 \text{ cm}^2 \text{ g}^{-1}$) são consistentes com as restrições para halos de massa $M_{200c} \lesssim 10^7 M_\odot$.

Significância e Alegações
O artigo afirma abrir uma "nova rota para conectar a estrutura de halos e objetos compactos recentemente relatados à microfísica do setor escuro". Ao estender o formalismo fSIDM para incluir a dissipação, os autores fornecem um arcabouço para testar modelos que apresentam mediadores leves e estados ligados. A principal significância reside em mostrar que a dissipação pode resolver a tensão entre o objeto compacto observado em JVAS B1938+666 e as grandes seções de choque exigidas pelos modelos elásticos. Os autores notam modestamente que suas conclusões dependem de condições iniciais específicas (halos NFW isolados) e que simulações cosmológicas e simulações que se estendam mais profundamente na fase de colapso (além da quebra numérica) são necessárias para restringir totalmente as propriedades dissipativas. Eles também reconhecem que a dependência angular da seção de choque pode se tornar mais relevante para satélites e fusões, um regime não coberto por este estudo de halo isolado.