Nilsson Diagram

While the spherical shell model brilliantly explains the stability of certain "magic" nuclei, the vast majority of atomic nuclei are not spherical. They are deformed, stretched or flattened in a way that the simple model cannot account for. This introduces a fundamental challenge: how can we understand the quantum states and properties of nucleons moving within a non-spherical potential? The Nilsson model provides a powerful and elegant answer to this question, serving as a cornerstone of modern nuclear structure physics. This article delves into this essential framework. In the first part, "Principles and Mechanisms", we will explore how the model is constructed, starting from a deformed harmonic oscillator and incorporating key physical ingredients like the spin-orbit force to generate the famous Nilsson diagram. Subsequently, in "Applications and Interdisciplinary Connections", we will witness the model's remarkable predictive power, seeing how it explains everything from the static shapes and moments of nuclei to their dynamic behavior in reactions, rotations, and even the cataclysmic process of fission.

Imagine you are trying to build a model of the solar system. A first good guess is that planets are spheres and their orbits are circles. This is simple, elegant, and captures the main idea. But soon you discover that planets are not perfect spheres, and their orbits are ellipses. To truly understand the dynamics, you need to account for these "deformations." The world of the atomic nucleus is much the same. The simple shell model, which pictures nucleons orbiting in a perfect spherical potential, is a brilliant starting point. It successfully explains the "magic numbers" where nuclei are exceptionally stable. But most nuclei are not spherical. They are stretched like a cigar (prolate) or squashed like a pancake (oblate). To understand this vast majority of nuclei, we need a model that embraces their deformed nature. This is the world of the Nilsson model.

