耦合簇方法

✍ dations ◷ 2025-10-23 08:41:05 #量子化学

耦合簇方法(coupled cluster, CC)是量子化学全始计算法中对多电子相关能的其中一种高精确计算方法。它从哈特里-福克分子轨道出发,通过指数形式的耦合算符运算得到真实体系的波函数。一些小分子和中等大小的分子精度最高的计算结果是通过 CC 方法得到的。

耦合簇方法提供了一种近似求解不含时薛定谔方程的方法:

这里 H ^ {\displaystyle {\hat {H}}} 来表示。耦合簇理论的其它变体,如运动方程耦合簇方法(英语:equation-of-motion coupled cluster) 和多参考态耦合簇方法(英语:multi-reference coupled cluster),则提供了求解体系激发态的方法。

体系的基态波函数可以用下面的拟设来表出:

式中 | Φ 0 {\displaystyle \vert {\Phi _{0}}\rangle } 表示占据轨道,而 表示空轨道。在耦合簇算符中的产生和湮没算符按照正规序排列。单粒子激发算符 T ^ 1 {\displaystyle {\hat {T}}_{1}} 和双粒子激发算符 T ^ 2 {\displaystyle {\hat {T}}_{2}} 分别把 | Φ 0 {\displaystyle \vert {\Phi _{0}}\rangle } 变为单激发和双激发斯莱特行列式的线性组合。为了最终得到体系的波函数,需要求解拟设中的待定系数 t i a {\displaystyle t_{i}^{a}} t i j a b {\displaystyle t_{ij}^{ab}} 等。

考虑到簇算符 T ^ {\displaystyle {\hat {T}}} 的结构后,指数耦合算符 e T ^ {\displaystyle e^{\hat {T}}} 可以展开成泰勒级数:

事实上,这一级数是有限的,因为分子轨道的数目与激发的数目都是有限的。为了简化求解系数 t {\displaystyle t} 的过程, T ^ {\displaystyle {\hat {T}}} 的展开式中一般在双激发或略高一点的激发处截断,很少有超过四激发的。这是因为是否包含五激发以上的算符 T ^ 5 {\displaystyle {\hat {T}}_{5}} T ^ 6 {\displaystyle {\hat {T}}_{6}} 等,对最终计算结果的影响很小。而且,即使只在簇算符的表达式中取前 n {\displaystyle n} 项:

那么由于耦合算符具有指数形式,高于 n {\displaystyle n} 激发的斯莱特行列式仍然会对最终的波函数有贡献。因此,在 T ^ n {\displaystyle {\hat {T}}_{n}} 处截断的 CC 方法通常能比激发数最高为 n {\displaystyle n} 的 CI 方法获得更多的电子相关能修正。

耦合簇方程就是展开系数 t {\displaystyle t} 所满足的方程。有多种方法来书写这一方程,其中标准的做法是会得到一个可以迭代求解的方程组。耦合簇方法的薛定谔方程可以写成:

假设现在共有 q {\displaystyle q} t {\displaystyle t} 系数需要求解。于是我们需要 q {\displaystyle q} 个方程。注意到每一个 t {\displaystyle t} 系数都与唯一的一个激发斯莱特行列式相关联: t i j k . . . a b c . . . {\displaystyle t_{ijk...}^{abc...}} 对应的是 | Φ 0 {\displaystyle \vert {\Phi _{0}}\rangle } 中处于 i , j , k , {\displaystyle i,j,k,\cdots } 轨道上的电子分别被激发到 a , b , c , {\displaystyle a,b,c,\cdots } 轨道上所得的行列式。上式两边向对应的行列式投影,就得到了我们所要的 q {\displaystyle q} 个方程。

式中 | Ψ {\displaystyle \vert {\Psi ^{*}}\rangle } 表示任意一个与待求的 t {\displaystyle t} 系数相关联的激发行列式。为了更好地利用这些方程之间的联系,我们可以把上面的方程改写成一种更方便的形式,将 e T ^ {\displaystyle e^{-{\hat {T}}}} 乘到耦合簇薛定谔方程两端,然后分别向 Ψ 0 {\displaystyle \Psi _{0}} Ψ {\displaystyle \Psi ^{*}} 投影,我们得到:

