计算化学公社

标题: cp2k进行AIMD如何处理降温过程 [打印本页]

作者
Author: nusiew 时间: 2021-6-30 19:23
标题: cp2k进行AIMD如何处理降温过程
各位老师好，
本人cp2k新手，做过一段时间的硅酸盐熔体经典分子动力学模拟。现在用cp2k做AIMD研究1823K条件下硅酸盐熔体结构。以前基于经典力场用lammps研究硅酸盐熔体结构时，由于建立的初始构型是原子随机均匀分布的构型，会首先将初始构型在5000K跑若干ps以消除初始构型记忆，再以恒定冷速经降温到1823K, 然后再跑若干皮秒使体系达到平衡。由于缺少某些体系经典力场，目前转向AIMD。
现在用cp2k研究硅酸盐熔体结构是否还需要在5000K的条件下弛豫以消除初始结构记忆呢？如果需要的话，后续从5000降温到1823K过程怎么来处理呢？我想的是有两个方案：
1. 使用cp2k自带关键词ANNEALING进行处理。
2. 分别在5000，4500，4000，3500，3000，2500，2000，1823K条件下依次进行AIMD模拟，从而实现从5000降温到1823这个过程。
请问老师们如何看待以上两个方案，有没有其他建议？
最后，1823k条件下硅酸盐体系已经处于熔体状态，能不能省略降温过程，对初始随机均匀分布构型直接在1823K条件下进行AIMD呢？
谢谢各位老师，望各位老师不吝赐教！

作者
Author: sobereva 时间: 2021-7-1 04:52
也需要消除初始结构记忆

用ANNEALING就可以，也可以用https://manual.cp2k.org/trunk/CP ... MPERATURE_ANNEALING

如果你要直接在1823K下跑，原理上需要跑更长时间才能消除初始结构记忆

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