Angular and Kinetic Properties of Scission Neutrons within Time-dependent Density Functional Theory

Este estudo utiliza a teoria do funcional da densidade dependente do tempo para demonstrar que os nêutrons de cisão constituem um componente significativo de alta energia do espectro de nêutrons de fissão prompt em fissões de $^{235}\mathrm{U}$, $^{239}\mathrm{Pu}$ e $^{252}\mathrm{Cf}$, fornecendo evidência direta de que sua inclusão resolve a subestimação sistemática dos rendimentos de alta energia observada em modelos tradicionais baseados apenas em evaporação.

O essencial

Imagine um balão gigante e instável (um núcleo atômico) que de repente se parte em dois. Isso é a fissão nuclear. Por décadas, os cientistas sabem que, quando isso acontece, os dois pedaços resultantes (chamados de fragmentos) voam para longe a velocidades incríveis e cospem minúsculas partículas chamadas nêutrons.

Por muito tempo, os físicos pensaram que todos esses nêutrons eram "evaporados" mais tarde, como o vapor subindo de uma xícara de café quente depois que a água ferveu. Eles assumiam que os fragmentos já estavam totalmente formados e se movendo de forma constante antes de começarem a soltar esses nêutrons.

Este novo artigo sugere que alguns nêutrons são na verdade "chutados para fora" exatamente no momento em que o balão se parte. Estes são chamados de nêutrons de cisão. Eles nascem no caos de um segundo — o caos da divisão — e não de um fragmento calmo e resfriando mais tarde.

Aqui está como os pesquisadores encontraram a prova desses nêutrons do "momento do estalo", explicada de forma simples:

1. A Simulação de Supercomputador

Para ver o que acontece durante a divisão, os cientistas não usaram um microscópio; eles usaram um supercomputador para rodar um filme do evento usando uma teoria chamada Teoria do Funcional de Densidade Dependente do Tempo (TDDFT).

Pense nisso como um videogame 3D de alta velocidade onde eles simulam os átomos dançando e se partindo. Em versões anteriores deste "jogo", o mundo virtual era pequeno demais. Os nêutrons batiam na borda da tela antes que os cientistas pudessem descobrir exatamente quão rápido estavam indo ou para qual direção estavam voando.

Neste estudo, eles construíram um mundo virtual muito maior (cerca de 3 vezes maior que antes). Isso deu aos nêutrons espaço suficiente para voar e se estabilizar, permitindo que os cientistas os medissem com precisão sem que as "paredes" da simulação atrapalhassem os dados.

2. A Descoberta do "Limite de Velocidade"

Uma vez que tiveram uma visão clara, eles observaram os nêutrons saindo em ângulos específicos (principalmente para os lados e levemente para trás em relação à divisão). Eles descobriram algo surpreendente:

• A "Zona Proibida": Não existem nenhum nêutrons de cisão com baixa energia (abaixo de cerca de 1,5 a 2 milhões de elétron-volts). É como se houvesse um limite de velocidade; nada lento é permitido ser um nêutron de "momento do estalo".
• A Multidão de Alta Velocidade: Em vez disso, todos esses nêutrons são rápidos. Eles se agrupam em torno de uma velocidade específica e alta (3–3,5 MeV) e depois se estendem em uma longa cauda de partículas ainda mais rápidas.

É como uma multidão de pessoas pulando de um trampolim. Os nêutrons "evaporados" são como pessoas caminhando casualmente para fora da borda da piscina mais tarde. Os nêutrons de "cisão" são como pessoas sendo violentamente arremessadas para fora do trampolim no exato segundo em que ele quebra. Aqueles arremessados do trampolim são sempre rápidos; você nunca vê um lento vindo desse evento específico.

3. Resolvendo o Mistério da Energia "Faltante"

Cientistas tentam há anos combinar seus modelos de computador com experimentos do mundo real. Eles tinham um problema:

