Double Chooz publica la primera medición cuantitativa del flujo residual de antineutrinos

Un reactor nuclear apagado no queda en silencio. Los productos de fisión acumulados en su núcleo y en sus piscinas siguen decayendo durante meses o años, y sus cadenas beta emiten antineutrinos. La colaboración Double Chooz publica la primera medición cuantitativa de este flujo residual, obtenida gracias a paradas simultáneas de los dos reactores de la central nuclear de Chooz.

Fuente: journals.aps.org

En pocas palabras

Cuando se apaga un fuego de leña, las brasas siguen al rojo durante mucho tiempo. Un reactor nuclear hace lo mismo, pero con un resplandor que ningún ojo puede ver: un flujo de partículas fantasmales, los antineutrinos, producido por las cadenas de desintegración beta de los productos de fisión presentes en el combustible irradiado que permanece en los núcleos y en los conjuntos gastados almacenados en piscinas.

Durante 17,2 días de tiempo vivo en los que los dos reactores de Chooz, en las Ardenas francesas, estuvieron apagados, el detector cercano —situado a unos 400 metros de los núcleos y protegido por 115 metros de equivalente en agua— aisló 106 ± 18 candidatos residuales después de sustraer los fondos. Son unos seis al día. El modelo detallado predice que 44 % de la señal procede de las piscinas donde se enfría el combustible gastado; el detector no separa directamente esta fracción.

La aplicación práctica considerada es el control de los materiales nucleares: los antineutrinos atraviesan la materia con una probabilidad de interacción muy baja, lo que permite observar un reactor sin acceder a su núcleo. Conviene, no obstante, poner una cautela: este resultado se apoya en varias decenas de metros cúbicos de centelleador y una protección subterránea de 115 metros de equivalente en agua. Valida un cálculo del flujo residual; no demuestra que un dispositivo transportable funcionaría en la superficie de una central nuclear.

Explicación técnica

1. Cómo se captura un antineutrino y por qué hacen falta tantos. La reacción utilizada es la desintegración beta inversa (IBD): $\bar{\nu}_e + p \rightarrow e^+ + n$νˉe​+p→e++n. Un antineutrino electrónico interactúa con un protón libre —un núcleo de hidrógeno del centelleador— y produce un neutrón y un positrón. La firma es doble, y eso permite la medición pese a una tasa de apenas unos sucesos al día. El positrón deposita de inmediato su energía cinética y después se aniquila: es la señal prompt. El neutrón se dispersa antes de ser capturado; el método «total neutron capture» de Double Chooz aprovecha las capturas en gadolinio, hidrógeno y carbono dentro de los volúmenes internos para formar la señal retardada. Exigir las dos señales en una ventana de 0,5 a 800 µs y a menos de 1,2 m una de otra elimina la mayor parte de la radiactividad ambiental. La energía del antineutrino se reconstruye a partir de la señal prompt mediante $E_{\bar{\nu}_e} \simeq E_{\text{prompt}} + 0{,}78\ \text{MeV}$Eνˉe​​≃Eprompt​+0,78 MeV, y la reacción posee un umbral cinemático de unos 1,8 MeV: por debajo, este canal no detecta el antineutrino.

2. De dónde procede el resplandor residual: tres cadenas, tres relojes. Un reactor en funcionamiento obtiene la mayor parte de su flujo de antineutrinos de productos de fisión de vida corta, cuya contribución se desploma en las horas posteriores a la parada. El remanente detectable procede sobre todo de tres parejas padre–hija: durante la irradiación y después el enfriamiento, un padre de vida larga alimenta continuamente a una hija de vida mucho más corta. Como la hija desaparece con rapidez después de cada producción, su actividad sigue principalmente la del padre; la parada elimina ante todo el flujo dominante que ocultaba estas cadenas persistentes.

   Pareja	Semivida del padre	Hija: semivida y energía máxima $Q_\beta$Qβ​	Papel
   $^{144}\mathrm{Ce} \rightarrow {}^{144}\mathrm{Pr}$144Ce→144Pr	285 días	17,3 min · 3,00 MeV	≈ 54 % de la señal de 1–3 MeV esperada por el modelo
   $^{106}\mathrm{Ru} \rightarrow {}^{106}\mathrm{Rh}$106Ru→106Rh	372 días	30,1 s · 3,54 MeV	≈ 38 % de la señal esperada por el modelo
   $^{90}\mathrm{Sr} \rightarrow {}^{90}\mathrm{Y}$90Sr→90Y	28,9 años	3,19 h · 2,28 MeV	se vuelve dominante en la predicción del modelo a largo plazo (> 90 % después de una década)

   Estas energías condicionan todo lo demás: con valores de $Q_\beta$Qβ​ de 2 a 3,5 MeV y un umbral IBD de 1,8 MeV, solo puede detectarse la cola de alta energía de cada espectro beta. En consecuencia, los autores esperan 98,7 % del flujo residual detectable por debajo de 3 MeV; esto justifica restringir el análisis a la ventana de energía prompt de 1–3 MeV y también lo dificulta, porque ahí es donde los fondos son más densos.

