多体微扰理论

✍ dations ◷ 2025-09-12 19:50:59 #量子化学,微扰理论

多体微扰理论是一种基于分子轨域理论的高级量子化学计算方法。这种方法以Hartree-Fock方程的自洽场解为基础,应用微扰理论,获得考虑了相关能的多电子体系近似解,其计算精度与组态相互作用方法的DCI接近,但计算量远小于DCI,是应用比较广泛的高级量子化学计算方法。

多体微扰理论是由量子化学家Møller和Plesset在1934年提出的,所以这一方法也经常以二人的名字所写MP表示,MPn表示的是多体微扰n级近似。

多体微扰理论是以Hartree-Fock方程为基础的,应用微扰法处理的计算方法。微扰法要求将复杂体系的哈密顿算子分解为可精确求解项和微扰项两部分,在多体微扰理论中,引入Hartree-Fock哈密顿算子的概念:

H 0 = i N f i {\displaystyle H_{0}=\sum _{i}^{N}f_{i}}

可以证明,由Hartree-Fock方程解得的单电子分子轨道波函数所构成的斯莱特行列式波函数是 H 0 {\displaystyle H_{0}} 的本征函数,构成斯莱特行列式的各分子轨道轨道能的代数和是 H 0 {\displaystyle H_{0}} 的本征值:

H 0 | Ψ 0 >= i N ϵ i | Ψ 0 > {\displaystyle H_{0}|\Psi _{0}>=\sum _{i}^{N}\epsilon _{i}|\Psi _{0}>}

将多电子体系哈密顿算子分解为Hartree-Fock哈密顿算子和微扰项的代数和:

H e l e = H 0 + V {\displaystyle H_{ele}=H_{0}+V}

在这个假设下,多电子体系电子哈密顿算子被分解为可精确求解的Hartree-Fock哈密顿算子和微扰算子,应用微扰方法进行近似处理。

在多体微扰理论下,基态零级能量就是构成基态斯莱特行列式的各分子轨道轨道能的代数和,零级波函数就是基态斯莱特行列式波函数。可以看出,多体微扰理论的零级能量精度甚至不如Hartree-Fock方程所得的能量。

根据微扰理论,能量的一级校正 E 0 ( 1 ) {\displaystyle E_{0}^{(1)}} 为:

E 0 ( 1 ) =< Ψ 0 ( 0 ) | V | Ψ 0 ( 0 ) > {\displaystyle E_{0}^{(1)}=<\Psi _{0}^{(0)}|V|\Psi _{0}^{(0)}>}

将微扰算子V的表达式代入得到:

可以看到,经过能量的一级校正后体系能量为:

根据微扰理论,体系基态能量的二级校正 E 0 ( 2 ) {\displaystyle E_{0}^{(2)}} 为:

E 0 ( 2 ) = n 0 | < Ψ 0 ( 0 ) | V | Ψ n ( 0 ) > | 2 E 0 ( 0 ) E n ( 0 ) {\displaystyle E_{0}^{(2)}=\sum _{n\neq \;0}{\frac {{\begin{vmatrix}<\Psi _{0}^{(0)}|V|\Psi _{n}^{(0)}>\end{vmatrix}}^{2}}{E_{0}^{(0)}-E_{n}^{(0)}}}}

其中 Ψ n ( 0 ) {\displaystyle \Psi _{n}^{(0)}} 为Hartree-Fock哈密顿算子本征能量为 E n ( 0 ) {\displaystyle E_{n}^{(0)}} 的波函数,其本质是体系激发态的斯莱特行列式

可以证明,只有对双激发的斯莱特行列式才有 | < Ψ 0 | V | Ψ a , b r , s > | 0 {\displaystyle {\begin{vmatrix}<\Psi _{0}|V|\Psi _{a,b}^{r,s}>\end{vmatrix}}\neq \;0} 所以体系能量的二级校正为:

E 0 ( 2 ) = a < b , r < s | < Ψ 0 | V | Ψ a , b r , s > | 2 ϵ a + ϵ b ϵ r ϵ s {\displaystyle E_{0}^{(2)}=\sum _{a<b,r<s}{\frac {{\begin{vmatrix}<\Psi _{0}|V|\Psi _{a,b}^{r,s}>\end{vmatrix}}^{2}}{\epsilon _{a}+\epsilon _{b}-\epsilon _{r}-\epsilon _{s}}}}

将分子项展开,得到:

E 0 ( 2 ) = a < b , r < s | < χ a χ b | | χ r χ s > | 2 ϵ a + ϵ b ϵ r ϵ s {\displaystyle E_{0}^{(2)}=\sum _{a<b,r<s}{\frac {{\begin{vmatrix}<\chi _{a}\chi _{b}||\chi _{r}\chi _{s}>\end{vmatrix}}^{2}}{\epsilon _{a}+\epsilon _{b}-\epsilon _{r}-\epsilon _{s}}}}

