![点击分享此内容可以赚币 分享](/master/images/share_but.png)
密度泛函理论与从头计算分子动力学.pdf
《密度泛函理论与从头计算分子动力学.pdf》由会员分享,可在线阅读,更多相关《密度泛函理论与从头计算分子动力学.pdf(49页珍藏版)》请在咨信网上搜索。
1、密度泛函理论与从头计算分子动力学$1引言自从上世纪60年代以来,密度泛函理论(DFT)建立.并在局域密度近似(LDA)下导出著名 的Ko hn-Sham(KS)方程以来,DFT一直是凝聚态物理领域计算电子结构及其特性的最有力 的工具.近几年来DFT与分子动力学相结合,在材料设计,合成,模拟计算和评价诸多方面有 明显进展,成为计算材料科学的重要基础和核心技术.特别在量子化学计算领域1987年以 前主要用Har t r ee-Fo c k(HF)方法.但近年来用DFT的工作以指数增加.以致于HF方法应用已 经相当少.W.Ko hn因提出DFT获得1998年诺贝尔化学奖.已经表明了DFT在计算化学领
2、域 的核心作用与应用的广泛性.分子动力学计算机模拟是研究复杂的凝聚态系统的有力工具.这一技术既能得到原 子的运动轨迹,还能像做实验一样作各种观察.对于平衡系统何以在一个分子动力学观 察(o bs er v at io n t ime)内作时间平均来计算一个物理量的统计平均值.对一个非平衡系统过 程,只要发生在一个分子动力学观察时间内(Ll Ops)的物理现象也可以用分子动力学计算进 行直接模拟.可见数值实验对理论与实验的有力补充,特别是许多与原子有关的微观细节,在 实验中无法获得而在计算机模拟中可以方便的得到.DFT与分子动力学(MD)相结合可以有大量不同类型的应用.如品格生长孙延生长,离 子
3、移植,缺陷运动,无序结构,表而与界面的重构,电离势的计算,振动谱的研究化学反应的 问题,生物分子的结构,催化活性位置的特性以及材料电子结构和几何结构.固,液体的相变 等.现在这些方法已发展成为成熟的计算方法.DFT的另一个特点是,它提供了第一性原 理(fir s t-pr inc ipal)或称为从头计算的理论框架在这个框架下可以发展各种各样的能带计 算方法.虽然在DFT的所有这些实际应用中,几乎都采用局域密度近似(LDA),这是一种不能 控制精度的近似,因而DFT方法的有效性在很大程度上要看结果与实验的一致性.人们没 有任何直接的方法可以改善LDA的精度.然而DFT允许发展别的方法加以补充.
4、例如广义 梯度近似(GGA)等方法,把密度分布n)的空间变化包括在方法之中,实现了比较大幅度减 少LDA误差的目的.相比较传统的量子化学方法,如组态相互作用(CI)方法QFT+MD方法显然可以用于数 百个原子的大分子体系但对于具有强关联相互作用的体系,似乎化学家更愿意应用CI方 法.而对于凝聚态物理领域QFT+MD方法可以相当精确的计算材料的电子结构和相应的 许多物性.在DFT获得巨大成功的同时,也有些不可忽视的弱点与困难.针对这些问题已经发展了许2Self-consistent All-electron full potential Harris functional-Ail-electi-
5、cn mufnn-tin-Pseudopotential-JelliumGW.-Time-dependent DFT Beyond LDA-Optimized effective potentials(OEP)LDA+U,LDA+一 Self-interacticn correcions(SIC)Fully-relativistic-Semi-relativistic-Non-relativistic-Local density approximation(LDA)4v2+K(r)+(r)iirf 二E:必Non-periodic-Periodic-Symmetry-Locality,O(N)
6、一 Non-spin-polarized Spin-polarized-Atomic orbitals Slater type(STO)I Gaussians(GTO)*-Numerical(Dmol,TB)-1-Plane waves(FLAPW)Spherical waves(LMTO,ASW)-Fully numerical-O(N)一吨Plane waves-Augmentation图1密度泛函理论的发展多不同的方法,这些方法可以用Ko hn-Sham方程的有效Hamil t o nian的各部分和波函数构造 上的考虑进行归类,如图一所示传统的DFT LDA是预言多电子体系基态性质的理
7、论。对于激发态性质的描述总是与 实验不符合,这固然因为(1)激发态本身存在着复杂性质。(2)DFT LDA理论存在 着对激发态描述困难。此外,对于具有强关联相互作用体系LDA理论描述不够好。因而 发展出LDA+U和LDA+等方法。此外关于大原子数的复杂体系,近几年来发展了各种 线性标度的方法,也称为0(N)算法。为研究复杂体系提供了有力的工具。$2基本理论首先,对包括诸多电子,离子,原子的凝聚态体系来说,面临的首要问题是,由于原 子核较重而且运动中无法跟随电子的运动。因而可以把电子与原子分开考虑,原子实作为 3带有正电荷的外场存在。这就是所谓的绝热近似或称B0近似(Bo r n-Opp6nhe
8、imer).因而 我们可以把电子的哈密顿量写成H=T+V+Vext(2.1)$2.1 Hohenberg-Kohn定理1964年,HK提出HK两个定理:1)多电子系统在外场势匕J乍用下,其基态的电子密度p(f%与基态下任意力学量。的 可观测量是一一对应的,且是严格的基态电子密度p(竹的泛函.=0p(2.2)2)若。为哈密顿量方,则基态的总能量函数Hp=Ev p具有如下形式EveXtp=+(2.3)而对任意多电子体系的普适量尸HK,为FHKp=EVext p=Fhk p+fp vext(r)dr(2.4)这里,我们不准备证明上述两个定理,只说明如下3点1)任何p与力学量或势之间存在一一对应性(o
9、 ne t o o ne c o r r es po ndenc e)2)普适函数尸可以很容易用密度算符表示。3)HK第二定理可以使我们利用变分定理,通过求能量极小值来获得基态的电子密度。$2.2 Kohn-Sham方程K S方程诞生于1965年,K S方程提出后,才真正使DFT成为实际计算可能。设总能 量函数乙同(包括电子间相关能)及在Har t r ee-Fo c k近似下的总能量函数为石印加(不包 括电子间相关能)贝(J Ee=T+V(2.5)Ehf=Tq+(VH+Vx)(2.6)VT,V为严格的动能及电子间的势能函数。毁为无相互作用的电子气动能,质为电子间 的库仑能,口为电子间交换能。
10、因而,不难得到上述两种体系下能量差别仅仅在于电子 相关能.则Vc=Ee-Ehf=T-Tg(2.7)则HK的普适函数外长可以写为:Fhk T+V+TqTq Tq+V+(T 7q)Tq+V-Vc-Vh Vh Tq+Vh+Vc+(V Vn)Vc Vx4=n+0+(“+%)(2.8)-一 Vxc 这里“被称为交换-关联势函数,这样总能量泛函可以写为EVextp=T0p0+(1/2)/而加。言索+/Vextpdr+Exc p。(2这里交换关联势可由下式给出SEXcp yxc=7-而对应的系统哈密顿量(通常称为KS哈密顿量)写为ks6i=-|v +=E(2.10)这里单粒子波函数满足N,(了)=。:&(2
11、.11)2=1Veff=Vextf+/*产2T+笔F(2.12)J f-rf M(/)(2.9)-(2.12)式就构成了KS方法的基本方程组。遗憾的是,它还是形式上的东西,因为包含电子多体效应交换关联效应的Exc的具体 表达式还不清楚,需作某种近似.现在广泛应用的是局域密度近似(LDA),即把非均匀 电子系统分割成一些小块。在这些小块中认为电子分布是均匀的。这样F处的子块中的交 换关联能密度”只取决于该点的p/),整个系统的交换关联能为Exc=/历)5p(2.13)至于交换关联势Exc的具体表达式不断有著者给出。目前用的最广泛的 是Cepel ey和Al der在八十年代用Mo nt e Ca
12、r l o方法导出的解析表(P hys.R ev.l et t(566-568)Vo l.45.No.7)o从方法推导近似来看,LDA只局限于那些电子密度缓变的系统。因 而,在1996年,P er dew,Bur kgEr nzer ho f把它推广到包括电子梯度密度在内的函数,并给 出了解析表示式。称为GGA方法。(Gener al ized Gr adient Appr o ximat io n.P hys.R ev.l et t Vo l.77.No.8(3895-3897),如包含电子自旋极化(LSD)情况在内,贝(J(2.13)可写为E射n。八=/d3r n嘤订n=八(2.14)而在
13、GGA下,则为E肥 n?,ni=I d3r f(n n Vn Vn(2.15)5现在GGA方法已经得到广泛的应用。总之,局域密度泛函理论,自上世纪60年代提出以 来,经过许多人的发展与完善,已广泛的应用于各领域多原子体系的电子结构计算。由于 在此理论近似下,单电子运动方程中的交换关联势形式简单,使计算工作量大大减少。因 此,现在大多数能带或者原子集团模型中的计算方法,都是在局域密度泛函理论近似的基 础上建立起来的。$2.3 Kohn-Sham方程的求解在传统的固体,分子的局域密度泛函电子结构计算中,KS自洽方程组一般可以这样求 解:先给定一个电荷p/竹,比如说可以选用各原子电荷密度的迭加,然后
14、代入(2.12)计算有 效势匕万,再通过对角化来求解(2.10),得到方ks的本征矢。即&。然后可以求得新的分 布电荷密度外(可,继续进行迭代直致自洽。由于对角化计算工作量以河3增加(M为KS轨 道展开所需基函数的个数),因而对有数百个原子的系统工作量非常庞大。KS方程组也可以通过其它方法解。为了求得能量泛函对电子密度变分的最小值,我们可以在电子自由度出的空间中引入一个虚拟的动力学过程,来使泛函达到极小值。最简单 动力学考虑是速降法。(2.16)上式中的点代表对时间的微分。变分患=HksM依赖于仇,求解上式时,为保证波 函数的正交性,还得加上一个约束条件,虚拟动力学方程是:6”二+工 血 17
15、)实际上,先用无约束的动力学方程(2.16)计算所,然后通过适当的正交化过程对仇进行 修正,并得到系数入户当系统达到极小值时,加=0,即为5弧=2他与把相应的约束 矩他对角化后,既可以得到KS方程的本征矢与本征值。经验表明,对于一定的体系坐 标E只有一个极小值。这样速降法也能得到能量极小值。这种方法与传统的方法相比省去 了大量的对角化计算。速降法的有效性取决于波函数达到收敛的有效步数,这个步数可能很多,特别是在低 对称下,人们也做了许多改善速降法的尝试,如有更复杂的积分方程,以及改善的速降 法,及共扼梯度法等等。象速降法及共扼梯度法这样的算法在哈密顿量方ks作用到&的每步计算中都需要适当 的正
16、交化操作。如果仇是用平而波展开的(设N个单电子轨道均用M个平面波展开),则解就需要N2次浮点计算。加上正交化需要N2M次操作,计算量很大。为了减少 计算量,可以把算符Hks分成动能项与势能项,动能项在倒格矢空间是对角的,作用到 6波函数而上只需O(NM)次操作。势能匕万是局域的,则作用到波函数仇上也可方便的用 快速付立叶变换技术(FFT)在实空间计算,也只要O(NMl nM)次操作。因而工作量以线性 或O(NMl n 次增长,这与传统的标准对角化技术相比,工作量有显著减小.$2.4局域密度近似下的离子间相互作用势的计算传统的,如果给定原子坐标%,由从头计算方法推导离子间相互作用势V必,然后 进
17、行第一性原理的分子动力学模拟,一般有三个基本步骤,即在每个MD过程中:对于给定 的%求解自洽的KS方程.根据Hel l man-Feymam的力定理计算每个离子受到的力.求解牛顿运动方程=-7RlV但需要同时对电子和离子的自由度进行优化(Bendt,Zunger).利用速降法,通过下面的方 程组就可以简单的实现这样的想法.口加=一菽+E 血(2.18q).3 E=-(2.186)oRj其中约化“质量”小和方z的引入可以调节与两套参数6和 R 相关的时间标度.方 程组(2.18)就可以确定从能量而上初始点对附近的极小值之间的一条轨迹.在 R/不变 时,通过变化&,就可以找到能量而上E的唯一的极小
18、值.相应的在BO势能面上也有与之 对应的极小值;V Ez=mm Ea,R,这样能量而E上的局域极小值也就成了势能 皿面V上的局域极小值.通过同时对&,R 进行变分,方程组(2.18)就可以用于复杂的原子 结构的计算.与传统的第一性原理方法直接应用于分子动力学模拟相比,这一方法在计算量 上大为减少.不过,方程组(2.18)决定的动力学过程还有很大的局限性.它所得到的结果只是局域优化 的结果.实际上,势能而V一般有许多个极小值,为了解决这个问题1985年。ar Far r inel l o提 出一个比较满意的全局优化算法.$3从头计算分子动力学(CP方法)$3.1CP方法原理Car far r i
19、nel l o提出的从头计算分子动力学方法,最重要的一点是在真实的物理系统中 引入一个虚拟的电子动力学系统,这样组成的新系统的势能面E是离子与电子的一个总泛 函.使得这个虚拟系统产生的轨迹与具有同样势能面V的实际物理系统的轨迹一致.并能同 时处理电子-离子系统.7实际物理系统的经典拉氏量为:LCL=|/Mi瑶-(对离子品格而言)(3.1)由拉氏量1=则可以把虚拟系统的广义经典拉氏量写为:(3.2)这里,L为两套自由度&,必的泛函,它本身虽不显含时t,但&园是均与时间有关心为可 调参数.(单位:质量x长度x长度),相当于电子的“质量”实际上起着调节电子时间运动标 度的作用.可令的=4与波函数无关
20、.其中E 6,%应为电子和离子耦合的能量函数,应 包括电子动能,电子相互作用能(Har t r ee能),电子的交换关联能.以及离子实对电子的作用 能,以及离子间的相互作用能.即等价于在计算中的所需的“总能量”(即除了离子动能 项之外)。拉格朗日乘子八刀是为了保证仇的正交性而引入的。在经典力学中,就是一个 完整约束。根据拉格朗日-欧拉方程:i=1.2.3-n.从(3.2)式不难得到:b i=菽+52 人力功(3.3).dEMiRi=一万元-(3.36)Orii我们可以通过改变离子“速度”A/和电子“速度”而,可以控制这一虚拟系统的温度。从而可以对它们进行各种热力学准静态过程处理,如退热,冷却等
21、。在这一热力学过程 中,可以同时处理离子和电子的自由度。特别可以设计一个缓慢“降温”的过程,使系 统T-OK时,达到平衡状态,这一方法被称为动力学退热模拟。用动力学退热模拟来做全局优化计算的优点十分明显.在有限温度下,各自由度的速度是 不为零的,而且由于系统的复杂性,会产生随机的热力学运动.由各态经历可知,只要初始温 度足够高,体系就能达到势能面上的各个点.通过降温,体系在势能面上随机运动慢慢的局 限于几个极小值之间的跳跃.而且跳跃的频率逐渐减小.只要温度冷却的足够慢,系统最终 一定会平衡在势能而的最低点.但有时这种全局优化方案在系统即将达到能量最低点时,其 冷却速度会慢的无法忍受.以至于在实
22、际操作中,计算机无法收敛,此町就需要其他因素来 进行调节.上述我们提到的温度只是离子的“温度”,并没有涉及与电子相关的温度.在整个动力学 模拟过程中,电子自由度与离子自由度之间并不是始终处于热力学平衡状态,还可能会使按8这种方法模拟的离子运动轨迹与实际物理系统中离子的运动轨迹不一致.离子的实际运动 方程为:访的=一驾a(3.4)0R1由于V Ri=mm E6,R,相关,所以要使(3.3b)与(3.4)产生离子轨道一致,要 依求在模拟过程中,泛函始终处于对济的极小值.我们可以通过参数4和初始条 件6,戒的选取,使电子运动时间标度远小于离子运动时间标度(B。条件),以使电子在 离子坐标每次变动之前
23、尽量趋于基态,结果电子自由度在能量而E的极小值附近做简谐运 动,从而使离子的轨迹在高温和低温下基本上都保持在B0而上.即若从和6,4的选取使 离子和电子的自由耦合度很弱,它们之间能量转移是较小,而使电子的运动绝热于离子的运 动,电子可以调整自己的运动状态来“追随”离子的运动.这样除了对初始的原子构型外不 必每一步自洽求解KS方程组,从而大大减少了工作量.真正使第一性原理用于分子动力学模 拟.现在,这一方法已有效的运用到如Cl us t er,液态金属,无序系统等许多方而.$3.2基矢与原胞最初CP提出的从头计算分子动力学方法是在震势和平面波的基础上具体实现的.(当然 现在也有用缀加平面波(au
24、gment ed pl ane w av e)的算法.如软件包WIN2K.但主流是基于度 势和平面波的,主要是因为平面波基能方便的采用FFT技术,使能量,力等的计算在实空间和 倒空间快速变换,这样计算尽可能在方便的空间中进行.如前面讲的哈密顿量中动能项的矩 阵元,在倒空间中只有对角元非零.就比实空间减少了计算量.此外,平而波基函数具体形式 并不依赖于离子坐标,这样,一方而,价电子对离子作用力可以直接应用Hel l man-F6ymann定 理得到解析表示,另一方而也使总能量的计算在不同原子位型下有基本相同的精度.有利 于分子动力学模拟.此外,平而波的计算的收敛性和精确性比较容易控制,这可以通过
25、截断 能(Ecut)来控制平而波基的数量.当然平面波基也有缺点.一般电子轨道都有一定局域性,而 平面波是空间均匀的.因此,电子轨道展开时,与其它原子轨道波函数相比,平而波基的数量 要多得多,这样为了尽量减少平而波基的数量,需要引入鹰势来描写离子实与电子间相互作 用.使电子轨道波函数在离子实内部变化尽量平缓.分子动力学模拟中需要一个满足周期性边界条件的原胞(MDBo x)对于理想晶体的 计算,还是很自然的,因为哈密顿量具有平移对称性,只取其一个单胞即可.但对于无序系 统(无定型结构的固体与液体,表而,界而缺陷及杂质材料等问题,就必须把原胞取的足够 大,Mako v与P ayne已经证明:对于非周
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 密度 理论 从头 计算 分子 动力学
![提示](https://www.zixin.com.cn/images/bang_tan.gif)
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【曲****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【曲****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。