Collapse versus Disruption: The Fate of Compact Stellar Systems in Ultralight Dark Matter Halos

Resumo Técnico: Colapso versus Disrupção em Halos de Matéria Escura Ultraleve

Enunciado do Problema
Modelos de matéria escura ultraleve (ULDM) preveem que a interferência das ondas de ULDM gera grânulos de densidade estocásticos, levando a flutuações temporais no potencial gravitacional. Essas flutuações induzem o aquecimento dinâmico em sistemas estelares incorporados, fazendo-os expandir. Embora este mecanismo tenha sido usado para estabelecer limites inferiores rigorosos na massa da partícula de ULDM ($m_a$) usando sistemas pequenos como aglomerados estelares nucleares e galáxias anãs ultra-fracas (UFDs), o comportamento evolutivo de sistemas significativamente menores que o comprimento de onda de de Broglie ($\lambda_{dB}$) da ULDM permanece mal compreendido. Especificamente, não está claro como o interplay entre o aquecimento induzido por ULDM e o relaxamento estelar interno (impulsionado pelo espalhamento gravitacional de dois corpos) afeta o destino de sistemas compactos. Estudos anteriores sugeriram que o aquecimento poderia ser ineficaz para sistemas muito pequenos, enquanto outros observaram que sistemas compactos possuem tempos de relaxamento mais curtos que a idade do Universo, potencialmente permitindo que o colapso do núcleo domine.

Metodologia
Os autores empregam uma estrutura de simulação numérica de dois estágios para resolver as escalas espaciais e temporais díspares envolvidas:

1. Estágio 1 (Evolução do Halo): O halo de ULDM é modelado como um todo, evoluindo dentro dos campos de maré da Via Láctea (MW) e da Grande Nuvem de Magalhães (LMC). O sistema estelar é tratado como uma massa pontual ($M_\star = 10^3 M_\odot$) para registrar sua trajetória e o tensor de maré dependente do tempo induzido pelas inhomogeneidades da ULDM. A simulação adota uma história de montagem realista onde o halo evolui fora da MW por $\sim 3$ Gyr antes de entrar em uma trajetória similar à da UFD Segue 1.
2. Estágio 2 (Evolução Interna): Usando o código de $N$-corpos direto PeTar, a evolução interna do sistema estelar é simulada sob o tensor de maré dependente do tempo derivado do Estágio 1. O sistema estelar é inicializado com um perfil Plummer em equilíbrio dinâmico aproximado.
   • Parâmetros: Seis sistemas representativos (S1–S6) são simulados com variando massa estelar total ($M_\star$), número de partículas ($N$) e raios de meia-massa ($R_h$).
   • Comparação: Cada sistema é evoluído tanto com o aquecimento induzido por ULDM quanto isoladamente (sem ULDM) para isolar o efeito de aquecimento.
   • Física: A aceleração da ULDM é aproximada via o tensor de maré atuando sobre as partículas estelares em relação ao centro de massa do sistema. As simulações rodam por 13 Gyr.

Contribuições Principais e Resultados
O estudo revela que a evolução de sistemas estelares compactos em halos de ULDM é governada por uma competição entre o relaxamento estelar interno e o aquecimento induzido por ULDM, levando a três regimes evolutivos distintos:

1. Regime Dominado pelo Relaxamento (Colapso do Núcleo):

   • Em sistemas suficientemente compactos (ex: S1, S3, S4), o relaxamento interno de dois corpos impulsiona o colapso do núcleo.
   • Achado Contraintuitivo: O aquecimento da ULDM não impede o colapso; em vez disso, ele o acelera. Ao remover estrelas externas, o aquecimento da ULDM reduz a massa estelar ligada e o raio de meia-massa. Essa redução encurta a escala de tempo de relaxamento ($t_{rh}$) e reduz a dispersão de velocidade, aumentando a eficiência do relaxamento de dois corpos. Consequentemente, esses sistemas permanecem ligados e tornam-se ainda mais concentrados do que estariam isoladamente.

2. Regime Dominado pelo Aquecimento (Disrupção):

   • Em sistemas mais estendidos (ex: S2, S5), a densidade estelar inicial é baixa o suficiente para que o tempo de relaxamento seja longo. Aqui, o aquecimento da ULDM domina, injetando energia que impulsiona o sistema em direção a um estado de dispersão.
   • Além de um limiar crítico, o sistema é completamente disrompido.

3. O Análogo da "Segue 1":

   • Os autores identificam um regime de quase-disrupção cujos remanescentes se assemelham às UFDs observadas, como a Segue 1.
   • A simulação S5, aproximando-se da disrupção em $t \approx 10$ Gyr, produz perfis de densidade superficial estelar e dispersão de velocidade que coincidem amplamente com os dados observacionais da Segue 1.
   • Um sistema com massas estelares mais realistas (S6) também reproduz propriedades semelhantes à Segue 1 em uma idade mais jovem ($\sim 6,5$ Gyr). Isso sugere que, para $m_a = 10^{-22}$ eV, a Segue 1 poderia ser interpretada como um sistema atualmente próximo da disrupção.

Diagrama de Fase e Quantificação
Para generalizar essas descobertas, os autores introduzem um parâmetro adimensional $\alpha$, definido como a razão entre o tempo de relaxamento de meia-massa ($t_{rh}$) e o tempo de disrupção característico da ULDM ($t_{dis}$):
$$\alpha \equiv \frac{t_{rh}}{t_{dis}}$$

• $\alpha < 0,05$: O relaxamento domina, levando ao colapso do núcleo.
• $\alpha > 0,05$: O aquecimento domina, levando à disrupção.
• $t_{rh} > 4,7$ Gyr (e $t_{dis} > 85$ Gyr): Ambos os efeitos são negligenciáveis durante a idade do Universo; os sistemas permanecem estáveis.

Usando simulações extensas, os autores constroem um diagrama de fase bidimensional no plano $(N, M_\star R_h^{-3})$. Este diagrama separa claramente o espaço de parâmetros em regiões de colapso do núcleo, disrupção e evolução negligenciável, demonstrando que o destino evolutivo é determinado pelos parâmetros iniciais do sistema ($M_\star, N, R_h$) e pela massa da partícula de ULDM.

Significância e Alegações
O artigo afirma que a dinâmica de sistemas estelares compactos em halos de ULDM é "rica e não trivial", desafiando a visão simplista de que o aquecimento da ULDM sempre leva à expansão ou que sistemas pequenos são imunes aos efeitos de aquecimento.

• Restrições Robustas: Os autores argumentam que a modelagem precisa do interplay entre relaxamento e aquecimento é essencial para derivar restrições robustas sobre a ULDM. Ignorar a aceleração do colapso do núcleo ou os limiares específicos de disrupção poderia levar a limites incorretos para a massa da partícula de ULDM.
• Interpretação Observacional: Os resultados fornecem um cenário evolutivo plausível para UFDs como a Segue 1, sugerindo que elas podem ser sistemas em um estado transitório próximo da disrupção, em vez de objetos em equilíbrio estático.
• Limitações: Os autores declaram explicitamente que sua análise é realizada em uma massa de partícula de ULDM e configuração de halo fixas. Eles observam que incorporar física adicional, como uma função de massa inicial (IMF) estelar realista e evolução estelar, complicaria ainda mais o quadro e é necessário para estabelecer conexões mais robustas com as observações.