Castep的基本原理和参数设置.ppt

单击此处编辑母版标题样式 单击此处编辑母版副标题样式 单击此处编辑母版标题样式 单击此处编辑母版文本样式 第二级 第三级 第四级 第五级 CASTEP的基本原理和参数设置 CASTEP将电子波函数通过平面波函数展开 赝势 pseudo wave function pseudopotential 有效的减少平面波数目 只适用于周期性体系 计算任务(Task) Energy, Geometry optimization, Elastic constant, TS search, Properties 精度控制(Quality) 交换相关函数(Functional) LDA: CA-PZ GGA: PBE, RPBE, PW91, WC, PBEsol B3LYP 0PBE HF, HF-LDA sX, sX-LDA 自旋极化(Spin polarized) 初始自旋值(Initial spin) 体系电荷 动能截断值(Energy cutoff) 自洽计算收敛精度(SCF tolerance) K点设置(k- point set) 赝势(Pseudopotential) Ultrasoft, Normconserving, On the fly 赝势的表达，实空间 or 倒空间 Core hole的选择 More…… Basis ——基组的相关参数设置 动能截断值的设定 FFT格子设定 基组的测试(Finite basis correction ) SCF ——自洽计算的相关参数设置 k-points Potential ——选择赝势 Basis ——基组的相关参数设置 SCF ——自洽计算的相关参数设置 收敛精度的具体设定 自洽计算允许的步数 Density mixing 轨道占据数 初始自旋值的优化 k-points Potential ——选择赝势 Density Mixing参数控制体系中，如何根据特征方程来如何构造新的电子密度。在整个体系中，通过加入阻尼振荡来确保整个体系的平滑收敛。对于一个最简单的阻尼方法，遵循如下方程： A就是阻尼系数。ρin是当前前一步SCF迭代计算得到的电子密度，根据这一项来构建特征矩阵。ρout是没有考虑阻尼系数，根据新的分子轨道系数构建得到的电子密度。而ρnew则是真正用于下一步迭代的电子密度。 DIIS(direct inversion in an iterative subspace) ： 也称为Pulay mixing方法，由Peter Pulay发展得到，是一种外推方法。用于Harter-Fock SCF迭代方法的加速并稳定体系的收敛过程中。 当选中时，能够加速整个体系的SCF迭代速度。 The DIIS method assumes that a good approximation to the final solution pf can be obtained as a linear combination of the previous guess vectors where m is the number of previous vectors (in practice, only the most recent few vectors are used). The coefficients ci are obtained by requiring that the associated residual vector  approximates the zero vector in a least-squares sense. Furthermore, the coefficients are required to add to one,  初始的情况 采用Smearing后 初始的情况 采用Smearing后 Final energy, E =