CASTEP-计算理论分析总结实例分析.doc

《CASTEP-计算理论分析总结实例分析.doc》由会员分享，可在线阅读，更多相关《CASTEP-计算理论分析总结实例分析.doc（28页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、-_CASTEPCASTEP 计算理论总结计算理论总结XBAPRSXBAPRSCASTEP 特点是适合于计算周期性结构，对于非周期性结构一般要将特定的部分作为周期性结构，建 立单位晶胞后方可进行计算。CASTEP 计算步骤可以概括为三步：首先建立周期性的目标物质的晶体；其 次对建立的结构进行优化，这包括体系电子能量的最小化和几何结构稳定化。最后是计算要求的性质， 如电子密度分布(Electron density distribution)，能带结构(Band structure)、状态密度分布 (Density of states)、声子能谱(Phonon spectrum)、声子状态密度分布

2、(DOS of phonon)，轨道群分 布(Orbital populations)以及光学性质(Optical properties)等。本文主要将就各个步骤中的计算原 理进行阐述，并结合作者对计算实践经验，在文章最后给出了几个计算事例，以备参考。 CASTEP 计算总体上是基于 DFT，但实现运算具体理论有： 离子实与价电子之间相互作用采用赝势来表示； 超晶胞的周期性边界条件； 平面波基组描述体系电子波函数； 广泛采用快速 fast Fourier transform (FFT) 对体系哈密顿量进行数值化计算； 体系电子自恰能量最小化采用迭带计算的方式； 采用最普遍使用的交换-相关泛函实

3、现 DFT 的计算，泛函含概了精确形式和屏蔽形式。 一，一，CASTEPCASTEP 中周期性结构计算优点中周期性结构计算优点 与 MS 中其他计算包不同，非周期性结构在 CASTEP 中不能进行计算。将晶面或非周期性结构置于一个有 限长度空间方盒中，按照周期性结构来处理，周期性空间方盒形状没有限制。之所以采用周期性结构原 因在于：依据 Bloch 定理，周期性结构中每个电子波函数可以表示为一个波函数与晶体周期部分乘积的 形式。他们可以用以晶体倒易点阵矢量为波矢一系列分离平面波函数来展开。这样每个电子波函数就是 平面波和，但最主要的是可以极大简化 Kohn-Sham 方程。这样动能是对角化的，

4、与各种势函数可以表示 为相应 Fourier 形式。2()()(),k GVGGVGGVGGCCionHxci i k GGGi k G采用周期性结构的另一个优点是可以方便计算出原子位移引起的整体能量的变化，在 CASTEP 中引入外 力或压强进行计算是很方便的，可以有效实施几何结构优化和分子动力学的模拟。平面波基组可以直接 达到有效的收敛。 计算采用超晶胞结构的一个缺点是对于某些有单点限缺陷结构建立模型时，体系中的单个缺陷将以 无限缺陷阵列形式出现，因此在建立人为缺陷时，它们之间的相互距离应该足够的远，避免缺陷之间相 互作用影响计算结果。在计算表面结构时，切片模型应当足够的薄，减小切片间的人

5、为相互作用。 CASTEP 中采用的交换-相关泛函有局域密度近似（LDA） （LDA） 、广义梯度近似（GGA）和非定域交换-相 关泛函。CASTEP 中提供的唯一定域泛函是 CA-PZ，Perdew and Zunger 将 Ceperley and Alder 数值化结果进 行了参数拟和。交换-相关泛函的定域表示形式是目前较为准确的一种描述。NameDescriptionReference PW91Perdew-Wang generalized-gradient approximation, PW91Perdew and Wang PBEPerdew-Burke-Ernzerhof fun

6、ctional, PBEPerdew et al.RPBERevised Perdew-Burke-Ernzerhof functional, RPBEHammer et al.-_采用梯度校正的非定域或广义梯度近似泛函与电子密度梯度和电子密度都有关，这样可以d dr同时提高能量和结构预测的准确性，但计算耗时。CASTEP 中提供的非定域泛函有三种：PBE 泛函与 PW91 泛函计算在本质上实际是相同的，但在电子密度变化迅速体系中 PBE 泛函实用性更好；RPBE 是 特别用来提高 DFT 描述金属表面吸附分子能量的泛函，White and Bird 描述了各种梯度校正泛函计算方 法，利用广义

