DFT 计算材料电荷密度

这张图是本项目用 DFT 算某种材料 (001) 晶面电荷密度的切片结果。颜色越亮表示电子密度越高，暗区则接近 0。图里最显眼的规律是：每个原子核周围都有一个亮黄的”芯”，往外迅速衰减到红紫相间的电子云，再过渡到晶面间隙的暗区。这种图不是拿来”好看”的，而是用来判断成键类型、电荷转移方向、表面态位置的。

这张图在说什么

亮黄色集中在原子核附近，说明芯电子被赝势或全电子波函数紧紧束缚；红橙过渡区是价电子云，相邻原子之间的红橙桥接意味着成键区有电子共享。如果两个原子之间是暗带，说明键合弱或离子性强、电子被拉向一侧。本案例中相邻原子间有连续的电子云桥，价电子分布比较均匀，说明体系偏向共价/金属键合，而不是离子键。

右上角颜色条把电子密度归一化到 0–1.0，这种归一化便于直观比较不同体系，但做定量分析时必须回到原始 CHGCAR 数值。

电荷密度分析有哪几种用法

本项目最常从 CHGCAR 里提取三类信息：

| 分析类型 | 做法 | 能回答什么 |

|———-|——|——|

| 差分电荷密度 Δρ | ρ(AB) − ρ(A) − ρ(B) | 电荷转移方向与成键极性 |

| 平面平均电荷密度 ρ̄(z) | 沿 z 方向平均 | 表面/界面电荷重新分布、偶极矩 |

| Bader / Mulliken 电荷 | 划分原子电荷 | 定量氧化态、电荷转移量 |

其中差分电荷密度最能讲故事：电子从哪边流向哪边，一目了然。

VASP 关键参数

| 参数 | 取值 | 说明 |

|——|——|——|

| LCHARG | .TRUE. | 输出 CHGCAR |

| NGXF / NGYF / NGZF | 加密 FFT 网格 | 电荷密度图光滑、积分误差小 |

| PREC | Accurate | 防止网格截断失真 |

| ISPIN | 2 | 磁性体系必须开，否则自旋密度丢失 |

NGXF 系列参数如果设太粗，画出来的等值面会有锯齿，Bader 电荷划分误差也会增大。本项目习惯在标准 PREC=Accurate 基础上再把 NGXF 提高 20–30%。

怎么做二维切片

本项目用 VESTA 或 Python（`pymatgen` / `ASE`）读取 CHGCAR，选定一个包含关键原子的晶面做 slice。要点：

– 切片厚度不要太薄：至少 1–2 个 FFT 网格点厚度，避免单点噪声。

– 等值面/色标要固定：不同体系用同一色标才能比较。

– 差分电荷密度做 slice 时，确保三个子体系的网格完全一致（同一 NGXF）。

从图里读什么

本图能看到：

1. 电子局域在原子核附近，但价区有一定延展；
2. 最近邻原子间有连续电子云，说明键合较强；
3. 晶面间隙区域暗，说明层间/间隙电荷少；

如果这是二维材料，间隙暗说明层间耦合弱；如果这是表面吸附体系，还要看吸附原子下方是否出现新的电荷堆积峰。

常见坑

– **归一化陷阱**：彩色图归一化到 1.0 会抹掉绝对密度差异，定量必须用 CHGCAR 原始值。

– **赝势影响**：PAW 赝势把芯电子冻结，芯区电荷不是真实全电子密度；差分密度分析应避开芯区。

– **自旋密度**：磁性体系只做电荷密度图不够，必须补自旋密度 ρ↑−ρ↓，否则磁性起源完全丢失。

DFT 计算材料电荷密度是判断”电子去哪了”最直接的证据。更多电荷密度分析与可视化脚本可参考科研学术网 [DFT计算材料电荷密度](https://www.keyanxueshu.com/category/dft/) 栏目，完整流程见首页 https://www.keyanxueshu.com/。

局限

电荷密度只是基态单粒子密度；激发态、动态电荷涨落、强关联体系里的 d/f 电子局域性需要更高级方法（GW、DMFT、BSE）。