第性原理和密度泛函理论教育课件.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,第一性原理和密度泛函理论PPT讲座,引 言,量力力学是反映,微观粒子,(,分子、原子、原子核、基本粒子,),运动规律的理论。,以量子力学为基础，结合高速发展的计算技术分别建立起来的计算材料科学、计算物理、量子化学等分枝学科，促进了物理学、化学和材料科学的发展，为发展和设计新型材料提供了理论基础和,新的,研究方法。,在理论上最具诱惑力，且在将来最有可能开展真正意义上的材料设计的计算就是解体系的,Schrdinger,方程，即为计算材料学中的第一原理计算。,多粒子体系的第一原理,第一原理,的基本思想：,将多原子构成的体系理解为由电子和原子核组成的多粒子系统，在解体系,Schrdinger,方程的过程中，最大限度地进行“非经验性”处理，即不涉及任何经验参数，所要输入的只是原子的核电荷数和一些模拟环境参量。,计算所求得的结果是体系,Schrdinger,方程的本征值和本征函数（波函数），有了这两项结果，就可研究体系的基本物理性质。,多粒子体系的第一原理,从微观角度看，一块物质材料是由大量（每立方厘米约,10,23,个）,原子核和游离于原子核之间的,电子,组成。因此，材料的性质（如硬度、电磁和光学性质）和发生在固体内的物理和化学过程是由它所包含的原子核及其电子的行为决定的。,？,20,世纪初量子力学的出现，原则上提出了像原子核和电子这样的微观粒子运动和交互作用的定律,。,理论上，给定一块固体化学成分（即所含原子核的电荷和质量），我们就可以计算这些固体的性质。因为一块固体实际上是一个,多粒子体系,。决定这个体系性质的波函数可以通过解薛定谔（,Schrdinger,）波动方程来获得。,多粒子体系的第一原理,多粒子系统的,Schrdinger,方程,其中,和,H,分别对应于多粒子系统的波函数和哈密顿量。,原则上只要对,上式进行求解即可得出所有物理性质，然而由于电子之间的相互作用的复杂性，要严格求出多电子体系的,Schrdinger,方程解是不可能的,，,必须在物理模型上进一步,作一系列的近似,。,多粒子体系的第一原理,换句话说，做第一原理计算（,ab initio calculation,）便可知道一块固体的性质。,可是，这个薛定谔波动方程有,3N,个变量（,N,是粒子总数），极其复杂，假使我们把目前世界上的所有电脑都用上，让它跑千年、万年都不可能算出来。正如,1929,年量子物理大师,狄拉克,（,Dirac,）所言：,处理大部分物理学和全部化学问题的基本定理已经完全知道。困难在於这些定律的应用所引出的数学方程（,Schrdinger,方程）太复杂以致於无法解决。,伤不起啊！,多粒子体系的第一原理,价电子近似,在,原子模型中，近核电子可被紧紧,束缚,在核的周围，这些束缚电子是定域的，比较稳定，因而对固体性质的贡献很小。而在外层轨道的价电子可以是离域的，当原子结合在一起组成固体时，这些电子的状态变化很大，对固体的电学和光学性质有决定性的影响。因此，可将固体看作是,由原子核,束缚电子构成的,离子核,和,价电子,组合而成,。,多粒子体系的第一原理,绝热,近似,由于,原子核的质量远大于电子的质量，原子核的运动速度要比电子慢很多，因此可以认为电子运动在固定不动的原子核的势场中，所以原子核的动能为零，而势能为一个常数,。,中子,/,质子的质量是电子质量的约,1835,倍，即电子的运动速率,比核,的运动速率要高,3,个数量级，因此可以实现电子运动方程,和核,运动方程的,近似脱耦,。这样，电子可以看作是在一组准,静态原子核,的平均势场下运动。,为了对问题进行简化，可以把多电子系统中的相互作用视为有效场下的无关联的单电子的运动。一个具有,N,个电子的系统的总波函数可以写成所有单电子的波函数的乘积,。,Hartree,方程,Hartree-Fock,近似,Hartree-Fock,近似,Hartree-Fock,近似,单粒子,Hartree-Fock,算符是自恰的，即它决定于所有其他单粒子,HF,方程的解，必须通过叠代计算来求解，具体操作如下：,猜测试探波函数,构造所有算符,求解单粒子赝薛定谔方程,对于解出的新的波函数，重新构造,Hartree-Fock,算符,重复以上循环，直到收敛（即前后叠代的结果相同）,自恰场（,SCF,）方法是求解材料电子结构问题的常用方法,对处理原子数较少的系统来说，,Hartree-Fock,近似是一种很方便的近似方法。,但用于,原子数大,的系统，问题就变得非常,复杂,，此计算方法的计算量随着电子数的增多呈,指数增加,，这种计算对计算机的内存大小和,CPU,的运算速度有着非常苛刻的要求，它使得对具有较多电子数的计算变得不可能。同时,Hartree-Fock,近似方法给出的一些金属费米能和半导体能带的计算结果和实验结果偏差较大。,这在很大程度上导致了,密度泛函理论,的产生。,Hartree-Fock,近似,Hartree-Fock,方法的主要缺陷,完全,忽略,电子关联作用,计算量,偏大，随系统尺度,4,次方关系增长,Density Functional Theory,(DFT 1964),一种用电子密度分布,n,(,r,),作为基本变量，研究多粒子体系基态性质的新理论,密度泛函理论,W.Kohn,W.Kohn,荣获,1998,年,Nobel,化学奖,密度泛函理论,1927,年，,Thomas-Fermi,理论首先运用,作为计算原子中电子结构的基本变量，从而给出一个简化处理方法。,如果能够证明，对于任意电子系统的基态电子密度的分布唯一地决定电子系统的情况，则基于的一个计算电子结构的新的理论框架就可以建立起来,，,这就是表述密度泛函理论,（,density functional theory,，简称,DFT,）,的,出发点,。,密度泛函理论,在密度泛函理论中，将电子密度作为描述体系状态的基本变量，可追溯到,Thomas,和,Fermi,用简并的非均匀电子气来描述单个原子的多电子结构。,直到,Hohenberg,和,Kohn,提出了两个基本定理才奠定了密度泛函理论的基石。,随后,Kohn,和,Sham,的工作使密度泛函理论成为实际可行的理论方法。,密度泛函理论,1964,年,Hohenberg,和,Kohn,推证（,P.Hoheuberg and W.Kohn,Physical Review 136B,864,（,1964,）,）,：一个多粒子（如电子）体系性质（如体系的能量）由粒子密度的空间分布（即空间任意一点的粒子平均个数）完全决定。换句话说：,我们无需去求解描述每个粒子运动的体系波函数，我们只要设法找出仅有,3,个变量的粒子密度的空间函数就行了。,1998,获奖诺贝尔化学奖,表彰,Walter Kohn,在,60,年代提出,密度泛函理论,及,John A.Pople,发明了测验化学结构和物质特性的,计算机技术,密度泛函理论,1965,年柯恩又和沈吕九证明,(,W.Kohn and L.J.Shan,Physical Review 140,All33,),：一个多粒子体系的粒子密度函数可以通过一个简单的单粒子波动方程获得。这个单粒子波动方程现在被称作柯恩沈（,Kohn-Sham,）方程。,Hohenberg,，,Kohn,和,Shan,的理论就是诺贝尔化学奖颁词所指的密度泛函理论，显然，密度泛函理论大大简化了应用量子力学探讨材料物理性质所涉及的数学问题。,W.Kohn,Lu.J.Sham,（沈吕九）,密度泛函理论,然而，当时用密度泛函理论计算一块固体材料的性质仍然有相当的困难，到十年过后，电脑运算能力的提高和有效数值方法的发明和发展，第一原理计算才成為现实。,20,世纪,70,年代中期，密度泛函理论首先在固体物理学领域的应用，定量地预测了的大量简单金属材料的结构参数、磁学性质等等。,密度泛函理论,密度泛函理论,物质电子结构的新理论,1.,氢原子,1,）,Bohr:,电子粒子,2,）,Schrodinger:,电子波,(,r,),.,3,）,DFT:,电子是电子云的密度分布。,n,(,r,).,密度泛函理论,3,）,DFT:,电子是电子云的密度分布。,2.DFT,中的氢分子。,由密度分布表示。,密度泛函理论,3.,大分子（例如,DNA,）,N,个原子,Schrodinger:(r,1,r,2,r,3,r,N,)3N,维空间,DFT:,n,(,r,),3,维空间,Hohenberg-Kohn,定理,I,1964,年，,P.Hohenberg,和,W.Kohn,在,非均匀电子气,理论的基础上，提出两个基本定理，奠定了密度泛函理论的基础。,定理,1,：,对于一个共同的外部势,v,(,r,),相互作用的多粒子系统的所有基态性质都由（,非简并）,基态的电子密度分布,n,(,r,),唯一地决定,。,或,:,对于非简併基态，粒子密度分布,n,(,r,),是系统的基本变量。,密度泛函理论,考虑一个多粒子系（电子体系、粒子数任意），在外部势和相互作用,Coulomb,势作用下，,Hamiltonian,为,密度泛函理论,电子密度算符,电子密度分布,n,(,r,),是 的期待值：,Hohenberg-Kohn,定理,II,定理,2,：,如果,n,(r),是体系正确的密度分布，则,E,n,(r),是最低的能量，即体系的基态能量,密度泛函理论,实际计算是利用,能量变分原理,，使系统能量达到最低（有一定精度要求）。由此求出体系的真正电荷密度,n,(,r,),进而计算体系的所有其它基态性质。如能带结构，晶格参数，体模量等等。,密度泛函理论,Hohenberg-Kohn,定理指出，多粒子系统的所有基态性质都是密度泛函的唯一泛函。因此系统的总能量,可表达为,其中,T,n,代表无相互作用的粒子的动能，,U,n,代表一般的库仑能，,E,xc,n,代表多体效应的交换关联能,。,假设,N,个电子的系统有,N,个单电子波函数,，则电荷密度为：,动能可以写为,密度泛函理论,库仑能可具体表达为：,此时总能量,E,t,的表达式可写为,交换关联能密度,密度泛函理论,由上式可得到一个本征方程,即,Kohn-Sham,方程,交换关联势,求解本征值,总能量的另一表达式,密度泛函理论,密度泛函理论,密度泛函理论,Kohn-Sham,近似的核心思想：,1,.,动能的大部分通过相同电子密度的无相互作用体系来计算,2.,电子相互作用中库仑作用占据了主要部分，而交换相关是相对次要的,3.,非经典的交换和相关作用，动能校正项，自相互作用折入交换相关泛函中,解,Kohn-Sham,方程的流程图,n(r)=,n,at,(r),求解,、,V,xc,、,V,eff,求解,Kohn-Sham,方程,得到,i,由,i,构造,n,out,(r),比较,n,in,与,n,out,(r),计算总能,E,tot,No,Yes,n,in,与,n,out,混合,原子计算,精度控制,No,Yes,输出结果：,E,tot,、,i,、,n(r),V,xc,、,V,eff,、,E,n,(k),、,N(E),33,虽然,K-S,方程十分简单，其计算量也只有,Hartree,方程的水平，但却包含着深刻得多的物理内容。其中一个重要的,概念性结果,是，,多体基态的解被准确地简化为基态密度分布之解，而这个密度是由单粒子的,Schrdinger,方程给出的,。,由此，方程中的,有效势,在原理上,包括了所有的相互作用效应，即,Hartree,势、交换势（由,Pauli,原理决定的相互作用所产生的势）和关联势（一个给定的电子对整个电荷分布的影响所产生的势）,。,在这个意义上，它比,Hartree-Fock,方程要优越得多。,密度泛函理论,34,Formally,equivalent,Electron,Interaction,External potential,Hard problem to solve,Schr,dinger equation,“Easy”problem,To Solve DFT,Properties of the,system,Non-interacting electron,(KS particle),Effective potential,LDA,GGA,etc,量子力学体系的性质可以通过,求解薛定格方程,(SE),进行计算,(,上图左边,),。,但更加容易的、形式上等价的方法是,求解,DFT,的,KS,方程,(,上图右边,),。,但是准确的,E,xc,n(r),并不知道。需要采用近似方法，如,LDA or GGA,。,这就会影响,KS,解的精度。,SE,DFT,电子运动服从量子力学规律,电子体系的性质由其状态波函数确定,波函数包含,3N,个变量,(N,为电子数目,),对于含很多电子的大体系,通过求出波函数来计算体系的性质其计算量非常大,很难实现,。,根据密度泛函理论,体系的性质由其,电子密度分布,唯一确定,电子密度分布是只含,3,个变量的函数,通过它研究体系的性质可以大大减少计算量,对大体系的量子力学计算就比较,容易进行,。,密度泛函理论研究的基本内容,寻找体系的性质,(,特别是,动能,和,交换关联能,),作为电子密度分布的泛函的精确或近似形式、相关的计算方法和程序以及在各科学领域的应用。,密度泛函理论,密度泛函理论,它提供了第一性原理或从头算的计算框架。在这个框架下可以发展各式各样的能带计算方法。,在凝聚态物理中，如：材料电子结构和几何结构，固体和液态金属中的相变等。,这些方法都可以发展成为用量子力学方法计算力的精确的分子动力学方法。,密度泛函理论,DFT,适应于大量不同类型的应用：,(1),电子基态能量与原子（核）位置之间的关系可以用来确定分子或晶体的结构；,(2),当原子不处在它的平衡位置时，,DFT,可以给出作用在原子,(,核,),位置上的力。,2.,因此，,DFT,可以解决原子分子物理中的许多问题，如,(1),电离势的计算，,(2),振动谱研究，,(3),化学反应问题，,(4),生物分子的结构，,(5),催化活性位置的特性等等。,3.,另一个重要优点是降低维数（,Kohn,的演讲）,密度泛函理论,要运用,密度泛函理论,来解决实际问题，还需要进一步确定交换关联泛函,E,xc,n,(,r,),的近似表达式。目前最流行的交换关联泛函是局域密度,近似（,LDA,）和广义梯度密度近似（,GGA,）,HK,定理已经建立了密度泛函理论（,DFT,）的框架，但在实际执行上遇到了严重困难。主要是,相互作用,电子体系的,交换关联能,E,xc,n,无法精确得到。,为了使,DFT,理论能够付诸实施，,Kohn-Sham,提出了,局域密度近似,(Local Density Approximation,LDA,),。,密度泛函理论,局域密度近似（,LDA,）,密度泛函理论,DFT,本身未做任何近似，对基态严格成立，但是，除了均匀电子气体以外，对于真实的体系并不知道交换关联势,xc,的具体形式。,局域密度近似（,LDA,）是所有的近似交换关联泛函,E,xc,的基础。考虑一个电子密度缓变的系统，把,E,xc,写成局域量,xc,n,(,r,),的积分形式：,其中,xc,n,(,r,),是每个粒子的交换关联能，可进一步分解成交换能,x,和关联能,c,两个部分,密度泛函理论,对于均匀电子气，,xc,n,(,r,),和,交换关联势,xc,可以通过忽略电子关联得到：,LDA,在自旋极化的系统中，,LDA,演化为局域自旋密度近似（简称,LSDA,）,e,e,交换关联泛函表达式,密度泛函理论,尽管,L(S)DA获得了巨大的成功，但是也有许多不足之处，如系统地高估结合能。在,L(S)DA,基础上的改进有广义梯度近似（,GGA,）,在,GGA,近似下，交换相关能是电子（极化）密度及其梯度的泛函。,密度泛函理论,广义梯度近似（,GGA,）,局域（自旋）密度近似是在均匀电子气或电子密度变化足够缓慢的系统中提出的，用于描述密度变化大的非均匀电子气系统并不适合。因此在,LDA,的基础上，引入电子密度的梯度展开因子,，提出了广义梯度密度近似（,GGA,），,迄今已经发展出很多种,GGA,的具体形式，常见的有,Langreth-Mehl,（,LM,）,Perdew-Wang 91,（,PW91,）,Perdew-Burke-Ernzerhof,（,PBE,）,广义梯度近似（,GGA),构造,GGA,交换相关泛函的方法分为两个流派。,Becke,为首的一派,认为“,一切都是合法的,”，所以人们可以以任何原因选择,任何可能的泛函形式,，而这种形式的好坏由实际计算来决定。通常，这样的泛函的参数由,拟合大量的计算数据,得到。,以,Perdew,为首的一派,他们认为发展交换相关泛函必须以一定的,物理规律为基础,，这些规律包括标度关系、渐进行为等。基于这种理念构造的一个著名的,GGA,泛函是,PBE,泛函，也是现在用的最广泛的,GGA,泛函之一,不同的,LDA,方案之间大同小异，但不同的,GGA,方案可能给出完全不同的结果。,密度泛函理论,LDA/GGA,近似下密度泛函理论的局限与改进,在交换关联能中考虑的是,电子电子,之间的相互作用，但要求出系统的总能量，还必须计算出,电子原子核,的库仑能。,运用,价电子近似法,，多电子,Schrdinger,方程将大为简化。此时可以忽略大多数的近核电子与原子核的库仑作用，而只需关心价电子和离子芯的相互作用。,Fermi,首先在,1934,年提出用一种有效势场，即,赝势,，来近似地描述价电子所处的势场。,密度泛函理论,赝势,密度泛函理论,赝势,赝势方法从量子力学证明了为什么大多数金属中价电子能够看作象自由电子一样在具有离子芯的非常强的晶体势场中运动,赝 势,构造赝势的具体方法有：,经验赝势（,Empirical Pseudo Potentials,，简称,EPP,）,模,方守恒赝势（,Norm Converving Pseudo Potential,，简称,NCPP,）,超,软赝势（,UltraSoft Pseudo Potentials,，简称,USPP,）,其中,NCPP,和,USPP,是基于密度泛函理论的,第一原理赝势,。,密度泛函理论,赝势,密度泛函理论,NCPP,的方法是在原子中指定一个芯区半径,r,c,，称之为截断半径，在,r,c,以外赝波函数,v,和真实的价电子波函数,v,相同，而在,r,c,之内,v,和,v,有相同的模方。,其中,R,l,r,),是波函数的径向部分，,PP,代表赝波函数，,AE,代表全电子（,All-Electron,）的真实波函数，下标,l,表示角动量，,n,表示价电子能级。此外还有两个强制性的条件：赝波函数不应出现节面；赝能量本征值和价电子的真实能量本征值相同,赝势,USPP,与,NCPP,不同，为了能够得到尽可能平滑的赝势并减少波函数的数目，由,Vanderbilt,提出的超软赝势（,USPP,）模型放弃了模方守恒条件。提出,USPP,的基本思路是：在大多数情况下，只有当紧束缚电子的大部分轨道定域在芯区时，才需要高的截断能。因此可以把这些电荷从芯区移去，来达到平滑赝波函数和降低截断能的目的。为了使总电荷数不变，引入缀加函数相应增加了芯区的电子密度。,赝势,赝势,厄密重叠算符,Kohn-Sham,方程可写成,通过自洽求解,方程可,得,赝波函数,和电子密度,n(r),其它物理性质,与密度泛函理论相关的研究主要有以下三方面的工作,:,密度泛函理论本身的研究,寻找基态体系性质,(,特别是,动能,和,交换相关能,),作为电子密度分布的泛函的精确形式或者尽可能精确的近似形式。,拓宽密度泛函理论的内涵,。,密度泛函理论,(2),密度泛函计算方法的研究,新算法和程序的优化,对于大的体系,计算很复杂,是能否用密度泛函理论方法进行研究的,瓶颈,。因此,发展高效率的计算方法和相关程序是很重要的工作。,目前的研究热点是实现对,大体系,的高精度计算,其中结合使用密度泛函理论的线性标度算法和分区算法特别受到重视,迄今对此也已经提出过很多算法,并且推出了相关的计算程序,发展对含重元素体系的相对论密度泛函计算方法也受到重视。,密度泛函理论,密度泛函理论,(3),用以近似能量密度泛函为基础建立的方法研究各种化学和物理问题,密度泛函方法由于其计算量比从头计算方法小得多,可以用来计算大的复杂体系,结果精度可以满足很多研究工作的要求,因此目前已经得到广泛应用。,随着更精确的密度泛函形式的发现和更高效率的计算方法和程序的推出,密度泛函理论方法将在,化学、物理学、材料科学,(,纳米科学,),、生命科学、药物化学,等领域的研究工作中发挥更大的作用。,密度泛函理论的应用实例,原子总能,LDA,计算结果相对较小，,GGA,与实验值更接近,密度泛函理论的应用实例,小分子的结合能,LDA,计算结果通常偏大，,GGA,计算结果与实验值较为接近,密度泛函理论的应用实例,固体晶格常数,与实验值相比，,LDA,结果偏大，,GGA,偏小，但,GGA,的误差更明显,密度泛函理论的应用实例,固体弹性模量,LDA,结果偏硬，,GGA,偏软,密度泛函理论的应用领域,随着密度泛函理论的发展，它的应用领域越来越广泛，在物理、化学和生物等多门学科中，密度泛函理论已成为强有力的研究工具。,密度泛函理论研究涉及的体系,零维（如小分子、团簇、量子点一,维（如纳米管,）,二,维（如固体,表面）,三维,（如高温,超导）,物理学：强相关体系,传统：,使用模型哈密顿量处理，通过量子,蒙特卡洛,模拟数值求解,计算能力增加,DFT+U,DFT-DMFT,强相关体系,密度泛函理论,随着计算能力的增加和,DFT+U,和,DFT-DMFT,等方法的出现，现在已经可以在基于密度泛函的理论框架下处理强关联体系。,化学：弱作用体系,范德瓦尔斯（,vdW,）力或者说长程色散力对松散堆积的软物质、惰性气体、生物分子和聚合物，以及对物理吸附等过程都非常关键,传统的密度泛函理论中对长程密度涨落效应考虑的不足，基于,LDA/GGA,的密度泛函理论被发现不能很好的处理这类体系，以前这种弱作用体系都是通过半经验的参数拟合分子力场的方法处理。,发展密度泛函的,vdW,理论的目标是得到一个既能产生,vdW,相互作用系数又能产生总关联能的非局域泛函。,Kohn,等人和,Lein,等人分别提出了一种对不同系统间距都适用的计算方案。,既可以计算范德华系数，又能准确的计算总关联能，被称为无缝的（,seamless,）方法。,Dobson,和,Wang,成功的将这种方法用于计算任意间隔的两个,jellium,平板之间的作用力,化学：弱作用体系,密度泛函加衰减色散（,DFdD,）方案,为了处理大的复杂的体系，需要通过较小的计算量得到比较精确的结果，所以人们希望对,DFT,加以简单的修正来计算,vdW,体系。,Wu,（,X.Wu,et al.,J.Chem.Phys.115,8748(2001).20,）等人提出了密度泛函加衰减色散（,DFdD,）方案。,DFdD,方案对,H,2,O,2,和,(C,6,H,6,),2,等体系的成功说明了它用于诸如水、表面物,理吸附和涉及氢键的化学反应的可能性。,生命科学：空间尺寸与时间尺度,密度泛函理论在生物体系中的应用的障碍,:,空间尺度,时间尺度,体系庞大、生物环境、生物过程（时间）,要完全理解生物体系,又必须在,原子尺度,考察其中,的化学物理,过程,为了克服生物体系的尺寸限制,，最,简单的办法就是把与要研究的,过程,关系最密切的部分拿出来进行,第一性原理,的研究,。,周围,结构的长程,作用可以,近似成,静电场,，而溶液环境可以通过,自洽反应场,（,SCRF,）近似成,介电,连续,媒质,。,QM/MM,方法，即,在一个用经典的分子力学描述的大系统中嵌入,量子力学,描述的,核心片断,。,量子和经典区域的长程相互作用可以用,静电场,来,描述,描述区域界面打断的共价键：,饱和,原子,法,断键外围加上,氢,原子来饱和悬,键,&,冻结,轨道,法,在每个界面原子中,包,括一个局域的杂化冻结轨道,生命科学：空间尺寸与时间尺度,时间尺度,通常可以通过一些简单势能面方法（如,LST,和,QST,）来研究反应能量和寻找过渡态。,通过基于过渡态理论的方法来扩展模拟的时间尺度。,并行复制（,parallel replica,）,动力学,对,势能面进行修正的超,动力学,（,hyperdynamics,）,赝,动力学（,metadynamics,）,温度,加速,的（,temperature-accelerated,）,动力学,分子电子学：输运性质,随着电子学器件越变越小，其尺寸将很快逼近,原子分子尺度,。,早期对介观输运性质研究使用简单的有效质量理论和,Boltzmann,方程，没有考虑器件的电子结构。,随着密度泛函理论与输运理论的发展和结合，人们开始可以对分子体,系输运性质进行第一性原理研究。,分子器件中存在电流和电压：非平衡态体系。,典型的处理方法：散射态方法和非平衡格林函数方法。,Kosov,通过利用流限制的,拉格朗日乘子变分法,得到在某个确,定的电流强度下的,薛定谔方程,和,KS,方程，将一个开放体系的模,拟变成一个等效的封闭体系。由于采用简单的拉格朗日表述形式,，这种方法有可能被用来进行载流体系的分子动力学模拟，从而,处理电流与分子振动耦合等复杂的问题,光谱学：激发态和外场,电子激发往往是各种谱测量的基础，如光电子谱、电子能损谱和吸收光谱等。,现在基于,DFT,激发态问题的处理方法有许多种，包括,系综密度泛函理论,，或者化学中用得较多的考虑系统对称性，用求和方法计算多重态激发能的方法。,更主流的两种方法：,多体微扰理论,基于一系列格林函数方程，从,单电子传播,子开始考虑电子空穴格林函数响应。其中的核心物理量是电子自能和电子空穴相互作用。电子自能可以用密度泛函计算结果为零阶近似。,含时密度泛函理论,比多体微扰理论的计算量要小得多,人们开始考虑用含时密度泛函理论计算激发态时的非绝热效应。,小 结,DFT,是当今处理相互作用多电子体系电子结构和几何结构最有力的工具。所谓从头算或第一性原理方法就是基于,DFT,框架建立起来的。它独立于实验，只需很少几个熟知的基本物理参数便可运作。,DFT,并不要求原子的周期性排列，它具有十分广泛的适应性。已经在计算凝聚态物理、计算材料科学、量子化学、量子生物学和许多工业技术部门获得成功的应用。,71,3,.DFT,最基本的应用依赖于,LDA,近似，,DFT-LDA,理论在原子、分子和（有机、无机）固体的基态性质研究中同样获得巨大成功。,4.DFT-LDA,理论也存在明显的,缺点,。在处理半导体材料的电子结构中，存在著名的,带隙偏小问题,；它在研究有电子,强关联,效应的材料，如高温超导材料、某些磁性材料和过渡金属氧化物等的电子性质时，也存在严重问题；对于,超大体系,（如纳米材料和生物体系）的研究仍然需要发展新的辅助方法和提高计算机能力。,5.,这些缺点的来源主要是,DFT-LDA,作为基态理论不能满意地描述激发态造成的。其中最重要的是必须解决激发态的交换关联能问题。,6.,最近已经提出了不少超越,LDA,的理论方法，它们正朝着解决上述各种困难深入发展。如准粒子,GW,方法，屏蔽交换局域密度近似（,sX-LDA,）和含时间密度泛函理论,(TDDFT),等。,近几年来，,DFT,同分子动力学方法相结合，有许多新发展,在材料设计、合成、模拟计算和评价诸多方面有明显的进展,已成为计算凝聚态物理、计算材料科学和计算量子化学的重要基础和核心技术,在工业技术领域的应用开始令人关注,THE END,