求教：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：使用的是结构优化后的文件（二维材料，没有固定原子坐标）

bobosiji

1.没看懂问题
2.MS中的Forcite模块（NPT）计算出来的能量是非常稳定的 而vasp计算
Forcite和vasp根本就不是一回事，没有可比性

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

agent99

ENCUT=150 是不是太小了点

飞翔的猪

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

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

飞翔的猪

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

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

sige

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

飞翔的猪

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

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

李朋

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

150太小了，至少300以上

dindindindindin

学习了

haitaier

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

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