Nuclear molecule of heavy nuclei

This paper proposes a nuclear molecule model for heavy nuclei composed of two interacting heavy nuclei, deriving a Hamiltonian to analytically and numerically describe roto-vibrational excitations in $^{240}$Pu and predict hyperdeformed states in $^{232}$Th, while also analyzing fission fragment angular distributions.

The Gist

Imagine the nucleus of an atom not as a single, solid ball of dough, but as a cosmic dance partnership between two heavy partners. This is the core idea of the paper you shared: the concept of a "nuclear molecule."

Here is a simple breakdown of what the authors, T. M. Shneidman and R. G. Nazmitdinov, are proposing, using everyday analogies.

1. The Big Idea: Two Nuclei Holding Hands

Usually, we think of an atomic nucleus as one big lump. But the authors suggest that under certain conditions, a heavy nucleus can split into two distinct parts that stay stuck together, like two people holding hands.

• The Partners: One partner is a perfect sphere (like a billiard ball), and the other is a squashed, egg-shaped ball (like a rugby ball).
• The Glue: They are held together by the "nuclear force," which acts like a very strong, sticky glue.
• The Tension: At the same time, they are pushing away from each other because they both have a positive electrical charge (Coulomb repulsion), like trying to push two north poles of magnets together.

The paper argues that when these two forces balance out, they form a stable "molecule" that can vibrate and spin.

2. How They Move: The Dance Floor

The authors created a mathematical model (a Hamiltonian) to describe how this "dance" works. They looked at two main ways these partners can move:

• The Pole-to-Pole Dance (The "Top" Position):
  Imagine the spherical partner sitting right on the "North Pole" or "South Pole" of the egg-shaped partner.

  • The Motion: The sphere can wiggle back and forth around the pole, like a child spinning a top that is slightly off-center. It can also vibrate up and down.
  • The Result: This creates specific energy levels (notes) that the molecule can sing. The authors found that if the egg-shaped partner is very squashed, the sphere gets "stuck" near the pole and can't easily jump to the other side.

• The Equatorial Dance (The "Waist" Position):
  Now, imagine the spherical partner moving to the "waist" or equator of the egg-shaped partner.

  • The Motion: This happens when the system spins very fast. The sphere starts orbiting around the waist of the egg.
  • The Wobble: As it orbits, the whole system starts to wobble or "nutate" (like a spinning top that is tilting and wobbling). The authors compare this to a specific type of instability in physics called an "Andronov-Hopf bifurcation"—basically, a smooth circle turning into a wobbly, precessing motion.

3. The "Phase Transition"

One of the paper's cool discoveries is that the dance changes depending on how fast the system spins.

• Slow Spin: The partners stay at the poles (the "Pole-to-Pole" mode).
• Fast Spin: Once the spin gets fast enough (reaching a "critical speed"), the partners suddenly switch. The sphere slides down to the waist and starts orbiting there (the "Equatorial" mode).
• The Analogy: Think of a spinning coin. When it spins slowly, it stands up. When it spins fast enough, it flattens out and spins on its edge. The nucleus does something similar with its shape.

4. Testing the Theory

The authors didn't just make up the math; they tested it against real-world data.

• Case Study 1: The "Hyperdeformed" Nucleus (232Th):
  They looked at a heavy nucleus called Thorium-232. They suggested that its most stretched-out, excited states look exactly like a molecule made of a Tin-132 nucleus and a Zirconium-100 nucleus.

  • The Result: Their mathematical predictions for the energy levels of this "molecule" matched the experimental data very well.

• Case Study 2: Fission (240Pu):
  They looked at Plutonium-240 just before it splits (fissions). They treated the moment right before the split as a nuclear molecule.

  • The Prediction: They calculated how the pieces would fly apart (the angular distribution).
  • The Result: Their model predicted the angles at which the fragments fly out, and it matched the experimental data very closely, especially at lower energies.

5. Why This Matters

The authors point out that previous models often ignored the complex interaction between the spinning of the whole molecule and the wobbling of the individual parts. By fixing this math, they get a more accurate picture of how heavy nuclei behave.

In summary: This paper proposes that heavy atomic nuclei can act like dancing couples. Depending on how fast they spin, they can either hold hands at the poles or orbit around the waist. The authors built a new set of rules to describe this dance and proved that these rules accurately predict how real nuclei (like Thorium and Plutonium) behave in experiments.

Technical Summary

Technical Summary: Nuclear Molecule of Heavy Nuclei

Problem Statement
The paper addresses the theoretical description of heavy nuclei as molecular systems composed of two interacting fragments (a dinuclear system, or DNS). While the concept of nuclear molecules has been successfully applied to light fragments (e.g., $\alpha$-clusters) and alpha decay, previous attempts to formulate a Hamiltonian for heavy fragments often relied on inconsistent approximations. Specifically, earlier models (Refs. [9, 10]) replaced the total kinetic energy Hamiltonian, including Coriolis terms, with an asymmetric rotator Hamiltonian. The authors argue that this approximation leads to an underestimated rotational moment of inertia and fails to correctly describe the evolution of heavy fragments, particularly when one fragment is strongly deformed. The goal is to derive a self-consistent Hamiltonian that accurately describes the collective dynamics (rotational and vibrational) of a di-nuclear system consisting of a heavy, axially-symmetric deformed fragment and a spherical fragment.

