第一性原理计算方法讲义

第一性原理计算原理和方法(总40页)-CAL-FENGHAI.-(YICAI)-Company One1-CAL-本页仅作为文档封面，使用请直接删除第二章 计算方法及其基本原理介绍化学反应的本质是旧键的断裂和新建的形成，参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。

因此阐述化学键的理论应当描写电子壳层的相互作用与重排，借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学，为解决这一问题提供了一般的方法，然而，对于一些实际的体系，不引入一些近似，就不可能求解其Schrodinger 方程。

这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。

这些近似和关于分子波函数的方程形成计算量子化学的数学基础。

2.1 SCF-MO 方法的基本原理分子轨道的自洽场计算方法(SCF-MO)是各种计算方法的理论基础和核心部分，因此在介绍本文计算工作所用方法之前，有必要对其关键的部分作一简要阐述。

2.1.1 Schrodinger 方程及一些基本近似 为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便，这里将主要阐明用于本文量子化学计算的一些重要的基本近似，给出SCF MO 方法的一些基本方程，并对这些方程作简略说明，因为在大量的文献和教材中对这些方程已有系统的推导和阐述[1-5]。

确定任何一个分子的可能稳定状态的电子结构和性质，在非相对论近似下，须求解定态Schrodinger 方程 ''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =⎥⎥⎦⎤⎢⎢⎣⎡-++∇-∇-∑∑∑∑∑∑≠≠ （2.1）其中分子波函数依赖于电子和原子核的坐标，Hamilton 算符包含了电子p 的动能和电子p 与q 的静电排斥算符，图2-1分子体系的坐标∑∑≠+∇-=p q p pq p e r H 12121ˆ2 (2.2) 以及原子核的动能∑∇-=A A AN M H 2121ˆ (2.3) 和电子与核的相互作用及核排斥能 ∑∑≠+-=p A B A AB B A pAA eN R Z Z r Z H ,21ˆ (2.4) 式中Z A 和M A 是原子核A 的电荷和质量，r pq =|r p -r q |，r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离（均以原子单位表示之）。

第一性原理计算方法讲义标准化管理部编码-[99968T-6889628-J68568-1689N]第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中，求解的过程中需要知道一些物理参量，如温度场方程中的热传导系数和浓度场方程中的扩散系数等，这些参量随着材料的不同而改变，需要通过实验或经验来确定，所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等，而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学，即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现，使得人们对于物质微观结构的认识日益深入。

原则上，量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学，使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题，引发了通讯技术和计算机技术的重大变革。

目前，结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有～1023数量级粒子的多粒子系统，具体应用量子理论时会导致物理方程过于复杂以至于无法求解，所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似（Born-Oppenheimei近似）将电子的运动和原子核的运动分开，从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数（分子轨道）为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量；对于研究分子体系，他可以作为一个很好的出发点，但是不适于研究固态体系。

1964年，Hohenberg和Kohn提出了严格的密度泛函理论（DensityFunctional Theory, DFT）。

第一性原理计算原理和方法精编Document number：WTT-LKK-GBB-08921-EIGG-22986第二章 计算方法及其基本原理介绍化学反应的本质是旧键的断裂和新建的形成，参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。

因此阐述化学键的理论应当描写电子壳层的相互作用与重排，借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学，为解决这一问题提供了一般的方法，然而，对于一些实际的体系，不引入一些近似，就不可能求解其Schrodinger 方程。

这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。

这些近似和关于分子波函数的方程形成计算量子化学的数学基础。

SCF-MO 方法的基本原理分子轨道的自洽场计算方法(SCF-MO)是各种计算方法的理论基础和核心部分，因此在介绍本文计算工作所用方法之前，有必要对其关键的部分作一简要阐述。

Schrodinger 方程及一些基本近似为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便，这里将主要阐明用于本文量子化学计算的一些重要的基本R AB =R 图2-1分子体系的坐标近似，给出SCF MO 方法的一些基本方程，并对这些方程作简略说明，因为在大量的文献和教材中对这些方程已有系统的推导和阐述[1-5]。