7、梯度近似计算总能量使用平面波基组与定域泛函相比并不直接。包含梯度近似的交换-相关 泛函计算时对电子密度数据的精度要求较高，对计算机内存占用会增大。通过采用与平面波基组总能量 计算中分裂交换-相关能量采用一系列空间网格相一致的方法来定义交换-相关势。平面波基组（平面波基组（Plane wave basis set）Bloch 理论表明每个 k 点处电子波函数都可以展开成离散的平面波基组形式，理论上讲这种展开形 式包含的平面波数量是无限多的。然而相对于动能较大的情况，动能|k k+G G|2很小时平面波系数 Ck k+G G更重 要。调节平面波基组，其中包含的 平面波动能小于某个设定的截止能 量，

8、如图所示（球体半径与截止能 量平方根成比例）：总能量计算会因为平面波特定能量 截止而产生误差，通过增加体系能 量截止数值就可以减小误差幅度。 理论上截止能量必须提高到总能量 计算结果达到设定的精确度为止， 如果你在进行关于相稳定性的研究， 而需要对比每个相能量的绝对值时， 这是一种推荐计算方法。不过，同 一个结构在低的截止能量下收敛引 起的差别要小于总体能量本身。因 此可以选用合适的平面波基组对几 何结构进行优化或进行分子动力学 研究。以上的方法对 Brillouin 区 取样收敛测试同样成立。 有限平面波基组的校正有限平面波基组的校正采用平面波基组的一个问题是截止能量与基组数量的变化是间断的

9、，一般而言在点基组（k-point set）中不同点对应不同能量截止（cutoffs）时就会产生这种不连续性。此外，在截止能量不变时， 晶胞形状和尺寸的变化都会引起平面波基组的间断。通过采用更加致密的点基组就可以解决这个问题， 与特定平面波基组相关的加权性也会消除。然而即使在基组取样很致密的情况下，这个问题依然存在， 对其近似的解决方法就是引入一个校正因子（correction factor），利用某个状态基组计算使用了无 限数量的点与实际采用的数量之间的差别来确定。晶体结构在进行几何优化时如果基组不能真正的达 到绝对的收敛，有限基组的纠正就很重要。比如硅的规范保守赝势很“软”，在平面波基组截

10、止能量 是时就已经可以得到准确的计算结果了。但如果计算状态方程时使用上述截止能量（比如体 积与总能量和压强都有关系），能量最小时对应的体积与体系内压为零时对应体积是不同的。在提高截 止能量和增加点取样基础上重复对状态方程的计算，这两个体积之间的差别会越来越小。此外截止能-_量低时计算得到的 E-V 曲线呈现锯齿状，提高截止能量计算的曲线连续而平滑。E-V 曲线中出现锯齿状 的原因在于平面波基组在相同的截止能量时由于晶体点阵常数不同引起的平面波基组数量的间断。对总 能量进行有限基组的校正，使得我们可以在一个恒定数量基组状态下进行计算，即使采用了恒定的截止 能量这个更强制条件也可以纠正计算结果。M

11、ilman 等详细的讨论了这种计算方法的细节。进行这种校正 所需要的唯一的参数就是 dEtot/d lnEcut，Etot是体系总能量，Ecut是截止能量。dEtot/d lnEcut的值很好 的表示了能量截止和点取样计算收敛性质。当它的数值（每个原子）小于 0.01 eV/atom 时，计算就 达到了良好的收敛精度，对于大多数计算 0.1 eV/atom 就足够了。非定域交换非定域交换-相关泛函相关泛函基于 LDA 或 GGA 的泛函的 Kohn-Sham 方程在计算能带带隙上存在低估。这对晶体或分子相关性质 以及能量的描述是没有影响的。然而要理解半导体和绝缘体性质，就必须得到关于电子能带结