最终体系经过二级校正的基态能量为:

E = E 0 H F + a < b , r < s | < χ a χ b | | χ r χ s > | 2 ϵ a + ϵ b ϵ r ϵ s {\displaystyle E=E_{0}^{HF}+\sum _{a<b,r<s}{\frac {{\begin{vmatrix}<\chi _{a}\chi _{b}||\chi _{r}\chi _{s}>\end{vmatrix}}^{2}}{\epsilon _{a}+\epsilon _{b}-\epsilon _{r}-\epsilon _{s}}}}

由于式中 ϵ r {\displaystyle \epsilon _{r}} ϵ s {\displaystyle \epsilon _{s}} 是体系未占据分子轨道的轨道能,在基态,其能量衡高于 ϵ a {\displaystyle \epsilon _{a}} ϵ b {\displaystyle \epsilon _{b}} 所以能量的二级微扰是一个负值,因而考虑二级微扰的体系能量低于Hartree-Fock方程得到的体系能量,这一差异来自电子相互作用。

考虑二级校正的多体微扰计算简称MP2

更高级的校正是以较低级校正为计算基础的,随着校正级别的提高,计算量也急剧增加,理论上讲,随着校正级别的提高最终的体系能量会逐渐逼近真实值。目前的计算方法最高可以进行MP5计算,即体系能量的五级校正。

多体微扰理论方法是一种量子化学高级计算方法,在考虑相关能的计算方法中,多体微扰理论方法是计算量相对最小的。MP1可以达到HF方程的计算精度,MP2一般可以达到60%的相关能,与DCI方法相当,但计算过程仅需要计算少量双电子积分,远远小于DCI;MP4一般可以达到85%的相关能。

MPn方法是一个大小一致的方法,即对电子数不同的体系,使用MPn计算的精度是相同的,这一特性使得MPn方法特别适合进行化学反应的模拟计算。但是由于MPn方法以HF方程为基础,因而受到HF方程的局限,对于那些应用HF方程不能很好处理的体系,如非限制性开壳层体系,MPn方法也不能很好处理。

相关

  • 生态区从生态学和地理学上讲,生态区小于生物区,而生物区又小于生物地理分布区。生态区涵盖了大范围的陆地或水域,有着独特的天然群落和物种的组合。植物、动物和生态系统的生物多样性
  • 江苏省侯集高级中学江苏省侯集高级中学始建于1956年,位于江苏省徐州市大庙镇,是江苏省四星级重点中学。校园占地300余亩,建筑面积15多万平方米。
  • BZFlag《BZFlag》 ( ) 是三维第一人称坦克射击游戏。Chris Schoeneman原为SGI系列计算机编写,基于经典游戏 Battlezone。Tim Riker和一个团队维护现在的版本,而且已经开源了。这个项
  • 宾夕法尼亚州同性婚姻宾夕法尼亚州同性婚姻于2014年5月20日合法化。20日美国联邦地区法院法官裁定1996年州禁止同性婚姻的法律是违宪的。宾夕法尼亚州此前不承认民事结合,是美国东北部地区最后一
  • 刘良佐刘良佐,字明辅,中国山西大同左卫人。明朝末年至南明时期的军事人物。刘良佐与高杰本来是李自成麾下的战将,崇祯十一年十月,中曹变蛟伏兵而投降明军,被封广昌伯。崇祯十四年,击败饥
  • 武玉笑武玉笑(1929年2月6日-2018年8月9日),男,陕西葭县人,中国剧作家,甘肃省文学艺术界联合会原副主席,甘肃省作家协会原主席。
  • 和穆熙和穆熙(1930年-),笔名巴紫、马伏,男,河南沁阳人,中国编辑、作家,曾任长江文艺出版社总编辑,第七、八届全国政协委员。
  • 桓济桓济,字仲道,谯国龙亢(今安徽怀远)人。东晋大司马桓温之次子。桓熙之弟,桓歆、桓祎、桓伟、桓玄之兄。桓济是东晋权臣桓温之次子,于晋穆帝升平四年(360年)时封为临贺县公,后来官至给
  • 郭符甲郭符甲(1605年-1648年),字辅伯,号介庵,福建省泉州府晋江县(今晋江市)人,明末官员,进士出身。郭符甲少年博学,弱冠时即为孝廉,然而此后六次春闱落第,十八载寒窗,直到崇祯十六年(1643年)方考中
  • 超越 (游戏)《超越》(日语:アウトラン,英语:Out Run)是一款由世嘉公司制作和发行的竞速游戏。本游戏最初于1986年在街机上发行。 后移植至多个平台《Out Run》是一款3D竞速游戏。游戏非常受