计算化学公社

标题: VASP计算AIMD,结构变化程度大 [打印本页]

作者
Author:
种花家的蓝兔    时间: 2022-8-23 16:40
标题: VASP计算AIMD,结构变化程度大
本帖最后由 种花家的蓝兔 于 2022-8-23 16:41 编辑

请问各位前辈,用VASP计算分子动力学,目的想看结构的稳定性。一个单层二维结构,参考文献上写明了独立的这个单层二维结构是稳定的,但是我计算将近1000步之后(1ps),结构严重变形,如下图。【用的和文献一样的系综、恒温器、时间步长】按理说有文献做出来它是稳定的,应该是可以证明稳定性的吧。所以我想请问各位前辈老师,做AIMD时,对分子结构影响程度大的参数有哪些

作者
Author:
Eudaimonia    时间: 2022-8-23 18:51
AIMD之前弛豫过么,弛豫一下会好点
作者
Author:
种花家的蓝兔    时间: 2022-8-23 19:04
Eudaimonia 发表于 2022-8-23 18:51
AIMD之前弛豫过么,弛豫一下会好点

弛豫过了
作者
Author:
waitingseven    时间: 2022-8-23 20:09
真空层厚度也跟文献一致吗
作者
Author:
种花家的蓝兔    时间: 2022-8-23 20:26
waitingseven 发表于 2022-8-23 20:09
真空层厚度也跟文献一致吗

是的,at least 20埃
作者
Author:
Eudaimonia    时间: 2022-8-23 21:47
按我不多的AIMD经验来看,体系大一点更稳定,我单胞跑5ps很容易结构跑乱,但是2x2超胞就不会
作者
Author:
waitingseven    时间: 2022-8-23 22:03
种花家的蓝兔 发表于 2022-8-23 20:26
是的,at least 20埃

是不是用的K点不一样,我看晶胞很小,如果直接用gamma点肯定不行
作者
Author:
种花家的蓝兔    时间: 2022-8-23 22:35
waitingseven 发表于 2022-8-23 22:03
是不是用的K点不一样,我看晶胞很小,如果直接用gamma点肯定不行

不是gamma点,MP那个
作者
Author:
种花家的蓝兔    时间: 2022-8-23 23:10
Eudaimonia 发表于 2022-8-23 21:47
按我不多的AIMD经验来看,体系大一点更稳定,我单胞跑5ps很容易结构跑乱,但是2x2超胞就不会

我这是3×3的。请问您的结构是单层的吗
作者
Author:
Eudaimonia    时间: 2022-8-24 08:34
种花家的蓝兔 发表于 2022-8-23 23:10
我这是3×3的。请问您的结构是单层的吗

单层石墨,单原子的结构,按原始的石墨晶胞来算我那个可能有14x14了
作者
Author:
种花家的蓝兔    时间: 2022-8-24 08:43
Eudaimonia 发表于 2022-8-24 08:34
单层石墨,单原子的结构,按原始的石墨晶胞来算我那个可能有14x14了

是用的NVT系综吗,如果是的话,方便说一下SMASS设置的值是多少吗
作者
Author:
Eudaimonia    时间: 2022-8-24 09:01
种花家的蓝兔 发表于 2022-8-24 08:43
是用的NVT系综吗,如果是的话,方便说一下SMASS设置的值是多少吗

不太确定,我用的cp2k计算,分子动力学相关的设置只设置了Nose-Hoover热浴,我也不确定算哪种
作者
Author:
种花家的蓝兔    时间: 2022-8-24 09:06
Eudaimonia 发表于 2022-8-24 09:01
不太确定,我用的cp2k计算,分子动力学相关的设置只设置了Nose-Hoover热浴,我也不确定算哪种

好的,谢谢您
作者
Author:
CHLllooll    时间: 2023-12-11 22:08
种花家的蓝兔 发表于 2022-8-24 09:06
好的,谢谢您

前辈你好,请问问题解决了吗




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