Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

技術要約：幾何学的に情報に基づいた変分量子固有値ソルバーを用いた、カゴメ格子上の反強磁性ハイゼンベルク模型の基底状態の探索

問題提起
本論文は、幾何学的フラストレーションと潜在的な量子スピン液体相で知られる系である、カゴメ格子上の反強磁性ハイゼンベルク模型の基底状態を見出すという課題に取り組んでいる。このような$k$-局所ハミルトニアンの基底状態を決定することはQMA困難な問題であり、ヒルベルト空間の指数関数的な増大により、大規模な系に対する古典的なシミュレーションは困難である。変分量子固有値ソルバー（VQE）は、ハイブリッドな量子・古典的アプローチを提供するものの、NISQ（Noisy Intermediate-Scale Quantum）デバイスに適している一方で、深い回路の必要性、ノイズ（コヒーレント誤差および非コヒーレント誤差）への感受性、勾配が消失するバレン・プラトー（barren plateaus）、そして表現力とハードウェア効率の両立が困難なアンサンブル設計といった大きな障壁に直面している。具体的には、著者らは、問題固有の複雑な簡約化や最適化ループ中の広範な誤差緩和に依存することなく、実機上でスピン・スピン相関や静的構造因子などの基底状態の特性を特徴付けることを目的としている。

手法
著者らは、カゴメ格子に特化した、幾何学的に情報に基づいたVQEアプローチを提案しており、単一の三角形とカゴメ・スター（6つの頂点を共有する三角形）という2つの基本的な幾何学構造に焦点を当てている。

1. 幾何学的に情報に基づいたアンサンブル（Ansatz）:
   核心となる革新は、カスタムメイドの浅いハードウェア効率的なアンザブルである。HVA（Hamiltonian Variational Ansatz）とは異なり、このアンザブルは、SWAPゲートを回避し、ネイティブなハードウェア接続性に制限されたパラメータ化された単一量子ビット回転（$R_y$）と絡み合いゲートのパターンを用いて構築されている。

   • ユークリッド・パラメータ空間: この設計の重要な特徴は、フビニ・ストディ・メトリック・テンソルが解析的に対角かつ一定であることが証明されている点である。これにより、パラメータ空間は自然にユークリッド的、特異点フリー、かつ過剰パラメータ化のないものとなる。その結果、量子自然勾配は標準的な勾配と等価になり、最適化のランドスケープが単純化される。
   • ハードウェア適合性: 回路の深さは最小限（1つの絡み合い層）に抑えられ、$R_y$ゲートの配置を交互に入れ替えることでパラメータ数を削減している。

2. 最適化戦略（AQNGD/AGD）:
   著者らは、適応型勾配降下法（AGD）オプティマイザを採用している。このアンザブルは対角かつ一定のフビニ・ストディ・メトリックを持つため、適応型量子自然勾配降下法（AQNGD）はAGDへと簡約される。学習率はアルモ（Armijo）のバックトラッキング規則を用いて動的に適応される。これにより、収束を加速させるための大きな初期ステップが可能になると同時に、オプティマイザが最小値に近づくにつれてステップサイズを自動的に縮小し、狭い幾何学的形状における停滞を回避できる。これは、同時摂動確率近似（SPSA）と比較されている。

3. 誤差緩和（ポストプロセッシング）:
   誤差緩和は、反復的なトレーニング中ではなく、VQE最適化の後に最終的なエネルギー推定値や構造的特性を精緻化するために適用される。

   • 読み出し誤差緩和（REM）: 分割応答行列アプローチを用い、SPAM（状態準備および測定）誤差を逆転させる。これは、疎な測定に対してフル行列反転よりも効率的にスケールする。
   • ゼロノイズ外挿（ZNE）: 著者らは、二次外挿（およびベイズ多項式回帰）を用いたZNEを適用し、結果をゼロノイズ極限へと外挿している。ZNEは厳密にはレイリー・リッツの変分原理に従わないが、ここでは精度向上のために使用されている。

