DFT计算材料电荷密度

结合案例图分析，我把这个项目的完整过程、参数设置、结果解读和踩坑经验整理如下。

项目过程：DFT计算材料电荷密度

客户给的是一个类石墨烯的二维蜂窝状材料，想搞清楚相邻原子之间是不是共价键、电荷在成键区域如何分布。这个需求来源于审稿人质疑”你们说的层内强相互作用有没有电荷密度证据”。

我用的软件是 VASP 6.3.2，泛函 PBE，赝势 PAW。体系建的是 4×4 超胞，真空层厚度 18 Å，避免周期性镜像相互作用。为了得到平滑的电荷密度图，我把 ENCUT 提到 520 eV，K 点取 6×6×1，并且增大了 FFT 网格（NGXF=NGYF=168）。计算分两步：先做标准自洽得到 CHGCAR，然后用 VESTA 做二维切片和 1D 线积分。

具体参数：PBE，ENCUT=520 eV，ISMEAR=0，SIGMA=0.05 eV，EDIFF=1E-6。后处理时沿最近邻原子连线方向提取电荷差分 Δρ，同时做了平面切片展示二维分布。

案例分析：从数据到结论

第一张案例图是二维总电荷密度切片。可以看到每个原子中心是红色高值区（核电荷+内层电子密集），原子之间由黄色/绿色桥连接，说明存在明显的电荷共享。特别值得注意的是六边形蜂窝格子的每条边上都有连续的电荷云，这是共价键存在的直接可视化证据。

第二张案例图是沿 Si–Si 键连线的 1D 差分电荷密度分布 Δρ。横轴用 r/d_Si–Si 归一化，纵轴是 Δρ。曲线在键中心 r/d≈0 处出现明显的正峰，标注为”covalent bond charge build-up”；在原子核位置 r/d≈±1 附近出现负峰，标注为”charge depletion”（蓝色负值区）。这明确说明：成键区域有电荷积聚，而原子外侧有电荷耗尽。这种”中间积聚、两侧耗尽”的图案是典型共价键特征，和离子键的单一转移图案完全不同。

总结教训：复盘必须警惕的三件事

第一，差分电荷密度一定要定义清楚参考态。我用的是”成键后总密度减去原子密度叠加”，也就是 Δρ = ρ(material) – Σρ(atom)。如果参考态搞错，图的物理解释会完全相反。

第二，FFT 网格不能太粗。第一次用默认 NGXF，切片出现锯齿状伪影；把 NGXF 从 112 提到 168 后，成键桥才变得平滑连续。

第三，真空层厚度要足够。我试过 12 Å，能看到上下层镜像电荷的微弱干扰；提到 18 Å 后镜像效应消失。

项目经验：下次再做我会怎么做

以后再出电荷密度图，我会固定这个后处理流程：CHGCAR → VESTA 切片 → 1D 线积分 → 和原子参考态对比。交付时同时给两张图：一张二维切片让读者”看见”成键，一张 1D 曲线让读者”量化”电荷转移量。另外，对于审稿人，我会在正文里明确写出 Δρ 的参考态定义，避免被质疑图是怎么做出来的。

DFT计算材料电荷密度

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验 需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。