计算化学公社

标题: 求助CP2K做AIMD研究高浓度水溶液体系的流程问题 [打印本页]

作者
Author:
花非花    时间: 5 day ago
标题: 求助CP2K做AIMD研究高浓度水溶液体系的流程问题
各位大神,新手刚开始学习cp2k的AIMD的模拟,摸索了一段时间,想学习AIMD和实验测试的进行结合,目的是得到平衡后的溶剂化结果,最多需要扩散系数。我先MS搭建模拟的高浓度(3~20 mol/kg)的不同电解质(如LiNO3,NaNO3,KNO3)体系的盒子(由于密度相差很大,搭建的时候考虑实验的密度,盒子一般大于15个埃,原子数500左右),然后选择COMPASIII力场进行MD得到初始结构: 运用此结构可以直接到CP2K固定晶胞选择OPT优化,得到合理的结构:1。直接进行NVT系综的模拟吗?还是2:需要先用经典的NPT预平衡后再NVT计算? 谢谢!~



作者
Author:
sobereva    时间: 5 day ago
如果怀疑力场下MD得到的密度不合理,CP2K做NVT MD前显然得先经过NPT的MD

更多信息参考
谈谈做分子动力学模拟什么时候适合用NPT和NVT系综
http://sobereva.com/761http://bbs.keinsci.com/thread-57835-1-1.html
作者
Author:
花非花    时间: 5 day ago
谢谢SOB大神!您的帖子我拜读过了。我的研究方向和主要内容还是在实验,AIMD只是补充。这样我是不是偷懒一点,采取MS去做NPT稳定晶胞,然后去CP2K-固定晶胞OPT----低温预热 ------300K正式开跑NVT?谢谢!~!
作者
Author:
sobereva    时间: 4 day ago
花非花 发表于 2026-4-8 17:47
谢谢SOB大神!您的帖子我拜读过了。我的研究方向和主要内容还是在实验,AIMD只是补充。这样我是不是偷懒一 ...

取决于审稿人是否是内行
内行审稿人可能会怀疑M$得到的稳定晶胞尺寸是否合理,导致CP2K的NVT其实是在比如很高压下跑的




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