Double Chooz publica a primeira medição quantitativa do fluxo residual de antineutrinos

Um reator nuclear desligado não fica silencioso. Produtos de fissão acumulados em seu núcleo e em suas piscinas continuam decaindo durante meses ou anos, e suas cadeias beta emitem antineutrinos. A colaboração Double Chooz publica a primeira medição quantitativa desse fluxo residual, obtida graças a desligamentos simultâneos dos dois reatores da usina de Chooz.

Fonte: journals.aps.org

Em termos simples

Quando se apaga uma fogueira, as brasas continuam incandescentes por muito tempo. Um reator nuclear faz algo semelhante, mas com um brilho que nenhum olho consegue ver: um fluxo de partículas fantasmagóricas, os antineutrinos, produzido pelas cadeias de decaimento beta dos produtos de fissão presentes no combustível irradiado que permanece nos núcleos e nos conjuntos usados armazenados em piscinas.

Durante 17,2 dias de tempo efetivo em que os dois reatores de Chooz, nas Ardenas, estavam desligados, o detector próximo — situado a cerca de 400 metros dos núcleos e protegido por 115 metros de equivalente de água — isolou 106 ± 18 candidatos residuais após a subtração dos fundos. Isso corresponde a cerca de seis por dia. O modelo detalhado prevê que 44 % do sinal venha das piscinas onde o combustível usado esfria; essa distribuição não é separada diretamente pelo detector.

A aplicação prática considerada é o controle de materiais nucleares: os antineutrinos atravessam a matéria com probabilidade muito baixa de interação, o que permite observar um reator sem acesso ao seu núcleo. Há, porém, uma ressalva: esse resultado depende de várias dezenas de metros cúbicos de cintilador e de uma blindagem subterrânea de 115 metros de equivalente de água. Ele valida um cálculo de fluxo residual; não demonstra que um dispositivo transportável funcionaria na superfície de uma usina.

Explicação técnica

1. Como se captura um antineutrino e por que são necessários muitos deles. A reação utilizada é o decaimento beta inverso (IBD): $\bar{\nu}_e + p \rightarrow e^+ + n$νˉe​+p→e++n. Um antineutrino eletrônico interage com um próton livre — um núcleo de hidrogênio do cintilador — e produz um nêutron e um pósitron. A assinatura é dupla, e é ela que torna a medição possível apesar de uma taxa de poucos eventos por dia. O pósitron deposita imediatamente sua energia cinética e depois se aniquila: esse é o sinal prompt. O nêutron se difunde antes de ser capturado; o método “total neutron capture” do Double Chooz explora as capturas no gadolínio, no hidrogênio e no carbono nos volumes internos para formar o sinal retardado. Exigir os dois sinais em uma janela de 0,5 a 800 µs e a menos de 1,2 m um do outro elimina a maior parte da radioatividade ambiente. A energia do antineutrino é reconstruída a partir do prompt por $E_{\bar{\nu}_e} \simeq E_{\text{prompt}} + 0{,}78\ \text{MeV}$Eνˉe​​≃Eprompt​+0,78 MeV, e a reação possui um limiar cinemático de cerca de 1,8 MeV: abaixo dele, esse canal não detecta o antineutrino.

2. De onde vem o brilho residual: três cadeias, três relógios. Um reator em operação obtém a maior parte de seu fluxo de antineutrinos de produtos de fissão de vida curta, cuja contribuição despenca nas horas seguintes ao desligamento. O restante detectável vem sobretudo de três pares pai-filho: durante a irradiação e depois durante o resfriamento, um pai de vida longa alimenta continuamente um filho de vida muito mais curta. Como o filho desaparece rapidamente após cada produção, sua atividade acompanha principalmente a do pai; o desligamento remove sobretudo o fluxo dominante que mascarava essas cadeias persistentes.

   Par	Meia-vida do pai	Filho: meia-vida e energia máxima $Q_\beta$Qβ​	Papel
   $^{144}\mathrm{Ce} \rightarrow {}^{144}\mathrm{Pr}$144Ce→144Pr	285 dias	17,3 min · 3,00 MeV	≈ 54 % do sinal de 1–3 MeV esperado pelo modelo
   $^{106}\mathrm{Ru} \rightarrow {}^{106}\mathrm{Rh}$106Ru→106Rh	372 dias	30,1 s · 3,54 MeV	≈ 38 % do sinal esperado pelo modelo
   $^{90}\mathrm{Sr} \rightarrow {}^{90}\mathrm{Y}$90Sr→90Y	28,9 anos	3,19 h · 2,28 MeV	torna-se dominante na previsão do modelo em longo prazo (> 90 % após uma década)

   Essas energias determinam todo o restante: com $Q_\beta$Qβ​ de 2 a 3,5 MeV e um limiar IBD de 1,8 MeV, apenas a cauda superior de cada espectro beta é detectável. Por isso, os autores esperam 98,7 % do fluxo residual detectável abaixo de 3 MeV — o que justifica restringir a análise à janela de 1–3 MeV em energia prompt e também explica sua dificuldade: é nessa região que os ruídos de fundo são mais densos.

