求助cp2k计算NPT_I控制压强压缩水盒子以及带真空层的水盒子动力学计算问题

liurui

各位老师，学生用cp2k在M06-2X with ADMM/DZVP-MOLOPT-SR-GTH-q1层次计算30个水以及60个水，密度1g/cm3构建的无定形盒子，用的csvr方法控制热浴，TIMECON =200，之前看博主老师在其他帖子回复NPT不用Nose-Hoover therm ;但是，学生感觉奇怪的点在于学生发现out文件里输出的pressure并没有维持在输入文件里设置的1bar，输入和输出文件如附件所示,请各位老师帮学生看看是不是哪里整的不对，谢谢各位老师！非常感谢！

sobereva

纯水体系跑MD没事甭用杂化泛函，尤其是M06-2X，又贵又没好处。一般目的PBE-D3(BJ)足矣
又是用了SCREEN_ON_INITIAL_P T而不读wfn文件，典型的致命错误，好好看
CP2K做杂化泛函计算的关键要点和简单例子
http://sobereva.com/690（http://bbs.keinsci.com/thread-41087-1-1.html）
绝对不要在没杂化泛函使用常识的情况下用杂化泛函

“NPT不用Nose-Hoover therm ”完全语义不明

带那么大真空区作甚。倘若就是要模拟气水界面显然应该用NVT

liurui

sobereva 发表于 2024-12-11 10:56
纯水体系跑MD没事甭用杂化泛函，尤其是M06-2X，又贵又没好处。一般目的PBE-D3(BJ)足矣
又是用了SCREEN_ON_ ...

老师，学生NVT和NPT都要跑, canonical ensemble [NVT] and isothermal-isobaric ensemble [NPT], temperature [270, 300, 330, 400, 500, 600 K], and pressure [0.2, 1, 5, 10 bar]),想与前人发表的一篇工作做的类似，前人的这篇工作用vasp进行的计算；老师，您之前在回复帖子回复变温过程用Nose-Hoover热浴是众所周知的禁忌，会出现震荡我建议一律用CSVR热浴。Nose-Hoover没必要考虑，想着这里用CSVR，是不是学生这里用Nose-Hoover也是合理的哇，谢谢老师！

Uus/pMeC6H4-/キ

怎么这帖也有种http://bbs.keinsci.com/thread-50088-1-1.html的既视感，也是没有太多CP2K做水溶液AIMD的经验却要急着复现文献用杂化泛函跑MD的工作……我就不重复一遍那楼里面我的两条评论了，请自行移步好好阅读。

回一楼“输出的pressure并没有维持在输入文件里设置的1bar”，看下述社长博文，特别是有关压力的部分。
谈谈怎么判断分子动力学模拟是否达到了平衡
http://sobereva.com/627（http://bbs.keinsci.com/thread-27122-1-1.html）

回三楼“是不是学生这里用Nose-Hoover也是合理的”，这取决于目的是要复现文献还是要自己产生数据，注意文献的设置未必合理。

liurui

Uus/pMeC6H4-/キ 发表于 2024-12-11 21:46
怎么这帖也有种http://bbs.keinsci.com/thread-50088-1-1.html的既视感，也是没有太多CP2K做水溶液AIMD的经 ...

收到，谢谢老师，我马上仔细阅读理解一下子

sobereva

liurui 发表于 2024-12-11 16:09
老师，学生NVT和NPT都要跑, canonical ensemble [NVT] and isothermal-isobaric ensemble [NPT], tempera ...

CP2K里若无特殊情况一律别用Hose-Hoover，没事用它干嘛，明明CSVR更普适

liurui

sobereva 发表于 2024-12-11 23:20
CP2K里若无特殊情况一律别用Hose-Hoover，没事用它干嘛，明明CSVR更普适

收到，谢谢老师