【DFT】ダブルハイブリッド汎関数の概要：数理構造と特性について

はじめに#

B3LYPなどのハイブリッド汎関数は、DFTの交換エネルギーの一部をハートリー・フォック（HF）交換で置き換える手法ですが、ダブルハイブリッド（Double-Hybrid, DH）汎関数は、さらに相関エネルギーの一部も波動関数理論（MP2）で置き換えるというアプローチを採ります。

2006年にStefan Grimmeによって提案された B2PLYP 以降、この手法は計算化学における「ヤコブの梯子」の第5段に位置づけられる高精度な手法として広く利用されています。本稿では、その数理的構造と特性について解説します。

1. 定義と数理構造#

通常のハイブリッド汎関数が「交換」のみを混合するのに対し、ダブルハイブリッド汎関数は 「交換」と「相関」の両方 にパラメータを持ちます。その一般式は以下の通りです。

E_{XC}^{\text{DH}} = \underbrace{(1 - a_x) E_X^{\text{DFT}} + a_x E_X^{\text{HF}}}_{\text{交換項}} + \underbrace{(1 - a_c) E_C^{\text{DFT}} + a_c E_C^{\text{PT2}}}_{\text{相関項}}

$a_x$ : HF交換の混合比。
$a_c$ : MP2（2次摂動論）相関の混合比。
$E_C^{\text{PT2}}$ : MP2法による相関エネルギー補正。

計算の手順#

DH-DFTの計算は一般に以下の2段階で行われます。

SCF計算: 通常のDFT（MP2項なし）を行い、Kohn-Sham軌道と軌道エネルギーを収束させます。
摂動計算: 得られた軌道を用いて、以下のMP2公式に基づき相関エネルギーを算出します。

E_C^{\text{PT2}} = \frac{1}{4} \sum_{ij}^{\text{occ}} \sum_{ab}^{\text{virt}} \frac{|\langle ij || ab \rangle|^2}{\epsilon_i + \epsilon_j - \epsilon_a - \epsilon_b}

このMP2項は空軌道（Virtual Orbitals）の情報を含むため、電子の動的相関（分散力など）を記述する上で寄与します。

2. パラメータ決定のアプローチ#

混合パラメータ $a_x$ と $a_c$ の決定方法には、大きく分けて2つのアプローチが存在します。

アプローチA：経験的フィッティング#

代表例：B2PLYP, mPW2-PLYP

手法: G3データセットなどの実験値（原子化エネルギー等）との誤差が最小になるように係数を決定します。
- 例：B2PLYP ( $a_x=0.53, a_c=0.27$ )
特徴: ターゲットとする化学プロパティに対して高い精度を示しますが、トレーニングセットに含まれない系での挙動については留意が必要です。

アプローチB：理論的導出#

代表例：PBE0-DH

手法: 「断熱接続（Adiabatic Connection）」の理論形式に基づき、係数間の関係式を導きます。
- 例： $a_c \approx a_x^2$ や $a_c \approx a_x^3$ といった制約を用いる。
特徴: パラメータ数が少なく理論的な背景を持ちますが、実用的な精度においては経験的な手法と同程度、あるいは系によっては劣る場合もあります。

3. Post-MP2への拡張とその限界#

「MP2（2次摂動）よりも高次の摂動項（MP3, MP4など）を用いれば精度が向上するか」という点について、Chan, Goerigk, Radomらによる検証が行われています。

検証結果#

報告によると、高次項の追加は必ずしも精度の向上にはつながらず、場合によっては精度が低下することも確認されています。これは、DH-DFTの精度が 「DFTの誤差とMP2の誤差の相殺（Error Cancellation）」 に依存している部分が大きいためと考えられています。摂動の次数を上げることでこのバランスが崩れる可能性があります。

スピン成分スケーリング (SCS)#

一方で、MP2項をスピン成分ごとにスケーリングする手法（DSD-DFT）は有効とされています。

a_c E_C^{\text{PT2}} \rightarrow c_O E_C^{\text{OS-PT2}} + c_S E_C^{\text{SS-PT2}}

平行スピンと反対スピンの寄与を個別に調整することで、分散力や反応障壁の記述が改善される傾向にあります。

4. 計算コストと留意点#

ダブルハイブリッド汎関数の利用にあたっては、以下のコストと制約を考慮する必要があります。

① 計算コスト#

MP2項の計算には $O(N^5)$ （原子数の5乗）のスケールでコストがかかります。通常のDFTは $O(N^3)$ 〜 $O(N^4)$ であるため、大規模な系への適用には計算リソースが必要です。

② 基底関数依存性#

波動関数理論に基づくMP2項が含まれるため、基底関数のサイズに対する依存性が強くなります。

DFT: 比較的小さな基底関数でも収束しやすい。
MP2: 収束が遅く、大きな基底関数が必要。

そのため、DH-DFTで十分な精度を得るには、一般的に Triple-Zeta (TZVPP) 以上 の基底関数の使用が推奨されます。不十分な基底関数を用いた場合、期待される精度が得られない可能性があります。

まとめ#

ダブルハイブリッド汎関数（DH-DFT）は、DFTと波動関数理論を組み合わせることで、相関エネルギーの記述を改善した手法です。

B2PLYP 等の経験的手法は実用的な選択肢として定着しています。
理論的なアプローチや、スピン成分を考慮した DSD-DFT などの改良も進められています。
利用の際は、計算コストの増加と、適切な基底関数の選択が必要である点を理解しておくことが重要です。

参考文献#

S. Grimme, J. Chem. Phys. 124, 034108 (2006).
D. Bousquet et al., J. Chem. Theory Comput. 9, 3444 (2013).
B. Chan et al., J. Comput. Chem. 37, 183 (2016).

【DFT】ダブルハイブリッド汎関数の概要：数理構造と特性について

https://ss0832.github.io/posts/20251228_dft_double_hybrid_functional/

Author

ss0832

Published at

2025-12-28

License

CC BY-NC-SA 4.0

【DFT】B97系汎関数の到達点：ωB97M-V の理論的背景と特徴

【DFT】物理学の標準原器：PBE汎関数の哲学と構造