PASPT2: a size-extensive and size-consistent partial-active-space multi-state multi-reference second-order perturbation theory for strongly correlated electrons

本論文は、IN-GMS-SU-CCSD に基づいて導出され、特殊な参照特異的な零次ハミルトニアンを通じて非連結項を排除することで厳密なサイズ拡張性とサイズ整合性を達成する、新規な部分活性空間多状態多参照第二順序摂動理論である PASPT2 を導入する。

要約

以下は、論文「PASPT2」を平易な言葉と創造的な比喩を用いて解説したものです。

全体像：電子の「渋滞」を解決する

分子を賑やかな都市だと想像してください。電子は車であり、軌道は道路です。ほとんどの単純な分子では、交通はスムーズに流れています。ほぼすべての車が通る一本の主要ルートがあるからです。これは予測しやすいものです。

しかし、強相関系（複雑な遷移金属錯体など）では、交通は完全に地獄です。同等に優れた道路が多数あり、車は絶えず車線を変更し、合流し、分断し、巨大で混沌とした重ね合わせ状態になっています。単一のルートが支配的ではありません。このような系のエネルギーを予測することは、すべての車の動きが他のすべての車の動きに依存する巨大な渋滞の正確な流れを予測しようとするようなものです。

何十年もの間、化学者たちは、数学が破綻したり、不可能なほど遅くなったりすることなく、これらの「渋滞」のエネルギーを正確に計算することに苦慮してきました。

従来の方法の問題点

この論文は、PASPT2と呼ばれる新しい手法を紹介しています。それがなぜ特別なのかを理解するには、従来の方法の問題点を検討する必要があります。

1. 「完全」アプローチ（CAS）: 街のすべての可能な交通パターンをマッピングしようとするようなものです。これは正確ですが、大規模な都市では計算が不可能です。パターンの数が指数関数的に増えるためです。砂浜のすべての砂粒を数えようとするようなものです。
2. 「部分的」アプローチ（sCI）: 時間を節約するために、科学者たちは「最も重要な」交通パターンだけを選ぶようになりました。これは速いですが、欠点があります。街を相互作用しない二つの独立した町に分割した場合、数学が正しく足し合わされないことがあるのです。二つの別々のパーティーのコストを計算し、数学が混乱したために、二つの合計よりも高かったり低かったりする総額が出てくるようなものです。
3. 「普遍的」アプローチ（IN-GMS-SU-CC）: これを修正しようとした洗練された理論（IN-GMS-SU-CCSD）がありましたが、この論文の著者たちは隠れた欠陥を発見しました。それは実際には「サイズ広義性（size-extensive）」を持っていなかったのです。簡単に言えば、系が大きくなるにつれて誤差が一定に留まらず、増大するため、大規模分子に対する結果が信頼できなくなることを意味します。

解決策：PASPT2

著者である劉春章、寧張、劉文建は、PASPT2を開発しました。これがどのように機能するかを、比喩を用いて説明します。

「部分的活性空間（PAS）」戦略
街全体をマッピングする代わりに、PASPT2 は交通が最も混沌としている特定の「ダウンタウン地区」（活性空間）に焦点を当てます。すべての可能性ではなく、最も重要な交通パターン（配置）の賢明な部分集合を選択して集中します。

「特別なゼロ次ハミルトニアン」（交通警官）
核心的な革新は、道路のルール（数学）を設定する新しい方法にあります。

• 従来の方法: 従来の数学には「非連結項」がありました。交通報告書で、「車 A が動いており、別々に車 B も動いている」と述べ、しかし車 A の動きが実際には車 B の経路を変えるという事実を報告書が考慮していなかったようなものです。これにより、都市が大きくなるにつれて誤差が積み重なりました。
• PASPT2 の方法: 著者たちは特別な「交通警官」（参照固有のゼロ次ハミルトニアン）を導入しました。この警官は、すべての計算が「連結」であることを保証します。それは、システムのすべての部分が相互に関連していることを数学に認識させるように強制します。これにより、以前の理論を悩ませてきた「非連結項」を排除しました。

結果：完璧にバランスの取れた秤
この新しい「交通警官」のおかげで、PASPT2 はサイズ広義性とサイズ整合性を持っています。

• サイズ広義性: 分子のサイズを倍にすると、エネルギー計算も完璧に倍になります。誤差は増大しません。
• サイズ整合性: 遠く離れた二つの独立した分子がある場合、それらを一緒に計算した総エネルギーは、個別に計算したエネルギーの合計と正確に一致します。距離によって数学が混乱することはありません。

検証方法

著者たちは単に数学を書き留めただけでなく、現実世界の「渋滞」でそれをテストしました。

