第一性原理计算方法论文.docx

第一性原理计算的理论方法随着科技的发展，计算机性能也得到了飞速的提高，人们对物理理论的认识也更 加的深入，利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手 段这主要是因为在许多情况下计算机模拟比实验更快、更省，还得意于计算机模拟可 以预测一些当前实验水平难以达到的情况 然而在众多的模拟方法中，第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员 的青睐，成为计算材料学的重要基础和核心计算本章将介绍第一性原理计算的理论基 础，研究方法和ABINIT软件包1.1第一性原理第一性原理计算（简称从头计算，the abinitio calculation），指从所要研究的 材料的原子组分出发，运用量子力学及其它物理规律，通过自洽计算来确定指定材料 的几何结构、电子结构、热力学性质和光学性质等材料物性的方法基本思想是将多原 子构成的实际体系理解成为只有电子和原子核组成的多粒子系统，运用量子力学等最基 本的物理原理最大限度的对问题进行”非经验”处理第一性原理计算就只需要用到五个最基本的物理常量即 （mo.ehckb）和元素周期表中各组分元素的电子结构，就可以合理地预测材料的许多物理性质。

用第一 性原理计算的晶胞大小和实验值相比误差只有几个百分点， 其他性质也和实验结果比较吻合，体现了该理论的正确性第一性原理计算按照如下三个基本假设把问题简化：1 .利用Born-Oppenheimer绝热近似把包含原子核和电子的多粒子问题转化为 多电子问题2 .利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单 电子方程3.利用自洽迭代法求解单电子方程得到系统基态和其他性质以下我将简单介绍这些第一性原理计算的理论基础和实现方法：绝热近似、密度 泛函理论、局域密度近似（LDA）和广义梯度近似（GGA）平面波及鹰势方法、密 度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINITo1 . 2量子力学与Born-Oppenheimer近似固体是由原子核和核外的电子组成的， 在原子核与电子之间，电子与电子之间，原子核与原子核之间都存在着相互作用从物理学的角度来看，固体是一个多体 的量子力学体系，相应的体系哈密顿量可以写成如下形式：H - （r,R）二（r,R） （1-1）其中r,R分别代表所有电子坐标的集合、所有原子核坐标的集合在不计外场作 用下，体系的哈密顿量日包括体系所有粒子（原子核和电子）的动能和粒 子之间的相 互作用能，即H 二 HeHN He』 （1-2）其中，以是电子部分的哈密顿量，形式为：卉2 2 1,,(r) = , rH；) d需J+好厂j卅I j一打上式的前一项代表电子的动能 用能，n是电子的质量。

II后一项表示电子.电子之间的库仑相互作原子核部分的哈密顿量Hn，可以写成:(1-4)Hn(R) = 7 _"Mjj1二、v「Rj-可j原子核与电子的相互作用项可以写成:(1-5)H (r, R)— 〃V (r )e JN — e JN i - rji ,j对于这样一个多粒子体系要对其实际精确求解是非常困难的， 因此对其进 行简化和近似是非常的必要考虑到电子的质量比原子核的质量小很多(约103个数量 级)，相对来说，电子的运动速度比核的运动速度要快近千倍当电子在做高速运动时，原子核只在平衡位置附近缓慢振动， 电子能够绝热于原子核的运动因此，可以将上面的多体问题分成两部分考虑：当考虑电子运动时，原子核要 处在它们的瞬时位置上； 当考虑原子核运动时，就不需要考虑不电子在空间的具体分布这就是波恩(M.Born)和奥本海默(J.E.Oppenheimer)提出的绝热 近似，或称波恩.奥本海默近似，即Born-Oppe nheimer绝热近似此时系统的哈密 顿量简化为：(1-6)A2 2 1 e2，「― -…VejN(ri- Rj)2m 2 A |「i ； " - i，j1.3 Hartree-Fock 轨道近似利用Born-Oppenheimer绝热近似就容易把包含原子核和电子的多粒子问题转 化为多电子问题。

