⚛️ nuclear theory

SU3HOB-cgvcs: harmonic-oscillator brackets in the SU(3) basis with isofactors from an external SU(3) $\,\supset\,$ SO(3) Clebsch--Gordan library

This paper introduces SU3HOB-cgvcs, a Fortran 2008 package that computes general Talmi–Moshinsky harmonic-oscillator transformation brackets in the SU(3) basis by leveraging isofactors from the external Su3cgvcs library, while implementing necessary phase corrections and analytical fixes to achieve machine-precision agreement with established codes.

Original authors: Augustinas Stepšys, Saulius Mickevičius, Darius Germanas

Published 2026-08-05

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Augustinas Stepšys, Saulius Mickevičius, Darius Germanas

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

The Cosmic Dance of Tiny Particles

Imagine trying to describe a dance floor where millions of tiny partners are spinning, jumping, and swapping places all at once. In the world of nuclear physics, this is exactly what happens inside an atom's nucleus. Scientists use a mathematical tool called the "harmonic oscillator" to model how these particles move, kind of like imagining them attached to invisible springs that bounce them around. But here's the tricky part: when you have two particles dancing together, it's often easier to describe their motion from the perspective of the whole pair (like watching a couple waltz) rather than looking at each dancer individually.

To switch between these two viewpoints—looking at the individual dancers versus the couple as a whole—physicists need a special set of instructions called "transformation brackets." Think of these as a universal translator that converts the language of "individual steps" into the language of "group choreography." These translations are crucial for building accurate models of atomic nuclei, which helps us understand everything from how stars burn to how we might one day harness nuclear energy. However, calculating these translations is notoriously difficult, especially when the particles have different masses or when you need to do the math for thousands of different scenarios. It's like trying to translate a novel into a different language while the author is still writing it, and you need to do it perfectly every single time.

The New Translator: SU3HOB-cgvcs

This paper introduces a new, highly efficient software package called SU3HOB-cgvcs that acts as a super-smart translator for these nuclear dance moves. The authors, a team from Lithuania, have built a tool that calculates these complex transformation brackets with incredible speed and precision. Instead of reinventing the wheel, they decided to build their translator on top of an existing, well-tested library of "coupling coefficients" (which are like a dictionary of how different quantum states fit together).

Here is the clever trick they used: They realized that the math for these brackets could be broken down into two parts. The first part is a set of "isofactors"—these are like the static rules of the dance that don't change no matter how heavy or light the particles are. The second part involves a specific parameter called the "mass-ratio" (denoted as d), which changes depending on the specific particles involved. The authors' method calculates the static rules just once and stores them. Then, whenever they need to solve a problem with a new mass ratio, they simply plug that new number into a quick formula. This "prepare once, evaluate many" approach is like baking a giant cake base once and then just adding different frostings later, rather than baking a whole new cake for every flavor.

The software works by taking these pre-calculated rules and combining them with a specific mathematical function (called a Wigner d-function) that handles the mass-ratio changes. To make sure their new translator speaks the same language as the old, trusted ones, the authors had to fix a small "accent" issue. The library they borrowed from uses a slightly different sign convention (a positive or negative sign) based on the "radial quantum number" of the state. The authors discovered they needed to flip these signs in a very specific way to match the physical reality of the particles. They also had to write a special rule for a "scalar" case (where one of the particles has zero energy), which the borrowed library didn't handle automatically.

The results are impressive. The team tested their new code against a gold-standard reference program and found that the numbers matched almost perfectly, with errors so tiny they are only about 0.0000000000001 (or $10^{-13}$ ). This level of precision is what you'd expect from a computer doing double-precision math. They also checked that the math holds together logically (orthonormality), ensuring that the total probability of all possible dance moves adds up to exactly one. The paper demonstrates that by combining a specialized method with a general-purpose library, they created a robust, fast, and accurate tool that can be used by other scientists to build better models of the atomic nucleus. It's a successful bridge between two different ways of doing the math, proving that you can build a specialized, high-speed engine using a reliable, general-purpose chassis.

Technical Summary of SU3HOB-cgvcs

