MS计算介电常数：Materials Studio CASTEP介电函数计算全流程

介电常数不是一个数字，它是一组矩阵——至少对于各向异性晶体是这样。本项目组在计算BaTiO₃的介电响应时，CASTEP输出的介电矩阵对角元分别是ε₁₁=5.2、ε₂₂=5.2、ε₃₃=120，非对角元为零。这个结果在组内引发了争论：ε₃₃=120是静态介电常数还是高频介电常数？为什么对角元之间相差20倍？要回答这些问题，必须从CASTEP的介电计算原理出发，理解每一个输出项的物理含义。

静态与高频：两个介电常数的故事

CASTEP计算介电常数有两条路径：密度泛函微扰理论（DFPT）和有限场法（Finite Electric Field）。两条路径的核心区别在于是否包含离子对极化的贡献。

DFPT路径直接计算Born有效电荷和介电矩阵的电子贡献部分，输出ε∞（高频/光学介电常数，仅含电子极化）。要得到ε₀（静态介电常数，含电子+离子极化），需要额外做声子计算获取离子贡献：

ε0=ε∞+4πV∑νZν∗⋅Zν∗ων2ε0​=ε∞​+V4π​∑ν​ων2​Zν∗​⋅Zν∗​​

其中Z*ν是Born有效电荷矩阵，ων是声子频率，V是晶胞体积。

有限场法则通过在体系上施加外部电场E，计算极化响应P=ε₀E，一步直接得到静态介电常数ε₀。但有限场法对晶胞形状敏感——施加电场后晶胞可能发生非物理畸变，需要固定晶胞形状（仅在内部做原子弛豫）。

本项目组的选择逻辑：对于需要分解电子和离子贡献的体系（如铁电体），DFPT路径更优，因为它天然分离了两种机制；对于只需要静态介电常数总值的体系（如绝缘涂层材料的介电性能评估），有限场法更省计算资源。

模型构建与对称性预处理

BaTiO₃在室温下是四方相（P4mm空间群），c/a≈1.01的微小四方畸变是铁电性的来源。CASTEP建模时最关键的步骤是确认对称性不被误判——如果CIF文件的对称性标注为立方相（Pm3m），CASTEP会按立方对称性计算介电矩阵，输出ε₁₁=ε₂₂=ε₃₃，直接抹掉了铁电畸变带来的巨量极化差异。

项目组的做法：从实验CIF导入结构后，先在CASTEP中做结构优化（允许原子位置和晶胞形状弛豫），确认优化后的对称群从Pm3m降为P4mm。只有确认了铁电畸变被正确保留，后续的介电计算才有意义。

参数设置的三个关键决策

ENCUT：介电计算比结构优化更敏感。

介电矩阵的计算需要精确的波函数展开，因为极化响应P∝⟨ψ|x|ψ⟩对基组完备性的敏感度高于总能量E∝⟨ψ|H|ψ⟩。项目组在BaTiO₃上做了ENCUT收敛测试：结构优化ENCUT=500eV时已收敛（ΔE<0.5meV/atom），但介电矩阵在ENCUT=500eV时ε₃₃=95，在ENCUT=600eV时ε₃₃=118——差异高达24%。介电计算的ENCUT至少要比结构优化高100eV，项目组默认ENCUT≥600eV。

K点密度：光学介电函数尤其依赖K点。

光学介电函数ε(ω)的计算需要在足够密的K点网格上做频率依赖的矩阵元积分。项目组对比了4×4×4和8×8×8的K点网格：在4×4×4网格下，ε∞的虚部在可见光区出现了明显的spurious peak（数值积分噪声导致的虚假吸收峰），8×8×8网格下这些峰消失。对于光学介电函数计算，K点密度至少要达到0.02Å⁻¹的分辨率，介电常数总值计算则6×6×6通常够用。

声子计算路径：DFPT vs 有限位移。

DFPT介电计算需要声子数据配合，CASTEP支持DFPT声子计算（线性响应法）和有限位移法（Phonopy风格的超胞法）。DFPT声子计算速度快（不需要超胞），但对低频声学支的精度取决于K点积分；有限位移法精度更可控（直接在Γ点附近做位移），但超胞尺寸需要足够大（至少2×2×2）。项目组在BaTiO₃上用两种方法对比：DFPT声子给出的软模频率ω=47cm⁻¹，有限位移法给出ω=42cm⁻¹——差异约12%，影响ε₀约5%。对于铁电体这类软模体系，推荐有限位移法配合2×2×2超胞。

结果分析：介电矩阵的物理解读

CASTEP输出介电矩阵后，解读工作才真正开始。BaTiO₃的结果：

ε₀与ε∞之间近100倍的差距，正是铁电体中离子极化（Ti离子沿c轴位移约0.12Å）贡献的巨量极化响应。离子贡献集中在软模频率ω≈42cm⁻¹附近——声子越软，离子对极化的贡献越大。

但ε₃₃=1200这个数值与实验值约1500-2000之间的偏差需要正视。CASTEP的PBE泛函对铁电体的极化响应系统性低估，原因在于PBE对Ti-O键的共价性描述不足——交换势的屏蔽效应让电荷转移偏小，Born有效电荷Z*偏低约10-15%。修正路径：用HSE06泛函重算，ε₃₃提升到约1400，与实验值更接近。

适用边界与注意事项

CASTEP介电常数计算有三个明确的适用边界：

第一，金属体系不能直接算介电常数。 CASTEP的介电计算假设体系是绝缘体（存在带隙），金属体系的Drude响应（自由电子贡献）需要额外的频率依赖处理，CASTEP标准流程不包含这一部分。如果误对金属体系运行介电计算，输出的ε∞会是一个无意义的极小负值。

第二，各向异性晶体的介电矩阵需要完整输出。 CASTEP默认输出对角化后的介电常数，对于非对角元不为零的体系（如单斜晶系），必须要求输出完整的3×3介电矩阵。项目组在计算LiNbO₃（三角晶系）时，默认对角化输出把ε₁₂≈0.3的非对角贡献抹掉了——虽然数值不大，但在相匹配器件设计中0.3的差异足以让折射率偏离预期。

第三，高频介电常数ε∞不能代替静态介电常数ε₀使用。 工程应用中需要的介电常数几乎都是静态值（ε₀），因为实际电场频率远低于晶格振动频率。ε∞仅对光学折射率有意义（n²≈ε∞），用ε∞代替ε₀会导致数值偏低数倍甚至数十倍。

CASTEP的介电常数计算模块封装了复杂的物理过程，但封装并不意味着可以跳过理解。从对称性保留到泛函精度，从声子配合到矩阵解读——每一步决策都决定了最终结果是否可以被信任。关于Materials Studio更多模块的计算指南，可在MS计算方法栏目中找到系统性整理。回到科研学术网首页探索更多计算材料学的实战路径。