计算化学公社

标题: cp2k跑分子动力学模拟时，设置的目标温度出现突然下降 [打印本页]

作者
Author: ljb874722957 时间: 2025-3-4 21:46
标题: cp2k跑分子动力学模拟时，设置的目标温度出现突然下降
本人cp2k初学者，我最近在用NVT系综对结构进行分子动力学模拟时，初步猜想是设置一个温度找到结构的相变温度，我跑到1600K时结构还是稳定的，因此我设置一个较大的温度3000K，但我发现设置3000K时温度会突降，降低到2000K左右，这是怎么回事？

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-3-4 22:04
本帖最后由 Uus/pMeC6H4-/キ于 2025-3-4 22:05 编辑

不知道体系结构长什么样、如何获得，没法具体评论。另外开OT却不用OUTER_SCF、仍把SCF迭代次数上限设到200并不合适，需按http://sobereva.com/665第3节的说明修改。热浴推荐用CSVR。创建CP2K输入文件（以及入门CP2K计算）的途径在http://sobereva.com/587有介绍。

另外确认一下，这个20多埃尺寸的盒子在3000 K下跑NVT动力学3 ps真能看到相变么？这要求没有高温分解的副反应、相变发生足够快、不改变晶胞体积且不产生超过20埃的长程有序结构。

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-3-4 22:51
本帖最后由 Uus/pMeC6H4-/キ于 2025-3-4 22:53 编辑

ljb874722957 发表于 2025-3-4 22:26
老师您好，这是我的结构文件，麻烦您看一下有什么问题没有，我用NVT的原因主要是我主用的软件是vasp，然 ...

原来是经典的g-C3N4单层啊，单看结构也没啥问题。1楼那个分子动力学初始阶段温度快速降低，应该是因为体系中原子受力还比较小（所以我问了怎么获得，就是估计做过几何优化），速度不快，热浴控温的效果也一时没跟上来。延长模拟时间、热浴充分作用了，温度自然就回升了（参考http://sobereva.com/627所说的平衡判断），如果观察每一步SCF均收敛、轨迹也正常的话就不用担心太多。注意CP2K暂时没有线性升温/控制速度退火的功能，其温度设置既用于初始化速度也作为热浴的参考目标。

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