用MS的Dmol3模块计算：SnO2 SnS SnS2 CuS Cu2S ZnS PbS CdS FeS2吸附In3+可行性？

1814

老师您好，我测了一下SnO2 SnS SnS2 CuS Cu2S ZnS PbS CdS FeS2九种材料的XRD，想根据最强的信号对应的晶面，使用Materials Studio的Dmol3模块对In3+离子的进行吸附能大小的计算，另外再计算一下差分电荷密度、Bader电荷、态密度，比较一下九种材料的区别。请问这种方法可行吗，需要注意一些什么吗？另外如果可行，请问需要如何改进呀？要是不可行，更应该如何改进呢？

kaikai

2020以前的版本里面Dmol3模块应该是没办法计算Bader电荷，VASP，QE，CP2K是可以的。另外Dmol3使用的数值型函数，大部分人使用的VSAP这种平面波基组软件，看你的问题你应该是新手，你要是想理论计算结果与实验数值相吻合，需要相当的功底。

卡开发发

kaikai 发表于 2025-10-12 10:08
2020以前的版本里面Dmol3模块应该是没办法计算Bader电荷，VASP，QE，CP2K是可以的。另外Dmol3使用的数值型 ...

真正存在问题的是In(3+)离子的吸附，带电表面吸附的问题可能一直是诸多程序难用处理的问题。

Bader在MS没办法直接算，但是因为DMol3能够提供基于标量相对论的全电子计算，可以按照如下步骤实现：
1、用赝势和全电子分别进行单点自洽计算得到的电子密度的grd，计算需要足够小的网格间距。
2、写脚本将赝势和全电子分别得到的grd文件转换为CHGCAR和CHGCAR_sum。
3、和vasp处理Bader电荷方式类似，bader CHGCAR -ref CHGCAR_sum。
其实CHGCAR的格式更换为cube格式也是可以的，其中可能需要花时间去了解grd和CHGCAR/cube的格式以及单位转换等诸多细节。

当然了，Bader电荷本身有些“缺陷”，论坛有些讨论这里不详细展开。