多体微扰理论

✍ dations ◷ 2025-06-09 03:46:32 #量子化学,微扰理论

多体微扰理论是一种基于分子轨域理论的高级量子化学计算方法。这种方法以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方法也不能很好处理。

相关

  • 沧龙科沧龙超科(Mosasauroidea),Mosa在拉丁语意为荷兰的默兹河(Meuse river),sauros在希腊文意为蜥蜴。沧龙类是种如蛇般弯曲的海生爬行动物。第一个沧龙类化石在1780年于默兹河流域的马
  • 财政保守主义财政保守主义(Fiscal conservatism)是用来形容主张避免赤字开支(Deficit spending)的保守主义财政政策。财政保守主义者往往支持减少整体政府开支、赤字和国债,以达到最重要的平
  • 磨损作用 (牙科)磨损作用(abrasion),牙科上称为牙磨损,指受外来元素形成的机械力(mechanical forces),导致牙齿结构脱落。外来元素包括牙刷、牙签、牙线等。
  • 双璧传说《双璧传说》,是一部台湾偶像剧,共有15集,由中国电视公司所拍摄。由张天霖一人分饰二位男主角。家财万贯的方耀德(张天霖 饰)是知名财团方氏企业董事长方建山(冯粹帆 饰)的独子;从小
  • 腓特烈·弗朗茨三世腓特烈·弗朗茨三世(Friedrich Franz III,1851年3月19日-1897年4月10日),全名腓特烈·弗朗茨·保罗·尼库劳斯·恩斯特·亨利(Friedrich Franz Paul Nikolaus Ernst Heinrich),梅克
  • 核素图核素图或核素表是一个二维图表。其中一维表示核素的中子数,另一维表示它的质子数。因此图中的每一点代表某个元素真实存在的或者假想的核素。相比于更关注化学性质的元素周期
  • 九州铁道纪念馆九州铁道纪念馆(日语:九州鉄道記念館)是位于日本福冈县北九州市门司区清滝二丁目3番29号的一座以铁路为主题的博物馆。九州铁道纪念馆开业于2003年8月9日,靠近门司港站。2014年,
  • GOT YA! 公园少女《 GOT YA! 公园少女》(韩语:GOT YA! 공원소녀,英语:GOT YA! Gong Won So Nyeo)为韩国Mnet电视台2018年推出的综艺节目,由金信英、Dindin(朝鲜语:딘딘)、孙东云(HIGHLIGHT)、素珍(Girl's
  • 印尼电影资料馆印尼电影资料馆(印尼语:Sinematek Indonesia,简称Sinematek)是一座位于印度尼西亚(印尼)雅加达的电影资料馆。该馆于1975年由米斯巴赫·尤萨·比兰(英语:Misbach Yusa Biran)和阿斯鲁
  • 黄友谋黄友谋(1910年-1988年),男,广东梅县人,中国教育家、物理学家,九三学社社员,曾任广东省人大常委会副主任、暨南大学副校长、中山大学副校长。