差分电荷密度matlab,CASTEP差分电荷密度操作流程
原标题：CASTEP 差分电荷密度操作流程
本⽰例以 CO 吸附在 Pt(111) 表⾯上为例
1、对吸附的结构进⾏结构优化，如图1所⽰：
图1 「优化后的结构」
2、如果要研究吸附分⼦ CO 和表⾯ Pt 的之间由于吸附引起的电荷转移，则要选中 CO 分⼦并定义 Set 。
【操作】：主标题栏 Edit-Edit ts ，在 Define New Set 对话框的 Name 参数设置处，输⼊的 t 名称要包含 DensityDifference 字段，如图2所⽰：
图2 「定义 t 」
包含 t 的结构如图3所⽰：
图3 「定义 t 的结构」
3、提交计算任务
【操作】：Task 中选择 Energy ，Properties 中选择Electron density difference ，如图4所⽰，提交计算。
图4 「参数设置」
图4中的参数说明：
Atomic densities ：软件将计算当前结构中的状态中的每个原⼦与其孤⽴存在时的电⼦密度差异。
Sets of atoms ：软件将计算各个t在当前结构中与t包含的分⼦孤⽴状态时的电⼦密度差异。
4、分析差分电荷密度
使⽤ CASTEP ⾃带的分析⼯具分析差分电荷密度，如图5所⽰。选择 Electron density difference ，当直接点击 import 时，引⼊基于t 的差分电荷密度。当勾选 U atomic densities 时，引⼊基于 atom 的差分电荷密度。
图5 「 CASTEP 分析⾯板」
使⽤如下⼯具条调整差分电荷密度的显⽰。
其中 Volumetric Selection 选项，可呈现结构中包含了哪类差分电荷密度。
通过⼯具条中的 Color Maps 调节差分电荷密度的显⽰
调整前的结构如图6所⽰：
图6 「调整前的结构和参数设置⾯板」
调整后的显⽰如图7所⽰：
图7 「差分电荷密度」
从差分电荷密度图的计算结果可以看出，C原⼦和与C原⼦结合的Pt原⼦失去了电⼦，周围差分电荷密度显⽰为蓝⾊(蓝⾊对应负值，失去电⼦)，O原⼦得到了电⼦(红⾊对应正值，得到电⼦)。
图7的显⽰对应的 Color Maps 对话框的参数设置如图8所⽰。不显⽰的⾊段是通过⿏标点选/划选删除的，这些⾊段对应的是差分电荷密度接近0的区域。
图8 「图7显⽰所需的参数设置」返回搜狐，查看更多
责任编辑：
