Exact chiral symmetry with quantum signal processing

技術要約：量子シグナルプロセッシングを用いた厳密なカイラル対称性

問題提起
格子量子色力学（QCD）における非摂動的、非平衡的な観測量の量子シミュレーションは、フェルミオンのダブリング問題と、ハミルトニアン定式化におけるカイラル対称性の保持に関して大きな課題に直面している。ユークリッド時空におけるギンスパルク・ウィルソン（GW）関係式やオーバーラップ演算子はエレガントな解決策を提供するが、ハミルトニアン定式化における標準的な対応物となるオーバーラップ・フェルミオンは欠けていた。ドメインウォール・フェルミオンのような既存のアプローチはハミルトニアンの設定に適合するが、明示的な余剰次元を必要とするため、量子ビットのコストが増大する。一方、オーバーラップ・フェルミオンは物理的な空間次元内で定式化されるが、ウィルソン・ディラック・ハミルトニアンの符号関数 $\varepsilon(h_W)$ を含む、高度に非局所的な全対全（all-to-all）相互作用を伴う。これは量子ハードウェア上で効率的に実装することが困難である。本研究が扱う中心的な問題は、メモリ（量子ビット）コストとゲート複雑性のバランスを取りつつ、厳密な（あるいは制御可能な形で破れた）カイラル対称性を維持する、効率的な量子アルゴリズムをいかに構築するかである。

手法
著者らは、オーバーラップ・フェルミオン・ハミルトニアンをシミュレートするための量子シグナルプロセッシング（QSP）アルゴリズムを提案している。その手法は以下の主要なステップで構成される：

1. ハミルトニアンの定式化： 本研究では、オーバーラップ・ハミルトニアン $H_{ov} = \psi^\dagger (\Gamma_0 + \varepsilon(h_W)) \psi$ を利用する。ここで $h_W$ はウィルソン・ディラック単一粒子ハミルトニアンであり、符号関数 $\varepsilon(X) = X/\sqrt{X^\dagger X}$ は次数 $M$ の多項式 $p_M(x)$ によって近似される。
2. ブロック符号化： 著者らは、単一粒子ウィルソン・ハミルトニアン $\hat{h}_W$（ゲージ場を含む）のブロック符号化を、「準備（prepare）」演算子 $P$ と「選択（select）」演算子 $S$ を用いて構築する。これにより、拡張されたヒルベルト空間（アンシラ量子ビットを含む）内にハミルトニアンを符号化する。
3. 量子シグナルプロセッシング（QSP）： 符号関数の近似を実装するために、著者らはQSP（具体的には量子特異値変換）を採用する。これにより、ブロック符号化された演算子に対して、次数 $M$ の多項式 $p_M(U_{\hat{h}_W})$ を適用することが可能となる。多項式は、スペクトルギャップ $[\kappa, 1]$ 内で誤差 $\epsilon_e$ を持つように選択される。
4. 時間発展： オーバーラップ・ハミルトニアンが一度ブロック符号化されると、時間発展演算子 $e^{-iHt}$ はQSPを用いて近似される。これには、ブロック符号化のコストと発展時間の両方に比例するゲート数が必要となる。

主な貢献と結果

• アルゴリズムの構築： 本論文は、制御可能な誤差 $\epsilon_e$ までGW関係式を保持する、オーバーラップ・ハミルトニアンのための具体的なQSPベースのアルゴリズムを提供している。修正されたカイラル演算子 $\hat{\gamma}_5$ は、誤差のオーダー $\epsilon_e$ まで近似ハミルトニアンと可換であることが示されている。
• 複雑度解析：
  • ゲート複雑度： オーバーラップ・ハミルトニアンをブロック符号化するコストは $O(Q \kappa^{-1} \log(1/\epsilon_e))$ でスケールする。ここで $Q$ は格子サイト数（内部自由度を含む）であり、$\kappa$ はスペクトルギャップに関連する。時間発展のための総ゲート複雑度は $O(Q \kappa^{-1} \log(1/\epsilon_e) [Qt + \log(1/\epsilon_t)])$ とスケールする。
  • 量子ビットコスト： アルゴリズムは $O(Q + \log Q)$ 量子ビットを必要とする。これは、余剰次元の広がり $L_5$ を必要とするドメインウォール・フェルミオンの $O(Q L_5)$ 量子ビットと比較して、大幅な削減となっている。
  • 比較： ドメインウォール・フェルミオンは幾何学的な局所性により恩恵を受ける（鈴木・トロッター法などにより、ほぼ線形なコスト $O(Q t \text{polylog})$ が可能）一方で、オーバーラップ定式化は符号関数の近似に伴う非局所性により、より高いゲート深さを招く。しかし、オーバーラップ・アプローチはメモリにおける優れた漸近的スケーリングを提供する。
• QSPの物理的解釈： 著者らは、誤差 $\epsilon_e$ を達成するために必要な多項式の次数 $M$ が $M = O(\kappa^{-1} \log(1/\epsilon_e))$ とスケールすることを示す。彼らはこのスケーリングを、ドメインウォール・フェルミオンにおける余剰次元 $L_5$ のサイズ（誤差が $e^{-c L_5}$ でスケールする）と同一視している。したがって、QSPによる実装は、回路の深さを通じて実質的に「余剰次元を構築」しており、これはオーバーラップ演算子とドメインウォール・フェルミオンの境界理論との間の物理的な対応関係を反映している。

意義と主張
本論文は、厳密なカイラル対称性を持つディラック・フェルミオンの量子シミュレーションは、誤差許容度に対する対数因子分だけ（ウィルソン・ディラック・ハミルトニアンと比較して）「ほぼ無料（nearly free）」であると主張している。

中心的な洞察は、QSPが、オーバーラップ演算子とドメインウォール・フェルミオンの余剰次元との間の既知の対応関係を、量子アルゴリズム的な実現として提供している点にある。そこには明確なトレードオフが存在する：

• ドメインウォール・フェルミオン： 量子ビットコストは高い（$O(Q L_5)$）が、局所性によりゲート深さは低い。
• オーバーラップ・フェルミオン（QSP経由）： 量子ビットコストは低い（$O(Q)$）が、非局所性によりゲート深さは高い（時間発展における最悪のスケールは $O(Q^2)$）。

著者らは、どちらの定式化を選択するかは、ハードウェアの制約（量子ビット数か、それとも回路の深さか）に依存すると結論付けている。また、古典的な格子QCDでは符号関数のために（Zolotarev法などの）有理近似がより効率的であるために用いられるが、これらは標準的なQSPには直接的な類似が存在しないことを指摘しており、これが今後の研究の方向性を示唆している。最後に、スケーリングの結果はゲージ場の符号化方法に対して頑健であるが、定数因子は変化し得ることを明確にしている。