计算化学公社

标题: 使用CP2K做NPT模拟时遇到slab严重变形问题 [打印本页]

作者
Author:
Stardust0831    时间: 2024-4-5 22:10
标题: 使用CP2K做NPT模拟时遇到slab严重变形问题
本帖最后由 Stardust0831 于 2024-4-6 00:00 编辑

老师们好,我在用cp2k跑NPT的aimd的时候发现z轴方向变得异常的大,体系在扩成超胞前用PBEsol优化过晶格常数,并经过2ps的NVT预平衡,压控参数为1.01325,时间常数最低尝试过50,也尝试提高CUTOFF到600,其余参数都是严格按照第三届第一性原理计算培训班中的TiO2与水界面的模拟部分设置,全程SCF都能正常收敛。现附上我正常结束的NVT预平衡的文件以及遇到异常变形的NPT部分文件,望各位老师不吝赐教。
4月6日更新:尝试去掉偶极矫正,异常变形问题仍然存在。
作者
Author:
sobereva    时间: 2024-4-6 00:17
你这个体系也没有水,本来就不应该用NPT。一直用NVT就完了
作者
Author:
Stardust0831    时间: 2024-4-6 09:23
sobereva 发表于 2024-4-6 00:17
你这个体系也没有水,本来就不应该用NPT。一直用NVT就完了

老师,我做NVT的时候材料分子只在真空层自己振动而不和晶体表面发生反应。请问既然NPT可以调节水层密度,能不能也调节材料分子层密度、帮助它接近晶体表面?如果继续做NVT,会面临材料分子不靠近晶体表面的问题,是否建议加大分子数量来解决?有其他更好的解决方法吗?
作者
Author:
sobereva    时间: 2024-4-7 05:22
Stardust0831 发表于 2024-4-6 09:23
老师,我做NVT的时候材料分子只在真空层自己振动而不和晶体表面发生反应。请问既然NPT可以调节水层密度, ...

原地振动说明一开始摆得离表面太远,靠非键作用都难以自发吸附到表面上




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