Growth of Structure in Multi-species Wave Dark Matter

Resumo Técnico: Crescimento de Estrutura em Matéria Escura Ondulatória Multi-espécie

Declaração do Problema
A natureza do setor escuro permanece desconhecida, com possibilidades que vão desde um único campo bosônico (por exemplo, o áxion da QCD) até um conjunto multi-espécie contendo componentes bosônicos e fermiônicos. Embora trabalhos anteriores tenham abordado o crescimento de estrutura em matéria escura ondulatória de espécie única ou em matéria escura de partículas multi-espécie, faltava um quadro abrangente para calcular o crescimento gravitacional de estrutura em um sistema multi-espécie que compreende tanto componentes ondulatórios (campos bosônicos) quanto componentes tipo partícula. Tais sistemas podem surgir em teorias de alta energia como o "universo de áxions das cordas" (string axiverse), onde diferentes campos possuem massas distintas, frações de densidade e condições iniciais (por exemplo, escalares globalmente desalinhados versus campos localmente desalinhados ou componentes de fóton escuro vetorial). O desafio reside em modelar como esses componentes diversos, potencialmente exibindo "calor" (movimento aleatório) e flutuações intrínsecas semelhantes a Poisson, interagem gravitacionalmente para moldar o espectro de potência da matéria.

Metodologia
Os autores desenvolvem um quadro teórico baseado no sistema Schrödinger–Poisson para $N$ campos escalares não relativísticos $\phi_s(t, \mathbf{x})$. A evolução é analisada usando transformadas de Fourier dos campos e seus conjugados complexos, tratando-os como variáveis aleatórias com estatísticas iniciais especificadas.

As etapas metodológicas-chave incluem:

1. Definição Bilinar: Os autores definem uma quantidade bilinear $\hat{f}^s_k(p)$ relacionada à distribuição no espaço de fase, permitindo que a perturbação fracional de densidade $\delta$ seja expressa como uma integral sobre o espaço de momento.
2. Equações de Evolução: Ao aplicar o sistema Schrödinger–Poisson, eles derivam equações de evolução para as funções de correlação dessas quantidades bilineares. Essas equações accountam para o acoplamento entre diferentes espécies via potencial gravitacional.
3. Decomposição de Perturbações: O espectro de potência total de densidade é decomposto em duas contribuições distintas:
   • Perturbações Adiabáticas: Originadas de perturbações de curvatura primordiais, evoluindo de forma semelhante à Matéria Escura Fria (CDM) padrão, mas modificadas pela dinâmica multi-espécie.
   • Perturbações Isocurvatura: Originadas de flutuações intrínsecas semelhantes a Poisson (ruído de disparo) inerentes à natureza ondulatória dos campos, particularmente relevantes para mecanismos de produção pós-inflacionária ou campos localmente desalinhados.
4. Equações de Volterra: A evolução temporal do espectro de potência é governada por um conjunto de equações integrais de Volterra acopladas. Essas equações envolvem "kernels de livre-escoamento" ($T^{fs}$) que dependem dos espectros de potência iniciais dos campos ($f^s_0(p)$) e da massa de cada espécie. Os kernels capturam efeitos como livre-escoamento, pressão efetiva decorrente do movimento aleatório ("calor") e "pressão quântica" dinâmica ondulatória.
5. Implementação Numérica: Os autores fornecem um algoritmo numérico para resolver essas equações integrais de Volterra acopladas de forma iterativa. Um código público é lançado para calcular espectros de potência de densidade para números arbitrários de espécies, frações de densidade e espectros de potência iniciais de campos.