确定任何一个分子的可能稳定状态的电子结构和性质，在非相对论近似下，须求解定态Schrodinger 方程 ''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =⎥⎥⎦⎤⎢⎢⎣⎡-++∇-∇-∑∑∑∑∑∑≠≠ （）其中分子波函数依赖于电子和原子核的坐标，Hamilton 算符包含了电子p 的动能和电子p 与q 的静电排斥算符，∑∑≠+∇-=p q p pq p e r H 12121ˆ2以及原子核的动能∑∇-=A A AN M H 2121ˆ 和电子与核的相互作用及核排斥能∑∑≠+-=p A B A AB B A pA A eN R Z Z r Z H ,21ˆ式中Z A 和M A 是原子核A 的电荷和质量，r pq =|r p -r q |，r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离（均以原子单位表示之）。

第三节第一性原理计算简介第一性原理计算简介在物理学中,第一性原理计算或称从头计算是指,基于构建物理学的基础定理,不作任何假设,例如:经验模型和拟合参数,所进行的计算研究。

特别地,在凝聚态物理中,指的是运用薛定愕方程在一定的近似情况下,但不包括拟合实验数据所得到的参数和模型,对物质的电子结构进行计算r 从而得到所研究物质的性质的一种研究方法。

近些年,随着计算机技术的飞速发展,其运算能力越来越强大,使得人们可以处理更庞大更繁杂的物质结构体系,同时也使得计算物理成为了现代物理学,尤其是在凝聚态物理领域的一个重要分支。

众所周知,固体是由相对重且带正电的粒子——原子核,以及相对轻且带负电的粒子——电子聚集在一起构成的。

如果有个原子,需要处理的问题是包含有N+ZN(Z 为原子核所含的质子的个数)个粒子的电磁相互作用,是一个多体问题。

另一方面,由于处理的是微观粒子的运动,所以需要运用量子力学来描述其基本的运动规律和相互作用。

对于该系统,精确的多粒子哈密顿量可以写作:i 2i ii 1122R H M ?=--∑∑Fuuuuuuuuj其中位于為处的原子核的质量为M,.,位于巧处的电子的质量为m 一第一项是原子核的动能算符,第二项是电子的动能算符。

后三项分别是描述电子与原子核,单个电子与其它电子以及单个原子核与其它原子核之间的库伦相互作用。

很显然,直接精确求解(1.64)式几乎是不可能的。

为了在合理的近似条件下得到体系的本征值,需要作不同层次的近似。

1.3.1波恩-奥本海默(Bom-Oppenheimer)近似由于原子核的质量远大于电子质量,所以,原子核的运动速度远小于电子。

因此,可以将原子“冻结”在固定的位置,并假设电子在瞬时与原子核是平衡的。

或者说,只有电子在这个多体问题中是考察对象,原子核仅仅被当作一个带正电的外源场,相对于电子云是外在独立的。

该近似被称为波恩-奥本海默(Bom-Oppenheimer)近似。

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中，求解的过程中需要知道一些物理参量，如温度场方程中的热传导系数和浓度场方程中的扩散系数等，这些参量随着材料的不同而改变，需要通过实验或经验来确定，所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等，而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学，即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现，使得人们对于物质微观结构的认识日益深入。

原则上，量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学，使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题，引发了通讯技术和计算机技术的重大变革。

目前，结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有〜1023数量级粒子的多粒子系统，具体应用量子理论时会导致物理方程过于复杂以至于无法求解，所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似（Born-Oppenheimei 近似）将电子的运动和原子核的运动分开，从而将多粒子系统简化为多电子系统。

Hartree-Fock 近似将多电子问题简化为仅与以单电子波函数（分子轨道）为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量；对于研究分子体系，他可以作为一个很好的出发点，但是不适于研究固态体系。

1964年，Hohenberg和Kohn提出了严格的密度泛函理论（Density Functional Theory, DFT ）。

它建立在非均匀电子气理论基础之上，以粒子数密度（『）作为基本变量。

第一性原理计算方法第一性原理计算方法是一种基于量子力学原理的计算方法，它可以用来研究原子和分子的结构、性质和反应。

与传统的经验性方法相比，第一性原理计算方法具有更高的精度和可靠性，能够提供更多的物理和化学信息。

本文将介绍第一性原理计算方法的基本原理和应用。

首先，第一性原理计算方法是建立在薛定谔方程的基础上的。

薛定谔方程描述了体系的波函数随时间的演化，通过求解薛定谔方程，我们可以得到体系的能量、波函数和其他物理性质。

在第一性原理计算中，我们通常采用密度泛函理论来近似求解薛定谔方程，通过求解库仑势和交换-相关势的作用，得到体系的基态能量和波函数。

