计算化学公社

标题: Cp2k 做NPT动力学模拟的几个问题 [打印本页]

作者
Author: wanlichuan 时间: 2024-3-6 11:11
标题: Cp2k 做NPT动力学模拟的几个问题
本帖最后由 wanlichuan 于 2024-3-6 11:13 编辑

大家好，我想做一个晶体材料的NPT系综动力学，目标温度300K，本来设定的是跑1000步，500fs，但是跑到802步就出错了。

前面802步的信息都保留了。

以下是温度的变化情况：

以下的超胞体积的变化情况：

inp文件（坐标省略）和Out文件（中间步骤重复内容删掉）见附件。

我的问题：

1）out文件中没有出错信息，怎么知道哪儿的问题呢？

2）为什么温度升到232.6K就基本不再上升了呢？我设的是300K，就算步数没跑够（差200步），升温时间常数设的有点大（300 fs，已经跑了400 fs），但是升温趋势不应该停止啊。

3）为什么温度基本不再上升的时候，体积却在一直增加呢？

谢谢。

作者
Author: 泰然汪汪汪 时间: 2024-3-6 12:38
1. 考虑你硬件的问题，再跑一遍看能不能复现错误。
2. 3. 结构AIMD前优化了吗？SCF都收敛了吗？估计结构跑散掉了已经

作者
Author: wanlichuan 时间: 2024-3-6 13:05
本帖最后由 wanlichuan 于 2024-3-6 13:32 编辑

谢谢回复。
1）我跑过几个类似的结构和相同的参数设定，都有类似的问题，应该不是偶然。区别在于有些跑完了，有些像此例一样，差一点没跑完，没跑完的也是没有错误信息；无论跑完与否，温度都无法达到设定温度，并且体积一直增加。
2）由于来源于实验晶体结构，并且是跑动力学，就没事先优化。最后没有跑散，只是分子间距离增大了，有规律的增大。
再次感谢。

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