CASTE计算理论总结
XBAPRS
CASTEP寺点是适合于计算周期性结构，对于非周期性结构一般要将特定的部分作为周期性结构，建
立单位晶胞后方可进行计算。
CASTEF十算步骤可以概括为三步：首先建立周期性的目标物质的晶体;
次对建立的结构进行优化，这包括体系电子能量的最小化和几何结构稳定化。最后是计算要求的性质, 如电子密度分布(Electron
density distribution) ,能带结构(Band structure)、状态密度分布
of states)、声子能谱(P ho non sp ectrum)、声子状态密度分布 (DOS of p ho non)，轨道群分布 populations) 以及光学性质(Optical properties)
等。本文主要将就各个步骤中的计算原理进行阐述,
并结合作者对计算实践经验，在文章最后给出了几个计算事例，以备参考。
CASTEP十算总体上是基于 DFT但实现运算具体理论有:
离子实与价电子之间相互作用采用赝势来表示; 超晶胞的周期性边界条件; 平面波基组描述体系电子波函数;
广泛采用快速 fast Fourier tran sform （FFT） 对体系哈密顿量进行数值化计算;
体系电子自恰能量最小化采用迭带计算的方式;
采用最普遍使用的交换-相关泛函实现 DFT的计算，泛函含概了精确形式和屏蔽形式。
CASTEF中周期性结构计算优点
与MS中其他计算包不同，非周期性结构在
CASTEF中不能进行计算。将晶面或非周期性结构置于一个有限长度空间方盒中，按照周期性结构来处理，周期性空间方盒形状没有限制。之所以采用周期性结构原 因在于：依据 Bloch定理，周期性结构中每个电子波函数可以表示为一个波函数与晶体周期部分乘积的 形式。他们可以用以晶体倒易点阵矢量为波矢一系列分离平面波函数来展开。这样每个电子波函数就是 平面波和，但最主要的是可以极大简化 Kohn-Sham方程。这样动能是对角化的，与各种势函数可以表示
为相应Fourier形式。
.2
Den sity Orbital
（（
[
l GG Von（ G） *（ G） %（
kG
G
G
G G）]C
i,kG
ii,kG
CASTEF中引入外力
C
采用周期性结构的另一个优点是可以方便计算出原子位移引起的整体能量的变化，在
或压强进行计算是很方便的，可以有效实施几何结构优化和分子动力学的模拟。平面波基组可以直接达 到有效的收敛。
计算采用超晶胞结构的一个缺点是对于某些有单点限缺陷结构建立模型时，体系中的单个缺陷将以 无限缺陷阵列形式出现，因此在建立人为缺陷时，它们之间的相互距离应该足够的远，避免缺陷之间相 互作用影响计算结果。在计算表面结构时，切片模型应当足够的薄，减小切片间的人为相互作用。
CASTEP中采用的交换-相关泛函有局域密度近似（LDA（ LDA、广义梯度近似（GGA和非定域交换-相关
泛函。CASTE冲提供的唯一定域泛函是
CA-PZ Perdew and Zunger 将Ceperley and Alder 数值化结果
进行了参数拟和。交换-相关泛函的定域表示形式是目前较为准确的一种描述。
Name Descri pti on Refere nee Perdew and Wang P erdew et al. Hammer et al.
P W91 P erdew-Wa ng gen eralized-gradie nt app roximati on, P W91 PBE Perdew-Burke-Ernzerhof functional, PBE
RPBE Revised Perdew-Burke-Ernzerhof functional, RPBE
d
采用梯度校正的非定域或广义梯度近似泛函与电子密度梯度
/dr和电子密度
都有关，这样可以
同时提高能量和结构预测的准确性，但计算耗时。
泛函计算在本质上实际是相同的，但在电子密度变化迅速体系中
CASTEF中提供的非定域泛函有三种：
PBE泛函实用性更好;
PBE泛函与PW91 RPBE是特别用来
提高DFT描述金属表面吸附分子能量的泛函， White and Bird 描述了各种梯度校正泛函计算方法，利用 广义梯度近似计算总能量使用平面波基组与定域泛函相比并不直接。包含梯度近似的交换 时对电子密度数据的精度要求较高，对计算机内存占用会增大。通过采用与平面波基组总能量计算中分
-相关泛函计算 裂交换-相关能量采用一系列空间网格相一致的方法来定义交换
-相关势。
平面波基组(Plane wave basis set )
CASTEP计算理论总结+实例分析
