第一性原理计算的理论方法随着科技的发展，计算机性能也得到了飞速的提高，人们对物理理论的认识也更加的深入，利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手段。
这主要是因为在许多情况下计算机模拟比实验更快、更省，还得意于计算机模拟可以预测一些当前实验水平难以达到的情况。
然而在众多的模拟方法中，第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员的青睐，成为计算材料学的重要基础和核心计算。
本章将介绍第一性原理计算的理论基础，研究方法和ABINIT软件包。
1.1 第一性原理第一性原理计算( 简称从头计算，the abinitio calculation)，指从所要研究的材料的原子组分出发，运用量子力学及其它物理规律，通过自洽计算来确定指定材料的几何结构、电子结构、热力学性质和光学性质等材料物性的方法。
基本思想是将多原子构成的实际体系理解成为只有电子和原子核组成的多粒子系统，运用量子力学等最基本的物理原理最大限度的对问题进行”非经验”处理。
【1】第一性原理计算就只需要用到五个最基本的物理常量即( m o.e.h.c.k b ) 和元素周期表中各组分元素的电子结构，就可以合理地预测材料的许多物理性质。
用第一性原理计算的晶胞大小和实验值相比误差只有几个百分点，其他性质也和实验结果比较吻合，体现了该理论的正确性。
第一性原理计算按照如下三个基本假设把问题简化：1．利用Born-Oppenheimer 绝热近似把包含原子核和电子的多粒子问题转化为多电子问题。
2．利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单电子方程。
3．利用自洽迭代法求解单电子方程得到系统基态和其他性质。
以下我将简单介绍这些第一性原理计算的理论基础和实现方法：绝热近似、密度泛函理论、局域密度近似(LDA)和广义梯度近似(GGA)、平面波及赝势方法、密度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINIT。
1．2量子力学与Born-Oppenheimer 近似固体是由原子核和核外的电子组成的，在原子核与电子之间，电子与电子之间，原子核与原子核之间都存在着相互作用。
从物理学的角度来看，固体是一个多体的量子力学体系【2】，相应的体系哈密顿量可以写成如下形式：H (r,R) E H(r ,R) (1-1)其中r,R 分别代表所有电子坐标的集合、所有原子核坐标的集合。
在不计外场作用下，体系的哈密顿量日包括体系所有粒子( 原子核和电子) 的动能和粒子之间的相互作用能，即H H e H N H e N (1-2) 其中，以是电子部分的哈密顿量，形式为：22 1 e2H e(r) r2i 1 e(1-3)i 2m 2i i,i i|r i r i |上式的前一项代表电子的动能，后一项表示电子．电子之间的库仑相互作用能，m是电子的质量。
原子核部分的哈密顿量H N ，可以写成：221H N (R) 2R j1V N(R j R j )N j2M j R j2 j,j N j jjj(1-4)原子核与电子的相互作用项可以写成：H e N (r,R) V e N (r i r j )i,j(1-5)对于这样一个多粒子体系要对其实际精确求解是非常困难的，因此对其进行简化和近似是非常的必要。
考虑到电子的质量比原子核的质量小很多( 约103个数量级) ，相对来说，电子的运动速度比核的运动速度要快近千倍。
当电子在做高速运动时，原子核只在平衡位置附近缓慢振动，电子能够绝热于原子核的运动。
因此，可以将上面的多体问题分成两部分考虑：当考虑电子运动时，原子核要处在它们的瞬时位置上；当考虑原子核运动时，就不需要考虑不电子在空间的具体分布。
这就是波恩(M.Born) 和奥本海默(J.E.Oppenheimer) 提出的绝热近似，或称波恩．奥本海默近似【2】，即Born-Oppenheimer 绝热近似。
此时系统的哈密顿量简化为：(1-6)1.3 Hartree-Fock 轨道近似利用 Born-Oppenheimer 绝热近似就容易把包含原子核和电子的多 粒子问题转化为多电子问题。
求解方程 (1-6) 的困难在于电子与电 子之间的库伦相互作 用项。
假设不考虑电子之间的相互作用，就容易得到相互独立的单电子近似哈密顿量。
为了把多电子问题简化成单电子问题【 3】， 如果把其他电子对所考虑电子的瞬时作用平均化和球对称化，则(1-7) 这样就可以把多电子问题转变成单单子问题。
这时，整个系统的波 函数就是每个电子波函数 i (r i ) 连乘积。
单电子波函数应该满足单电 子的 Hartree 方程：H i 2V(r i )d r i| i (ri )|2m ei (i i) i|r i r r |(1-8)其中V(r) 是该电子所受到的核的作用势。
Hartree 方程描述了每个坐 标r 处单电子在核作用势和其它电子的平均势中的运动， E 是单电子 的能量，简化后就可以从假设的一组 i (r i ) 出发，求解波函数时引入 自治场方法，则整个系统的能量可以写为： E |H i (r) H i (r) E iHr 2i 1i 2m2i i,ii|r i e2V e N (r i R j )r i | i,jV i (r i )i (i i)dr i| i (r i ) |2|r r |(1-9) 上式并没有考虑到波函数是电子交换反对称的，于是需要考虑尸口 础不相容原理，即把波函数写成 ( 斯莱特)Slater 行列式。
此时体系 的总能要增加一个由电子交换引起的交换项，体系的总能可改写 成：1 i (r i ) i (r i , ) i , (r i ,) i , (r i )E | Hdr i i (r i )H i i (r i ) dr i dr i . i i i i i i i ii 2 i,i ,i|r i t i , |(1-10)对应的单电子方程为：2 2 | (r ) |2(r ) (r )[ 2 V(r i )] i (r i ) d r ii ii (r i ) dr ii i i ii(r i ) ii i (r i ) 2m i (i i) i|r i r r |i (i i)|r i r i | i(1-11) 这就是 Hartree-Fock 方程【 4】。
2．1密度泛函的理论基础密度泛函理论( Density Functional Theoty, 简称 DFT )【 5】是从量子力学的基本原理出发，考虑电子结构，用体系的粒子数密度函 数替代电子波函数来描述体系的理论。
也就是说，假定固体、 原子、 分子等系统的基态能量和物理性质可以用电子密度函数唯一的确 定。
密度泛函理论是由于考虑了电子相关作用的 Thomas-Fermi 模型 【6、7】，并在 Hobenberg 以及 Kohn 等人的工作【 8】后发展成的， 在经过 Kohn 和Sham(沈吕九 ) 改进得到的电子密度泛函理论中的单 电子方程，即Kohn-Sham 方程【 9】，最终才使密度泛函理论得到实际的应用。
密度泛函理论是研究多粒子系统基态的重要方法之一，它不但成功将多电子问题转化为简单的单电子方程理论，而且也成为计算分子、固体等的电子结构和总能的有效手段。
2．2Thomas-Fermi-Dirac 近似在1927年，H.Thomas和E.Fermi 就已经提出来建立在均匀电子气基础上的Thomas-Fermi模型【6、7】。
在这个均匀的电子气模型中，电子不受外力，电子与电子之间也没有相互作用，经过求解电子运动的波动方程和简单的推导，就能看出，体系的能量仅与电子密度的函数有关。
在1930年，Dirac 考虑了电子的交换相互作用并推导出来在外势V ext (r)中的电子的能量泛函的表达式如下：E TF (n) C1 d 3rn(r)53d3rV ext(r)n(r) C2 d3rn(r)43 1d3rd3rn(r)n(r )TF 1 ext 22 | r r | (2-12)上式从左到右各项表达式分别表示：动能的局域近似、外力能作用、交换关联相互作用、经典的经典作用能。
由于Thomas-Fermi-Dirac 近似太粗略简单，没有考虑到物理、化学中的一些本质现象而没用得到广泛的应用f 鲫。
2． 3 Hobenberg-Kohn 定理密度泛函理论的基本理论基础是Hobenberg和Kohn提出的非均匀电子气理论的第一、第二定理。
第一定理：处于外势V ext ( r )中的不计自旋的电子体系，不可能存在另外一个外势V ext (r )也有相同的密度函数，即其外势V ext (r ) 可由电子密度唯一决定。
此时系统的哈密顿量H=T+V+U，这里T表示电子动能，V是外势，U为电子相互作用势。
在不同体系的哈密顿量H中，外势V是不一样的，而电子动能T和电子相互作用势U的表达式是相同的。
因此只要外势确定，体系的哈密顿量H也就确定了。
根据公式H E ，只要H是确定的，系统的波函数也确定，也可以说电子密度决定了系统波函数的所有性质。
第二定理：对于已定的外势，体系基态能量能于基态能量泛函E(n(r)) 的极小值。
对于不计自旋的全同电子体系，其能量泛函E(n(r)) 可写为：2E(n (r)) V(r)n (r)dr T[n (r)] e Cn (r)drdr E xc[n(r)]2 |r r |(2-13) 其中，第一项是电子在外势场中的势能，第二项表示无相互作用电子气的动能，第三项是电子间的库伦作用能，第四项是电子间的交换关联能。
第二定理的基本点是在粒子数不变条件下求能量对密度函数的变分，就可以得到体系基态的能量E(n) 。
但是Hobenberg-Kohn定理中还存在一些不足之处：(1)电子密度分布函数n (r) 的具体形式不明确。
(2)无相互作用电子气的动能泛函T[ n (r) ]不知道。
(3)电子间的交换关联能泛函E xc[n(r)] 不清楚。
针对前两个问题可以用Kohn-Sham方程解决。
第三个问题，通常是采用各种近似得到电子间的交换关联能。
·2． 4有效单电子近似： Kohn-Sham 方程1965年， Kohn 和Sham 提出了这样一个假设：体系的电荷密度可以用 电子波函数构造。
此时电荷密度N2n(r) | i (r)|2i1(2-14)这样前面遇到的问题就可以顺利解决。
将 i (r) 代到(2 ．13)变形成； E[n(r)] T o [n(r)] n(r)V ext (r)dr E h [n(r)] E xc [n(r)] ( 2-15 )2 N2其中T o [n(r)]i 2 i，2m e i 1( 2-16 )E [n(r)]1 n(r)n(r )drdrE h [n(r)]drdr2 r r(2-17)虽然 E xc [n(r)] 与电子密度 n(r) 之间的函数表达式不知道，但是Kohn 和Sham 成功的将多电子体系的薛定谔方程问题简单的归结为单电 子在周期性势场中的运动的单电子方程。