Symmetry-adapted qubit encoding with complete active space and Bravyi--Kitaev mapping for quantum chemistry on a quantum computer

本論文は、近似的なZ対称性とBravyi-Kitaevマッピングを統合した完全活性空間（SAE-CAS）による対称適応型量子ビットエンコーディングを紹介するものであり、これにより量子化学シミュレーションにおける量子ビット数と回路の複雑さを大幅に削減し、近未来型およびフォールトトレラントな量子プロセッサの両方において、標準的な手法よりも優れた収束性とリソース効率を実証している。

要約

あなたは、分子がどのように機能するかを理解するために、非常に巨大で複雑なジグソーパズルを解こうとしていると想像してください。量子化学の世界では、このパズルは分子の「電子構造」と呼ばれます。

これを量子コンピュータで解くには、通常、電子が存在しうるすべての場所に、コンピュータの極めて小さな断片（「量子ビット」）を割り当てる必要があります。

問題は、たとえ小さな分子であっても、これには数千個のパズルピース（量子ビット）が必要になり、それらを組み立てるための指示書（回路）があまりに長く、複雑に絡み合ってしまうため、現在のコンピュータでは処理できず、将来のコンピュータでさえ苦労することになるという点です。

この論文は、このパズルを解くための巧妙な新しい方法であるSAE-CASを紹介しています。その仕組みを、簡単な比喩を使って説明しましょう。

1. 「凍結と無視」戦略（CASの部分）

分子を、忙しいオフィスビルだと考えてみてください。

• 凍結コア（Frozen Core）: 地下や最上階には、決して外に出ることもなく、他のフロアとも関わることがない人々が常に満席の状態です。量子論的には、これらは「凍結コア」電子です。彼らは退屈で、予測可能です。
• 仮想軌道（Virtuals）: 屋根裏部屋は完全に空っぽで、おそらく空のままです。これらは「仮想」軌道です。
• 活性空間（Active Space）: 中間のフロアこそが、本当の活動が行われる場所です。人々が動き回り、話し、変化を起こしています。これが「活性空間」です。

従来の方法では、退屈な地下や空っぽの屋根裏を含め、あらゆるフロアに量子ビットを割り当てようとします。SAE-CASはこう言います。「地下や屋根裏は無視しましょう。」私たちは、興味深い化学反応が起きる中間のフロアにだけ、量子ビットを割り当てます。これにより、解くべきパズルのサイズが即座に縮小されます。

2. 「対称性のショートカット」（SAEの部分）

この忙しい中間のフロアの中にも、ルールが存在します。例えば、水分子の場合、左側は右側の鏡像になっています。もし左側で何が起きているかを知っていれば、右側で何が起きているかも自動的に分かります。

通常、コンピュータは両方のサイドを別々に計算するため、時間とリソースを無駄にします。SAE-CASは、数学的な「魔法のトリック」（アフィン・クリフォード変換と呼ばれます）を使用して、鏡のルールがあるため、両側に別々の量子ビットを用意する必要はないということを突き止めます。私たちはパズルを数学的に「半分に折りたたむ」ことができるのです。これにより、さらに多くの量子ビットが削減され、パズルはより小さく、解きやすくなります。

3. 「翻訳」（Bravyi–Kitaevの部分）

縮小され、折りたたまれたパズルができたら、次はそれを量子コンピュータが理解できる言語に翻訳する必要があります。そこには主に2つの翻訳者がいます。

• ジョルダン・ウィグナー（JW）: 標準的な翻訳者です。単純ですが、指示書が非常に長くなります（すべての単語が繰り返される本を読んでいるようなものです）。
• ブラヴィ・キタエフ（BK）: より賢い翻訳者です。情報をより効率的に整理するため、指示書がより短く、絡まりが少なくなります。

著者たちは、この「折りたたまれたパズル」の手法（SAE-CAS）が、どちらの翻訳者とも併用できることを示しました。彼らは、よりスマートな翻訳者を使用したSAE-CAS-BKというバージョンを作成しました。これは最終的な答えを変えるものではありませんが、そこに至るまでの経路をよりスムーズで高速にします。

何を見出したのか？

著者たちは、分子のエネルギーを見つけるために、2つの異なる「探索戦略（アルゴリズム）」を用いて、9つの小さな分子（水、酸素、窒素など）に対してこの手法をテストしました。

1. UCCSD: 化学的に精密ですが、複雑な探索。
2. HE-SCA: より単純で、ハードウェアに適した探索。

