🔬 physics

A Discrete Occupation-Flow Diagnostic for Orbital-Filling Crossovers in sd-Shell Nuclei: A Multi-Interaction Benchmark in Ne, Mg, and Si

This paper introduces a discrete orbital occupation-flow diagnostic, specifically the curvature $\kappa(N)$ , to map sequential orbital crossovers in Ne, Mg, and Si isotopes using USD-family interactions, revealing that this metric traces shell evolution distinct from collective coherence and highlighting the model's quantitative reliability except in the well-known island-of-inversion nuclei where the sd-shell space fails.

Original authors: Luay F. Sultan

Published 2026-08-28

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Luay F. Sultan

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Inside the heart of every atom lies a tiny, bustling city of protons and neutrons, arranged in layers much like the floors of a building. Physicists call this arrangement the nuclear shell model. Just as people fill up the rooms on a floor before moving to the next, neutrons fill specific energy levels, or "orbits," as they are added to an atom. For decades, scientists have been able to count how many neutrons sit in each orbit for any given atom, creating a static snapshot of the atom's interior. However, these snapshots tell us where the neutrons are, but they do not show how the atoms change as they grow. They miss the story of the transition: how the neutrons rearrange themselves as the atom adds new members, shifting from one orbit to another. Understanding these shifts is crucial because they reveal the fundamental rules that govern how matter holds together, especially in the region of the periodic table where atoms are known to behave strangely, such as in the "island of inversion," a zone where standard rules seem to break down.

In a recent study, a researcher named Luay F. Sultan sought to capture this missing story of change. Instead of just looking at the final count of neutrons in each orbit, he developed a new way to watch the flow of neutrons as they move from one atom to the next in a family of elements known as neon, magnesium, and silicon. He treated the addition of neutrons not as a series of isolated snapshots, but as a continuous journey. By calculating the difference in neutron numbers between atoms that are two steps apart, he created a vector, or a directional arrow, that points toward the specific orbit where the new neutrons are settling. He then measured the angle between these arrows as the atoms grew larger. If the neutrons kept filling the same orbit, the arrows would point in the same direction. But if the neutrons suddenly switched to a different orbit, the arrows would turn sharply. This turning angle became a new diagnostic tool, a precise measure of exactly when and how the internal structure of the atom reorganizes itself.

The researcher tested this method using three different, highly respected mathematical models that describe how neutrons interact with one another. He applied these models to the entire range of neon, magnesium, and silicon isotopes, from their lightest forms up to the point where the neutron shell is completely full. The results were strikingly consistent. No matter which of the three models was used, the turning angle pointed to the exact same moments in the growth of these atoms. The neutrons consistently filled the lowest orbit first, then switched to a middle orbit, and finally moved to the highest orbit. The new tool successfully mapped these switches with a high degree of precision, identifying the specific number of neutrons where the dominant filling orbit changed. This confirmed that the tool was not an artifact of a single mathematical choice but a robust feature of how these atomic shells actually behave.

Perhaps the most surprising finding was what this new tool did not correlate with. For a long time, scientists have used the strength of electromagnetic transitions—how easily an atom can change its shape and emit energy—as a sign of how much the nucleus is changing its structure. It was widely expected that the sharp turning of the neutron flow would happen at the same time these electromagnetic signals changed, suggesting a link between the rearrangement of orbits and the collective behavior of the nucleus. However, the study found no such connection. The sharp turns in the neutron flow happened independently of the electromagnetic signals. This means that the way neutrons shuffle between orbits is a distinct physical process from the way the nucleus as a whole vibrates or deforms. The two phenomena, while happening in the same atom, are not driving each other in the way previously assumed.

The study also revealed a specific point of instability in these atomic families. Around the time the nucleus contained sixteen neutrons, the different mathematical models began to disagree with each other regarding exactly how the neutrons were distributed. This suggests that sixteen neutrons is a critical threshold where the internal structure is most sensitive to the details of the forces holding the nucleus together. Furthermore, when the researcher compared his calculations to real-world experimental data, the tool worked perfectly for most atoms, accurately predicting their behavior. However, for two specific heavy isotopes of neon and magnesium, the standard model failed to match reality. This was not a failure of the new tool, but a confirmation of its limits: these two specific atoms are known to be "islands of inversion," where neutrons jump into a completely different set of orbits that the standard model cannot see. The tool correctly identified that the standard rules were breaking down in these cases, highlighting exactly where our current understanding needs to expand to include more complex interactions.

Ultimately, this work provides a new, clear lens for watching the internal life of atoms. By turning static counts into a dynamic map of flow and direction, the researcher has shown that we can pinpoint exactly when an atom's internal architecture shifts. The study demonstrates that the rearrangement of individual neutrons is a precise, measurable event that follows its own rules, separate from the broader, collective movements of the nucleus. While the tool works beautifully within the known boundaries of the atomic shell, it also clearly marks the edges where our current knowledge ends and where the strange, complex world of exotic nuclei begins.

Technical Summary: A Discrete Occupation-Flow Diagnostic for Orbital-Filling Crossovers in sd-Shell Nuclei

Problem Statement
In the nuclear shell model, orbital occupation numbers are traditionally reported as static snapshots (tables of $n(1d_{5/2})$ , $n(2s_{1/2})$ , $n(1d_{3/2})$ ) for specific isotopes. While accurate, this representation fails to explicitly characterize the dynamics of how a nucleus reorganizes its valence configuration as neutrons are added. Existing diagnostics for shell evolution—such as excitation energies $E(2_1^+)$ , transition strengths $B(E2)$ , and effective single-particle energies (ESPE)—either integrate over the full many-body wavefunction or provide single-N snapshots. There is a lack of a diagnostic tool that explicitly tracks the directional, step-to-step reorganization of orbital filling along an isotopic chain without requiring additional computational machinery beyond standard shell-model outputs.

