什么是电荷密度？第一性原理计算定义、来源、可视化图谱及分析范围解析

说明：本文华算科技主要介绍电荷密度在第一性原理计算中的定义、计算来源、常见图谱和判断范围。

电荷密度的基础概念是什么？

在 DFT计算里，电荷密度 ρ(r) 描述电子在实空间位置 r 附近的分布。它描述的对象并非某个原子名下的固定电荷数；连续场铺在整个晶胞、分子或界面模型中。电荷密度的对象是空间分布，它把电子云、键区、空位邻域和界面区域放到同一个三维网格里。

Kohn-Sham 方程给出一组轨道 ψ_i(r) 和占据数 f_i。软件把占据轨道的模平方求和，得到 ρ(r) = Σ_i f_i|ψ_i(r)|²。轨道属于能量和量子态描述，电荷密度属于实空间描述，两者来自同一电子结构计算，但读出的物理对象不同。

平面波程序通常把密度写成网格文件，例如 VASP 的 CHGCAR 或 PARCHG。总电荷密度可包含核心电子和价电子贡献，价电子密度更贴近赝势计算中显式处理的电子。总密度、价电子密度、自旋密度和部分电荷密度在文件口径上不同，混用时会改变图谱含义。

bcc Fe-H 体系里的密度图把 Fe 原子骨架和 H 引入后的电子分布放在同一晶胞内。颜色或等值面并非“某原子得失电子”的直接标签；密度高的区域通常靠近原子核和成键区域，相邻模型的差异才支撑电子重新排布的判断，这类比较优先于单张图的颜色深浅。

若计算采用自旋极化，ρ_↑(r) 与 ρ_↓(r) 还可相减得到自旋密度。若只选取某条能带、某个能量窗口或某个缺陷态，得到的是部分电荷密度。同一个电子体系可以导出多种密度图，图名和生成口径必须一一对应。

电荷密度图由哪些计算结果构成？

电荷密度可直接来自收敛后的电子基态，也可通过数学组合得到新的图谱。总密度 ρ_system 描述完整体系；差分密度 Δρ 描述重排量；平面平均密度把三维场压缩成沿界面法向的曲线。图谱名称变了，参与相减或积分的密度文件也随之改变，这决定了它回答的是空间分布、重排方向还是层向变化。

界面或吸附体系常用 Δρ = ρ_interface – ρ_slab – ρ_adsorbate。片段密度通常要在吸附后几何、相同晶胞、相同网格和相同计算设置下生成。参考片段几何不同，Δρ 的正负区域会混入结构弛豫带来的差异，电子重排信号就会偏离原本的吸附或界面问题。

Fe(111)/Ti₃SiC₂(0001)界面图把电荷密度和差分电荷密度并列给出。前者描绘界面附近电子主要分布在哪里，后者突出形成界面后哪些区域出现积累或耗尽。总密度看电子云本身，差分密度看相对于参考片段的变化量，二者在文章中应分开使用。

把 Δρ(x,y,z) 沿界面平面积分，可得到 Δρ(z) 曲线。层状材料、异质结和金属-半导体接触中，这类曲线能定位电荷重排主要发生在界面层、吸附层还是远离界面的基底层。平面平均曲线牺牲横向细节，换来法向位置和符号的定量比较，适合处理周期性 slab 模型。

部分电荷密度来自选定能带或能量窗口。它只提取某个带边、缺陷态、界面态或反应相关轨道在实空间中的分布，已经脱离全体占据电子的总和。能量窗口的上下限会直接决定图中出现哪些电子态，PARCHG、LDOS 和带分解电荷密度应和能带位置一起标明。

PtSe₂/PtTe₂ 异质双层的部分电荷密度把选定能量窗口、三维等值面和层向密度曲线放到同一组图内。若某个带边态主要分布在一侧层，界面电荷分离、带偏移和层间耦合就有了空间依据；若态密度跨越两层，界面杂化会在部分密度中留下连续的空间分布。

电荷密度能回答哪些材料问题？

在成键分析中，电荷密度常用于观察电子是否集中在两个原子之间、是否偏向某个配位原子，或者是否沿界面形成连续分布。键区电子分布可支撑共价成分、离子成分和界面极化的讨论，但键强还要由键长、COHP/COOP、振动频率或能量差校正。

ELF 使用电子密度及其梯度构造局域化指标，数值接近 1 的区域通常表示电子更局域，接近 0.5 的区域更接近均匀电子气。电荷密度回答电子多不多，ELF 回答电子是否局域，它们都在实空间显示电子结构，对应的物理量不同。

Ti₃C₂/Cu 复合界面中的 ELF 截面把局域化区域放在原子层之间比较。若界面处出现连续的局域化带，通常对应更强的方向性相互作用；若 ELF 只在各自原子附近增强，界面成键贡献较弱。这类判断应和 PDOS 或 COHP 对同一原子对的贡献同步核对，避免把局域色块误写成键能数值。

在离子插入、扩散和缺陷筛选中，电荷密度还可作为空间场参与位点搜索。局域势给出带电离子所感受到的静电环境，电荷密度给出电子云和分子骨架的空间占据。对 Li⁺、Na⁺ 等插入离子而言，低排斥区域和可配位原子环境共同限定候选位点，随后再用 DFT 弛豫比较总能。

有机晶体离子插入位点搜索把局域势、电荷密度和几何空隙放到同一流程中。电荷密度在这里并非装饰图；离子靠近哪些官能团、避开哪些高电子密度区域，都要参考这一路径上的空间输入。电池材料、离子导体和有机电极里的候选位点筛选，经常要把密度场和结构弛豫能量共同使用。

在催化表面和缺陷体系中，电荷密度可指出电子在吸附物、金属中心、空位邻域或载体之间的重排位置。它能支撑“界面极化增强”“缺陷态局域”“吸附物轨道参与相互作用”等判断范围；若要谈反应活化，仍要回到 E_ads、ΔG、NEB 能垒和 PDOS/COHP。

电荷密度判断应限定哪些计算条件？

密度图能否比较，取决于计算口径是否一致。赝势决定核心电子如何处理，网格密度决定零通量面和等值面的分辨率，泛函与 U 值会改变局域 d/f 电子。同一系列模型应使用同一泛函、赝势、ENCUT、FFT 网格和自旋设置，否则密度差可能来自数值口径，并非材料结构变化。

等值面大小、切片平面和颜色符号会改变读者看到的形状。黄色、蓝色、红色或绿色没有通用含义，必须以原图图注和绘图设置为准。等值面过低会把远处弱信号连成片，等值面过高会漏掉界面弱极化，同一组模型最好固定 isovalue 和色标范围。

原子电荷数来自分区或投影方法，并非连续密度图自动给出的项目。Bader、Hirshfeld、Mulliken 和 DDEC 使用不同的空间分割或基组投影口径，得到的绝对值常有差异。电荷密度图负责空间位置，原子电荷负责积分结果，二者在文章中应分工明确。

BaTiO₃ 掺杂模型中的 Mulliken 电荷分布把连续电子密度投影到原子或轨道基函数上。它适合比较同一方法下不同掺杂原子周围的电荷偏移，但不宜把数值直接写成形式价态。电荷数值的符号、参考电子数和基组口径，决定了它能支撑哪一种材料判断。

把电荷密度用于论文图谱解释时，可核对六项内容：密度文件类型、参考片段几何、网格和赝势、颜色正负号、等值面或色标范围、选定能量窗口。对于 CO₂ 吸附、H^* 中间体、氧空位、异质结和二维材料表面，同一套检查能区分总电子云、界面重排、局域态分布和原子电荷积分。