Algebraic Nilsson cranking model and its prediction for 20Ne

This paper presents an algebraic solution to the self-consistent Nilsson cranking model for the 20Ne nucleus, which yields ground-state rotational-band energies in significantly better agreement with experimental data than previous numerical methods and offers insights into rotation type transitions and weak pairing correlations.

The Gist

The Big Picture: Spinning a Nucleus

Imagine an atomic nucleus (like the one in Neon-20) not as a tiny, static ball, but as a spinning, squishy water balloon.

Physicists want to understand how this balloon spins. Does it spin like a perfect top (one axis)? Does it wobble like a spinning coin (two axes)? Or does it do something weird and complex (three axes)?

For decades, scientists used a standard tool called the "Cranking Model" to predict this. Think of this model like a manual car. You (the physicist) have to guess how hard to press the gas pedal (the "angular velocity") to make the car go a certain speed. If you guess wrong, the prediction doesn't match reality. It's a bit of a "best guess" approach.

The Problem with the Old Way

In this paper, the authors look at a specific nucleus: Neon-20.

• The Old Prediction: Using the "manual car" method (a numerical solution), previous scientists predicted the energy levels of Neon-20. They were okay, but not perfect. They missed some details, especially at higher spin speeds (when the nucleus is spinning very fast).
• The Mystery: When you look at real data, the energy required to spin the nucleus at certain speeds (specifically when the spin is 4 or 8) is lower than the old models predicted. It's like the nucleus suddenly becomes "easier" to spin at those specific moments.

The New Solution: The "Self-Driving" Car

The authors (Gulshani and Lahbas) developed a new way to solve the equations. Instead of guessing the gas pedal (angular velocity), they built a self-driving car.

In their new model, the "speed" of the spin isn't a guess; it is calculated automatically based on how the nucleus is shaped and how the particles inside are moving. It's a "microscopic" and "self-consistent" approach. This means the model adjusts itself until the shape of the nucleus and the speed of its spin perfectly match each other.

They used a special algebraic method (a set of mathematical shortcuts) to solve this, rather than just brute-forcing the numbers with a computer.

The "Magic" of Iteration: The Bouncing Ball

Here is the most interesting part of their discovery, explained with an analogy:

Imagine you are trying to balance a ball on a hill.

1. For some spins (I=2 and I=6): The ball rolls down and settles in a nice, quiet valley. The math is stable. The prediction is smooth.
2. For other spins (I=4 and I=8): The ball doesn't settle. It starts bouncing back and forth between two different valleys.
   • In the math, this looks like an oscillation. The model keeps switching between two slightly different shapes for the nucleus.
   • Why? It's because the internal particles (the "nucleons") are crossing paths, like cars changing lanes on a highway. When they switch lanes, the energy of the system jumps up and down.

The authors realized that this "bouncing" isn't a mistake; it's a feature!

• At Spin 4, the energy bounces between two values. By picking the lower of the two, they found the true energy, which matched the real-world measurements perfectly.
• At Spin 8, the bouncing eventually settles down into a new, lower energy state. This explains why the nucleus is "easier" to spin at this speed in real life.

The "Quenching" Effect

The paper mentions a transition from "planar rotation" to "uniaxial rotation."

• Planar Rotation: Imagine a Frisbee spinning flat on the ground. It's a wide, flat spin.
• Uniaxial Rotation: Imagine a spinning top standing straight up.

The authors suggest that at high speeds (Spin 8), the nucleus decides to stop spinning like a flat Frisbee and starts spinning like a top. This change in shape (from flat to tall) releases energy, making the nucleus more stable. This is what they call "quenching" (damping) the planar rotation.

Why Does This Matter?

1. Better Accuracy: Their new "algebraic" method predicted the energy levels of Neon-20 much better than the old "numerical" method. It matched the real-world data almost perfectly.
2. No Pairing Needed: In heavier nuclei, particles like to pair up (like dance partners), which complicates the math. But in Neon-20, these pairs don't really matter. The authors proved that you can get perfect results without worrying about these pairs, simplifying the physics.
3. Solving the Mystery: They explained why the energy drops at Spin 4 and 8. It's because the nucleus is constantly reshaping itself (oscillating) to find the most comfortable, lowest-energy position.

Summary

Think of the nucleus as a dancer.

• Old Model: The choreographer guessed the steps. Sometimes the dancer looked a bit stiff or off-rhythm.
• New Model: The dancer figures out their own steps based on their own body mechanics.
• The Result: The dancer (Neon-20) performs a complex routine where they wobble and change shape (oscillate) at specific beats (Spin 4 and 8). This wobble actually saves energy, making the performance smoother and more efficient than anyone expected.

The authors have shown that by letting the math "self-correct" and watching how the nucleus wobbles, we can understand the secrets of how atoms spin.

Technical Summary

1. Problem Statement

The paper addresses limitations in the Conventional Cranking Model (CCRM3) used to describe triaxial collective rotation in deformed nuclei.

• Semi-classical Nature: The standard CCRM3 relies on a constant, adjustable angular velocity parameter ($\vec{\Omega}$), making it semi-classical and phenomenological rather than fully microscopic.
• Symmetry Breaking: The model breaks time-reversal and $D_2$ symmetries, leading to conceptual issues when comparing predicted spectra with experimental data (specifically regarding the superposition of even and odd angular momentum eigenstates).
• Previous Discrepancies: Previous numerical solutions of the Nilsson-CCRM3 model (Nilsson-Ragnarsson) for the $^{20}\text{Ne}$ nucleus predicted ground-state rotational band energies ($I=2, 4, 6, 8$) that did not fully align with measured experimental values, particularly failing to explain the lower-than-expected excitation energies at $I=4$ and $I=8$.
• Goal: The authors aim to solve the self-consistent Nilsson-CCRM3 Schrödinger equation using a new algebraic method to improve agreement with experimental data and provide a physical explanation for the energy anomalies in $^{20}\text{Ne}$.