第一式提供了求解 CC 能量的方法,第二式则是用来求解 t {\displaystyle t} 系数的方程。以标准的 CCSD 方法为例,方程组中包括下面三组方程:

上式中经相似变换后的哈密顿量(用 H ¯ {\displaystyle {\bar {H}}} 表示)可以通过BCH 公式(英语:BCH formula) 求出:

H ¯ {\displaystyle {\bar {H}}} 不是厄米的。

传统上耦合簇方法依照 T ^ {\displaystyle {\hat {T}}} 中包含哪些 T ^ n {\displaystyle {\hat {T}}_{n}} 算符来进行分类。相应的方法名称则由 CC 后面加上相应的字母构成:

例如,CCSDT 方法里面簇算符 T ^ {\displaystyle {\hat {T}}} 的表达式如下:

在圆括号里面的项则表示它们是通过微扰理论求得的。例如 CCSD(T) 表示:

相关

  • 铜合金铜合金是以铜为主的合金,它们能有效地抵抗腐蚀。最多人知道的种类是青铜(铜为主要,锡为次要)和黄铜(铜为主要,锌为次要)。在古代欧洲“青铜”和“黄铜”术语常混用,所以现在欧洲(尤其
  • 朗伯·阿道夫·雅克·凯特勒凯特勒(Lambert Adolphe Jacques Quetelet,1796年-1874年),19世纪比利时的通才,他既是统计学家、又是数学家和天文学家。他从统计学角度出发看人,认为人的成长是会依从一套既定的法
  • 塔拉迪加塔拉迪加国家森林(英语:Talladega National Forest)位于美国亚拉巴马州阿巴拉契亚山脉南端,面积392,567英亩(158,866公顷)。在被联邦政府收购前,森林曾一度遭滥砍滥伐。目前,这片土
  • 2019冠状病毒病丹麦疫情2月2月3月3月4月4月过去15日过去15日2019冠状病毒病丹麦疫情,介绍在2019新型冠状病毒疫情中,在丹麦发生的情况。2020年2月27日,丹麦出现了首例确诊新冠肺炎病例。1月底到2月初
  • 詹姆斯·菲尼莫尔·库珀詹姆士·菲尼莫尔·库柏(英语:James Fenimore Cooper,1789年9月15日-1851年9月14日)是十九世纪最早赢得国际声誉的美国作家之一。也是浪漫主义的代表作家。他的主要作品是以猎人
  • 关心亚洲学者委员会关心亚洲学者委员会(英语:Committee of Concerned Asian Scholars)是一个由研究生和年轻大学教员于1968年创立的组织,旨在反对美国参与越南战争。该会于当年的亚洲研究协会费城
  • 兴威兴威(英语:Hsenwi),又称赛尼(英语:Hseni)或登尼(英语:Theinni,掸语:သဵၼ်ႈဝီ,缅甸语:သိန္နီ,缅甸语发音:,泰语:แสนหวี,皇家音译:,泰语发音:)缅甸掸邦腊戍县兴威镇区下辖的一个镇
  • 王子晏王子晏(1892年-1951年),名咸熙,浙江嘉兴人,中国近代围棋名家。他与顾水如并称为“南王北顾”。王子晏早年曾跟随方金题、吴祥麐学棋。1920年代初他移居上海。因受张澹如赏识,供职于
  • 种谔种谔(1017年-1083年),字子正,洛阳人,北宋大将。北宋名将种世衡第五子,世代守宋夏边境。以父亲的荫蔽入官,升到左藏库副使(从七品武官阶),延帅(鄜延路安抚使)陆诜推荐种谔任青涧城(城为设在
  • Q阶乘幂q阶乘幂是阶乘幂的Q-模拟。与阶乘幂在广义超几何函数中的作用类似,q阶乘幂也是定义基本超几何函数的基础。因为很多关于q阶乘幂的等式都含有多个q阶乘幂相乘,因此在标准写法中