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

Este estudio utiliza la teoría del funcional de la densidad dependiente del tiempo para demostrar que los neutrones de escisión constituyen un componente significativo de alta energía del espectro de neutrones de fisión prompt en la fisión de $^{235}\mathrm{U}$, $^{239}\mathrm{Pu}$ y $^{252}\mathrm{Cf}$, proporcionando evidencia directa de que su inclusión resuelve la subestimación sistemática de los rendimientos de alta energía observada en los modelos tradicionales de solo evaporación.

Lo esencial

Imagina un globo gigante e inestable (un núcleo atómico) que de repente se rompe en dos. Esto es la fisión nuclear. Durante décadas, los científicos han sabido que cuando esto sucede, las dos piezas resultantes (llamadas fragmentos) salen disparadas a velocidades increíbles y escupen diminutas partículas llamadas neutrones.

Durante mucho tiempo, los físicos pensaron que todos estos neutrones se "evaporaban" más tarde, como el vapor que sube de una taza de café caliente después de que el agua ha hervido. Asumían que los fragmentos ya estaban completamente formados y moviéndose de manera constante antes de empezar a soltar estos neutrones.

Sin embargo, este nuevo artículo sugiere que algunos neutrones son en realidad "expulsados" justo en el momento en que el globo se rompe. Estos son los neutrones de escisión. Nacen en el caos de un segundo en el momento de la ruptura, no de un fragmento tranquilo y enfriándose más tarde.

Aquí es cómo los investigadores encontraron la prueba de estos neutrones del "momento del chasquido", explicado de forma sencilla:

1. La simulación de supercomputadora

Para ver qué sucede durante la división, los científicos no usaron un microscopio; usaron una supercomputadora para ejecutar una película del evento utilizando una teoría llamada Teoría de la Funcional de la Densidad Dependiente del Tiempo (TDDFT).

Piensa en esto como un videojuego 3D de alta velocidad donde simulan los átomos bailando y rompiéndose. En versiones anteriores de este "juego", el mundo virtual era demasiado pequeño. Los neutrones chocaban contra el borde de la pantalla antes de que los científicos pudieran determinar exactamente qué tan rápido iban o hacia qué dirección volaban.

En este estudio, construyeron un mundo virtual mucho más grande (unas 3 veces más grande que antes). Esto les dio a los neutrones suficiente espacio para salir volando y asentarse, de modo que pudieran medir sus velocidades y direcciones con precisión sin que las "paredes" de la simulación alteraran los datos.

2. El descubrimiento del "límite de velocidad"

Una vez que tuvieron una visión clara, observaron los neutrones que salían en ángulos específicos (principalmente hacia los lados y ligeramente hacia atrás respecto a la división). Encontraron algo sorprendente:

• La zona de "prohibido el paso": No hay ningún neutrón de escisión con baja energía (por debajo de aproximadamente 1.5 a 2 millones de electronvoltios). Es como si hubiera un límite de velocidad; nada lento puede ser un neutrón de "momento de escisión".
• La multitud de alta velocidad: En cambio, estos neutrones son todos rápidos. Se agrupan alrededor de una velocidad específica alta (3–3.5 MeV) y luego se desvanecen en una larga cola de partículas aún más rápidas.

Es como una multitud de personas saltando desde un trampolín. Los neutrones "evaporados" son como personas que caminan tranquilamente fuera de la plataforma de la piscina más tarde. Los neutrones de "escisión" son como personas que son lanzadas violentamente desde el trampolín en el segundo exacto en que este se rompe. Los que son lanzados desde el trampolín siempre se mueven rápido; nunca verás a uno lento proveniente de ese evento específico.

3. Resolviendo el misterio de la energía "faltante"

Los científicos han intentado durante años hacer coincidir sus modelos computacionales con experimentos del mundo real. Tenían un problema:

• El modelo antiguo: Si solo cuentas el "vapor" (neutrones evaporados), tu modelo computacional predice muy pocos neutrones de alta energía. Es como intentar llenar un cubo con una taza pequeña, pero el cubo necesita más agua de la que la taza puede proporcionar.
• El nuevo modelo: Cuando los investigadores añadieron los neutrones del "momante del chasquido" (los que encontraron en su gran simulación) a los neutrones de "vapor", las matemáticas finalmente funcionaron. El modelo combinado coincidió perfectamente con los datos de alta energía medidos en experimentos reales para el Uranio y el Californio.

