材料物理1现代电子理论

1、12化学反应位能面早在早在1950 1950 年代，波普尔年代，波普尔( (PoplePople J A) J A)就对半经验就对半经验量子化学方法作出重要量子化学方法作出重要贡献，他是贡献，他是Pariser-Pariser-Parr-PopleParr-Pople 方法方法(PPP)(PPP)的创始人之一，的创始人之一，PPP PPP 可可用来研究不饱和有机分用来研究不饱和有机分子的电子结构。子的电子结构。以后他将这一以后他将这一电子理论扩展为所有价电子电子理论扩展为所有价电子的理论称为的理论称为CNDOCNDO。开始时，波普尔并不相信。开始时，波普尔并不相信非经验的从头计算法能够和半经验

2、方法竞争，非经验的从头计算法能够和半经验方法竞争，后来他自己改变了这种状况。后来他自己改变了这种状况。3哈特里哈特里- -福克方法福克方法，虽然引入了基函数虽然引入了基函数和高斯型基函数，离实用仍有很大距和高斯型基函数，离实用仍有很大距离，瓶颈在于如离，瓶颈在于如H-FH-F表达式那样的至少表达式那样的至少是六维的多中心积分，这种积分在计是六维的多中心积分，这种积分在计算中可达算中可达106106至至109 109 个。波普尔发明了个。波普尔发明了一种方法，可使计算工作量减少一到一种方法，可使计算工作量减少一到两个数量级。他构造了一系列质量逐两个数量级。他构造了一系列质量逐步改进的基函数，并进

3、行了大量的计步改进的基函数，并进行了大量的计算加以检验，这是一个决定性的成就，算加以检验，这是一个决定性的成就，由于这种改进才使哈特里由于这种改进才使哈特里- -福克方法可福克方法可以应用于实际的化学问题。以应用于实际的化学问题。4为了进一步得到分子的结构，人为了进一步得到分子的结构，人们不仅要计算一定构型的能量，们不仅要计算一定构型的能量，为了求得平衡的构型，还要知道为了求得平衡的构型，还要知道能量的一阶乃至二阶导数，波普能量的一阶乃至二阶导数，波普尔在早期尔在早期PulayPulay P P 工作的基础上，工作的基础上，开发了有效计算这些导数的方法。开发了有效计算这些导数的方法。5这种方法

4、不仅能研究平衡态，还这种方法不仅能研究平衡态，还可以沿着反应的途径，确定能垒，可以沿着反应的途径，确定能垒，以及过渡态的结构。波普尔在以及过渡态的结构。波普尔在1970 1970 年将这些进展开发成为一年将这些进展开发成为一个应用软件个应用软件GUASSIANGUASSIAN，并不断做，并不断做出改进，出改进，1990 1990 年代又将年代又将密度泛密度泛函方法函方法引入。这一软件已为化学引入。这一软件已为化学界广泛采用。界广泛采用。6“泛函泛函”这个名词是由法国数学家阿达这个名词是由法国数学家阿达马（马（HadamardHadamard，1865-19631865-1963）在）在1897

5、1897年研年研究变分问题时引进的。究变分问题时引进的。“泛函泛函”也称泛也称泛函数，它是对实（复）值函数概念的拓函数，它是对实（复）值函数概念的拓广或发展，通俗地说，泛函就是以函数广或发展，通俗地说，泛函就是以函数为变元的函数，其基本思想是把函数为变元的函数，其基本思想是把函数（或曲线等）看作空间的元素或点，而（或曲线等）看作空间的元素或点，而函数的集合构成了空间，函数的集合构成了空间，“泛函分析泛函分析”是研究无限维线性空间的拓扑性质及其是研究无限维线性空间的拓扑性质及其“泛函泛函”与与“算子算子”的一般性质的一个的一般性质的一个现代数学分支，它是无限维分析学的一现代数学分支，它是无限维分