3. O que a janela de 1–3 MeV mede exatamente e o controle que a valida. Restringir a análise a uma janela na qual o sinal é esperado expõe uma objeção clássica: encontra-se aquilo que se procura exatamente onde se procura. Os autores respondem com um controle negativo. Acima de 3 MeV, as taxas após a subtração do ruído de fundo são compatíveis com zero. O modelo prevê, sim, uma pequena cauda residual, mas abaixo de 0,1 evento por dia e, portanto, aquém da sensibilidade atual. A ausência de um excesso mensurável nessa região dá suporte ao modelo de fundo; não prova que o fluxo físico seja estritamente nulo ali.

4. Os fundos e por que a proteção subterrânea importa. Com seis eventos por dia, tudo conta. Três famílias são tratadas separadamente. Os acidentais — dois depósitos de energia sem relação, um tomado como prompt e o outro como retardado — são eliminados por uma rede neural cujo poder de rejeição supera 400; o que resta é então medido in situ a partir da taxa de depósitos isolados. Os correlacionados vêm dos múons cósmicos: podem ser nêutrons rápidos produzidos por espalação na rocha, que imitam um prompt e depois são capturados, ou isótopos β-n como o $^{9}\mathrm{Li}$9Li, gerado por espalação muônica no $^{12}\mathrm{C}$12C do cintilador, cujo decaimento fornece sozinho o par elétron-nêutron. Estes últimos são identificados pela correlação temporal e espacial com o múon progenitor. Os 115 m de equivalente de água de proteção do detector próximo reduzem o fluxo muônico e, portanto, a produção desses isótopos, que nenhum veto elimina perfeitamente.

5. A previsão: um cálculo de inventário, conjunto por conjunto. O número esperado não é uma ordem de grandeza, mas uma simulação detalhada. O espectro IBD previsto é escrito como
   $\frac{dN_{\bar{\nu}_e}}{dE_{\text{vis}}} = \sum_i \frac{N_p}{4\pi L_i^2}\iint \varphi_i(E,t)\,\sigma_{\text{IBD}}(E)\,P_{ee}(E,L_i)\,\mathcal{D}(E_{\text{vis}}|E)\,dE\,dt,$dEvis​dNνˉe​​​=i∑​4πLi2​Np​​∬φi​(E,t)σIBD​(E)Pee​(E,Li​)D(Evis​∣E)dEdt,
   em que $N_p$Np​ é o número de prótons-alvo, $L_i$Li​ é a distância até a fonte $i$i, $\varphi_i$φi​ é o espectro de emissão dependente do tempo, $\sigma_{\text{IBD}}$σIBD​ é a seção de choque, $P_{ee}$Pee​ é a probabilidade de sobrevivência por oscilação e $\mathcal{D}$D é a resposta do detector. As atividades dos produtos de fissão são obtidas pelos códigos APOLLO-2.8.4 e DARWIN-3 com a biblioteca de dados nucleares JEFF-3.1.1; os espectros beta vêm da biblioteca BESTIOLE, que aplica a teoria de Fermi com correções eletromagnéticas, tamanho finito do núcleo, blindagem atômica e fatores de forma para transições permitidas e proibidas. Cada conjunto de combustível teve seu histórico de irradiação simulado individualmente, com o acompanhamento de sua potência térmica durante toda a vida, até o início de cada período de desligamento.

6. Por que o modelo atribui 44 % às piscinas. As duas piscinas de armazenamento ficam próximas dos núcleos — portanto, vistas de um detector a cerca de 400 m, seu fator geométrico $1/L^2$1/L2 continua comparável. Cada conjunto chega ali depois de vários ciclos de irradiação e esfria durante anos. Um núcleo desligado contém um inventário recente que decai rapidamente; uma piscina acumula conjuntos de idades variadas, alguns com isótopos que persistem por muito tempo. Durante um desligamento prolongado, a soma prevista das contribuições das piscinas torna-se comparável à dos núcleos. O cálculo integrado fornece 56 % para os núcleos e 44 % para as piscinas; a geometria atual não permite medir separadamente essas duas frações.