Methodology
The authors develop a quantum mechanical model based on the following steps:

1. Coordinate Systems and Kinetic Energy:

   • The system is defined by a laboratory frame and intrinsic frames for the two fragments. The relative distance vector $\mathbf{R}$ connects the centers of mass, and the orientation of the deformed fragment is defined by Euler angles.
   • A classical kinetic energy expression is derived by separating the motion of the center of mass, the relative motion of the fragments, and the internal rotation of the deformed fragment.
   • Using the Pauli quantization procedure, the classical kinetic energy is converted into a quantum operator. The metric tensor for the generalized coordinates (Euler angles and relative distance) is explicitly calculated.

2. Potential Energy:

   • The interaction potential is modeled as the sum of a Coulomb term and a nuclear term.
   • The nuclear interaction is approximated using a Morse-type potential derived from folding nucleon-nucleon interactions (Migdal forces) with fragment densities.
   • The effective potential is analyzed to identify a "pocket" (local minimum) where the fragments are bound in a touching configuration. The stability of this molecular state is shown to depend on the interplay between Coulomb repulsion and nuclear attraction, leading to a stability criterion ($Z^2/A^{4/3} \lesssim 6$) analogous to the liquid drop fissility parameter.
   • The potential is expanded in a multipole series, resulting in an angular dependence approximated by a double-well potential ($\sin^2 \epsilon$), where $\epsilon$ is the angle between the symmetry axis of the deformed fragment and the inter-fragment vector.

3. Hamiltonian and Approximations:

   • The Born-Oppenheimer approximation is employed, separating the fast radial motion (relative distance $R$) from the slow angular motion. The radial motion is treated as creating an effective potential for the angular degrees of freedom.
   • The total Hamiltonian is diagonalized in the basis of bipolar spherical functions.
   • Two limiting analytical cases are solved:
     • Pole-to-Pole Configuration: Occurs when the spherical fragment is near the poles of the deformed fragment ($\epsilon \approx 0, \pi$). This regime is characterized by low $K$ values (projection of angular momentum on the symmetry axis) and is treated as a harmonic oscillator in the bending mode.
     • Equatorial Configuration: Occurs when the angular momentum projection $K$ exceeds a critical value ($K > K_{cr}$). In this regime, the centrifugal term overcomes the potential barrier, creating a single minimum at the equator ($\epsilon = \pi/2$). This leads to a tilted rotation and precessional motion (nutation) of the system.

Key Contributions and Results

• Derivation of a Consistent Hamiltonian: The paper provides a rigorous derivation of the Hamiltonian for a heavy nuclear molecule, explicitly retaining the coupling between rotational and vibrational modes (Coriolis coupling) which was neglected or approximated in previous works.
• Analytical Solutions:
  • For the pole-to-pole limit, the authors derive energy spectra resembling a rotor plus harmonic vibrational bands, with corrections for the coupling between rotation and vibration calculated via second-order perturbation theory.
  • For the equatorial limit ($K > K_{cr}$), the authors identify a phase transition where the potential minimum shifts to the equator. They derive analytical expressions for the energy levels, describing a system undergoing tilted rotation and nutation. The ratio of rotational to nutation frequencies is shown to depend on the mass asymmetry of the system.
• Numerical Validation: The analytical results are compared with full numerical diagonalization of the Hamiltonian. The agreement is found to be remarkable, particularly for the description of roto-vibrational excitations in $^{240}\text{Pu}$ (modeled as $^{236}\text{U} + \alpha$).
• Application to Hyperdeformed (HD) States:
  • The model is applied to the HD states of $^{232}\text{Th}$, treated as a $^{132}\text{Sn} + ^{100}\text{Zr}$ molecular system.
  • The calculated excitation spectrum and rotational constants show good agreement with experimental data for HD states in $^{232,234,236}\text{U}$.
  • The authors provide predictions for the multipole moments ($Q_2, Q_3$) of these molecular states.
• Fission Fragment Angular Anisotropy:
  • The model is applied to the angular distribution of fission fragments in neutron-induced fission of $^{240}\text{Pu}$.
  • By calculating the probability distribution of $K$-values at the scission point, the authors reproduce the experimental angular anisotropy at low incident neutron energies ($E_n < 4$ MeV).
  • The study suggests that the angular anisotropy is strongly influenced by the distribution of $K$-values in the low-energy excitation spectrum of the fissile nucleus.

Significance and Claims
The authors claim that their unified model successfully resolves inconsistencies found in previous treatments of heavy nuclear molecules, specifically regarding the moment of inertia and the coupling of rotational and vibrational modes. The approach is presented as a consistent framework that:

1. Correctly describes the collective dynamics of heavy fragments where previous approximations failed.
2. Identifies distinct structural regimes (pole-to-pole vs. equatorial) based on the angular momentum projection $K$.
3. Provides a theoretical basis for interpreting hyperdeformed states as molecular configurations.
4. Offers a mechanism for understanding the angular momentum generation and angular distribution of fission fragments through the dynamics of the dinuclear system at scission.

The paper emphasizes that while the model is consistent with existing data for specific cases (like $^{240}\text{Pu}$ and actinide HD states), its primary contribution is the development of a self-consistent theoretical tool that can be applied to various nuclear molecular phenomena without relying on ad-hoc approximations.