计算化学公社

标题: 求助：VASP跑AIMD怎么记录每一步的bader电荷 [打印本页]

作者
Author: we000000 时间: 2021-11-3 10:30
标题: 求助：VASP跑AIMD怎么记录每一步的bader电荷
本帖最后由 we000000 于 2021-11-3 10:32 编辑

看到了一篇用VASP跑AIMD的文章，如图： (, 下载次数 Times of downloads: 18) 请问这种每一步的电荷是怎么计算的。我自己在INCAR中开启LCHARG = .TRUE. 和 LAECHG = .TRUE.跑完5000步之后用bader脚本处理的得到的ACF.dat文件里面只有一步的。他是这每一步都自己用脚本循环提交bader任务跑过吗.... 感觉这样文件目录都爆炸了，感谢大家解答

作者
Author: 卡开发发 时间: 2021-11-3 11:08
我不确定文献用的是否为Bader电荷，如果基于WS划分我以前给过一些办法可以尝试（http://bbs.keinsci.com/thread-5920-1-1.html）
对Bader电荷，理论上能做，但咋改程序需要另外想办法，实际可能也不难。思路上应该是做一步离子步分析一次电荷，然后更新结构再自洽刷新当前的chgcar。

作者
Author: we000000 时间: 2021-11-3 11:21

卡开发发发表于 2021-11-3 11:08
我不确定文献用的是否为Bader电荷，如果基于WS划分我以前给过一些办法可以尝试（http://bbs.keinsci.com/th ...

十分感谢老师

作者
Author: sobereva 时间: 2021-11-3 23:09
如果这是分子，用ORCA+Multiwfn获得原子电荷随时间的变化是最好的做法，最简单而且速度最快

ORCA做AIMD模拟很容易
使用ORCA做从头算动力学(AIMD)的简单例子
http://sobereva.com/576（http://bbs.keinsci.com/thread-20800-1-1.html）

此文里有获得每一帧自旋布居现成的脚本，改一一行就能得到每帧的原子电荷
详谈Multiwfn的命令行方式运行和批量运行的方法
http://sobereva.com/612（http://bbs.keinsci.com/thread-24929-1-1.html）

作者
Author: we000000 时间: 2021-11-4 10:29

sobereva 发表于 2021-11-3 23:09
如果这是分子，用ORCA+Multiwfn获得原子电荷随时间的变化是最好的做法，最简单而且速度最快

ORCA做AIMD ...

谢谢sober老师

作者
Author: reedfy 时间: 2021-11-4 10:45
你好，这篇文章题目能提供下嘛，也想学习学习

作者
Author: we000000 时间: 2021-11-4 10:57

reedfy 发表于 2021-11-4 10:45
你好，这篇文章题目能提供下嘛，也想学习学习

J. Phys. Chem. C 2015, 119, 26828−26839

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