12、构的准确 的描述。DFT 能带带隙计算误差可以通过引入经验“剪刀” 校正，相对于价带而言导带产生了一个刚 性的变化。当实验提供的能带带隙准确时，光学性质计算得到了较为准确的结果。电子结构实验数据缺 乏时采用“剪刀”工具进行预测性研究或对能带带隙调整是不可靠的。关于 DFT 计算中能带带隙问题已 经发展许多技术，但这些技术大多复杂而且很耗时，实际计算中最常用的是屏蔽交换（Sx-LDA），建 立在广义 Kohn-Sham 方法基础上。广义 Kohn-Sham 泛函允许我们将总能量交换分布泛函分离为非定域、 定域以及屏蔽密度组元。在 CASTEP 计算中采用的广义 Kohn-Sham 方法有：HF:

13、 exact exchange, no correlation HF-LDA: exact exchange, LDA correlation sX: screened exchange, no correlation sX-LDA: screened exchange, LDA correlation 与 LDA 和 GGA 相比 No local functionals 也有一些缺陷。在屏蔽交换泛函中不存在已知形式应力张量表 达方式，因此没有完全的非定域势可以用于单位晶胞结构优化或进行 NPT/NPH 动力学。这样利用这些 泛函计算的光学性质很有可能是不准确的。在哈密顿量中引入一个完全非定

14、域组元就可以解决这个问题， 这个额外的矩阵元破坏了光学矩阵元素由位置算符转换为动量算符常用表达形式，使得哈密顿量对易很 复杂。规范保守赝势和超软赝势规范保守赝势和超软赝势赝势是利用平面波基组计算体系总能量中关键的一个概念，价电子与离子实之间强烈的库仑势用全 势表示时由于力的长程作用很难准确的用少量的 Fourier 变换组元表示。解决这个问题的另一种方法从 体系电子的波函数入手，我们将固体看作价电子和离子实的集合体。离子实部分由原子核和紧密结合的 芯电子组成。价电子波函数与离子实波函数满足正交化条件，全电子 DFT 理论处理价电子和芯电子时采 取等同对待，而在赝势中离子芯电子是被冻结的，因此采

15、用赝势计算固体或分子性质时认为芯电子是不 参与化学成键的，在体系结构进行调整时也不涉及到离子的芯电子。为了满足正交化条件全电子波函数 中的价电子波函数在芯区剧烈的振荡，这样的波函数很难采用一个合适的波矢来表达。在赝势近似中芯 电子和强烈库仑势替代为一个较弱的赝势作用于一系列赝波函数。赝势可以用少量的 Fourier 变换系数 来表示。理想的赝势在芯电子区域是没有驻点的，因此需要平面波矢数量很少。众所周知的是现在将赝 势与平面波矢相结合对描述化学键是很有用的。全离子势的散射性质可以通过构筑赝势得到重现，价电 子波函数相位变化与芯电子角动量成分有关，因此赝势的散射性质就与轨道角动量是相关的。赝势最

16、普 遍表达方式是：-_VNL = |lm Vl are the spherical harmonics and Vl is the Pseudopotential for angular momentum l.在不同角动量通道均采用同一个赝势值称为定域赝势（Local Pseudopotential），定域赝势计算 效率更高，一些元素采用定域赝势就可以达到准确描述。赝势的硬度（hardness）在赝势的应用中是一 个重要的概念，当一个赝势可以用很少的 Fourier 变换组元就可以准确描述时称为“软赝势”，硬赝势 与此相反。早期发展的准确规范保守赝势很快就发现在过渡元素和第一周期元素（C、N、

17、O，等）中的 描述十分“硬”，提高规范保守赝势收敛性质的各种方法都已经被提出，在 CASTEP 中采用了由 Lin 等 提出的动能优化而来的规范保守赝势。 Vanderbilt 提出了另一种更基本的方法，放宽规范保守赝势的 要求，从而生成更软的赝势。在超软赝势方法中，芯电子区的赝平面波函数可以尽可能的“软”，这样 截止能量就可以大幅度的减少。超软赝势与规范保守赝势相比除了“更软”以外还有其它的优点，在一 系列预先设定的能量范围内遗传算法确保了良好的散射性质，从而使赝势获得更好变换性和准确性。超 软赝势通常将外部芯区按照价层处理，每个角动量通道中的占据态都包含了复合矢。这样就增加了赝势 的变换性