Problem Statement
The paper addresses the computational challenge of evaluating general Talmi–Moshinsky harmonic-oscillator brackets (HOBs). These brackets describe the transformation between two-particle oscillator product states in one set of Jacobi coordinates and a rotated set, parameterized by an arbitrary mass-ratio parameter $d$ . HOBs are fundamental building blocks for constructing translationally invariant few-body and ab initio no-core shell-model matrix elements. Existing requirements demand a compact, numerically robust evaluation scheme capable of handling many two-oscillator blocks (defined by total oscillator quanta $E$ and total orbital angular momentum $L$ ) across a range of mass-ratio parameters $d$ .

Methodology
The authors implement the method previously reformulated by Kalinauskas et al. [1], which expresses HOBs entirely within the language of the $SU(3)$ group. The core methodology reduces the general bracket to a single sum over small Wigner $d$ -functions, weighted by products of $SU(3) \supset SO(3)$ isofactors (SO(3)-reduced Wigner coefficients) belonging to the chain $U(6) \supset U(3) \times U(2)$ .

Key aspects of the implementation include:

Separation of Variables: The isofactors are independent of the mass-ratio parameter $d$ . Consequently, the "isofactor towers" (the matrix of coefficients) are computed once per $(E, L)$ block and reused for any number of $d$ values. This "prepare once, evaluate many" factorization significantly improves efficiency compared to element-wise evaluation.
External Library Integration: The package, SU3HOB-cgvcs, obtains the necessary $SU(3) \supset SO(3)$ isofactors from the established external library Su3cgvcs [2], which calculates coupling coefficients using the vector-coherent-state (VCS) method.
Convention Bridging: A critical technical contribution is the bridging of the VCS convention used by Su3cgvcs to the harmonic-oscillator (Moshinsky) convention required for physical HOBs. This involves:
- Phase Correction: Applying a state-dependent sign factor, $(-1)^{n_1+n_2}$ , where $n_i = (e_i - l_i)/2$ represents the radial quantum number. While omitting this phase preserves matrix orthonormality, it yields incorrect signs for individual physical brackets.
- Scalar Factor Handling: The general-purpose coupling routines in Su3cgvcs are not guarded against scalar $(0,0)$ representations (occurring when $e_1=0$ or $e_2=0$ ). The package analytically handles these trivial cases, setting the isofactor to the appropriate phase factor and skipping the library call to avoid degenerate eigen-decomposition issues.

Key Contributions

Software Package: The release of SU3HOB-cgvcs, a Fortran 2008 package that computes general Talmi–Moshinsky brackets for arbitrary mass-ratio parameters.
Interface Implementation: A successful demonstration of how an $SU(3)$-scheme HOB code can be built atop a general-purpose $SU(3)$ Clebsch–Gordan library (Su3cgvcs).
Validation Protocol: A rigorous validation strategy that checks both the orthonormality of the resulting matrices ( $HOB \cdot HOBT = I$ ) and element-by-element agreement with an independent reference code (the HOTB tool of Kamuntavičius et al. [3]).

Results
The package was validated against the general Talmi–Moshinsky bracket code of Kamuntavičius et al. [3] and the closed-form expression of Buck and Merchant [5].

Numerical Precision: The computed brackets reproduce the reference values to machine precision ( $\sim 10^{-13}$ ) for all tested $(E, L)$ blocks and mass ratios $d$ .
Orthonormality: The orthonormality error remains at the machine-precision level ( $\sim 10^{-13}$ ) across the tested range ( $E \le 12$ ).
Comparison Limits: For $E=12$ , deviations rise to $\sim 10^{-11}$ ; the authors attribute this to the double-precision floor of the reference code's raw binomial tables, not a failure of the $SU(3)$ method.
Performance: The "prepare + evaluate" time for full sets of $L$ blocks scales efficiently, with times ranging from $7.3 \times 10^{-5}$ s for $E=4$ to $9.2 \times 10^{-3}$ s for $E=12$ on a single core.

Significance
The paper claims that SU3HOB-cgvcs provides a library-based route to the $SU(3)$-scheme HOB, demonstrating that a self-contained $SU(3)$-scheme code can be constructed using existing, separately validated coupling libraries. This approach offers a convenient interface for codes that already utilize Su3cgvcs for symmetry-adapted bases. Furthermore, building the brackets from an independently developed library serves as a stringent check of the underlying $SU(3)$ reformulation method, confirming its equivalence to traditional coordinate-space evaluations.

The Cosmic Dance of Tiny Particles

The New Translator: SU3HOB-cgvcs

Technical Summary of SU3HOB-cgvcs

More like this