多体微扰理论

✍ dations ◷ 2025-07-07 04:31:05 #量子化学,微扰理论

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

相关

  • 概率论概率论(英语:Probability theory)是集中研究概率及随机现象的数学分支,是研究随机性或不确定性等现象的数学。概率论主要研究对象为随机事件、随机变量以及随机过程。对于随机事
  • 努瓦克肖特努瓦克肖特(阿拉伯语:نواكشوط‎),毛里塔尼亚首都和最大城市,政治、文化、商业、金融中心。人口611,883人(2000年),约占全国人口的四分之一。位于撒哈拉沙漠西部,距大西洋海岸5
  • 机械电子工业部机械电子部是中华人民共和国国务院曾经存在的部门,1988年,原有的电子工业部与国家机械工业委员会合并成为机械电子工业部,1993年3月,将机械电子工业部又分成电子工业部和机械工
  • 高碘酸钠高碘酸钠是高碘酸的钠盐,有两种形式,偏高碘酸钠(化学式:NaIO4)和原高碘酸钠(化学式:Na5IO6),都可用作有机合成中的氧化剂。高碘酸钠密度为3.865g/cm3,可溶于水,加热时分解为碘酸钠(NaIO3
  • 袋貂属13种,见内文袋貂属()是属于袋貂科的有袋类动物。本属包括下列各物种:
  • 弗雷德里克王子 (荷兰)弗雷德里克王子(荷兰语:,1797年2月28日-1881年9月8日),奥兰治亲王威廉六世(1815年成为荷兰国王,称威廉一世)的次子,威廉二世的弟弟。1813年,随着拿破仑在莱比锡战役的失败,弗雷德里克随
  • 绿米绿米(越南语:Cốm)是一道越南甜品,味甜并附有坚果味,呈绿色。绿米的颜色不是染上去的,而是通过低温烘培未成熟的米粒,并用杵臼捣碎压平所致。绿米作为一道秋季的时令菜肴,在越南红河
  • 朝代 (2017年电视剧)《朝代》(英语:)是美国黄金时段肥皂剧,自2017年10月11日起于CW电视网首播。
  • 希维娅·邦戈·翁丁巴希维娅·邦戈·翁丁巴(英语:Sylvia Bongo Ondimba,1963年3月11日-),加蓬总统阿里·邦戈·翁丁巴的夫人。两人在1989年结为夫妻。他们育有四名子女。希维娅出生于法国巴黎,大约两岁
  • 史提夫·卡根史提夫·卡根(Steve Kagen;1949年12月12日-)是美国的一位政治人物。在2007年至2011年期间,他是威斯康星州第8选举区选出的美国众议院议员。他的党籍是民主党。卡根也是一位医师,主