結果:

• より少ない量子ビット: 退屈な部分を無視し、対称な部分を折りたたむことで、量子ビットの数を大幅に削減できました（時には半分以下に抑えました）。
• より短い回路: シミュレーションを実行するための指示書が、より短く、絡まりも少なくなりました。
• 高い成功率: より単純な探索戦略（HE-SCA）を用いた際、彼らの手法はテストされたすべての分子において正しい答えを見つけ出しました。一方、従来の方法（JW-CAS）は、酸素や一酸化炭素において、制限時間内に答えを見つけられず失敗しました。
• 精度の損失なし: 退屈な電子を無視し、パズルを折りたたんだとしても、最終的なエネルギーの値は、標準的な大規模計算と同じ精度でした。

結論

著者たちは、「リソース効率の高い」ツールキットを構築しました。彼らは、変化しない部分を捨て、対称な部分を折りたたんでも、正しい答えを失うことなく、これらが安全に行えることを証明しました。これにより、以前考えられていたよりもはるかに小さく、低性能な量子コンピュータでも、これらの複雑な化学シミュレーションを実行することが可能になります。

彼らは、この「魔法のトリック」のコードを（QuantumSymmetryというパッケージとして）無料で公開しており、他の人々が自分自身の量子コンピュータで分子をシミュレートできるようにしています。

技術概要

技術要約：完全活性空間およびBravyi–Kitaevマッピングを用いた対称適応型量子ビットエンコーディング

問題提起
近未来（NISQ）およびフォールトトレラント量子プロセッサにおける量子化学シミュレーションは、分子の電子構造をエンコードするために必要な量子ビット数と、結果として得られるハミルトニアンの複雑さ（パウリ重みおよび回路深さ）に起因する重大なリソース制約に直面している。完全活性空間（CAS）近似は、コア軌道を凍結し仮想軌道を破棄することで問題のサイズを縮小するが、標準的な実装では、量子ビット数をさらに最小化するために対称性を十分に活用できていないことが多い。さらに、Jordan–Wigner（JW）やBravyi–Kitaev（BK）といった異なるフェルミオン・量子ビットマッピングは、局所性と量子ビット要件の間でトレードオフを提供しているが、これらが対称性削減技術と最適に統合されているとは限らない。

手法
著者らは、Symmetry-Adapted Encoding with Complete Active Space (SAE-CAS) と呼ばれる統一フレームワークを提案している。このアプローチは、CAS近似を単なるヒルベルト空間の切り捨てとしてではなく、凍結コアおよび仮想軌道に対応する近似的な $Z$ 対称性の賦課として扱うものである。

1. 理論的枠組み:

   • 近似対称性: 本手法は、正確な対称適応型エンコーディング（SAE）を、固定された軌道占有数（凍結コアおよび仮想軌道）に対応する近似的な対称性へと拡張する。これらは、$Z$ 対称性としてモデル化され、特定の量子ビットが $-1$（凍結）または$+1$（仮想）の固有値に拘束される。
   • アフィン・クリフォード変換: 削減は、アフィン・クリフォード写像を介して行われる。著者らは、対称性を表すバイナリ軌道特性行列 $S$ を定義する。$\mathbb{F}_2$ 上の可逆なアフィン写像 $T$ を構築することで、計算基底を変換し、対称性制約が特定の量子ビット上の固定値（通常は0）に対応するようにする。
   • 量子ビットの除去: クリフォード演算子によってハミルトニアンが共役変換された後、$X$ または $Y$ を含む項は破棄され、拘束された $Z$ 演算子はスカラーに置き換えられる。これにより、冗長な量子ビットが永久に除去され、システムは $n$ 個のスピン軌道から $n-k$ 個の量子ビットへと減少する（ここで $k$ は独立した対称性の数である）。
   • 等価性の証明: 著者らは、得られた活性空間量子ビットハミルトニアンが、分子電子構造論における凍結コアおよび仮想軌道の射影を介して導出される標準的なCASハミルトニアンと完全に等価であることを証明している。近似誤差は、エンコーディングではなく、活性空間の選択にのみ起因する。

