REVIEW 3 major objections 5 minor 7 references
SU3HOB-cgvcs: harmonic-oscillator brackets in the SU(3) basis with isofactors from an external SU(3)$\,\supset\,$SO(3) Clebsch--Gordan library
T0 review · 3 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read A Fortran 2008 package computes Talmi–Moshinsky harmonic-oscillator brackets in an SU(3) basis, taking the needed isofactors from a general-purpose SU(3) Clebsch–Gordan library and bridging conventions with a single state-dependent sign; th
desk verdict Useful SU(3)-scheme HOB package with a new phase bridge to Su3cgvcs; the abstract's unqualified machine-precision claim outruns the actual validation window. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the SU(3)-basis expression of Eq. (1): the bracket is a single sum over J of small Wigner d-functions of the reflection-containing Talmi–Moshinsky matrix, weighted by products of SU(3)⊃SO(3) isofactors of U(6)⊃U(3)×U(2). The load-bearing bridge is Eq. (5), a state-dependent sign (−1)^(n1+n2) that maps the external library's vector-coherent-state convention coefficients to the harmonic-oscillator isofactor convention, plus an analytic override for the (0,0) representation. Because the isofactor towers in C are d-independent, they are calculated once per (E,L) and reused for all mass ratios d—that factorization is what makes the "prepare once, evaluate many" design work.
What would settle it
Take a block beyond the validated range, e.g. E=16, L=8 with d=2, and compare every bracket from this package against an arbitrary-precision direct evaluation of the Talmi–Moshinsky integral (or an independent high-precision HOB code). If the restored-phase brackets do not agree to the package's own precision, the claimed universal sign bridge—and with it the general claim of machine-precision agreement—is falsified. A second check: compare the external library's isofactors for a high (λ,μ) such as (16,0) with those from an independently derived SU(3) CG algorithm.
Extended reading notes
Core claim
The paper's central claim is that the general Talmi–Moshinsky bracket factorizes as H(d)=C D(d) C^T per (E,L) block: a d-independent matrix C of SU(3)⊃SO(3) isofactors, and a d-dependent diagonal-in-J matrix D of small Wigner d-functions of the reflection-containing Talmi–Moshinsky transformation. The isofactors can be read directly from a general-purpose SU(3) Clebsch–Gordan library, provided the user reinserts one sign per state, (−1)^(n1+n2), and supplies the scalar (0,0) factor analytically. With these conventions restored, the computed brackets agree with the coordinate-space Talmi–Moshinsky code of Ref. [3] to machine precision. The paper further notes that the assembled bracket is inv
Load-bearing premise
The bridge relies on the claim that one state-dependent sign, (−1)^(n1+n2), completely converts the external library's coefficient convention to the harmonic-oscillator convention for every (E,L) block; this is verified only up to E=12 because the independent reference code loses numerical reliability above that.
Editorial extensions
If this is right
- For any (E,L) block, all brackets for all mass ratios d are obtained cheaply after one isofactor preparation; this makes wide scans over d practical in few-body calculations.
- The multiplicity-basis invariance means a code using any SU(3) coupling library, regardless of its convention for degenerate α-states, can produce the same physical HOBs without re-orthogonalizing.
- The validated agreement for E≤10 (and near-machine-precision at E=12) supports use of SU(3)-scheme HOBs as a drop-in replacement for coordinate-space Talmi–Moshinsky brackets in existing shell-model workflows.
- The package provides a call-through template for embedding SU(3)-scheme HOBs in codes that already link the external SU(3) coupling library.
Reading between the lines
- The same phase bridge likely applies to other SU(3)⊃SO(3) recoupling coefficients drawn from a vector-coherent-state library, so quantities such as SU(3) Racah factors may need only the same sign fix; this is testable but not claimed in the paper.
- Since the phase depends only on radial parity, the method may tolerate convention changes inside the external library that do not alter that parity; conversely, if a future library version flips the VCS phase convention, the bridge would silently break—a version-pinning concern.
- The restriction to E≤12 is a validation-range limit, not a method limit; extending the comparison with arbitrary-precision reference brackets for E>12 would directly test the universality of the sign mapping.
- The factorization H=CDC^T suggests a route to generalized brackets where each Jacobi coordinate has its own mass ratio (multiple d parameters), since C is unchanged and only D would gain additional structure.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents SU3HOB-cgvcs, a Fortran 2008 package for computing general Talmi–Moshinsky harmonic-oscillator transformation brackets in an SU(3)-coupled basis with an arbitrary mass-ratio parameter d. The method follows the authors' earlier SU(3) reformulation, reducing the bracket to a single sum over Wigner d-functions weighted by products of SU(3)⊃SO(3) isofactors. The new contribution is the interface to the external Su3cgvcs Clebsch–Gordan library, which requires a state-dependent sign fix (Eq. 5) and an analytic treatment of the (0,0) scalar oscillator factor. The package is validated by orthonormality tests and by element-by-element comparison with the independent HOTB code of Kamuntavičius et al. for selected (E,L) blocks up to E=12. The central claim is that the package reproduces the reference brackets to machine precision (~1e-13) and provides a template for building SU(3)-scheme HOB codes on top of a general-purpose SU(3) CG library.
Significance. If the central claim is fully supported, this is a useful contribution to the computational nuclear-structure toolbox. The 'prepare once, evaluate many' factorization is attractive: isofactor towers are d-independent and built once per (E,L) block, while subsequent evaluations only require small Wigner d-functions and matrix products. The use of an external, separately published CG library rather than a new bespoke coefficient generator is a commendable design choice, and the comparison against an independent coordinate-space code is the right kind of validation. The paper also supplies explicit multiplicity checks, a public repository, and a build-and-check workflow, which enhance reproducibility. However, the validity of the phase relation in Eq. (5) is load-bearing for the universal 'bracket-for-bracket' claim, and the reported validation leaves gaps that need to be closed before the claim is fully credible.
major comments (3)
- [Sec. 3.2, Eq. (5)] The phase relation (−1)^{n1+n2} is asserted, not derived or traced to the explicit phase conventions of Refs. [2] and [3]. This is the only nontrivial bridge between the external VCS-based library and the harmonic-oscillator isofactor convention, and it is exactly the kind of step that can fail in a state-dependent way. As the manuscript itself notes, orthonormality is invariant under per-state sign flips, so the orthonormality tests cannot detect an error in Eq. (5). The absolute element-by-element test in Table 2 covers only L=E/2 and even E up to 12; low-L blocks and odd-E blocks are not reported. The universal statement 'makes Eq. (1) reproduce the physical Talmi–Moshinsky bracket exactly' is therefore stronger than the evidence. I recommend either deriving Eq. (5) from the phase conventions of the two libraries or extending the absolute comparison to all L values (and, where feasibl
- [Sec. 6, Table 2] Table 2 shows max|SU3HOB − HOTB| = 1.4×10^{-11} at E=12, while the abstract and summary claim agreement to ~10^{-13}. This discrepancy is attributed to 'the double-precision floor of the reference's raw binomial tables,' but the stated reason is not convincing: the text says those tables are exact only while C(n,m)<2^53 (n≲57), whereas E=12 involves n values far below that threshold. If the reference is still reliable at E=12, then 1.4×10^{-11} is a genuine disagreement with the machine-precision claim; if it is not reliable, the authors need to provide a concrete error analysis (e.g., cancellation in the HOTB summation) rather than the current binomial-overflow argument. This point directly affects the central numerical claim and should be resolved.
- [Sec. 1 and Sec. 6] The introduction states that the computed brackets agree with the reference 'for every (E,L) block and mass ratio d,' but the quantitative validation presented in Table 2 is limited to the representative mid-L blocks L=E/2 for E=2,4,...,12, with d=1 tabulated and d=1/2,2 only described as differing in the last digit. No table or supplementary output is provided for all L, all odd-E blocks, or the full d range. The 'PASS' statement in the Makefile check may cover more, but the manuscript should show or cite the complete validation data. At minimum, the claims in the abstract and introduction should match the range actually demonstrated, or the full sweep should be reported.
minor comments (5)
- [Abstract] Typos: 'These brackets arte constructed' should be 'are constructed'; the phrase 'In this library they are computed there' is redundant.
- [Sec. 1] 'ab initiono-core' should be 'ab initio no-core'.
- [Sec. 3.2, Eq. (5)] The notation n_i = (e_i − l_i)/2 implicitly assumes e_i and l_i have the same parity; this is standard for HO states but should be stated explicitly for readers outside the field.
- [Sec. 6, Table 2] The column layout of Table 2 is hard to read because the entries are packed together; adding explicit column separators or aligning the numbers would improve clarity. Also, the table caption should state the ranges of e1,e2 and l1,l2 included in 'every bracket' for each block.
- [Sec. 3.3] The phrase 'not guarded against it' is informal; suggest a more precise statement, e.g., 'the routine does not handle the degenerate scalar-factor case.'
Circularity Check
No significant circularity: the derivation chain is an interface among a prior SU(3)-method paper [1], an external CG library [2], and an independent coordinate-space benchmark [3]; the phase fix is asserted and only partially validated, but it is not fitted to the benchmark and the central claim does not reduce to its inputs by construction.
full rationale
The claimed derivation chain is: Eq. (1) from Ref. [1] (a method paper by the same group, hence a self-citation); isofactors from the separately authored and separately validated Su3cgvcs library [2] (Bahri, Rowe, Draayer); and validation against the general Talmi–Moshinsky bracket code of Kamuntavičius et al. [3], which shares two authors with the present paper. These citations are real, independent evidence: [1] and [3] are published, separately implemented derivations, [2] is an external library, and the comparison is element-by-element against a coordinate-space formula (Buck–Merchant). The only nontrivial bridge, the convention-phase relation (−1)^{n1+n2} in Eq. (5), is asserted as the difference between the Draayer VCS and HO conventions and is then tested absolutely; it is not a parameter fitted to the benchmark. The paper itself flags the limitation at Sec. 6, Table 2: the absolute bracket-level comparison is restricted to E=2,...,12 and to mid-L blocks, with the remark 'We restrict the quantitative comparison to E≲12: the reference bracket relies on raw factorial/binomial tables that are exact only while C(n,m) < 2^53 (n≲57)'. The paper also explicitly notes in Sec. 3.2 that orthonormality is phase-blind: 'As noted in Sec. 2, the orthonormality of the assembled matrix is insensitive to this phase.' This means the universal validity of Eq. (5) at larger E and low-l/high-n states is an unverified extrapolation, not a demonstrable failure or a circular reduction. Because the load-bearing inputs (prior method, external CG library, independent coordinate-space reference) are not defined in terms of the output brackets, and the validation is an absolute test against an independent implementation, no step in the paper reduces to its own conclusion by construction. Score 1 reflects only the minor self-citation overlap in Refs. [1] and [3], which does not make the central claim circular.
Assumptions & free parameters
assumptions (5)
- domain assumption Eq. (1) from Ref. [1] correctly expresses the general Talmi-Moshinsky bracket in the SU(3) basis.
- domain assumption The Su3cgvcs library returns correct SO(3)-reduced SU(3) Clebsch-Gordan coefficients for all (e1,0) x (e2,0) couplings needed.
- ad hoc to paper The phase relation (−1)^{n1+n2} (Eq. 5) maps the VCS convention to the HO convention universally for all E and L.
- standard math The Racah multiplicity formula in Eq. (3) correctly counts the number of isofactor columns.
- standard math det(TM) = -1, giving the phase factor Δ0^{E2} in Eq. (1).
Cite this review
Pith. "Pith review of SU3HOB-cgvcs: harmonic-oscillator brackets in the SU(3) basis with isofactors from an external SU(3)$\,\supset\,$SO(3) Clebsch--Gordan library." pith.science (2026). https://pith.science/paper/7TARFGSX
@misc{pith2026260802640,
author = {Pith},
title = {Pith review of: SU3HOB-cgvcs: harmonic-oscillator brackets in the SU(3) basis with isofactors from an external SU(3)$\,\supset\,$SO(3) Clebsch--Gordan library},
year = {2026},
howpublished = {\url{https://pith.science/paper/7TARFGSX}},
note = {Machine review of arXiv:2608.02640}
}
abstract
\DeclareRobustCommand{\Su}{\code{Su3cgvcs}} We present SU3HOB-cgvcs, a Fortran~2008 package that computes the general Talmi--Moshinsky harmonic-oscillator transformation brackets, for an arbitrary mass-ratio parameter $d$. These brackets arte constructed in the $\su{SU}(3)$ basis by the method of Kalinauskas et al. . The transformation is reduced to a single sum over Wigner $d$-functions of the reflection-containing Talmi--Moshinsky matrix weighted by products of $\su{SU}(3)\supset\su{SO}(3)$ isofactors of the chain $\su{U}(6)\supset\su{U}(3)\times\su{U}(2)$. SU3HOB-cgvcs obtains these isofactors from the established external Clebsch--Gordan library \Su{} . In this library they are computed there by the vector-coherent-state (VCS) method. We bridge \Su{} formalism to ours and this requires two fixes: reinserting a sign/phase factor that depends on each state's radial-quantum number and handling the trivial $(0,0)$ representation analytically. The package reproduces, bracket for bracket to machine precision ($\sim\!10^{-13}$), the general Talmi--Moshinsky bracket code of Kamuntavi\v{c}ius et al. , and demonstrates that a self-contained SU(3)-scheme HOB code can be built on top of a general-purpose SU(3) Clebsch--Gordan library.
Reference graph
Works this paper leans on
-
[2]
Computer Physics Communications , volume =
Programs for generating Clebsch--Gordan coefficients of. Computer Physics Communications , volume =. 2004 , issn =. doi:https://doi.org/10.1016/j.cpc.2004.01.005 , author =
- [3]
-
[1]
Calculation of harmonic oscillator brackets in SU(3) basis , journal =
Ramutis Kazys Kalinauskas and Augustinas Stepšys and Darius Germanas and Saulius Mickevičius , keywords =. Calculation of harmonic oscillator brackets in SU(3) basis , journal =. 2025 , issn =. doi:https://doi.org/10.1016/j.nuclphysa.2025.123033 , url =
-
[4]
The general harmonic-oscillator brackets: compact expression, symmetries, sums and
Kamuntavi. The general harmonic-oscillator brackets: compact expression, symmetries, sums and. Nucl. Phys. A , volume =
- [5]
-
[6]
Bahri, C. and Rowe, D. J. and Draayer, J. P. , title =. Comput. Phys. Commun. , volume =
-
[7]
Elliott, J. P. , title =. Proc. R. Soc. Lond. A , volume =
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.