使用CP2K做NPT模拟时遇到slab严重变形问题

Stardust0831 · 发表于 Post on 2024-4-5 22:10:41

本帖最后由 Stardust0831 于 2024-4-6 00:00 编辑

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

sobereva · 发表于 Post on 2024-4-6 00:17:02

你这个体系也没有水，本来就不应该用NPT。一直用NVT就完了

Stardust0831 · 发表于 Post on 2024-4-6 09:23:35

sobereva 发表于 2024-4-6 00:17
你这个体系也没有水，本来就不应该用NPT。一直用NVT就完了

老师，我做NVT的时候材料分子只在真空层自己振动而不和晶体表面发生反应。请问既然NPT可以调节水层密度，能不能也调节材料分子层密度、帮助它接近晶体表面？如果继续做NVT，会面临材料分子不靠近晶体表面的问题，是否建议加大分子数量来解决？有其他更好的解决方法吗？

sobereva · 发表于 Post on 2024-4-7 05:22:08

Stardust0831 发表于 2024-4-6 09:23
老师，我做NVT的时候材料分子只在真空层自己振动而不和晶体表面发生反应。请问既然NPT可以调节水层密度， ...

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 使用CP2K做NPT模拟时遇到slab严重变形问题