第一性原理

第二章 第一性原理计算方法与软件介绍

19世纪末，科学家们发现经典力学和经典电动力学在描述物质的微观系统时存在明显不足，对实验中的许多现象不能做出真正合理的解释。鉴于此，20世纪初物理学家们在旧量子论的基础上建立了量子力学，主要研究原子、分子、凝聚态物质等内部微观粒子的结构、运动规律等性质，目前已广泛应用于物理、化学、材料等学科领域。随着量子力学理论的不断完善，并结合日趋成熟的计算机技术，量子计算模拟成为了现代科学中必不可少的研究手段之一。第一性原理计算(First -principles calculation)，亦称为从头算(Ab -initio calculation)。该计算方法可根据量子力学基本原理，基于密度泛函理论对材料微观体系的状态和性质进行理论上的预测，且计算过程中不需要使用任何经验参数，只需要一些基本物理量(电子电荷质量e 、电子静止质量m 0、光速c 、普朗克常数h 、波尔兹曼常数k B )。本工作所选用的计算程序为Materials Studio 软件中的CASTEP 量子力学模块，该模块是基于密度泛函理论的从头算量子力学程序。本章节将简要的介绍密度泛函理论和CASTEP 计算模块。

2.1密度泛函理论概述

第一性原理主要的研究对象是多原子体系。它依据量子力学原理，且在无任何实验参数引入的情况下，将多原子体系当作由自由电子和原子核组成的多粒子体系进行处理。然而，关于量子力学中多粒子体系处理的出发点则为著名的薛定谔方程(Schrödinger Equation)。Schrödinger 方程是量子力学的一个基本方程，也是第一性原理计算方法的核心，它是由奥地利物理学家薛定谔(Schrödinger)于1926年提出的。该方程可用于描述微观粒子的运动规律，故亦被称为薛定谔波动方程(Schrödinger Wave Equation)，其定态方程描述如下：

2

2[()]()(,)2V r r,t i r t t ψψμ∂-∇+=∂

(2-1) 式中ħ为约化普朗克(Plank)常数；μ和V(r)分别表示粒子质量和势场；r 和t 则为体系中所有电子与原子核的位置坐标；Ψ(r,t)是系统波函数，即运动的微观粒子

在V(r)势场下的波函数。但Schrödinger 方程在描述真实的复杂系统时求解过程非常困难，只能处理氢原子等简单的电子体系。因此，研究者们采用一些方法对Schrödinger 方程进行了科学合理地近似和简化，其中最经典的方法是密度泛函理论(Density Functional Theory ，DFT)。DFT 主要是将多电子的问题转化为单电子进行处理，通过基态电子密度来描述体系的基态物理性质。

2.1.1早期Thomas -Fermi 模型

1927年，托马斯(Thomas)和费米(Fermi)建立了Thomas -Fermi(T -F)模型，该模型提出了一个用电子密度来描述体系动能的表达式。但是，T -F 模型认为电子彼此之间互不影响，忽略了电子-电子的交换相互作用，故其计算精度大大降低，并未得到广泛应用。1930年，Dirac 在T -F 模型上引入了描述电子之间的交换相互作用的局域近似，即Thomas -Fermi -Dirac(T -F -D)理论，该理论表达式描述如下： ')'()('21)()()()()(3334334531r r r n r n r rd d r rn d r n r rV d r rn d C n E ext TF -+++=⎰⎰⎰⎰

(2-2)

式中第一、二、三、四项分别为系统的电子动能、外场势能、电子-电子交换相互作用能和电子间静电作用能的表达式。T -F -D 理论可以通过最低能量的寻找，来获得基态电子密度。但T -F -D 理论模型仍存在较多的漏洞，并且过于简单化，对于分子、原子和凝聚态物质等微观性质的准确描述还需要进一步完善。

2.1.2 Hohenberg -Kohn 定理

对于密度泛函理论，其坚实的理论依据始于Hohenberg -Kohn(H -K)定理[92]。该定理是霍恩伯格(Hohenberg)和科恩(Kohn)于1964年在T -F 模型的基础上提出的，它可以总结归纳为两个定理：定理①证明了体系的基态能量仅仅是电子密度ρ(r)的泛函；定理②证明了以基态密度ρ(r)为变量，将体系能量最小化之后就能够得到基态能量。

根据H -K 定理，系统的基态能量泛函可表示为：

()][()]()()G E r =E r +V r r dr ρρρ[⎰ (2-3)

其中，第二项代表外场中电子的势能作用；而泛函E[ρ(r )]可表示为：

()()[()][()][()]xc 1r r E r =T r +drdr E r 2r r ρρρρρ'''+∣-∣

⎰⎰ (2-4) 式(2-4)中，第一、二、三项分别代表系统的动能、电子间库仑作用能和电子间交换-关联能。虽然H -K 定理证明了通过基态的电子密度分布函数能够得到体系的总能，但仍不能确定体系总能与基态电子密度分布函数之间所存在的具体泛函形式，并且对于如何利用泛函极值性质求解体系总能的问题尚未得到解决。

2.1.3 Kohn -Sham 方程

密度泛函理论真正意义上的广泛应用是通过Kohn -Sham(K -S)方程实现的。由于上述H -K 定理无法确定电子密度函数ρ(r )、动能泛函T [ρ(r )]和交换-关联能泛函E xc [ρ(r )]。因此，在1965年科恩(Kohn)和沈呂九(Lu -Jeu Sham)在电子气近似的基础上对H -K 定理做了一些创新性的改进，他们通过证明认为一个多粒子体系的粒子密度数可以由一个简单的单粒子波动方程获得，这对密度泛函理论的实际应用起到了关键作用，该方程即为Kohn -Sham 方程，详细描述如下：

体系的电子密度分布是其组成各电子的波函数的平方和，即：

()()N

2i i 1r r ρψ==||∑ (2-5)

而无相互作用体系的动能泛函T s [ρ(r )]则简单认为是N 个单电子的动能之和：

[()]=()()()N 2s i i i 21T r r r dr 2m *ρ-ψψ=∇∑ħ (2-6) 故可得到： [()]2KS i i i V r r r 2m

ρψEψ{-∇+}()=()ħ (2-7) 在DFT 中，体系能量是电子密度ρ(r )的泛函：

[)][)]())[)][)]ext Coul xc E r =r V r r dr r r ρρρE ρE ρ((+(+(+(⎰s T (2-8) 式中E Coul [ρ(r )]为库仑能，而库仑势为： [)]Coul r V r dr r r '''ρρ()

(=|-|⎰ (2-9) 对于H -K 定理中未确定的交换-关联能泛函E xc [ρ(r )]则用电子体系的交换-关联能密度εxc [ρ(r )]来代替。