量子化学复习提纲

量子化学复习提纲

1 量子力学QM与分子力学MM

研究化合物的结构和性能，可以从两个层次：宏观层次和微观层次。宏观层次就是我们熟悉的实验观测；微观则是计算分子结构。计算分子结构基本上可以分为两种方法：量子力学（Quantum Mechanics）和分子力学(Molecular Mechanics)

量子力学是计算电子波函数的，计算电子波函数就要解薛定谔方程，简称S方程。通过核和电子的相互作用原理和规律等，运用量子力学原理，经过近似处理直接求解S方程得到电子波函数，从而求得电子总能和分子结构，我们称为第一性原理（First Principle），狭义的第一性原理计算就是从头算（ab initio）,它不采用经验参数，只用电子质量，光速，质子中子质量等少数实验数据去做量子计算。所以说量子力学是纯净无污染的（自己种的苹果）。

分子力学则从另一个方面，其直接计算原子和分子，不计算电子波函数，又叫力场方法（force field method）。它应用的原理是能量最小值方法。即原子间相互作用势下，改变原子（分子）的几何分布，以能量最小为判据，求得体系最佳构型。分子立场方法来源不清，假象的（超市卖的苹果）。

2 薛定谔方程

薛定谔方程HΨ=EΨ。薛定谔方程可以由驻波方程去推导。因为电子本身既是粒子又是波，而驻波的方程恰好可以反应粒子的性质。将驻波方程求二阶微导并与波粒二象性公式联系即可得到S方程。S方程不含自旋，而一个完整的电子波函数应该包括空间波函数和自旋波函数。S方程得到的波函数Ψ是不含自旋的空间波函数。既包括空间波函数也包括自旋波函数的方程是狄拉克方程D方程。

3 能量E

能量包括两部分：动能和势能。S方程中H代表哈密顿算符。H与E对应，所以H应该包括两个算符动能算符+势能算符。动能算符就是拉布拉斯算符（倒三角）作用于(-h2/8π2m)

一个完整的哈密顿算符在具体的电子计算中包括五个部分：电子动能（-）+核动能（-）+电子核吸引势能（-）+电子之间排斥能（+）+核之间排斥能（+）,括号内代表值的正负

实际计算中我们都采用B-O近似，也就是玻恩奥本海默近似(绝热近似)。与电子相比，核质量大，基本不动，所以我们把核动能+核之间排斥势视为常数，所以B-O近似下的哈密顿算符就只包括三项。这样就把薛定谔方程分为了两部分乘积：核S方程X 电子S方程。

4 结构优化

结构优化就是dE/dx=0的过程

5 全同粒子

所谓全同粒子是指质量电荷自旋等固有属性都相同的粒子。如所有的电子是全同粒子。全同粒子构成的体系特点是粒子无法区分，因为属性都一样。全同粒子可以分为两类：费米子和波色子。电子、质子、中子等是费米子，特点是自旋是半整数（如电子是

1/2）,费米子构成的全同粒子体系，如果交换任意两个粒子的位置，其波函数变号。

6 Slater行列式

基于费米子全同粒子体系交换位置变号的特点，可以用矩阵的形式表示出来。因为交换矩阵任意两行，矩阵值变负。倘若有任意两行相同，矩阵值为零（Pauli不相容原理，因为矩阵值为零，所以不允许两个完全一样属性的粒子占据同一个轨道）。这样的一个矩阵我们称之为Slater 行列式。矩阵前系数是1/√N！。矩阵中的每一项都是一个分子轨道。

7 周期性边界条件（Periodic Boundary Conditions,PBC）

一种边界性条件，在这种边界性条件下，当有一个粒子跑出我们设定的盒子时，必有一个粒子（称之为镜像）跑进来，从而维持系统粒子的稳定。周期性边界条件是分子动力学模拟的基本条件。整个体系就是这个盒子的无限三维扩展。（Born-von-Karman边界条件，从头算分子动力学模拟AIMD）

8 计算的基本理念

要根据问题选软件，而不是根据软件选择问题。

计算结果不是越接近实验值就是正确的。一个错误的方法可以得到更接近实验值的数值，但方法是错误的，值是没意义的。

计算的好不好与理论水平不一定成正比。比如有时候用高水平的理论结算得到的结果还不如低水平的理论计算得到的结果。但一般来说选择大基组会提供精确度。

精确度与cost关系。一般来说，牺牲精确度会提高效率。算得快的方法大部分都是以牺牲精确度为前提的。

9 原子单位

量化计算采用的单位制，me=1,e=1,h/2π=1,距离采用波尔半径a0=0.529埃，光速C=137，能量单位为Hartree,1Hartree=27.21ev=627.5095kcal/mol=2625.51kJ/mol

10 势能面（Potential energy surface ,PEC）

以能量为纵坐标，分子键长、键角等为横坐标，绘图，得到势能与分子结构参数的函数曲面，曲面上一阶导数为零的点我们称之为驻点（stationary point），最小的驻点对应于分子的平衡构型（equilibrium geometry），1级鞍点(saddle point)对应于过渡态。势能曲面上的最小能量途径(minimum energy path)对应于最可能的反应路径。(反应路径很多，最小能量路径是反应发生的几率最大的路径，但不是唯一路径)

11 基组(basis set)

首先介绍LCAO（linear combination of atomic orbitals）,原子轨道线性组合。我们知道薛定谔方程只能对一些简单体系，如类氢粒子等有精确解，当处理多电子原子时，我们不得不采用LCAO的方法。我们用一部分原子轨道的线性组合来表示分子轨道。选用的这些原子轨道，我们称之为基函数。基函数分为两部分：径向部分和角度部分。我们熟知的STO和GTO在径向部分表达式方面是不同的。

STO，称为Slater 函数，GTO称为高斯型函数。STO数值逼近真实值，但是收敛慢，所以计算速度慢；GTO函数收敛快，计算速度快，但是结果偏离真实值。因此人们结合上面两个函数特点，用多个GTO去模拟STO，同时达到计算精度较高和计算速度较快的目的

基组分类很多，王老师讲课中主要涉及到了三种基组分类。Pople基组和相关一致性CC（correlation consistent）基组、赝势基。

Pople基组。最小基组：就是填满电子所需的最少轨道。STO-3G是最小基组。分裂基组：处理价层轨道时，我们对价层轨道进行分层，分成两层，叫做双ζ基组。比如6-31G，表示对内层核心轨道采用6个GTO，对价层轨道分成两层，内层用3个GTO，外层用1个GTO。如果分成三层，叫做三ζ基组，比如6-311G，表示外层轨道分成三层，分别用3、1、1个GTO表示。极化：6-31G（d,p），d表示对重原子加d极化，p表示对氢原子加p极化。弥散：6-31++G,第一个+表示对重原子加弥散，第二个+表示对氢原子+弥散

CC基组：相关一致性基组考虑了价电子的相关能量，CC-pVDZ、CC-pVTZ、基组可以通过添加弥散函数的方式增大，符号表示为aug-，以aug-CC-pVDZ 为例，需要在原函数基础上加上1s,1p和1d函数

ECP(effective core potential):又称赝势。ECP是经验函数，哪个更精确需通过实验验证，，赝势没有结点。为什么引入赝势呢？我们知道在计算的时候，很多内核电子是不参与化学反应的，这些内核电子又要耗去我们大量的基函数描述。所以引入赝势的概念，ECP代表所有的内核电子，这是一种半经验的方法。内核电子通过合适的函数模型化，从而只考虑价层电子。

12 理论水平