使用CP2K跑Al-Mg-Cu合金的MD预平衡过程，体系温度和压强都一直增加怎么解决(已解决)

xyl · 发表于 Post on 2023-6-30 12:28:04

本帖最后由 xyl 于 2023-7-3 10:18 编辑

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

sobereva · 发表于 Post on 2023-7-1 07:44:46

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

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

xyl · 发表于 Post on 2023-7-1 09:55:31

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

sobereva老师用的是“MSST方法模拟石墨在冲击压缩下变成金刚石”例子，scf是收敛的

sobereva · 发表于 Post on 2023-7-2 11:18:28

用popt版跑，当前速度明显慢得要死
结合轨迹判断情况

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

xyl · 发表于 Post on 2023-7-2 11:27:23

好的老师

xyl · 发表于 Post on 2023-7-28 09:33:35

sobereva 发表于 2023-7-2 11:18
用popt版跑，当前速度明显慢得要死
结合轨迹判断情况

sob老师你好我想请问一下还是这个例子。在MSST过程中我为了把压缩时间控制在1ps以内跳小了晶胞虚拟质量和把步长增大到0.5fs。现在导致晶体不是呈现压缩过程而是一个膨胀过程是什么原因。

sobereva · 发表于 Post on 2023-7-28 14:04:53

xyl 发表于 2023-7-28 09:33
sob老师你好我想请问一下还是这个例子。在MSST过程中我为了把压缩时间控制在1ps以内跳小了晶胞虚拟质 ...

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

xyl · 发表于 Post on 2023-7-28 14:50:18

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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