3. Qué mide exactamente la ventana de 1–3 MeV y el control que la valida. Restringir el análisis a una ventana donde se espera la señal expone a una objeción clásica: se encuentra lo que se busca allí donde se ha decidido buscar. Los autores responden con un control negativo. Por encima de 3 MeV, los rendimientos tras sustraer el fondo son compatibles con cero. El modelo sí predice una pequeña cola residual, pero inferior a 0,1 suceso al día y, por tanto, por debajo de la sensibilidad actual. La ausencia de un exceso medible en esta zona respalda el modelo de fondo; no demuestra que el flujo físico sea allí estrictamente nulo.

4. Los fondos y por qué importa la protección subterránea. Con seis sucesos al día, todo cuenta. Se tratan por separado tres familias. Los accidentales —dos depósitos de energía no relacionados, uno confundido con una señal prompt y el otro con una retardada— se suprimen mediante una red neuronal cuyo poder de rechazo supera 400; lo que queda se mide después in situ a partir de la tasa de depósitos aislados. Los fondos correlacionados proceden de los muones cósmicos: bien de neutrones rápidos producidos por espalación en la roca, que imitan una señal prompt antes de ser capturados, bien de isótopos β-n como el $^{9}\mathrm{Li}$9Li, generado por espalación muónica sobre el $^{12}\mathrm{C}$12C del centelleador, cuya desintegración aporta por sí sola la pareja electrón–neutrón. Estos últimos se etiquetan por su correlación temporal y espacial con el muón progenitor. Los 115 m de equivalente en agua que protegen el detector cercano reducen el flujo muónico y, con él, la producción de esos isótopos, que ningún veto elimina por completo.

5. La predicción: un cálculo de inventario, conjunto por conjunto. El número esperado no es una estimación de orden de magnitud, sino una simulación detallada. El espectro IBD predicho se escribe
   $\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,
   donde $N_p$Np​ es el número de protones blanco, $L_i$Li​ la distancia a la fuente $i$i, $\varphi_i$φi​ el espectro de emisión dependiente del tiempo, $\sigma_{\text{IBD}}$σIBD​ la sección eficaz, $P_{ee}$Pee​ la probabilidad de supervivencia por oscilación y $\mathcal{D}$D la respuesta del detector. Las actividades de los productos de fisión proceden de los códigos APOLLO-2.8.4 y DARWIN-3 con la biblioteca de datos nucleares JEFF-3.1.1; los espectros beta provienen de la biblioteca BESTIOLE, que aplica la teoría de Fermi con correcciones electromagnéticas, de tamaño nuclear finito y apantallamiento atómico, además de factores de forma para las transiciones permitidas y prohibidas. La historia de irradiación de cada conjunto de combustible se simuló individualmente, siguiendo su potencia térmica durante toda su vida hasta el comienzo de cada periodo de parada.

6. Por qué el modelo atribuye 44 % a las piscinas. Las dos piscinas de almacenamiento están cerca de los núcleos, de modo que, vistas desde un detector situado a unos 400 m, su factor geométrico $1/L^2$1/L2 sigue siendo comparable. Cada conjunto llega allí tras varios ciclos de irradiación y se enfría durante años. Un núcleo apagado contiene un inventario reciente que decae deprisa; una piscina acumula conjuntos de edades escalonadas, algunos de cuyos isótopos persisten mucho tiempo. Durante una parada prolongada, la suma predicha de las contribuciones de las piscinas se vuelve comparable a la de los núcleos. El cálculo integrado da 56 % para los núcleos y 44 % para las piscinas; la geometría actual no permite medir por separado estas dos fracciones.

