Impact of neutron captures on the r-process: I. Propagation of statistical model uncertainties

O universo é uma vasta fábrica química, mas durante décadas, os cientistas lutaram para explicar como aproximadamente metade de todos os elementos mais pesados que o ferro foram forjados. Esses elementos pesados, que incluem desde o ouro até o urânio, são criados nos eventos cósmicos mais violentos imagináveis, como a colisão de duas estrelas de nêutrons. Nesses ambientes extremos, núcleos atômicos são bombardeados por uma enxurrada de nêutrons, construindo rapidamente átomos cada vez mais pesados em um processo conhecido como processo de captura de nêutrons rápida. Embora os astrônomos tenham confirmado que essas colisões ocorrem e produzem os tipos certos de elementos, a física precisa que governa como os átomos mudam durante a explosão permanece em grande parte um mistério. Como não podemos recriar essas condições em um laboratório, os pesquisadores devem confiar em modelos computacionais complexos para simular a jornada da matéria, de um átomo simples a um elemento pesado. A precisão dessas simulações depende inteiramente dos dados inseridos nelas, especificamente de como entendemos os estados de energia interna de átomos instáveis e ricos em nêutrons que existem apenas por um momento fugaz.

Uma equipe de pesquisadores liderada por Yu An, da Universidade Técnica de Darmstadt, mergulhou profundamente em uma das peças mais críticas deste quebra-cabeça: os modelos estatísticos usados para prever como esses átomos instáveis capturam nêutrons. Em seu estudo, eles não tentaram inventar uma nova física, mas sim testar o quanto os resultados finais de uma explosão cósmica mudam quando as suposições subjacentes sobre a estrutura atômica são ligeiramente alteradas. Eles se concentraram em dois ingredientes teóricos específicos: a densidade de níveis nucleares, que descreve de quantas maneiras diferentes um núcleo atômico pode vibrar ou armazenar energia, e a função de força de raios gama, que mede a probabilidade de o núcleo emitir energia luminosa ao se estabilizar. Ao executar milhares de simulações com diferentes combinações desses modelos teóricos, a equipe mapeou exatamente onde nossa falta de conhecimento cria as maiores lacunas em nossa compreensão da história química do universo.

Os pesquisadores começaram simulando o rescaldo de uma fusão de estrelas de nêutrons, um cenário onde duas estrelas mortas colidem e ejetam uma nuvem de material superdenso. Eles executaram essa simulação usando sessenta combinações diferentes dos modelos teóricos disponíveis, criando uma ampla gama de possíveis resultados para a mistura final de elementos. Os resultados revelaram uma hierarquia clara de incerteza. Os modelos que descrevem como os núcleos armazenam energia, conhecidos como densidades de níveis nucleares, tiveram um impacto muito maior na composição química final do que os modelos que descrevem como eles emitem luz. Enquanto os grandes picos de elementos pesados, como os correspondentes ao ouro e à platina, permaneceram relativamente estáveis independentemente do modelo usado, as regiões entre esses picos mostraram diferenças dramáticas. Em algumas simulações, a abundância de elementos em uma faixa de massa específica em torno de 172 a 178 unidades atômicas triplicaria ou desapareceria inteiramente dependendo de qual modelo de energia fosse escolhido.

Para entender por que isso acontece, a equipe observou de perto o tempo dos eventos dentro da nuvem de detritos em expansão. Após a explosão inicial, a temperatura cai rapidamente. No início, os átomos estão tão quentes que os nêutrons são constantemente capturados e depois expelidos em uma dança equilibrada. Eventualmente, a temperatura cai o suficiente para que esse equilíbrio se quebre, e os átomos comecem a se estabelecer em suas formas finais. Os pesquisadores descobriram que a escolha do modelo de energia torna-se crítica na fase "fria" deste processo, que ocorre quando a temperatura cai para cerca de 0,45 bilhão de graus. Nesta janela, se o modelo teórico prevê que os átomos capturam nêutrons lentamente, o material permanece em uma região específica da tabela periódica, criando um subpico distinto de elementos. Se o modelo prevê uma captura mais rápida, esse material passa rapidamente para se tornar elementos ainda mais pesados, deixando um vazio para trás. Essa sensibilidade não é apenas uma peculiaridade de uma única simulação; a equipe descobriu que mesmo ao calcular a média dos resultados de mais de 11.000 diferentes cenários de explosão para imitar uma mistura galáctica realista, esses gaps e picos específicos persistiam.

O estudo também explorou uma vasta grade de condições hipotéticas, variando a densidade inicial e a razão entre prótons e elétrons no material para ver como diferentes ambientes reagem a essas incertezas. Eles descobriram que o impacto dos modelos de energia não é uniforme em todos os eventos cósmicos. Em ambientes onde o suprimento de nêutrons é escasso, a escolha do modelo importa muito pouco porque os átomos simplesmente ficam sem combustível antes que possam ser significativamente alterados. No entanto, em ambientes com um suprimento rico de nêutrons e condições termodinâmicas específicas, a diferença entre os modelos pode determinar se um determinado elemento é produzido em abundância ou quase não é produzido. Os pesquisadores identificaram que as incertezas mais significativas surgem em cenários "frios", onde os átomos continuam a capturar nêutrons muito depois que o equilíbrio inicial foi quebrado. Nessas condições, pequenas mudanças na velocidade prevista de captura de nêutrons podem deslocar todo o fluxo de matéria, decidindo se o processo para no segundo pico de elementos pesados ou se avança para o terceiro.

Em última análise, este trabalho serve como um mapa crucial para experimentos futuros. Ele diz aos físicos nucleares exatamente quais regiões da tabela atômica são mais sensíveis às nossas lacunas teóricas atuais. O estudo confirma que, embora tenhamos uma imagem geral de como os elementos pesados são feitos, os detalhes finos da paisagem química dependem fortemente de quão bem entendemos os estados de energia interna de núcleos instáveis. Os pesquisadores concluem que, para refinar nossa compreensão da evolução química do universo, precisamos de melhores dados experimentais e descrições teóricas mais precisas de densidades de níveis nucleares, particularmente para os isótopos específicos que se situam entre os principais picos de elementos pesados. Sem essas melhorias, nossas simulações do cosmos continuarão a carregar uma sombra de incerteza, incapazes de explicar plenamente a abundância precisa dos elementos pesados que compõem o nosso mundo.