2. Bravyi–Kitaevとの統合 (SAE-CAS-BK):

   • 本フレームワークは、Bravyi–Kitaev (BK) マッピングへと拡張される。BK変換自体が、計算ビット列に対するアフィン・クリフォード基底変換であるため、SAE-CASのクリフォード写像と合成することができる。
   • 得られる SAE-CAS-BK エンコーディングは、SAE-CASとユニタリ等価である。これは量子ビット数、パウリ項の数、および変分パラメータを保持するが、局所構造（項あたりのパウリ重み）および結果としての回路深さ／CNOTカウントを変更する。

3. 数値ベンチマーク:

   • 本手法は、最小の STO-3G 基底を用いた9つの小さな分子（H$_2$O、C$_2$H$_4$、O$_2$、N$_2$、C$_4$H$_6$ を含む）に対してテストされた。
   • 2つの変分アンサンブルが用いられた：一電子および二電子励起を含むユニタリ結合クラスター（UCCSD）、およびハードウェア効率的なシフト・サーキュラー・オルタネイティング（HE-SCA）アンサンブルである。
   • 標準的な JW、JW-CAS（JWを用い、かつ対称性削減を行わないCAS）、および回路レベルでの対称性フィルタリング（SF）を伴う JW-CAS と比較を行った。

主な結果

• リソース削減: SAE-CASは、JW-CASと比較して、量子ビット数、パウリ演算子数、回路深さ、CNOTカウント、および変分パラメータを一貫して削減した。例えば、C$_4$H$_4$ (4,4) CASにおいて、SAE-CASは量子ビット数を8個（JW-CAS）から4個へ、パウリ項の数を97個から46個へと削減した。
• 収束性と精度:
  • UCCSD: SAE-CASは、JW-CASおよび JW-CAS (SF) と同等の収束エネルギーに到達したが、より少ない最適化イテレーションを必要とした（例：N$_2$ では 109 回に対し 21 回）。特筆すべきは、C$_4$H$_4$ および C$_4$H$_6$ において、JW-CAS および JW-CAS (SF) の UCCSD アンサンブルは、励起リストの膨大さによるメモリ不足のため実行に失敗したが、SAE-CAS はマッピング段階で励起リストをフィルタリングすることで成功した。
  • HE-SCA: SAE-CASは優れた収束特性を示した。テストされたすべての分子において、SAE-CASはテストされたレイヤー予算内でCAS参照エネルギーに収束した。対照的に、JW-CASは O$_2$ および CO において200レイヤーまで収束に失敗した。著者らはこれを、JW-CASにおける「バレン・プラトー（平坦な台地）」効果とセクターリークに起因すると考えている。これは、JW-CASの変分多様体が誤ったスピン・パリティ・セクターを跨いでしまう現象である。一方、SAE-CASは、ネイティブに正しい対称性セクター内に探索を制限している。
  • SAE-CAS vs. SAE-CAS-BK: これら2つのエンコーディングは、同一の収束エネルギーとイテレーション回数を示し、それらのユニタリ等価性を裏付けた。差異は回路深さとCNOTカウントにおいてのみ観察され、ほとんどの分子で6%未満であった（COについては+20%の増加が見られた）。著者らは、これらの小さな活性空間サイズにおいては、BKの漸近的な $O(\log n)$ の局所性の利点がまだ実現されていないと指摘している。

意義と主張
本論文は、SAE-CASがリソース効率の高い分子シミュレーションのための系統的かつ構成可能な経路を提供すると主張している。CAS近似を近似的な対称性として扱うことで、本手法は厳密な対称性削減（点群およびスピン・パリティ）と活性空間の切り捨てを統合する。主な意義は以下の通りである：

1. ネイティブな対称性の強制: 対称性をエンコーディングレベルで強制することにより（アンサンブル内での事後的な処理ではなく）、冗長な探索方向を取り除き、最適化が誤った粒子数やスピンセクターへドリフトすることを防ぐ。
2. スケーラビリティ: 本手法は、アンサンブルの構築におけるメモリ制約（C$_4$H$_4$ や C$_4$H$_6$ のケースで見られた）によって、本来であれば実行不可能となるような大きな活性空間のシミュレーションを可能にする。
3. 柔軟性: 本フレームワークは、アフィン・クリフォード基底変換として表現可能なあらゆるフェルミオン・量子ビットマッピングと互換性があり、これは JW と BK の両方を用いた実験によって示されている。

著者らは、Pythonパッケージ QuantumSymmetry にオープンソースの実装を提供しており、SAE-CASが、近未来およびフォールトトレラント量子プロセッサの両方において、量子ビット数および演算の複雑さを削減するための実行可能な手段であることを結論付けている。