• O Modelo Antigo: Se você contar apenas o "vapor" (nêutrons evaporados), seu modelo de computador prevê poucos nêutrons de alta energia. É como tentar encher um balde com uma caneca pequena, mas o balde continua precisando de mais água do que a caneca pode fornecer.
• O Novo Modelo: Quando os pesquisadores adicionaram os nêutrons do "momento do estalo" (aqueles que encontraram em sua grande simulação) aos nêutrons de "vapor", a matemática finalmente funcionou. O modelo combinado combinou perfeitamente com os dados de alta energia medidos em experimentos reais para Urânio e Califórnio.

4. Por Que Isso Importa

Isso é um grande feito porque é a primeira vez que uma teoria puramente microscópica (uma que não apenas supõe ou assume que coisas existem) previu esses nêutrons de "momento do estalo" e provou que eles são reais.

• Antes: Os cientistas tinham que supor que esses nêutrons existiam porque a matemática não batia.
• Agora: A simulação de computador produziu naturalmente esses nêutrons sem que lhe dissessem para fazê-lo. É como prever uma tempestade observando o movimento das nuvens, em vez de apenas assumir que uma tempestade está chegando porque o relatório meteorológico diz que sim.

A Conclusão

O artigo conclui que, quando um átomo se divide, uma pequena mas importante parte dos nêutrons (cerca de 6% a 10% do total) nasce no momento violento da quebra. Esses nêutrons são distintos porque são sempre rápidos e nunca lentos dentro de certos ângulos.

Ao encontrar essa "impressão digital" nos dados, os pesquisadores finalmente separaram os nêutrons do "momento do estalo" dos nêutrons de "vapor", dando-nos uma imagem mais clara e precisa de como a fissão nuclear realmente funciona. Isso ajuda a refinar nossa compreensão das forças fundamentais que mantêm a matéria unida e a despedaça.

Resumo técnico

Resumo Técnico: Propriedades Angulares e Cinéticas de Nêutrons de Escisão dentro da Teoria do Funcional de Densidade Dependente do Tempo

Definição do Problema
Apesar da descoberta da fissão nuclear há quase um século, uma descrição microscópica completa do processo permanece elusiva, particularmente em relação à transição do ponto de sela externo para a escisão e a subsequente formação de fragmentos de fissão (FFs) primários. Códigos estatísticos de evaporação atuais assumem que todos os nêutrons prontos são emitidos isotropicamente a partir de fragmentos totalmente acelerados. No entanto, esse arcabouço não consegue explicar nêutrons emitidos durante o processo dinâmico altamente não-equilibrado de ruptura do colo, conhecido como nêutrons de escisão (SNs). Embora os SNs tenham sido conjecturados desde 1939 e estudados via diversos modelos, eles permanecem experimentalmente não resolvidos porque observáveis de nêutrons prontos sozinhos não conseguem distinguir entre nêutrons emitidos na escisão e aqueles evaporados posteriormente. Tentativas anteriores de prever SNs usando a teoria do funcional de densidade dependente do tempo (TDDFT) foram limitadas por domínios de simulação pequenos, impedindo a determinação simultânea das propriedades angulares e cinéticas antes que os nêutrons atingissem as fronteiras computacionais.

Metodologia
Este estudo emprega um arcabouço totalmente microscópico baseado em TDDFT estendida para sistemas superfluidos para investigar a emissão de SN em $^{235}\text{U}(n_{th}, f)$, $^{239}\text{Pu}(n_{th}, f)$ e $^{252}\text{Cf}(sf)$.

