(完整版)2-第一性原理与密度泛函理论课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《(完整版)2-第一性原理与密度泛函理论课件.ppt》由用户(三亚风情)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 完整版 第一性 原理 密度 理论 课件
- 资源描述:
-
1、材料设计与计算机模拟第一性原理与密度泛函理论王 锐Material designing&Computer simulation引 言l 量力力学是反映微观粒子(分子、原子、原子核、基本粒子)运动规律的理论。l 以量子力学为基础,结合高速发展的计算技术分别建立起来的计算材料科学、计算物理、量子化学等分枝学科,促进了物理学、化学和材料科学的发展,为发展和设计新型材料提供了理论基础和新的研究方法。l 在理论上最具诱惑力,且在将来最有可能开展真正意义上的材料设计的计算就是解体系的Schrdinger方程,即为计算材料学中的第一原理计算。多粒子体系的第一原理 第一原理的基本思想:将多原子构成的体系理解为
2、由电子和原子核组成的多粒子系统,在解体系Schrdinger方程的过程中,最大限度地进行“非经验性”处理,即不涉及任何经验参数,所要输入的只是原子的核电荷数和一些模拟环境参量。计算所求得的结果是体系Schrdinger方程的本征值和本征函数(波函数),有了这两项结果,就可研究体系的基本物理性质。多粒子体系的第一原理 从微观角度看,一块物质材料是由大量(每立方厘米约1023个)原子核和游离于原子核之间的电子组成。因此,材料的性质(如硬度、电磁和光学性质)和发生在固体内的物理和化学过程是由它所包含的原子核及其电子的行为决定的。20世纪初量子力学的出现,原则上提出了像原子核和电子这样的微观粒子运动和
3、交互作用的定律。理论上,给定一块固体化学成分(即所含原子核的电荷和质量),我们就可以计算这些固体的性质。因为一块固体实际上是一个多粒子体系。决定这个体系性 质 的 波 函 数 可 以 通 过 解 薛 定 谔(Schrdinger)波动方程来获得。多粒子体系的第一原理 多粒子系统的Schrdinger方程其中 和H分别对应于多粒子系统的波函数和哈密顿量。原则上只要对上式进行求解即可得出所有物理性质,然而由于电子之间的相互作用的复杂性,要严格求出多电子体系的Schrdinger方程解是不可能的,必须在物理模型上进一步作一系列的近似。多粒子体系的第一原理换句话说,做第一原理计算(ab initio
4、calculation)便可知道一块固体的性质。可是,这个薛定谔波动方程有3N个变量(N是粒子总数),极其复杂,假使我们把目前世界上的所有电脑都用上,让它跑千年、万年都不可能算出来。正如1929年量子物理大师狄拉克(Dirac)所言:处理大部分物理学和全部化学问题的基本定理已经完全知道。困难在於这些定律的应用所引出的数学方程(Schrdinger方程)太复杂以致於无法解决。多粒子体系的第一原理l 价电子近似 在原子模型中,近核电子可被紧紧束缚在核的周围,这些束缚电子是定域的,比较稳定,因而对固体性质的贡献很小。而在外层轨道的价电子可以是离域的,当原子结合在一起组成固体时,这些电子的状态变化很大
5、,对固体的电学和光学性质有决定性的影响。因此,可将固体看作是由原子核束缚电子构成的离子核和价电子组合而成。多粒子体系的第一原理l 绝热近似 由于原子核的质量远大于电子的质量,原子核的运动速度要比电子慢很多,因此可以认为电子运动在固定不动的原子核的势场中,所以原子核的动能为零,而势能为一个常数。中子/质子的质量是电子质量的约1835倍,即电子的运动速率比核的运动速率要高3个数量级,因此可以实现电子运动方程和核运动方程的近似脱耦。这样,电子可以看作是在一组准静态原子核的平均势场下运动。为了对问题进行简化,可以把多电子系统中的相互作用视为有效场下的无关联的单电子的运动。一个具有N个电子的系统的总波函
6、数可以写成所有单电子的波函数的乘积。Hartree方程Hartree-Fock近似Hartree-Fock近似Hartree-Fock近似单粒子Hartree-Fock算符是自恰的,即它决定于所有其他单粒子HF方程的解,必须通过叠代计算来求解,具体操作如下:猜测试探波函数 构造所有算符 求解单粒子赝薛定谔方程 对于解出的新的波函数,重新构造Hartree-Fock算符 重复以上循环,直到收敛(即前后叠代的结果相同)自恰场(SCF)方法是求解材料电子结构问题的常用方法对处理原子数较少的系统来说,Hartree-Fock近似是一种很方便的近似方法。但用于原子数大的系统,问题就变得非常复杂,此计算方
7、法的计算量随着电子数的增多呈指数增加,这种计算对计算机的内存大小和CPU的运算速度有着非常苛刻的要求,它使得对具有较多电子数的计算变得不可能。同时Hartree-Fock近似方法给出的一些金属费米能和半导体能带的计算结果和实验结果偏差较大。这在很大程度上导致了密度泛函理论的产生。Hartree-Fock近似 Hartree-Fock方法的主要缺陷 完全忽略电子关联作用 计算量偏大,随系统尺度4次方关系增长 Density Functional Theory(DFT 1964)一种用电子密度分布n(r)作为基本变量,研究多粒子体系基态性质的新理论密度泛函理论W.KohnW.Kohn 荣获1998
8、年Nobel 化学奖密度泛函理论 1927年,Thomas-Fermi理论首先运用 作为计算原子中电子结构的基本变量,从而给出一个简化处理方法。如果能够证明,对于任意电子系统的基态电子密度的分布唯一地决定电子系统的情况,则基于的一个计算电子结构的新的理论框架就可以建立起来,这就是表述密度泛函理论(density functional theory,简称DFT)的出发点。密度泛函理论在密度泛函理论中,将电子密度作为描述体系状态的基本变量,可追溯到Thomas和Fermi用简并的非均匀电子气来描述单个原子的多电子结构。直到Hohenberg和Kohn提出了两个基本定理才奠定了密度泛函理论的基石。随
9、后Kohn和Sham的工作使密度泛函理论成为实际可行的理论方法。密度泛函理论1964年年Hohenberg和和Kohn推证(推证(P.Hoheuberg and W.Kohn,Physical Review 136B,864(1964):一个多粒子(如电子)体系性质(如体:一个多粒子(如电子)体系性质(如体系的能量)由粒子密度的空间分布(即空间任意一点的粒子平均个数)系的能量)由粒子密度的空间分布(即空间任意一点的粒子平均个数)完全决定。换句话说:完全决定。换句话说:我们无需去求解描述每个粒子运动的体系波函我们无需去求解描述每个粒子运动的体系波函数,我们只要设法找出仅有数,我们只要设法找出仅有
10、3个变量的粒子密度的空间函数就行了。个变量的粒子密度的空间函数就行了。1998获奖诺贝尔化学奖获奖诺贝尔化学奖表彰表彰Walter Kohn在在60年代提出年代提出密度泛函理论密度泛函理论及及John A.Pople发明了测验化学结构和物质特性的发明了测验化学结构和物质特性的计算机技术计算机技术密度泛函理论1965年柯恩又和沈吕九证明(W.Kohn and L.J.Shan,Physical Review 140,All33):一个多粒子体系的粒子密度函数可以通过一个简单的单粒子波动方程获得。这个单粒子波动方程现在被称作柯恩沈(Kohn-Sham)方程。Hohenberg,Kohn和Shan的
11、理论就是诺贝尔化学奖颁词所指的密度泛函理论,显然,密度泛函理论大大简化了应用量子力学探讨材料物理性质所涉及的数学问题。W.KohnLu.J.Sham(沈吕九)密度泛函理论密度泛函理论物质电子结构的新理论1.氢原子氢原子1)Bohr:电子粒子电子粒子2)Schrodinger:电子波电子波(r).3)DFT:电子是电子云的密度分布。电子是电子云的密度分布。n(r).密度泛函理论3)DFT:电子是电子云的密度分布。电子是电子云的密度分布。2.DFT中的氢分子。中的氢分子。由密度分布表示。由密度分布表示。密度泛函理论 3.大分子(例如DNA)N个原子Schrodinger:(r1,r2,r3,rN)
12、,3N维空间DFT:n(r)3维空间 Hohenberg-Kohn定理I1964年,P.Hohenberg和W.Kohn在非均匀电子气理论的基础上,提出两个基本定理,奠定了密度泛函理论的基础。l 定理1:对于一个共同的外部势v(r),相互作用的多粒子系统的所有基态性质都由(非简并)基态的电子密度分布n(r)唯一地决定。或:对于非简併基态,粒子密度分布n(r)是系统的基本变量。密度泛函理论 考虑一个多粒子系(电子体系、粒子数任意),在外部势和相互作用Coulomb势作用下,Hamiltonian为密度泛函理论r rHTVUTrr drVv rrr drUrrrr drdr 12112()()()
13、()()()()()()电子密度算符电子密度算符电子密度分布电子密度分布n(r)是是 的期待值:的期待值:)()()(rrrn)(rn)(,()(rnrn Hohenberg-Kohn定理IIl 定理2:如果n(r)是体系的密度分布,则En(r)是最低的能量,即体系的基态能量密度泛函理论实际计算是利用能量变分原理,使系统能量达到最低(有一定精度要求)。由此求出体系的真正电荷密度n(r),进而计算体系的所有其它基态性质。如能带结构,晶格参数,体模量等等。密度泛函理论 Hohenberg-Kohn定理指出,多粒子系统的所有基态性质都是密度泛函的唯一泛函。因此系统的总能量 可表达为其中 Tn代表无相
14、互作用的粒子的动能,Un代表一般的库仑能,Excn 代表多体效应的交换关联能。假设N个电子的系统有N个单电子波函数,则电荷密度为:密度泛函理论密度泛函理论密度泛函理论Kohn-Sham近似的核心思想:近似的核心思想:1.动能的大部分通过相同电子密度的无相互作用体系来计算动能的大部分通过相同电子密度的无相互作用体系来计算2.电子相互作用中库仑作用占据了主要部分,而交换相关是电子相互作用中库仑作用占据了主要部分,而交换相关是相对次要的相对次要的3.非经典的交换和相关作用,动能校正项,自相互作用折入非经典的交换和相关作用,动能校正项,自相互作用折入交换相关泛函中交换相关泛函中解解Kohn-ShamK
15、ohn-Sham方程的流程图方程的流程图nin(r)n(r)=nat(r)求解、Vxc、Veff求解Kohn-Sham方程得到i由i构造nout(r)比较nin与 nout(r)计算总能EtotNoYesnin与nout混合原子计算精度控制NoYes输出结果:Etot、i、n(r)Vxc、Veff、En(k)、N(E)30虽然K-S方程十分简单,其计算量也只有Hartree方程的水平,但却包含着深刻得多的物理内容。其中一个重要的概念性结果是,多体基态的解被准确地简化为基态密度分布之解,而这个密度是由单粒子的Schrdinger方程给出的。由此,方程中的有效势在原理上包括了所有的相互作用效应,即
展开阅读全文