Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

This paper extends spherical coupled-cluster theory to open-shell nuclei with two nucleons removed, validating the method against experimental data for oxygen and calcium isotopes while demonstrating high accuracy for binding energies and excited states but noting an underestimation of electric dipole polarizabilities.

The Gist

Imagine the atomic nucleus as a bustling city made of tiny citizens called protons and neutrons. In some cities, the population is perfectly balanced, with every street (energy level) either completely full or completely empty. These are "closed-shell" nuclei, and scientists have been very good at mapping them out.

But many nuclei are "open-shell," meaning they have a few extra citizens or a few missing ones, leaving the streets partially empty or partially full. This makes them much harder to study because the citizens interact in messy, unpredictable ways.

This paper is about a new, clever way to map these messy, open-shell cities using a method called Coupled-Cluster Theory. Here is how the authors did it, explained simply:

1. The "Neighbor" Trick

Instead of trying to solve the messy open-shell city directly, the authors decided to look at it as a "neighbor" of a perfect, closed-shell city.

• The Analogy: Imagine you want to understand a house with two missing bricks (an open-shell nucleus). Instead of analyzing the broken house from scratch, you start with the perfect, intact house next door (the closed-shell nucleus).
• The Method: They use a mathematical "excitation operator" to simulate removing two bricks (two particles) from the perfect house. This allows them to describe the broken house as an "excited state" of the perfect one. This is called the Two-Particle-Removed (2PR) method.

2. Building the Map (Ground State Energies)

First, they tested if this "neighbor trick" could accurately predict how heavy (or tightly bound) these nuclei are.

• The Result: They looked at Oxygen and Calcium isotopes (different versions of these elements). When they included more complex interactions (like accounting for triplets of particles moving together, not just pairs), their predictions became incredibly accurate.
• The Takeaway: For the basic structure and weight of these nuclei, their new method works just as well as the established methods used for perfect, closed-shell nuclei. It matches experimental data very closely.

3. Predicting the "Vibe" (Excited States)

Next, they tried to predict what happens when these nuclei get "excited" (like when a city lights up or vibrates).

• The Challenge: Some states are easy to predict (like a simple vibration), but others are tricky because they involve complex cross-talk between different energy levels.
• The Result:
  • For simple states (like in Carbon-14 or Oxygen-22), the method worked beautifully, correctly predicting the order and energy of the excited states.
  • For very complex, "negative parity" states (a specific type of quantum vibration), the method struggled a bit, overestimating the energy. This suggests that for these specific, messy states, they might need to add even more layers of complexity to their math in the future.

4. The "Sponge" Test (Electric Dipole Polarizability)

Finally, they tested how these nuclei react to an external electric field. Think of this as seeing how much a sponge squishes when you squeeze it. In physics, this is called Electric Dipole Polarizability.

• The Setup: They used a technique called the Lorentz Integral Transform (LIT), which is like a special filter that helps them see the "squishiness" of the nucleus without getting lost in the infinite possibilities of breaking it apart.
• The Result: Here is where they hit a snag. While their method worked great for the weight and structure of the nuclei, it consistently underestimated how "squishy" the Calcium isotopes were compared to real-world experiments.
• Why? The math showed that their method was missing some of the low-energy "wiggles" or "soft modes" that happen in these nuclei. It's as if their map showed the city as being stiffer than it actually is. They suspect they need to include even higher-order interactions (more complex particle groupings) to fix this.

Summary

The authors successfully built a new mathematical tool to study "imperfect" atomic nuclei by treating them as slightly modified versions of "perfect" ones.

• What worked: They can now predict the weight and basic energy levels of these nuclei with high accuracy, rivaling the best existing methods.
• What needs work: When predicting how these nuclei react to electric fields (specifically in Calcium), the method is a bit too "stiff" and misses some of the softer, low-energy behaviors seen in real life.

The paper concludes that this approach is a powerful, unified way to study open-shell nuclei, but to get the electric reaction perfect, they will need to add even more detailed layers of complexity to their calculations in the future.

Technical Summary

Technical Summary: Structure and Dynamics of Open-Shell Nuclei from Spherical Coupled-Cluster Theory

Problem Statement
While ab initio nuclear theory has achieved significant success with closed-shell nuclei using the Coupled-Cluster (CC) method, extending these capabilities to open-shell systems remains a challenge. Traditional approaches often rely on symmetry-breaking techniques (e.g., Bogoliubov CC or deformed CC) followed by symmetry restoration, which significantly increases computational complexity. Alternatively, Equation-of-Motion (EOM) CC offers a symmetry-preserving framework where open-shell nuclei are treated as excited states of a closed-shell core. While the two-particle-attached (2PA) EOM-CC method has been established for nuclei with two extra nucleons, and recently coupled with the Lorentz Integral Transform (LIT) for response functions, a corresponding two-particle-removed (2PR) formalism for nuclei with two fewer nucleons was lacking. This paper addresses the need to extend the spherical coupled-cluster framework to describe the structure and dipole response of open-shell nuclei formed by removing two nucleons from a shell closure.

