手机版

castep计算静电势

发布时间：2026-08-27 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

静电势（Electrostatic Potential, ESP）是固体表界面电子结构分析的核心物理量。沿表面法线方向的静电势分布直接决定了功函数、电子亲和势、界面能带对齐以及吸附分子与表面的静电相互作用。在光催化（如 TiO₂ 光解水）、钙钛矿太阳能电池和半导体器件设计中，准确获得表面的静电势剖面与真空能级位置，是预测界面电荷转移与能带弯曲的关键前提。

本项目以锐钛矿 TiO₂(101) 表面为对象，利用 CASTEP 模块计算其平面平均（planar average）与宏观平均（macroscopic average）静电势，提取真空能级 E_vac、费米能级 E_F 与功函数 Φ，并分析表面态密度与能带结构。

二、计算方法与理论背景

2.1 静电势的第一性原理计算

DFT 自洽计算得到电子密度 n(r) 后，静电势由两部分组成：离子实产生的局域势与电子密度产生的 Hartree 势，二者通过求解泊松方程获得。周期性 slab 计算中，静电势在真空区域趋于平台值，该平台值即为真空能级 E_vac 的参考。

2.2 平面平均与宏观平均

平面平均势定义为 V(z) = (1/A)∫V(r)dxdy，即在平行表面方向取平均，得到沿法线 z 的一维分布。宏观平均势进一步对 V(z) 作周期平滑（如 Hanning 窗卷积），消除原子层内的振荡，凸显层间势垒与表面区域的平均势特征，是功函数与能带对齐分析的标准工具。

2.3 功函数提取

Φ = E_vac − E_F。其中 E_F 由自洽计算的费米能级给出，E_vac 取真空区宏观平均势的平台值。表面偶极层会使 E_vac 相对于体相参考发生偏移，因此功函数直接反映表面电荷重排效应。

三、计算设置与模型构建

软件：Materials Studio 的 CASTEP 模块。
表面：锐钛矿 TiO₂(101) 面，采用周期性 slab，厚度 6 个 TiO₂ 层（约 10 Å），真空层 15–20 Å。
泛函与基组：GGA-PBE，超软赝势，截断能 500 eV。
弛豫：顶层原子弛豫至原子力 < 0.02 eV/Å，底层固定模拟体相。
计算输出：LOCPOT（或 CASTEP 电位文件），DOS 与能带结构。

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

图 6 的左面板给出静电势沿表面法线（c 轴）的分布。蓝色细线为平面平均势 V(z)，在原子层区域呈现周期振荡，对应各原子层的离子势与电子势；红色粗线为宏观平均势，平滑后清晰显示三个区域：真空区平台（约 6.9 eV）、表面区与体相区。真空平台值对应 E_vac，红色虚线标注的费米能级约 2.3 eV。图中橙色双向箭头标出 Φ = E_vac − E_F ≈ 4.6 eV，与文献报道的锐钛矿 TiO₂(101) 功函数（4.3–4.9 eV）一致，验证了计算方案的可靠性。

图 6 的右面板给出表面附近的总态密度（DOS）示意：价带（VB）态在 −0.6 eV 以下，导带（CB）态在 2.2 eV 以上，费米能级位于禁带内。标注的带隙 E_g = 2.8 eV 为 PBE 计算值（文献 PBE 值约 2.2–2.9 eV，低于实验值 3.2 eV，属已知的低估效应）。该面板将静电势分析（左）与电子结构分析（右）关联，说明真空能级与带边的相对位置共同决定了光催化氧化还原能力。

左面板中，真空区与体相区之间的宏观平均势降幅即表面功函数；而表面偶极层（表面几个原子层内宏观平均势的倾斜段）正是 TiO₂ 表面发生电荷重排、形成表面态的证据。这一静电势剖面对论文中”表面电子性质”章节极具说服力。