Methodology
The author introduces a discrete dynamical quantity called the orbital occupation-flow vector, $\vec{F}(N)$ , defined as the finite difference in neutron occupations of the three sd-shell orbitals ( $1d_{5/2}$ , $2s_{1/2}$ , $1d_{3/2}$ ) between isotopes with neutron numbers $N$ and $N+2$ :
$\vec{F}(N) = (\Delta n_{d5/2}, \Delta n_{s1/2}, \Delta n_{d3/2})$
where $\Delta n_j = n_j(N+2) - n_j(N)$ .

From this vector, two scalar indicators are derived:

Curvature ( $\kappa(N)$ ): The angle (in degrees) between consecutive flow vectors $\vec{F}(N)$ and $\vec{F}(N+2)$ . This measures the abruptness of the change in the dominant filling orbital.
Divergence ( $\nabla \cdot \vec{F}(N)$ ): The change in the magnitude of the flow vector between steps, indicating acceleration or deceleration in the reorganization rate.

The study utilizes the NuShellX code to perform shell-model calculations in the sd valence space (outside a $^{16}$ O core) for the even-even isotopic chains of Neon (Ne), Magnesium (Mg), and Silicon (Si) up to the $N=20$ shell closure. Three independent USD-family interactions (USDA, USDB, and USDC) are employed to test the robustness of the diagnostic. The author also computes monopole-corrected effective single-particle energies (ESPE) derived directly from the interaction Hamiltonians to ground the flow behavior in the underlying physics.

Key Results

Robustness of Curvature: The curvature $\kappa(N)$ exhibits a highly reproducible pattern across all three interactions (USDA, USDB, USDC), with coefficients of variation (CV) typically below 9% and below 2% at dominant peaks. This confirms that the sequential crossover of filling orbitals is a robust feature of the sd-shell structure, not an artifact of a specific interaction parametrization.
Identification of Crossovers: $\kappa(N)$ identifies specific neutron numbers where the dominant filling orbital switches. In Ne and Si, a dominant peak occurs at $N=14$ (marking the $1d_{5/2} \to 2s_{1/2}$ crossover), while in Mg, a peak occurs at $N=12$ . A secondary crossover ( $2s_{1/2} \to 1d_{3/2}$ ) is observed near $N=16-18$ .
Independence from Collectivity: Contrary to the initial hypothesis, statistical tests (Pearson and Spearman correlations) reveal no systematic correlation between $\kappa(N)$ and $B(E2; 0_1^+ \to 2_1^+)$ . The author concludes that $\kappa$ tracks the identity of the orbital receiving neutrons (single-particle flow), whereas $B(E2)$ tracks the coherence of the resulting quadrupole deformation (collective coherence). These are distinct physical aspects of shell evolution.
Interaction Sensitivity: An interaction-sensitivity map, based on the summed standard deviation of occupation numbers across the three interactions, identifies $N=16$ as the point of maximal interaction dependence in all three chains. This corresponds to the region where the $2s_{1/2} \to 1d_{3/2}$ crossover is actively progressing and where the relative single-particle energies of the interactions differ most significantly.
Validation and Limitations: Comparison with experimental data shows the sd-shell calculations are quantitatively reliable for $E(2_1^+)$ and $B(E2)$ away from $N=20$ . However, as expected, the model fails qualitatively for $^{30}$ Ne and $^{32}$ Mg (the "island of inversion"), where it cannot reproduce the low $E(2_1^+)$ and high $B(E2)$ caused by $pf$-shell intruder configurations. The model correctly fails here because the sd-shell space excludes these intruder states.
Hamiltonian Connection: Monopole-based ESPE calculations confirm that the peaks in $\kappa(N)$ correspond to the narrowing of the $1d_{3/2} - 2s_{1/2}$ energy gap driven by the neutron-neutron monopole interaction, grounding the flow diagnostic in the Hamiltonian.

Significance and Claims
The paper claims that the occupation-flow vector $\vec{F}(N)$ and its curvature $\kappa(N)$ provide a novel, low-cost diagnostic for shell evolution. Its primary significance lies in:

Dynamic Insight: It transforms static occupation tables into a directional, step-by-step map of orbital reorganization, revealing "trajectories" of filling that static data obscures.
Physical Distinction: It establishes that orbital-filling crossovers (tracked by $\kappa$ ) and collective structural changes (tracked by $B(E2)$ ) are not necessarily coupled, challenging the assumption that a sharp change in filling implies a simultaneous change in collectivity.
Diagnostic Utility: It identifies $N=16$ as a critical region of model sensitivity in the sd-shell, where predictions for occupation composition are most dependent on the choice of interaction, even if derived observables like $B(E2)$ remain stable.
Generality: The framework is presented as a general methodology applicable to any isotopic chain within a fixed shell-model space, requiring no additional computational cost beyond standard diagonalization outputs.

The author explicitly states that $\kappa$ is not a parameterization of occupation numbers but a distinct quantity that captures trajectory information. They also note that while $\kappa$ is robust within the sd-shell, its application to the island of inversion requires a cross-shell model space (sd-pf), which is beyond the scope of the current benchmark. Future work is proposed to extend this to static quadrupole moments ( $Q(2_1^+)$ ) and to include proton-neutron monopole terms in ESPE calculations for a fully predictive test.

Technical Summary: A Discrete Occupation-Flow Diagnostic for Orbital-Filling Crossovers in sd-Shell Nuclei

More like this