GroupFunctions.jl: computing individual entries of the irreducible representations of the unitary group U(d)

GroupFunctions.jl に関する技術要約

問題提起
ユニタリ群 $U(d)$ の表現は、量子光学や代数的組合せ論といった様々な物理学および数学の分野において極めて重要です。これらの分野における中心的な計算課題は、既約表現（irreps）の個々の行列要素、すなわち「群関数」を評価することです。これらの要素は $D^{(\lambda)}_{\text{out,init}}(U) = \langle \text{out} | D^{(\lambda)}(U) | \text{init} \rangle$ と表記され、線形光学ネットワークにおける光子の遷移振幅の計算から、量子デバイスの特性評価、対称性の解析に至るまで、幅広いタスクで必要とされます。

多くの用途では数値的な推定で十分な場合もありますが、遷移振幅や指標（character）の厳密な記号表現を導出する必要がある文脈では、精密な記号計算が求められます。既存のソフトウェアソリューションは、関連してはいるものの異なる問題に対処しています。例えば、クレブシュ・ゴルダン係数を計算するもの（SUNRepresentations.jl）、群関数の間接的な数値アクセスを提供するもの（RepLAB）、あるいは群上の記号積分に焦点を当てたもの（IntegrateUnitary.jl）などがあります。量子光学に特化したライブラリ（BosonSampling.jl や The Walrus など）は、通常、デバイスを数値的にモデル化したり、パーマネントのような特定の振幅を計算したりしますが、個々の表現行列の成分を記号的に計算するように設計されたものはありません。さらに、単一の成分のみが必要な場合に、表現行列全体を組み立てることはしばしば計算効率が悪くなります。

手法
本論文では、個々の行列要素を記号的または数値的に計算するために設計された Julia ライブラリである GroupFunctions.jl を紹介します。このライブラリは、ゲルファント・ツェトリン（GT）パターンによって列挙されるキャリア空間の基底内で動作します。

中核となるアルゴリズム的アプローチは、行列要素を入力行列 $U$ の成分の単項式の和として表す Grabmeier-Kerber 公式に基づいています。これらの単項式は、既約表現および特定の基底状態によって重み付けされます。実装においては、以下のプロセスによって最適化が行われています：

1. 二重剰余類（Double Cosets）の列挙： 対称群 $S_N$（ここで $N$ はヤング図形のボックスの数）のすべての置換を反復する代わりに、同一の単項式に寄与する置型をグループ化します。
2. 効率的な総和： 二重剰余類を構築し、その代表的な置換要素を反復することで、素朴な置換ごとのアプローチと比較して計算オーバーヘッドを大幅に削減します。
3. 内部表現： 基底状態は内部的に半標準ヤング表（semistandard Young tableaux）として表現されますが、ユーザーインターフェースからは GTPattern データ構造を介して等価な GT パターンが公開されます。

本ライブラリは、記号的な入力行列（SymEngine.Basic を使用）をサポートしており、GT パターンを占有数ケット（occupation-number kets）に変換したり、一般的な量子光学のパラメータ化から入力ユニタリを構成したり、関連するシュー関数を計算したりするユーティリティ関数を提供します。

主な貢献

• 個々の成分の記号計算： 作者の知る限り、GroupFunctions.jl は $U(d)$ 表現行列の個々の成分を記号的に計算できる最初のパッケージです。
• 統一されたフレームワーク： 長さが最大 $d$ の整数分割によって指定される既約表現に対して、記号的および数値的な入力を扱う統一された手法を提供します。
• 量子光学への特化： ボゾンサンプリングのための遷移振幅の評価や、対称な既約表現（完全に区別不可能なボゾンをモデル化するもの）の扱いなど、量子光学のシナリオに特化したツールが含まれています。
• エクスポート機能： 結果を Mathematica と互換性のある形式でエクスポートでき、さらなる記号操作を容易にします。

結果と性能
本論文では、$SU(2)$ 群関数（Wigner D-関数への簡約）や、7個のボゾンを含む $U(10)$ 対称既約表現の計算を含む例を通じて、ライブラリの能力を実証しています。

• 性能： AMD Ryzen 7 PRO 4750U CPU 上で、特定の $U(10)$ 行列要素の記号計算は、コンパイル後のウォームアップなしで、キャッシュされた結果を再利用せずに約4秒で完了しました。
• スケーラビリティ： 計算コストは、群の次元 $m$、ゲルファント・ツェトリン空間の次元 $d_\lambda(m)$、および対称群の既約表現の次元 $f^\lambda$（標準ヤング表の数）という3つの要因によって支配されます。著者らは、基礎となる作業には $N!$ に近づく置換が含まれる可能性があるものの、実用上のボトルネックは通常、末尾のゼロによる $m$ の増加ではなく、$N$ または $f^\lambda$ の増大であると指摘しています。
• 応用： 本ライブラリは、ランダムベンチマーキングにおける $SU(d)$ 群の指標を数値的に評価するために既に活用されています。

意義
GroupFunctions.jl の意義は、個々の表現行列の成分に対して厳密な記号解を提供できる能力にあります。これは、既存のソフトウェアエコシステムではこれまで利用できなかった機能です。表現行列全体を組み立てる必要を避けることで、本ライブラリは、代数的組合せ論や量子物理学で必要とされる特定の計算に対して、より効率的な経路を提供します。著者らは、本ツールを、純粋な数値的光学シミュレータと、直接的な記号的成分計算機能を欠いた汎用的な群論ライブラリとの間のギャップを埋める、遷移振幅の厳密な解析解や指標の評価を必要とする研究者のための専門的なリソースとして位置付けています。