1、第一性原理计算简介在物理学中,第一性原理计算或称从头计算是指,基于构建物理学的基础定理,不作任何假设,例如:经验模型和拟合参数 ,所进行的计算研究。特别地 ,在凝聚态物理中,指的是运用薛定愕方程在一定的近似情况下,但不包括拟合实验数据所得到的参数和模型,对物质的电子结构进行计算 r 从而得到所研究物质的性质的一种研究方法。近些年,随着计算机技术的飞速发展,其运算能力越来越强大,使得人们可以处理更庞大更繁杂的物质结构体系,同时也使得计算物理成为了现代物理学,尤其是在凝聚态物理领域的一个重要分支。众所周知,固体是由相对重且带正电的粒子原子核,以及相对轻且带负电的粒子电子聚集在一起构成的。如果有个原
2、子,需要处理的问题是包含有 N+ZN(Z 为原子核所含的质子的个数) 个粒子的电磁相互作用,是一个多体问题。另一方面,由于处理的是微观粒子的运动,所以需要运用量子力学来描述其基本的运动规律和相互作用。对于该系统,精确的多粒子哈密顿量可以写作: i2ii1RHMFuuuuuuuuj其中位于為处的原子核的质量为 M,.,位于巧处的电子的质量为 m 一第一项是原子核的动能算符,第二项是电子的动能算符。后三项分别是描述电子与原子核,单个电子与其它电子以及单个原子核与其它原子核之间的库伦相互作用。很显然,直接精确求解(1.64)式几乎是不可能的。为了在合理的近似条件下得到体系的本征值,需要作不同层次的近
3、似。1.3.1 波恩-奥本海默(Bom-Oppenheimer)近似由于原子核的质量远大于电子质量,所以,原子核的运动速度远小于电子。因此,可以将原子“冻结”在固定的位置,并假设电子在瞬时与原子核是平衡的。或者说,只有电子在这个多体问题中是考察对象,原子核仅仅被当作一个带正电的外源场,相对于电子云是外在独立的。该近似被称为波恩-奥本海默(Bom-Oppenheimer)近似。原来的多体问题被简化成在原子的静电势下,瓜个带负电的粒子的相互作用。波恩-奥本海默认为,原子核不再运动,其动能为零,因此,(1.64)式的第一项被消除,最后一项退化为常数。(1.64)式简化为只含有电子气的动能,电子与电子
4、之间的相互作用所产生的势能,以及电子在可看作外源的原子核的势中的势能。(1.64)式可重写为:H = f + V + V,值得注意的是,(1.65)式中的动能以及电子与电子间的相互作用只取决于所处理的是系统是多电子系统,而不是多质子系统中强的原子内部作用力,并不依赖于特定的多电子系统本身,例如,Br2 或者水分子,Cu 还是 Fe, bcc-Fe 还是 fcc-Fe,等等。因此,前两项是普适的,包含特定系统信息的部分均在第三项中。1.3.2 密度泛函理论(Density functional theory)在波恩-奥本海默近似后,该量子多体问题得到了极大的简化,但是,依然很难直接求解。存在许多
5、方法将方程(1.65)进一步近似变为易于处理的形式,历史上非常重要的是 Hartree-Fock 方法。该方法在处理原子以及分子时效果很好,因此在量子化学中被广泛使用。但对于处理固体问题,其精度不够高。本文中使用的是更为现代且可能更强大的方法:密度泛函理论。密度泛函理论的建立可以追溯到 1964 年 Hohenberg 和 Kohn7提出的两条定理。1.3.2.1 Hohenberg-Kohn 定理两条定理的原始表述如下:第一定理:多电子体系(原子,分子,固体)基态时的电荷密度 pOO 与外源的势之间存在着一一对应。直接的结果是导致任意可观察量的算符 d 的基态期望值与基态的电子密度存在唯一泛
6、函:第二定理:对于为哈密顿量片的情况 ,基态总能的泛函 M 的形式如下:p=叫 f + v力 +叫 q 平 (1.67)其中 Hohenberg-Kohn 密度泛函p对任何多电子体系都适用。Eyp在厂所对应的基态密度下达到最小值(等于基态的总能) 。以上定理思想的三个关键词为可逆性(一一对应,p,00 所决定。因此,陷入了一个自洽问题: 等式的解(1,从可以得到/?,。很可能 P,与猜测的 P。不同。然后继续从 P,构造/fn,求得广 2,如此循环下去,最后会收敛到密度 P,由它构造的方程再次求得 P,那么,该密度则为所要求解的值。1.3.2.3 交换关联泛函为了求解多电子问题,首先采用了波恩
7、-奥本海默近似,然后通过 Kohn-Sham 方法,将相互稱合的多体问题变成了求解虚拟系统中无相互作用单粒子在外源势中的运动,其粒子的密度与真实的多电子系统的密度一致,并且,可以通过自洽方法求得该密度。然而,在这一过程中,并不知道交换关联泛函的具体形式,因此,必须对其形式进行猜测,给出进一步的近似。其中,被广泛使用的交换关联泛函近似叫做局域密度近似(Local DensityApproximation, LDA)。该近似认为交换关联泛函的形式如下:ET=pr)s,cpr)dr(1.80)其中,自由电子气函数的数值形式是可以求得的。该近似认为特定密度对交换关联能的贡献可以被分成许多具有恒定密度的
8、小空间。每一个小空间对总的交换关联能的贡献等于同等体积下均匀电子气的贡献,且均匀电子气的密度总体上与原材料在该体积的密度一致。局域密度近似被认为适用于电荷密度变化缓慢的系统,事实上,计算发现对于其它真实系统也足够精确。在局域密度近似下,各个极小体积对交换关联能的贡献只取决于该体积内的局域密度,因此,自然地想到可以改进该近似,使得对交换关联能的贡献还取决于邻近体积。也就是说,可以考虑密度的梯度。如此得到的近似被称为广义梯度近似(Generalized Gradient Approximation, GGA)=I p7)slr (pf)F, (p(r),|Vp(F)|)c/r其中 Fa 是无量纲的
9、标量泛函,是均匀电子气的交换关联能。通过对以各阶梯度为参数进行多项式展开,将交换关联能与密度的梯度产生联系。因此,该近似更适用于处理电子密度不均匀的原子、分子体系。目前,已经构造出许多形式的 GGA 泛函,包括 B88 (Becke) 9,PW91 (Perdew 和 Wang)10和 PBE (Perdew, Burke 和 Enzerhof)“。1.3.3 方程求解按照密度泛函理论,多电子问题被最终近似为单电子问题,通过求解如下的 Kohn-Sham 方程即可得到体系的基态密度二士+ + r 士 (1-82)方程的求解,使用的数学方法是寻找一组基函数以及系数将上式中的展开波函数属于一个无穷
10、维的函数空间,因此,P 在理论上应该是无限的。实际中,会将波函数用有限个数的基函数展开。虽然,有限组的基函数永远不能确切描述但是,可以尽量找寻接近一组?的基函数。一般地,存在三种选择基函数的方法,分别是:平面波和格点方法(Plane waves and grids)。这种方法采用平面波作为基函数展开,是解微分方程的通用方法,比如薛定谔方程和泊松方程。平面波尤其适用于周期性的晶体,釆用该方法可以对所求解的问题提供很多直观的理解,以及一些利于实际计算的简单算法。另一方面,格点方法极为适合实空间中的有限系统,在许多科学和工程领域均有应用。现代电 4结构的计算算法用到平面波和包含快速傅立叶变化的格点算
11、法。局域轨道方法(Localized orbitals)。这种方法采用局域的原子轨道作为基底函数,在化学中被广泛应用。相对于平面波方法,其基函数的数目会极大减少,因此,特别适用于较大体系的计算模拟,比如大分子,以及纳米器件设计。缀加函数方法(Augmentedfunctions) 。这种方法将空间分为两个区域 ,一个为原子位置处的球形区域内,另一个则是各个球形之间的区域。球形区域之间的区域,波函数变化平缓可以用平面波或者变化平缓的函数表示,在球形区域内部波函数变化剧烈可以用径向函数乘以球谐函数表示,两种函数应在球形边界处连续。原问题转化为寻找满足边界条件的函数。假定基函数己经确定,则可以处理(1.82)式的本征值问题,将(1.83)式带入(1.82)式中,并且左乘|(丨=1,., 户),得到? ? ? ? ?通过对角化哈密顿量对应的矩阵,可以得到尸个本征值以及对应的 P 个用所选基函数表示的本征函数,从而完成对方程(1.82)的求解。
