Development of an accurate formalism to predict properties of two-neutron halo nuclei: case study of $^{22}$C

This paper develops and validates an accurate, efficient three-body hyperspherical harmonics formalism combined with the R-matrix method to predict the properties of the two-neutron halo nucleus $^{22}$C, demonstrating that the projection method outperforms the supersymmetric approach in enforcing the Pauli principle while introducing algorithmic optimizations that significantly reduce computational costs.

The Gist

The Big Picture: The "Ghostly" Nucleus

Imagine an atom's nucleus not as a solid marble, but as a tiny solar system. Usually, the planets (protons and neutrons) hug the sun (the core) tightly. But in some very unstable, exotic atoms, the outer planets are so loosely held that they drift far away, creating a giant, fuzzy cloud around the core.

This is called a Halo Nucleus. Specifically, this paper studies Carbon-22 ($^{22}$C), which has a core of Carbon-20 and two "ghostly" neutrons floating far out in a halo.

The problem? These two neutrons are so far out that they act like a single unit (a few-body system), but they are also made of individual particles that must obey strict quantum rules (the Pauli Exclusion Principle). It's like trying to describe a dance where the dancers are both a single couple and individual people who can't occupy the same space.

The Goal: Building a Better Calculator

The authors wanted to build a super-accurate computer model to predict how this Carbon-22 nucleus behaves. They had two main goals:

1. Fix the Rules: Make sure the model correctly enforces the "no two particles in the same seat" rule (Pauli Principle) without breaking the math.
2. Speed It Up: Make the calculations fast enough so scientists can run them thousands of times to understand the uncertainty in their predictions.

The Two Methods: The "Magic Eraser" vs. The "Bouncer"