18、和准确性，但同时是以消耗计算效率为代价的。可转移性是赝势的主要优点。赝势是通过孤立 的原子或离子特定的电子排部状态下构建的，因此可以准确的描述原子在那些特定排部下芯区的散射性 质。在相应条件下产生的赝势可以用于各种原子电子排部状态以及各种各样的固体中，同样也确保了在 不同的能量范围内具有正确的散射状态。Milman 给出了不同化学环境和一系列结构中采用赝势描述准确 性事例。非定域赝势即使在最有效离散表示情况下，体系能量赝势计算依然占用了大量计算时间。此外， 在倒易点阵空间采用非定域赝势会因原子数目增多而耗时以原子立方数增大，因此对于大体系是很适用 的。赝势非定域性是指只有在超过原子芯区时它才会

19、扩展，由于芯区是很小的，特别是当体系包含有许 多的真空腔体时，在实空间采用赝势来计算就有很大的优势。这时计算量随体系中原子数目平方增长， 因此是很适合大体系计算的。将电子划分为芯电子和价电子在处理交换-相关相互作用时会产生新问题， 在原子芯区两个亚体系叠加在赝势产生过程中很难完全去屏蔽。在赝势能量算符中与电子密度存在非线 形关系的项就是交换-相关能。Louie 等采用了一种简单的方法来处理芯电子和价电子密度之间非线性的 交换-相关能。这种方法在很大程度上提高了赝势的可变换性，特别是自旋极化的计算更为准确。当准 芯区电子不能简单处理为价电子时非线性核校正就很重要。另一方面将他们简单地包含在价层亚

20、体系中 从本质上可以避免 NLCC 处理的必要性。规范保守赝势规范保守赝势：采用赝势计算关键在于可以有效的对化学键的价电子进行可再现的近似，赝势与全势在超过离子实半径 以后具有完全相同的函数形式。Figure 1. Schematic representation of the all-electron and pseudized wave functions and potentials两个函数平方幅度的积分数值应该 是相同的，这等同于要求赝势波函数具 有规范-保守性，比如每个赝波函数只 能描述一个电子的行为。这样的条件就 确保了赝势可以再现正确的散射 (Scattering Propert

21、ies)性质。 -_生成赝势的典型方法如下所述：选择某个特定的电子排部状态（不一定就是基态）全部电子计算在 一个孤立的原子中进行。从而得到原子价电子能量本征值和价电子波函数。选择一个离子赝势或赝波函 数参数形式，通过对参数的调节，使得赝原子计算和全电子原子赝势计算采用相同的交换-相关势，在 超过截止半径后与价电子波函数形式相同，赝势的本征值等于价电子的本征值。如果电子波函数和赝势 波函数满足正交归一，两者在截止半径以外的匹配性决定了规范-保守条件自动成立。离子赝势的截止 半径是实际物理芯区的二到三倍。截止半径越小，赝势越“硬”而适用性（transferability）好。计算精 度和效率决定了

22、实际中采用的截止半径的大小。 规范规范- -保守赝势优化保守赝势优化在固体计算中依据能量的截止存在一系列优化赝势的方法，Lin 基于 Rappe 早期工作提出了下列赝势产 生方法：在截止半径（cutoff radius）内，赝势波函数可以表达为： 4()()( )(),()()1PSjq RRli clcra jq ri liljq RRici clc是球形 Bessel 函数，在 r=0 和 r=Rc 之间有（i-1）个零点。为保证赝势的实用性，截止半径越大()jq rli越好。超过截止矢量 qc对动能最小化可得到系数。i2,22(,)( )( )( ) 00qcPSPSPSEqqdrrrd

23、qqqilllkc 在第一个方程中让qc等于q4。其他的三个限制条件使得赝波函数在进行 Lagrange 连乘（Lagrange multipliers）时保持正交化（normalization） ，并且使赝波函数在 Rc 处的第一个二介偏微分是连续的。半 径相关 Kohn-Sham 方程反转标准步骤产生的一个具有理想收敛性质的平滑赝势函数。Lee 提出了进一 步改进的方法，在 CASTEP 数据库中固体规范保守赝势就是采用他的思想设计的。这种通用的方法消除 了在特定的截止半径处赝波函数的二介偏微分必须是连续的条件，因为它是自动满足这个条件的。这 样对于特定截止半径 Rc 允许我们通过调节 q