4. Por qué esto es importante

Esto es algo importante porque es la primera vez que una teoría puramente microscópica (una que no solo supone o asume que las cosas existen) ha predicho estos neutrones de "momento de escisión" y ha demostrado que son reales.

• Antes: Los científicos tenían que suponer que estos neutrones existían porque las matemáticas no cuadraban.
• Ahora: La simulación por computadora produjo naturalmente estos neutrones sin que se le ordenara que existieran. Es como predecir una tormenta observando el movimiento de las nubes, en lugar de simplemente asumir que viene una tormenta porque el reporte meteorológico lo dice.

La conclusión

El artículo conclifica que cuando un átomo se divide, una pequeña pero importante parte de los neutrones (aproximadamente del 6% al 10% del total) nacen en el violento momento de la ruptura. Estos neutrones son distintos porque siempre son rápidos y nunca lentos dentro de ciertos ángulos.

Al encontrar esta "huella digital" en los datos, los investigadores han logrado finalmente separar los neutrones del "momento del chasquido" de los neutrones de "vapor", dándonos una imagen más clara y precisa de cómo funciona realmente la fisión nuclear. Esto ayuda a refinar nuestra comprensión de las fuerzas fundamentales que mantienen la materia unida y la desgarran.

Resumen técnico

Resumen Técnico: Propiedades Angulares y Cinéticas de los Neutrones de Escisión dentro de la Teoría de la Funcional de la Densidad Dependiente del Tiempo

Planteamiento del Problema
A pesar del descubrimiento de la fisión nuclear hace casi un siglo, la descripción microscópica completa del proceso sigue siendo esquiva, particularmente en lo que respecta a la transición desde el punto de silla exterior hacia la escisión y la formación subsiguiente de los fragmentos de fisión (FF) primarios. Los códigos de evaporación estadística actuales asumen que todos los neutrones rápidos son emitidos isotrópicamente desde fragmentos totalmente acelerados. Sin embargo, este marco no puede explicar los neutrones emitidos durante el proceso dinámico altamente fuera de equilibrio de la ruptura del cuello, conocidos como neutrones de escisión (SN). Aunque los SN han sido conjeturados desde 1939 y estudiados mediante diversos modelos, permanecen sin resolverse experimentalmente porque las observables de neutrones rápidos por sí solas no pueden distinguir entre los neutrones emitidos en la escisión y aquellos evaporados posteriormente. Los intentos previos para predecir los SN utilizando la teoría de la funcional de la densidad dependiente del tiempo (TDDFT) se vieron limitados por dominios de simulación pequeños, lo que impedía la determinación simultánea de las propiedades angulares y cinéticas antes de que los neutrones alcanzaran los límites computacionales.

Metodología
Este estudio emplea un marco plenamente microscópico basado en TDDFT extendido a sistemas superfluidos para investigar la emisión de SN en $^{235}\text{U}(n_{th}, f)$, $^{239}\text{Pu}(n_{th}, f)$ y $^{252}\text{Cf}(sf)$.

