【做计算 找华算】理论计算助攻顶刊，50000+成功案例，全职海归技术团队、正版商业软件版权！
考虑一个结构稳不稳定，或者通过第一性原理分子动力学获得合理初始构型，需要使用VASP进行第一性原理分子动力学计算。
首先，结构的晶胞应该足够大，包含几十个原子，这样温度才会波动比较小。
其次，KPOINTS采用比较少的撒点，因为晶胞比较大。采用1 1 1的采样就可以。
KPOINTS(Gamma-only)
0
Gamma
1 1 1
0 0 0
第三，在INCAR里设置参数。降低EDIFF的精度，一般10-4, 10-5 eV/atom就很高了。加入分子动力学的参数：
#md
IBRION=0
NSW=100
NWRITE=0
TEBEG=300
TEEND=1300
SMASS=3
NBLOCK=50
POTIM=3
第四，注意INCAR里电子结构计算的一些参数设置，有些系统计算出来OSZICAR里温度T一直为0，经过摸索才发现与ISMEAR, SIGMA等有关：
ISMEAR = 0 ; SIGMA = 0.05
ALGO=Fast
#POTIM = 0.10
#Wavefunction and charge
LWAVE = FALSE
LCHARG = FALSE
#Target Pressure
#PSTRESS = 0.01
LREAL=Auto
NCORE=4
LVDW=.TRUE.
ISYM=0
第五，从OSZICAR里提取温度和能量。
grep "T= " OSZICAR | awk '{print $1 " " $3 " " $5}' > T_E.txt
第六，一些摸索的经验。POTIM不宜取得太大，否则会导致非常大的温度T，使得程序停掉。尝试过POTIM=5，温度infinite。POTIM取1.5，2都是不错的。
本文转自胡涛科学网博客，转载目的在于知识分享,如有侵权问题，请联系我们删除。https://blog.sciencenet.cn/blog-907836-1237593.html
想了解更多VASP相关理论与实操？关注华算科技系列培训课程，让你的入门更轻松！
做计算，找华算 🏅 我们提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。 🎯我们的理论计算服务，累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果，计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏
👉 点击阅读原文加我，探索更多优惠💖
