From closed shells to open shells: Coupled-cluster calculations of atomic nuclei

This study provides a comprehensive comparison of coupled-cluster formulations based on symmetry-broken reference states and equation-of-motion techniques, demonstrating that both approaches yield consistent descriptions of bulk nuclear properties across calcium and nickel isotopic chains when using chiral effective field theory interactions.

The Gist

Imagine the atomic nucleus as a bustling, crowded dance floor inside a tiny ballroom. The dancers are protons and neutrons, and their movements are governed by the "music" of the strong nuclear force. For decades, physicists have been trying to predict exactly how these dancers move, how tightly they hold hands, and what shapes the dance floor takes.

This paper is a report card on a specific set of mathematical tools called Coupled-Cluster (CC) theory. The authors are asking: "Can we use these tools to predict the behavior of nuclei that aren't perfectly organized, or do we need new tricks?"

Here is the breakdown of their journey from "closed shells" to "open shells," explained simply.

1. The Setup: The Perfect vs. The Messy Dance Floor

In the world of atoms, some nuclei are like a perfectly choreographed ballet troupe. Every dancer has a specific spot, and the formation is a perfect sphere. These are called closed-shell nuclei (like Calcium-40 or Nickel-56).

• The Easy Mode: Calculating these is relatively easy because the symmetry (the perfect order) makes the math simpler. It's like solving a puzzle where all the pieces are the same shape.

However, most nuclei are open-shell. Imagine a dance floor where the music changes, the dancers get crowded, and they start pairing up or forming irregular shapes (deformation).

• The Hard Mode: These are messy, dynamic systems. The old "perfect order" math breaks down. To solve this, physicists have developed three different "strategies" to handle the chaos.

2. The Three Strategies (The Tools)

The paper compares three different ways to calculate the energy of these messy nuclei:

• Strategy A: The "Neighborly Excitation" (EOM-CC)

  • The Analogy: Imagine you want to know what happens if you remove two dancers from a perfect troupe. Instead of re-calculating the whole dance from scratch, you take the perfect troupe (a closed-shell nucleus) and say, "Okay, let's pretend two dancers just left." You calculate the "excitation" or the ripple effect of that missing pair.
  • Pros: Very fast and accurate for nuclei right next to the perfect ones.
  • Cons: It struggles if the nucleus is too far from the "perfect" state. It's like trying to describe a chaotic mosh pit by only looking at a quiet line dance.

• Strategy B: The "Pairing Breaker" (Bogoliubov CC)

  • The Analogy: In some nuclei, neutrons like to pair up (superfluidity), kind of like couples holding hands and spinning together. This strategy breaks the rule that says "you must have an exact number of dancers." Instead, it allows for a "fuzzy" number of dancers to account for these pairs.
  • Pros: Great for nuclei where pairing is the main event.
  • Cons: It's computationally expensive (requires more computer power) because it has to track all those fuzzy pairs.

• Strategy C: The "Shape Shifter" (Deformed CC)

  • The Analogy: Sometimes the dance floor isn't a sphere; it's an egg or a football. This strategy admits that the nucleus is deformed from the start. It calculates the energy based on this squashed shape rather than trying to force it into a sphere.
  • Pros: Excellent for nuclei that are naturally weird shapes.
  • Cons: Also computationally heavy because you lose the mathematical shortcuts that come with perfect spheres.

3. The Experiment: Calcium and Nickel

The authors tested these three strategies on two families of elements: Calcium and Nickel. They used two different "rulebooks" (nuclear interactions) to see which strategy worked best.

The Results:

• They all agreed: Surprisingly, all three strategies gave almost the same answer for the total energy of the nucleus. Whether they used the "Neighborly Excitation," the "Pairing Breaker," or the "Shape Shifter," the results were consistent.
• The "Bulk" Properties: For big-picture things like how heavy the nucleus is or how tightly it holds together (binding energy), all methods work well.
• The "Drip Line": They looked at how many neutrons a nucleus can hold before it falls apart (the drip line). Their calculations suggest that Nickel can hold more neutrons than we thought, potentially extending the "island of stability."

4. The Big Takeaway

