Nuclear molecule of heavy nuclei

Resumen Técnico: Molécula Nuclear de Núcleos Pesados

Planteamiento del Problema
El artículo aborda la descripción teórica de núcleos pesados como sistemas moleculares compuestos por dos fragmentos que interactúan (un sistema dinuclear, o DNS). Si bien el concepto de moléculas nucleares se ha aplicado con éxito a fragmentos ligeros (p. ej., cúmulos $\alpha$) y a la desintegración alfa, los intentos previos de formular un Hamiltoniano para fragmentos pesados a menudo dependieron de aproximaciones inconsistentes. Específicamente, modelos anteriores (Refs. [9, 10]) reemplazaron el Hamiltoniano de la energía cinética total, incluyendo los términos de Coriolis, por un Hamiltoniano de rotador asimétrico. Los autores argumentan que esta aproximación conduce a un momento de inercia rotacional subestimado y no describe correctamente la evolución de los fragmentos pesados, particularmente cuando uno de los fragmentos está fuertemente deformado. El objetivo es derivar un Hamiltoniano autoconsistente que describa con precisión la dinámica colectiva (rotacional y vibracional) de un sistema dinuclear compuesto por un fragmento pesado, axialmente simétrico y deformado, y un fragmento esférico.

Metodología
Los autores desarrollan un modelo de mecánica cuántica basado en los siguientes pasos:

1. Sistemas de Coordenadas y Energía Cinética:

   • El sistema se define mediante un marco de laboratorio y marcos intrínsecos para los dos fragmentos. El vector de distancia relativa $\mathbf{R}$ conecta los centros de masa, y la orientación del fragmento deformado se define mediante ángulos de Euler.
   • Se deriva una expresión clásica de la energía cinética separando el movimiento del centro de masa, el movimiento relativo de los fragmentos y la rotación interna del fragmento deformado.
   • Utilizando el procedimiento de cuantización de Pauli, la energía cinética clásica se convierte en un operador cuántico. El tensor métrico para las coordenadas generalizadas (ángulos de Euler y distancia relativa) se calcula explícitamente.

2. Energía Potencial:

   • El potencial de interacción se modela como la suma de un término de Coulomb y un término nuclear.
   • La interacción nuclear se aproxima mediante un potencial de tipo Morse derivado de la plegadura (folding) de las interacciones nucleon-nucleón (fuerzas de Migdal) con las densidades de los fragmentos.
   • El potencial efectivo se analiza para identificar un "bolsillo" (mínimo local) donde los fragmentos están ligados en una configuración de contacto. La estabilidad de este estado molecular depende de la interacción entre la repulsión de Coulomb y la atracción nuclear, lo que conduce a un criterio de estabilidad ($Z^2/A^{4/3} \lesssim 6$) análogo al parámetro de fisilidad de la gota líquida.
   • El potencial se expande en una serie multipolar, lo que resulta en una dependencia angular aproximada por un potencial de doble pozo ($\sin^2 \epsilon$), donde $\epsilon$ es el ángulo entre el eje de simetría del fragmento deformado y el vector interfragmento.

3. Hamiltoniano y Aproximaciones:

   • Se emplea la aproximación de Born-Oppenheimer, separando el movimiento radial rápido (distancia relativa $R$) del movimiento angular lento. El movimiento radial se trata como un creador de un potencial efectivo para los grados de libertad angulares.
   • El Hamiltoniano total se diagonaliza en la base de funciones esféricas bipolares.
   • Se resuelven dos casos analíticos límite:
     • Configuración Polo a Polo: Ocurre cuando el fragmento esférico se encuentra cerca de los polos del fragmento deformado ($\epsilon \approx 0, \pi$). Este régimen se caracteriza por valores bajos de $K$ (proyección del momento angular sobre el eje de simetría) y se trata como un oscilador armónico en el modo de flexión (bending mode).
     • Configuración Ecuatorial: Ocurre cuando la proyección del momento angular $K$ excede un valor crítico ($K > K_{cr}$). En este régimen, el término centrífugo supera la barrera de potencial, creando un único mínimo en el ecuador ($\epsilon = \pi/2$). Esto conduce a una rotación inclinada y al movimiento de precesión (nutación) del sistema.