How do we describe a deformed shape? We can start with a sphere of radius $R_0$ and stretch it along one axis (let's call it the $z$-axis) while shrinking it in the perpendicular directions to keep the total volume constant—a reasonable assumption for the nearly incompressible nuclear matter. We can define a deformation parameter, say $\delta$, to quantify this. For a prolate shape, the length along the symmetry axis might be $R_z = R_0 (1 + \frac{2}{3}\delta)$ and the radius in the perpendicular plane would be $R_x = R_y = R_0 (1 - \frac{1}{3}\delta)$. Notice how a positive $\delta$ gives you a cigar shape. This is a common way to parameterize the nuclear shape. There are other ways, too, like using spherical harmonics, which connect to the collective motion of the nucleus as a liquid drop. For small deformations, these different descriptions are directly related, providing a consistent language to talk about nuclear shapes.

Now, the crucial step taken by Sven Gösta Nilsson was to assume that the potential well that a single nucleon feels has the same shape as the nucleus itself. If the nucleus is a prolate spheroid, the potential is also a prolate spheroidal well. The simplest way to model this is with an ​​anisotropic harmonic oscillator​​. Instead of one spring constant (or frequency, $\omega_0$), we have two: one for motion along the symmetry axis ($\omega_z$) and another for motion perpendicular to it ($\omega_\perp$).

A fundamental assumption connects the shape of the nucleus to the shape of the potential. If we define the geometric shape by a parameter $\delta$ and the potential's anisotropy by a parameter $\epsilon$, which relates $\omega_\perp$ and $\omega_z$, we can derive a direct relationship between them. This elegant link ensures that our model of a single particle moving in a potential is consistent with the collective shape of the entire nucleus. For a prolate shape ($\delta > 0$), the potential well is shallower and wider along the $z$-axis, meaning $\omega_z \omega_\perp$. A nucleon can travel further along this axis before being pulled back.

In the perfectly spherical shell model, energy levels are highly degenerate. For instance, in a $1f_{7/2}$ shell, a nucleon can have a total angular momentum of $j=7/2$. Its projection on any axis, $m_j$, can take any of the $2j+1=8$ values from $-7/2$ to $+7/2$, and all these states have the exact same energy.

But the moment we deform the nucleus, this beautiful symmetry is broken. The $z$-axis is now special—it's the symmetry axis of our spheroid. The total angular momentum $j$ is no longer a conserved quantity because a torque is exerted on the orbiting nucleon. However, the projection of the total angular momentum onto this special axis, which we call $\Omega$, is conserved.

What happens to our degenerate energy levels? They split. States with different absolute values of $\Omega$ now have different energies. Imagine our nucleon's orbit: if it orbits mostly in the plane perpendicular to the symmetry axis (high $|m_j|$), it feels a different potential than if it orbits mostly along the symmetry axis (low $|m_j|$). For a prolate (cigar-shaped) nucleus, orbits along the long axis are energetically favored. This leads to a splitting of the original $j$-shell. For example, a simple deformation potential proportional to $3\hat{j}_z^2 - \hat{j}^2$ shows that the states with the smallest $|\Omega|$ (orbits along the long axis) decrease in energy, while states with the largest $|\Omega|$ (orbits in the short plane) increase in energy. The original single energy level of the $1f_{7/2}$ shell fractures into four distinct levels, each corresponding to a pair of states $\pm\Omega$. This splitting is the central feature of the Nilsson diagram.

A simple anisotropic oscillator is a good start, but it's not realistic enough. The true nuclear potential is not a perfect parabola; its bottom is flatter, more like a wine bottle. Furthermore, there is a powerful force inside the nucleus that couples a nucleon's orbital motion to its intrinsic spin—the ​​spin-orbit interaction​​.

The full ​​Nilsson Hamiltonian​​ is a masterful blend of these ingredients: $H_{Nilsson} = H_{osc} + C \, \mathbf{l} \cdot \mathbf{s} + D \, \mathbf{l}^2$

Let's break this down.

1. $H_{osc}$ is our anisotropic harmonic oscillator, which accounts for the deformed shape.
2. $D \, \mathbf{l}^2$ is a correction term. With $D0$, it lowers the energy of states with higher orbital angular momentum $l$, making the potential well more like that realistic, flat-bottomed shape.
3. $C \, \mathbf{l} \cdot \mathbf{s}$ is the crucial spin-orbit term. With $C0$, it dramatically lowers the energy of states where the spin and orbital angular momentum are aligned ($j = l + 1/2$) and raises the energy of states where they are anti-aligned ($j = l - 1/2$).

It is this combination of the $\mathbf{l} \cdot \mathbf{s}$ and $\mathbf{l}^2$ terms that correctly reproduces the magic numbers of the spherical shell model before we even introduce deformation. Interestingly, for specific ratios of the strengths $C$ and $D$, certain pairs of orbitals can become degenerate, giving rise to "pseudo-spin symmetries" that are a fascinating feature of modern nuclear theory.

The full magic of the Nilsson model comes from the competition between these terms. The deformation part of $H_{osc}$ wants to sort states by their spatial orientation ($\Omega$). The spin-orbit term, $\mathbf{l} \cdot \mathbf{s}$, wants to sort states by the coupling of spin and orbit ($j$). As we "dial up" the deformation, we see a complex interplay where states of the same $\Omega$ but different $j$ (and $l$) can mix. A state is no longer a pure $|l, j\rangle$ configuration but becomes a superposition of several basis states, all sharing the same conserved quantum number $\Omega$.

The Nilsson diagram is the graphical summary of this entire story. It's a plot with deformation on the horizontal axis and single-particle energy on the vertical axis. Each line on the diagram represents the energy of a single Nilsson orbital as the nucleus deforms.

At zero deformation on the far left, you see the familiar, highly degenerate energy levels of the spherical shell model. As you move to the right (increasing prolate deformation), you see these levels split dramatically, fanning out into a beautiful and complex web of crisscrossing lines.

For large deformations, the harmonic oscillator part of the potential dominates, and the states can be labeled by a new set of "asymptotic" quantum numbers: $[N, n_z, \Lambda]\Omega$.

• $N$ is the total oscillator quantum number, just like in the spherical case.
• $n_z$ is the number of energy quanta along the symmetry axis. It tells you how much the nucleon's motion is concentrated along the long axis.
• $\Lambda$ is the projection of the orbital angular momentum on the symmetry axis.
• $\Omega$, our conserved quantity, is the projection of the total angular momentum. The spin projection $\Sigma$ is simply $\Omega - \Lambda$.

These labels give us a deep physical intuition for each orbital. For example, a state with a low $n_z$ is a "down-sloping" orbital; its energy decreases with prolate deformation because its motion avoids the "squeezed" perpendicular directions. The diagram is filled with ​​level crossings​​, points where two orbitals meet. These crossings are profoundly important. If the ground state of a nucleus corresponds to filling orbitals up to a certain level, and that level crosses another as deformation changes, the nucleus might suddenly find it energetically favorable to snap into a different shape.

So far, we've only considered axially symmetric shapes—cigars and pancakes. But what if the nucleus is a bit more like a potato, with three different axis lengths? This is called ​​triaxiality​​, described by a parameter $\gamma$.

Introducing triaxiality breaks the axial symmetry. Now, not even $\Omega$ is a conserved quantity! The Hamiltonian becomes more complex. What happens to our beautiful level crossings? If the two crossing orbitals have the correct quantum numbers to interact, they no longer cross. As they approach each other, a residual interaction, perhaps induced by the triaxiality itself, causes them to repel. This is a classic quantum mechanical phenomenon known as an ​​avoided crossing​​ or ​​level repulsion​​. The energy splitting at the point of closest approach depends on the strength of this interaction. Instead of crossing, the wave functions of the two states mix strongly and the orbitals exchange their character. This phenomenon is critical for understanding the fine details of nuclear spectra and the transition between different nuclear shapes.

The principles behind the Nilsson diagram are a testament to the power of quantum mechanics in describing complex systems. Starting from simple ideas of shape and symmetry, and adding the essential physical ingredients of spin-orbit coupling, we arrive at a model of breathtaking predictive power. The diagram is not just a collection of lines; it is a roadmap of the nucleus, guiding our understanding of its structure, stability, and the rich dynamics hidden within the heart of the atom.

In our previous discussion, we dismantled the machinery of the Nilsson model, laying bare its gears and cogs—the deformed potential, the spin-orbit force, and the resulting energy levels. We saw how this elegant model organizes the complex world of single-particle states within a non-spherical nucleus. But a model, no matter how elegant, is only as good as its ability to describe reality. Now, we embark on a more exciting journey: to see how the Nilsson diagram transforms from a theoretical map into a powerful predictive tool, a veritable Rosetta Stone for deciphering the language of the nucleus. We will find that this single framework brings a remarkable unity to a vast landscape of nuclear phenomena, from the static shape of a nucleus to the violent dynamics of its fission.

The most immediate success of the Nilsson model is its justification for the very thing it assumes: that many nuclei are not spherical. By examining the energy levels as a function of deformation $\delta$, we see that for certain numbers of protons or neutrons, the total energy of the system is minimized when the nucleus is shaped like a football (prolate) or a doorknob (oblate). But how do we confirm this? We can measure a nucleus’s electric quadrupole moment, $Q_0$, which is the definitive measure of its deviation from a sphere. The Nilsson model provides a direct path to calculating it. The total quadrupole moment is simply the sum of the contributions from each individual nucleon, and the contribution of a nucleon in a specific Nilsson orbital $|k\rangle$ is given by the expectation value $q_k = \langle k|\hat{q}_0|k\rangle$.

However, the reality is a bit more subtle and beautiful. Nucleons are not just independent particles; they feel a residual "pairing" force, which encourages them to form spin-zero pairs, much like electrons in a superconductor. The Bardeen-Cooper-Schrieffer (BCS) theory accounts for this by telling us that Nilsson orbitals are not simply "full" or "empty." Instead, each level has a certain occupation probability, $v_k^2$. A more refined calculation of the total quadrupole moment therefore involves summing over all valence orbitals, with each single-particle moment $q_k$ weighted by its occupation probability $v_k^2$. The agreement between these calculated values and experimental measurements provides stunning confirmation that this combined picture—of independent particles moving in a deformed field, modified by pairing correlations—is fundamentally correct.

Beyond shape, the model makes sharp predictions about a nucleus's spin and parity, $J^\pi$. For an odd-A nucleus, these properties are determined by the last, unpaired nucleon. By simply identifying which Nilsson orbital this nucleon occupies on the diagram, we can read off its spin projection $\Omega$ and its parity $\pi=(-1)^N$, which, for many nuclei, correspond directly to the ground-state spin and parity of the entire nucleus.

A far more sensitive test, however, is the nuclear magnetic dipole moment, $\mu$. This property depends not just on the nucleon's total angular momentum, but on the delicate balance of its orbital and intrinsic spin contributions, which the Nilsson wavefunction explicitly details. For a rotational band built on an orbital with projection $K$, the magnetic moment reveals the interplay between the intrinsic magnetism of the odd nucleon (with g-factor $g_K$) and the collective rotation of the nuclear core ($g_R$). A particularly fascinating case arises for bands with $K=1/2$. Here, the Coriolis force—the same "fictitious" force you feel on a merry-go-round—becomes so strong that it significantly perturbs the motion, introducing a "decoupling parameter" $a$ into the formula for the magnetic moment. The success of this modified formula is a beautiful testament to the model's ability to incorporate the subtle dynamics of particle-core coupling.

The model's power is not limited to odd-A systems. For odd-odd nuclei, where we have both an unpaired proton and an unpaired neutron, the situation might seem hopelessly complex. Yet, the Gallagher-Moszkowski coupling rules bring order to the chaos. They tell us how the two individual projections, $\Omega_p$ and $\Omega_n$, will couple to form the total $K$ of the ground state. With this, and the Nilsson wavefunctions for the two odd nucleons, we can predict the magnetic moment of the entire odd-odd nucleus, as demonstrated for a nucleus like $^{176}\text{Lu}$.

Having characterized the nucleus at rest, we now turn to its dynamics—how it transforms and interacts. Here, the Nilsson model shines as a bridge between nuclear structure and nuclear reactions.

Consider beta decay, the process by which a neutron transforms into a proton (or vice versa), emitting an electron and a neutrino. The rate of this decay depends critically on the overlap between the initial nucleon's wavefunction and the final nucleon's wavefunction. If the parent and daughter nuclei are deformed, then these wavefunctions are Nilsson states. The model gives us these states not as simple entities, but as precise superpositions of spherical shell model states. This detailed knowledge allows us to directly calculate the Gamow-Teller matrix element, a key factor governing the decay probability. The Nilsson model thus provides the essential structural input needed to understand the dynamics of the weak nuclear force.

We can also probe nuclear structure by instigating reactions in a particle accelerator. In a direct reaction, like the deuteron stripping reaction $(d,p)$, a deuteron grazes a target nucleus, dropping off its neutron, which is then captured into a specific orbital. If the final nucleus is deformed, the neutron will land in a Nilsson orbital, and the states populated will be the members of the rotational band built upon that orbital. The theory of these reactions predicts that the probability (or cross-section) of populating a final state with spin $I$ is directly proportional to a spectroscopic factor, $S_I$. In the context of the Nilsson model, these spectroscopic factors are related to the coefficients $C_{jl}$ that define the expansion of the Nilsson orbital in a spherical basis. This provides an extraordinary tool: by measuring the relative intensities of the outgoing protons corresponding to different final states, we create a "fingerprint" of the Nilsson orbital and directly test the predicted composition of its wavefunction.

The Nilsson model is not just a descriptor of single-particle behavior; it is a cornerstone for understanding the collective motion of the entire nucleus.

A deformed nucleus is not a rigid object. It can vibrate and rotate. The stiffness against vibrations in its shape—so-called beta-vibrations—is determined by the curvature of its potential energy surface. This surface is not the smooth bowl predicted by the classical Liquid Drop Model. Instead, the quantum mechanics of the Nilsson levels introduces wiggles, or "shell corrections." The total energy is a sum of the macroscopic liquid-drop energy and the microscopic shell-correction energy. For a deformed nucleus, this sum results in an energy minimum at a non-zero deformation $\delta_0$. The stiffness against beta-vibrations, $C_{\delta\delta}$, is simply the second derivative of this total energy curve at its minimum. Thus, the collective vibrational properties of the nucleus are dictated by the underlying single-particle level structure.

Even more dramatic is what happens when we spin a deformed nucleus to extreme rotational frequencies. Initially, the nucleus spins faster as a whole, like a rigid rotating body. But as the Coriolis forces grow, a new, more energy-efficient mechanism becomes available. It can become energetically cheaper to break a pair of high-j nucleons and align their individual angular momenta with the axis of rotation. This sudden change in strategy, from collective to single-particle alignment, causes a dramatic kink, or "backbending," in the plot of the nucleus's moment of inertia versus rotational frequency. Using the cranked Nilsson model, we can model this by comparing the energy in the rotating frame (the Routhian) of the ground-state band with that of the aligned "super-band." The crossing of these two bands predicts the critical frequency $\omega_c$ at which backbending will occur.

The Coriolis force is also responsible for another subtle effect: it can mix states from different rotational bands. This explains why some electromagnetic transitions that should be strictly "forbidden" by the selection rule on the quantum number $K$ are nevertheless observed, albeit weakly. For instance, an E2 transition that requires $|\Delta K|=3$ is normally forbidden. However, if the initial state is Coriolis-mixed with another state for which the transition is allowed ($|\Delta K| \le 2$), the transition can proceed by "borrowing" strength from the allowed pathway. The observed transition rate becomes a sensitive measure of the degree of mixing, providing a window into the complex interplay of forces inside the nucleus.

Perhaps the most profound application of the ideas underpinning the Nilsson model lies in understanding nuclear fission and the very existence of the heaviest elements. The classical Liquid Drop Model predicts a single, smooth barrier that a nucleus must overcome to split apart. This picture, however, cannot explain phenomena like fission isomers—metastable states that exist in a highly deformed configuration.

The resolution comes from shell corrections. As a nucleus deforms towards fission, its single-particle levels, as mapped by a Nilsson-like diagram, shift dramatically in energy. Summing these energies provides a quantum correction to the smooth liquid-drop potential. This correction is not smooth; it creates hills and valleys on the path to fission. A "valley" at large deformation can create a second minimum in the potential, trapping the nucleus in an isomeric state. The height and shape of the multi-humped fission barrier are dictated entirely by these microscopic shell effects.

This very same physics is the foundation for the predicted "island of stability" for superheavy elements. These elements, with over 100 protons, should not exist according to the simple liquid drop model; they would instantly fly apart. Their predicted (and in some cases, observed) existence is a pure quantum mechanical effect. The specific arrangements of protons and neutrons create a strong, negative shell-correction energy that digs a deep hole in the potential energy surface, providing an extra barrier against spontaneous fission. The Nilsson model, in its generalized forms, is the primary tool used to map this unknown territory and guide the experimental search for these exotic forms of matter.

From the subtle shift in a magnetic moment to the cataclysmic process of fission, the Nilsson model serves as our indispensable guide. It reveals a world of profound beauty and unity, where the shape of a nucleus, the spin of its ground state, the rates of its decays, its response to being struck or spun, and its ultimate stability are all intricately connected, all governed by the quantum dance of nucleons in a deforming field.