The paper concludes that Coupled-Cluster theory is a Swiss Army Knife.

• If you are near a perfect nucleus, use the Neighborly Excitation (it's fast).
• If you are in the middle of a messy, pairing-heavy nucleus, use the Pairing Breaker or Shape Shifter (they are more flexible).

The most exciting part? It doesn't matter which tool you pick for the big picture. They all paint the same consistent portrait of the nucleus. This gives physicists confidence that they can now use these tools to explore even heavier, stranger nuclei that we haven't been able to study before.

In a Nutshell

Think of this paper as a group of engineers testing three different blueprints to build a bridge over a turbulent river.

• Blueprint A uses the calm water upstream as a reference.
• Blueprint B accounts for the swirling eddies.
• Blueprint C accounts for the rocks in the middle.

They found that all three blueprints result in a bridge that stands just as strong. This means we can now build bridges (predict nuclear properties) across the most turbulent parts of the atomic river with confidence, paving the way to discover new elements and understand the universe's building blocks.

Technical Summary

1. Problem Statement

Ab initio nuclear structure calculations have successfully described closed-shell (magic) nuclei using Coupled-Cluster (CC) theory. However, the majority of atomic nuclei are open-shell, characterized by partially filled valence shells, deformation, and pairing correlations. Standard single-reference CC theory, which relies on a spherical, symmetry-preserving reference state (Hartree-Fock), struggles with these systems due to the breakdown of spherical symmetry and the presence of static correlations.

While extensions exist to handle open-shell nuclei, there is a need for a comprehensive benchmarking of different formulations to determine their consistency, accuracy, and computational efficiency across medium-mass isotopic chains. Specifically, the paper addresses the challenge of describing doubly open-shell nuclei (where both proton and neutron shells are open) and singly open-shell nuclei (where one shell is open) using a unified framework.

2. Methodology

The authors benchmark three distinct formulations of Coupled-Cluster theory applied to even-even Calcium ($Z=20$) and Nickel ($Z=28$) isotopes using chiral effective field theory ($\chi$EFT) interactions (specifically EM 1.8/2.0 and $\Delta$NNLOGO(394)).

The three approaches compared are:

1. Equation-of-Motion CC (EOM-CC):

   • Strategy: Treats open-shell nuclei as excitations (particle-attached or particle-removed) from a neighboring closed-shell (or semi-magic) reference nucleus.
   • Implementation: Uses spherical CCSD for the reference state. The target nucleus is accessed via 2-particle-attached (2PA) or 2-particle-removed (2PR) excitation operators (including up to 3p-1h or 1p-3h configurations).
   • Limitation: Restricted to nuclei near shell closures; less effective for mid-shell nuclei with strong deformation.

2. Bogoliubov CC (BCC):

   • Strategy: Employs a reference state that breaks particle-number symmetry (U(1) gauge symmetry) to capture pairing correlations.
   • Implementation: Uses a spherical Hartree-Fock-Bogoliubov (HFB) vacuum. The cluster operator is expanded in terms of quasi-particles.
   • Target: Primarily suited for singly open-shell nuclei where pairing is dominant but deformation is weak.

3. Deformed CC (CC on Deformed Reference):

   • Strategy: Employs a reference state that breaks rotational symmetry (SO(3)) to capture static deformation.
   • Implementation: Uses an axially symmetric deformed Hartree-Fock (HF) Slater determinant. The calculation proceeds without angular momentum projection (unprojected) to maintain computational feasibility, though projection is noted as a future step for spectroscopy.
   • Target: Suited for doubly open-shell nuclei with significant deformation.

Computational Details:

• Basis: Spherical harmonic oscillator basis with $e_{max}=12$ (13 major shells).
• Truncation: All calculations are performed at the CCSD (Singles and Doubles) level.
• Three-Body Forces: Treated via the Normal-Ordered Two-Body (NO2B) approximation.
• Uncertainty Estimation: Theoretical uncertainties are estimated based on model-space dependence and the omission of triple excitations (approximated as 10% of the BCCSD correlation energy).

3. Key Contributions

• Comprehensive Benchmarking: The first large-scale comparison of EOM-CC, BCC, and Deformed CC across the same isotopic chains (Ca and Ni) using identical interactions and truncation levels.
• Validation of Symmetry-Broken Approaches: Demonstrates that symmetry-broken reference states (deformed HF and spherical HFB) provide consistent descriptions of bulk properties (ground-state energies, separation energies) comparable to EOM-CC near shell closures.
• Extension to Mid-Shell: Highlights that while EOM-CC is efficient near magic numbers, symmetry-broken methods are essential for describing mid-shell isotopes where EOM ansatzes fail due to the complexity of the excitation space.
• Quantification of Triple Contributions: Provides estimates for the impact of missing triple excitations, showing they are significant for absolute energies but largely cancel out in differential observables.

4. Key Results

Ground-State Energies

• Consistency: For both Ca and Ni chains, the different CC methods (CCSD, BCCSD, 2PA-EOM, 2PR-EOM) yield ground-state energies that agree within the estimated theoretical uncertainties (gray bands in figures).
• Interaction Dependence: The EM 1.8/2.0 interaction generally yields results closer to experimental data than the harder $\Delta$NNLOGO(394) potential. The discrepancy with experiment is largely attributed to the missing triple excitations (estimated at ~20 MeV for $^{48}$Ca with the harder potential).
• Symmetry Breaking: In Calcium ($Z=20$), deformation is minimal ($\beta \lesssim 0.1$), so both BCC and Deformed CC perform similarly. In Nickel ($Z=28$), discrepancies between CCSD and BCCSD reach 2–3 MeV in mid-shell isotopes ($^{62,64}$Ni), reflecting the absence of a subshell closure at $N=34$.

Two-Neutron Separation Energies ($S_{2n}$)

• Error Cancellation: As differential quantities, $S_{2n}$ values show excellent agreement with experiment even at the CCSD level, as systematic errors (model space truncation, missing triples) cancel out.
• Shell Closures: All methods correctly reproduce the sharp drops in $S_{2n}$ at magic numbers ($N=20, 28, 40, 50$).
• Ni Chain: The methods accurately describe the trend in the Ni chain, confirming the magicity of $^{56}$Ni, $^{68}$Ni, and $^{78}$Ni. The BCCSD method predicts $^{80}$Ni to be bound, suggesting the neutron drip line lies beyond $N=52$.

Two-Neutron Shell Gaps ($\Delta_{2n}$)

• Magic Numbers: Both CCSD and BCCSD reproduce the strong shell gaps at $N=28$ and $N=50$.
• Discrepancies:
  • At $N=28$ ($^{56}$Ni), both methods overpredict the experimental shell gap by ~30%.
  • At $N=40$ ($^{68}$Ni), the calculations predict a strong magic character (enhanced gap), whereas experimental data does not show a clear magic signature. This suggests CCSD tends to over-stabilize closed shells in this region.
  • The methods agree well for $A \leq 56$ and $A \geq 70$, but show larger deviations in the open-shell region between $N=28$ and $N=40$.

5. Significance and Outlook

• Robustness of CC Theory: The study confirms that Coupled-Cluster theory is a robust tool for describing medium-mass nuclei, capable of handling both spherical and deformed systems through different symmetry-breaking strategies.
• Method Selection:
  • EOM-CC is optimal for nuclei near closed shells due to its conceptual simplicity and lower computational cost (using $j$-scheme).
  • Symmetry-Broken CC (BCC/Deformed) is necessary for mid-shell nuclei and systems with strong static correlations, offering a scalable path to heavier nuclei where EOM-CC becomes intractable.
• Future Directions: The authors emphasize the need for a unified framework that incorporates both deformation and pairing correlations within a single reference state. Additionally, the restoration of broken symmetries (angular momentum and particle number projection) is identified as a critical next step for improving spectroscopic predictions, although it is less critical for bulk properties like binding energies.

In conclusion, the paper establishes that different CC formulations offer consistent descriptions of nuclear bulk properties across open-shell isotopic chains, validating the use of symmetry-breaking reference states as a powerful extension of ab initio nuclear structure theory.