Chebyshev interpolation in Einstein-Boltzmann codes

Para compreender o universo, os cosmólogos constroem modelos matemáticos elaborados que descrevem como a matéria e a luz se comportaram desde o Big Bang. Esses modelos baseiam-se num conjunto de equações complexas que acompanham a evolução do cosmos ao longo do tempo, através de diferentes escalas de distância e em todas as direções do céu. Para transformar estas equações em previsões que possam ser comparadas com dados reais de telescópios, os cientistas utilizam programas de computador poderosos conhecidos como códigos Einstein–Boltzmann. Estes programas atuam como laboratórios virtuais, simulando o nascimento da radiação cósmica de fundo — o brilho residual do Big Bang — e a distribuição das galáxias. No entanto, executar estas simulações é incrivelmente difícil porque as equações são rígidas e mudam rapidamente em certos momentos, como quando o universo arrefeceu o suficiente para os átomos se formarem. Para tornar os cálculos geríveis, os programas tradicionalmente resolvem as equações apenas em alguns pontos específicos e depois preenchem as lacunas entre eles utilizando um método chamado interpolação, que é essencialmente desenhar uma curva suave através de pontos de dados conhecidos para adivinhar os valores intermédios. A precisão de toda a simulação depende fortemente de quão bem este jogo de adivinhação funciona.

Um investigador da Universidade de Oslo, Herman Sletmoen, encontrou uma forma de tornar este jogo de adivinhação significativamente mais preciso e rápido. Num novo estudo, Sletmoen demonstra que, ao alterar a técnica matemática utilizada para preencher essas lacunas, o computador pode alcançar uma precisão muito superior com muito menos cálculos. O método tradicional, que utiliza splines cúbicas — uma técnica que liga pontos com pequenas curvas suaves — funciona bem, mas requer um grande número de pontos de dados para atingir altos níveis de precisão. Sletmoen substituiu isto por uma abordagem diferente baseada em polinómios de Chebyshev, uma ferramenta matemática que é excecionalmente boa a ajustar curvas suaves quando os pontos de dados são colocados em intervalos específicos e não uniformes. Ao resolver as difíceis equações da física apenas nestes pontos especiais e depois utilizar o novo método para preencher o resto, o computador pode produzir resultados com erros de interpolação até dez mil vezes menores do que o método antigo para a mesma quantidade de trabalho.

O cerne do problema reside na natureza dos dados que estes códigos geram. O universo evolui suavemente ao longo do tempo e do espaço, exceto em alguns breves momentos de mudança rápida. Quando um computador tenta simular isto, deve calcular o comportamento da luz e da matéria para biliões de diferentes distâncias e ângulos. Fazer isto para cada possibilidade individual é impossível, por isso o código calcula alguns milhares de exemplos e depois interpola o resto. Durante décadas, a ferramenta padrão para este trabalho foram as splines cúbicas. Embora fiáveis, este método é um pouco rudimentar; para obter uma resposta muito precisa, precisa de ser alimentado com uma grelha densa de pontos de dados, o que abranda a simulação. O trabalho de Sletmoen mostra que, como a física subjacente é suave, uma abordagem matemática mais sofisticada pode fazer o trabalho com uma fração dos dados. Ao colocar os pontos de cálculo onde o novo método mais precisa, em vez de os espalhar uniformemente, o computador pode capturar os detalhes finos da evolução do universo de forma muito mais eficiente.

No estudo, Sletmoen testou esta nova abordagem nas duas variáveis mais críticas da simulação: a escala das ondas no universo primitivo e o ângulo através do qual observamos o céu. Quando o código utilizou o novo método para preencher as lacunas para as escalas de onda, o erro na previsão final diminuiu drasticamente. Enquanto o método antigo precisava de cerca de duzentos pontos de dados para atingir um determinado nível de precisão, o novo método atingiu a mesma ou melhor precisão com apenas cinquenta a oitenta pontos. Esta redução nos cálculos necessários traduz-se diretamente em velocidade. Para uma simulação típica da temperatura da radiação cósmica de fundo, o novo método tornou todo o processo de duas vezes e meia a quatro vezes mais rápido, mantendo o mesmo elevado padrão de precisão. Este é um ganho significativo, pois estas simulações são frequentemente executadas milhares de vezes para testar diferentes teorias sobre o universo.

O desafio não era apenas matemático, mas também na forma como o computador lida com os dados. O novo método funciona melhor quando os pontos são colocados em intervalos específicos e desiguais, mas as quantidades físicas sendo medidas, como o ângulo no céu, só existem como números inteiros. Não se pode ter uma medição a meio grau da mesma forma que se pode ter uma medição a meio metro. Para resolver isto, Sletmoen desenvolveu um contorno inteligente. Em vez de forçar o computador a calcular valores em ângulos não inteiros impossíveis, o método arredonda os pontos de cálculo especiais para o número inteiro mais próximo. Este pequeno ajuste permite que a poderosa nova matemática seja utilizada sem quebrar as regras da simulação física. O estudo descobriu que esta abordagem de "arredondamento" funcionou tão bem quanto a versão matemática ideal, provando que a técnica é robusta o suficiente para a cosmologia do mundo real.

Os resultados sugerem uma mudança na forma como estas simulações complexas são construídas. Os códigos de computador mais antigos utilizavam frequentemente atalhos e aproximações para acelerar as coisas, alternando entre diferentes equações dependendo da situação. Embora rápidos, estes atalhos podem introduzir pequenos relevos ou irregularidades nos dados que fazem com os métodos de interpolação tradicionais terem dificuldades. A nova abordagem brilha mais quando utilizada com uma nova geração de códigos que evitam estes atalhos e resolvem as equações continuamente, mantendo os dados perfeitamente suaves. Ao combinar estes códigos suaves e livres de aproximações com o novo método de interpolação, os cientistas podem obter o melhor de dois mundos: a velocidade de menos cálculos e a precisão de um modelo altamente exato. Esta combinação permite aos investigadores explorar o universo com uma clareza que era anteriormente demasiado dispendiosa para computar, abrindo a porta a testes mais detalhados de teorias cosmológicas sem a necessidade de esperar por computadores mais rápidos.