Contribuciones Clave y Resultados

• Derivación de un Hamiltoniano Consistente: El artículo proporciona una derivación rigurosa del Hamiltoniano para una molécula nuclear pesada, reteniendo explícitamente el acoplamiento entre los modos rotacionales y vibracionales (acoplamiento de Coriolis) que fue omitido o aproximado en trabajos previos.
• Soluciones Analíticas:
  • Para el límite polo a polo, los autores derivan espectros de energía que asemejan bandas de un rotador más bandas vibracionales armónicas, con correcciones para el acoplamiento entre rotación y vibración calculadas mediante la teoría de perturbaciones de segundo orden.
  • Para el límite ecuatorial ($K > K_{cr}$), los autores identifican una transición de fase donde el mínimo de potencial se desplaza al ecuador. Derivan expresiones analíticas para los niveles de energía, describiendo un sistema que experimenta rotación inclinada y nutación. Se muestra que la relación entre las frecuencias de rotación y nutación depende de la asimetría de masa del sistema.
• Validación Numérica: Los resultados analíticos se comparan con la diagonalización numérica completa del Hamiltoniano. La concordancia es notable, particularmente para la descripción de las excitaciones roto-vibracionales en $^{240}\text{Pu}$ (modelado como $^{236}\text{U} + \alpha$).
• Aplicación a Estados Hiperdeformados (HD):
  • El modelo se aplica a los estados HD de $^{232}\text{Th}$, tratado como un sistema molecular $^{132}\text{Sn} + ^{100}\text{Zr}$.
  • El espectro de excitación calculado y las constantes rotacionales muestran una buena concordancia con los datos experimentales para los estados HD en $^{232,234,236}\text{U}$.
  • Los autores proporcionan predicciones para los momentos multipolares ($Q_2, Q_3$) de estos estados moleculares.
• Anisotropía Angular de Fragmentos de Fisión:
  • El modelo se aplica a la distribución angular de los fragmentos de fisión en la fisión inducida por neutrones de $^{240}\text{Pu}$.
  • Al calcular la distribución de probabilidad de los valores de $K$ en el punto de escisión, los autores reproducen la anisotropía angular experimental a bajas energías de neutrones incidentes ($E_n < 4$ MeV).
  • El estudio sugiere que la anisotropía angular está fuertemente influenciada por la distribución de los valores de $K$ en el espectro de excitación de baja energía del núcleo fisible.

Significancia y Reivindicaciones
Los autores afirman que su modelo unificado resuelve con éxito las inconsistencias encontradas en tratamientos previos de moléculas nucleares pesadas, específicamente respecto al momento de inercia y al acoplamiento de los modos rotacionales y vibracionales. El enfoque se presenta como un marco consistente que:

1. Describe correctamente la dinámica colectiva de fragmentos pesados donde las aproximaciones previas fallaron.
2. Identifica distintos regímenes estructurales (polo a polo frente a ecuatorial) basados en la proyección del momento angular $K$.
3. Proporciona una base teórica para interpretar los estados hiperdeformados como configuraciones moleculares.
4. Ofrece un mecanismo para comprender la generación del momento angular y la distribución angular de los fragmentos de fisión a través de la dinámica del sistema dinuclear en la escisión.

El artículo enfatiza que, si bien el modelo es consistente con los datos existentes para casos específicos (como $^{240}\text{Pu}$ y estados HD de actínidos), su principal contribución es el desarrollo de una herramienta teórica autoconsistente que puede aplicarse a diversos fenómenos moleculares nucleares sin depender de aproximaciones ad hoc.