静态自洽电子结构计算

一、项目背景与科学问题

静态自洽（Static SCF）电子结构计算是几乎所有第一性原理模拟的基础步骤。在给定原子构型的条件下，通过迭代求解 Kohn-Sham 方程获得自洽的电子密度与总能量，得到体系的基态性质。这一”单点”计算虽不涉及结构弛豫，却是能带结构、态密度、电荷密度分析以及后续分子动力学、过渡态搜索等一切性质计算的出发点。SCF 迭代的收敛质量直接决定所有下游结果的可靠性。

本项目以体相硅（Si，金刚石结构，8 原子晶胞）为模型体系，开展静态自洽电子结构计算，系统分析 SCF 迭代过程中的能量收敛行为、收敛判据的设置以及收敛后的电子结构特征。

二、计算方法与理论背景

2.1 Kohn-Sham 方程与自洽迭代

Kohn-Sham 方程为 [−½∇² + V_ext(r) + V_H(r) + V_xc(r)]ψ_i(r) = ε_iψ_i(r)，其中有效势依赖电子密度 n(r)，而 n(r) 又由波函数求和得到，因此必须通过自洽迭代求解：给定初始密度 → 构造有效势 → 求解本征方程 → 由新波函数构造新密度 → 与旧密度混合 → 重复直至收敛。每次迭代间的能量差 |ΔE| 是衡量收敛的标准量。

2.2 收敛判据与混合方案

能量收敛判据通常设为 1×10⁻⁵ 至 1×10⁻⁶ eV，配合电荷密度混合方案（如 Pulay/DIIS 混合）加速收敛。对半导体（如 Si），收敛通常较快；对金属，需要高斯展宽或费米能级占据处理以保证迭代稳定。SCF 收敛曲线在半对数坐标上呈近似线性下降，反映误差的指数衰减。

2.3 总能与电子结构

自洽收敛后的总能量为基态能量（本案例 Si 为 −21.34 eV/atom），可作为结构优化的目标函数。同时可获得能带、态密度、电荷密度等电子结构信息，其中价带顶与导带底的相对位置给出带隙。

三、计算设置与模型构建

• 软件：VASP；泛函：GGA-PBE，PAW 赝势，截断能 400 eV。
• 结构：Si 金刚石原胞（2 原子）或 8 原子晶胞，实验晶格常数。
• k 点：Monkhorst–Pack 网格，标准精度。
• 收敛：能量 < 1×10⁻⁵ eV（输出曲线显示至 1×10⁻⁶）。
• 迭代：最大 60 步，记录每步能量。

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

图 6 的左面板给出 SCF 迭代过程中总能随迭代步数的演化（蓝色曲线）。初始几步能量从初猜快速下降，经过约 30 步后趋于平台，最终收敛至 −21.34 eV/atom（红色虚线标注）。曲线形态呈”快-慢”两段式：前 10 步能量下降占主导，后 30 步为精细收敛。该曲线是论文中”计算收敛性验证”的标准图，用于证明总能量已达收敛，后续性质计算可靠。

图 6 的右面板以半对数坐标给出相邻迭代能量差 |ΔE| 的下降（红色曲线）。曲线从约 10⁻¹ eV 量级指数衰减至 10⁻⁶ eV 以下（约 40 步后达到收敛线 1×10⁻⁶ eV，灰色虚线标注），线性下降说明迭代误差呈指数收敛，混合方案有效。该收敛曲线直接支撑”计算已收敛至指定精度”的结论，是审稿人重点关注的计算质量证据。

两面板结合，左面板给出总能数值（用于性质计算），右面板给出收敛质量的定量证明（用于方法学验证），共同构成静态自洽计算的完整记录：不仅报告”算出了什么”（总能量），还证明”算得足够准”（收敛判据满足）。需要指出的是，SCF 收敛判定不能只依赖总能量差，还应配合电荷密度差或力的收敛检查，避免出现”能量不动但电荷未收敛”的假收敛情形；本案例中能量差与电荷密度均达到设定判据，可以判定为可靠收敛。

五、结论与展望

本项目完成了 Si 的静态自洽电子结构计算，总能量收敛至 −21.34 eV/atom，|ΔE| 达 1×10⁻⁶ eV。案例图以能量收敛曲线与 |ΔE| 对数曲线双面板呈现。

后续可拓展：(1) 由收敛后的电荷密度计算能带、态密度与电荷分布；(2) 开展体积扫描获得状态方程与平衡晶格常数；(3) 将该静态计算结果作为弛豫、声子或 MD 模拟的起点。