计算化学公社

标题: 关于QE做第一性原理MD收敛的问题 [打印本页]

作者
Author: 茶味果冻 时间: 2015-8-25 09:27
标题: 关于QE做第一性原理MD收敛的问题
最近在用QE中的cp.x做cpmd的计算，按照教程文件中三步走，第一步电子基态，收敛标准的话看out文件中ekinc的值即可，这个好判断。紧接着是第二步，即原子结构的弛豫（relax the system,check the forces on atoms），不然的话d第三步MD计算原子会全部跑散。那问题就来了啊，这里看判定标准的话教程文件中已经说明了怎么去估算自己建立的结构需不需要进行这一步的弛豫，但是没有说明做弛豫的话，结果out文件去看哪个值来判定我的结构到底有没有达到力最小，下图贴出了out文件中唯一找到forces acting on atoms(au)，这里求教啊，out文件中判定的标准是什么？还有下图中为啥每个原子会给出3个值，这三个值是指xyz三个方向上的值么，判定的话是不是每个值都需要满足in文件中设置的forc conv thr.值？求解答，谢谢了。
@万里云

作者
Author: 茶味果冻 时间: 2015-8-25 09:29
@jiahh

作者
Author: 茶味果冻 时间: 2015-8-25 09:51
@万里云

作者
Author: 万里云 时间: 2015-8-25 10:39

这个……我对分子动力学完全不懂啊。

作者
Author: 茶味果冻 时间: 2015-8-25 13:52

万里云发表于 2015-8-25 10:39
这个……我对分子动力学完全不懂啊。

好吧，同样感谢！我做静态计算的话也用QE，刚学一年多，以后有问题希望可以多多交流学习

作者
Author: yoya 时间: 2019-11-23 11:56
教程文件中讲到relax the system分两种情况，一种是你自己手动搭建的结构，第二种是用之前AIMD的方法的得到的体系原子位置，给出了两种解决方案。但我之前按照教程跑的时候，原子全部散掉了。所以后面我直接跳过了第二步再算，就没有问题了。

作者
Author: YANGQING 时间: 2022-3-25 09:35
大佬你好，我最近也在做cpmd,但是在结构优化第二步离子能量最小化时能量一致波动达不到收敛，我是否可以用pw.x进行结构优化，然后cp.x中忽略结构优化直接跑md呢？还有cp.x中ensemble-dft是只针对金属系统结构优化过程还是md过程也要这么设置呢？因为对于这条命令一直没懂，有报错，所以想跳过它。

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