24、c提高赝势的精度和计算效率。超软赝势（超软赝势（ULTRASOFT PSEDUPOTENTIAL）为了能够使平面波基组计算中所采用的截止能量尽可能的小，Vanderbilt 提出了超软赝势方法。众 所周知规范-保守赝势在收敛优化中存在本身缺陷，所以就设计了另一种方法。超软赝势基础是在大多 数情况下只有当紧密结合原子价轨道加权性分数大部分在芯区时，利用平面波基组计算才要求较高的截 止能量。在这种情况下，减少平面波基组的唯一方法就是解除（violate）规范-保守赝势成立条件，将这 些轨道中的电子从芯区移去。芯区的赝势就可以尽可能的“软”，从而使截止能量降低达到要求。从技 术上讲，通过引入一个广义

25、的正交归一化条件就可以完成。为了覆盖全部电子电荷，在芯区对由电子波-_函数模平方产生的电子密度进行适度放大（augmented） 。电子密度划分成两部分：扩展在整个晶体中 “软”部分和定域在芯区的“硬”部分。固体中超软赝势公式固体中超软赝势公式超软赝势中总能量与采用其他赝势平面波方法时相同，非定域势 VNL表达如下：(0) ,IIVDnmnmNlnm I投影算符 和系数 D（0）分别表征赝势和原子种类的差别，指数 I 对应于一个原子位置 。总能量用电 子密度可以表示为：2( )( )( ) ,IIIn rrQrinminmiinm I 是波函数， Q(r r) 是严格位于芯区的附加函数(Aug

26、ment function) 。超软赝势完全由定域部分， Vlocion(r r) 和系数 D(0), Q and 确定，这些变量计算方法在下文中将做介绍。 引入一个广义正交归一条件来解除规范-保守赝势的限制条件：S 是哈密顿重叠算符（Hermitian overlap operator） Sijij1,IISqnmnmnm I系数 q 是通过对 Q（r）积分得到，超软赝势的 Kohn-Sham 方程可以写为：HSiiiH 代表了动能和定域势能之和，如下所示：-_,IIIHTVDnmnmeffnm I在 Veff中包含离子定域势 Vlocion(r r)，Hartree 势和交换-相关势等项。

27、通过定义一些新参数就可以将因附 加（augmented）电子密度而产生所有项全部包含在赝势的非定域部分。(0) ( )( )II DDdrVr Qreffnmnmnm与规范-保守赝势对比，不同之处在于在超软赝势中存在重叠算符 S，波函数与 D 有关而且事实上投影算 符函数 （projector function）数量要比规范-保守赝势中大两倍多。与附加（augmented）电荷相关的一 系列计算可以在实空间(real space)中进行，这与函数中定域势的性质有关。多余的步骤不会对计算效 率产生较大的影响。在 Laasonen 文献中提供了超软赝势计算的详细方法以及总能量微分表达式。赝势生成：

28、与规范-保守赝势情况一样，在自由原子上对所有的电子进行计算，得到屏蔽原子势 VAE(r) （screened atomic potential） 。每个角动量选择一系列的参考能量l，一般两个能量参考点就足够了。这 些能量参考范围必须包含良好散射性质，在每个参考能量处求解与半径相关的 Kohn-Sham 方程，得到规 则初始点。选择截止半径，对上面产生的每个全电子波函数构筑一个赝势，唯一的限制条件是它必须 在 Rcl处与平滑相交。定义一个比所有芯区半径稍大的辅助半径 R。最后就形成了定域轨道（超过 R 时 就消失）：()TVnnnloc以及它们矩阵内积（inner products）： Bnmn

29、m这样就可以定义用于固体计算的变量 (Vlocion(r r), D(0), Q and ): ( )( )( )( )( );( )QrrrrrqdrQrnmnmnmnmnm -_1()Bnnmmm采用去屏蔽（descreening procedure）方法计算 Vlocion(r r), D(0)系数: ()()()()(0)()()ionVrVrVrVrloclocHxcDBqdrVr n rnmnmm nmloc在去屏蔽方法中可以引入非线性核校正方法（The nonlinear core correction (NLCC)） ，这与规范-保守赝 势中所采用的方法完全一致。在以下情况下超

