计算化学公社

标题: cp2k电势分布后处理怎么做 [打印本页]

作者
Author: zzk 时间: 2022-5-15 20:16
标题: cp2k电势分布后处理怎么做
本帖最后由 zzk 于 2022-5-15 20:35 编辑

请问cp2k得到结构后怎么进行电势分布的后处理得到如下图的结果，multiwfn可以做吗

作者
Author: zzk 时间: 2022-5-15 20:27
类似这种结果

作者
Author: chengsq 时间: 2022-5-15 20:55
我也遇到了同样的问题。。。求老师解答！

作者
Author: sobereva 时间: 2022-5-15 21:37

zzk 发表于 2022-5-15 20:27
类似这种结果

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。
原文贴出来

作者
Author: zzk 时间: 2022-5-15 21:50
本帖最后由 zzk 于 2022-5-15 22:56 编辑

这个是论文的方法部分，先通过gmx和cp2k得到结构，然后分析dos和电势这是论文链接

作者
Author: sobereva 时间: 2022-5-15 23:22
6式里的rho大抵对应平面平均的电子密度，用下文的方法就可以得到rho在不同z位置的值（文中的密度差改为电子密度）
使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线
http://sobereva.com/638（http://bbs.keinsci.com/thread-28225-1-1.html）

然后自己写个程序做梯形积分就完了

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