实例分享：MS castep的参数设置与计算电子性质的操作步骤（免费、新手向）

目录

 理论基础：第一性原理与DFT概述

参数设置

建模优化、计算性质的操作步骤

---

基于密度泛函理论和赝势理论的第一性原理的计算方法，采用MS castep模块，构建模型及结构优化，计算并分析了能带结构（Band Structure）、态密度（Density of State，DOS）、电荷密度（Electron density）、布居分布（Population Analysis）等电子性质。研究技术路线为：文献调研→建模优化→性质计算→数据分析。基于密度泛函理论和赝势理论的第一性原理的计算方法：

参数设置

软件中关键参数需手动设定，如泛函设定、赝势选择、截断能（Ecut）、布里渊区k点取样、是否包含自旋极化、以及计算精度等级等，其它参数直接使用系统值。

常见的交换关联势：局域密度近似(LDA)[34]和广义梯度近似(GGA)，LDA：均匀电子气密度n(r)；GGA：还兼顾了周围区域的电子密度，从而提高了准确性。GGA、LDA计算禁带宽度偏小，HSE弥补了缺陷，但计算成本增加。因为只研究禁带宽度的变化趋势，使用GGA PBE。

赝势：模守恒假势、超软赝势(ultrasoft pseudopotentials)。模守恒赝势要求在一定区域内，电子对于赝势和全电子函数保持一致，以增加赝势的精确度，但需要较高的Ecut；超软赝势可显著降低运算负担，却不必牺牲计算的精准程度，但是需依依靠一系列经验参数。

计算精度等级：根据计算能力、计算耗时、体系的大小复杂程度，因为本体系原子数极少，故采用fine。结构优化时，建议阶梯推进计算精度：先使用coarse，在其结果文件夹重新使用medium迭代计算，在medium结果文件夹使用fine作为最终结构优化的结果。

收敛性测试/查阅前人文献：得到截断能(Ecut)和k点的参数设置。在计算时，Ecut设置越大，计算结果也更精确，可计算时耗更严重；k点取样越密，虽然计算越准确，可计算时间剧烈增加。故进行收敛性测试，来寻找精度良好࿰