30、软赝势是很适用的：赝本征值与所有电子本征值相同，在芯半径截止区以外赝轨道波函数与价电子波函数匹配一致；对于每赝本征值与所有电子本征值相同，在芯半径截止区以外赝轨道波函数与价电子波函数匹配一致；对于每 个参考能量散射性质都是正确的，这样通过增加参考能量点数目就可以系统的提高赝势的适用性；在参个参考能量散射性质都是正确的，这样通过增加参考能量点数目就可以系统的提高赝势的适用性；在参 考电子排部情况下，赝势价电子密度与全价电子密度相同。考电子排部情况下，赝势价电子密度与全价电子密度相同。关于非线性核校正关于非线性核校正Louie 等人第一次提出了非线性芯校正，使得赝势对磁系统的描述更准确。然而，对于

31、非自旋极化 体系中准芯区电子，NLCC 也具有同样的作用。DFT 总能量准确表达需要 NLCC，如下： totioneexcETEEE在赝势的计算式中，电子密度分别来自于芯区电子和价电子。将芯区能量假设为一常数并切不计入计算。 用一个价电子密度和由赝势计算得到的离子定域势 Eion来代替总电子密度，这样芯区电子与价电子之间 所有的相互作用全部转移到赝势上。由此可以推断电子密度线性化只是对动能和简单非线性交换-相关 能的一个近似，很明显当芯区电子和价电子在空间很好分离时是一个良好的近似。但如果两个区域电子 密度的叠加密切时，计算体系本身就会产生错误，进一步减弱赝势实用性。解决 NLCC 问题的方

32、法就是 调节赝势生成方法以及在固体中计算方法。在产生赝势时每个角动量通道对应一个屏蔽势，并且满足一 定的条件，比如规范-保守，赝波函数本征值与全电子波函数本征值相同等。这些屏蔽势（screened potentials）对应的原子赝波函数（atomic pseudowavefunctions）仅表示价电子。从这些波函数可以 得到价层赝电子密度（pseudo charge density），通过对势的屏蔽得到“光秃”离子势（bare）：( )( )( )( )ll ioneexcVrVrVrVr由于交换-相关势泛函是电子密度的非线性函数，对自旋极化体系采用这种方法产生的离子势与价电子 排列有关。

33、Louie 等提出了将上面方程替换为如下表达： -_( )( )( )( )( )ll ioneexccVrVrVrVrr在屏蔽原子势中减去总交换-相关势。此外，在计算交换-相关势时芯区电荷必须加到价电子中去，这个 额外原子状态信息传递给 CASTEP，在所有计算中芯区电荷认为（deemed）是相同的，这种做法的一个缺 点是在利用赝势计算时芯区电荷很难准确的用 Fourier 网格表示。而且通常芯区电子密度比价电子密度 大，这很容易将与价电子密度有关的影响掩盖掉。以下部分将对部分芯区校正方程建立做介绍，该方法 充分的认识到价电子与芯电子密度重叠的区域才是我们感兴趣的。靠近原子核的芯电子密度不会

34、产生物 理结果，虽然有如上所述的一些问题。部分 NLCC 采用一个在特定半径以外与c一致的函数替代全芯电 子密度，在原子核周围这个函数起伏是平滑的。在 CASTEP 中对一些特定元素在赝势中采用的部分芯区 校正使用了数值化的芯区电子密度。在规范保守赝势中虽然有相关的内容，但在计算中并没有采用这个 方法。A A IntroductionIntroduction toto DFTDFT第一性原理（The first principle）计算也称为从头算起（ab-initial calculation） ，由于固体的许多基本 的物理性质是由其微观的电子结构决定的，因此通过求解多粒子系统的 Schod

35、inger 方程，来获取固体全 部的微观信息从而预测宏观的性质。利用这个思想建立的能量的哈密顿量非相对论形式可以表示如下：22212(, .,.,)(, .,.,)22123, 1, 2,123, 1, 2,2221 2()(,12 ,1,21 212 1 21 2rrrr rrrrrrr rrrimk eeeei nnnjjk eeeei nnnjemijnjz ez zeejjjrrk eerrrri ii jj jei ejei njrrnnjjjji i ,.,.,)3, 1, 2,(,.,.,)123, 1, 2,rr rrreei nnnjErrrr rrrk k eeeei nn

