真空层厚度如何影响 DFT 计算结果？

说明：本文华算科技主要介绍真空层厚度对中性薄层能量、带电缺陷、二维介电响应及 GW/BSE 收敛行为的影响。

真空层厚度改变的是哪个计算盒子？

薄膜与二维材料仍在三维周期条件下求解。一个含单层材料的超胞会沿 z 方向无限复制，真空层把相邻周期映像的原子与电子密度尾部隔开。若薄层几何厚度为 d_2D，相邻层之间的空隙为 d_vac，超胞高度可写为 L_z=d_2D+d_vac。这里的 d_vac应按最上方与最下方原子之间的实际间距定义，不能把整个晶格常数 c 直接称为真空厚度。

真空区的电荷密度很低，却仍属于数值计算的盒子。增大 L_z会扩大体积 Ω，使倒格矢 G_z 的间隔变为更密的 2π/L_z，也会增加固定平面波截断能下的基函数与 FFT 网格点。长程库仑核、宏观介电量的体积归一化和带电超胞的补偿背景都会随之改变。只把 k_z 取成 1，并不会取消 z 方向的周期复制。

图 1 把这三个长度分开后，真空层收敛的对象也更明确：它检验的是周期映像相互作用、长程核和体积归一化误差，而薄层自身厚度检验的是内部原子层数是否足以代表目标表面或体相。两种收敛测试改变的物理模型不同，一条总能曲线无法兼顾。

中性体系、准粒子能量和带电缺陷为何表现不同？

对无净电荷、无强垂直偶极的薄层，半局域 DFT 总能、原子力和 Kohn–Sham 带边常在相邻密度尾部不再明显重叠后趋于稳定。吸附能或表面能属于能量差，参与相减的各个模型还要保持相同的 L_z、面内晶胞和数值精度，否则基函数数量与归一化方式的变化也会进入差值。范德华作用、精确交换、自能和屏蔽库仑作用具有更长的空间范围，不能沿用半局域总能已经收敛的真空间距。

图 2 中，未截断库仑作用的 GW 带隙随着周期层间距增大仍持续移动；同一层间距下，9×9、15×15、21×21 与 33×33 网格也没有汇到同一点。二维屏蔽在 q→0 附近变化很陡，增大真空会让这一窄区更难由粗网格表示。真空厚度与面内 k/q 点密度在 GW 中属于耦合收敛量，单独加宽真空可能把采样误差放大。

图 3 的等值线没有平行于任一坐标轴，说明只沿一个方向测试容易遇到偶然平坦区。某个粗 k 网格下，两组真空间距得到相近带隙，可能来自层间屏蔽误差与布里渊区积分误差相互抵消；换用更密网格后，所谓“收敛值”仍会移动。检查二维材料的 GW 结果时，至少要在两个真空间距上比较多组面内采样，或采用与超胞周期长度相匹配的二维库仑截断。

带净电荷的周期超胞通常引入均匀补偿背景，以维持泊松方程可解。局域缺陷电荷与分布在整个超胞中的背景形成虚假的静电电容；背景体积随 L_z 增大，未校正形成能便可能近似线性漂移。此时更厚的真空并不自动接近孤立带电层极限。带电缺陷需要针对二维几何的有限尺寸校正或库仑截断，还要确认缺陷电荷确实局域，价带或导带中的离域载流子不能硬套局域电荷模型。

图 4 的中性硫空位对 10、15 和 20 Å 真空间距不敏感，q=+1 与 q=−1 的空心点却明显分开；加入二维静电校正后，三组实心点大幅靠拢。相同原子结构因为电荷态不同，真空依赖便从密度尾部问题转成泊松方程的周期约束与补偿背景问题，不能用中性缺陷的收敛曲线替代。

二维介电/光学强度为何随真空层“变小”

常规平面波程序输出的宏观介电函数按三维超胞体积 Ω 归一化。薄层本身的面内极化响应固定时，真空增大等于把同一份诱导偶极矩摊进更大的体积。静态面内分量可写成 ε_∥^SC−1∝α_2D/L_z，其中α_2D 是单位面积二维极化率。ε^SC 随 L_z 增大趋近 1，是超胞归一化的结果，不表示材料的电子极化能力被真空削弱。

图 5 中各条曲线都向 1 靠近，较大的面内波矢 q 具有较强屏蔽，却仍保留明显的 L 依赖。若两篇工作使用 15 Å 和 30 Å 真空，直接比较它们的 ε_∥^SC 会把超胞高度差误当成材料差。可比对象应换成 α_2D、片层电导或明确采用物理厚度换算的有效介电常数，并在文中给出厚度定义。

电子结构已对真空间距收敛后，Im[ε(ω)] 以及 Re[ε(ω)]−1 通常与 Ω⁻¹ 成正比。真空从 15 Å 增至 30 Å 后，吸收峰位置可以几乎不动，峰高却接近减半；把这种下降解释成跃迁概率变小，物理对象就错了。峰位收敛与谱强度归一化必须分开检查。

图 6 的虚线保留超胞体积，实线按二维响应恢复后再与实验比较。恢复操作只处理真空稀释因子，不会修正泛函、激子核、展宽参数或实验环境造成的峰位偏差。谱强度经 L_z 缩放后仍与实验不符时，应继续检查电子结构和多体效应，不能再归因于真空厚度。

怎样为目标结果确定真空层收敛点？

真空层没有适用于所有任务的固定厚度。几何优化关心最大残余力与层间伪相互作用，吸附和表面反应关心能量差，功函数关心真空静电势平台，带电缺陷关心校正后的形成能，介电和光学计算关心二维化后的响应，GW/BSE 还关心库仑截断与小 q 区域的积分。停止标准应写成目标量允许变化多少，不能只写“使用了 20 Å 真空”。

收敛序列要固定薄层结构、面内晶胞、赝势、截断能、占据与展宽方式，再逐级改变 d_vac。非对称薄膜应分别检查两侧真空势是否出现平坦区，并区分偶极校正前后的能量与势能参照。能量差计算的各个构型使用同一 L_z，可减少纯数值体积变化混入吸附能或反应能。

图 7 中，体相 cBN 的 ε₀₀⁻¹(q) 在 Γ 点附近平缓，单层 hBN 却出现急剧上翘。普通均匀网格可能用很少的点代表整块 Γ 点微区，真空增大或启用二维库仑截断后，小 q 极限还会改变。GW 的收敛记录需列出L_z、截断方式与面内网格；其中任一项变化，都要重新判断另外两项的误差。

图 8 上部的光学带隙随层间距变化较小，下部的激子结合能却仍大幅移动。光学带隙由准粒子带隙减去激子结合能得到，两项同向误差可能在相减时抵消。若研究目标是激子束缚，光学峰看似稳定不能作为收敛证据；若只讨论光学跃迁位置，也要报告组成它的两项是否各自稳定，防止结论依赖误差抵消。

最终记录至少包含实际真空间距的定义、测试过的 L_z 序列、被监测的物理量及容差。带电模型补充电荷校正与定位检验，二维响应补充由 ε^SC 转成 α_2D 或片层量的公式，GW/BSE 补充库仑截断和 k/q 网格。读者据此才能分辨结果已接近孤立薄层极限，还是恰好停在一处数值平台。