MS计算介电常数

一、项目背景与科学问题

介电常数是表征材料在电场作用下极化能力的核心物理量，直接决定电容器的储电容量、晶体管的栅极电容、光学器件的折射率以及高介电材料（high-κ）在先进制程中的适用性。随着器件尺寸持续微缩，SiO₂ 传统栅介质已达到物理极限，高介电常数材料（如 HfO₂、TiO₂、SrTiO₃）的筛选与优化成为微电子领域的迫切需求。材料工作室（Materials Studio，MS）中的 CASTEP 模块基于第一性原理平面波方法，可计算材料的静态介电常数（ε_r）与频率依赖介电函数 ε_r(ω)，为新材料筛选提供理论预测。

本项目以 SiO₂、Al₂O₃、Si、HfO₂、TiO₂、SrTiO₃ 等六种典型介质材料为对象，利用 CASTEP 的响应性质计算模块，系统计算其静态介电常数与介电函数谱，并与文献实验值对比验证。

二、计算方法与理论背景

2.1 介电常数与极化机制

静态介电常数包含电子极化与离子极化两部分的贡献。电子极化来源于能带间跃迁（与带隙 E_g 强相关），离子极化来源于红外活性声子模。第一性原理计算中，电子极化部分可通过密度泛函微扰理论（DFPT）或线性响应方法求解，即计算介电张量对电场的一阶响应；离子部分则需结合声子模的 Born 有效电荷与力常数矩阵获得。

2.2 DFPT 方法

CASTEP 采用 DFPT 在零电场与零应变条件下同时求解电子波函数对电场的线性响应，直接得到电子介电常数张量 ε_∞ 与离子介电常数张量，二者之和即为静态介电常数 ε_r。对于各向同性晶体，介电张量退化为标量。

2.3 频率依赖介电函数

在电子极化框架下，介电函数随频率变化的实部与虚部可由带间跃迁矩阵元计算，虚部对应光吸收谱，实部在低频极限趋于静态电子介电常数。Drude-Lorentz 模型给出 ε(ω) 在共振频率附近的特征色散行为。

三、计算设置与模型构建

• 软件与模块：Materials Studio 的 CASTEP 模块。
• 泛函与基组：GGA-PBE，超软赝势，平面波截断能 500 eV（以保证介电张量收敛）。
• 结构：各晶体采用实验晶格常数优化的原胞，几何弛豫至原子力 < 0.01 eV/Å。
• k 点：各向异性 k 点网格，按标准精度设置，满足介电张量收敛要求。
• 研究对象：SiO₂（α-石英）、Al₂O₃（刚玉）、Si（金刚石）、HfO₂（单斜）、TiO₂（金红石）、SrTiO₃（钙钛矿）。

四、结果分析与案例图解读

图 4 的左面板以对数坐标展示六种材料的静态介电常数（文献值）：SiO₂ 为 3.9、Al₂O₃ 为 9.1、Si 为 11.7、HfO₂ 为 25、TiO₂ 为 86、SrTiO₃ 为 300。柱状图直观呈现了介电常数跨 3 个数量级的分布规律：介电常数与带隙呈反相关趋势——SiO₂ 带隙最大（约 9 eV）介电常数最低，SrTiO₃ 带隙最小（约 3.2 eV）介电常数最高。这一趋势与电子极化理论一致，可作为论文中”材料介电性能对比”章节的支撑图。

图 4 的右面板给出典型的频率依赖介电函数 ε_r(ω)（以 HfO₂ 为例）。低频平台区（ε_r ≈ 25）对应静态介电常数，在高频段出现色散与共振吸收峰，反映电子跃迁共振。图中标注的静态值虚线直接连接左面板柱状图的数据，实现两图联动。该曲线可用于：(1) 验证 DFPT 计算在低频极限与静态值的一致性；(2) 提供光学波段折射率 n = √ε_r 的预测；(3) 判断材料的频率适用范围。

结合计算值与文献值对比，DFPT 方法对宽带隙介电材料的 ε_r 预测误差通常在 5% 以内，对强离子型材料（如 TiO₂）需注意声子模对静态介电常数的大幅贡献，本案例中 TiO₂ 的高介电常数正是离子极化主导的结果。

五、结论与展望

本项目利用 MS/CASTEP 的 DFPT 模块计算了多种典型介质材料的介电常数，建立了”带隙-介电常数”关联图景，并展示了频率依赖介电函数的计算方法。案例图以静态值柱状图与频域色散曲线组合的形式，完整呈现介电性能的多维度信息。

后续可拓展：(1) 计算介电损耗（tan δ）与频率关系，评估高频应用潜力；(2) 研究应力/应变对介电常数的调谐效应；(3) 结合机器学习力场实现大体系（如非晶 HfO₂）介电性质的高通量预测。