36、nj考虑到原子核与核外电子质量差别以及电子驰豫时间比原子核驰豫时间要小三个数量级，因此利用 Born-Oppenheimer 近似将原子核运动和电子的运动分离，从而将体系波函数划分为电子波函数和原子核波函数两个部分，分别用和表示：(, .,.,)(, .,) (.,)123, 1, 2,1231, 2,rrrrrrrrrrrrrrkeeeei nnnjeeeeinnnj能量的哈密顿量可以分解为如下的两个方程：-_2222()(, .,(, .,2123)123),1, 2 1 2ez eejrrrrErrrrimeeeeikeeeeeieii ii jrrrreiejeinj ii 22221

37、12() (.,)(.,)21,2,1,2,1 21212nz ezzejjjrrrErrrjnnnjknnnnjmji jjjrrnjeinjrrnnjjjj 第一性原理严格求解仅在氢分子中实现了，对于多粒子体系的计算几乎是不可能的。目前均采用不同的 近似方法来实现计算，主要方法有 Hartree-Fock 近似和 DFT 近似。在 Hartree-Fock 近似中体系的哈密顿 量表示如下：()EJKHFijijTotalij ii为第 i 个电子的 Hartree-Fock 的轨道能，是库仑积分，表示电子静电互斥能，为交换积分。HFiJijKij交换积分所代表的交换能指电子由于自旋平行而引

38、起的电子轨道库仑能量减少的部分。 密度泛函理论(Density Functional Theory)建立了将多电子体系化为单电子方程的理论基础，并且给出了 有效势计算方法，是目前研究多粒子体系性质的一种普遍使用的重要方法。 该理论认为对于处于外势场 V（r）中相互作用的多电子系统，电子密度分布函数 (r)是决定该系统基 态物理性质的基本变量。密度泛函理论中体系的能量泛函表示如下：()()()()ETUEtxc：Kinetic energy; :classical electrostatic energy; :exchange and correlation energy()T()U()Exc由

39、上表达式可见体系能量是电子密度的泛函，因此可以进一步将上式表达为：2( )()( )( )( )( )( )2errErV rr dr TrdrdrErtxcr r在上式中第一项为电子在外场中的势能，第二项为电子的动能，第三项为电子相互之间的库仑能，第四 项为交换相关能，最后一项形式是未知的。 系统的电子密度分布可以表达如下：2( )( )rrii-_利用上式可以将动能项表示为：2 ( )2nTriii表达为：()U 11 ( )( )( )( )()()121)221 211( )( ) ()112211 2( )( ) 2nNNz zzaUrrrrrrrriijijRrr riiaRRaa

40、jaaNNz zzaarrr Rrr rRRaaaaaV rer VrVNNN 形式确定有两种方法：局域密度近似（LDA,Local Density Approximation）和广义梯度近似()Exc（GGA, General Gradient Approximation） 。在局域密度近似（LDA）中采用了均匀电子气的分布函数推倒出了非均匀电子气的交换-相关能泛函，从而得到的具体形式。从近期计算结果相关报道来()Exc看采用局域密度近似（LDA）计算在绝缘体中会产生较大的误差，而且对带隙宽的半导体等得到不正 确的结果。采用局域密度近似（LDA）主要的缺陷现归纳如下： 1 对光学跃迁带隙预测

41、很差（一般是过低估计带隙宽度） 。这虽然对基态性质如电荷密度，总能量以及 力影响不大，但在导带状态计算中却是个大问题，如关于光学性质，运输性质等的计算。在诸如光 伏装置等领域的研究中，带隙就是个很重要的问题。采用“剪刀”(Scissors)工具在固体带隙计算中 很有用，但对我们未获得实验结果的物质，是不能采用这个方法的。 2 对类似于二氧化硅这样的电子气分布极不均匀体系，基本假设中关于电子密度分布在空间是缓慢变 化的条件是不满足的，这样的体系采用 LDA 处理就存在难题。 3 LDA 简单的认为计算体系是顺磁性(Paramagnetic)的，对于包含未配对(Unpaired)自旋体系采用局域

