Non-linear Evolution of Dark Plasma Subhalos

Resumo Técnico: Evolução Não Linear de Subhalos de Plasma Escuro

Definição do Problema
Este artigo investiga a evolução dinâmica de subhalos de matéria escura (DM) orbitando dentro do potencial da Via Láctea (MW) sob a premissa de que a DM constitui um "plasma escuro". Em modelos onde a DM interage consigo mesma através de forças de longo alcance mediadas por um fóton oculto (um setor $U(1)$ escuro), a dinâmica assemelha-se à dos plasmas do Modelo Padrão. Especificamente, quando um subhalo de DM se move através do plasma escuro ambiente da MW, a velocidade de deriva relativa entre as duas populações pode exceder suas dispersões de velocidade internas. Esta configuração é suscetível a instabilidades de plasma (por exemplo, instabilidades de dois fluxos e de feixe-plasma), que geram campos eletrostáticos e eletromagnéticos escuros. Esses campos espalham partículas, modificam a distribuição no espaço de fase e conduzem o sistema ao equilíbrio térmico. O problema central abordado é a evolução não linear dessas instabilidades e seu impacto na perda de massa dos subhalos, comparando especificamente o aquecimento induzido pelo plasma contra a remoção por maré gravitacional padrão.

Metodologia
Os autores empregam uma abordagem híbrida semianalítica e numérica para modelar a evolução do subhalo:

1. Estrutura Semi-Analítica:

   • Condições Iniciais: Os subhalos são modelados como plasmas esféricos de neutralidade líquida com perfis de densidade NFW e distribuições de velocidade de Maxwell-Boltzmann. Duas massas de referência são consideradas: $10^7 M_\odot$ e $10^9 M_\odot$.
   • Dinâmica Orbital: As órbitas dos subhalos são integradas usando o pacote galpy dentro de um potencial da MW (MWPotential2014). As órbitas são inicializadas no raio virial com diferentes excentricidades ($\epsilon = 0,1$ a $0,9$).
   • Acoplamento de Perda de Massa: O perfil de densidade do subhalo é discretizado em camadas radiais concêntricas. A evolução de cada camada é rastreada ao longo de passos de tempo ($\delta t = 7$ Myr) até $1,2 \times t_{\text{peri}}$ (tempo da primeira passagem pelo pericentro).
   • Mecanismos: A perda de massa é calculada via dois canais distintos:
     • Remoção por Maré (Tidal Stripping): Modelada usando um parâmetro de eficiência padrão ($\alpha = 0,55$) aplicado às camadas fora do raio de maré instantâneo.
     • Aquecimento por Plasma: A fração de partículas ligadas ($f_{\text{bound}}$) em cada camada é determinada pela interpolação de resultados de um conjunto de simulações de Particle-in-Cell (PIC). Partículas aquecidas acima da velocidade de escape local ($v_{\text{esc}} \approx 4\sigma_{\text{SH}}$) são removidas, e a massa restante ligada se revirializa.

2. Simulações de Particle-in-Cell (PIC):

   • Código: O estudo utiliza o código TRISTAN-MP, adaptado para plasmas escuros via invariância de escalonamento de coordenadas das equações de Vlasov-Poisson.
   • Configuração: Simulações eletrostáticas unidimensionais (1D) são realizadas com condições de contorno periódicas. A configuração modela feixes de contraviajem representando o subhalo e o fundo da MW.
   • Parâmetros: As simulações variam a razão de densidade ($\rho_{\text{SH}}/\rho_{\text{MW}} \in [10^{-3}, 10^2]$) e a velocidade de streaming relativa ($v_{\text{SH}}/\sigma_{\text{MW}}$). A razão carga-massa de referência é $q_\chi/m_\chi = 10^{-14} e/m_p$.
   • Dimensionalidade: Embora os resultados primários sejam 1D, os autores realizam uma simulação 2D-2V limitada (duas dimensões espaciais, duas de velocidade) para verificar que a eficiência de aquecimento aumenta em dimensões superiores, sugerindo que os resultados 1D são conservadores.

Principais Contribuições