7. Qué limita la predicción y qué impide concluir. La incertidumbre de normalización del número predicho alcanza 7,4 %, desglosada en 6,0 % de modelización espectral, 2,9 % de geometría de las líneas de base, 2,1 % de inventario de productos de fisión, 2,0 % de inventario de los conjuntos en las piscinas y menos de 0,8 % para el detector. El término dominante se debe a la forma de la transición beta del $^{144}\mathrm{Pr}$144Pr, cuya incertidumbre se toma como la diferencia entre un cálculo detallado de estructura nuclear y la llamada aproximación $\xi$ξ. Estos 7,4 % no miden ni la precisión del acuerdo entre datos y modelo ni el rendimiento de una inferencia de inventario: este estudio compara un recuento con una predicción, sin estimar a la inversa una masa o una composición de combustible.

Notación del esquema: $N_{\rm obs}$Nobs​ es el recuento residual medido, $N_{\rm pred}$Npred​ la predicción del modelo y $H_0(B)$H0​(B) la hipótesis de solo fondo. Los $5{,}9\,\sigma$5,9σ comparan $N_{\rm obs}$Nobs​ con $H_0(B)$H0​(B), no con $N_{\rm pred}$Npred​.

Por qué funcionó

La medición depende de una circunstancia y una geometría. La circunstancia: en 2017, los dos reactores de Chooz estuvieron apagados simultáneamente cuatro veces, una de ellas durante un periodo de 20,8 días, condición indispensable porque el detector cercano registra unos 900 sucesos IBD al día a plena potencia, frente a 6,2 ± 1,1 sucesos residuales diarios aquí. La geometría: un detector ya construido, calibrado y bien conocido, instalado a 400 m, cuyos efectos sistemáticos se habían establecido durante años de mediciones de oscilación.

La comparación numérica es directa. Un análisis anterior, con los datos de 2011–2012, solo había reunido una veintena de candidatos, demasiado pocos para extraer un espectro residual. Aquí, después de sustraer los fondos, se comparan 106 ± 18 candidatos residuales con 88 ± 7 esperados: la diferencia es de $\Delta = 18 \pm 19$Δ=18±19 sucesos, es decir, hay acuerdo. En el detector lejano, 27 ± 13 frente a 14 ± 1 esperados, con una barra de error que deja una conclusión mucho menos restrictiva.

Los límites son explícitos. El análisis de la señal abarca la ventana de 1–3 MeV. Por encima de 3 MeV, los rendimientos sustraídos son compatibles con cero, mientras que el modelo conserva una pequeña cola residual hasta alrededor de 4,5 MeV, inferior a 0,1 suceso al día y no resuelta en esta muestra. La incertidumbre de normalización de 7,4 % caracteriza la predicción; el experimento no demuestra la precisión de una inferencia de inventario.

Merece explicitarse la distancia entre lo demostrado y lo considerado. Lo demostrado: un flujo residual se distingue del fondo a $5{,}9\,\sigma$5,9σ, y el recuento $106\pm18$106±18 es compatible con la predicción $88\pm7$88±7 porque $\Delta=18\pm19$Δ=18±19. Lo no demostrado: una vigilancia operativa del combustible. La demostración emplea más de 30 m³ de centelleador bajo 115 m de equivalente en agua, para unos seis sucesos al día a 400 m. Un instrumento compacto en superficie aún tendría que controlar el fondo cósmico y demostrar su sensibilidad a cambios de inventario.

Cadena causal

Fisión en el núcleo → acumulación, durante la irradiación, de padres de vida larga ($^{144}\mathrm{Ce}$144Ce, $^{106}\mathrm{Ru}$106Ru, $^{90}\mathrm{Sr}$90Sr) que alimentan hijas beta de vida corta → parada de los reactores y rápida desaparición del flujo dominante → persistencia de las cadenas ya establecidas en los núcleos y en los conjuntos transferidos a las piscinas → cuatro paradas simultáneas que sumaron 24,4 días en 2017 → 17,2 días de tiempo vivo en el detector cercano → $106\pm18$106±18 candidatos residuales entre 1 y 3 MeV, incompatibles con el solo fondo a $5{,}9\,\sigma$5,9σ → pequeña cola predicha por encima de 3 MeV, inferior a 0,1 suceso al día y no resuelta → compatibilidad con el modelo conjunto por conjunto ($88\pm7$88±7, $\Delta=18\pm19$Δ=18±19) → primera referencia cuantitativa para probar los cálculos de flujo residual → futuros estudios para separar núcleos y piscinas y para usar detectores dedicados.

