An Enhanced Isothermal Jeans Approach to Constraining Dark Matter Self-Interactions from Galactic Kinematics

Resumo Técnico: Uma Abordagem Isotérmica de Jeans Aperfeiçoada para Restringir Interações de Auto-espalhamento da Matéria Escura a partir da Cinemática Galáctica

Definição do Problema
O paradigma cosmológico padrão $\Lambda$CDM enfrenta desafios em escalas subgalácticas, particularmente em relação à diversidade estrutural de galáxias anãs. Embora o feedback bariônico possa explicar algumas discrepâncias, a Matéria Escura com Auto-interação (SIDM) oferece um mecanismo alternativo onde as partículas de matéria escura sofrem espalhamento próprio frequente, impulsionando uma evolução de dois estágios: a formação inicial de um núcleo seguida pelo potencial colapso gravotérmico. Tentativas anteriores de restringir a seção de choque de auto-espalhamento da SIDM usando curvas de rotação galáctica foram limitadas pelos custos computacionais (dependendo de simulações completas de $N$-corpos ou de fluidos 1D) ou pela incapacidade de modelos semi-analíticos de capturar a fase de colapso do núcleo. Além disso, implementações anteriores do modelo de Jeans isotérmico frequentemente excluíam sistemas com centros dominados por bárions devido a problemas de robustez numérica.

Metodologia
Os autores apresentam um modelo semi-analítico melhorado baseado na abordagem de Jeans isotérmica, construído sobre estruturas estabelecidas por Jiang et al. (2023) e Yang et al. (2024). A metodologia envolve três atualizações técnicas fundamentais:

1. Seções de Choque Dependentes da Velocidade: O modelo incorpora o espalhamento de Rutherford elástico via um potencial de Yukawa, substituindo a suposição de seção de choque constante por uma seção de choque efetiva ($\sigma_{\text{eff}}$) que depende da velocidade circular máxima ($V_{\text{max}}$) do halo. Isso é parametrizado por uma normalização de baixa velocidade ($\sigma_0$) e uma escala de velocidade ($\omega$).
2. Tratamento Empírico do Colapso do Núcleo: Para abordar a quebra da suposição isotérmica durante o colapso gravotérmico, os autores adotam um contorno empírico. Eles utilizam a solução de alta densidade das equações de Jeans, que é tipicamente descartada, e a espelham em torno do tempo de fusão ($t_{\text{merge}}$) onde as soluções de baixa e alta densidade convergem. Isso permite que o modelo aproxime perfis de densidade pós-colapso.
3. Robustez Numérica Aperfeiçoada: Os autores refinam o procedimento de "costura" (combinando o núcleo isotérmico com o halo externo do tipo CDM) implementando limiares de convergência rigorosos ($\delta^2 = 10^{-5}$) e particionamento de intervalos de busca. Isso permite a identificação confiável de ambas as soluções de baixa densidade (crescimento do núcleo) e de alta densidade (colapso do núcleo), mesmo para galáxias com centros dominados por bárions.

O modelo é aplicado a uma amostra de 68 galáxias do banco de dados SPARC. O processo de ajuste envolve:

• Modelagem das distribuições bariônicas (gás e estrelas) usando um perfil de Hernquist, com razões massa-luminosidade estelar derivadas de ajustes de múltiplos perfis anteriores.
• Uso de métodos de Cadeia de Markov Monte Carlo (MCMC) para inferir a massa virial do halo ($M_{\text{vir}}$), concentração ($c$) e o produto degenerado $p \equiv \sigma_{\text{eff}} t_{\text{age}}$.
• Conversão de $p$ para $\sigma_{\text{eff}}$ estimando as idades dos halos via teoria de Press–Schechter estendida.
• Realização de uma análise conjunta para restringir os parâmetros globais $\sigma_0$ e $\omega$ minimizando o desacordo entre as previsões do modelo e as curvas de rotação observadas.

Principais Contribuições e Resultados

• Expansão da Amostra: A estabilidade numérica melhorada permite a inclusão de galáxias com centros dominados por bárions, expandindo significativamente o tamanho da amostra em comparação com estudos anteriores que utilizaram métodos semi-analíticos similares.
• Bimodalidade e Estados Evolutivos: Aproximadamente 1/6 da amostra (11 galáxias) exibe posteriors bimodais, onde tanto uma solução de crescimento de núcleo quanto uma solução de colapso de núcleo fornecem ajustes igualmente bons. Essa degenerescência surge porque halos em diferentes estágios evolutivos podem possuir densidades centrais quase idênticas. As galáxias restantes estão predominantemente na fase de crescimento do núcleo ou próximas ao estágio de núcleo máximo.
• Restrições Dependentes da Velocidade: Restrições conjuntas através da amostra revelam uma clara dependência de velocidade no espaço de parâmetros permitido. Os resultados formam uma degenerescência em forma de "L" no plano $(\sigma_0, \omega)$, indicando que dois regimes distintos são viáveis:
  • Seções de choque baixas, quase constantes: $\sigma_0 \sim 2 \, \text{cm}^2/\text{g}$ com $\omega \gtrsim 500 \, \text{km/s}$.
  • Modelos fortemente dependentes da velocidade: $\sigma_0 \sim 100 \, \text{cm}^2/\text{g}$ com $\omega \sim 60 \, \text{km/s}$.
  • Adotando a interpretação de crescimento do núcleo, os valores de melhor ajuste são $\sigma_0 \simeq 5 \, \text{cm}^2/\text{g}$ e $\omega \simeq 250 \, \text{km/s}$.
• Comparação com CDM: Quando comparado aos modelos de Matéria Escura Fria (CDM), o modelo SIDM supera perfis simples de contração adiabática e rivaliza com perfis de CDM baseados em feedback empírico (especificamente o modelo de Di Cintio et al. 2014) no ajuste de curvas de rotação. No entanto, ao contrário dos modelos de feedback bariônico, as seções de choque de SIDM inferidas não mostram correlação significativa com a razão massa estelar-halo, sugerindo que a SIDM fornece um mecanismo distinto para a diversidade estrutural.
• Relações de Escala: O estudo descobre que os modelos de massa SIDM resultam em concentrações de halo sistematicamente menores do que as relações padrão de CDM (ex: Dutton & Macciò 2014). Além disso, embora algumas galáxias nos ajustes de SIDM exibam eficiências de formação estelar elevadas (aproximando-se da unidade), essa tendência também é observada em estudos anteriores baseados em CDM, indicando que a SIDM não resolve unicamente as tensões na relação massa estelar-halo.

Significância
O artigo afirma que o paradigma SIDM permanece um quadro convincente para explicar a diversidade da estrutura em pequena escala. As restrições derivadas são notavelmente consistentes com sondas independentes, incluindo cinemática estelar em anãs ultra-fracas, lentes fortes e as funções de massa de buracos negros supermassivos no universo primordial (especificamente os "Little Red Dots"). Os autores enfatizam que, embora a degenerescência em forma de "L" impeça um modelo único de seção de choque que se ajuste a toda a amostra de galáxias, os resultados destacam uma região específica do espaço de parâmetros que satisfaz simultaneamente restrições ortogonais de cinemática galáctica local e demografia de buracos negros de alto redshift. O estudo conclui que futuros dados cinemáticos de baixa massa, particularmente para anãs ultra-fracas ($V_{\text{max}} \lesssim 30 \, \text{km/s}$), serão cruciais para quebrar as atuais degenerescências nos modelos de seção de choque dependentes da velocidade.