计算化学公社

标题: 通过Multiwfn绘制电子密度差图时，电子密度区域偏离晶体晶体结构 [打印本页]

作者
Author: 故1996 时间: 2026-3-21 10:22
标题: 通过Multiwfn绘制电子密度差图时，电子密度区域偏离晶体晶体结构
老师你好，我所研究的体系是缓蚀剂在金属氧化物（Fe2O3)表面上的吸附行为。
我通过构建缓蚀剂加slabFe2O3模型，通过cp2k计算了两者的吸附能，然后想绘制电子密度差图。
但是绘制时发现，电子密度差图中，电子密度相对晶格有一定位移。并且其中的平面平均密度差（plane-averaged EDD）曲线，也是在界面处是锯齿状线条。
我怀疑是在计算单点能过程中k点取值过小导致，还有就是体系中的cutoff取值是600，60。因为体系中含有Na离子。
不知该如何改进。
还有个问题如果发生化学吸附等值面差值是不是应该选的交大一些（原本看博文，物理吸附是0.0005）
cp2k计算中是采用了全电子基组pob-TZVP，尝试多次只有这个可以收敛，开启OT，采用色散校正。XY周期性，自旋磁矩按照讲义设置，4和-4。

作者
Author: sobereva 时间: 2026-3-21 22:14
Multiwfn屏幕上有明确提示，无法显示非正交晶胞的等值面。需要导出cube文件用VMD、VESTA等第三方程序看等值面图。先确保等值面图正确再说其它的
并且要确保求差的时候原子坐标相互对应（更具体来说，被冻结的基底原子的坐标在两次计算后得到的给Multiwfn的cube或molden文件中完全一样）

欢迎光临计算化学公社 (http://bbs.keinsci.com/)