6、析学的一个重要组成部分。个重要组成部分。 7虽然上述方法只计及电子交换，而忽略了电子间的相关，但对能量计算来说，误差已不足1%。为了达到化学准确性的要求，以更好地预测键能、活化能、反应焓等，他又进一步研究电子相关，采用微扰理论达到四阶，还进行波函数的集团展开，进行了更精确的计算。由于上述成就，波普尔在1998 年荣获Nobel 化学奖。89用密度泛函用密度泛函理论计算在理论计算在钯催化剂上钯催化剂上合成醋酸乙合成醋酸乙烯酯，上部烯酯，上部是结构变化，是结构变化，下部是能量下部是能量变化变化1011沃尔特沃尔特库恩库恩的密度泛函理论对化学作出的密度泛函理论对化学作出了巨大的贡献。量子化学理论和计

7、算的丰了巨大的贡献。量子化学理论和计算的丰硕成果被认为正在引起整个化学的革命。硕成果被认为正在引起整个化学的革命。量子化学家几十年的辛勤耕耘得到了充分量子化学家几十年的辛勤耕耘得到了充分的肯定。的肯定。这标志着古老的化学已发展成为这标志着古老的化学已发展成为理论和实验紧密结合的科学。沃尔特理论和实验紧密结合的科学。沃尔特库库恩的密度泛函理论构成了简化以数学处理恩的密度泛函理论构成了简化以数学处理原子间成键问题的理论基础原子间成键问题的理论基础, ,是目前许多是目前许多计算得以实现的先决条件。传统的分子性计算得以实现的先决条件。传统的分子性质计算基于每个单电子运动的描写质计算基于每个单电子运动的

8、描写, ,使得使得计算本身在数学上非常复杂。计算本身在数学上非常复杂。12沃尔特库恩指出,知道分布在空间任意一点上的平均电子数已经足够了,没有必要考虑每一个单电子的运动行为。这一思想带来了一种十分简便的计算方法密度泛函理论。方法上的简化使大分子系统的研究成为可能,酶反应机制的理论计算就是其中典型的实例,而这种理论计算的成功凝聚着无数理论工作者30余年的心血。如今,密度泛函方法已经成为量子化学中应用最广泛的计算方法。College of Technology, Portsmouth, England（1955-1957）B. Sc., University of London（Imperial

9、College）, England（1957-1960）Ph. D., University of Cambridge, England（1960-1963）University of California, San Diego. Research Associate（1963-1966）、Associate Professor（1968-1974）、Professor（1975 ）. University of California, Irvine（1966-1967）. Assistant Professor. Reader. University of London（Queen Mary

10、 College）（1967-1968）.Research Physicist. IBM Research Center, Yorktown Heigts（1974-1975）. Dean. Division of Natural Sciences（1985-1989）. Director. Institute of Pure and Applied Physical Sciences（1991-1995）. Chairman. Department of Physics, UCSD（1995-1998）14 DFTDFT的优点的优点 HohenbergHohenberg-Kohn-Kohn定

11、理定理 能量泛函公式能量泛函公式 局域密度近似（局域密度近似（LDALDA） Kohn-Sham Kohn-Sham方程方程 总能总能E Etottot表达式表达式 DFT DFT的意义的意义概述概述DFTDFT即即Density Density Functional Theory Functional Theory 一种用电子密度分布一种用电子密度分布n n( r)( r)作为基本变量，研究多粒子作为基本变量，研究多粒子体系基态性质的新理论。体系基态性质的新理论。 W. Kohn W. Kohn 荣获荣获19981998年年Nobel Nobel 化学化学奖奖DFT- LDA- KSDFT-

