STEPSIC: Initial condition generator for stereographic cosmological simulations

Resumo Técnico do stepsic: Gerador de condições iniciais para simulações cosmológicas estereográficas

Declaração do Problema
Os geradores convencionais de condições iniciais (CI) cosmológicas são projetados exclusivamente para domínios cúbicos totalmente periódicos (topologia $T^3$). Embora bem-sucedidos para estudos padrão de formação de estruturas, a topologia $T^3$ introduz correlações artificiais em escalas comparáveis ao tamanho da caixa e falha em conservar o momento angular devido à falta de invariância sob o grupo completo de rotações $SO(3)$. Essas limitações tornam as ferramentas padrão inadequadas para simulações não periódicas centradas no observador ou para fenômenos que exigem invariância rotacional. Especificamente, o código StePS (STEreographically Projected cosmological Simulations), que utiliza projeção estereográfica inversa para mapear a variedade $R^3$ em domínios esféricos ($R^3$) e cilíndricos ($S^1 \times R^2$) finitos com resolução de partículas variando radialmente, carece de um gerador de CI público compatível. Códigos existentes (por exemplo, 2lptic, music, monofonic) não conseguem produzir as geometrias não padrão necessárias ou lidar com os gradientes intrínsecos de resolução de massa das projeções estereográficas.

Metodologia
O artigo apresenta o stepsic, um gerador de condições iniciais de código aberto projetado para preencher essa lacuna. O código estende a Teoria de Perturbação Lagrangiana (LPT) para suportar domínios esféricos, cilíndricos e de paralelepípedos com razão de aspecto arbitrária. O fluxo de trabalho técnico envolve:

1. Geração de Carga de Partículas: O stepsic suporta quatro estratégias de carregamento:

   • Rede Regular: Vóxeis cúbicos para domínios periódicos $T^3$.
   • Distribuição de Poisson: Amostragem aleatória homogênea para testes ou geração de vidro.
   • Cascas/Anéis Concêntricos: Específico para geometrias StePS, preenchendo o volume com cascas (esféricas) ou anéis (cilíndricos) contendo contagens fixas de partículas. Isso cria a resolução de massa variando radialmente intrínseca às projeções estereográficas. Dois modos de binagem são suportados: ângulo estereográfico constante ($\Delta\omega$) e volume de espaço compacto constante. Uma zona interna de resolução constante é definida para mitigar a degradação da precisão da força causada pela interação de partículas com massas drasticamente diferentes.
   • Entradas de Vidro: Capacidade de ler e reescalar configurações de vidro geradas externamente.

2. Construção de Campos:

   • Grade de Fourier: Grades de Fourier livres de anisotropia com vóxeis cúbicos são construídas. Para domínios não cúbicos, a grade é condicionada pela geometria da caixa, o que modifica o conjunto de modos acessíveis (suprimindo a variância ao longo dos eixos curtos) em vez de simplesmente recortar uma realização isotrópica maior.
   • Campos de Densidade: Campos aleatórios gaussianos são gerados usando ruído branco, reescalados para o espectro de potência da matéria linear alvo $P(k, z)$ no desvio para o vermelho inicial. O código suporta realizações padrão, realizações emparelhadas (com deslocamento de fase) e realizações de amplitude fixa para reduzir a variância cósmica.
   • Implementação LPT: O código implementa a teoria de perturbação lagrangiana de primeira ordem (aproximação de Zel'dovich) e de segunda ordem (2LPT). Campos de deslocamento e velocidade são computados no espaço de Fourier e transformados para o espaço físico. A implementação segue convenções estabelecidas pelo 2lptic e monofonic, garantindo compatibilidade com fluxos de trabalho padrão enquanto se estende a topologias não padrão.

3. Interpolação e Multirresolução:

   • Os campos de deslocamento e velocidade são interpolados nas posições das partículas usando kernels de atribuição de massa B-spline (NGP, CIC, TSC) com deconvolução no espaço de Fourier para corrigir o alisamento.
   • Um esquema de multirresolução mapeia os campos de deslocamento através da resolução de massa variando radialmente. As partículas recebem campos de deslocamento de grades de resolução variável correspondentes à sua massa, garantindo que partículas de alta massa (representando volumes maiores) não recebam flutuações de pequena escala não resolvidas que atuariam como ruído.

Principais Resultados e Validação
Os autores validam o stepsic através de três níveis de teste: reprodutibilidade, correção numérica e concordância estatística com a teoria e outros códigos.

• Recuperação do Espectro de Potência: Em caixas cúbicas periódicas, o espectro de potência da matéria recuperado concorda com a previsão da teoria linear de entrada em melhor que 0,5% até metade do número de onda de Nyquist. Essa precisão mantém-se independente da razão de aspecto da caixa (testada até 10:1).
• Validação Cruzada: Quando comparado ao código estabelecido monofonic usando campos de ruído branco e parâmetros cosmológicos idênticos, o stepsic alcança concordância de espectro de potência subpercentual. Um pequeno deslocamento residual ($\approx 0,6\%$) é atribuído a diferenças na implementação dos cálculos do fator de crescimento linear (usando a biblioteca colossus no stepsic versus um solucionador interno no monofonic).
• Supressão de Transientes: A evolução completa N-corpos de execuções StePS cilíndricas correspondentes confirma que a LPT de segunda ordem suprime corretamente o excesso de potência transitória de 2–3% presente nas condições iniciais de primeira ordem em escalas quase lineares.
• Consistência Geométrica: Testes em geometrias de placas demonstram que a construção de grade de Fourier com vóxeis cúbicos não introduz viés sistemático dependente da geometria para razões de aspecto até 10:1, além da redução esperada na variância amostral para volumes menores.

Significado e Alegações
O artigo posiciona o stepsic como o primeiro gerador de condições iniciais público capaz de suportar as geometrias não padrão ($R^3$ e $S^1 \times R^2$) e as resoluções de massa variando radialmente exigidas pelo código StePS. Ao estender com sucesso a LPT a esses domínios e validar a saída contra expectativas teóricas e códigos existentes, os autores demonstram que simulações cosmológicas precisas e centradas no observador podem ser realizadas sem as correlações artificiais inerentes às topologias periódicas $T^3$. O código é apresentado como uma ferramenta necessária para estudar estruturas filamentares e modelos cosmológicos anisotrópicos dentro do framework StePS. Os autores observam modestamente que, embora a LPT de terceira ordem e a não gaussianidade primordial estejam planejadas para futuras versões, a implementação atual foca em fornecer uma base robusta de segunda ordem para simulações estereográficas.