⚛️ nuclear theory

Constrained coupled-cluster computations yield negative-parity rotational bands

This paper introduces octupole-constrained coupled-cluster theory to compute negative-parity rotational bands in deformed nuclei, demonstrating through $^{20}$ Ne calculations that while positive-parity band energies are robust, negative-parity band energies are more sensitive to the constrained octupole moment and reference state.

Original authors: Gaute Hagen, Benjamin Bally, Thomas Papenbrock

Published 2026-10-12

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Gaute Hagen, Benjamin Bally, Thomas Papenbrock

Original paper licensed under CC BY 4.0 (http://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 nucleus, a dense cluster of protons and neutrons that behaves less like a solid marble and more like a fluid drop capable of twisting, stretching, and spinning. These tiny systems possess a dual nature: they contain individual particles moving in specific orbits, yet they also move together in large, coordinated patterns. When a nucleus spins, it often forms a rotational band, a sequence of energy levels that rise in a predictable way, much like the steps of a ladder. For many nuclei, these bands are straightforward, but some exhibit a peculiar quirk: they possess a second set of rotational bands that spin in a mirror-image fashion, known as negative parity. These bands are crucial for understanding the fundamental shape of matter, particularly in the search for rare physical phenomena that could explain why the universe is made of matter rather than antimatter. However, calculating these mirror-image states is notoriously difficult for physicists because the standard mathematical tools they use often fail to see them unless the nucleus is already distorted in a very specific way.

A team of researchers has developed a new computational method to force these hidden states into view, allowing them to map out the energy levels of a specific nucleus, neon-20, with greater precision than before. In their work, published in a recent study, the scientists tackled a problem where the nucleus naturally prefers a symmetric shape, hiding the very asymmetry needed to create the negative-parity bands. By introducing a controlled distortion into their calculations and then carefully removing the artificial forces used to create that distortion, they were able to reveal the full spectrum of the nucleus, including the elusive mirror-image bands. Their findings show that while their method successfully reproduces the known energy levels of the standard rotational band, it still struggles slightly to get the exact energy of the mirror-image band right, suggesting that even more complex interactions between the particles are needed to fully describe the nucleus.

The challenge the researchers faced stems from how nuclear shapes are usually calculated. In many cases, a nucleus naturally settles into a shape that breaks its own symmetry, such as a pear shape, which makes it easy to calculate its mirror-image states. But in the case of neon-20, the nucleus naturally prefers a shape that is perfectly symmetric, like a football, which preserves its mirror symmetry. In this symmetric state, the mathematical equations that describe the nucleus cannot spontaneously generate the mirror-image bands. To solve this, the researchers had to manually break the symmetry. They did this by applying a mathematical constraint that forced the nucleus to adopt a pear-like distortion, even though it didn't want to be that way naturally. This is similar to holding a flexible object in a twisted position to study how it would behave if it stayed that way, rather than letting it snap back to its comfortable shape.

Once they forced the nucleus into this distorted state, they used a sophisticated technique called coupled-cluster theory to calculate the energy of the system. This method accounts for the complex interactions between all the protons and neutrons, treating them as a collective whole rather than just a sum of individual parts. However, simply calculating the energy of a distorted nucleus isn't enough, because the real nucleus in nature is not permanently distorted; it fluctuates and rotates. To get the true energy levels, the researchers had to perform a second step called symmetry projection. This process mathematically "un-twists" the forced distortion, averaging out the shape to restore the natural symmetry of the nucleus while keeping the information about the rotational bands. It is a delicate operation: if the initial distortion is too small, the calculation fails to find the mirror-image band; if it is too large, the results become unreliable. The team had to test a wide range of distortions to find the sweet spot where the calculation was both stable and accurate.

The results of these simulations provided a detailed map of the energy levels for neon-20. The researchers found that their method worked exceptionally well for the standard, positive-parity rotational band. The calculated energy levels for this band matched experimental data very closely, and the results were robust, meaning they did not change much regardless of how much distortion was initially applied. This confirmed that the method was sound for describing the familiar behavior of the nucleus. However, the story was different for the negative-parity band. The energy levels for this mirror-image band were more sensitive to the details of the calculation. While the method successfully located the band, the predicted energy was slightly higher than what is observed in experiments, and the results varied more depending on the specific starting conditions of the simulation.

This discrepancy highlights a remaining gap in our understanding of nuclear forces. The researchers noted that the negative-parity band is more sensitive to the specific way the nucleus is modeled and the forces used in the calculation. They suspect that the current models might be missing some subtle interactions between the particles, or that the mathematical approximation used to restore symmetry might need refinement. Despite this, the study represents a significant step forward. It demonstrates that it is possible to compute these difficult mirror-image states even when the nucleus does not naturally break its symmetry, a feat that was previously out of reach for this type of high-precision calculation. By successfully generating these states, the team has opened a new door for studying octupole collectivity, a type of nuclear vibration that is essential for understanding the fundamental symmetries of the universe.

The work also revealed that simply mixing different distorted states together, a technique often used to improve accuracy, did not significantly change the results in this specific case. The calculations showed that the different distorted states were nearly identical in their mathematical behavior, suggesting that the method had already captured the essential physics of the system. The researchers concluded that while their current approach is powerful, future improvements will likely require even larger computational models and more advanced ways of handling the symmetry restoration. For now, the study stands as a proof of concept, showing that by carefully guiding the mathematics through a forced distortion and then restoring the natural order, scientists can uncover the hidden rotational bands of nuclei that were previously invisible to their equations. This approach offers a promising path toward a more complete understanding of the quantum world inside the atom.

Technical Summary: Constrained Coupled-Cluster Computations Yield Negative-Parity Rotational Bands

Problem Statement
Many deformed atomic nuclei exhibit low-lying rotational bands with negative parity. In standard mean-field calculations, such bands typically arise from the restoration of parity when the variational solution spontaneously breaks reflection symmetry. However, in certain nuclei, such as $^{20}\text{Ne}$ , the Hartree–Fock (HF) minimum preserves parity (exhibiting quadrupole deformation but maintaining reflection symmetry). In these cases, negative-parity states cannot be accessed via spontaneous symmetry breaking. Instead, reflection symmetry must be broken explicitly, for instance, by imposing an octupole constraint in the mean-field calculation. While this strategy is established in mean-field theory and projected generator coordinate methods (PGCM), extending it to ab initio single-reference methods like coupled-cluster (CC) theory presents challenges. Specifically, if an octupole constraint is applied only to the reference state, the coupled-cluster singles ( $T_1$ ) amplitudes tend to relax the orbitals back toward the lower-energy, parity-conserving HF minimum, effectively undoing the symmetry breaking.

Methodology
The authors introduce octupole-constrained coupled-cluster theory to address these challenges. The approach involves the following key steps:

Constrained Hamiltonian: An octupole constraint is enforced via a Lagrange multiplier ( $\lambda_3$ ) in the Hamiltonian (Routhian), $H' = H - \lambda_1 Q_{10} - \lambda_3 Q_{30}$ . To prevent the center of mass from shifting due to the laboratory-coordinate definition of the octupole operator, a dipole constraint ( $\lambda_1$ ) is simultaneously enforced to fix the center of mass at the origin.
CCSD Implementation: The authors employ the coupled-cluster singles-and-doubles (CCSD) approximation. The right eigenstate is expanded as $|\Psi\rangle = e^T |\Phi_0\rangle$ and the left eigenstate as $\langle \tilde{\Psi}| = \langle \Phi_0|(1 + \Lambda)e^{-T}$ . The amplitudes $T$ and $\Lambda$ are determined by solving the bivariational equations for the constrained Routhian. Crucially, the octupole constraint is maintained throughout the correlated calculation to prevent the $T_1$ operator from relaxing the system back to a parity-conserving state.
Symmetry Projection: After variation, the intrinsic states (which have broken parity) are projected onto good angular momentum ( $J$ ) and parity ( $\pi$ ) quantum numbers. The authors utilize the disentangled projected coupled-cluster formulation (PCCSD), employing a projected linear-response approach where both the right and left states enter the energy functional. The projection operators are implemented using a unitary transformation that includes parity projection via an exponential form.
Model Space and Interactions: Calculations are performed for $^{20}\text{Ne}$ using chiral effective field theory interactions (two- and three-nucleon forces, specifically the 1.8/2.0 EM interaction). The model space includes spherical harmonic-oscillator states up to $e_{\text{max}} = 8$ with an oscillator frequency of $\hbar\omega = 16$ MeV.

Key Contributions

Extension of Constrained Techniques: The paper extends the technique of enforcing symmetry breaking via constraints from mean-field theory to coupled-cluster theory.
Octupole-Constrained CCSD: The authors demonstrate that maintaining the octupole constraint within the CCSD equations is necessary to preserve the broken symmetry in the correlated state, preventing the relaxation of orbitals toward the parity-conserving minimum.
Dipole Constraint Necessity: The work highlights the necessity of a dipole constraint to fix the center of mass when using laboratory-coordinate octupole operators, preventing spurious dipole moments from collapsing the octupole deformation.
Projection After Variation: The study applies combined angular-momentum and parity projection to constrained CCSD states, utilizing a disentangled expansion truncated at the two-particle–two-hole level ( $W_2$ ) for the rotated states.

Results
The authors computed the rotational bands of $^{20}\text{Ne}$ using four different Hartree–Fock references constrained to octupole moments ( $Q_{30}^{\text{HF}}$ ) of 75, 100, 125, and 150 fm $^3$ .

Positive-Parity Band: The energies of the positive-parity ground-state rotational band ( $0^+, 2^+, 4^+$ ) are close to experimental data. These energies depend only weakly on the constrained octupole moment and the specific HF reference state. The ratio $R_{42} = E(4^+)/E(2^+)$ is approximately 3.11, closer to the rigid-rotor value (10/3) than the experimental value (2.61), a result consistent with previous projected mean-field studies.
Negative-Parity Band: In contrast, the negative-parity band ( $1^-, 3^-, 5^-$ ) shows significant sensitivity to the auxiliary quantities (the octupole constraint and the HF reference). The energy minima for negative-parity states occur at higher octupole moments than for positive-parity states. The calculated excitation energy of the $1^-$ band head is approximately 6.93 MeV (relative to the $0^+$ ground state), which is about 1.1 MeV higher than the experimental value. The rotational constant is also too large (moment of inertia too small) compared to experiment.
Dependence Analysis: The study reveals that while positive-parity states are robust against variations in the constraint, negative-parity states vary by more than 1 MeV in excitation energy across the different references.
Mixing of States: Attempts to mix projected constrained CCSD states (analogous to PGCM) along a single trajectory yielded negligible improvements (less than 0.2 MeV) because the projected states were found to be nearly linearly dependent.

Significance and Claims
The paper claims that constrained coupled-cluster theory provides a viable method to access collective negative-parity states in nuclei where reflection asymmetry does not emerge spontaneously at the mean-field minimum. The method successfully generates a low-lying negative-parity band using chiral two- and three-nucleon forces.

The authors note that while the positive-parity band is well-described, the negative-parity band energies are too high and exhibit stronger dependence on the constraint and reference state. They attribute potential discrepancies to:

Limitations of the current nuclear Hamiltonian.
Lack of full convergence with respect to the model space size ( $e_{\text{max}}$ ).
Shortcomings of the "projection after variation" approach (as opposed to "variation after projection").
The omission of triaxial deformations (both quadrupole and octupole), which recent PGCM calculations suggest are important for $^{20}\text{Ne}$ .

The work concludes that constrained coupled-cluster theory makes octupole collective states accessible in a wider class of nuclei, though further improvements in model spaces, projection techniques, and the inclusion of triaxial degrees of freedom are required for quantitative precision.

More like this