Ab initio mapping of the boundary of the $N=20$ island of inversion

Using an ab initio approach combining the multi-reference in-medium similarity renormalization group and the quantum-number projected generator coordinate method with chiral two- plus three-nucleon interactions, this study systematically reproduces the low-lying properties of neutron-rich nuclei around $N=20$ and delineates the specific boundaries of the $N=20$ island of inversion.

The Gist

Imagine the atomic nucleus not as a static, hard marble, but as a squishy, dynamic blob of dough. In the world of nuclear physics, most of these "dough balls" have a preferred shape: a perfect sphere. This happens when the protons and neutrons inside are neatly packed into specific energy layers, much like electrons in an atom. When these layers are full, the nucleus is stable and spherical. These are the "magic" nuclei.

However, there is a strange region on the nuclear map, around the element Magnesium (specifically where there are 20 neutrons), where the rules change. Here, the dough doesn't want to be a sphere; it wants to stretch out into a football shape. This region is called the "Island of Inversion" (IOI). It's called an "island" because it's surrounded by "normal" spherical nuclei, and an "inversion" because the usual rules of stability are flipped upside down.

For decades, scientists have been trying to draw the exact coastline of this island: Which specific atoms belong to this squishy, deformed world, and which ones are still normal spheres?

This paper, written by a team of physicists, uses a powerful new computational tool to map out that coastline with high precision. Here is how they did it, explained simply:

The Problem: The "Shape-Shifting" Puzzle

In this region of the nuclear chart, things get messy. Some nuclei are spherical, some are football-shaped, and some are a weird mix of both (like a nucleus that is simultaneously a sphere and a football). This is called shape coexistence.

Previous methods were like trying to take a photo of a spinning, shape-shifting dancer with a camera that only takes blurry snapshots. They could guess the shape, but they often got the energy levels wrong or missed the subtle mix of shapes entirely.

The Solution: The "Super-Telescope" (IM-GCM)

The authors used a method called IM-GCM (In-Medium Generator Coordinate Method). Think of this as a high-tech, multi-lens telescope that can see the nucleus in three different ways at once:

1. The Microscope (IMSRG): First, they use a technique called the Similarity Renormalization Group to zoom in on the tiny, individual interactions between protons and neutrons. It's like taking a complex recipe and simplifying the ingredients so you can understand how they interact without getting overwhelmed by the noise.
2. The Shape-Shifter (PGCM): Next, they use the Generator Coordinate Method. Imagine the nucleus is a piece of clay. This method asks the computer to mold the clay into every possible shape (flat, round, long, short) and calculate the energy of each shape.
3. The Blender (Mixing): Finally, they don't just pick the "best" shape. They realize the nucleus is actually a quantum superposition—a blend of all those shapes at once. The IM-GCM "blends" these different shapes together to find the true, final state of the nucleus.

Crucially, they did this starting from the most fundamental laws of physics (chiral interactions) rather than making up rules to fit the data. This is an ab initio approach, meaning "from the beginning."

The Discovery: Mapping the Coastline

By running these massive simulations on supercomputers, the team mapped out the "Island of Inversion" for the first time with such clarity.

Who lives on the Island? (The Deformed Nuclei)
These nuclei are the "football players." They are strongly stretched out and unstable in the traditional sense. The paper confirms that the following are definitely inside the island:

• Neon-30
• Sodium-29, 31, 33
• Magnesium-31, 32, 33, 34
• Aluminum-35

Who lives on the Shore? (The Spherical Nuclei)
These are the "spherical" neighbors that stay round and don't join the deformation party. They fall outside the island:

• Fluorine-29
• Neon-29
• Magnesium-30
• Aluminum-31, 33
• Silicon-34, 35
• Phosphorus-35

Why Does This Matter?

Think of the "Island of Inversion" as a special zone where the laws of nuclear structure are rewritten. By accurately mapping the boundary, scientists can:

• Test our fundamental theories: It proves that our current understanding of how protons and neutrons stick together (the strong force) is correct, even in these extreme, exotic conditions.
• Predict the unknown: There are many nuclei in this region that we haven't discovered yet because they are hard to make in a lab. This map tells us what to expect when we finally find them.
• Understand the Universe: These exotic nuclei play a role in how stars explode (supernovae) and how heavy elements are forged in the cosmos. Knowing their properties helps us understand the history of the universe.