7. O que limita a previsão e o que isso impede concluir. A incerteza de normalização do número previsto chega a 7,4 %, decomposta em 6,0 % de modelagem espectral, 2,9 % de geometria das linhas de base, 2,1 % de inventário dos produtos de fissão, 2,0 % de inventário dos conjuntos nas piscinas e menos de 0,8 % para o detector. O termo dominante decorre da forma da transição beta do $^{144}\mathrm{Pr}$144Pr, cuja incerteza é tomada como a diferença entre um cálculo detalhado de estrutura nuclear e a aproximação chamada $\xi$ξ. Esses 7,4 % não medem nem a precisão da concordância entre dados e modelo nem um desempenho de inferência de inventário: este estudo compara uma contagem a uma previsão, sem estimar inversamente uma massa ou uma composição de combustível.

Notação do diagrama: $N_{\rm obs}$Nobs​ é a contagem residual medida, $N_{\rm pred}$Npred​ é a previsão do modelo e $H_0(B)$H0​(B) é a hipótese de apenas fundo. Os $5{,}9\,\sigma$5,9σ comparam $N_{\rm obs}$Nobs​ a $H_0(B)$H0​(B), e não a $N_{\rm pred}$Npred​.

Por que funcionou

A medição depende de uma circunstância e de uma geometria. A circunstância: em 2017, os dois reatores de Chooz estiveram simultaneamente desligados quatro vezes, incluindo um período de 20,8 dias — condição indispensável, pois o detector próximo conta cerca de 900 eventos IBD por dia em potência máxima, contra 6,2 ± 1,1 eventos residuais por dia aqui. A geometria: um detector já construído, calibrado e compreendido, instalado a 400 m, cujas incertezas sistemáticas haviam sido estabelecidas por anos de medições de oscilação.

A comparação numérica é direta. Uma análise anterior, com os dados de 2011–2012, havia reunido apenas cerca de vinte candidatos — poucos demais para extrair um espectro residual. Aqui, após a subtração dos fundos, 106 ± 18 candidatos residuais são comparados aos 88 ± 7 esperados: a diferença é $\Delta = 18 \pm 19$Δ=18±19 eventos, isto é, uma concordância. No detector distante, são 27 ± 13 para 14 ± 1 esperados, com uma barra de erro que torna a conclusão muito menos restritiva.

Os limites são explícitos. A análise do sinal cobre a janela de 1–3 MeV. Acima de 3 MeV, as taxas subtraídas são compatíveis com zero, enquanto o modelo mantém uma pequena cauda residual até cerca de 4,5 MeV, inferior a 0,1 evento por dia e não resolvida por esta amostra. A incerteza de normalização de 7,4 % caracteriza a previsão; o experimento não demonstra a precisão de uma inferência de inventário.

A diferença entre o que foi demonstrado e o que é considerado precisa ficar clara. O que foi demonstrado: um fluxo residual é distinguido do fundo a $5{,}9\,\sigma$5,9σ, e a contagem $106\pm18$106±18 é compatível com a previsão $88\pm7$88±7, pois $\Delta=18\pm19$Δ=18±19. O que não foi demonstrado: uma vigilância operacional do combustível. A demonstração mobiliza mais de 30 m³ de cintilador sob 115 m de equivalente de água para cerca de seis eventos por dia a 400 m. Um instrumento compacto na superfície ainda teria de controlar o ruído cósmico e demonstrar sensibilidade a mudanças de inventário.

Cadeia causal

Fissão no núcleo → acumulação, durante a irradiação, de pais de vida longa ($^{144}\mathrm{Ce}$144Ce, $^{106}\mathrm{Ru}$106Ru, $^{90}\mathrm{Sr}$90Sr) que alimentam filhas beta de vida curta → desligamento dos reatores e rápido desaparecimento do fluxo dominante → persistência das cadeias já estabelecidas nos núcleos e nos conjuntos transferidos para as piscinas → quatro desligamentos simultâneos totalizando 24,4 dias em 2017 → 17,2 dias de tempo efetivo no detector próximo → $106\pm18$106±18 candidatos residuais entre 1 e 3 MeV, incompatíveis com apenas fundo a $5{,}9\,\sigma$5,9σ → pequena cauda prevista acima de 3 MeV, inferior a 0,1 evento por dia e não resolvida → compatibilidade com o modelo conjunto por conjunto ($88\pm7$88±7, $\Delta=18\pm19$Δ=18±19) → primeira referência quantitativa para testar os cálculos de fluxo residual → estudos futuros de separação entre núcleo e piscina e de detectores dedicados.

