MS计算介电常数的CASTEP模块操作详解

为什么介电常数计算容易出错

介电常数是材料光学和电子性质的重要参数，但用Materials Studio（以下简称MS）的CASTEP模块计算时，结果对参数极为敏感。本项目近期完成了一组钙钛矿材料的介电常数计算任务，同样的BaTiO₃结构，不同参数设置下静态介电常数从15到120不等——差了8倍。问题出在K点密度、SCF收敛精度和声子贡献三个方面。

介电常数的完整计算包含两部分：

• 电子贡献（高频介电常数ε∞）：来自电子极化，CASTEP直接通过DFPT计算
• 离子贡献（晶格极化）：来自原子位移极化，需要声子频率和Born有效电荷

总静态介电常数 ε_static = ε∞ + ε_ionic

很多研究者只算了电子贡献就当最终结果报告，这在离子晶体中会严重低估。以NaCl为例，ε∞≈2.3，但ε_static≈5.9，离子贡献占了一大半。

CASTEP介电常数计算的前置条件

在MS中计算介电常数之前，必须完成以下准备：

1. 高精度结构优化

结构优化的精度直接决定介电常数计算的可靠性。本项目推荐的参数：

这套参数比普通结构优化严格得多，计算量也大3-5倍，但对于介电常数计算是必要的。

2. 能带间隙检查

介电常数与能带间隙存在反比关系（Penn模型：ε∞ ≈ 1 + (ℏω_p/E_g)²），因此先检查GGA-PBE计算的带隙是否合理。如果PBE严重低估带隙（常见于半导体），需要考虑用HSE06或GW修正。本项目在TiO₂体系中实测：

• PBE带隙：1.86 eV（实验值3.0 eV，严重低估）
• HSE06带隙：2.95 eV（接近实验值）
• PBE ε∞：8.7（实验值7.4，偏高）
• HSE06 ε∞：7.1（更接近实验）

所以对于带隙敏感的体系，用HSE06算介电常数更靠谱。关于DFT计算中泛函选择的影响，本项目在其他文章中有系统讨论。

CASTEP操作流程

步骤一：结构优化

在MS界面中：CASTEP → Calculation → Task: Geometry Optimization → 设置上述参数 → 运行。

关键勾选项：

• ✅ Optimize cell（晶胞参数优化）
• ✅ Quality: Fine 或 Ultra-fine
• ✅ More… → 电子最小化方法：All Bands/EDFT（对金属体系）

步骤二：介电常数计算

在优化后的结构基础上：

CASTEP → Calculation → Task: Properties → 勾选 Dielectric function

参数设置：

注意：对于非立方晶系，介电常数是张量，有多个独立分量。MS会输出完整的3×3介电张量，需要根据晶体对称性取平均值。

步骤三：声子计算（离子贡献）

静态介电常数的离子贡献需要声子频率和Born有效电荷：

ε_ionic = (4πe²/V) × Σ_m (Z*_m)² / (M_m × ω_m²)

其中Z*是Born有效电荷张量，ω是Γ点声子频率，V是晶胞体积。

CASTEP在介电计算中会自动计算这些量，但需要确保：

• 声子计算在Γ点完成
• 没有虚频（虚频意味着结构不稳定）
• Born有效电荷张量的对称性正确

TiO₂金红石实例

本项目以金红石TiO₂为例展示完整计算流程：

结构参数：a=4.594Å, c=2.959Å, 空间群P4₂/mnm

优化参数：PBE, 600 eV, 6×6×9 K-mesh, SCF=1E-8

结果：

金红石TiO₂的静态介电常数极大（~100），主要来自离子贡献——Ti⁴⁺和O²⁻的Born有效电荷远大于形式电荷（Z*_Ti ≈ 7.1 vs +4），表明Ti-O键有显著的共价性。这个结果在科研学术网的材料数据库中已有收录。

常见问题排查

问题1：介电常数远大于实验值

最常见原因是K点不够密或SCF没收敛。介电常数对波函数的导数敏感，任何数值噪声都会被放大。解决方案：提高K点密度到0.02/Å，SCF tolerance到1E-9。

问题2：介电张量不对称

如果晶体有对称性但计算出的介电张量不满足对称性约束，说明计算精度不够。检查：结构优化是否充分（应力<0.05GPa）、K点网格是否满足对称性。

问题3：声子有虚频

虚频意味着当前结构不是真正的能量极小值。可能原因：结构优化未收敛、泛函选择不当、声子计算的超胞太小。解决方案：重新优化（更严参数）或换泛函。

问题4：CASTEP报内存不足

介电计算内存消耗大，特别是高K点密度时。解决方案：减少并行进程数（增加每进程内存）、降低K点密度（先粗算验证流程）。

与实验对比的注意事项

1. 单晶 vs 多晶：CASTEP计算的是单晶介电张量，多晶样品的介电常数是各方向的平均值。对四方晶系：ε_avg = (2ε⊥ + ε∥)/3
2. 频率依赖：CASTEP给出的是特定频率下的介电常数（静态和光学），实验测量可能在不同频率下进行
3. 温度效应：DFT计算是0K结果，室温下的热膨胀和声子软化会改变介电常数
4. 缺陷效应：实际晶体中的缺陷和杂质会显著影响介电性能，DFT计算的完美晶体结果只能作为理想上限

MS计算介电常数是一项精度要求很高的工作，从结构优化到参数设置再到结果分析，每个环节都需要严格把关。上述参数和流程是本项目在多个材料体系中反复验证后的推荐配置，实际应用中可根据具体材料类型做适当调整。