第一性原理计算方法

第一性原理计算1基本概念利用自洽场法求解薛定额方程，得到系统的各种性质根据量子力学基本原理最大限度对问题进行非经验处理输入普朗克常数，电子电量，电子质量，光速等基本物理常数分子团簇、晶体表面、体材料，各种原子、分子体计算与电子结构有关的物理、化学以及力学性能2发展简况1量子力学基础单个粒子时间有关的薛定额方程ttritrVzyxm),(),(})(2{2222222外加势场不依赖于时间t)()(),(tTrtr)()(222rErVm2222222zyxVmH222EH为本征函数，E为本征值微分本征方程解释ryydxd)(算符为d/dx本征值方程：算符作用于函数（本征函数），得到本征函数与一个量（本征值）的乘积的结果一个本征函数为y=ex本征值r为1.1算符量子值，如能量，位置，动量都可以用算符来得到。能量算符－哈密顿算符ddHE**•哈密顿算符由势能和动能两部分组成动能算符222m势能算符rZeV024沿x轴的动量算符xi这个量的期望值ddxipx**1.2原子单位1单位电荷＝电子的绝对电量，e的绝对值＝1.6021910-19C1质量单位＝电子的绝对质量，91059310-31kg1单位长度=波尔半径,1单位能量=1Hartreememeae1122201029177.54JaeEa1800221035981.441.3薛定额方程的精确的解只有一部分的薛定额方程可以精确求解,箱体中粒子,简谐振子,环中粒子共同特点是必须对可能的解加入限制条件(常称为边界条件).在无限高势垒中的粒子波函数在边界处必须为0环中的粒子必须具有2的周期性波函数的解的特点：1*d正交)(0*nmdnm2单电子原子rZemH022242球坐标原子具有球状结构波函数可以写为径向函数与角度函数（球谐函数）乘积)()(lmnlnlmYrRn主量子数，0,1,2l角量子数0，1，2（n-1）m磁量子数－l,-(l-1),0.(l-1),l径向函数部分)()2exp(])![(2)!1()2()(1212/1330lnlnlLlnnlnnaZrR0/2naZr•a0为波尔半径•方括号内为标准化因子)(121lnLLaguerre多项式轨道系数，＝z/nnlRnl(r)1023/2exp(-r)2023/2(1-r)exp(-r)21(4/3)1/25/2rexp(-r)30(2/3)1/23/2(3-6r+22r2)exp(-r)31(8/9)1/25/2(2-r)rexp(-r)32(8/45)1/27/2r2exp(-r)径向分布函数与主量子数的关系)()(),(mlmlmY)exp(21)(imm)(cos!(!(2)12()(2/1mllmPmlmll)(cosmlP连带Legendre多项式轨道的普通图形表示3多电子原子和分子多电子原子和分子的薛定额方程求解复杂化薛定额方程不能精确求解波函数可以取多种形式电子自旋量子数s1/2和－1/2自旋角动量z轴的投影+h/2和－h/2电子波函数为依靠于空间坐标的空间函数和依赖于自旋的自旋函数乘积空间函数描述了电子密度在空间的分布；自旋部分定义了电子的自旋部分，分别为（1/2）=1,（-1/2）=0,(1/2)=0,(-1/2)=1（电子自旋，可参看量子化学上册261页）电子是不可区分的（费米子）交换一对电子，电子密度的分布保持不变反对称性电子互换的时候，波函数改变符号（波利不相容原理）3.1Born-Oppenheime近似原子核的质量远远大于电子的质量根据原子核的运动，电子可以瞬时进行调整电子从核子的运动中分离开来tot=(电子)（核子）3．2氦原子假设：赝原子，两个电子与核相互作用，电子之间不存在相互作用)2,1()2,1(4242202222102212rrErrrZemrZem)2,1()2,1(}21{rrErrHH波函数可以写为两个单电子波函数乘积的形式，)2,1()2(2)1(1)2,1(rrErrrr)2(2)1(1)2(2)1(121rrErrHH两边乘以12，对整个空间积分有)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(121ss)]1()2()2()1()][1(2)2(1)2(2)1(1[21ssss)2()1()]1(2)2(1)2(2)1(1[21ssss)2()1()]1(2)2(1)2(2)1(1[21ssss)]1()2()2()1()][1(2)2(1)2(2)1(1[21ssss3．3一般的多电子系统和Slater行列式)()2()1(),3,2,1(111NN这种形式的波函数称为Hartree方程•系统的能量等于单个电子自旋轨道能量的和（H=H1+….+HN）•在空间某一特殊点找到一个电子的几率并不依赖于在空间中一点找到其他电子的几率•不符合反对称性原则•电子的运动时关联的•系统具有N个电子，并且具有反对称性假设：粒子之间没有相互作用。（多个电子之间）低能状态下可以接受的函数)]1()2()2()1()[2(1)1(121ss)2()2(1)2()2(1)1()1(1)1()1(1ssss两个自旋轨道)1()1(11s)2()2(12s行列式是描述允许的多电子波函数符合反对称性条件的最方便的方法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(1ssssssss221;221'2'1ssss2)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('2'1'2'1ssssssss4分子轨道计算4．1氢原子：从波函数中计算能量分子自旋轨道可以表达为原子轨道的线性组合。（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(1g)1()1(1)1(2g)2()2(1)2(1g)2()2(1)2(2g哈密顿算符：•每个电子的动能算符；•两个电子与两个原子核之间由于库仑作用•两个电子之间的排斥作用122211222112121rrZrZrZrZHBBAABBAA)]2()1()2()1([]12121[)]2()1()2()1([2112211222112221122121rrZrZrZrZddEBBAABBAA)]2()1()2()1(][121[)]2()1()2()1([21122112122121rHHddEBBAArZrZH1121211BBrZH122212)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(12122121211221211222121212212112121212112121rddrddHddHddHddHddE对上式中的第一项)1()21()1()2()2()2()1()2()1(11121112222112121BBAArZrZddHdd)1()1()1(1)21)(1(1)1()21)(1(1112111112111drZrZdrZrZdgBBAAgBBAA对第二项进行积分)1()21()1()2()2()2()1()2()1(11121111221212121BBAArZrZddHdd0)2()2(122d•上式积分得0•电子－原子核积分，只有4项非0。•每一项等于一个单电子在两个氢原子核场中能量个剩余的4项为电子与电子的相互作用)2()1(1)2()1()2()1(1)2()1()2()1(1)2()1()2()1(1)2()1(21121221121221211212122121122121