|
|
[求助]
求助VASP进行MD计算时INCAR设置-为何原子跑到了一起?
大家好,我想作H2在Ni(111)表面的吸附-分解的MD模拟,用VASP。我先优化了下H2在表面的吸附,H2距离表面约4A,然后开始跑动力学,VASP大概跑了五步,就停了,这时候力变的非常到一下子到6百多了。我看了下结构,发现两个H原子竟然重叠在一起了(请见图)……下面是我的INCAR,请大家帮忙看下是否是INCAR设置不对?NVT系综。非常感谢。
Ni(111),2层,底下一层固定,上面一层T T T,真空层10A
INCAR内容:
#PREC=Accurate
ISTART=0 ; ICHARG=2
ENCUT=400
ISMEAR=1 ; SIGMA=0.1
ALGO=V
#NELMDL=5
NSW=1000
POTIM=3
IBRION=0
ISIF=2
ISPIN=2
ISYM=0
NBLOCK=5
KBLOCK=10
#APACO=10
#NPACO=200
LCHARG=F
TEBEG=300
TEEND=300
SMASS=2
#EDIFF=1E-4;
#EDIFFG=-0.05
LPLANE=.TRUE.
NPAR=4
NSIM=4
LREAL=Auto
![]()
初始结构.jpg
![]()
中断结构-HH原子重叠.jpg |
» 猜你喜欢
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有5人回复
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有10人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有11人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
|