计算化学公社

标题: 求教：vasp计算分子动力学的困惑 [打印本页]

作者
Author: 飞翔的猪 时间: 2019-10-25 20:38
标题: 求教：vasp计算分子动力学的困惑
大家好，最近刚接触使用vasp计算分子动力学有如下困惑不是很理解：
1）根据官网手册 ISIF=3对应的是NPT系综恒温恒压过程，但是我看OSZICAR中信息温度并不恒定，波动很大并不恒温https://cms.mpi.univie.ac.at/wiki/index.php/Molecular_dynamics_calculations 但是查看到小*虫上的帖子：说vasp计算主要使用的是NVT和NVE系综 http://mu****g.com/t-1744011-1 有的说vasp计算动力学先NVE再NVT 。自己对于这些信息不是很能够判断信息的正误，想跟大家交流一下

2）对于同一个结构优化文件，MS中的Forcite模块（NPT）计算出来的能量是非常稳定的而vasp计算出来的能量是一直增大的，是因为vasp中步数计算的不够而导致的吗？vasp计算分子动力学一般需要计算多少步

3）用以下输入文件计算出现以下警示可能是什么原因造成的The initial velocities result in a center of mass drift but          |
|    there must be no drift.                                              |
|    The drift will be removed !

以下是我的INCAR
ENCUT=150
NSW=300
IBRION=0
ISIF=3
IALGO=48
LREAL=Auto
NELMIN=8
BMIX=2.0
MAXMIN=50
ISYM=0
POTIM=3
TEBEG=298
LWAVE=.FALSE
LCHARG=.FALSE
ISMEAR=0
SIGMA=0.1

KPOINTS：G 1 1 1

POSCAR：使用的是结构优化后的文件（二维材料，没有固定原子坐标）

作者
Author: bobosiji 时间: 2019-10-26 17:43
1.没看懂问题
2.MS中的Forcite模块（NPT）计算出来的能量是非常稳定的而vasp计算
Forcite和vasp根本就不是一回事，没有可比性

你的vasp是自学成才还是课题组培训呀？

作者
Author: agent99 时间: 2019-10-27 07:31
ENCUT=150 是不是太小了点

作者
Author: 飞翔的猪 时间: 2019-10-27 14:33

bobosiji 发表于 2019-10-26 17:43
1.没看懂问题
2.MS中的Forcite模块（NPT）计算出来的能量是非常稳定的而vasp计算
Forcite和vasp根本就不 ...

一般自己先尝试解决不清楚的再去跟老师沟通

作者
Author: 飞翔的猪 时间: 2019-10-27 14:35

agent99 发表于 2019-10-27 07:31
ENCUT=150 是不是太小了点

ENCUT设置主要参考：https://cndaqiang.github.io/2018/12/07/MD/

作者
Author: sige 时间: 2019-10-28 10:49
VASP是AIMD，是第一性原理动力学，forcite是分子动力学，二者没有什么可比性。VASP跑动力学控温很差，能量震荡也正常，时间尺度至少要到ps级别

作者
Author: 飞翔的猪 时间: 2019-10-28 16:49

sige 发表于 2019-10-28 10:49
VASP是AIMD，是第一性原理动力学，forcite是分子动力学，二者没有什么可比性。VASP跑动力学控温很差，能量 ...

非常感谢回复这两个的区别之前是没有弄清楚谢谢

作者
Author: 李朋 时间: 2019-11-12 19:49

飞翔的猪发表于 2019-10-27 14:35
ENCUT设置主要参考：https://cndaqiang.github.io/2018/12/07/MD/

150太小了，至少300以上

作者
Author: dindindindindin 时间: 2020-4-27 06:40
学习了

作者
Author: haitaier 时间: 2024-6-3 13:17

飞翔的猪发表于 2019-10-27 14:35
ENCUT设置主要参考：https://cndaqiang.github.io/2018/12/07/MD/

ENCUT的设置要根据自己的体系来修改，不能看别人是啥就是啥啊，看看类似体系别人算的时候选的大概多少

欢迎光临计算化学公社 (http://ccc.keinsci.com/)