Anécdota

Double Chooz ya había intentado esta medición con sus datos de 2011–2012. La colaboración solo había reunido entonces una veintena de candidatos residuales: el detector cercano todavía no estaba disponible y el tiempo de parada simultánea era demasiado corto para extraer un espectro. En 2017, el detector cercano funcionaba y cuatro paradas simultáneas, una de ellas de 20,8 días, aportaron 17,2 días de tiempo vivo; esos dos cambios hicieron posible la medición cuantitativa.

Legado y datos actuales

La medición ofrece el primer punto de comparación cuantitativo publicado entre un cálculo de sumación del flujo residual y un recuento real. La predicción tiene una incertidumbre de normalización de 7,4 %; la comparación observada es, por su parte, $\Delta=18\pm19$Δ=18±19 sucesos. Estas dos cantidades no deben confundirse.

En cuanto a las aplicaciones, el reparto predicho es la información nueva con utilidad más directa: el modelo atribuye 44 % de la señal del emplazamiento apagado a las piscinas, no a los núcleos. Por tanto, un método de vigilancia que solo modelizara el combustible todavía presente en los núcleos subestimaría considerablemente la señal esperada. Lo que la publicación no permite cuantificar, en cambio, es la sensibilidad que alcanzaría un detector compacto en superficie: el artículo no contiene datos de este tipo y los autores no extrapolan a partir de los volúmenes activos del dispositivo subterráneo.

La mirada del investigador — preguntas abiertas

Estas líneas son posibles prolongaciones, no resultados del estudio.

• La prueba espectral decisiva. Con más estadística, seguir durante varios meses la evolución relativa del $^{144}\mathrm{Ce}$144Ce y del $^{106}\mathrm{Ru}$106Ru pondría a prueba la forma temporal del modelo. Observar el dominio de la pareja $^{90}\mathrm{Sr}$90Sr–$^{90}\mathrm{Y}$90Y exigiría años: la predicción solo supera 90 % después de más de diez años.
• Separar las piscinas de los núcleos mediante la geometría. Los 56 %/44 % son aquí una predicción del modelo, no una separación medida: a 400 m, una diferencia de 38 m no puede resolverse. Un detector situado mucho más cerca, o dos detectores en ángulos diferentes, podría disociar experimentalmente las dos contribuciones en lugar de repartirlas por cálculo.
• El cuello de botella de la estructura nuclear. La incertidumbre dominante no procede del detector, sino de la forma de la transición beta del $^{144}\mathrm{Pr}$144Pr. Una medición específica de este espectro beta podría reducir el término de 6,0 %, pero el estudio no cuantifica la precisión final que alcanzaría entonces el modelo.

Fuentes

Referencias verificadas durante la auditoría de comprobación de hechos del 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 — artículo revisado por pares. DOI: 10.1103/dr26-j19g
• Versión preprint del mismo trabajo: arXiv:2510.04869 — texto completo consultado (resumen, métodos y tablas de incertidumbres sistemáticas).

Referencias de contexto

1. P. Vogel y J. F. Beacom — Angular distribution of neutron inverse beta decay — cinemática y sección eficaz de referencia para la desintegración beta inversa.

Declaración de confianza

Confianza alta en los recuentos, las significancias, las ventanas de energía, la geometría y los componentes de incertidumbre reproducidos a partir del artículo revisado por pares. Confianza media en el alcance práctico: la publicación aborda las salvaguardias nucleares, pero no prueba ni un detector compacto en superficie ni la inferencia de inventario. La fracción de 44 % atribuida a las piscinas es una salida del modelo, no una separación medida. Este primer análisis cuantitativo aún no cuenta con una réplica independiente. Control de unicidad: ningún DOI, identificador arXiv ni resultado central equivalente aparece en el registro A001–A018; veredicto UNIQUE. Se comprobaron el DOI, las fechas, las fórmulas, los enlaces y la ausencia de duplicado.