Controlling artificial surface heating in neutron star simulations: Application to hybrid equations of state

As estrelas de nêutrons são os objetos mais densos do universo, os núcleos colapsados de estrelas massivas que esgotaram seu combustível. Imagine uma esfera do tamanho de uma cidade contendo mais massa do que todo o nosso Sol, compactada tão fortemente que uma única colher de chá de seu material pesaria bilhões de toneladas. Esses laboratórios cósmicos permitem que cientistas estudem a matéria sob condições de pressão e densidade que não podem ser recreadas na Terra. Para entender como essas estrelas se comportam, como se movem e como podem colidir umas com as outras, pesquisadores dependem de simulações computacionais poderosas. Esses modelos digitais atuam como laboratórios virtuais, resolvendo equações complexas para prever a evolução das estrelas e as ondas gravitacionais que elas emitem quando colidem. No entanto, para que essas previsões sejam úteis, os modelos de computador devem ser incrivelmente precisos, especialmente na borda da estrela, onde a matéria densa encontra o vácuo vazio do espaço.

A fronteira entre uma estrela de nêutrons e o vazio circundante é um lugar de extrema dificuldade para simulações computacionais. No mundo real, essa transição é abrupta, mas em um computador, a estrela é representada por uma grade de pontos, como uma malha digital. À medida que a estrela se move ou muda de forma dentro dessa grade, o computador tem dificuldade em lidar com o salto repentino do material denso para o nada. Esse esforço cria um erro sutil, mas persistente: o computador acidentalmente adiciona calor à superfície da estrela onde não deveria existir calor físico. Esse aquecimento artificial não é causado por fricção ou reações nucleares; é uma falha na matemática. Com o tempo, esse calor falso se espalha para o interior, fazendo com que a estrela se expanda ligeiramente e altere sua energia interna. Se não for controlado, esse erro numérico pode distorcer a simulação, levando os cientistas a tirar conclusões incorretas sobre a estrutura da estrela ou os sinais que ela envia ao universo.

Em um estudo recente, um pesquisador da Universidade Goethe de Frankfurt investigou como interromper esse aquecimento artificial. O trabalho focou em um tipo específico de código de computador usado para simular estrelas de nêutrons, conhecido como código BAM. Este código é amplamente utilizado para modelar as colisões violentas de estrelas de nêutrons binárias, onde duas dessas estrelas giram uma em direção à outra e eventualmente se fundem. O pesquisador testou uma nova estratégia matemática chamada limitação de fluxo baseada em entropia. Em termos simples, este método atua como um filtro inteligente para os cálculos do computador. Ele monitora o fluxo de energia e matéria na superfície da estrela e ajusta automaticamente o método de cálculo quando detecta o tipo de mudanças acentuadas que geralmente causam erros. Ao fazer isso, evita que o computador gere esse calor espúrio logo de início.

Para verificar se este método funcionou, o pesquisador executou dois conjuntos de simulações. O primeiro conjunto envolveu estrelas de nêutrons isoladas, sentadas sozinhas no espaço, enquanto o segundo conjunto simulou pares de estrelas girando um em direção ao outro. Em ambos os cenários, o computador teve uma escolha: usar o método de cálculo padrão ou usar o novo filtro baseado em entropia. Os resultados foram claros e consistentes. Quando o método padrão era usado, a energia interna das estrelas aumentava constantemente ao longo do tempo, e as próprias estrelas cresciam, como se estivessem inflamando devido ao calor falso. Essa expansão não era uma realidade física, mas um sintoma do erro numérico. Em contraste, quando o novo filtro era aplicado, as estrelas permaneciam estáveis. Sua energia interna permanecia próxima ao seu valor inicial e seu tamanho não mudava artificialmente. O novo método conseguiu manter as estrelas sem aquecer e expandir, mantendo uma representação muito mais fiel de seu estado físico.

O estudo examinou uma ampla variedade de modelos de estrelas de nêutrons, cada um baseado em diferentes teorias sobre como a matéria se comporta em densidades tão extremas. Alguns modelos assumiam que as estrelas eram feitas de um tipo de partícula, enquanto outros incluíam misturas mais complexas. Independentemente de qual modelo fosse usado, o novo método de filtragem desempenhou igualmente bem. Ele reduziu o aquecimento artificial por uma margem significativa, mantendo o crescimento da energia interna em uma fração minúscula do que era visto com o método antigo. Nas simulações de estrelas binárias, onde a força gravitacional e as forças de maré são intensas, o novo método continuou a suprimir o aquecimento durante toda a longa dança da aproximação, até o momento da colisão. As estrelas simuladas com o novo método permaneceram menores e mais frias do que aquelas simuladas com o método antigo, sugerindo que os artefatos digitais estavam sendo efetivamente removidos.

Este trabalho não afirma ter resolvido todos os problemas de simulação de estrelas de nêutrons, nem sugere que o novo método seja uma solução mágica para todos os desafios numéricos. Em vez disso, fornece evidências fortes de que uma fonte específica de erro — o aquecimento artificial na superfície — pode ser controlada. O pesquisador observa que, embora esses testes tenham sido feitos usando uma forma simplificada de modelar o calor da estrela, o sucesso do método sugere que ele funcionará tão bem quanto ao ser aplicado a modelos mais complexos e realistas que rastreiam a temperatura com maior detalhe. As descobertas oferecem um caminho para a criação de simulações mais precisas de colisões de estrelas de nêutrons. Ao remover essas distorções digitais, os cientistas podem confiar mais profundamente em seus modelos, garantindo que os sinais que eles preveem para os detectores de ondas gravitacionais sejam baseados na verdadeira física do universo, e não nas limitações do código de computador.