Curiosidade

O Double Chooz já havia tentado essa medição com seus dados de 2011–2012. Na ocasião, a colaboração reuniu apenas cerca de vinte candidatos residuais: o detector próximo ainda não estava disponível, e o tempo de desligamento simultâneo era curto demais para extrair um espectro. Em 2017, o detector próximo estava em operação, e quatro desligamentos simultâneos, incluindo um de 20,8 dias, forneceram 17,2 dias de tempo efetivo; essas duas mudanças tornaram possível a medição quantitativa.

Legado e dados atuais

A medição fornece o primeiro ponto de comparação quantitativo publicado entre um cálculo de soma do fluxo residual e uma contagem real. A previsão tem uma incerteza de normalização de 7,4 %; a comparação observada é $\Delta=18\pm19$Δ=18±19 eventos. Essas duas grandezas não devem ser confundidas.

Na vertente aplicada, a distribuição prevista é a informação nova mais diretamente utilizável: o modelo atribui 44 % do sinal do local desligado às piscinas, não aos núcleos. Um método de vigilância que modelasse apenas o combustível ainda presente nos núcleos subestimaria fortemente o sinal esperado. O que a publicação não permite quantificar, porém, é a sensibilidade que um detector compacto na superfície alcançaria: o artigo não contém dados desse tipo, e os autores não fazem a extrapolação a partir dos volumes ativos do dispositivo subterrâneo.

O olhar do pesquisador — questões em aberto

Estas possibilidades são desdobramentos potenciais, não resultados do estudo.

• O teste espectral decisivo. Com mais estatística, acompanhar durante vários meses a evolução relativa do $^{144}\mathrm{Ce}$144Ce e do $^{106}\mathrm{Ru}$106Ru testaria a forma temporal do modelo. Observar o domínio do par $^{90}\mathrm{Sr}$90Sr–$^{90}\mathrm{Y}$90Y exigiria anos: a previsão só ultrapassa 90 % depois de mais de dez anos.
• Separar as piscinas dos núcleos pela geometria. Os 56 %/44 % são aqui uma previsão do modelo, não uma separação medida: a 400 m, uma diferença de 38 m não pode ser resolvida. Um detector colocado bem mais perto, ou dois detectores em ângulos diferentes, poderia separar experimentalmente as duas contribuições em vez de distribuí-las por cálculo.
• O obstáculo da estrutura nuclear. A incerteza dominante não vem do detector, mas da forma da transição beta do $^{144}\mathrm{Pr}$144Pr. Uma medição dedicada desse espectro beta poderia reduzir o termo de 6,0 %, mas o estudo não quantifica a precisão final que o modelo alcançaria nesse caso.

Fontes

Referências verificadas durante a auditoria de checagem de fatos de 16 de agosto de 2026.

• Double Chooz Collaboration, “First Measurement of Neutrino Emissions from Spent Nuclear Fuel by the Double Chooz Experiment”, Physical Review Letters 137, 061803, 4 de agosto de 2026 — artigo revisado por pares. DOI: 10.1103/dr26-j19g
• Versão preprint do mesmo trabalho: arXiv:2510.04869 — texto integral consultado (resumo, métodos, tabelas de incertezas sistemáticas).

Referências de base

1. P. Vogel e J. F. Beacom — Angular distribution of neutron inverse beta decay — cinemática e seção de choque de referência para o decaimento beta inverso.

Declaração de confiança

Confiança alta nas contagens, significâncias, janelas de energia, geometria e componentes de incerteza reproduzidos do artigo revisado por pares. Confiança média no alcance aplicado: a publicação menciona as salvaguardas nucleares, mas não testa um detector compacto na superfície nem a inferência de inventário. A parcela de 44 % atribuída às piscinas é uma saída do modelo, não uma separação medida. Esta primeira análise quantitativa ainda não teve replicação independente. Controle de unicidade: nenhum DOI, identificador arXiv ou resultado central equivalente aparece no registro A001–A018; veredito UNIQUE. DOI, datas, fórmulas, links e ausência de duplicata foram verificados.