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

SU3HOB-cgvcs に関する技術要約

問題提起
本論文は、一般的なタルミ・モシンスキー調和振動子ブラケット（HOB）の計算における課題に対処している。これらのブラケットは、任意の質量比パラメータ $d$ によってパラメータ化された、1組のヤコビ座標における2粒子振動子積状態と、回転したセットとの間の変換を記述するものである。HOBは、並進不変な数体および ab initio 無核殻模型行列要素を構成するための基本的な構成要素である。既存の要件では、多くの2振動子ブロック（全振動子量子数 $E$ および全軌道角運動量 $L$ で定義される）に対し、広範な質量比パラメータ $d$ の範囲で、コンパクトかつ数値的に堅牢な評価スキームが求められている。

手法
著者らは、Kalinauskasら [1] によって以前に再定式化された手法を実装している。この手法は、HOBを完全に $SU(3)$群の言語を用いて表現するものである。核心となる手法は、一般のブラケットを、小さなウィグナー$d$関数の単一の和へと還元することであり、その和は$U(6) \supset U(3) \times U(2)$の鎖に属する$SU(3) \supset SO(3)$アイソファクター（$SO(3)$ 縮退ウィグナー係数）の積によって重み付けされている。

実装の主な側面は以下の通りである：

1. 変数の分離： アイソファクターは質量比パラメータ $d$ に依存しない。したがって、「アイソファクターの塔（係数の行列）」は、各 $(E, L)$ ブロックごとに一度だけ計算され、任意の数の $d$ 値に対して再利用される。この「一度準備して、何度も評価する」という因数分解により、要素ごとの評価と比較して効率が大幅に向上する。
2. 外部ライブラリの統合： パッケージ SU3HOB-cgvcs は、ベクトル・コヒーレント・ステート（VCS）法を用いて結合係数を計算する確立された外部ライブラリ Su3cgvcs [2] から、必要な $SU(3) \supset SO(3)$ アイソファクターを取得する。
3. 慣習の架け橋： 技術的な重要な貢献は、Su3cgvcs で使用されている VCS の慣習と、物理的な HOB に必要な調和振動子（モシンスキー）の慣習との間の架け橋となることである。これには以下が含まれる：
   • 位相補正： $n_i = (e_i - l_i)/2$ （ここで $e_i$ は径方向量子数）を表す状態依存の符号因子 $(-1)^{n_1+n_2}$ を適用する。この位相を省略すると行列の直交性は保たれるが、個々の物理的ブラケットの符号が誤ってしまう。
   • スカラー因子の取り扱い： Su3cgvcs の汎用的な結合ルーチンは、スカラー $(0,0)$ 表現（$e_1=0$ または $e_2=0$ の場合に発生）に対してガードされていない。本パッケージは、これらの自明なケースを解析的に処理し、適切な位相因子を設定した上で、縮退した固有値分解の問題を避けるためにライブラリの呼び出しをスキップする。

主要な貢献

• ソフトウェアパッケージ： 任意の質量比パラメータに対する一般的なタルミ・モシンスキー・ブラケットを計算する Fortran 2008 パッケージ SU3HOB-cgvcs のリリース。
• インターフェースの実装： 一般的な $SU(3)$ 結合ライブラリ（Su3cgvcs）の上に、$SU(3)$ スキームによる HOB コードを構築できることの成功した実証。
• 検証プロトコル： 結果として得られる行列の直交性（$HOB \cdot HOB^T = I$）と、独立した参照コードである Kamuntavičiusら [3] の HOTB ツールとの要素ごとの一致の両方をチェックする厳格な検証戦略。

結果
本パッケージは、Kamuntavičiusら [3] の一般タルミ・モシンスキー・ブラケット・コード、および Buck と Merchant [5] の閉形式の式を用いて検証された。

• 数値精度： 計算されたブラケットは、テストされたすべての $(E, L)$ ブロックおよび質量比 $d$ に対して、参照値とマシン精度（$\sim 10^{-13}$）で一致している。
• 直交性： 直交性の誤差は、テストされた範囲（$E \le 12$）においてマシン精度レベル（$\sim 10^{-13}$）に留まっている。
• 比較限界： $E=12$ において、偏差は $\sim 10^{-11}$ まで上昇する。著者らは、これを $SU(3)$ 手法の失敗ではなく、参照コードの生の二項係数テーブルの倍精度浮動小数点数の限界によるものとしている。
• パフォーマンス： フルセットの $L$ ブロックに対する「準備＋評価」の時間は、単一コア上で $E=4$ の場合の $7.3 \times 10^{-5}$ 秒から $E=12$ の場合の $9.2 \times 10^{-3}$ 秒までの範囲で、効率的にスケールする。

意義
本論文は、SU3HOB-cgvcs が $SU(3)$ スキームによる HOB へのライブラリベースの経路を提供していることを主張しており、既存の、個別に検証された結合ライブラリを使用して、自己完結型の $SU(3)$ スキーム・コードを構築できることを示している。このアプローチは、対称適応基底のために既に Su3cgvcs を利用しているコードにとって便利なインターフェースを提供する。さらに、独立して開発されたライブラリを用いてブラケットを構築することは、基礎となる $SU(3)$ 再定式化手法の厳格なチェックとして機能し、それが従来の座標空間での評価と等価であることを裏付けている。