• Configuração da Simulação: Os autores utilizaram o solver de oscilador harmônico axialmente simétrico AxialHOHFB com o funcional de densidade de energia (EDF) SkM*. Um avanço crítico foi o uso de um domínio de simulação substancialmente maior (um elipsoide com semieixos $Z_{max} = 70$ fm e $R_{max} = 55$ fm), aproximadamente 3,2 vezes maior que em estudos anteriores.
• Integração em Duas Etapas: Para gerenciar restrições computacionais mantendo a precisão, um procedimento de integração temporal de duas etapas foi empregado. O sistema foi evoluído da solução estática de Hartree-Fock-Bogoliubov (HFB) com restrição para o momento de escisão em uma base menor ($B_1$), sendo então transformado para a base maior ($B_2$) para continuar a evolução até que os nêutrons se aproximassem da fronteira.
• Extração de Propriedades de SN: Os SNs foram identificados definindo uma região $V^{(SN)}_\eta$ onde a razão entre a densidade de energia de fluxo coletivo e a densidade de energia cinética intrínseca excedia um limiar $\eta = 3$. Este critério isola nêutrons dominados pelo fluxo coletivo externo.
• Análise Angular e Cinética: Dentro desta região, campos de velocidade local e energia de fluxo coletivo por nêutron foram calculados. O estudo focou no intervalo angular $110^\circ \leq \theta \leq 150^\circ$ (medido em relação à direção do fragmento leve) onde o rendimento de SN satura antes das reflexões de fronteira ocorrerem.
• Comparação com Experimentos: O espectro de energia cinética de SN calculado foi combinado com um modelo de evaporação Maxwelliana para o componente de nêutron evaporado (EN). Os parâmetros do modelo de evaporação (temperatura do fragmento $T$ e razão de multiplicidade $\bar{\nu}_L/\bar{\nu}_H$) foram restringidos exclusivamente por dados experimentais de baixa energia ($E \leq E_{thr}$). Este modelo combinado foi então comparado contra dados de espectro de nêutrons de fissão pronta (PFNS) de alta energia para $^{239}\text{Pu}$ e $^{252}\text{Cf}$.

Principais Resultados

• Características do Espectro de SN: Os espectros de SN extraídos para os três sistemas exibem um comportamento de limiar distinto. Os SNs estão ausentes abaixo de um limiar de energia específico ($E_{thr}$), que é aproximadamente 1,5–2 MeV (especificamente $\approx 2$ MeV para fissão induzida e $\approx 1,5$ MeV para fissão espontânea). Os espectros apresentam um pico próximo a 3–3,5 MeV e possuem uma cauda de alta energia exponencial que se estende até cerca de 18 MeV.
• Estimativas de Rendimento: Para $^{239}\text{Pu}(n_{th}, f)$, o rendimento total de SN é calculado como $0,19 \pm 0,01$ (aproximadamente 6,6% do total de nêutrons prontos). Para $^{252}\text{Cf}(sf)$, o rendimento é $0,39 \pm 0,02$ (aproximadamente 10,4%).
• Validação do Modelo: Quando o componente de SN calculado é adicionado ao modelo de evaporação (restringido apenas por dados de baixa energia), o PFNS total resultante reproduz com precisão o rendimento de alta energia medido para $^{239}\text{Pu}$ e $^{252}\text{Cf}$. Em contraste, o modelo de apenas evaporação subestima sistematicamente o rendimento de alta energia, mesmo dentro dos intervalos de confiança de 95%.
• Estabilidade Angular: O estudo confirma que propriedades assintóticas confiáveis só podem ser extraídas dentro do intervalo angular $110^\circ \leq \theta \leq 150^\circ$. Fora deste intervalo, efeitos de interferência transitórios e reflexões de fronteira impedem que o rendimento sature dentro do atual domínio de simulação.

Significância e Alegações
O artigo afirma fornecer a primeira identificação de uma assinatura de nêutron de escisão dentro dos espectros medidos de nêutrons de fissão pronta existentes utilizando um arcabouço totalmente microscópico. Ao demonstrar que a TDDFT produz naturalmente SNs sem suposições a priori, e que esses SNs resolvem a discrepância entre modelos de apenas evaporação e dados experimentais de alta energia, o estudo oferece evidência direta para um componente de nêutron de escisão não negligenciável na fissão pronta.

Os autores postulam que este trabalho estabelece uma fundação microscópica para interpretar espectros de nêutrons de fissão além do padrão estatístico-evaporação. Além disso, identificam uma assinatura experimental específica: a separação das contribuições de SN e EN abaixo do limiar de energia ($E_{thr}$). Isso implica que medições de baixa energia no intervalo angular especificado fornecem as restrições mais robustas para o componente de evaporação, enquanto a cauda de alta energia serve como o principal indicador da contribuição de SN. O estudo conclui que, embora as simulações atuais sejam limitadas a intervalos angulares específicos, a metodologia isola com sucesso uma assinatura mensurável de nêutrons de escisão que foi teoricamente conjecturada por quase um século.