Contribuições e Resultados Principais
O artigo deriva expressões gerais para o espectro de potência do contraste de densidade total $P_\delta(k)$ e os espectros cruzados entre diferentes espécies ($P^{ss'}_\delta(k)$).

• Formalismo Geral: O espectro de potência total é expresso como a soma de um termo adiabático e um termo de isocurvatura, cada um evoluído via funções de transferência específicas ($T^{(ad)}_k$ e $T^{(iso)}_k$). Essas funções de transferência são soluções de equações de Volterra que incorporam a soma ponderada dos kernels de livre-escoamento de todas as espécies.
• Escala Relevantes: Os autores identificam e analisam várias escalas críticas que suprimem o crescimento de estrutura:
  • Escala de Livre-escoamento ($k_{fs}$): Apaga perturbações adiabáticas iniciais.
  • Escala de Jeans "Quente" ($k_J$): Suprime o crescimento devido à pressão efetiva decorrente do movimento aleatório.
  • Escala de Jeans "Fuzzy" ($k_j$): Suprime o crescimento devido à pressão quântica dinâmica ondulatória.
  • Escala de Isocurvatura Inicial ($k_*$): A escala na qual o espectro de potência inicial semelhante a Poisson atinge o pico.
• Interação de Componentes: O quadro demonstra que um componente ondulatório quente subdominante pode influenciar significativamente o espectro de potência total.
  • Potencialização de Isocurvatura: Mesmo que um componente seja subdominante em massa, suas flutuações intrínsecas de isocurvatura podem semear perturbações no componente dominante, potencialmente aumentando o espectro de potência total em pequenas escalas.
  • Mecanismos de Supressão: A combinação de livre-escoamento, supressão de Jeans quente e supressão de Jeans fuzzy leva a características complexas e dependentes da escala no espectro de potência. Por exemplo, em cenários com um escalar globalmente desalinhado dominante e um componente ondulatório quente subdominante, o espectro total mostra supressão em alto $k$ devido à pressão fuzzy do componente dominante, enquanto as sementes de isocurvatura do componente subdominante geram crescimento que é subsequentemente suprimido pelo seu próprio calor.
• Limites: O formalismo recupera naturalmente limites conhecidos:
  • Matéria Escura Ondulatória Globalmente Desalinhada: Corresponde a $k_* \to 0$, resultando em isocurvatura e calor negligenciáveis, comportando-se como matéria escura fuzzy padrão.
  • Matéria Escura de Partículas: Corresponde a $k_* \to \infty$ com dispersão de velocidade fixa, recuperando os resultados padrão de matéria escura de partículas quente/fria.
• Exemplos Numéricos: O artigo apresenta resultados numéricos para cenários de dois componentes (por exemplo, CDM + matéria escura ondulatória quente, ou matéria escura ondulatória globalmente desalinhada + matéria escura ondulatória quente). Esses exemplos ilustram como variar a fração de massa do componente quente altera o espectro de potência, mostrando que componentes subdominantes podem deixar marcas distintas no espectro de potência da matéria em escalas acessíveis a florestas de Lyman-$\alpha$ e levantamentos de galáxias.

Significado
O artigo afirma fornecer um quadro necessário e geral para estudar a formação de estrutura em cenários complexos de matéria escura que não podem ser reduzidos a um único fluido efetivo. Ao unificar cálculos dinâmicos ondulatórios com interações multi-espécie, permite a investigação de:

• Cenários onde a matéria escura consiste em múltiplos campos com diferentes massas e mecanismos de produção (por exemplo, o universo de áxions das cordas).
• O impacto de flutuações intrínsecas de Poisson (isocurvatura) no crescimento de estrutura, o que pode ser significativo mesmo para componentes subdominantes.
• A interação entre diferentes mecanismos de supressão (livre-escoamento, pressão térmica, pressão quântica) em sistemas multi-componentes.

Os autores enfatizam que seus resultados são gerais o suficiente para incluir vários candidatos a matéria escura, desde escalares ultraleves até campos vetoriais e partículas clássicas. Eles fornecem um código publicamente disponível para facilitar o cálculo de espectros de potência para esses modelos complexos, permitindo comparações com dados observacionais para restringir a natureza do setor escuro. O trabalho serve como base para futuros estudos sobre formação de estrutura não linear e as assinaturas observacionais da matéria escura multi-espécie.