其次，第一性原理计算方法的应用非常广泛。

它可以用来研究固体、液体和气体的结构和性质，预测材料的稳定相和晶体结构，计算分子的几何构型和振动频率，分析化学反应的动力学过程等。

同时，第一性原理计算方法还可以用来设计新型的功能材料，优化催化剂的性能，预测分子的电子结构和光学性质，研究纳米材料的电子输运行为等。

在第一性原理计算方法的发展过程中，科学家们提出了许多不同的计算框架和方法，如密度泛函理论、量子蒙特卡洛方法、格林函数方法等。

这些方法在不同的体系和问题上都有各自的优势和局限性，需要根据具体的研究目的来选择合适的方法。

总的来说，第一性原理计算方法是一种强大的工具，它在材料科学、物理化学、生物化学等领域都有重要的应用价值。

随着计算机硬件和软件的不断发展，第一性原理计算方法将会变得更加高效和精确，为科学研究和工程应用提供更多的支持和帮助。

通过以上介绍，我们可以看到第一性原理计算方法在材料科学和化学领域的重要性和广泛应用。

它不仅可以帮助我们理解物质的基本性质，还可以指导新材料的设计和合成，促进科学技术的发展和进步。

因此，掌握和应用第一性原理计算方法对于科研工作者和工程技术人员来说都是非常重要的。

希望本文的介绍能够为读者提供一些有益的信息，引起对第一性原理计算方法的兴趣和关注。

第一性原理计算方法讲义 This manuscript was revised on November 28, 2020第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中，求解的过程中需要知道一些物理参量，如温度场方程中的热传导系数和浓度场方程中的扩散系数等，这些参量随着材料的不同而改变，需要通过实验或经验来确定，所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等，而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学，即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现，使得人们对于物质微观结构的认识日益深入。

原则上，量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学，使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题，引发了通讯技术和计算机技术的重大变革。

目前，结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有～1023数量级粒子的多粒子系统，具体应用量子理论时会导致物理方程过于复杂以至于无法求解，所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似（Born-Oppenheimei近似）将电子的运动和原子核的运动分开，从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数（分子轨道）为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量；对于研究分子体系，他可以作为一个很好的出发点，但是不适于研究固态体系。

1964年，Hohenberg和Kohn提出了严格的密度泛函理论（DensityFunctional Theory, DFT）。

第一性原理计算的原理和应用随着计算机技术的不断发展和物理化学科学的深入研究，人们发现可以使用计算机模拟复杂的现象和过程，这就是第一性原理计算。

本文将介绍第一性原理计算的原理和应用。

一、第一性原理计算的原理所谓第一性原理计算，是指基于量子力学的原理和公式推导出固体、液体和气体内部物理化学现象的计算方法。

其中最基本的公式是薛定谔方程式：HΨ = EΨ其中H是系统的哈密顿算符，Ψ是波函数，E是系统状态的能量。

这个方程可用来计算电子运动的态函数和能量。

但这个方程式无法直接解出来，因为它涉及到太多的变量。

因此，研究者们发明了一种数值算法，称为密度泛函理论（DFT）。

密度泛函理论中的密度泛函表述的是体系中全部粒子的费米分布函数，它是电子密度的函数。

通过求解密度泛函，就可以推算出化学反应、材料表面的反应、气态中的自由基反应等等。

二、第一性原理计算的应用第一性原理计算是基于量子力学的计算方法，也可以称为第一原理分析计算。

它可以帮助我们理解物理和化学的基本原理，对于材料和化学的设计也有很大帮助。

1、材料设计组成纳米和宏观物质的原子是复杂的物理系统，它们的内部结构和外部特性带有很多未知因素。

第一性原理计算可以让我们更好地理解原子和分子之间的物理作用原理，通过模拟构建物质结构，预测材料的性质，帮助科学家们设计新的材料。

2、化学反应在化学反应中，基本的机理是原子之间的结构、强度和电性互相作用并且相互作用引入新的物质。

为了利用化学反应进行新的合成，我们需要在原子和分子层面上理解化学反应机理。

第一性原理计算可以揭示反应的原则，为我们提供了在计算机上模拟和预测化学反应的能力。

3、超导研究超导指的是电流在特定材料中不受电阻的限制传导。