求解方程(1-6)的困难在于电子与电子之间的库伦相互 作用项假设不考虑电子之间的相互作用，就容易得到相互独立的单电子近似哈 密顿量为了把多电子问题简化成单电子问题，如果把其他电子对所考虑电子的瞬 时作用平均化和球对称化，则vi(rj 八 a 占 i(i 分 IT-rr 1这样就可以把多电子问题转变成单单子问题 (1-7)这时，整个系统的波函数就是每个电子波函数'■ \(ri)连乘积单电子波函数应该满足单电子的 Hartree方程:-2V(r>r" |nr-r Iri (i %其中V (r)是该电子所受到的核的作用势Hartree方程描述了每个坐标r处单电子在 核作用势和其它电子的平均势中的运动， E是单电子的能量，简化后就可以从假设的一组'■ v(ri)出发，求解波函数时引入自治场方法，则整个系统的能量可以写为：e_,HH)=瓦化i (r) H性(「二瓦巳 (1-9)匚一 i上式并没有考虑到波函数是电子交换反对称的，于是需要考虑尸□础不相容原理，即把 波函数写成(斯莱特)Slater行列式此时体系的总能要增加一个由电子交换引起的交换 项，体系的总能可改写成：E=(空|H|空)二瓦 蚀出百)比空心)一上无Jdrprj 〃田⑴忆下「叫⑴\ i、 2i,j 山 一t「|(1-10)对应的单电子方程为：[ 2 丫山川匕⑴)^二 idr.M(―可(r) _vdr'' 〃〃 '川10；)2m .(寸厂 | r. ■ -r』.(.-./ | r-r|(1-11)这就是Hartree-Fock方程。

2.1密度泛函的理论基础密度泛函理论(Dens i ty Fu ncti on al Theoty,简称DFT是从量子力学的基本原理出发，考虑电子结构，用体系的粒子数密度函数替代电子波函数来描述体系的理论 也就是说，假定固体、原子、分子等系统的基态能量和物理性质可以用电子密度函数唯一的确 定密度泛函理论是由于考虑了电子相关作用的Thomas-Ferm模型，并在Hobenberg以 及Kohr等人的工作后发展成的，在经过Kohr和Sham沈吕九)改进得到的电子密度泛函理 论中的单电子方程，即Koh n-Shan方程，最终才使密度泛函理论得到实际的应用密度泛函理论是研究多粒 子系统基态的重要方法之一， 它不但成功将多电子问题转化为简单的单电子方程理论，而且也成为计算分子、固体等的电子结构和总能的有效手段2. 2Thomas-Fermi-Dirac 近似在1927年，H.Thoma和E.Fermi就已经提出来建立在均匀电子气基础上的Thomas-Fer m模型在这个均匀的电子气模型中，电子不受外力，电子与电子之间也没有相互作用， 经过求解电子运动的波动方程和简单的推导， 就能看出，体系的能量仅与电子密度的函数有关。

在 1930年，Dirac考虑了电子的交换相互作用并推导出来在外势Vext (r)中的电子的能量泛函的表达式如下：r -rETF (n) Jd3rn (r) 53 + Jd3rVext (r) n (r) +C?Jd3rn (r) 43 2d'rd3r上式从左到右各项表达式分别表示： 动能的局域近似、外力能作用、交换关联相互作用、经典的经典作用能由于Thomas-Fermi-Dirac近似太粗略简单，没有考 虑到物理、化学中的一些本质现象而没用得到广泛的应用 f鲫2. 3 Hobenberg-Kohn 定理密度泛函理论的基本理论基础是Hobe nberg和Kohn提出的非均匀电子气理 论的第 一、第二定理第一定理：处于外势Vext(r)中的不计自旋的电子体系，不可能存在另外一个外势Vext(r)也 有相同的密度函数，即其外势Vext(r)可由电子密度唯一决定此时系统的哈密顿量H=T+V+U这里T表示电子动能，V是外势，U为电子相互作用势在 不同体系的哈密顿量H中，外势V是不一样的，而电子动能T和电子相 互作用势的表达式 是相同的因此只要外势确定，体系的哈密顿量H也就确定了根据公式H，，只要H是确 定的，系统的波函数也确定，也可以说电子密度决定了系统波函数的所有性质。

第二定理：对于已定的外势，体系基态能量能于基态能量泛函 E( n(r))的极小值对于不计自旋的全同电子体系，其能量泛函 E(n(r))可写为：2E(n (r))「V(r)n(r)dr T[n (r)]斗；CAdrdr E』n(r)](2-13)其中，第一项是电子在外势场中的势能，第二项表示无相互作用电子气的动能，第三项是电子间的库伦作用能，第四项是电子间的交换关联能第二定理的基本点是在 粒子数不变条件下求能量对密度函数的变分， 就可以得到体系基态的能量E(n)但是Hobenberg-Kohn定理中还存在一些不足之处：(1)电子密度分布函数n (r)的具体形式不明确⑵无相互作用电子气的动能泛函T[ n (r)]不知道3)电子间的交换关联能泛函Exc[n(r)]不清楚针对前两个问题可以用Kohn-Shan方程解决第三个问题，通常是采用各种近似得到 电子间的交换关联能2. 4有效单电子近似：Kohn-Shan方程1965年，Kohn和Shan提出了这样一个假设：体系的电荷密度可以用电子波函数构造此时电荷密度n(r)A |T(r)|2 (2-14)i 4这样前面遇到的问题就可以顺利解决将\(「)代到(2. 13)E[n(r)] =T°[n(r)] . n(rgjdrEh[n(r)] EJn(r)](2-15)其中，To[n(耽五泸乃切-2 NEh[ n(r)]1 n (r) n(r)drd r2 ' r -r(2-17)虽然Exc[n(r)]与电子密度n(r)之间的函数表达式不知道，但是Kohn和Sham 成功的将多电子体系的薛定谔方程问题简单的归结为单电子在周期性势场中 的 运动的单电子方程。

此时，只要求解在周期性势场方程：N个无相互作用的单电子[—厂、2 Vks[ n(r)]?i(r) = ; ih(r) (2-18)2m其中，论®n(r)]旦但旦凹] (219)6n(r) On(r) *根据Kohn-Sham勺本征值「，体系的总能量可写成：二 1 _n(r)n(r ) , 「E = ； i drdr - VXC[n(r)]n(r)dr EXC[n(r)]i2 r -r XC XC(2-20)需要注意的是Kohn-ShamT程中本征值没有实际的物理意义唯一的例外是体系 的最高占据轨道，它的本征值对应于体系的离子化能2. 5交换关联能近似电子间的交换关联能泛函 Exc[n(r)]表示的是所有其它多体项对总能的贡献它的物理意思是：当单电子在一个多电子体系运动中，由于考虑电子之间的库伦排 斥，电子与体系之间就有交互关联作用换句话说，就是在同一时刻两个电子不可能占据同一个位置，也就产生了交换关联能 Exc[n(r)]在HoBenerg-Kohn-Shant勺理论框架下，多电子体系基态的薛定谔方程问题转化成了有 效的单电子方程问题，这种形式的描述比胁舰 P粕出方程更严密更简洁但前提是要处理好交换关联能后这个理论才有实际的应用价值。

所以交换关联能泛函在密度泛函理论中占有非常重。