The Bottom Line

This paper is like drawing the first accurate map of a mysterious, shape-shifting continent. The researchers built a digital microscope that could see the nucleus not just as a static object, but as a dynamic, quantum dance of shapes. Their map confirms that the "Island of Inversion" is a real, distinct place where nuclear matter behaves differently, and they have successfully drawn the line between the "normal" world and this "inverted" one.

Technical Summary

Here is a detailed technical summary of the paper "Ab initio mapping of the boundary of the N = 20 island of inversion" by E. F. Zhou et al.

1. Problem Statement

The "Island of Inversion" (IOI) around $N=20$ is a region in the nuclear chart where neutron-rich nuclei exhibit large quadrupole deformations and low-lying $2^+_1$states, deviating significantly from the magic number behavior expected at$N=20$. While the existence of the IOI in nuclei like$^{32}\text{Mg}$ is well-established, the precise boundary of this region remains a subject of debate.

• The Challenge: Accurately describing these nuclei requires capturing both collective correlations (shape coexistence, shape mixing) and dynamical correlations (multi-particle multi-hole excitations across the shell gap).
• Limitations of Existing Methods:
  • Standard Valence-Space IMSRG (VS-IMSRG) often overestimates excitation energies and underestimates $B(E2)$ values, partly due to the lack of angular momentum projection in deformed frames.
  • Projected Generator Coordinate Method (PGCM) based on phenomenological functionals or limited valence spaces struggles to predict the correct energy ordering of competing weakly and strongly deformed bands (e.g., in $^{30}\text{Ne}$).
  • Existing ab initio approaches often struggle to simultaneously describe shape coexistence and the correct ground-state spin-parity in odd-mass nuclei within this region.

2. Methodology: The IM-GCM Framework

The authors employ the In-Medium Generator Coordinate Method (IM-GCM), a fully ab initio approach that combines two advanced many-body techniques:

1. Multi-Reference In-Medium Similarity Renormalization Group (MR-IMSRG):
   • Starts from a chiral two- plus three-nucleon interaction (EM1.8/2.0).
   • Uses an ensemble of symmetry-restored Hartree-Fock-Bogoliubov (HFB) reference states (determined via Variation After Particle-Number Projection, VAPNP) to capture static correlations and shape coexistence.
   • Evolves the Hamiltonian and all operators (including transition operators) via a flow equation to decouple the reference state from excited states, incorporating dynamic correlations.
   • Employs the NO2B approximation (truncating to normal-ordered two-body terms) to make calculations feasible while retaining essential physics.
2. Quantum-Number Projected Generator Coordinate Method (PGCM):
   • Constructs the final wavefunctions as linear combinations of projected HFB states.
   • Projections: Explicitly restores broken symmetries: Particle Number ($N, Z$) and Angular Momentum ($J$).
   • Collective Coordinates: Uses the quadrupole deformation parameter $\beta_2$ as the generator coordinate.
   • Solves the Hill-Wheeler-Griffin (HWG) equation to obtain the energy spectrum and mixing weights of states with different intrinsic shapes.
3. Treatment of Odd-Mass Nuclei:
   • Extends the framework to odd-mass systems by constructing one-quasiparticle excitations on top of the even-even core.
   • Handles time-reversal symmetry breaking in odd-mass HFB calculations and uses the Pfaffian formula for overlap integrals.

Key Technical Innovation: The authors consistently evolve tensor operators (for electromagnetic transitions) using the Baker-Campbell-Hausdorff (BCH) expansion within the MR-IMSRG flow. This ensures that induced two-body currents are included, avoiding the need for empirical effective charges.

3. Key Contributions

• Systematic Ab Initio Mapping: The study provides a comprehensive, parameter-free (using only chiral forces) mapping of the $N=20$ region, covering both even-even and odd-mass nuclei.
• Resolution of Shape Coexistence: The method successfully reproduces the coexistence of weakly and strongly deformed states at comparable energies, a phenomenon difficult to capture with mean-field or single-reference methods.
• Accurate Prediction of Ground States: The framework correctly predicts the spin-parity and deformation of ground states in odd-mass nuclei (e.g., $^{33}\text{Mg}$, $^{31}\text{Na}$), resolving previous discrepancies in shell-model and coupled-cluster predictions.
• Operator Evolution: Demonstrates the necessity of evolving transition operators (including induced two-body terms) to reproduce experimental $B(E2)$ and magnetic moments without effective charges.

