Core-halo scaling relations in self-interacting scalar field dark matter

Declaração do Problema
Matéria Escura de Campo Escalar (SFDM), frequentemente referida como Matéria Escura Fuzzy (FDM) quando as auto-interações são negligenciadas, postula que a matéria escura consiste em bósons ultraleves ($m \sim 10^{-22}$ eV/$c^2$). Embora os modelos padrão de FDM explicem com sucesso a formação de núcleos solitônicos em halos galácticos, o impacto das auto-interações (SI) na estrutura e evolução desses halos permanece uma questão em aberto. A literatura anterior estabeleceu relações de escala para FDM não interagentes (por exemplo, relações massa-núcleo), mas não está claro como auto-interações repulsivas ou atrativas modificam essas relações, particularmente no que diz respeito à estabilidade dos núcleos, à relação massa-raio e ao potencial de colapso gravitacional em buracos negros supermassivos. Além disso, embora o regime de Thomas-Fermi (TF) de SFDM fortemente repulsivo tenha sido estudado usando aproximações de fluido, simulações numéricas diretas de fusões de solitons nesse regime sem aproximações de fluido têm sido escassas.

Metodologia
Os autores realizam simulações numéricas tridimensionais de múltiplas fusões de solitons para investigar halos de SFDM em uma ampla gama de intensidades e sinais de auto-interação.

• Equações Governantes: O estudo resolve as equações de Gross-Pitaevskii-Poisson (GPP) não relativísticas, que descrevem a evolução da função de onda do campo escalar $\psi$ acoplada à gravidade. As equações incluem um termo de auto-interação cúbico (proporcional ao comprimento de espalhamento $a_s$) e um termo quártico de ordem superior.
• Esquema Numérico: Os autores utilizam um método pseudo-espectral de divisão temporal otimizado para GPUs. Essa abordagem trata os termos de auto-interação como um potencial efetivo, permitindo o uso de um esquema de chute-deriva-chute (KDK). Para interações atrativas, um método KDK de ordem superior (O(4)) é implementado para garantir estabilidade.
• Configuração da Simulação: As simulações envolvem a fusão de $N_{sol}$ solitons (variando de 10 a 50) dentro de uma caixa de simulação estática e periódica de $L=100$ kpc. A massa da partícula é fixada em $m = 10^{-22}$ eV/$c^2$.
• Espaço de Parâmetros: O estudo explora três regimes distintos:
  1. Livre (FDM): $a_s = 0$.
  2. SI Repulsiva: $a_s > 0$, incluindo o regime de Thomas-Fermi fortemente repulsivo ($\Lambda \gg 1$).
  3. SI Atrativa: $a_s < 0$, até valores onde o sistema se aproxima do colapso gravitacional.
• Análise: Os autores analisam perfis de densidade, massa do núcleo ($M_c$), raio do núcleo ($r_c$) e energia total ($E$) para derivar relações de escala. Eles comparam os resultados das simulações com previsões analíticas derivadas da minimização de energia e soluções auto-semelhantes.

Principais Contribuições e Resultados

1. Estrutura e Densidade do Núcleo:

   • SI Repulsiva: Leva à formação de núcleos mais massivos e estendidos com densidades centrais mais baixas em comparação com o caso de FDM livre. A interação repulsiva fornece suporte de pressão adicional, impedindo que o núcleo se torne tão compacto quanto no caso não interagentes.
   • SI Atrativa: Aumenta as densidades centrais. Para interações atrativas suficientemente fortes, a densidade do núcleo diverge, indicando um colapso gravitacional. Os autores identificam uma escala de massa crítica ($M_{max}$) para modelos atrativos; se a massa do núcleo exceder esse limiar, o equilíbrio é perdido, potencialmente levando à formação de buracos negros supermassivos.
   • Evolução: Núcleos repulsivos permanecem estáveis e estendidos ao longo da simulação ($t = 10\tau_{dyn}$), enquanto núcleos atrativos mostram um aumento monótono na densidade central até o colapso.

2. Relações de Escala Massa-Raio e Núcleo-Halo:

   • Relação Massa-Raio: Os resultados da simulação para a relação massa-raio do núcleo ($M_c$ vs. $r_c$) alinham-se bem com previsões analíticas (Chavanis 2011; Padilla et al. 2021) mesmo na presença de fortes auto-interações. Interações repulsivas deslocam a relação para cima (permitindo mais massa para um dado raio), enquanto interações atrativas deslocam-na para baixo.
   • Relação Massa-Núcleo-Halo: Os autores investigam a relação entre a razão massa-núcleo/massa-halo ($M_c/M$) e a quantidade invariante $\Xi = |E| / [M^3(Gm/\hbar^2)]$. Eles constatam que, embora a inclinação dessa relação permaneça relativamente universal (semelhante ao caso FDM $M_c/M \propto \Xi^{1/3}$), a normalização depende da intensidade da auto-interação. Isso sugere que o mecanismo que alimenta massa ao núcleo é primordialmente gravitacional e universal uma vez que os halos são virializados, enquanto a SI determina localmente o suporte de pressão e o raio do núcleo.
   • Relação Tamanho-Energia: A relação entre o tamanho do núcleo ($r_c$) e a energia total ($E$) é sensível à SI. Interações repulsivas resultam em inclinações mais suaves (núcleos maiores para uma dada energia), enquanto interações atrativas resultam em inclinações mais íngremes.

3. Decomposição de Energia:

   • Nos casos repulsivos, a energia de auto-interação domina a região do núcleo, explicando a menor densidade e o maior tamanho.
   • Nos casos atrativos e livres, a energia potencial gravitacional domina o núcleo.
   • Os autores confirmam que as relações de escala não são estritamente universais, mas dependem do estágio evolutivo do halo e dos parâmetros específicos de SI.

Significado e Afirmações
O artigo afirma que as auto-interações fornecem um mecanismo natural para regular as propriedades dos núcleos de halos de SFDM. O estudo demonstra que:

• As relações de escala núcleo-halo derivadas para FDM livre são casos limites dentro de uma família mais ampla de modelos de SFDM que incluem auto-interações.
• Interações repulsivas podem estabilizar núcleos contra o colapso e produzir estruturas maiores e mais difusas, o que pode ter implicações para o problema "núcleo-cúspide" na formação de galáxias.
• Interações atrativas podem levar núcleos à instabilidade e colapso, oferecendo um canal potencial de formação para sementes de buracos negros supermassivos em centros galácticos.
• O método pseudo-espectral é validado como uma ferramenta robusta para simular SFDM no regime de Thomas-Fermi sem recorrer a aproximações de fluido, capturando com sucesso a transição do comportamento solitônico para o comportamento fluido.

Os autores concluem que, embora a dinâmica gravitacional fundamental do crescimento do núcleo seja universal, as propriedades estruturais específicas (densidade, tamanho, estabilidade) são altamente sensíveis ao sinal e magnitude das auto-interações e à história evolutiva do halo. Essas descobertas sugerem que restrições observacionais em núcleos galácticos poderiam ser usadas para restringir o comprimento de espalhamento e a massa da partícula de SFDM.