1. ヘリウム鎖: ヘリウム原子を家の一列のように並べました。より多くの家を追加するにつれて、エネルギー計算が完全に直線的に増大することを示し、この手法が大規模系に対して信頼性があることを証明しました。
2. 水分子（H2O）: 電子をより高いエネルギー準位にジャンプさせるために必要なエネルギー（励起）を計算しました。彼らの結果を「ほぼ正確な」基準値と比較しました。彼らの手法は非常に優れていましたが、精度は初期の「交通パターン」（モデル空間）をどの程度よく選んだかに大きく依存していることが分かりました。より良い出発点を選べば、結果はほぼ完璧になりました。
3. 窒素分子（N2）: これは化学結合の切断をテストする古典的なケースです。窒素原子が離れるにつれて、「交通」は非常に混沌とします。PASPT2 は、原子が遠く離れていても、最も正確な利用可能な基準値と一致するよう、エネルギー曲線をスムーズに追跡することに成功しました。

結論

この論文は、PASPT2が以下の点においてその種類としては初の手法であるため、画期的であると主張しています。

1. 部分的活性空間に基づいている（複雑な分子に対して十分な速度を持つ）。
2. 多状態である（複数のエネルギー準位を同時に処理できる）。
3. 厳密にサイズ広義性とサイズ整合性を持つ（大規模系に対して数学的に信頼できる）。
4. イントラダーフリーである（エネルギー準位が近くなりすぎたときに発生する数学的な「不具合」を回避する）。

著者たちは、この手法は現在、大きな前進であるが、次の課題は「スピン適応」化（電子の磁気スピンをさらに完璧に処理すること）であり、それは今後の作業で計画されていると結論付けています。現時点では、PASPT2 は化学における最も困難な電子系を理解するための、堅牢で正確かつ拡張可能なツールを提供しています。

技術概要

以下は、論文「PASPT2: A Size-Extensive and Size-Consistent Partial-Active-Space Multistate Multireference Second-Order Perturbation Theory for Strongly Correlated Electrons」の詳細な技術的概要です。

1. 問題提起

強相関系（多核遷移金属錯体など）の正確な量子化学記述は、依然として重大な課題です。これらの系では、単一の配置が支配的ではない電子配置の巨大な重ね合わせの扱いが必要です。

• 既存手法の限界:
  • CASCI/CASSCF: 基礎的ではあるものの、活性軌道の数に対して組合せ的にスケーリングするため、約 12 軌道までに制限されます。また、静的相関は厳密に扱いますが動的相関は近似で扱うため、不均衡な扱いや、高エネルギー状態における「イントレーダー状態」（摂動論における特異点）を引き起こします。
  • Selected CI + PT2 (sCIPT2): 効果的ですが、sCIPT2 で用いられる標準的な Epstein-Nesbet 摂動論は、完全ヒルベルト空間ではない選択された空間に適用された場合、サイズ広大性（size-extensivity）を持ちません。
  • IN-GMS-SU-CCSD: Li と Paldus が提案した親理論（Intermediate Normalization-based General-Model-Space State-Universal Coupled Cluster with Singles and Doubles）は汎用的を目指していますが、その振幅方程式に**非連結項（disconnected terms）**が含まれるため、以前の主張とは異なり、サイズ広大性を欠くという問題があります。
• ギャップ: 部分活性空間（PAS）に基づき、厳密にサイズ広大かつサイズ一貫性（size-consistent）を持ち、かつイントレーダー状態を回避する多状態（MS）・多参照（MR）・**第二階摂動論（PT2）**手法の必要性があります。

2. 手法：PASPT2

著者らは、IN-GMS-SU-CCSD 理論を線形化することでPASPT2を定式化しました。核心的な革新は、ゼロ次ハミルトニアンの特定の構築とモデル空間の扱いにあります。

A. 理論的枠組み

• 波動演算子: Jeziorski-Monkhorst (JM) アンサッツに基づき、部分活性空間（PAS, $M_0$）を正確な多様体へ写像します。
• 参照特異的な分割: 単一のグローバルハミルトニアンを使用する標準的なアプローチとは異なり、PASPT2 は各参照行列式 $|\alpha\rangle$ に対して**参照特異的なゼロ次ハミルトニアン（$H_0^\alpha$）**を採用します。
  • $H_0^\alpha = H_{0,I}^\alpha + H_{0,A}^\alpha$
  • 不活性部分（$H_{0,I}^\alpha$）: 不活性軌道に作用する対角型（フォック型）演算子。
  • 活性部分（$H_{0,A}^\alpha$）: 活性軌道に作用する非エルミート演算子。重要なのは、振幅方程式における非連結項を生成する特定の非対角結合を排除するように構築されている点です。
• 振幅方程式: 一次クラスタ振幅（$t$）は、線形系 $At = V$を解くことで得られます。$H_0^\alpha$ の特定の選択により、これらの方程式の右辺は**項ごとに連結（termwise connected）**となり、IN-GMS-SU-CCSD に見られた非連結項が排除されます。

B. イントレーダー状態の扱いと閉包

• 中間ハミルトニアン: 有効ハミルトニアン（$\bar{H}_{eff}^{[2]}$）は連結かつ閉じたものとして構築されます。
• バッファ空間（$\bar{M}_0$）: ハミルトニアンを閉じたもの（参照状態を外部空間 $Q$ へ励起しない）にするために、モデル空間は自動的に $M_0$ から中間空間 $M_X = M_0 \cup \bar{M}_0$ に拡張されます。
  • $\bar{M}_0$ 内の状態は「バッファ」として機能します。これらはハミルトニアンを閉じるために含まれますが、摂動されません（それらの $T$-クラスタはゼロに設定されます）。
  • このメカニズムは、アドホックなレベルシフトなしにイントレーダー状態の問題を実質的に解決します。