42、自旋密度近似（LSDA） （对自旋向上(spin up)和向下(spin down)的电子分别采用密度泛函计算）是很 有用的，比如费米能级(Fermi level)处半填充的系统。 4 最后一个很少关注的领域就是玻璃陶瓷工业，LDA 对弱的结合键（如偶极涨落）很难描述，氢键 (Hydrogen bond)在 LDA 中也无法获得准确的计算结果。 GGA 近似则改进了 L(S)DA,将相关交换能确定为电子密度极其梯度的函数，在 GGA 学派中以 Perdew 等 人认为交换相关能的泛函形式应该以一定的物理规律为基础，构造了著名的 PBE 泛函。 将电子密度分布函数带入体系能量电子密度泛函中，对泛

43、函变分求极小值，可以得到 Kohn-Sham 方程：-_()()FrEriii 2 ( )FVrKS ()()() ()ExcriVrVrd rKSrirr 交换-相关能可以按照下式计算：( )( )Er Erdrxcxc:number of particles; :exchange-correlation energy per particles in an uniform electron gas ; ( )Erxc:distribution function of electron density.( )r称为交换-相关势和，表示为：()Excr( ) ()Exc rxcr 在 Cast

44、ep 计算中采用了周期性边界条件，单电子的轨道波函数满足 Bloch 定理，采用平面波展开式有：( )( )iK R iirer周期性边界条件下的波函数扩展为一系列分离的平面波波矢，这些波矢与晶体的倒易点阵矢量相联系。,( )iG R ii GrC e2.2 晶体光学性质的计算基于以下原理：-_电磁波在真空以及某种材料介质中传播时差别可以用一个复数式的折射指数来表示：Nnik在真空中 N 为实数，而且其大小为 1；在其他介质中时若材料对于光是透明的则是一个纯实数，虚部 对应材料的吸收系数（Adsorption Coefficient） 。它们之间的关系方程 2 所示： 2 .2k c 吸收系数

45、表示的是电磁波通过单位厚度的材料时能量的衰减分数，通常可以用材料焦耳热的产生来衡量。反射系数（Reflection Coefficient）可以简单通过将垂直光束照射材料的表面引起222(1)1.322(1)1nkNRnkN在计算光学性质时一般先计算虚部的介电常数，其他的性质与介电常数之间建立关系。虚部介电常数计 算式由下方程确定：222 .512nkN这样折射指数的实部和虚部以及介电常数之间的关系可以写为：22 2.612nknk光导率（Optical conductivity）也是一个普遍用来描述材料光学性质的物理量。光导率的表达式为方程 7：.712i这个参数用来描述金属的光学性质，但在

46、 CASTEP 中将计算范围扩大到了绝缘体和半导体。计算过程的 主要的区别在于前者的光学谱中 IR 部分与内部能带之间的转变密切相关，而者则在计算内时并没有完 全考虑到这些因素。 从虚部介电常数可以进一步得到材料电子的能量损失函数(Energy Loss Function)，它描述了电子通过均 匀的电介质时能量的损失情况，计算式如下所示：-_1 Im().8() 在实验中我们可以测定的光学性质参数有吸收系数和反射系数。从理论上而言，得到这( ) ( )R些参数以后可以将方程 2、3、4 表示为复数的形式之后得到表达式 1 中的实数部和虚数部。但在实际情 况下由于入射光源的复杂性，而且晶体结构中

47、极化效应使得材料介电常数并非是各向同性的。此外材料 表面几何结构也不是理想的平滑表面。这些因素就限制了对其光学参数的预测。在 CASTEP 中提供的光 学性质的计算支持体系极化，但状态只能在同种自旋间相互转化。 晶体中声子和电子之间的相互作用可以用电子基态波函数中包含的含时微扰项来表示，声子电场扰动引 起了电子函数占据态和未占据态间的转变（磁场引起的效应要弱一个因数 V/C） ，这些激发态（激子）聚 集态称为等离波子。单独的态激发称为单粒子激子，这些激子对光谱产生的结果是导带和价带的状态密 度之间的连接可以通过选择合适的加权性矩阵元素来实现。在 CASTEP 虚部介电常数的计算按照方程 9 进行：222 (,),().920, ,eccqOu rEEEkkkkkcu