Lattice Quantum Chromodynamics for Quantum Simulations

技術要約：量子シミュレーションのための格子量子色力学

問題提起
量子色力学（QCD）は強い相互作用の基本理論であるが、その非摂動領域（ハドロン構造や高密度物質の理解に不可欠な領域）は計算上極めて困難である。標準的な格子QCDのアプローチは、ユークリッド経路積分と重要サンプリングに依存しているが、これらは実時間ダイナミクスや$\theta$角が含まれる場合のような「符号問題」が存在する状況下では機能しない。ハミルトニアン定式化は実時間シミュレーションへの道筋を提供するが、格子体積とともにヒルベルト空間が指数関数的に増大するため、大規模な系における古典的なシミュレーションは困難となる。したがって、ゲージおよびフェルミオンの自由度のエンコーディング、物理的（ゲージ不変な）状態の構築、および効率的な時間発展アルゴリズムの開発に特化した、量子コンピュータへの格子ゲージ理論のマッピング手法が切実に求められている。

手法
著者らは、$N_c=3$（QCD）の2次元および3次元空間における、クォークを含む$SU(N_c)$格子ゲージ理論の量子シミュレーションのための包括的なフレームワークを開発した。本手法は以下の柱に基づいている：

1. ハミルトニアン定式化と離散化：
   本研究では、Kogut-Susskindハミルトニアン定式化を利用している。著者らは、フェルミオンのダブリング問題を解決するために、スタガード・フェルミオンとウィルソン・フェルミオンの両方を検討している。理論は連続的な時間を持つ立方格子上に離散化されている。ハミルトニアンは、電場項（$H_E$）、磁場項（$H_B$）、$\theta$項（$H_\theta$）、運動項（$H_{kin}$）、および質量項（$H_{mass}$）に分解される。$\theta$項は、トポロジカル電荷密度への結合として含まれており、ハミルトニアンに複素係数を導入する。

2. 表現基底と物理ヒルベルト空間：
   中心的な貢献は、**表現基底（representation basis）**における物理ヒルベルト空間の定式化である。群の要素を直接量子ビットにマッピングする代わりに、リンク上のゲージ自由度とサイト上のフェルミオン自由度は、$SU(N_c)$の既約表現（irreps）を用いてエンコードされる。

   • リンク・ヒルベルト空間： リンクは状態 $|(R, r_R, r_L)\rangle$ によって記述される。ここで $R$ は既約表現であり、$r_{R,L}$ は右側および左側ハーフリンクの基底ベクトルである。この基底は二次カシミール演算子を対角化する。
   • サイト・ヒルベルト空間： フェルミオンは占有数基底によって扱われ、その後、クレブシュ・ゴルダン係数（CGC）を用いて表現基底へと変換される。これにより、粒子および反粒子の状態を$SU(N_c)$の下での変換特性によって分類することが可能になる。
   • ゲージ不変性： 物理状態は、格子サイトで出会うすべての既約表現（リンクおよびフェルミオン）のテンソル積によって形成されるシングレット（自明な表現）として構成される。著者らは、これらのシングレット状態を曖昧さなく構成するために必要な「サイト順序（site order）」の規約を明示的に定義し、ガウスの法則の制約が満たされることを保証している。

3. 打ち切りとエンコーディング：
   $SU(N_c)$のヒルベルト空間は無限次元であるため、著者らは打ち切りスキームを実装している。これには、ゲージリンクに対する最高ウェイト（$n_1$）または二次カシミール固有値（$C_2$）のカットオフ、およびフェルミオンのカラー占有数の制限が含まれる。物理状態は、特定の既約表現とシングレット多重度指数を表す**局所エンコーディング（local encoding）**スキームを用いて量子ビットにエンコードされる。この手法は、大きなヒルベルト空間に対する直接的なエンコーディングよりも量子ビット効率が高いことが指摘されている。

4. 量子シミュレーション・アルゴリズム：
   著者らは、時間発展のための**トロッター化（Trotterization）**アルゴリズムを提案している。時間発展演算子 $e^{-iHt}$ は、分解されたハミルトニアン項に対応するユニタリ演算子の積として近似される。

   • ゲート構成： プラケットおよびリンク演算子の発展は、ギブス回転（2レベル・ユニタリ）および位相シフトへとマップされる。著者らは、これらの操作のための量子回路の構築を詳述しており、フェルミオンの反交換関係に必要なジョルダン・ウィグナー・ストリング（位相反転）の扱いについても述べている。
   • 行列要素： 本論文は、すべてのハミルトニアン演算子の行列要素を、CGCを含む「サイト因子」の積として、物理基底において明示的な公式で提供している。これらの公式は、量子ゲートに必要な回転角や位相を計算するために不可欠である。

