多体微扰理论是以Hartree-Fock方程为基础的,应用微扰法处理的计算方法。微扰法要求将复杂体系的哈密顿算子分解为可精确求解项和微扰项两部分,在多体微扰理论中,引入Hartree-Fock哈密顿算子的概念:
可以证明,由Hartree-Fock方程解得的单电子分子轨道波函数所构成的斯莱特行列式波函数是 的本征函数,构成斯莱特行列式的各分子轨道轨道能的代数和是 的本征值:
将多电子体系哈密顿算子分解为Hartree-Fock哈密顿算子和微扰项的代数和:
|
|
|
|
在这个假设下,多电子体系电子哈密顿算子被分解为可精确求解的Hartree-Fock哈密顿算子和微扰算子,应用微扰方法进行近似处理。
零级能量和波函数
在多体微扰理论下,基态零级能量就是构成基态斯莱特行列式的各分子轨道轨道能的代数和,零级波函数就是基态斯莱特行列式波函数。可以看出,多体微扰理论的零级能量精度甚至不如Hartree-Fock方程所得的能量。
一级校正
根据微扰理论,能量的一级校正 为:
将微扰算子V的表达式代入得到:
|
|
|
|
可以看到,经过能量的一级校正后体系能量为:
|
|
|
|
|
|
二级校正
根据微扰理论,体系基态能量的二级校正 为:
其中 为Hartree-Fock哈密顿算子本征能量为 的波函数,其本质是体系激发态的斯莱特行列式
可以证明,只有对双激发的斯莱特行列式才有 所以体系能量的二级校正为:
将分子项展开,得到:
最终体系经过二级校正的基态能量为:
由于式中 是体系未占据分子轨道的轨道能,在基态,其能量衡高于 和 所以能量的二级微扰是一个负值,因而考虑二级微扰的体系能量低于Hartree-Fock方程得到的体系能量,这一差异来自电子相互作用。
考虑二级校正的多体微扰计算简称MP2
高级校正
更高级的校正是以较低级校正为计算基础的,随着校正级别的提高,计算量也急剧增加,理论上讲,随着校正级别的提高最终的体系能量会逐渐逼近真实值。目前的计算方法最高可以进行MP5计算,即体系能量的五级校正。