2. Methodology

The authors developed an algebraic iterative solution for the Nilsson-CCRM3 equations, incorporating spin-orbit interaction but neglecting pairing correlations (justified as insignificant in the light $s-d$ shell nucleus $^{20}\text{Ne}$).

• Hamiltonian Formulation:

  • The single-particle Hamiltonian includes a deformed harmonic oscillator potential and a spin-orbit term ($l \cdot s$).
  • The cranking term ($-\vec{\Omega} \cdot \vec{j}$) is included, where $\vec{j} = \vec{l} + \vec{s}$.
  • The spin-orbit term is treated by re-expressing the Hamiltonian to combine the cranking angular velocity and the spin expectation value into an effective angular velocity vector.

• Algebraic Solution Technique:

  • Instead of numerical diagonalization, the authors utilize Heisenberg equations of motion for the particle-coordinate radius vector.
  • They assume a periodic motion solution with frequency $\alpha$, leading to a cubic characteristic equation for the normal mode frequencies.
  • The eigenstates are expressed in terms of decoupled harmonic oscillator phonon quantum numbers ($n_1, n_2, n_3$).

• Self-Consistency and Constraints:

  • Constant Volume: The product of oscillator frequencies is kept constant ($\omega_1 \omega_2 \omega_3 = \omega_0^3$).
  • Shape Minimization: The intrinsic energy is minimized with respect to deformation parameters ($\epsilon, \gamma$) and the angular velocity vector direction.
  • Iterative Process: The model iterates on deformation parameters until convergence is achieved. For specific angular momenta ($I$), the magnitude of the cranking frequency is adjusted until the calculated expectation value of angular momentum matches the target $I\hbar$.

3. Key Contributions

• Algebraic vs. Numerical: The primary contribution is the derivation of an algebraic method to solve the self-consistent Nilsson-CCRM3 equations, contrasting with previous numerical approaches.
• Microscopic Insight into Oscillations: The study reveals that the iterative process for certain angular momenta ($I=4, 8$) does not converge to a single static value but exhibits periodic/cyclic oscillations in energy and deformation parameters.
• Level Crossing Mechanism: The authors identify that these oscillations are caused by the crossing of single-particle energy levels during the iteration. The system alternates between two coupled states (lower and higher energy) due to self-consistency conditions.
• Decoupling Strategy: A novel approach is proposed where the "yrast" (lowest energy) state is isolated by selecting the deformation parameters corresponding to the lower-energy component of the oscillating system, effectively decoupling the states.

4. Results

The model was applied to the $^{20}\text{Ne}$ nucleus for angular momenta $I = 2, 4, 6, 8$.

• Improved Agreement with Experiment:

  • The algebraic model predicts excitation energies that show much better agreement with measured experimental values compared to the previous numerical solutions by Nilsson-Ragnarsson.
  • Specifically, the model successfully reproduces the reduction in excitation energy at $I=4$ and $I=8$ relative to $I=2$ and $I=6$, a feature previously difficult to explain.

• Behavior by Angular Momentum:

  • $I=2$ and $I=6$: The energy and deformation parameters converge steadily to stable values.
  • $I=4$: The energy exhibits persistent oscillations between two values (approx. 51 and 52.5 in units of $\hbar\omega_0$) corresponding to two different deformation sets. By isolating the lower-energy branch, the predicted energy is lowered, matching the experimental trend.
  • $I=8$: The energy and deformation undergo cyclic variations that eventually decay to a steady, lower-energy state. This is attributed to the crossing of different single-particle levels.

• Rotation Types and Deformations:

  • $I=2$: Triaxial prolate shape, nearly planar rotation in the $x-y$ plane.
  • $I=4$: Highly prolate, nearly axially symmetric about the $z$-axis, with triaxial rotation mostly about the $x$-axis.
  • $I=6$: Triaxial highly prolate shape, triaxial rotation about the $x$-axis.
  • $I=8$: Almost oblate axially symmetric about the $x$-axis, with triaxial rotation.

5. Significance and Conclusions

• Explanation of Energy Quenching: The paper provides a mechanism for the observed "quenching" or reduction of yrast energies at $I=4$ and $I=8$. The cyclic variations and level crossings inherent in the self-consistent algebraic solution naturally lead to lower energy states, supporting suggestions by Bohr and Mottelson regarding the impact of oscillations on nuclear spectra.
• Validation of Algebraic Approach: The results validate the algebraic method as a powerful tool for solving cranking models, offering physical insights (such as the nature of level crossings) that might be obscured in purely numerical black-box solutions.
• Implications for Pairing: The success of the model without pairing correlations reinforces the understanding that pairing is weak in light $s-d$ shell nuclei like $^{20}\text{Ne}$.
• Future Directions: The authors plan to integrate this algebraic Nilsson-CCRM3 framework with their previously developed Microscopic Self-Consistent Cranking Model (MSCRM3), where the angular velocity is derived microscopically rather than being an adjustable parameter, to further refine predictions for rotational bands.

In summary, this work demonstrates that an algebraic treatment of the self-consistent Nilsson cranking model can resolve discrepancies in predicted rotational spectra for $^{20}\text{Ne}$ by explicitly accounting for the dynamic interplay between single-particle level crossings and nuclear deformation.