计算化学公社

标题: 使用CP2K跑Al-Mg-Cu合金的MD预平衡过程,体系温度和压强都一直增加怎么解决(已解决) [打印本页]

作者
Author:
xyl    时间: 2023-6-30 12:28
标题: 使用CP2K跑Al-Mg-Cu合金的MD预平衡过程,体系温度和压强都一直增加怎么解决(已解决)
本帖最后由 xyl 于 2023-7-3 10:18 编辑

求助使用CP2K跑Al-Mg-Cu合金的MD预平衡过程,采用的NPT_I,但是体系温度和压强上升的很快,且上升到几千K,是否合理,不合理的话需要怎么设置条件。我参照冲击石墨变成金刚石的算例,设置的inp文件。

作者
Author:
sobereva    时间: 2023-7-1 07:44
不要用“冲击石墨变成金刚石的算例”来描述,别人都不知道你准确说的是什么东西,仔细看下文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620http://bbs.keinsci.com/thread-25787-1-1.html
如果是指北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)里第一性原理动力学部分我提供的“用MSST方法模拟石墨在冲击压缩下变成金刚石”例子,交代清楚

没输出文件没法判断。注意看SCF收敛情况,并且观察轨迹。


作者
Author:
xyl    时间: 2023-7-1 09:55
sobereva 发表于 2023-7-1 07:44
不要用“冲击石墨变成金刚石的算例”来描述,别人都不知道你准确说的是什么东西,仔细看下文
在网上求助计 ...

sobereva老师用的是“MSST方法模拟石墨在冲击压缩下变成金刚石”例子,scf是收敛的
作者
Author:
sobereva    时间: 2023-7-2 11:18
用popt版跑,当前速度明显慢得要死
结合轨迹判断情况

如果怀疑是热浴问题,减小timecon

作者
Author:
xyl    时间: 2023-7-2 11:27
好的老师
作者
Author:
xyl    时间: 2023-7-28 09:33
sobereva 发表于 2023-7-2 11:18
用popt版跑,当前速度明显慢得要死
结合轨迹判断情况

sob老师你好 我想请问一下还是这个例子。  在MSST过程中我为了把压缩时间控制在1ps以内跳小了晶胞虚拟质量和把步长增大到0.5fs。现在导致晶体不是呈现压缩过程而是一个膨胀过程是什么原因。
作者
Author:
sobereva    时间: 2023-7-28 14:04
xyl 发表于 2023-7-28 09:33
sob老师你好 我想请问一下还是这个例子。  在MSST过程中我为了把压缩时间控制在1ps以内跳小了晶胞虚拟质 ...

至少得在比较短的晶胞方向上考虑k点,否则就算不跑MSST,优化晶胞结构也很烂
另外,Mg用DZVP-MOLOPT-SR-GTH对CUTOFF要求很高,300根本不够用。最起码用800。如果嫌太贵,Mg也可以用诸如pob-DZVP-rev2

作者
Author:
xyl    时间: 2023-7-28 14:50
sobereva 发表于 2023-7-28 14:04
至少得在比较短的晶胞方向上考虑k点,否则就算不跑MSST,优化晶胞结构也很烂
另外,Mg用DZVP-MOLOPT-SR- ...

谢谢sob老师答疑,我去调整一下




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3