第一节第一性原理计算方法

第一性原理计算的理论方法

随着科技的发展，计算机性能也得到了飞速的提高，人们对物理理论的认识也更加的深入，利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手段。这主要是因为在许多情况下计算机模拟比实验更快、更省，还得意于计算机模拟可以预测一些当前实验水平难以达到的情况。然而在众多的模拟方法中，第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员的青睐，成为计算材料学的重要基础和核心计算。本章将介绍第一性原理计算的理论基础，研究方法和ABINIT软件包。

1.1第一性原理

第一性原理计算(简称从头计算，the abinitio calculation)，指从所要研究的材料的原子组分出发，运用量子力学及其它物理规律，通过自洽计算来确定指定材料的几何结构、电子结构、热力学性质和光学性质等材料物性的方法。基本思想是将多原子构成的实际体系理解成为只有电子和原子核组成的多粒子系统，运用量子力学等最基本的物理原理最大限度的对问题进行”非经验”处理。【1】第一性原理计算就只需要用到五个最基本的物理常量即(mo.e.h.c.kb)和元素周期表中各组分元素的电子结构，就可以合理地预测材料的许多物理性质。用第一性原理计算的晶胞大小和实验值相比误差只有几个百分点，其他性质也和实验结果比较吻合，体现了该理论的正确性。

第一性原理计算按照如下三个基本假设把问题简化：

1．利用Born-Oppenheimer绝热近似把包含原子核和电子的多粒子问题转化为多电子问题。

2．利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单电子方程。

3．利用自洽迭代法求解单电子方程得到系统基态和其他性质。 以下我将简单介绍这些第一性原理计算的理论基础和实现方法：绝热近似、密度泛函理论、局域密度近似(LDA)和广义梯度近似(GGA)、平面波及赝势方法、密度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINIT。

1．2量子力学与Born-Oppenheimer近似

固体是由原子核和核外的电子组成的，在原子核与电子之间，电子与电子之间，原子核与原子核之间都存在着相互作用。从物理学的角度来看，固体是一个多体的量子力学体系【2】，相应的体系哈密顿量可以写成如下形式：

H?(r,R)?EH?(r,R) （1-1）

其中r,R分别代表所有电子坐标的集合、所有原子核坐标的集合。在不计外场作用下，体系的哈密顿量日包括体系所有粒子(原子核和电子)的动能和粒子之间的相互作用能，即

H?He?HN?He?N （1-2）其中，以是电子部分的哈密顿量，形式为：

?221e2 （1-3） He(r)????ri???2m2ii,i?|r?r|iii?i?上式的前一项代表电子的动能，后一项表示电子．电子之间的库仑相互作用

能，m是电子的质量。

原子核部分的哈密顿量HN，可以写成：

?221HN(R)????Rj??VN(Rj?Rj?)2j,j?j2Mjj?j?

(1-4)

原子核与电子的相互作用项可以写成：

He?N(r,R)???Ve?N(ri?rj)i,j

(1-5)

对于这样一个多粒子体系要对其实际精确求解是非常困难的，因此对其进行简化和近似是非常的必要。考虑到电子的质量比原子核的质量小很多(约103个数量级)，相对来说，电子的运动速度比核的运动速度要快近千倍。当电子在做高速运动时，原子核只在平衡位置附近缓慢振动，电子能够绝热于原子核的运动。因此，可以将上面的多体问题分成两部分考虑：当考虑电子运动时，原子核要处在它们的瞬时位置上；当考虑原子核运动时，就不需要考虑不电子在空间的具体

分布。这就是波恩(M.Born)和奥本海默(J.E.Oppenheimer)提出的绝热近似，或称波恩．奥本海默近似【2】，即Born-Oppenheimer绝热近似。此时系统的哈密顿量简化为：

?221e2H????ri????Ve?N(ri?Rj)?2m2ii,i?|r?r|i,jiii?i?

(1-6)

1.3 Hartree-Fock轨道近似

利用Born-Oppenheimer绝热近似就容易把包含原子核和电子的多粒子问题转化为多电子问题。求解方程(1-6)的困难在于电子与电子之间的库伦相互作

用项。假设不考虑电子之间的相互作用，就容易得到相互独立的单电子近似哈密顿量。为了把多电子问题简化成单电子问题【3】，如果把其他电子对所考虑电子的瞬时作用平均化和球对称化，则

|?i?(ri?)|2Vi(ri)???dri?|ri??rr|i?(i??i)

(1-7)

这样就可以把多电子问题转变成单单子问题。这时，整个系统的波函数就是每个电子波函数?i(ri)连乘积。单电子波函数应该满足单电子的Hartree方程：

|?i?(ri?)|2?22Hi????V(ri)???dri?2me|ri??rr|i?(i??i)

(1-8)

其中V(r)是该电子所受到的核的作用势。Hartree方程描述了每个坐标r处单电子在核作用势和其它电子的平均势中的运动，E是单电子的能量，简化后就可以从假设的一组?i(ri)出发，求解波函数时引入自治场方法，则整个系统的能量可以写为：

E??|H????i(r)H?i(r)??Eiii

(1-9)

上式并没有考虑到波函数是电子交换反对称的，于是需要考虑尸口础不相容原理，即把波函数写成(斯莱特)Slater行列式。此时体系的总能要增加一个由电子交换引起的交换项，体系的总能可改写成：

?i?(ri)?i(ri,)?i?,(r,)?,(ri)1iiE??|H????dri?(ri)Hi?i(ri)???dridri.2i,i,|ri?ti,|i?i（1-10）

对应的单电子方程为：

|?i?(ri?)|2?i?(r)?(r)?22[???V(ri)]?i(ri)???dri??i(ri)???dri??i?ii??i?(ri)???ii??i?(ri)2m|ri??rr||ri?ri?|i?(i??i)i?(i??i)i?（1-11）

这就是Hartree-Fock方程【4】。

2．1密度泛函的理论基础

密度泛函理论（Density Functional Theoty,简称DFT）【5】是从量子力学的基本原理出发，考虑电子结构，用体系的粒子数密度函数替代电子波函数来描述体系的理论。也就是说，假定固体、原子、分子等系统的基态能量和物理性质可以用电子密度函数唯一的确定。密度泛函理论是由于考虑了电子相关作用的Thomas-Fermi模型【6、7】，并在Hobenberg以及Kohn等人的工作【8】后发展成的，在经过Kohn和Sham(沈吕九)改进得到的电子密度泛函理论中的单电子方程，即

本文来源：https://www.bwwdw.com/article/pxy5.html