C. サイズ広大性とサイズ一貫性

• サイズ広大性: $t$-振幅と有効ハミルトニアンが厳密に連結であることを示すことで証明されます。$H_{0,A}^\alpha$ の特定の形式は、親理論である IN-GMS-SU-CCSD には欠けていた特徴である、共通軌道項が正しく相殺されることを保証します。
• サイズ一貫性: 非相互作用フラグメント $A$ と $B$ の PAS の直積である PAS を持つ超分子 $AB$に対して証明されます。$H_0^\alpha$の加法的分離可能性により、$AB$のエネルギーは$E_A + E_B$ に等しくなります。

3. 主要な貢献

1. 親理論の修正: 著者らは、以前提案された IN-GMS-SU-CCSD が非連結項によりサイズ広大性を欠くことを厳密に実証し、この欠陥が PASPT2 で修正されたことを示しました。
2. 新規ゼロ次ハミルトニアン: 振幅方程式における項ごとの連結性を強制する、参照特異的な非エルミート活性部分（$H_{0,A}^\alpha$）の導入。
3. 自動イントレーダー状態解決: 自動的に決定されるバッファ空間（$\bar{M}_0$）を介した閉じた中間ハミルトニアン（$\bar{H}_{eff,X}^{[2]}$）の開発により、手動のレベルシフトの必要性を排除しました。
4. 初の PAS ベースの MS-MRPT2: PASPT2 は、中間規格化の下で部分活性空間に基づいた、初のサイズ広大・サイズ一貫・イントレーダーフリーな多状態 MRPT2 手法として提示されました。

4. 結果と検証

この手法はBDFプログラムパッケージに実装され、3 つの代表的な系でテストされました。

• 直鎖状ヘリウム鎖（$He_n$）:
  • テスト: サイズ広大性。
  • 結果: 相関エネルギーは原子数に対して線形にスケーリングし、MP2 および準正確な iCIPT2 の結果と一致しました。
• 水（$H_2O$）の垂直励起:
  • テスト: iCIPT2 ベンチマークに対する励起エネルギー（VEE）の精度。
  • 知見:
    • 基底状態の HF 軌道を使用すると中程度の誤差が生じました。
    • 半標準化 SA-CASSCF 軌道を使用すると低エネルギー状態は改善されましたが、高エネルギー状態は悪化しました。
    • 最適戦略: 手動で選択されたモデル空間（$M_0$）を$\bar{H}_{eff,X}^{[2]}$バリアントを用いて閉じた空間（$M_X$）に拡張することが、最良のバランスを提供しました。
    • この手法は、モデル空間が良質であれば、43 の励起状態を成功裡に記述し、平均絶対誤差（MAE）は IN-GMS-SU-CCSD や CASSCF-SDSPT2 と同等かそれ以上でした。
• 窒素（$N_2$）のポテンシャルエネルギー曲線（PEC）:
  • テスト: $^1\Sigma_g^+$ および $^3\Sigma_u^+$ 状態の解離挙動と分光定数（$R_e, \omega_e, D_e$）。
  • 課題: $N_2$ の解離には CAS(6,6) が必要であり、選択を行わない標準的な MRPT2 ではしばしば大きすぎます。
  • 知見:
    • PASPT2 の結果は iCIPT2 ベンチマークとよく一致しました。
    • 改良: 結合領域において、いくつかの非結合振幅が物理的に不合理に大きくなり（ブリルアン条件に違反）、PEC の滑らかさや非平行性誤差（NPE）に問題がありました。振幅 $|t| > 0.3$ の外部励起を除去することで、PEC の滑らかさと NPE が大幅に改善されました。
    • 両状態の分光定数は非常に高精度であり、標準的な CASSCF を凌駕し、CASSCF-SDSPT2 と同等の性能を示しました。

5. 意義

• 理論的厳密性: PASPT2 は、一般モデル空間多参照摂動論における長年の課題であるサイズ広大性の問題を解決しました。
• 実用的適用性: 完全活性空間（CAS）ではなく**部分活性空間（PAS）**を利用することで、適切な初期配置選択がなされていれば、従来の CASSCF ベースの手法よりも大きな活性空間を持つ系を扱えるようになります。
• 堅牢性: バッファ空間（$\bar{M}_0$）の自動構築により、MRPT2 で一般的な失敗要因であるイントレーダー状態に対して堅牢になります。
• 将来展望: 著者らは、現在の実装がスピン軌道（対角型不活性部分）を使用しているものの、将来的にはスピン適応に焦点を当て、開殻系の精度をさらに向上させることを目指すと述べています。

要約すると、PASPT2 は電子構造理論における重要な進歩であり、強相関系に対して、多参照手法の精度と摂動論の効率性を両立させた、厳密かつサイズ広大で計算的に実行可能なアプローチを提供します。