第一性原理计算方法ppt课件.ppt
第一性原理计算1基本概念利用自洽场法求解薛定额方程,得到系统的各种性质根据量子力学基本原理最大限度对问题进行非经验处理输入普朗克常数,电子电量,电子质量,光速等基本物理常数 分子团簇、晶体表面、体材料,各种原子、分子体计算与电子结构有关的物理、化学以及力学性能2发展简况1量子力学基础单个粒子时间有关的薛定额方程 ttritrVzyxm),(),()(22222222外加势场不依赖于时间t )()(),(tTrtr)()(222rErVm2222222zyxVmH222EH为本征函数,E为本征值微分本征方程 解释ryydxd)(算符为d/dx本征值方程:算符作用于函数(本征函数),得到本征函数与一个量(本征值)的乘积的结果 一个本征函数为y=ex 本征值r为 1.1算符 量子值,如能量,位置,动量都可以用算符来得到。能量算符哈密顿算符 ddHE*哈密顿算符由势能和动能两部分组成动能算符 222m势能算符 rZeV024沿x轴的动量算符 xi 这个量的期望值ddxipx*1.2 原子单位 1单位电荷电子的绝对电量,e的绝对值1.6021910-19C1质量单位电子的绝对质量,91059310-31kg1单位长度=波尔半径,1单位能量=1Hartree memeae1122201029177. 54JaeEa1800221035981. 441.3 薛定额方程的精确的解 只有一部分的薛定额方程可以精确求解, 箱体中粒子,简谐振子,环中粒子共同特点是必须对可能的解加入限制条件(常称为边界条件). 在无限高势垒中的粒子波函数在边界处必须为0 环中的粒子必须具有2的周期性 波函数的解的特点:1*d正交)(0*nmdnm2 单电子原子 rZemH022242球坐标 原子具有球状结构 波函数可以写为径向函数与角度函数(球谐函数)乘积 )()(lmnlnlmYrRn 主量子数,0,1,2l 角量子数 0,1,2(n-1)m 磁量子数 l,-(l-1), 0. (l-1),l径向函数部分 )()2exp()!(2)!1()2()(1212/1330lnlnlLlnnlnnaZrR0/2naZra0为波尔半径 方括号内为标准化因子 )(121lnLLaguerre多项式 轨道系数,z/nnlRnl(r)1023/2exp(-r)2023/2(1-r)exp(-r)21(4/3)1/25/2rexp(-r)30(2/3)1/23/2(3-6r+22r2)exp(-r)31(8/9)1/25/2(2-r)rexp(-r)32(8/45)1/27/2r2exp(-r)径向分布函数与主量子数的关系 )()(),(mlmlmY)exp(21)(imm)(cos!(!(2) 12()(2/1mllmPmlmll)(cosmlP连带Legendre多项式 轨道的普通图形表示 3 多电子原子和分子 多电子原子和分子的薛定额方程求解复杂化 薛定额方程不能精确求解 波函数可以取多种形式 电子自旋 量子数s 1/2和1/2 自旋角动量z轴的投影+h/2和h/2 电子波函数为依靠于空间坐标的空间函数和依赖于自旋的自旋函数乘积 空间函数描述了电子密度在空间的分布;自旋部分定义了电子的自旋部分,分别为 (1/2)=1, (-1/2)=0,(1/2)=0, (-1/2)=1(电子自旋,可参看量子化学上册261页) 电子是不可区分的 (费米子) 交换一对电子,电子密度的分布保持不变 反对称性 电子互换的时候,波函数改变符号 (波利不相容原理)3.1 Born-Oppenheime近似 原子核的质量远远大于电子的质量 根据原子核的运动,电子可以瞬时进行调整 电子从核子的运动中分离开来 tot=(电子) (核子)32 氦原子 假设:赝原子,两个电子与核相互作用,电子之间不存在相互作用 )2, 1()2, 1(4242202222102212rrErrrZemrZem)2, 1()2, 1(21rrErrHH波函数可以写为两个单电子波函数乘积的形式,)2, 1()2(2) 1( 1)2, 1(rrErrrr)2(2) 1( 1)2(2) 1( 121rrErrHH两边乘以12,对整个空间积分有) 2() 1() 2() 1() 2() 1( 21)2() 1(212121212121rrrrddErrHHrrdd)2() 1()2() 1()2() 1(2)2() 1()2() 1(1)2() 1(212121212121212121rrrrddErrHrrddrrHrrdd波函数是归一化的,那么总的能量E可以写为E1以及E2氦原子中两个电子可能波函数的一般形式电子交换不依赖于电子标签,亦不影响电子密度。假设氦原子的每一个波函数是每个电子解的乘积 低能状态的波函数具有1s轨道的两个电子 1s(1)1s(2)函数满足不可区分原则交换电子的时候,1s(1)1s(2)等于1s(2)1s(1) 第一激发态,一个电子被激发到2s轨道 1s(1)2s(2)1s(2)2s(1)函数不满足不可区分原则线性组合 )1 (2)2(1)2(2) 1 (1 21ssss)1 (2)2(1)2(2) 1 (1 21ssss对称 反对称 1s(1)1s(2)对称 电子的自旋 (1)(2) 对称(1)(2) 对称 (1),(2),(1),(2) )1 ()2()2() 1 (21)1 ()2()2() 1 (21对称反对称多电子体系的总状态波函数一定是反对称的。(反对称原理)这是泊利原理的量子力学表达形式。电子交换的时候必须是反对称 联合一个对称空间函数和反对称的自旋函数 反对称的空间函数以及对称的自旋函数 氦原子基态第一激发态的可取函数形式 )1 ()2()2() 1 ()2(1 ) 1 (121ss)1 ()2()2() 1 ()1 (2)2(1)2(2) 1 (1 21ssss)2() 1 ()1 (2)2(1)2(2) 1 (1 21ssss)2()1()1(2)2(1)2(2)1(121ssss)1 ()2()2() 1 ()1 (2)2(1)2(2) 1 (1 21ssss33 一般的多电子系统和Slater行列式 )()2() 1 (), 3 , 2 , 1 (111NN这种形式的波函数称为Hartree方程 系统的能量等于单个电子自旋轨道能量的和 (H=H1+.+HN)在空间某一特殊点找到一个电子的几率并不依赖于在空间中一点找到其他电子的几率 不符合反对称性原则 电子的运动时关联的 系统具有N个电子,并且具有反对称性假设:粒子之间没有相互作用。 (多个电子之间)低能状态下可以接受的函数 )1 ()2()2() 1 ()2(1 ) 1 (121ss)2()2(1)2()2(1) 1 () 1 (1) 1 () 1 (1ssss两个自旋轨道 ) 1 () 1 (11s)2()2(12s行列式是描述允许的多电子波函数符合反对称性条件的最方便的方法n个电子具有自旋轨道1,n,每一轨道为一空间函数与自旋函数的乘积 )()() 1 () 1 () 1 (!1121NNNNN为Slater行列式交换行列式的任意两行,相当于交换两个电子。改变了行列式的符号,即相当于满足了反对称性的要求如果行列式的两行是相同的,或者说同一轨道上具有两个电子,行列式变为0。这符合了Pauli原则 电子交换奇数次,波函数改变符号;电子交换偶数次,最后仍得到原来的波函数。任意一列加到另一列上,而不改变行列式的值。这意味着自旋轨道并不是唯一的。其他的线性组合也具有相同的能量。 氦原子的低于激发态1s22s2 ) 1 () 1 (2)2()2(1)2()2(2) 1 () 1 (1)2()2(2)2()2(1) 1 () 1 (2) 1 () 1 (1ssssssss221;22121ssss2)2() 1 ()2(2)2(1)1 (2) 1 (1 2)2() 1 ()2(2)2(1)1 (2) 1 (1 )2()2() 1 () 1 (2121ssssssss4 分子轨道计算 41 氢原子:从波函数中计算能量 分子自旋轨道可以表达为原子轨道的线性组合。(LCAO) kiiic1)11 (1BAgssA为了解决分子轨道计算困难,把分子轨道按某个选定的安全基函数集合(基组)展开。这样就可以把对分子轨道的变分转化为对展开系数的变分。Hartree-Fock方程就从一组非线性的积分微分方程转化为一组数目有限的代数方程Hartree-Fock-Roothaan 方程。这组方程仍然是非线性方程,只能用迭代方法求解,但是比微分方程的求解简单了。这是一种近似逼近方法。把在选定的有限基组下满足Hartree-Fock-Roothaan方程的解称为自洽场分子轨道。自洽场分子轨道的极限精确值就是Hartree-Fock轨道。将分子轨道表达为原子轨道线性组合的方法称为LCAO-MO方法。H2低能状态简单的LCAO ) 1 ()2()2() 1 ()2()2() 1 () 1 (11212121) 1 () 1 (1) 1 (1g) 1 () 1 (1) 1 (2g)2()2(1)2(1g)2()2(1)2(2g哈密顿算符:每个电子的动能算符;两个电子与两个原子核之间由于库仑作用两个电子之间的排斥作用 122211222112121rrZrZrZrZHBBAABBAA)2 () 1 () 2 () 1 (12121)2 () 1 () 2 () 1 (2112211222112221122121rrZrZrZrZddEBBAABBAA)2 () 1 () 2 () 1 (121)2 () 1 () 2 () 1 (21122112122121rHHddEBBAArZrZH1121211BBrZH122212)2() 1 (1)2() 1 ()2() 1 (1)2() 1 ()2() 1 ()2() 1 ()2() 1 ()2() 1 ()2() 1 ()2() 1 ()2() 1 ()2() 1 (12122121211221211222121212212112121212112121rddrddHddHddHddHddE对上式中的第一项 ) 1 ()21( ) 1 ()2()2()2() 1 ()2() 1 (11121112222112121BBAArZrZddHdd) 1 () 1 () 1 (1 )21)(1 (1) 1 ()21)(1 (1112111112111drZrZdrZrZdgBBAAgBBAA对第二项进行积分 ) 1 ()21( ) 1 ()2()2()2() 1 ()2() 1 (11121111221212121BBAArZrZddHdd0)2()2(122d上式积分得0电子原子核积分,只有4项非0。每一项等于一个单电子在两个氢原子核场中能量 个剩余的4项为电子与电子的相互作用 )2() 1 (1)2() 1 ()2() 1 (1)2() 1 ()2() 1 (1)2() 1 ()2() 1 (1)2() 1 (21121221121221211212122121122121rddrddrddrdd上式的前两项 )2(1 ) 1 (11)2(1 ) 1 (1) 1 () 1 () 1 () 1 ()2(1 ) 1 (11)2(1 ) 1 (1)2() 1 ()1)(2() 1 (122111122121122121ggggggggrddddrddrdd对应两个轨道间的库仑作用 其他两项为 )2()2() 1 () 1 ()2(1 ) 1 (11)2(1 ) 1 (1)2() 1 ()1)(2() 1 (21122112122121ddrddrddgggg)2()2() 1 () 1 ()2(1 ) 1 (11)2(1 ) 1 (1)2() 1 ()1)(2() 1 (211221211212121ddrddrddgggg为0H2的三个激发态可以通过把电子激发到高能状态得到。这种高能轨道写为1u, 1u=A(1sA+1SB) )2(1)2(1) 1 (1) 1 (1ugug均具有spin() )2()2() 1 () 1 ()2(1 ) 1 (11)2(1 ) 1 (1)2() 1 ()1)(2() 1 (21122112122121ddrddrddugug交叉项并不象基态一样为0 这种作用为交换作用 使相同的自旋电子相互避免。每一个这样的电子好象具有一个洞,这种洞称为交换洞或费米孔 42 多电子系统的能量 N个电子的哈密顿量 )111121(11312112NiBAirrrrH)()2() 1 ()2()2()2() 1 () 1 () 1 (!12121NNNNNNN)3()2() 1 (kji ) 3 () 2 () 1 ()111121() 3 () 2 () 1 (*1131211221kjiNiBAikjiNrrrrdddH ) 3 () 2 () 1 () 3 () 2 () 1 (*21kjikjiNddd在电子在核子的作用场中运动的动能和势能 ) 1 ()21)(1 (121iMAiAAiicoreiirZdHNiNicoreiiiMAiAAiicoretotalHrZdE11121) 1 ()21)(1 (库仑排斥作用 NijjjiiNijjijicoulombirddrddE) 2 () 2 (1) 1 () 1 () 2 () 1 (1) 2 () 1 (12211221NiNijijNijjjiiNicoulombiJrddE11112211) 2() 2(1) 1 () 1 (交换作用 具有相同自旋的电子倾向于相互避免,它们承受了一个低的库仑排斥作用,给出了一个低的能量 ) 1 ()2(1)2() 1 (1221jijiijrddk) 1 ()2(1)2() 1 (1221jijiNijexchangeirddE) 1 ()2(1)2() 1 (122111jijiNijNiexchangetotalrddE库仑作用简写jijir12*1ijrij121交换作用简写jirij121(iijj) (ijji)