To stop the model from predicting impossible states (where neutrons sit in spots they shouldn't), the authors tested two different mathematical tricks. Think of these as two ways to manage a crowded party:

1. The Supersymmetric Method (The "Magic Eraser")

• How it works: This method takes the math describing the particles and subtly tweaks the "walls" of the room. It adds a tiny, invisible repulsive force that pushes the forbidden particles out of the wrong spots.
• The Analogy: Imagine you have a crowded dance floor. Instead of kicking people out, you magically make the floor slippery in the wrong spots so no one can stand there. It's computationally cheap and easy to do, but it changes the texture of the floor.
• The Result: It worked okay, but it made the "dance floor" slightly different than reality, leading to predictions that were a bit off.

2. The Projection Method (The "Bouncer")

• How it works: This method is more direct. It calculates exactly which states are forbidden and then mathematically "projects" them out of the equation, like a bouncer explicitly checking IDs and removing the wrong people before the party starts.
• The Analogy: You have a list of who can't come in. You check the list and physically remove those people from the room. It's harder work and takes more time, but the room is exactly as it should be.
• The Result: This turned out to be the winner. It gave a more accurate picture of the nucleus.

The Findings: Why the "Bouncer" Won

When they compared the two methods using Carbon-22, they found some surprising differences:

• Size Matters: The "Magic Eraser" method predicted the nucleus was slightly larger than it actually is. The "Bouncer" method gave the correct size.
• The Shape of the Dance: The two methods predicted the neutrons would arrange themselves differently.
  • The Projection (Bouncer) method showed the neutrons could be in a "cigar" shape (stretched out) or a "di-neutron" shape (hugging each other).
  • The Supersymmetric (Eraser) method forced the neutrons into a tight hug (di-neutron) and missed the other possibilities.
• The "Nodes": The most technical proof was looking at the "waves" of the particles. The Projection method showed extra "bumps" (nodes) in the waves, which is exactly what you expect when you correctly enforce the rule that particles can't overlap. The Eraser method smoothed these bumps out, losing important physical details.

The Verdict: If you want to know the truth about these exotic nuclei, you must use the Projection Method. It's harder to code, but it's the only one that tells the whole story.

The Speed Hack: Cutting the Fat

The second part of the paper was about making the math run faster. Calculating these nuclei is like trying to solve a puzzle with millions of pieces.

• The Problem: To get a perfect answer, you need to consider every possible way the particles can move. This takes forever.
• The Solution: The authors realized that some of those "pieces" of the puzzle don't matter much. They found a way to cut out the "unimportant" pieces (specifically, certain high-energy movements) without losing accuracy.
• The Result: They managed to cut the computer time by 20%. This might not sound like much, but in the world of super-computing, that's huge. It means scientists can run more tests, check more uncertainties, and get better answers faster.

Why This Matters

This paper isn't just about Carbon-22. It's about building a better toolkit for the future.

1. Better Physics: By proving the "Bouncer" method is superior, they are telling other scientists: "Stop using the easy shortcut; use the accurate one."
2. Future Discoveries: As we build better particle accelerators (like FRIB), we will discover even stranger, heavier halo nuclei. This new, fast, and accurate code is the foundation we need to understand them.
3. Uncertainty Quantification: Science isn't just about getting an answer; it's about knowing how sure we are of that answer. By making the code faster, scientists can run thousands of simulations to say, "We are 95% sure the nucleus is this big," rather than just guessing.

In short: The authors built a faster, more accurate microscope for looking at the weird, fuzzy edges of the atomic world, and they proved that doing the hard work (Projection) gives you a clearer picture than taking the easy shortcut (Supersymmetry).

Technical Summary

1. Problem Statement

The study addresses the challenge of accurately modeling two-neutron halo nuclei (specifically $^{22}$C), which are loosely bound systems exhibiting a transition from many-body to few-body dynamics.

• The Core Conflict: Standard microscopic models (using harmonic oscillator bases) struggle to describe the large spatial extension of halo nuclei and the continuum states. Conversely, few-body models (treating the nucleus as a core plus valence neutrons) naturally handle the spatial extension but fail to inherently reproduce shell structure and the Pauli exclusion principle.
• The Specific Challenge: In few-body models, the internal structure of the core is ignored, leading to "Pauli-forbidden" states where valence neutrons occupy orbitals already filled by core nucleons. These must be removed to avoid spurious bound states. Two primary methods exist for this removal:
  1. Supersymmetric (SUSY) Transformation: Modifies the potential to remove bound states without altering phase shifts.
  2. Projection Method: Explicitly projects out forbidden states from the Hilbert space.
• Knowledge Gap: Previous comparisons of these methods were limited to light systems (e.g., $^6$He) with only one forbidden state per subsystem. It was unclear which method is more accurate for heavier systems like $^{22}$C (where valence neutrons are in the $sd$-shell and multiple forbidden states exist) or how to implement these methods efficiently for uncertainty quantification.

2. Methodology

The authors developed a new computational framework, hyperboromir, based on the Hyperspherical Harmonics (HH) formalism combined with the R-matrix method.

• Formalism:
  • The three-body system (Core + $n$ + $n$) is described using Jacobi coordinates and hyperspherical coordinates ($\rho, \alpha, \Omega$).
  • The wavefunction is expanded in hyperspherical harmonics up to a cutoff $K_{max}$.
  • R-matrix Method: Used to solve the hyperradial Schrödinger equation. This allows for the consistent treatment of bound states and scattering states without discretizing the continuum, enforcing appropriate boundary conditions at a channel radius $\rho_{max}$.
• Pauli Exclusion Implementation:
  • SUSY Method: Adds a short-range repulsive potential derived from the logarithmic derivative of the forbidden bound state wavefunction.
  • Projection Method: Constructs a projection operator $P$ to explicitly remove forbidden states. The Hamiltonian is transformed to $H_{proj} = (1-P)H(1-P) + \lambda P$, where $\lambda$ is a large energy shift pushing forbidden states out of the physical spectrum.
• Interactions:
  • $n$-$n$ interaction: Minnesota potential in the $s$-wave.
  • Core-$n$ interaction: Woods-Saxon potential (with spin-orbit) tuned to reproduce $^{21}$C scattering data and bound states.
  • Three-body force: A local Gaussian force adjusted iteratively to reproduce the experimental two-neutron separation energy ($S_{2n} = 0.5$ MeV).
• Technical Developments:
  • Implementation of Raynal-Revai transformations to switch between T-basis and Y-basis for interaction matrix elements.
  • Derivation of a general expression for Electric Multipole Strengths ($B(E\lambda)$) in the HH formalism.
  • Development of algorithms to automate the adjustment of the three-body force strength.
  • Implementation of MPI/OpenMP parallelization and model space truncations (cutoffs on $n_{max}$ and $l_{max}$) to reduce computational cost.

3. Key Contributions

1. Comparative Analysis of Pauli Removal: A rigorous comparison between the SUSY and Projection methods for a system with multiple forbidden states ($^{22}$C).
2. New Code (hyperboromir): A robust, parallelized code capable of handling three-body halo nuclei with consistent bound/scattering state treatment.
3. Efficiency Optimization: Identification of specific truncations (specifically $n_{max}$) that reduce computational time by ~20% without sacrificing accuracy.
4. Theoretical Derivations: General formulas for electric multipole operators and automated force adjustment procedures within the HH-R-matrix framework.

4. Results

The study focuses on the ground state ($0^+$) and dipole response ($B(E1)$) of $^{22}$C.

• Convergence: Both bound and scattering states converge at $K_{max} \approx 40$ and $\rho_{max} = 60$ fm.
• Method Comparison ($^{22}$C Properties):
  • Binding Energy & Radius: Both methods yield similar binding energies ($\sim -0.5$ MeV) and root-mean-square hyperradii ($\sim 8.8$ fm) at convergence.
  • Wavefunction Structure: Significant differences emerge in the partial-wave decomposition.
    • The Projection method correctly enforces the Pauli principle, resulting in wavefunctions with extra nodes (anti-symmetrization with the core). It shows a dominant $d$-wave component in the T-basis.
    • The SUSY method populates higher-$l$ partial waves due to the artificial short-range repulsion, leading to a larger $s$-wave component and a slightly larger predicted radius.
  • Dalitz Plots: The probability distributions differ drastically. Projection favors a "cigar" or "di-neutron" mix, while SUSY is dominated by a "di-neutron" configuration.
  • Dipole Strength ($B(E1)$): The magnitude of the dipole strength differs by a factor of $\sim 2$ between the two methods. The Projection method yields a suppressed strength due to the correct occupation of partial waves.
• Conclusion on Accuracy: The Projection method is deemed more accurate. The SUSY method, while computationally simpler, fails to correctly anti-symmetrize the wavefunction in systems with complex shell structures (like the $sd$-shell in $^{22}$C), leading to incorrect structural predictions.
• Computational Efficiency: Truncating the number of channels using $n_{max}$ (limiting the number of nodes in the hyperradial function) reduces the dimension of the Hamiltonian significantly, achieving a 20% reduction in CPU time while maintaining convergence.

5. Significance

• Validation of Formalism: The paper establishes that for medium-mass halo nuclei with multiple forbidden states, the explicit Projection method is superior to the SUSY transformation. This suggests that previous predictions for halo nuclei relying solely on SUSY may need re-evaluation.
• Uncertainty Quantification: By demonstrating efficient truncations and parallelization, the authors enable the running of large statistical ensembles required for robust uncertainty quantification in few-body models.
• Future Applications: The developed formalism and code provide a foundation for studying more complex systems, such as three-body charged systems and other observables beyond dipole strength, facilitating the interpretation of data from next-generation facilities like FRIB.

In summary, this work bridges the gap between accurate shell-model constraints and few-body dynamics, providing a computationally efficient and physically rigorous tool for predicting the properties of exotic halo nuclei.