4. Key Results

The study analyzed 17 nuclei around $N=20$ ($Z=9-15$).

A. Even-Even Nuclei:

• $^{30}\text{Ne}$: Predicted to have a strongly deformed ground state ($\beta_2 \approx 0.5$) and a weakly deformed excited $0^+_2$ state, confirming its location inside the IOI.
• $^{32}\text{Mg}$: Reproduces the low-lying $2^+_1$state and large$B(E2)$ values, characteristic of the IOI.
• $^{34}\text{Si}$: Predicted to have a weakly deformed ground state (spherical/oblate), placing it outside the IOI, consistent with recent experimental findings of a "magic" behavior.
• Shape Coexistence: Successfully reproduced the inversion of energy ordering between spherical and deformed $0^+$states in$^{30}\text{Ne}$and$^{34}\text{Si}$.

B. Odd-Mass Nuclei:

• $^{33}\text{Mg}$: Correctly predicted a $3/2^-$ground state arising from a strongly deformed configuration ($K^\pi = 3/2^-$), driven by the crossing of$\nu f_{7/2}$and$\nu d_{3/2}$ orbitals.
• Sodium Isotopes ($^{29,31,33}\text{Na}$):
  • $^{29}\text{Na}$: Weakly deformed ($3/2^+$).
  • $^{31}\text{Na}$: Shows a shape transition; the ground state is strongly deformed ($3/2^+$), while low-lying states exhibit shape mixing.
  • $^{33}\text{Na}$: Strongly deformed ground state.
• Aluminum and Silicon: $^{31,33}\text{Al}$ and $^{34,35}\text{Si}$ are predicted to be weakly deformed (outside IOI), while $^{35}\text{Al}$ shows moderate deformation.

C. Electromagnetic Observables:

• Calculated magnetic dipole moments ($\mu$) and spectroscopic quadrupole moments ($Q_s$) agree well with experimental data.
• The study highlights that the Schmidt formula works for weakly deformed nuclei but fails for strongly deformed ones, where configuration mixing is crucial.
• $B(E2)$ values are reproduced without effective charges, validating the consistency of the evolved operators.

5. Determination of the IOI Boundary

Based on the mean quadrupole deformation ($\bar{\beta}_2$) of the ground states, the authors define the boundary as follows:

• Inside the Island of Inversion (Strongly/Moderately Deformed, $\bar{\beta}_2 > 0.2$):

  • $^{30}\text{Ne}$
  • $^{29,31,33}\text{Na}$
  • $^{31,32,33,34}\text{Mg}$
  • $^{35}\text{Al}$

• Outside the Island of Inversion (Weakly Deformed/Spherical, $\bar{\beta}_2 < 0.2$):

  • $^{29}\text{F}$
  • $^{29}\text{Ne}$
  • $^{30}\text{Mg}$
  • $^{31,33}\text{Al}$
  • $^{34,35}\text{Si}$
  • $^{35}\text{P}$

6. Significance

This work represents a significant advancement in ab initio nuclear structure physics:

1. Validation of IM-GCM: It demonstrates that the IM-GCM, starting from chiral interactions, can successfully describe complex phenomena like shape coexistence and shape mixing in medium-mass nuclei without phenomenological adjustments.
2. Definitive Boundary: It provides a theoretically robust and experimentally consistent definition of the $N=20$ IOI boundary, resolving ambiguities regarding specific isotopes like $^{30}\text{Ne}$ and $^{34}\text{Si}$.
3. Methodological Benchmark: The consistent treatment of operator evolution and the extension to odd-mass nuclei set a new standard for future ab initio studies of exotic nuclei.
4. Future Directions: The authors note that while the current results are robust, extending the method to odd-odd nuclei and incorporating continuum effects (crucial for nuclei near the drip line) are necessary for a complete understanding of the IOI.