MS计算介电函数：从光学性质模块到频率依赖介电常数的完整路径

介电函数是材料光学性质的门户参数，MS计算介电函数的流程看似简单——几何优化后勾选光学性质计算即可——但真正拿到一组与实验吻合的频率依赖介电常数，远比点击几个按钮复杂。团队在钙钛矿体系上曾因k点网格选取不当，介电函数峰值位置偏移了0.8 eV，两周迭代后才意识到问题根源不在后处理，而在前置设置的物理合理性。

几何优化：介电函数计算的地基

介电函数对结构高度敏感。以BaTiO₃为例，铁电相与顺电相的静态介电常数差异可达两个数量级。团队的做法是先用PBE泛函完成粗优化，再精细调整晶格参数收敛标准到0.01 Å——这个数值基于对结构弛豫对介电响应影响的评估：当晶格常数变化超过0.02 Å时，静态介电常数的变化已超过5%。

对于带隙偏小的体系（PBE计算带隙常偏低30%-50%），几何优化阶段不必追求带隙精确，因为后续光学性质计算中有专门的修正机制。

光学性质模块设置：参数的物理逻辑

进入CASTEP的光学性质计算模块，三个核心参数需要逐一审视：

k点网格密度。介电函数的计算需要足够密的k点来描述布里渊区中光学跃迁的细节。团队在BaTiO₃体系中采用8×8×8的网格，对比6×6×6的结果发现介电函数虚部ε₂(ω)在4.2 eV处的峰值强度差异达15%——差距不会说谎，k点不足直接导致跃迁矩阵元采样不完整。

scissors operator。这是CASTEP处理带隙偏低问题的内置方案——将计算带隙人为”剪开”到实验值，所有光学跃迁能量同步上移。团队选取scissors值为1.12 eV（实验带隙3.2 eV减去计算带隙2.08 eV），而非随意填一个近似值。选取依据是实验光学吸收边的精确位置。

频率范围与步长。介电函数的频率范围应覆盖材料所有重要光学跃迁。对氧化物体系，设定0-30 eV范围、0.1 eV步长，确保低频红外区域和高频电子跃迁都不遗漏。

频率依赖介电函数的物理含义

拿到结果后，介电函数实部ε₁(ω)和虚部ε₂(ω)各有明确的物理指向。ε₂(ω)直接反映光学吸收谱的轮廓，每一个峰值对应一次带间跃迁；ε₁(ω)通过Kramers-Kronig关系与ε₂(ω)耦合，静态介电常数ε₁(0)是评估材料极化能力的关键指标。

团队计算得到的BaTiO₃铁电相ε₁(0)≈340，与实验值约400的差异在15%以内。偏差的来源值得剖析：PBE泛函对电子局域化的描述不足是主因，scissors operator只修正跃迁能量不修正跃迁矩阵元——两者叠加造就了这个系统性偏低。认定当前方案适合定性趋势分析和相对比较，精确数值仍需更高精度泛函校准。

关于CASTEP光学性质计算方法的理论基础，可参考Segall等人对CASTEP方法的系统阐述以及Levine对介电函数Kramers-Kronig变换的推导。

方法局限与适用边界

MS计算介电函数的适用边界需要坦诚面对：金属体系由于费米面处的非零占据，CASTEP的光学性质模块计算结果可靠性显著下降；强关联体系（含Cu²⁺的过渡金属氧化物）中DFT+U的参数选取本身就是一个未完全解决的问题，介电函数的计算更增添了不确定性。

回过头看这两周的反复迭代，核心教训是：参数设置不应是”试错-看结果-再调”的循环，而应从物理含义出发做有依据的选择。实践讨论，也可以访问https://www.keyanxueshu.com/浏览完整方法专题体系。