• Configuración de la Simulación: Los autores utilizaron el resolvedor de oscilador armónico axialmente simétrico AxialHOHFB con el funcional de densidad de energía (EDF) SkM*. Un avance crítico fue el uso de un dominio de simulación sustancialmente mayor (un elipsoide con semiejes $Z_{max} = 70$ fm y $R_{max} = 55$ fm), aproximadamente 3.2 veces más grande que en estudios previos.
• Integración en Dos Pasos: Para gestionar las restricciones computacionales manteniendo la precisión, se empleó un procedimiento de integración temporal de dos pasos. El sistema se evolucionó desde la solución estática de Hartree-Fock-Bogoliubov (HFB) con restricción hasta el momento de la escisión en una base más pequeña ($B_1$), y luego se transformó a la base más grande ($B_2$) para continuar la evolución hasta que los neutrones se aproximaran al límite.
• Extracción de las Propiedades de los SN: Los SN se identificaron definiendo una región $V^{(SN)}_\eta$ donde la relación entre la densidad de energía del flujo colectivo y la densidad de energía cinética intrínseca excedía un umbral $\eta = 3$. Este criterio aísla los neutrones dominados por el flujo colectivo hacia el exterior.
• Análisis Angular y Cinético: Dentro de esta región, se calcularon los campos de velocidad local y la energía de flujo colectivo por neutrón. El estudio se centró en el intervalo angular $110^\circ \leq \theta \leq 150^\circ$ (medido respecto a la dirección del fragmento ligero), donde el rendimiento de los SN se satura antes de que ocurran las reflexiones de los límites.
• Comparación con el Experimento: El espectro de energía cinética de los SN calculado se combinó con un modelo de evaporación de Maxwell para el componente de neutrón evaporado (EN). Los parámetros del modelo de evaporación (temperatura del fragmento $T$ y la relación de multiplicidad $\bar{\nu}_L/\bar{\nu}_H$) fueron restringidos exclusivamente por datos experimentales de baja energía ($E \leq E_{thr}$). Este modelo combinado se comparó contra los datos del espectro de neutrones de fisión rápida (PFNS) de alta energía para $^{239}\text{Pu}$ y $^{252}\text{Cf}$.

Resultados Clave

• Características del Espectro de los SN: Los espectros de los SN extraídos para los tres sistemas exhiben un comportamiento de umbral distintivo. Los SN están ausentes por debajo de un umbral de energía específico ($E_{thr}$), el cual es aproximadamente de 1.5–2 MeV (específicamente $\approx 2$ MeV para la fisión inducida y $\approx 1.5$ MeV para la fisión espontánea). Los espectros alcanzan su máximo cerca de 3–3.5 MeV y presentan una cola exponencial de alta energía que se extiende hasta aproximadamente 18 MeV.
• Estimaciones de Rendimiento: Para $^{239}\text{Pu}(n_{th}, f)$, el rendimiento total de los SN se calcula como $0.19 \pm 0.01$ (aproximadamente el 6.6% de los neutrones rápidos totales). Para $^{252}\text{Cf}(sf)$, el rendimiento es de $0.39 \pm 0.02$ (aproximadamente el 10.4%).
• Validación del Modelo: Cuando se añade el componente calculado de los SN al modelo de evaporación (restringido solo por datos de baja energía), el PFNS total resultante reproduce con precisión el rendimiento medido de alta energía tanto para $^{239}\text{Pu}$ como para $^{252}\text{Cf}$. En contraste, el modelo de solo evaporación subestima sistemáticamente el rendimiento de alta energía, incluso dentro de los intervalos de confianza del 95%.
• Estabilidad Angular: El estudio confirma que las propiedades asintóticas fiables solo pueden extraerse dentro del rango angular $110^\circ \leq \theta \leq 150^\circ$. Fuera de este rango, los efectos de interferencia transitorios y las reflexiones de los límites impiden que el rendimiento se sature dentro del dominio actual de la simulación.

Significación y Reivindicaciones
El artículo afirma proporcionar la primera identificación de una firma de neutrones de escisión dentro de los espectros medidos de neutrones de fisión rápida existentes mediante un marco plenamente microscópico. Al demostrar que la TDDFT produce naturalmente los SN sin supuestos a priori, y que estos SN resuelven la discrepancia entre los modelos de solo evaporación y los datos experimentales de alta energía, el estudio ofrece evidencia directa de un componente de neutrones de escisión no despreciable en la fisión rápida.

Los autores postulan que este trabajo establece una base microscópica para interpretar los espectros de neutrones de fisión más allá del esquema estándar de evaporación estadística. Además, identifican una firma experimental específica: la separación de las contribuciones de SN y EN por debajo del umbral de energía ($E_{thr}$). Esto implica que las mediciones de baja energía en el rango angular especificado proporcionan las restricciones más robustas para el componente de evaporación, mientras que la cola de alta energía sirve como el indicador primario de la contribución de los SN. El estudio concluye que, si bien las simulaciones actuales están limitadas a rangos angulares específicos, la metodología logra aislar una firma medible de los neutrones de escisión que ha sido teóricamente conjeturada durante casi un siglo.