Methodology
The authors develop and implement the two-particle-removed (2PR) EOM-CC formalism within the spherical coupled-cluster framework, coupled with the Lorentz Integral Transform (LIT) technique to compute response functions.

1. Formalism:

   • Ground State: The method utilizes a single-reference CC approach (CCSD and CCSDT-1) to determine the ground state of a closed-shell reference nucleus (mass $A$). A mass shift is applied such that the reference energy corresponds to the target nucleus (mass $A-2$).
   • 2PR-EOM: The open-shell target nucleus is described by an excitation operator that removes two fermions from the reference. The ansatz includes $0p-2h$ (two-hole) configurations at leading order, complemented by $1p-3h$ contributions. The right ($\hat{R}$) and left ($\hat{L}$) excitation operators are constructed to solve the non-Hermitian eigenvalue problem for excitation energies relative to the mass-shifted reference.
   • Response Functions: To calculate electromagnetic responses, the 2PR-EOM ansatz is combined with the LIT technique. This involves solving inhomogeneous Schrödinger-like equations for auxiliary states using the Lanczos method. The electric dipole polarizability ($\alpha_D$) is extracted from the LIT moments.
   • Interactions: Calculations employ chiral effective field theory interactions ($\Delta$NNLOGO(394)) including two- and three-nucleon forces.

2. Truncation Schemes:

   • Ground state calculations are performed at the CCSD (denoted D) and CCSDT-1 (denoted T-1) levels.
   • EOM calculations vary the truncation of the excitation operators (e.g., $0p-2h$ vs. $1p-3h$ for 2PR).
   • Uncertainties are estimated by comparing different truncation schemes (e.g., D vs. T-1 for ground states; different EOM truncations for open-shell systems).

Key Contributions

• Development of 2PR-EOM-CC: The paper presents the detailed formalism for the two-particle-removed EOM-CC method, enabling the study of nuclei with two holes relative to a shell closure.
• Integration with LIT: The authors extend the 2PA-LIT-CC technique to the 2PR sector, allowing for the ab initio calculation of electric dipole response functions and polarizabilities in open-shell nuclei.
• Comprehensive Validation: The method is validated across the oxygen and calcium isotopic chains, covering ground-state energies, selected excited states, and electric dipole polarizabilities.

Results

• Ground-State Energies:
  • For $^{22}$O and $^{38}$Ca, the 2PR-EOM approach with $1p-3h$ excitations (D/1p-3h) recovers approximately 85% of the experimental binding energy.
  • Including triples in the reference ground state (T-1/1p-3h) significantly improves accuracy, bringing results within $\sim$1% of experimental values and matching the precision of established closed-shell CCSDT-1 calculations.
  • Consistency checks between 2PA and 2PR schemes for $^{38}$Ca show convergence to similar energies when higher-order excitations are included.
• Excited States:
  • In $^{14}$C, the 2PR-EOM method accurately describes positive-parity states (e.g., $2^+$) but struggles with negative-parity states ($1^-$), which require cross-shell configurations and higher-order correlations not fully captured by the current truncation.
  • In $^{22}$O, the 2PR-EOM approach successfully reproduces the ordering of the five lowest-lying excited states, outperforming the conventional closed-shell ansatz in this specific case.
• Electric Dipole Polarizability ($\alpha_D$):
  • Calculations for $^{22}$O and $^{38}$Ca show that 2PR-LIT-CC results generally underestimate experimental $\alpha_D$ values compared to closed-shell CC predictions.
  • The underestimation is attributed to a shift of the dipole strength to higher energies in the 2PR calculation, missing the low-energy "soft" dipole mode contributions seen in closed-shell and SCGF calculations.
  • For the calcium chain, both 2PA and 2PR schemes tend to underestimate $\alpha_D$, suggesting that higher-order many-body contributions are necessary for accurate response function predictions.

Significance and Claims
The paper claims to provide a unified, symmetry-preserving framework for studying open-shell nuclei near shell subclosures using the spherical coupled-cluster method. The authors assert that the 2PR-EOM-CC approach achieves accuracy comparable to closed-shell CCSDT-1 for binding energies and low-lying excited states when appropriate truncations are used. However, they modestly note that while the method is robust for structural properties, it currently underestimates electric dipole polarizabilities, particularly in calcium isotopes. This discrepancy highlights the need for higher-order many-body correlations in response function calculations. The work establishes the 2PR technique as a viable tool for expanding the reach of ab initio nuclear theory to open-shell systems, with future work planned to include higher-order excitations (up to 2p-4h) and to apply one-particle-attached/removed schemes to odd nuclei.