12、 LDA- KS 自从自从2020世纪世纪6060年代（年代（19641964）密度泛函理论（）密度泛函理论（DFTDFT）建立）建立并在局域密度近似（并在局域密度近似（LDALDA）下导出著名的）下导出著名的KohnKohnSham Sham ( (沈呂九沈呂九)(KS)(KS)方程以来，方程以来，DFTDFT一直是凝聚态物理领域计一直是凝聚态物理领域计算电子结构及其特性最有力的工具算电子结构及其特性最有力的工具。15近几年来，近几年来，DFTDFT同分子动力学方法相结同分子动力学方法相结合，有许多新发展；合，有许多新发展；在材料设计、合成、模拟计算和评价在材料设计、合成、模拟计算和评价诸多

13、方面有明显的进展；诸多方面有明显的进展；已成为计算凝聚态物理、计算材料科已成为计算凝聚态物理、计算材料科学和计算量子化学的重要基础和核心学和计算量子化学的重要基础和核心技术；技术； 在工业技术领域的应用开始令人关注。在工业技术领域的应用开始令人关注。16它提供了第一性原理或从头算的计算它提供了第一性原理或从头算的计算框架。在这个框架下可以发展各式各框架。在这个框架下可以发展各式各样的能带计算方法。样的能带计算方法。在凝聚态物理中，如：在凝聚态物理中，如： 材料电子结构和几何结构，材料电子结构和几何结构，固体和液态金属中的相变等。固体和液态金属中的相变等。171.1.DFTDFT适应于大量不同类

14、型的应用：适应于大量不同类型的应用： (1)(1)电子基态能量与原子（核）位置之间的关系可以电子基态能量与原子（核）位置之间的关系可以用来确定分子或晶体的结构；用来确定分子或晶体的结构； (2)(2)当原子不处在它的平衡位置时，当原子不处在它的平衡位置时，DFTDFT可以给出作可以给出作用在原子用在原子( (核核) )位置上的力。位置上的力。2. 2. 因此，因此，DFTDFT可以解决原子分子物理中的许多问题，如可以解决原子分子物理中的许多问题，如 (1)(1)电离势的计算，电离势的计算， (2)(2)振动谱研究，振动谱研究， (3)(3)化学反应问题，化学反应问题， (4)(4)生物分子的结

15、构，生物分子的结构， (5)(5)催化活性位置的特性等等。催化活性位置的特性等等。3. 3. 另一个重要优点是降低维数（另一个重要优点是降低维数（KohnKohn的演讲）的演讲）1819密度泛函理论密度泛函理论物质电子结构的新理论物质电子结构的新理论1 1。氢原子。氢原子1 1）Bohr: Bohr: 电子粒子电子粒子2 2）SchrodingerSchrodinger: : 电子波电子波 (r(r) .) .3 3）DFT: DFT: 电子是电子云的电子是电子云的密度分布。密度分布。 n n（r r）. .20 2 2。DFTDFT中的氢分子。中的氢分子。 由密度分布表示。由密度分布表示。2

16、13 3。大分子（例如。大分子（例如DNADNA）; ; N N个原子。个原子。SchrodingerSchrodinger: : (r(r1 1,r,r2 2,r,r3 3,r rN N) )，3N3N维空间。维空间。DFT: DFT: n(rn(r) )， 3 3维空间。维空间。在有机化学、生物技术、合在有机化学、生物技术、合金物理、表面科学、磁性等金物理、表面科学、磁性等领域领域DFTDFT最为重要。最为重要。1. 1. 定理定理1 1：对于一个共同的外部势对于一个共同的外部势v v( (r r), ), 相互作用的多粒子系统的相互作用的多粒子系统的所有基态性质都由（非简併）基态的电子密

17、度分布所有基态性质都由（非简併）基态的电子密度分布n n( (r r) )唯一地决唯一地决定。或定。或: : 对于非简併基态，粒子密度分布对于非简併基态，粒子密度分布n n( (r r) )是系统的基本变量。是系统的基本变量。2. 2. 考虑一个多粒子系（电子体系、粒子数任意），在外部势和相考虑一个多粒子系（电子体系、粒子数任意），在外部势和相互作用互作用CoulombCoulomb势作用下，势作用下，HamiltonianHamiltonian为为22r rHTVUTrr drVv rrr drUrrrr drdr 12112( )( )( )( ) ( )( )( ) ( ) ( )Har

18、tree单位外部势)( )()( rrrn电子密度算符电子密度算符电子密度分布电子密度分布n(r)是是 的期待值：的期待值：)( ,()(rnrn)( rn(1)(2)(3）(4)(5)(6)( rn（即 ）HKHK定理的证明定理的证明：外部势外部势v v( (r r) )是是n n( (r r) )的唯一泛函。即由的唯一泛函。即由n n( (r r) )唯一决唯一决定定。换句话说，如果有另一个。换句话说，如果有另一个vv( (r r) )，则不可能产生同样的，则不可能产生同样的n n( (r r).). 反证法：设有另一个反证法：设有另一个vv( (r r) ) ，其基态，其基态也会产生相同

19、的也会产生相同的n n( (r r).). v v( (r r) )vv( (r r) ) ， （除非（除非vv( (r r) )- -v v ( (r r)=const)=const）. . 与与 满足不同的满足不同的SchrSchr dinger dinger 方程：方程： H H = = E E HH = = E E 利用利用基态能量最小原理基态能量最小原理，有，有23UVTHVVHUVTH(, )(, )(,()(,)(,() ( )( ) ( )EHHHVVHVVEv rv r n r dr (7)(8)(9)HohenbergHohenberg-Kohn-Kohn定理的证明定理的证

20、明drrnrvrvEE)()()(即即同时，把带撇的与不带撇的交换得同时，把带撇的与不带撇的交换得drrnrvrvEE)()()(或者或者drrnrvrvEE)()()(10)(11)可见可见(10)与与(11)相互矛盾。表明相互矛盾。表明v(r) 不可能产生同样的不可能产生同样的n(r) .所以所以v(r) 是是n(r) 的唯一泛函。由于的唯一泛函。由于v(r) 决定整个决定整个H, 即系统的即系统的基态能量是基态能量是n(r) 的唯一泛函。的唯一泛函。 同理，同理，T和和U也是也是n(r) 的唯一泛函。可定义：的唯一泛函。可定义：)( ,()(UTrnF(12)式式(12)是一个普适函数，

21、适于任何粒子系和任何外部势。于是一个普适函数，适于任何粒子系和任何外部势。于是整个系统的基态能量泛函可写为：是整个系统的基态能量泛函可写为：)()()()(rnFdrrnrvrnE(13)定理定理2 2：如果：如果n(r) n(r) 是体系正确的密度分布，则是体系正确的密度分布，则En(r)En(r)是最低的是最低的能能 量，即体系的基态能量。量，即体系的基态能量。证明：设有另一个证明：设有另一个n n(r) ,(r) ,粒子数与粒子数与n(r) n(r) 相同为相同为N. N. 则则 实际实际计算是利用计算是利用能量变分原理能量变分原理，使系统能量达到最低（有一，使系统能量达到最低（有一定精

22、度要求）。由此求出体系的真正电荷密度定精度要求）。由此求出体系的真正电荷密度n n( (r r) ,) ,进而计算进而计算体系的所有其它基态性质。如，能带结构，晶格参数，体模量体系的所有其它基态性质。如，能带结构，晶格参数，体模量等等。等等。25)()()()( ,(),()( ,(),()()()()(rnErnErnEUTVUTVrnFdrrnrvrnE(14) 能量泛函公式能量泛函公式系统的基态能量泛函系统的基态能量泛函中，普适函数中，普适函数FnFn 可以把其包含的经典可以把其包含的经典CoulombCoulomb能部分写出：能部分写出：)()()()(rnFdrrnrvrnE（15）

23、rdrdnGnFrrrnrn)()(21)()()()(21nGrdrddrrnrvnErrrnrn其中其中Gn包括三部分：包括三部分：（16）（17）nEnEnTnGenergyselfxcsTsn=密度为密度为n(r) 的的非相互作用非相互作用电子体系的电子体系的动能动能。Excn=密度为密度为n(r) 的的相互作用相互作用电子体系的电子体系的交换关联能交换关联能。Eself-energyn=单个粒子的单个粒子的自能自能。应当扣除自能修正，下面。应当扣除自能修正，下面暂时忽略这一修正。暂时忽略这一修正。（18）局域密度近似局域密度近似(LDA)(LDA) HKHK定理已经建立了密度泛函理论

24、定理已经建立了密度泛函理论（DFTDFT）的框架，但在实际执行上）的框架，但在实际执行上遇到了严重困难。主要是遇到了严重困难。主要是相互作相互作用用电子体系的电子体系的交换关联能交换关联能E Excxcnn无法精确得到。无法精确得到。为了使为了使DFTDFT理论能理论能够付诸实施，够付诸实施，Kohn-ShamKohn-Sham提出了提出了局局域密度近似域密度近似(Local Density (Local Density Approximation, Approximation, LDALDA) )。本次只。本次只直接引用以便建立直接引用以便建立Kohn-ShamKohn-Sham方程。方程。

25、27局域密度近似（局域密度近似（LDALDA）LDA: LDA: 对于缓变的对于缓变的n(rn(r) ) 或或/ /和高电子密度情况，可采用如下近似：和高电子密度情况，可采用如下近似：28r)r () r (dnnnExcxc)r (nxc是是交换关联能密度交换关联能密度。它可以从均匀自由电子气的理。它可以从均匀自由电子气的理论结果得到。对于不同的论结果得到。对于不同的r, 有不同的有不同的n(r) .相应的有相应的有不同的不同的 。)r (nxc)r (nxc一种计算一种计算 的近似公式为（在的近似公式为（在Hartree单位下）：单位下）：0.4583341033211.4111230.0

26、333(1)ln(1)()( )ssxcrsnxrGGr axxxx rs是自由电子气的电子是自由电子气的电子”半径半径”。(19)(20)(21)利用利用LDALDA式式(19), (19), 能量泛函写为：能量泛函写为：drrnrndrdrdrrnrvnTnExcrrrnrns)( )( )( )( ) ( )( 21(22)上式考虑另一个电子密度上式考虑另一个电子密度n(r)。然后求。然后求En对对n的变分的变分En /n为最小。相当于改变为最小。相当于改变n(r) 使使En En。先求先求Tsn:为写出为写出Tsn，考虑，考虑v(r) 为一个试验的单电子势。可由为一个试验的单电子势。可

27、由v(r) 满足的单粒子方程，解出满足的单粒子方程，解出n(r) 。 21221( )( )( )( )( )iiiNiiv rrrn rr (23)(24)drrnrvnTdrrnrvnTrvrvNiissNiiiNiiiNiiiNii)( )( )( )( )( ,(),()( ( ,(11122112211(26)(25)于是能量泛函为于是能量泛函为 )( )()( )( )()(211nErdrddrrnrvdrrnrvnExcrrrnrnNii(27)求求 ，可得：，可得：0 nnE 由此得到：由此得到：31 ( ()()()( )( )0)xcEnn rnr rv rv rnnv

28、rv rdrn r drn r drconstrVconstdrrvconstdrrvrveffrVrVnnErrrrnnnErrrneffxcxcxc)()()()( )()()() ( ) ( 或或(28)(29).32由此得到由此得到Kohn-Sham方程：方程：( )21122( )( )( )( )( )( )( )( )( )( )()n rHreffiiiNiieffxcxcxcrrv rdrv rvEnVn rVrrVrrrrVrnrr i=Kohn-Sham本征值本征值称有效势称有效势经典经典Coulomb势势交换关联势交换关联势电子密度分布电子密度分布(30)Kohn-Sh

29、am方程是一个自洽方程组。先提供初始电子密度分布方程是一个自洽方程组。先提供初始电子密度分布n(r) , 它一般可由原子的它一般可由原子的nat(r) 叠加而成。依次求出经典叠加而成。依次求出经典Coulomb势、交换关联势、有效势。再求解势、交换关联势、有效势。再求解KS方程。再由方程。再由KS波函数构造新波函数构造新的电子密度分布。比较输入与输出的电子密度分布。如已自洽，的电子密度分布。比较输入与输出的电子密度分布。如已自洽，便计算总能，输出所有结果。便计算总能，输出所有结果。.33nin(r)n(r)=nat(r)求解、Vxc、Veff求解Kohn-Sham方程得到i由i构造nout(r

30、)比较nin与 nout(r)计算总能EtotNoYesnin与nout混合原子计算精度控制NoYes输出结果： Etot、 i、 n(r)Vxc、Veff、En(k)、N(E) 总能总能E Etottot表达式表达式.)()()()()()()()() ()(21;,211221drrndrrnrVrnEdrvrvrnrrExcxcirrrnrnixcmnmnRRZZHNiiitotmnmnHartree总能（不作详细推导，只了解物理意义）（不作详细推导，只了解物理意义）nxcnnExcxcxcnrV)(31)(32)第一项为动能，第二和第三项是总静电势能，最后一项是第一项为动能，第二和第三

31、项是总静电势能，最后一项是交换关联能。交换关联能。Zm是位于是位于Rm处的原子的核电荷。如果忽略处的原子的核电荷。如果忽略交换关联项，交换关联项，K-S方程的结果将与方程的结果将与Hartree近似一样。近似一样。虽然虽然K-SK-S方程十分简单，其计算量也只有方程十分简单，其计算量也只有HartreeHartree方程的水平，但却包含着深刻得多的物方程的水平，但却包含着深刻得多的物理内容。其中一个重要的概念性结果是，多体理内容。其中一个重要的概念性结果是，多体基态的解被准确地简化为基态密度分布之解，基态的解被准确地简化为基态密度分布之解，而这个密度是由单粒子的而这个密度是由单粒子的SchrS

32、chrdingerdinger方程给出方程给出的。的。由此，方程中的有效势在原理上包括了所有的由此，方程中的有效势在原理上包括了所有的相互作用效应，即相互作用效应，即HartreeHartree势、交换势（由势、交换势（由PauliPauli原理决定的相互作用所产生的势）和关联原理决定的相互作用所产生的势）和关联势（一个给定的电子对整个电荷分布的影响所势（一个给定的电子对整个电荷分布的影响所产生的势）。在这个意义上，它比产生的势）。在这个意义上，它比Hartree-Hartree-FockFock方程要优越得多。方程要优越得多。FormallyequivalentElectronInterac

33、tionExternal potentialHard problem to solveSchrdinger equation“Easy” problem To Solve DFTProperties of the systemNon-interacting electron(KS particle)Effective potentialxcxcEnn rnr rVrffreVv rrdr ( )( )( )()( xcVr ( )LDAGGAetc 量子力学体系的性质可以通过求解薛定格方程量子力学体系的性质可以通过求解薛定格方程(SE)进行计算进行计算(上图上图左边左边)。但更加容易的、形式上

34、等价的方法是求解。但更加容易的、形式上等价的方法是求解DFT的的KS方程方程(上图上图右边右边)。 但是准确的但是准确的 Excn(r) 并不知道。需要采用近似方法，如并不知道。需要采用近似方法，如 LDA or GGA。这就会影响。这就会影响 KS 解的精度。解的精度。SEDFT37LDFLDF近似下的近似下的电子电子相互作用电子电子相互作用示于图示于图C C，表明两种自旋的电子都有相同的交换，表明两种自旋的电子都有相同的交换关联空穴。如果进一步考虑不同自旋的关联空穴。如果进一步考虑不同自旋的电子有不同的分布，即所谓电子有不同的分布，即所谓局域自旋密局域自旋密度近似度近似（LSDLSD），则

35、不同自旋电子的交换），则不同自旋电子的交换关联空穴将有不同的形状，如图关联空穴将有不同的形状，如图d d所示。所示。电子电子相互作用图示电子电子相互作用图示38P(r)P(r)P(r)P(r)(a)(b)(c)(d)rrrr(a) Hartree (b) Hartree-Fock(c) DFT (d) SDFTP(r) =其余其余N-1个电子的几率分布个电子的几率分布r = 与固定电子的距离与固定电子的距离（固定电子（固定电子 r = 0）交换空穴交换空穴交换空穴交换空穴交换空穴交换空穴N-N-电子系统中电子电子相互作用电子系统中电子电子相互作用 上页给出了上页给出了N-N-电子系统中电子电子

36、相互作用的示电子系统中电子电子相互作用的示意图。考虑意图。考虑N N个电子中的一个电子（假定其自旋向上）个电子中的一个电子（假定其自旋向上）位于位于r r = 0 = 0处，横坐标表示与这一固定电子的距离，处，横坐标表示与这一固定电子的距离，纵坐标是其余纵坐标是其余N-1N-1个电子的几率分布个电子的几率分布p(rp(r) )。 a a）表示在）表示在HartreeHartree近似下，所有的电子都是独立的。近似下，所有的电子都是独立的。不管不管N-1N-1个电子的自旋是向上（实线）或向下（虚个电子的自旋是向上（实线）或向下（虚线），线），p(rp(r) )是均匀的并等于是均匀的并等于1 1，

37、没有结构，没有结构；39P(r)(a)固定电子固定电子r = 0r b b）说明在）说明在Hartree-FockHartree-Fock近似下，反对称的多电子波函近似下，反对称的多电子波函数反映了数反映了PauliPauli不相容原理，在不相容原理，在r r = 0= 0的固定电子周围的固定电子周围可以看到交换空穴，即自旋向上的电子被排斥，电子可以看到交换空穴，即自旋向上的电子被排斥，电子密度（实线）减少。但自旋相反的电子密度（虚线）密度（实线）减少。但自旋相反的电子密度（虚线）不受影响，也就是说，这些电子间的关联效应被忽略不受影响，也就是说，这些电子间的关联效应被忽略了了。事实上，事实上，

38、Hartree-FockHartree-Fock近似存在着一个严重的缺陷，近似存在着一个严重的缺陷，用它处理金属的电子结构时，用它处理金属的电子结构时，FermiFermi能级处的电子态密能级处的电子态密度为度为0 0，而且在实际计算上是如此的复杂，以至于很少，而且在实际计算上是如此的复杂，以至于很少有成功的计算结果。有成功的计算结果。40P(r)(b)交换空穴交换空穴r c c）LDFLDF近似下的电子电子相互作用，近似下的电子电子相互作用，表明两种自旋的电子都有相同的交换关联表明两种自旋的电子都有相同的交换关联空穴。空穴。 d d）如果进一步考虑不同自旋的电子有不）如果进一步考虑不同自旋的

39、电子有不同的分布，即所谓局域自旋密度近似同的分布，即所谓局域自旋密度近似（LSDLSD），则不同自旋电子的交换空穴将），则不同自旋电子的交换空穴将有不同的形状，如有不同的形状，如d d）所示。）所示。41P(r)P(r)(c)(d)r交换空穴交换空穴交换空穴交换空穴r电子电子相互作用图示电子电子相互作用图示42P(r)P(r)P(r)P(r)(a)(b)(c)(d)rrrr(a) Hartree (b) Hartree-Fock(c) DFT (d) SDFTP(r) =其余其余N-1个电子的几率分布个电子的几率分布r = 与固定电子的距离与固定电子的距离固定电子固定电子r = 0交换空穴交换

40、空穴交换空穴交换空穴交换空穴交换空穴SiSi中的对关联函数中的对关联函数g g43The pair correlation function g in the (110) plane, with one electron at the bond center. The atoms and bonds are schematically represented for bond chains along the 111 direction.(a), (b), and (c) show g with electron position r fixed on the bond center and

41、r ranging over the (110) plane, for parallel and antiparallel spins in VMC, and the spin averaged form in the LDA, respectively. The largest features are confined mainly to the bonding region where the first electron is located.parallel spin (VMC) opposite spin (VMC) spin averaged (LDA) SiSi中的交换关联空穴

42、中的交换关联空穴44 (a) spin-averaged pair correlation function (VMC) (b) exchange-correlation hole (VMC) (c) exchange correlation hole (LDA) One electron fixed at the tetrahedral interstitial site in the (110) plane. The atoms and bonds are schematically represented for bond chains along the 111 direction.R