多体微扰理论

✍ dations ◷ 2025-05-19 09:02:01 #量子化学,微扰理论

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

相关

  • 后转录修饰转录后修饰(RNA修饰,或称修饰RNA)是真核细胞中,将初级转录RNA转化为成熟RNA的加工过程。一个很好的例子就是前mRNA转化为成熟的mRNA,其中包括剪接,并发生在蛋白质生物合成之前。这
  • 伊斯梅尔·本·艾哈迈德阿布·易卜拉欣·伊斯梅尔·本·艾哈迈德(波斯语:ابو ابراهیم اسماعیل بن احمد سامانی‎,849年-907年11月),通常被称为伊斯梅尔·本·艾哈迈德、伊
  • 博恩哈德三世 (萨克森-迈宁根)博恩哈德三世·弗雷德里克·威廉·阿尔布雷希特·格奥尔格(Bernhard III Frederick Wilhelm Albrecht Georg,1851年4月1日-1928年1月16日),出生于迈宁根。萨克森-迈宁根末代公爵
  • 亚历山大·伊万诺维奇·乌斯宾斯基亚历山大·伊万诺维奇·乌斯宾斯基(俄语:Алекса́ндр Ива́нович Успе́нский,1902年2月27日-1940年1月27日)是西西伯利亚地区内务人民委员部副部长、
  • 王宗炎王宗炎(1755年-1826年),字以除,号谷塍,又号晚闻居士,浙江萧山人,进士出身。乾隆四十五年(1780年)登庚子科进士,未授官而归。学问渊博,通籍后,杜门不出,筑有十万卷楼藏书,以文史自娱,尤好诱掖
  • 月宫一号月宫一号(英语:Lunar Palace 1、Moon Palace 1或Yuegong-1)是中华人民共和国一个为建设月球基地而设立的生命保障研究实验装置。月宫一号位于北京航空航天大学学院路校区内,占地
  • 阿利·贝恩阿利·贝恩(英语:Ari Behn,1972年9月30日-2019年12月25日),挪威作家。他1999年出版的短篇小说集Trist som faen卖出100,000本书并得到不少好评。2002年5月24日,贝恩和挪威国王哈拉
  • 奥耳波特计量表奥耳波特计量表,亦称作“奥耳波特偏见计量表”,是一种用来量度社会上偏见的程度的记量方式,由美国心理学家高尔顿·威拉德·奥尔波特的著作《偏见的本质》中衍生而来 。奥耳波
  • 老雷斯的故事《老雷斯的故事》(英语: 或 )是一部2012年美国计算机动画3D音乐喜剧片。根据苏斯博士发表于1971年的同名儿童文学书籍 改编而成。是照明娱乐出品的首部IMAX 3D电影,环球影业负
  • 黄雨蒙黄雨蒙,女,又名黄木七,四川人。据黄雨蒙亲属称,其于2020年7月5日从南京乘火车到格尔木,7日和一朋友联系后失踪,次日学校向警方报案。据警方通报,她曾乘出租车到达可可西里索南达杰