MS计算差分电荷密度：CASTEP中参考态选取与界面电荷重分布的实战分析

MS计算差分电荷密度这个操作，在CASTEP里看起来只是几步——先算完整体系，再算各部分单独的电荷密度，最后做差值。但真正做过界面体系差分电荷分析的人都知道，差分电荷密度的可靠性取决于一个看似简单的选择：参考态怎么取。这个选择做得不对，后面所有的等值面图、积分分析、电荷转移量计算，全部建立在错误的根基上。

项目背景：CO₂电催化还原中的界面电荷转移

团队处理过一个Cu单原子催化剂上的CO₂电催化还原项目。催化活性位点是一个Cu原子锚定在石墨烯基底上的N掺杂位（Cu-N₄构型），客户需要知道CO₂吸附后电荷从催化剂向CO₂分子转移了多少、转移集中在哪些原子轨道上——这直接关系到CO₂活化程度和后续还原路径的选择。

差分电荷密度的定义是 Δρ = ρ(total) – ρ(substrate) – ρ(CO₂)，其中ρ(total)是吸附体系的总电荷密度，ρ(substrate)和ρ(CO₂)是各自在”参考态”下的电荷密度。关键问题就在于：参考态中的原子位置应该是什么？

参考态选取的两种策略及其物理含义

第一种做法是把参考态中各部分的原子放在和吸附体系完全相同的坐标上，只是去掉对方部分。这叫”冻结几何”参考态——substrate保持吸附后的弛豫构型，CO₂也保持吸附后弯曲的构型。这种做法的差分电荷密度反映的是”电子在吸附形成后如何重新分配”，不含几何弛豫的贡献。

第二种做法是让参考态中各部分分别优化到各自的平衡构型——substrate优化到无吸附时的构型，CO₂优化到自由分子的直线构型。这种做法的差分电荷密度同时包含了电子重分布和几何弛豫的影响。

两种策略各有适用的场景。如果目标是理解”吸附后电子结构的改变对催化活性的影响”，冻结几何参考态更纯粹——它只剥离了电子效应。但如果需要评估”整个吸附过程的电荷转移总量”，包含几何弛豫的参考态给出的数字更接近实验可观测的量。

团队在这个项目上认定冻结几何参考态更适合——因为客户关注的是”CO₂在Cu位点上的活化程度”，这本质上是一个电子结构问题，几何弛豫的贡献应该单独分析而非混入差分电荷中。

CASTEP中的实操步骤与密度网格精度

MS计算差分电荷密度在CASTEP中的操作流程：先完成吸附体系的自洽计算，保存.chg文件；然后分别计算substrate和CO₂在冻结几何下的自洽电荷密度（用相同截断能和k点网格）；最后在Materials Studio的Analysis模块中选择”Charge Density Difference”，指定两个参考态文件。

密度网格精度是容易被忽视的参数。CASTEP默认的fine grid cutoff为截断能的1.5倍，对于差分电荷密度来说这个精度可能不够——因为差分操作会在界面区域产生精细的电荷重分布特征（偶极层、局域极化），这些特征的空间尺度可能只有0.5-1 Å，需要更密的网格才能解析。团队将fine grid cutoff提高到截断能的2.5倍（从默认的750 eV提升到1250 eV），差分电荷在Cu-N₄位点附近的等值面细节明显更丰富——0.01 e/ų级别的电荷聚集区从模糊变得清晰。

这个精度提升的代价是计算量增加约60%，但对于差分电荷分析来说，这60%的投入是值得的——粗糙网格给出的差分电荷等值面在界面区域会出现虚假的平滑化，把物理上的局域电荷转移抹成大面积的模糊分布，分析结论的可信度直接打折。

可视化阈值：0.005 e/ų不是随意选的

差分电荷密度的等值面图需要选择阈值（isovalue）。在CASTEP的可视化模块中，默认阈值通常是0.01 e/ų，但团队在Cu-CO₂体系上发现，0.01的阈值漏掉了相当多的弱电荷转移特征——特别是CO₂的O端与Cu之间的π反馈作用产生的电荷富集区，其峰值只有约0.008 e/ų。

团队最终选取了0.005 e/ų作为主要展示阈值，0.02 e/ų作为强电荷转移区的辅助阈值。0.005的选择有物理依据：它大约是自由原子电荷密度在核区1 Å处的1/1000，低于这个阈值的信号更可能是数值噪声而非物理效应。

从差分电荷到电荷转移量的定量积分

差分电荷等值面图直观，但定量分析需要对差分电荷做空间积分。团队沿垂直于石墨烯平面的方向对Δρ做了面积分，得到从基底到CO₂分子方向的电荷转移分布曲线。积分结果显示，CO₂吸附后从Cu位点向CO₂方向转移了约0.32 e的电子——其中0.21 e集中在Cu的d轨道向CO₂的π*轨道的反馈作用上，0.11 e来自基底N原子向Cu的再分配。

这个0.32 e的转移量与Bader电荷分析的0.28 e差异约14%——差距不会说谎，两种方法对电荷边界的定义不同导致了这个偏差。Bader分析基于电荷密度的零通量面划分原子区域，差分积分则基于空间体积的直接求和。在界面体系中，两种方法的一致性在15%以内就足以支撑定性结论。

回过头看这个项目，MS计算差分电荷密度最关键的决策出现在第一步——参考态的选取方式。冻结几何还是自由优化，这个选择不是”哪种更方便”的问题，是”你想分析什么物理效应”的问题。密度网格精度和可视化阈值也有物理依据可循，不是随意调参数凑视觉效果。

更多内容请访问 https://www.keyanxueshu.com/