Multiconfiguration Pair-Density Functional Theory Calculations of Low-lying States of Complex Chemical Systems with Quantum Computers

以下は、提示された論文「Multiconfiguration Pair-Density Functional Theory Calculations of Low-lying States of Complex Chemical Systems with Quantum Computers（量子コンピュータを用いた複雑化学系の低励起状態に対する多配置ペア密度汎関数理論計算）」の技術的な要約です。

1. 研究の背景と課題 (Problem)

量子化学計算において、複雑な分子系の電子構造を正確に記述する際の最大の課題は、「静的相関（static correlation）」と「動的相関（dynamic correlation）」の両方を同時に適切に扱うことです。

• 従来の古典計算の限界: 単一参照法（HF, DFT, CC など）は動的相関を扱えますが、多重参照性が必要な系（結合解離や遷移金属など）では失敗します。一方、多参照法（CASSCF など）は静的相関を扱えますが、動的相関を欠くため精度が低く、完全な記述には計算コストが指数関数的に増大します。
• 量子計算の現状: 変分量子固有値ソルバー（VQE）は量子コンピュータの強みを活かす枠組みですが、現在の NISQ（ノイズあり中規模量子）デバイスでは、コヒーレンス時間の制約や回路の深さ（ゲート数）の限界により、動的相関を完全に記述するための大規模な回路（例：UCCSD）を実行することが困難です。特に、強相関系（例：Cr2 二原子分子）を扱うには、非常に大きな活性空間が必要となり、既存の VQE 手法では実用的ではありません。

2. 提案手法：VQE-MC-PDFT (Methodology)

本研究では、VQE と多配置ペア密度汎関数理論（MC-PDFT）を統合したハイブリッド量子・古典アルゴリズム「VQE-MC-PDFT」を提案しました。この手法は、相関効果を効率的に分離・処理する戦略を採用しています。

• 静的相関の量子処理: 化学的に重要な「活性空間（Active Space）」内の静的相関を、量子コンピュータ上でパラメータ化された多配置量子状態（CASCI 型のアンサッツ）を用いて VQE で解きます。これにより、電子の多重参照性を量子リソースで直接扱います。
• 動的相関の古典処理: 残りの動的相関（短距離相関）は、量子計算で得られた縮約密度行列（RDM: 1-RDM と 2-RDM）を古典コンピュータに渡し、古典的な「オン・トップ・ペア密度汎関数（on-top pair-density functional）」を用いて回復させます。本研究では ftPBE 汎関数を使用しました。
• 自己無撞着な軌道最適化: 従来の VQE-PDFT 手法の延長として、MC-PDFT エネルギーの軌道回転勾配（orbital-rotation gradient）を計算し、これを量子回路の「軌道回転ユニタリ演算子」としてエンコードします。これにより、配置係数と分子軌道の両方を自己無撞着に変分的に最適化するループを実現し、物理的な厳密性を保ちつつ量子リソースのオーバーヘッドを大幅に削減しました。
• 拡張技術:
  • 回路分割戦略: 弱くエンタングルしたクラスターを分解し、古典的な再構成を通じてより大きな活性空間を扱うことを可能にしました。
  • 密度情報に基づく誤差軽減: 測定ノイズに対する誤差軽減プロトコルとして、電子密度情報を利用した FEM（有限要素法）に着想を得た反復的なノイズモデル補正を採用し、RDM とエネルギーの精度を向上させました。

3. 主要な貢献と成果 (Key Contributions & Results)

本研究は、以下のベンチマーク分子において化学的精度（chemical accuracy）を達成し、NISQ デバイスでの実用性を示しました。

• 二炭素分子（C2）:

  • 強相関を持つ C2 の基底状態および励起状態のポテンシャルエネルギー曲線（PEC）を計算。
  • 平衡結合長において、実験値との平均絶対誤差（MAE）が 0.006 Å、励起エネルギーの MAE が 0.048 eV となり、高品質な古典計算（SHCI）や実験値と極めて良く一致しました。
  • VQE-SA-CASSCF と比較し、動的相関の回復により結合エネルギーの過小評価を解消し、より滑らかで物理的に正しい曲線を得ました。

• 二クロム分子（Cr2）:

  • 化学結合の「六重結合」として知られ、非常に強い静的・動的相関を持つ Cr2 を対象にしました。
  • 非常に大きな活性空間 CAS(48e, 42o)（実質的に 84 量子ビット相当）を、回路分割とクラスター再構成技術を用いて処理しました。
  • 従来の VQE-CASSCF は結合曲線が反発的（結合しない）になるなど失敗しましたが、VQE-MC-PDFT は束縛されたポテンシャルエネルギー曲線を再現し、解離挙動を定性的に正しく記述しました。
  • 84 量子ビット相当のスケールで得られた結合エネルギー（-0.465 eV）は、高レベルな古典ベンチマーク（MPS-LCC や SHCI）と良好に一致し、量子リソースの増加に伴う系統的な収束を確認しました。

• ベンゼン（Benzene）:

  • 芳香族系の励起状態を計算し、QUEST データベースの理論値との MAE を 0.048 eV としました。これは最先端の NES-VMC 法と同等の精度であり、標準的な TDDFT（MAE 0.37 eV）を大幅に上回ります。

4. 意義と結論 (Significance)

• NISQ 時代の新たなパラダイム: 本研究は、量子コンピュータに「静的相関の処理」という最も本質的な役割を割り当て、残りの「動的相関」を効率的な古典汎関数で補うという分離戦略の有効性を実証しました。これにより、回路の深さを抑えつつ、強相関系に対する高精度な予測が可能になりました。
• スケーラビリティの実証: 回路分割技術と誤差軽減手法を組み合わせることで、単一の量子ハードウェアの物理量子ビット数を超えた大規模な活性空間（Cr2 の 48 電子 42 軌道）をシミュレーション可能であることを示しました。
• 将来展望: 半芯電子（semi-core）の相関や、より複雑な励起状態への汎関数の転送性など、今後の課題は残されていますが、VQE-MC-PDFT は、近未来の量子ハードウェア上で信頼性の高い化学予測を行うための実用的なテンプレートとして確立されました。

要約すれば、この論文は「量子計算の限界を克服し、複雑な化学系を高精度にシミュレートするための、量子と古典の最適な役割分担（ハイブリッドアプローチ）」を提案し、その有効性を C2、Cr2、ベンゼンなどの難問系で実証した画期的な研究です。