第一性原理计算课件.ppt

《第一性原理计算课件.ppt》由会员分享，可在线阅读，更多相关《第一性原理计算课件.ppt（49页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、第一性原理计算第1页，此课件共49页哦 量子化学从头计算方法仅仅利用 Planck 常量、电子质量、电量三个基本物理常数以及元素的原子序数，不再借助于任何经验参数，计算体系全部电子的分子积分，求解 Schrodinger 方程。一般所说的量子化学从头计算是建立在三个基本近似基础上的计算方法，对于三个基本近似进行各种校正的计算方法是高级从头计算。从头计算法中的三个基本近似为：第2页，此课件共49页哦a）非相对论近似 认为电子质量等于其静止质量，即 me（电子质量）=me，0（静止质），并认为光速接近无穷大。第3页，此课件共49页哦b）BornOppenheimer 近似（也叫绝热近似）即将核运动

2、和电子运动分离开来处理。由于原子核质量一般比电子的质量约大 103105 倍，分子中核的运动比电子的运动要慢近千倍。因此在电子运动时，可以把核近似看作不动。第4页，此课件共49页哦c c）轨道近似（又叫单电子近似，由 Hartree 提出）轨道一词是从经典力学中借用来的概念，在量子化学中指单电子波函数，原子的单电子波函数称为原子轨道，分子的单电子波函数称为分子轨道。轨道近似是把 N个电子体系的总波函数写成 N 个单电子函数的乘积：(1 1,2 2,.,.,N N)=)=1 1(1 1)2 2(2 2).).N N (N N)其中每一个电子波函数i(i)只与一个电子坐标 i有关。第5页，此课件共

3、49页哦 从 头 计 算 方 法 基 于 上 述 三 个 基 本 近 似 后，不 再 作 任 何 近 似 来 求 解Schrodinger 方程，这三个基本近似即为从头计算方法的头，而 Schrodinger方程在引入三个近似后的具体表达形式为 HartreeFockRoothaan 方程。量子化学从头计算对于分子基态性质的研究一般是可靠的，但对于晶体材料来说计算量过大，必须采用进一步的近似。第6页，此课件共49页哦 针对实际材料体系和所研究的问题不依赖实验数据及经验公式，完全由最基本的物理定律出发。对材料科学研究而言，第一原理，是量子力学。用完备的形式理论体系统一理论解释和预测不同科学领域的

4、实验结果。量子力学的“第一原理”(First Principle)计算（从头算）只采用几个基本物理常数而不依赖任何经验参数即可正确预测微观体系的状态和性质。第7页，此课件共49页哦1量子力学基础单个粒子时间有关的薛定额方程 外加势场不依赖于时间t 第8页，此课件共49页哦为本征函数，E为本征值微分本征方程 第9页，此课件共49页哦解释算符为d/dx本征值方程：算符作用于函数（本征函数），得到本征函数与一个量（本征值）的乘积的结果 一个本征函数为y=ex 本征值r为 第10页，此课件共49页哦1.1算符 量子值，如能量，位置，动量都可以用算符来得到。能量算符哈密顿算符 哈密顿算符由势能和动能两部

5、分组成动能算符 势能算符 第11页，此课件共49页哦沿x轴的动量算符 这个量的期望值第12页，此课件共49页哦1.2 原子单位 1单位电荷电子的绝对电量，e的绝对值1.6021910-19C1质量单位电子的绝对质量，91059310-31kg1单位长度=波尔半径,1单位能量=1Hartree 第13页，此课件共49页哦1.3 薛定额方程的精确的解 只有一部分的薛定额方程可以精确求解,箱体中粒子,简谐振子,环中粒子共同特点是必须对可能的解加入限制条件(常称为边界条件).在无限高势垒中的粒子波函数在边界处必须为0 环中的粒子必须具有2的周期性 波函数的解的特点：正交第14页，此课件共49页哦2 单

6、电子原子 球坐标 第15页，此课件共49页哦原子具有球状结构 波函数可以写为径向函数与角度函数（球谐函数）乘积 n 主量子数，0,1,2l 角量子数 0，1，2（n-1）m 磁量子数 l,-(l-1),0.(l-1),l第16页，此课件共49页哦径向函数部分 a0为波尔半径 方括号内为标准化因子 Laguerre多项式 轨道系数，z/n第17页，此课件共49页哦nlRnl(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)3

