DSWIM:Efficient and Stable Deterministic Computation of Warm Inflation Perturbations

Resumo Técnico: DSWIM – Computação Determinística Eficiente e Estável de Perturbações de Inflação Quente

Problema
As perturbações da inflação quente (WI - Warm Inflation) são impulsionadas tanto por flutuações quânticas quanto térmicas, exigindo a solução de equações diferenciais estocásticas acopladas. Embora o código SWIM, disponível publicamente, implemente uma abordagem estocástica robusta, ele requer a média de um conjunto (ensemble averaging) de muitas realizações para atingir a convergência estatística, o que o torna computacionalmente caro para inferência de parâmetros e avaliações de verossimilhança (likelihood).

Uma alternativa determinística baseada na abordagem de Fokker-Planck, que evolui a matriz de correlação das variáveis de perturbação via equações diferenciais ordinárias (ODEs), oferece uma redução significativa no custo computacional. No entanto, os autores identificam que a evolução determinística direta sofre de um severo mau condicionamento numérico em muitos modelos de WI. Essa instabilidade surge porque as variáveis de perturbação evoluem em escalas amplamente distintas, fazendo com que os elementos da matriz de correlação abrancem muitos ordens de magnitude (ex: 13 a 90 ordens nos modelos testados). Essa hierarquia leva a erros de arredondamento de ponto flutuante, distorções artificiais no fator de correção $G(Q)$ e o eventual colapso da evolução numérica, particularmente no regime dissipativo forte. Além disso, implementações determinísticas anteriores (ex: WI2easy) exibiram discrepâncias sistemáticas com resultados estocásticos no regime dissipativo moderado devido ao tratamento incorreto do ruído térmico correlacionado.

Metodologia
Para abordar essas questões, os autores introduzem o DSWIM (Deterministic SWIM), um novo módulo dentro do framework SWIM que implementa um framework de matriz de escala para as equações de evolução determinística.

1. Transformação de Escala: Os autores definem um vetor de estado escalonado $\tilde{\Phi} = S\Phi$, onde $S$ é uma matriz de escala diagonal composta por potências do parâmetro de Hubble $H$. A escala é motivada fisamente pelo comportamento de escala efetivo das variáveis de perturbação durante a inflação:

   • A perturbação métrica $\psi$ (adimensional) permanece não escalonada.
   • As flutuações do inflaton $\delta\phi$ e suas derivadas são escalonadas por $H^{-1}$.
   • A perturbação da densidade de energia da radiação $\delta\rho_r$ é escalonada por $H^{-2}$.
   • A perturbação do momento da radiação $\delta q_r$ é escalonada por $H^{-1}$ com base em sua aparição nas equações de perturbação.
     Esta transformação normaliza as variáveis, garantindo que os elementos da matriz de correlação evoluída permaneçam numericamente comparáveis ao longo da integração, preservando assim a faixa dinâmica e evitando a perda de precisão.

2. Preservação de Observáveis Físicos: Os autores demonstram que esta transformação de escala deixa a perturbação de curvatura comóvel $R$ e, consequentemente, o espectro de potência escalar primordial $P_R$ exatamente invariantes. A matriz de deriva $\tilde{A}$ e a matriz de difusão $\tilde{D}$ transformadas são derivadas para manter a consistência física do sistema.

3. Correção de Correlações de Ruído: O formalismo determinístico é atualizado para incorporar explicitamente as correlações fora da diagonal na matriz de difusão $D$. Os autores mostram que os termos de ruído térmico nas equações do inflaton e da radiação não são independentes, mas sim correlacionados devido à conservação de energia-momento. A incorporação desses termos fora da diagonal não nulos resolve as discrepâncias sistemáticas observadas anteriormente entre códigos estocásticos e determinísticos no regime dissipativo moderado.

Resultos Principais

• Estabilidade Numérica: O framework determinístico escalonado suprime com sucesso os artefatos numéricos observados em implementações não escalonadas. Para um modelo de potencial exponencial de fuga (runaway exponential potential), a transformação de escala reduz a hierarquia da matriz de covariância em aproximadamente 40–50 ordens de magnitude, estendendo o alcance de evolução estável do solver determinístico por quase uma ordem de magnitude na razão de dissipação $Q$ em comparação com métodos não escalonados.
• Concordância com Métodos Estocásticos: Os resultados do DSWIM mostram excelente concordância com a implementação estocástica do SWIM através dos regimes dissipativos fraco, intermediário e forte, validando a precisão física da abordagem determinística escalonada.
• Eficiência Computacional: A abordagem determinística produz acelerações substanciais. Para modelos representativos, o DSWIM reduz o tempo de execução para o cálculo do fator de correção $G(Q)$ em aproximadamente duas ordens de magnitude em comparação com o SWIM estocástico e supera significativamente o código determinístico existente, WI2easy. Por exemplo, o tempo de execução para um modelo de potencial cúrtico caiu de ~116 segundos (estocástico) para ~1,56 segundos (determinístico).
• Resolução de Discrepâncias: Ao implementar corretamente a estrutura de ruído térmico correlacionado na matriz de difusão, o DSWIM elimina as desviações sistemáticas entre resultados estocásticos e determinísticos no regime dissipativo moderado ($0,01 \lesssim Q \lesssim 10$).

Significância e Alegações
O artigo estabelece o DSWIM como um framework numericamente robusto e computacionalmente eficiente para calcular perturbações escalares de inflação quente. Os autores alegam que o DSWIM é atualmente o único solver de perturbações de WI disponível publicamente que fornece um formalismo determinístico numericamente estabilizado diretamente interfaceado com o framework de análise de verossimilhança Cobaya. Esta capacidade permite estudos eficientes de inferência de parâmetros de modelos de inflação quente que eram anteriormente proibitivos do ponto de vista computacional devido ao custo da média estocástica ou da instabilidade numérica de métodos determinísticos não escalonados. O trabalho não propõe novos modelos físicos, mas sim fornece uma ferramenta computacional necessária para explorar com precisão e eficiência o espaço de parâmetros de cenários de inflação quente existentes.