从头计算法(ab initio method)在量子化学计算中指基于量子力学基本原理直接求解 薛定谔方程 的量子化学计算方法。从头计算法的特点是没有经验参数,并且对体系不作过多的简化。对各种不同的 化学 体系采用基本相同的方法进行计算。目前的从头计算法包括基于哈特里-福克 方程 的哈特里–福克方法、在哈特里–福克基础上引入电子相关作用校正而发展起来的后哈特里–福克方法,以及多组态多参考态方法等。与半经验方法相比,从头计算法精度高,但耗时长。
算法介绍
从头计算法中引入的近似很少,主要包括
• 单电子近似
• 非 相对论 近似
对于具体的方法而言,可能用到其中的一种或几种近似。
玻恩– 奥本海默近似
玻恩–奥本海默近似 (Born- Oppenheimer approximation,简称 BO近似 ,又称 绝热近似 )是一种普遍使用的解包含电子与 原子核 的体系的量子力学 方程 的近似方法。
在用量子力学处理分子或其他体系时,需要通过解 埃尔温·薛定谔 方程或其他类似的 偏微分方程 获得体系 波函数 。这个过程往往由于体系自由度过多而非常困难,甚至无法进行。据玻恩-奥本海默近似中,考虑到原子核的质量要比电子大很多,一般要大3-4个数量级,因而在同样的相互作用下,电子的移动速度会较原子核快很多,这一速度的差异的结果是使得电子在每一时刻仿佛运动在静止原子核构成的势场中,而原子核则感受不到电子的具体位置,而只能受到平均作用力。由此,可以实现 原子核 坐标与电子坐标的近似变量分离,将求解整个体系的 波函数 的复杂过程分解为求解电子波函数和求解原子核波函数两个相对简单得多的过程。
在玻恩-奥本海默近似下,体系 波函数 可以被写为电子波函数与原子核波函数的乘积
玻恩-奥本海默近似由于在大多数情况下非常精确,又极大地降低了量子力学处理的难度,被广泛应用于分子结构研究、凝聚体物理学、量子化学、化学反应动力学等领域。
马克斯·玻恩 – 罗伯特·奥本海默 近似是由物理学家奥本海默与其导师玻恩共同提出的。
首先将体系的哈密顿 算符 分解为 原子核 动能算符与电子哈密顿算符两项
其中电子哈密顿算符包含所有 静电 相互作用(电子-电子,电子-核及核-核相互作用),以及电子的动能算符。
接下来,解以下的电子 薛定谔方程 以获得电子 波函数 和绝热势能面:
这一步一般被称为电子结构计算,由于电子的数量往往远大于 原子核 ,这一步往往是整个计算最昂贵的一步。根据计算能力,体系复杂程度和精度的需要,可以采用各种从头计算法直接近似求解,也可以采用 密度泛函理论 ,半经验方法,分子力学势场等一系列方法获得势能面。
获得势能面之后,可以用其求解 原子核 的 定态 薛定谔方程 得到原子核波的定态函数和体系能量
或通过求解含时 薛定谔方程 模拟 原子核 波函数 随时间的演化
这两种方法获得的结果可以互相转化。原子核定态波函数描述了分子的振动,一般用来获得振动零振动点能,进行红外、 微波 光谱或电子光谱振动与转动结构的计算;而波函数随时间的演变则描述了化学反应的 动力学 过程。
注意到第二步 原子核 薛定谔方程 的求解中,需要电子态能量以原子核位置为坐标的函数,因而往往需进行多次,乃至数千万次电子结构的计算,造成计算严重的困难。即使在小分子体系中,往往也无法采用完全量子理学的方法求解 波函数 。由于 原子核 的质量比电子大得多,量子效应不如电子明显,因此这一步往往采用比电子结构步骤更高程度的近似方法,以减轻计算负担。例如常采用经典力学或者半经典近似处理 原子核 运动;在振动幅度较小的情况下,可将势能面近似为二次曲面,而原子核的运动则转化为量子 谐振子 问题而可解析解出。某些情况下亦可采用平衡态的近似,运用统计力学或过渡态理论等统计方法则,只需要在少量关键位置进行电子结构和振动计算就得到关于分子的分布或反应 速率 等信息。
适用性姐介绍
玻恩-奥本海默近似只有在所在电子态与其他电子态能量都足够分离的情况下才有效。当电子态出现交叉或者接近时,玻恩-奥本海默近似即失效。
最低值原理
对于 基态 ,不加其他近似的情况下通过玻恩-奥本海默近似得出的体系总能量一定小于体系真实能量,因而给出了真实能量的下限。与此相对,另一绝热近似方法 马克斯·玻恩 黄近似则给出体系真实能量的上限。
哈特里-福克 方程
哈特里-福克 方程( 英语 :Hartree–Fock equation),又称为 HF方程 ,是一个应用 变分法 计算多电子系统 波函数 的方程,是量子物理、凝聚态物理学、量子化学中最重要的方程之一。HF方程形式上是单电子本征方程,求得的本征态是单电子波函数,即 分子轨道 。以HF方程为核心的数值计算方法称为“哈特里-福克方法”(Hartree–Fock method)。
基于 分子轨道理论 的所有量子化学计算方法都是以HF方法为基础的。鉴于分子轨道理论在现代量子化学中的广泛应用,HF 方程 被视为现代量子化学的基石。
哈特里-福克近似
哈特里-福克近似 也称为 分子轨道近似 或 单 行列式 近似 ,认为多电子体系 波函数 可以由体系分子轨道波函数构造的单个斯莱特行列式表示:
这里 为电子数, 是能量最低的 个 分子轨道 波函数 。而分子轨道波函数是由体系中所有 原子轨道 波函数经过线性组合构成的,
通常原子轨道波函数 由计算时指定的基组给出。
在此近似下,由 变分法 ,体系的多电子 薛定谔方程 可化为单电子方程
其中 为福克 算符 , 为体系的 分子轨道 , 是分子轨道 的轨道能。这就是 哈特里-福克方程 。求解该方程,即求解使方程成立的一组分子轨道;这组分子轨道使得在HF近似下的体系总能量最低。
发展简史
1927年,物理学家沃特·海特勒和 弗里茨·伦敦 完成了 氢气 分子的量子力学计算之后,开启了量子化学的时代。从此,人们便开始尝试使用量子力学理论来解释 化学 物质结构和化学现象。
为了解决多电子体系 薛定谔方程 近似求解的问题,量子化学家道格拉斯·哈特里在1928年提出了哈特里假设,将每个电子看做是在其他所有电子构成的平均势场中运动的粒子,并且首先提出了迭代法的思路。哈特里根据他的假设,将体系电子哈密顿算子分解为若干个单电子哈密顿算子的简单代数和,每个单电子哈密顿算子中只包含一个电子的坐标,因而体系多电子 波函数 可以表示为单电子波函数的简单乘积,这就是 哈特里 方程 。但是由于哈特里没有考虑电子波函数的反对称要求,他的哈特里方程实际上是非常不成功的。
1930年,哈特里的学生弗拉基米尔·福克和约翰·斯莱特分别提出了考虑泡利原理的自洽场迭代方程和单行列式型多电子体系波函数,这就是 哈特里-福克方程 。但是由于计算上的困难,HF方程诞生后整整沉寂了二十年,在1950年,量子化学家克莱门斯·罗特汉想到将 分子轨道 用 原子轨道 的线性组合来近似展开,而得到了闭壳层结构的罗特汉方程。
1953年, 美国 的R.帕里瑟、R.帕尔和 英国 的 约翰·波普 花费两年时间使用手摇计算器分别独立地实现了对 氮气 分子的RHF自洽场计算,这是人类首次通过求解HF 方程 获得对化学结构的量子力学解释,也是量子化学计算方法第一次实际完成。在第一次成功之后,伴随着电脑技术的迅猛发展,HF方程与量子化学一道获得长足发展,在HF方程的基础上,人们发展出了高级量子化学计算方法,使得计算精度进一步提高,通过对HF方程电子积分的简化和参数化,人们大大缩减了量子化学的计算量,使得对超过1000个原子的中等大小分子的计算成为可能。
相关条目
• 非绝热耦合
• 薛定谔方程
• 马克斯·玻恩 黄近似
• 玻恩 定则