量子化学-HFR方程

Hartree-Fock-Roothaan方程•Hartree方程•Hartree-Fock方程•Hartree-Fock-Roothaan方程•总结•Koopmans定理及若干概念Hartree方程1.n电子原子体系的Hamiltonian算符为：单电子算符：双电子算符：总电子Hamiltonian算符：ninjijniiniirrZH11112012121ˆkikirZih221ˆ22211,ˆjijijiijzzyyxxrjigjiijigihH,ˆˆˆ0(k为原子核编号)当体系只有一个电子时(即氢原子和类氢离子)，Schrödinger方程为：该方程可以严格求解。当体系有两个或两个以上电子时，Schrödinger方程为：由于双电子算符ĝ(i,j)中含有两个电子的坐标，难以分离变量，故此Schrödinger方程难以求解。ErZ221ErrZnininijiijii1112121将电子间的库仑排斥作用平均化，每个电子均视为在核库仑场与其它电子对该电子作用的平均势相叠加而成的势场中运动，从而单个电子的运动特性只取决于其它电子的平均密度分布（即电子云）而与后者的瞬时位置无关。于是，每个电子的状态可分别用一个单电子波函数描述。又因各单电子波函数的自变量彼此独立，N电子体系的波函数可写成N个单电子波函数的乘积非相对论近似Born-Oppenheimer(玻恩-奥本海默)近似单电子近似（轨道近似）nnnxxxxxx221121,,,自洽场(SCF)方法三个近似单电子近似的结果是将原来需要求解含有n个电子坐标的体系总波函数的问题拆分为求解n个单电子波函数的问题。后者由于自变量数目大幅度减少，所以在数学上的求解难度也相应地大为降低。Hartree自洽场方法在n电子原子中，电子i的势能算符为：j(j)是描述电子j的波函数，对电子j的所有空间可能的位置取平均可以得到电子i和电子j的平均排斥能：由于上式右边是对电子j的坐标积分，所以只与电子i坐标有关，这样其他(n-1)个电子对电子i的排斥算符（即Coulomb算符）则可表示为：nijijiirrZnV1,,2,1ˆjjijjijdjrjijV1*ijjjijjidjrjJ1ˆ*ijVij单电子Hamiltonian算符可改写为：电子i的单电子Schrödinger方程为：这就是著名的Hartree方程。要求解i(i)，必须先知道j(j)，这相对于要解这个方程必须先知道方程的解。为解决这个困难，可采用迭代法进行求解。iidjrjrZiiinijjjijjii121*22*11ˆ2niijjjjiiijZhjjdrr•先假定一组N个函数作为波函数的初值：•将这组初值代入Hartree方程，得到一组N个新的波函数作为一级近似波函数：•再用这组一级近似波函数代入Hartree方程，可得到二级近似波函数，重复这个过程，直到最后求解的单电子波函数、对应的轨道能量或总能量收敛为止。000201,,,,,Ni111211,,,,,Ni此时，求解Hatree方程得到的是一系列的本征值Ei和本征函数i，如：E1，E2，E3…(单电子能级）1，2，3...（单电子轨道，分子中为分子轨道）缺陷•1.未考虑电子的自旋•2.电子作为费米子，没有将其波函数的反对称性表现出来波函数的反对称性和Slater行列式a)电子属于费米子(Fermions，自旋量子数为半整数的微观粒子)。由费米子构成全同粒子体系服从费米统计规律，体系总波函数必须满足两个必要条件：反对称性和全同性。b)Hartree方程得到的波函数为：该波函数不满足多电子波函数的反对称性条件，并且没有考虑电子自旋问题。nnn21,,2,121NjiNijNjiPij,,,,,,1,,,,,,1,,,,,,1ˆ22ˆ1,,,,,,1,,,,,,ijPijNjiN波函数反对称性：任意交换两个粒子的全坐标后波函数的绝对值相同，符号相反。即：全同性：任意交换两个粒子的全坐标后粒子的总几率分布不变。即：NNNNNNN212121!1222111c)在轨道近似条件下，多电子体系反对称波函数的最简单表达形式是由分子轨道波函数构成的单Slater行列式：注：行列式的行号：标注轨道列号：标注电子d)用单Slater行列式表示的波函数能够满足费米子特性的要求：反对称性：交换两个电子的全坐标，相当于交换行列式的两列，行列式的数值不变但正负符号相反；全同性：每个电子在每个MO上出现的机会均等。e)人们将采用Slater行列式(考虑了自旋波函数)后的Hartree方程称为Hartree-Fock(HF)方程。解方程思路要求得到的分子轨道是正交归一化的！Hartree-Fock算符和方程：1.电子i的Hartree-Fock方程为：jiijjidjriK1ˆ**1ˆijjjijJjjdr2.单电子Fock算符表示为：iiKJhijnijji)]ˆˆ(ˆ[)ˆˆ(ˆˆˆˆjnijjiijiiKJhghF能量表达1.先求解单电子Fock算符的期望值：]ˆˆ[ˆˆ****iKiiJiihiiFiiiiiiijiiiiiiii其中*****1ˆ1iiiijjjiijijjijiijijiJiijjdirijjiddrJJij为电子i和电子j之间的Coulomb排斥双电子积分。ihihiiiiiˆ*ijijiiijjiijiijjiiiiKddijriiidjriiiKi11ˆ*****其中Kij为电子i和电子j之间的交换双电子积分。所以单电子Fock算符的期望值为：njijijiiiKJh1ninjijijniininjijijniiiKJKJhE11111102121ninjijijniininjijijniiiKJKJhE11111102121记Gi=0iii1111()()22nniiiiEhGGε对n个电子的能量求和，所得的能量并不是体系的总能量，因为这里对双电子积分部分重复计算了一次，所以体系总能量为：Hartree-Fock-Roothaan方程当每个分子轨道按某个基组集合展开,用有限展开项,按一定精度逼近分子轨道。这样,对分子轨道的变分就转变为对展开系数的变分。H-F方程就从一组非线性的微积分方程转化为一组数目有限的代数方程,只需迭代求解分子轨道组合系数，称为Hartree-Fock-Roothaan方程。哈特里－福克－罗特汉方程(一)、Roothaan方程的简单推导由前已得到H-F方程：iiiFˆ分子轨道i用基组（或原子轨道）进行展开misissc(s=1,2,…,m；m为原子轨道个数)代入方程得ˆsisisisssFcc左乘*r对整个空间积分得**ˆsirsisirssscFdcd令*ˆˆHFrsrsFFd;*rsrsSdjjjiKJhF2ˆ代入前式得ˆ()0HFsirsirsscFS（s=1,…,m;i=1,…,m）该方程实际是一组代数方程，把它具体写出来：1111111211211211111221221222222222221112221c()()...()02c()()...()0......c()()...()0mmmmmmmmmmmmmmmmmmmmmiFScFScFSiFScFScFSimFScFScFS时，εεε时，εεε时，εεε将上面线性方程可写成矩阵形式SCCFˆRoothaan方程（指线性方程）非零解的条件是下列久期行列式值为零：ˆdet0rsirsFS该方程为广义本征方程。其中ε相当于算符F的本征值,C相当于算符F属于本征值ε的本征向量。])1()1(ˆ)1()1()1(ˆ)1(2[)1()1(ˆ)1()1()1(ˆ)1(2/1sjrsjrnjsrsrrsKJhFF212**212)2()2(*)1()2()2(*)1()1()1(ˆdrccdrJuttuujtjsjjssj)|)2()2(*)1()1(*)1()1(ˆ)1(**2112**tursccddrccJtuujtjutsrtuujtjsjr双电子排斥积分同理：)|()1()1(ˆ)1(*tsruccKtuujtjsjr双电子交换积分ˆ(1)(1)(1)rsrshh单电子积分单电子哈密顿算符jnjc由于矩阵元Frs中包含分子轨道，即含未知的展开系数，所以只能用自洽场迭代的方法解Roothaan方程。2/1*2njujtjtuccP/21*[2(|)(|)]nrsrstjujtujFhccrsturuts定义密度矩阵元1||2rsrstutuFhrsturutsPFrs又可表达为：rsrsrsFhG1||2rstutuGrsturutsPn为电子数，n/2为占据轨道数为Hamilton矩阵，为电子排斥矩阵rshrsG（注：P为Hermite矩阵）FCSCεFCC'ACC总能获得酉阵A，使得I=A+SA如果基组是正交归一的如果基组不是正交归一的εSAC'AFAC'AεSAC'FAC''''CCF注：F(Fock矩阵)、S（重叠矩阵）均为为Hermite矩阵。对于Hermite矩阵，总能找到酉矩阵作酉变换，使其对角化解HFR方程：1.选择基组；2.求解hrs，Srs和(rs|tu)；3.利用重叠积分Srs，用正交化过程计算矩阵A；4.初始猜测系数cni：，并得到密度矩阵P；5.计算Fock矩阵元，Frs；6.计算；7.矩阵对角化，得到，C’；8.计算系数阵C=AC’；9.计算新的密度矩阵P=2C*C；10.检查是否收敛：是，结束计算；否，返回步骤5。inisncFAAF'一个具体的例子—氨的自洽场计算：输入文件在每次迭代后输出核哈密顿矩阵、MO系数、本征值、密度矩阵、Fock矩阵等一大堆矩阵输出所有的分子轨道信息输出文件先求单电子积分、双电子积分以及重叠积分初猜一组系数，求得密度矩阵由密度矩阵，可求得Fock矩阵简单的单电子乘积是Hartree方程原子体系Slater行列式是Hartree-Fock方程原子和分子体系LCAO-MO方法Hartree-Fock-Roothaan方程（HFR方程）分子体系总结A:•轨道近似忽略了电子相关，所以Hartree方程两种相关都未考虑进去，而HF方程或HFR方程只考虑了自旋相关，并未考虑库仑相关。•在H-F方程中的迭代运算是对分子轨道本身Ψk的迭代，在迭代过程中每次都要计算库仑势和交换势，它们都是很不好算的积分，因此困难很大。而在HFR方程中的迭代是对分子轨道系数cn