7、2(8/45)1/27/2r2exp(-r)第18页，此课件共49页哦径向分布函数与主量子数的关系 第19页，此课件共49页哦连带Legendre多项式 第20页，此课件共49页哦轨道的普通图形表示 第21页，此课件共49页哦3 多电子原子和分子 多电子原子和分子的薛定额方程求解复杂化 薛定额方程不能精确求解 波函数可以取多种形式 电子自旋 量子数s 1/2和1/2 自旋角动量z轴的投影+h/2和h/2 电子波函数为依靠于空间坐标的空间函数和依赖于自旋的自旋函数乘积 空间函数描述了电子密度在空间的分布；自旋部分定义了电子的自旋部分，分别为（1/2）=1,（-1/2）=0,(1/2)=0,(-1

8、/2)=1 第22页，此课件共49页哦电子是不可区分的 交换一对电子，电子密度的分布保持不变 反对称性 电子互换的时候，波函数改变符号 第23页，此课件共49页哦3.1 Born-Oppenheime近似 原子核的质量远远大于电子的质量 根据原子核的运动，电子可以瞬时进行调整 电子从核子的运动中分离开来 tot=(电子)（核子)第24页，此课件共49页哦32 氦原子 假设：赝原子，两个电子与核相互作用，电子之间不存在相互作用 波函数可以写为两个单电子波函数乘积的形式，第25页，此课件共49页哦两边乘以12，对整个空间积分有波函数是归一化的，那么总的能量E可以写为E1以及E2第26页，此课件共4

9、9页哦氦原子中两个电子可能波函数的一般形式电子交换不依赖于电子标签，亦不影响电子密度。假设氦原子的每一个波函数是每个电子解的乘积 低能状态的波函数具有1s轨道的两个电子 1s(1)1s(2)函数满足不可区分原则交换电子的时候，1s(1)1s(2)等于1s(2)1s(1)第一激发态，一个电子被激发到2s轨道 1s(1)2s(2)1s(2)2s(1)函数不满足不可区分原则第27页，此课件共49页哦线性组合 对称 反对称 1s(1)1s(2)对称 第28页，此课件共49页哦电子的自旋（1）（2）对称（1）（2）对称（1）,（2）,（1）,（2）对称反对称第29页，此课件共49页哦电子交换的时候必须是

10、反对称 联合一个对称空间函数和反对称的自旋函数 反对称的空间函数以及对称的自旋函数 氦原子基态第一激发态的可取函数形式 第30页，此课件共49页哦33 一般的多电子系统和Slater行列式 这种形式的波函数称为Hartree方程 系统的能量等于单个电子自旋轨道能量的和 在空间某一特殊点找到一个电子的几率并不依赖于在空间中一点找到其他电子的几率 不符合反对称性原则 电子的运动时关联的 系统具有N个电子，并且具有反对称性第31页，此课件共49页哦低能状态下可以接受的函数 两个自旋轨道 第32页，此课件共49页哦行列式是描述允许的多电子波函数符合反对称性条件的最方便的方法n个电子具有自旋轨道1,n,

11、每一轨道为一空间函数与自旋函数的乘积 第33页，此课件共49页哦为Slater行列式交换行列式的任意两行，相当于交换两个电子。改变了行列式的符号，即相当于满足了反对称性的要求如果行列式的两行是相同的，或者说同一轨道上具有两个电子，行列式变为0。这符合了Pauli原则 电子交换奇数次，波函数改变符号；电子交换偶数次，最后仍得到原来的波函数。任意一列加到另一列上，而不改变行列式的值。这意味着自旋轨道并不是唯一的。其他的线性组合也具有相同的能量。第34页，此课件共49页哦氦原子的第一激发态1s22s2 第35页，此课件共49页哦4 分子轨道计算 41 氢原子：从波函数中计算能量 分子自旋轨道可以表达

12、为原子轨道的线性组合。（LCAO）H2低能状态简单的LCAO 第36页，此课件共49页哦哈密顿算符：每个电子的动能算符；两个电子与两个原子核之间由于库仑作用两个电子之间的排斥作用 第37页，此课件共49页哦第38页，此课件共49页哦对上式中的第一项 第39页，此课件共49页哦对第二项进行积分 上式积分得0电子原子核积分，只有4项非0。每一项等于一个单电子在两个氢原子核场中能量 第40页，此课件共49页哦剩余的4项为电子与电子的相互作用 上式的前两项 对应两个轨道间的库仑作用 第41页，此课件共49页哦其他两项为 为0第42页，此课件共49页哦H2的三个激发态可以通过把电子激发到高能状态得到。这

13、种高能轨道写为1u，1u=A(1sA-1SB)均具有spin()交叉项并不象基态一样为0 第43页，此课件共49页哦这种作用为交换作用 使相同的自旋电子相互避免。每一个这样的电子好象具有一个洞，这种洞称为交换洞或费米孔 第44页，此课件共49页哦42 多电子系统的能量 N个电子的哈密顿量 第45页，此课件共49页哦第46页，此课件共49页哦在电子在核子的作用场中运动的动能和势能 库仑排斥作用 第47页，此课件共49页哦交换作用 具有相同自旋的电子倾向于相互避免，它们承受了一个低的库仑排斥作用，给出了一个低的能量 第48页，此课件共49页哦库仑作用简写交换作用简写(iijj)(ijji)第49页，此课件共49页哦