Cylindrical cosmological simulations with StePS

Resumo Técnico: Simulações Cosmológicas Cilíndricas com StePS

Enunciado do Problema
As simulações de $N$-corpos cosmológicas padrão tipicamente evoluem a matéria dentro de uma caixa cúbica periódica, impondo uma topologia de toro triplo ($T^3 = S^1 \times S^1 \times S^1$). Embora numericamente convenientes e preservando a homogeneidade estatística, esta topologia quebra a isotropia estatística porque a geometria cúbica é invariante apenas sob o grupo octaédrico ($O_h$), não sob o grupo de rotação contínuo $SO(3)$. Consequentemente, o momento angular geralmente não é conservado em caixas periódicas, e a simetria cúbica pode introduzir anisotropias sutis no campo gravitacional, potencialmente enviesando as estatísticas de agrupamento em grande escala, particularmente quando o volume da simulação não é significativamente maior que as escalas de interesse. Além disso, as fronteiras periódicas padrão são inadequadas para simular sistemas com simetrias cilíndricas inerentes, como filamentos cósmicos ou modelos cosmológicos anisotrópicos específicos.

Metodologia
Os autores introduzem um arcabouço de simulação compactificado baseado na variedade topológica $S^1 \times \mathbb{R}^2$ (laje/slab). Nesta geometria, o domínio da simulação é periódico ao longo de um único eixo (o eixo $z$, com período $L_z$) mas aberto na direção radial ($\mathbb{R}^2$).

• Compactação Geométrica: Os autores empregam uma projeção estereográfica inversa ao longo da direção radial para compactificar o plano $\mathbb{R}^2$ não compacto. Isso mapeia a coordenada radial infinita para um domínio finito, permitindo uma estratégia de "zoom-in" onde a resolução espacial e de massa é alta próximo ao eixo central e diminui suavemente em direção à fronteira externa. Isso suprime artefatos de borda na fronteira externa livre.
• Equações de Movimento: As equações de movimento Newtonianas são resolvidas em coordenadas comóveis. O cálculo da força gravitacional leva em conta a topologia $S^1 \times \mathbb{R}^2$ ao somar imagens periódicas ao longo do eixo $z$. Para lidar com a convergência lenta da soma direta, os autores implementam um método de soma de Ewald adaptado para geometria cilíndrica, dividindo a força em componentes de espaço real e espaço recíproco.
• Algoritmos de Cálculo de Força: A implementação dentro do arcabouço StePS (STEreographically Projected cosmological Simulations) suporta dois métodos de cálculo de força:
  1. Direto $O(N^2)$: Mapeia eficientemente para GPUs, adequado para menores quantidades de partículas ou volumes cilíndricos menores.
  2. Octree $O(N \log N)$: Implementado para clusters de CPU. Para mitigar erros sistemáticos decorrentes da natureza cúbica da grade octree em um domínio cilíndrico, o código aplica deslocamentos e rotações aleatórios ao domínio em cada passo de tempo (timestep). Uma etapa crítica de calibração corrige erros de força radial para evitar colapso ou expansão radial espúria, garantindo a estabilidade da densidade de fundo.
• Condições Iniciais: As condições iniciais são geradas usando uma distribuição de partículas do tipo "vidro" (glass) criada via gravidade repulsiva na topologia $S^1 \times \mathbb{R}^2$. Perturbações são impressas usando a aproximação de Zel'dovich ou a teoria de perturbação lagrangiana de segunda ordem (2LPT), adaptadas para a geometria cilíndrica e a amostragem não uniforme da configuração de zoom-in. O espectro de potência é estimado usando o estimador de Feldman–Kaiser–Peacock (FKP) para levar em conta a função de janela não trivial.

Principais Contribuições

1. Implementação de Nova Topologia: O artigo apresenta a primeira implementação de simulações de $N$-corpos cosmológicas na variedade $S^1 \times \mathbb{R}^2$, estendendo as capacidades anteriores de $T^3$ e esféricas $\mathbb{R}^3$ do código StePS.
2. Leis de Conservação: Os autores demonstram que, nesta topologia, o Lagrangiano é invariante sob translação ao longo do eixo do cilindro e rotação no plano transversal ($SO(2)$). Consequentemente, o momento linear ao longo do eixo $z$ e o momento angular ao longo do mesmo eixo são conservados, abordando os problemas de quebra de simetria da topologia $T^3$.
3. Adaptações Algorítmicas: O trabalho detalha as modificações necessárias para o cálculo de força (soma de Ewald cilíndrica), geração de condições iniciais (vidro cilíndrico) e busca de halos (um buscador de sobre-densidade esférica customizado, adaptado para a periodicidade $S^1$ e variação de resolução radial).
4. Validação: Uma simulação $\Lambda$CDM completa é realizada para validar o método, demonstrando a formação de uma teia cósmica consistente com as expectativas cosmológicas padrão.

Resultos

• Execução da Simulação: Uma simulação foi executada com $2,4 \times 10^7$ partículas em um cilindro de altura $L_z = 1,0$ Gpc e raio $R_{sim} = 500$ Mpc, utilizando parâmetros $\Lambda$CDM do Planck 2018.
• Formação de Estrutura: A simulação reproduz com sucesso a teia cósmica não linear (vazios, paredes, filamentos e halos). A estratégia de zoom-in permite a resolução de alta fidelidade de halos de baixa massa próximos ao eixo central, enquanto volumes maiores resolvem halos mais massivos em raios maiores.
• Espectros de Potência: Os espectros de potência de matéria estimados concordam bem com as previsões lineares do CAMB em grandes escalas e com o emulador "halofit" não linear de Mead et al. (2021) em pequenas escalas, confirmando a precisão do método em capturar o crescimento das estruturas.
• Propriedades dos Halos: Os autores identificaram $\sim 1,1 \times 10^5$ halos de matéria escura em $z=0$ e calcularam 68 parâmetros para cada um, incluindo massa, concentração e momento angular.
• Teste de Anisotropia: Motivados por reivindicações recentes de anisotropia no spin das galáxias, os autores testaram a distribuição dos eixos de spin dos halos em relação ao eixo periódico $z$. Usando testes de Kolmogorov-Smirnov, descobriram que para $L_z = 1,0$ Gpc, a distribuição de spin é consistente com a isotropia estatística. Eles concluem que esta topologia, nesta escala, não explica naturalmente as assimetrias relatadas nos spins das galáxias.

Significância e Alegações
O artigo afirma que a topologia $S^1 \times \mathbb{R}^2$ oferece uma alternativa especializada às simulações cúbicas padrão, particularmente para o estudo de sistemas com simetrias cilíndricas (ex: filamentos) ou modelos onde a conservação do momento angular ao longo de um eixo específico é fisicamente relevante. Os autores enfatizam que, embora o método introduza sistemas distintos — como variação radial de resolução e efeitos de seleção — ele fornece um arcabouço autoconsistente livre da anisotropia cúbica das caixas $T^3$.

Os autores são modestos em suas alegações, afirmando que não esperam que este algoritmo substitua as bem testadas simulações cúbicas $T^3$ para cosmologia de propósito geral. Em vez disso, serve como uma ferramenta complementar para questões científicas específicas onde a simetria cúbica é um fator limitante. O trabalho destaca que o campo de força anisotrópico da topologia pode imprimir artefatos se o comprimento periódico for muito pequeno ($L_z \lesssim 200$ Mpc), um tópico reservado para caracterização detalhada futura. O código é disponibilizado como código aberto (StePS v2) para facilitar a exploração adicional destas geometrias.