主要な貢献

• **$\theta$角を伴う$SU(3)$のフレームワーク：** 本論文は、量子コンピュータ上でのハミルトニアン・アプローチを用いた、3次元空間における$\theta$角を伴う格子QCDの最初の定式化およびシミュレーションを提示している。
• 表現基底の定式化： 既約表現とシングレットCGCを用いた、物理ヒルベルト空間の厳密な構成を提供している。これは、群要素基底に代わるスケーラブルな選択肢となる。
• アルゴリズムの詳細： トロッター化された発展のための具体的な量子回路、および複雑なハミルトニアン項（$\theta$項やウィルソン運動項など）を基本ゲートへ分解する詳細を記述している。
• 行列要素の計算： 著者らは、プラケット、$\theta$、およびフェルミオン運動項の行列要素を導出し、表にまとめた。これらのデータは、シミュレーション回路の実装に必要となる。

結果
著者らは、小規模な格子（最大$2\times2$周期格子および$2\times2\times2$開境界立方体）を用い、最大32量子ビットの古典シミュレータ（Qiskit AerおよびCUDA-Q）上でノイズレス・シミュレーションを実施した。シミュレーションにより以下のことが示された：

• ヒルベルト空間の打ち切り効果： ゲージおよびフェルミオンの異なる打ち切りが、平均電場（$\langle \Pi^2 \rangle$）やカラー占有数（$\langle N_q \rangle$）などの観測量に与える影響の分析。結果は、強結合が高占有状態を抑制することを示しており、スタガード・フェルミオンに対しては控えめな打ち切りでも効果的である一方、ウィルソン・フェルミオンは打ち切りに対してより高い感度を示すことを示した。
• パラメータ依存性： 観測量の発展を、フェルミオン質量およびウィルソンパラメータ $r$ の関数として調査した。質量が重いほど、フェルミオンに富んだ状態への振動が抑制されることが判明した。
• $\theta$角のダイナミクス： $\theta$角を変化させたシミュレーションにより、強結合領域における非周期的挙動が明らかになり、正および負の$\theta$で異なるダイナミクスが観察された。
• ハドロン状態とストリング・ダイナミクス： 「荷電パイ中間子」状態の準備、およびストリング状態（ループおよびフェルミ・ストリング）の発展を実証した。著者らは、フラックスチューブがフェルミオン・反フェルミオン対へと崩壊するストリング破断（string breaking）のダイナミクスを観察し、破断していない構成と破断した構成の確率を追跡した。
• バリオン化学ポテンシャル： 単一の立方体を用いた厳密対角化により、バリオンを添加する際のエネルギーコストを計算した。これは、高密度物質の研究に向けたプルーフ・オブ・コンセプト（概念実証）として機能する。

意義および主張
本論文は、基礎物理学のための実用的な量子シミュレーションの開発における基礎的な貢献として位置づけられている。著者らは、自らの研究が、ハミルトニアン定式化における格子QCDの量子シミュレーションを設計するための必要な「材料」を提供すると主張している。

• 概念実証： 主要な意義は、$SU(3)$格子ゲージ理論の複雑な機構（$\theta$角およびフェルミオン物質を含む）が、量子ハードウェア（ここではシミュレートされたもの）上でエンコードされ、発展させられることを示した点にある。
• スケーラビリティへの道筋： 局所エンコーディングと表現基底は、より大きなシミュレーションへの実行可能な道筋として提示されており、グローバルなエンコーディングよりも少ない量子ビットを必要とする可能性がある。
• 今後の展望： 著者らは、現在の結果は探索的なものであり、格子サイズ、強結合領域、および用いられた積極的な打ち切りを考慮すると、まだ真の物理的な主張を行う段階には至っていないと謙虚に述べている。彼らは、連続極限に到達し、物理的に関連のある領域をシミュレートするために、改良されたハミルトニアン（カウンター項）、より優れた状態準備アルゴリズム、およびフォールトトレラントなハードウェアが必要であると指摘している。本研究は、将来的な高密度物質（例：中性子星）やアクシオンQCDのダイナミクスをシミュレートするための踏み台として位置づけられている。