主要な結果
研究は実機（IBM Oslo, IBM Kyoto, IBM Torino, および IBM Algiers）で実施され、状態ベクトルシミュレーションによって検証された。

• 収束と最適化:

  • 3量子ビットの三角形系において、AGD最適化VQEは、わずか19回のイテレーションで高い忠実度（シミュレーションで99.89%）で基底状態エネルギー（$E \approx -3$）へスムーズに収束した。最終的なパラメータは理論的な厳密値と密接に一致した。
  • 12量子ビットのカゴメ・スターにおいて、AGDはSPSAと比較して優れた収束性を示し、基底状態エネルギー（$E \approx -18$）に到達するために必要なイテレーション数が少なかった。シミュレーションにおける最終的な忠実度は96.37%であった。
  • アルモ規則によるステップサイズの適応は、大きなステップから開始し、収束が近づくにつれて精緻化するという、コストランドスケープのナビゲーションに効果的であることが示された。

• 構造的特性のノイズ耐性:

  • 重要な発見は、スピン・スピン相関と静的構造因子（$S(q)$）がノイズに対して強靭であることである。誤差緩和を行わない場合でも、実験的な$S(q)$マップは、基底状態の定性的なシグネチャ（例：強磁性秩序の不在を示すゼロ運動量での抑制）を保持していた。
  • ノイズを含む実験的な$S(q)$と厳密な理論的$S(q)$の間のピアソン相関係数は、極めて高い値（>99.9%）を維持した。
  • シミュレーション実験により、静的構造因子のノイズレベルに対する「ステップ状」の依存性が明らかになった。相関構造はある臨界ノイズ閾値までは頑健に維持されるが、その後急速に劣化し、相転移に似た挙動を示す。これは、エネルギーがノイズに対して単調に減衰する挙動とは対照的である。

• 誤差緩和の有効性:

  • REMおよびZNEを最適化後に適用することで、エネルギー推定値が大幅に改善された。カゴメ・スターの場合、生のエネルギーは約$-14.9$であったが、REMのみで約$-17.3$に改善され、二次外社のZNEを用いることで厳密値の$-18$にさらに近づいた。
  • 著者らは、REMとZNEを組み合わせると、基底エネルギーをオーバーシュート（超えてしまう）することがあることを観察した。これは、ZNEの非変分的な性質と一致する挙動である。ZNE単独の使用が、テストされた構成の中で最も高い精度を提供することが多かった。

意義と主張
本論文は、幾何学的に情報に基づいたアンザブル設計が、自然にユークリッド的なパラメータ空間を保証することで、複雑な問題固有の簡約化を必要とせずに、現在のNISQハードウェア上で効率的なVQEのトレーニングを可能にすることを主張している。

• 堅牢性: 本研究は、ノイズを含む量子ハードウェアから、誤差緩和中の最適化ループなしでも、基底状態を特徴付けるのに十分な精度で、スピン・スピン相関や静的構造因子といった意味のある物理的特性を抽出できることを実証した。
• 最適化効率: 適応型ステップサイズを持つAGDの使用は、この特定のハードウェア効率的なアンザブルに対してSPSAよりも効率的であり、イテレーション数においてより速い収束を実現した。
• 拡張性: 深い回路や重い誤差緩和をトレーニング中に避けることで、本アプローチはフラストレイテッド量子磁性体を調査するためのスケーラブルな経路を提供する。著者らは、このフレームワークをより大きなカゴメ格子や他のフラストレイテッド量子系に拡張できると考えている。
• 謙虚な姿勢: 著者らは、本手法が小さな基本セル（三角形とスター）のキャラクタライズには成功しているものの、無限カゴメ格子の一般的な基底状態は依然として文献における未解決問題であることを認めている。彼らは無限格子の問題を解決したと主張しているのではなく、現在のデバイスでこれらのシステムを調査するための実行可能な方法論を提示したものである。また、ZNEはエネルギー推定を改善するものの、変分原理を保持しないこと、およびREMとZNEの組み合わせには慎重な解釈が必要であることも述べている。