探索超导的机制和原理，以及发现可以用此技术制造的材料，可以为能源和电子技术领域带来重大发展机会。

第一性原理计算是超导研究中必不可少的工具，可以预测和评估新材料的超导行为。

三、结论第一性原理计算是一种计算复杂物理化学现象的方法。

模拟电子结构的第一性原理计算方法在物理学和化学领域，模拟电子结构的第一性原理计算方法已经成为研究原子、分子和固体材料的重要工具。

使用这种方法，科学家们能够计算出材料的电子结构，并从中推导出许多性质和行为。

本文将介绍第一性原理计算方法的基本概念、历史和应用，并探讨它的未来发展方向。

第一性原理计算方法是基于量子力学原理的计算方法，它可以预测原子、分子和材料的性质和行为。

与传统的计算方法不同，第一性原理计算方法只依赖于物理定律和区域性泛函理论，而不需要经验参数。

因此，它被认为是一种可靠的计算方法，能够提供最精确的结果。

第一性原理计算方法的历史可以追溯到20世纪50年代初期，当时物理学家们正在尝试计算原子和分子的性质。

他们首先解决了氢原子的波函数，并尝试将该方法扩展到更复杂的分子中。

随着计算机技术的发展，这种方法得到了更广泛的应用，并在化学和材料科学中取得了显著成果。

第一性原理计算方法的基础是库仑势能和基态波函数。

库仑势能是原子和分子之间的物理相互作用力，它可以用经典电磁学的定律计算。

基态波函数表示能量最低的电子状态，可以用量子力学的理论来计算。

使用这些基本概念，科学家们可以计算出原子、分子和材料的能带结构、电子密度等物理量，并用这些结果来预测材料的性质和行为。

第一性原理计算方法的应用非常广泛，涉及到许多领域。

在材料科学中，它被用于研究半导体、金属、陶瓷、催化剂等材料，探索材料的电学、磁学、光学、热学等性质。

在化学中，它被用于研究化学反应、催化剂设计、化学键合等问题。

它还被用于研究生物分子、纳米材料、太阳能电池等领域。

不过，第一性原理计算方法也存在一些问题和挑战。

首先，计算的复杂度非常高，耗时非常长。

其次，它对计算资源的要求非常高，需要大量的计算机处理能力。

第三，它对初始参数的依赖性非常强，需要对模型做出合理的假设才能得到可靠的结果。

因此，第一性原理计算方法的进一步发展需要依靠新的算法、理论和技术。

总的来说，第一性原理计算方法是一种非常强大的工具，能够帮助科学家们解决许多问题。

第一性原理计算方法引言前面讲述的有限元和有限差分等数值计算方法中，求解的过程中需要知道一些物理参量，如温度场方程中的热传导系数和浓度场方程中的扩散系数等，这些参量随着材料的不同而改变，需要通过实验或经验来确定，所以这些方法也叫做经验或者半经验方法。

而第一性原理计算方法只需要知道几个基本的物理参量如电子质量、电子的电量、原子的质量、原子的核电荷数、布朗克常数、波尔半径等，而不需要知道那些经验或半经验的参数。

第一性原理计算方法的理论基础是量子力学，即对体系薛定额方程的求解。

量子力学是反映微观粒子运动规律的理论。

量子力学的出现，使得人们对于物质微观结构的认识日益深入。

原则上，量子力学完全可以解释原子之间是如何相互作用从而构成固体的。

量子力学在物理、化学、材料、生物以及许多现代技术中得到了广泛的应用。

以量子力学为基础而发展起来的固体物理学，使人们搞清了“为什么物质有半导体、导体、绝缘体的区别”等一系列基本问题，引发了通讯技术和计算机技术的重大变革。

目前，结合高速发展的计算机技术建立起来的计算材料科学已经在材料设计、物性研究方面发挥着越来越重要的作用。

但是固体是具有～1023数量级粒子的多粒子系统，具体应用量子理论时会导致物理方程过于复杂以至于无法求解，所以将量子理论应用于固体系统必须采用一些近似和简化。

绝热近似（Born-Oppenheimei近似）将电子的运动和原子核的运动分开，从而将多粒子系统简化为多电子系统。

Hartree-Fock近似将多电子问题简化为仅与以单电子波函数（分子轨道）为基本变量的单粒子问题。

但是其中波函数的行列式表示使得求解需要非常大的计算量；对于研究分子体系，他可以作为一个很好的出发点，但是不适于研究固态体系。