• Primeiro Estudo Não Linear: Este trabalho apresenta o primeiro estudo da evolução não linear de instabilidades eletrostáticas em subhalos de DM utilizando simulações cinéticas PIC completas.
• Quantificação do Aquecimento por Plasma: Os autores quantificam a eficiência da perda de massa induzida pelo plasma, demonstrando que esta pode exceder significativamente a remoção por maré, particularmente para subhalos de baixa massa em órbitas excêntricas.
• Modelo Semi-Analítico: Um novo framework é desenvolvido que acopla o crescimento da instabilidade de plasma local (derivado de dados PIC) com a evolução orbital global para prever a evolução da massa e do perfil de densidade do subhalo.

Resultados

• Eficiência de Perda de Massa: Para subhalos em órbitas altamente excêntricas ($\epsilon = 0,9$), o aquecimento por plasma é o mecanismo dominante de perda de massa. No momento do primeiro pericentro ($t = 1,2 t_{\text{peri}}$):
  • Um subhalo de $10^7 M_\odot$ perde $\sim 97\%$ de sua massa inicial.
  • Um subhalo de $10^9 M_\odot$ perde $\sim 84\%$ de sua massa inicial.
  • Em contraste, a remoção por maré sozinha remove apenas $\sim 30\%$ da massa para ambos os casos.
• Dependência Radial ("Remoção de Fora para Dentro"): A perda de massa procede de fora para dentro. As instabilidades crescem mais eficientemente onde o contraste de densidade entre o subhalo e a MW é modesto ($\rho_{\text{SH}}/\rho_{\text{MW}} \sim 1$) e a velocidade relativa é alta. Consequentemente, as camadas externas são aquecidas e removidas antes do núcleo interno. Para órbitas altamente excêntricas, o núcleo é eventualmente afetado próximo ao pericentro; para órbitas circulares, a perda de massa é confinada às regiões externas.
• Implicações na Função de Massa: A remoção preferencial de subhalos de baixa massa (devido às suas menores dispersões de velocidade interna e velocidades de escape) sugere uma supressão da extremidade de baixa massa da Função de Massa de Subhalos (SHMF). Os autores estimam que, para os parâmetros de referência, a SHMF poderia exibir uma inclinação mais rasa ou um corte abrupto abaixo de $M_{\text{SH}} \sim 10^7 M_\odot$.
• Excentricidade Orbital: O efeito é altamente sensível à excentricidade. Para órbitas circulares ($\epsilon = 0,1$), o aquecimento por plasma remove apenas $\sim 50\%$ ($10^7 M_\odot$) e $\sim 30\%$ ($10^9 M_\odot$) da massa, valores comparáveis ou ligeiramente superiores aos efeitos de maré.

Significância e Alegações
O artigo afirma que a física do plasma escuro oferece um mecanismo distinto para a evolução de subhalos que pode deixar assinaturas observáveis na população de satélites da MW. Especificamente, os autores argumentam que o aquecimento e a remoção induzidos pelo plasma podem:

1. Remodelar a SHMF: Ao destruir preferencialmente subhalos de baixa massa e "frios", modelos de plasma escuro podem resolver ou alterar o problema dos "satélites faltantes" ou o problema do "demasiado grande para falhar" (too big to fail), potencialmente levando a uma inclinação suprimida na baixa massa da população observada de subhalos.
2. Fornecer Restrições Observacionais: A supressão de subhalos de baixa massa e a resultante modificação da SHMF fornecem uma nova via para restringir parâmetros do setor escuro (especificamente a razão carga-massa $q_\chi/m_\chi$) usando dados de levantamentos como ELVES, SPARC, SAGA e Merian.
3. Destacar a Dinâmica Não Linear: O estudo enfatiza que efeitos coletivos de plasma, impulsionados por instabilidades de streaming, podem operar em escalas de tempo comparáveis aos tempos dinâmicos galácticos, levando a uma evaporação rápida que simulações padrão de CDM sem colisões não capturam.

Os autores mantêm um tom modesto quanto ao escopo de seus achados, observando que seus resultados baseiam-se em simulações 1D (embora testes 2D sugiram que o 1D é conservador) e que uma previsão completa ao nível populacional requer trabalhos futuros incorporando gravidade autoconsistente e efeitos de plasma em dimensões superiores.