Evolution of Compact Stellar Systems in Ultralight Dark Matter Halos: Dependence on Stellar and Dark Matter Parameters

Resumo Técnico: Evolução de Sistemas Estelares Compactos em Halos de Matéria Escura Ultraleve

Enunciado do Problema
Sistemas estelares compactos, como galáxias anãs ultra-fracas (UFDs), são frequentemente utilizados para restringir a massa de partículas de matéria escura ultraleve (ULDM). A base teórica para essas restrições baseia-se no "efeito de aquecimento", onde a interferência de ondas no halo de ULDM transfere energia para sistemas estelares incorporados, fazendo com que eles se expandam e potencialmente contradigam dados observacionais. No entanto, a literatura recente apresenta conclusões conflitantes sobre a viabilidade dos modelos de ULDM. Alguns estudos sugerem que o arrancamento de maré (tidal stripping) pela Via Láctea (MW) suprime esse aquecimento, permitindo massas de ULDM mais baixas, enquanto outros argumentam que o aquecimento permanece significativo para massas mais altas. Além disso, o papel da dinâmica estelar interna — especificamente a relaxação de dois corpos e a evolução estelar (SE) — tem sido amplamente negligenciado. Este trabalho investiga o destino de longo prazo de populações estelares compactas em halos de ULDM, focando na interação entre a metalicidade estelar, o campo de maré da MW e a massa de partícula da ULDM, particularmente no regime onde o tamanho do sistema ($R_h$) é muito menor que o comprimento de onda de de Broglie da ULDM ($\lambda_{dB}$).

Metodologia
Os autores empregam uma estratégia de simulação de dois estágios usando simulações $N$-body diretas acopladas à evolução da função de onda da ULDM:

1. Evolução do Halo de ULDM: O halo de ULDM é modelado usando o método de decomposição de autovetores e evoluído via método pseudo-espectral (pacote PyUltraLight) em uma caixa cúbica de 40 kpc. O sistema estelar é inicialmente tratado como uma partícula pontual para registrar o tensor de maré da ULDM, $\nabla\nabla V_{\text{ULDM}}$, em sua posição.
2. Evolução do Sistema Estelar: A evolução interna do sistema estelar é simulada usando o código $N$-body direto PeTar. O sistema segue um perfil Plummer com $N=5 \times 10^4$ partículas e um raio de luz média inicial $R_h = 5$ pc.
3. Mecanismo de Acoplamento: A influência da ULDM é incorporada via uma aproximação de primeira ordem do tensor de maré. A aceleração de maré $\Delta a(r, r_c) \simeq -(r - r_c) \cdot \nabla\nabla V_{\text{ULDM}}(r_c)$ é adicionada à aceleração de cada partícula estelar em relação ao centro de massa do sistema.
4. Conjuntos de Simulação: O estudo conduz 11 execuções de simulação variando:
   • Metalicidade Estelar ($Z$): $10^{-4}$ vs. $10^{-2}$.
   • Evolução Estelar (SE) & IMF: SE completa com IMF, SE sem IMF, sem SE, e partículas de massa igual.
   • Contexto Ambiental: Evolução isolada vs. evolução dentro do campo de maré da MW (incluindo o arrancamento de maré do halo externo de ULDM).
   • Massa de Partícula de ULDM ($m_{22} \equiv m_a/10^{-22}$ eV): Valores variando de 1 a 3.
   • Dinâmica Orbital: Centro de massa fixo vs. movimento orbital livre.

Contribuições Principais e Resultados

• Impacto da Metalicidade e Evolução Estelar:
  O estudo descobre que a metalicidade estelar altera significativamente a resistência do sistema à desrupção. Sistemas com maior metalicidade ($Z=10^{-2}$) produzem buracos negros (BHs) de menor massa devido a ventos estelares mais fortes. Isso reduz a eficiência da segregação de massa, retardando a transferência de estrelas de baixa massa para as periferias do sistema e atrasando a desrupção. Inversamente, sistemas onde a SE é desligada exibem uma segregação de massa intensificada e desrupção rápida. Sistemas com partículas de massa igual (sem SE, sem IMF) não sofrem desrupção, mas em vez disso passam por colapso de núcleo.

• Papel do Campo de Maré da Via Láctea e Movimento Orbital:
  Contrariamente aos achados no regime $\lambda_{dB} \lesssim R_h$ (onde o arrancamento de maré da MW suprime o aquecimento ao remover estados de alta energia da ULDM), os autores descobrem que, no regime $\lambda_{dB} \gg R_h$, o campo de maré da MW aumenta a desrupção. O campo de maré da MW remove o halo externo de ULDM, alterando o movimento orbital do sistema estelar dentro do halo restante. Essa modificação aumenta a frequência de passagens pelo pericentro (aproximação mais próxima do sóliton), levando a eventos de aquecimento de maré mais frequentes e intensos. Sistemas que evoluem em isolamento (S4) permanecem estáveis por mais tempo do que aqueles sujeitos ao arrancamento de maré da MW (S0).

• Dominância do Aquecimento Induzido pelo Sóliton:
  Ao comparar simulações com flutuações completas de ULDM contra aquelas usando apenas um potencial de sóliton estático (S0-soliton, S4-soliton), os autores demonstram que o efeito de aquecimento é impulsionado principalmente pela influência de maré do sóliton, e não pelas estruturas granulares flutuantes do halo de ULDM. Uma perda significativa de partículas estelares ocorre quando o sistema estelar é deslocado além do raio do sóliton ($d > r_c$), onde ele experimenta o potencial gravitacional fortemente inhomogêneo do sóliton.

• Dependência da Massa de Partícula de ULDM:
  No regime $\lambda_{dB} \gg R_h$, o efeito de aquecimento torna-se mais forte com o aumento da massa da partícula de ULDM. Esta tendência é oposta à observada no regime $\lambda_{dB} \lesssim R_h$. Os autores atribuem isso à relação núcleo-halo do sóliton: conforme a massa da partícula aumenta, a massa do sóliton aumenta e seu raio diminui, levando a um campo de maré mais forte ($GM_c/r_c^3 \propto m_{22}^2$). Consequentemente, sistemas estelares dissolvem-se mais rapidamente para massas de partículas maiores neste regime.

Significância
O artigo conclui que a evolução de sistemas estelares compactos em halos de ULDM é altamente complexa e qualitativamente diferente dependendo da razão entre o tamanho do sistema e o comprimento de onda de de Broglie. As descobertas sugerem que as restrições existentes sobre a massa de partícula de ULDM derivadas de sistemas compactos (como Segue 1) podem exigir uma revisão cuidadosa. Especificamente, a interação entre a metalicidade estelar, processos de relaxação interna e a história orbital específica dentro do campo de maré da MW pode alterar significativamente os tempos de sobrevivência desses sistemas. O estudo enfatiza que os regimes $\lambda_{dB} \gg R_h$ e $\lambda_{dB} \lesssim R_h$ devem ser tratados separadamente ao derivar restrições astrofísicas sobre ULDM, pois os mecanismos físicos que impulsionam o aquecimento e a desrupção diferem fundamentalmente entre eles.