多体微扰理论

✍ dations ◷ 2025-08-03 20:17:52 #量子化学,微扰理论

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

相关

  • 免疫分析免疫分析是一种利用抗体或抗原对某种溶液中的大分子或小分子进行测量的生物化学实验技术。免疫分析既可以定性确认某种大分子或小分子是否存在于受测溶液中,也可以定量测量某
  • 寄生曲霉寄生曲霉(学名:Aspergillus parasiticus)是属于散囊菌目发菌科曲霉属的一种真菌,可生长在土壤、稻谷等基物上。该种分布于中国、美国、阿根廷、巴西、荷兰、印度、印度尼西亚、
  • 国际复兴开发银行国际复兴开发银行(International Bank for Reconstruction and Development,简称IBRD)是一所国际金融机构。它给与中等收入的发展中国家提供贷款。国际复兴开发银行是世界银行
  • 菲利普·弗兰兹·冯·西博尔德菲利普·弗兰兹·冯·西博尔德(Philipp Franz von Siebold,1796年2月17日-1866年10月18日),德国内科医生、植物学家、旅行家、日本学家和日本器物收藏家。他是日本第一个女医楠本
  • 南康市南康区是中华人民共和国江西省的一个市辖区,现是隶属于赣州市。南康区地处于江西省南部,距离赣州市中心城区约半小时车程,因“地接岭南,人安物阜”而得名。原南康市总面积1800平
  • 别墅喜剧别墅喜剧 (意大利语:Villa Commedia)是小普林尼拥有的罗马别墅,它在意大利北部科摩湖的岸边。普林尼在科莫湖有几栋别墅,但是他在给朋友小普林尼 (立陶宛语:Plinio)的信中提到,“别
  • 热莉卡·茨维亚诺维奇热莉卡·茨维亚诺维奇(Жељка Цвијановић,1967年8月2日-),波黑政治人物,现任塞族共和国总统。曾在2013年3月12日至2018年11月19日,担任波黑塞族共和国总理。1967年,茨
  • 卢克博士卢卡什·塞巴斯蒂安·戈特瓦尔德(英语:Lukasz Sebastian Gottwald;1973年9月26日-),以艺名卢克博士(Dr. Luke)知名,是一名美国词曲创作者和唱片制作人。他的音乐职业生涯始于1997年作
  • 埃布尔·P·厄普舍埃布尔·帕克·厄普舍(英语:Abel Paker Upshur,1790年6月17日 - 1844年2月28日),美国律师、政治家,曾任美国海军部长和美国国务卿。
  • 分凝分凝,又称偏析,是指含有杂质的晶态物质熔化后再结晶时,杂质在结晶的固体和为洁净的液体中的浓度不同的现象。分凝是区熔提纯的理论依据,是提纯半导体的方法。