The Mass Function of Neutron Stars from Core-Collapse Supernova Simulations

No coração profundo de estrelas massivas, um drama final e violento se desenrola. Quando uma estrela muito mais pesada que o nosso Sol esgota seu combustível, seu núcleo colapsa sob seu próprio peso, desencadeando uma explosão de supernova que pode ofuscar uma galáxia inteira. No rescaldo, o núcleo frequentemente sobrevive como uma estrela de nêutrons, um objeto tão denso que uma única colher de chá de seu material pesaria tanto quanto uma montanha. Por décadas, os astrônomos tentaram compreender o nascimento desses remanescentes cósmicos. Eles mediram as massas de estrelas de nêutrons encontradas em sistemas binários, onde duas estrelas orbitam uma à outra, mas as razões por trás da distribuição específica desses pesos permaneceram elusivas. Existe um peso padrão para uma estrela de nêutrons recém-nascida, ou elas nascem com uma ampla variedade de massas? E como o processo caótico da morte de uma estrela determina se ela se tornará uma estrela de nêutrons ou colapsará completamente em um buraco negro?

Um novo estudo realizado por uma equipe de pesquisadores reúne as simulações de computador mais avançadas dessas explosões estelares com as observações mais recentes de estrelas de nêutrons para responder a essas perguntas. Ao mapear a estrutura interna da estrela moribunda para o peso do remanescente que ela deixa para trás, a equipe construiu um quadro teórico de como as estrelas de nêutrons nascem. O trabalho deles sugere que a variedade de massas que vemos hoje não é apenas o resultado de estrelas ganhando peso mais tarde em suas vidas, mas está escrita no exato momento de sua criação. As descobertas oferecem uma explicação coerente para o motivo pelo qual a maioria das estrelas de nêutrons se agrupa em torno de um peso específico, enquanto um grupo menor se estende para massas muito mais pesadas, e fornecem uma nova estimativa de quão frequentemente estrelas massivas falham em se tornar estrelas de nêutrons, colapsando, em vez disso, em buracos negros.

Os pesquisadores começaram analisando a intrincada relação entre o tamanho inicial de uma estrela massiva e a estrutura de seu núcleo pouco antes de explodir. Eles utilizaram modelos de computador tridimensionais sofisticados que simulam a física complexa de uma supernova, rastreando como o núcleo da estrela se comporta nos momentos finais de sua vida. Essas simulações revelaram um padrão claro: a "compactação" interna do núcleo da estrela — uma medida de quão densamente compacta é a matéria em um ponto específico em seu interior — determina diretamente a massa da estrela de nêutrons que se forma. Um núcleo menos compacto tende a deixar para trás uma estrela de nêutrons mais leve, enquanto um núcleo mais compacto resulta em uma mais pesada. Essa conexão permitiu que a equipe pegasse a distribuição conhecida de estrelas massivas no universo e previsse a distribuição das massas das estrelas de nêutrons no momento de seu nascimento.

Quando compararam essa distribuição de nascimento teórica com as massas reais de estrelas de nêutrons observadas no céu, surgiu uma discrepância. As simulações previam que as estrelas de nêutrons recém-nascidas mais comuns deveriam ser ligeiramente mais leves que a massa média observada, com um grupo secundário de estrelas mais pesadas. No entanto, os dados observados mostraram um pico em um peso ligeiramente superior. Os pesquisadores perceberam que essa diferença poderia ser explicada por dois fatores principais: as estrelas de nêutrons ganhando massa após o nascimento e a maneira como o nascimento delas afeta se elas permanecem em um sistema binário.

Estrelas de nêutrons em sistemas binários frequentemente puxam material de suas estrelas companheiras, um processo chamado acreção, que adiciona peso à estrela de nêutrons ao longo do tempo. A equipe aplicou correções aos dados observados para estimar quais teriam sido as massas no momento do nascimento, removendo o peso ganho de seus companheiros. Esse ajuste deslocou o pico observado para baixo, alinhando-o muito mais proximamente com a previsão teórica. Além disso, o estudo destacou o papel dos "chutes" dados às estrelas de nêutrons no nascimento. Assim como um foguete é impulsionado para frente pelo seu escape, uma estrela de nêutrons recebe um empurrão repentino na direção oposta à da explosão. As simulações sugerem que estrelas de nêutrons mais pesadas recebem chutes mais fortes, em média. Esses chutes poderosos podem quebrar a ligação gravitacional que mantém o sistema binário unido, ejetando as estrelas de nêutrons mais pesadas para o vazio, onde são mais difíceis de encontrar. Essa remoção seletiva significa que os sistemas binários que observamos são enviesados para estrelas de nêutrons mais leves, que receberam chutes mais suaves e permaneceram em órbita.

Ao combinar os efeitos da formação de buracos negros, da acreção de massa e desses chutes de nascimento, a equipe criou um modelo refinado que combina com os dados observados de forma notável. Seu modelo sugere que o universo produz um amplo espectro de massas de estrelas de nêutrons no nascimento, variando de cerca de 1,2 a 2,0 vezes a massa do nosso Sol. A massa de nascimento mais comum é em torno de 1,35 vezes a massa do Sol, o que explica o pico proeminente visto nas observações. A cauda mais pesada da distribuição, estendendo-se até 2,0 massas solares, é um resultado natural das estruturas variáveis das estrelas progenitoras, e não apenas o resultado de acreção posterior.

O estudo também forneceu uma nova estimativa para a fração de estrelas massivas que terminam suas vidas como buracos negros em vez de estrelas de nêutrons. Ao analisar quais tipos de estrelas em suas simulações falharam em produzir uma estrela de nêutrons, os pesquisadores calcularam que aproximadamente 21 por cento das estrelas massivas na faixa de massa relevante colapsam diretamente em buracos negros. Esse número não é uma regra simples baseada no tamanho inicial da estrela; em vez disso, depende da estrutura complexa e não linear do núcleo da estrela, o que significa que estrelas de tamanhos semelhantes podem ter destinos muito diferentes.

Este trabalho representa um passo significativo para conectar a teoria de como as estrelas morrem com a realidade do que observamos. Ele vai além da antiga ideia de que todas as estrelas de nêutrons nascem com quase a mesma massa, substituindo-a por uma imagem dinâmica onde a estrutura inicial da estrela dita uma ampla gama de resultados. O acordo entre as simulações complexas e os dados observados sugere que nossa compreensão das explosões de supernova está amadurecendo. Embora ainda existam incertezas, como a influência precisa da composição química da estrela ou os detalhes das interações binárias, o estudo demonstra que a teoria moderna de supernovas pode agora explicar a população de objetos compactos que vemos no céu. Ele abre uma nova era onde o nascimento desses objetos extremos pode ser rastreado de volta às condições específicas das estrelas que os criaram, transformando uma coleção de medições em uma história coerente de morte e renascimento estelar.