1964年，Hohenberg和Kohn提出了严格的密度泛函理论（Density Functional Theory, DFT）。

它建立在非均匀电子气理论基础之上， 作为基本变量。

材料科学中的第一性原理计算第一性原理计算是材料科学研究中一种重要的计算方法。

它是基于量子力学理论和电子结构理论的计算模型，通过求解薛定谔方程，从基本粒子（原子、离子、电子）的特性出发，利用数学方法预测和描述材料的结构、能量、性质等基本信息。

本文将对第一性原理计算的原理、方法和应用进行详细介绍。

第一性原理计算的核心是量子力学。

量子力学是描述微观粒子行为的理论，它认为微观粒子的运动和相互作用需要用波函数描述，而波函数可以通过薛定谔方程求解。

在材料科学中，我们关心的是材料中电子的结构和性质。

通过解薛定谔方程，可以得到材料中电子的轨道分布、能带结构和电子密度等信息，进而预测和研究材料的各种性质。

第一性原理计算分为两个主要步骤：构建模型和求解薛定谔方程。

首先，需要确定材料的晶胞结构，即原子的排列方式和间距。

其次，需要选择合适的计算方法，如密度泛函理论（DFT）等。

DFT是一种基于电子密度的近似方法，它将材料中的电子相互作用简化为一个电子密度函数。

然后，需要选取计算所需的参数，包括平面波基组、能量截断和k点网格等。

最后，通过求解薛定谔方程，可以得到材料中电子的波函数和能量等信息。

第一性原理计算在材料科学中有广泛的应用。

首先，它可以用于材料的结构预测和优化。

通过计算不同原子和离子的结合能、晶格参数和局域构型能等信息，可以预测新材料的结构和稳定性，为材料设计和合成提供指导。

其次，第一性原理计算可以用于研究材料的电子性质。

通过计算材料的能带结构、禁带宽度和电子态密度等信息，可以预测材料的导电性和光学性质。

此外，第一性原理计算还可以用于模拟材料的机械性质、热学性质和磁学性质等。

尽管第一性原理计算有广泛的应用，但其存在一些限制。

首先，求解薛定谔方程是一项复杂且计算量大的任务，需要高性能计算机和大量的计算时间。

其次，第一性原理计算通常采用一些近似方法，如DFT等，会带来一定的误差。

此外，由于计算的复杂性，第一性原理计算通常只能研究小尺寸的体系，难以模拟大尺寸和复杂的材料。

第二章 计算方法及其基本原理介绍化学反应的本质是旧键的断裂和新建的形成，参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。

因此阐述化学键的理论应当描写电子壳层的相互作用与重排，借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学，为解决这一问题提供了一般的方法，然而，对于一些实际的体系，不引入一些近似，就不可能求解其Schrodinger 方程。

这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。

这些近似和关于分子波函数的方程形成计算量子化学的数学基础。

2.1 SCF-MO 方法的基本原理分子轨道的自洽场计算方法(SCF-MO)是各种计算方法的理论基础和核心部分，因此在介绍本文计算工作所用方法之前，有必要对其关键的部分作一简要阐述。

2.1.1 Schrodinger 方程及一些基本近似 为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便，这里将主要阐明用于本文量子化学计算的一些重要的基本近似，给出SCF MO 方法的一些基本方程，并对这些方程作简略说明，因为在大量的文献和教材中对这些方程已有系统的推导和阐述[1-5]。

确定任何一个分子的可能稳定状态的电子结构和性质，在非相对论近似下，须求R AB =R 图2-1分子体系的坐标解定态Schrodinger 方程''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =⎥⎥⎦⎤⎢⎢⎣⎡-++∇-∇-∑∑∑∑∑∑≠≠ （2.1） 其中分子波函数依赖于电子和原子核的坐标，Hamilton 算符包含了电子p 的动能和电子p 与q 的静电排斥算符，∑∑≠+∇-=p q p pqp e r H 12121ˆ2 (2.2) 以及原子核的动能∑∇-=A A AN M H 2121ˆ (2.3) 和电子与核的相互作用及核排斥能∑∑≠+-=p A B A AB B A pAA eN R Z Z r Z H ,21ˆ (2.4) 式中Z A 和M A 是原子核A 的电荷和质量，r pq =|r p -r q |，r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离（均以原子单位表示之）。