📆分子动力学计算模版：
ISTART = 0
ICHARG = 2
ISPIN = 1
LCHARG = .F.
LWAVE = .F.
PREC = N
#ALGO = N
NELM = 100
#EDIFF = 1E-6
ENCUT = 500
LREAL = A
#IBRION = 2
ISIF = 2
#NSW = 500
#EDIFFG = -0.05
IVDW = 11 #DFT-D3
ISMEAR = 0
SIGMA = 0.05
#Molecular Dynamics
IBRION = 0
NSW = 2000
EDIFF = 1E-5
POTIM = 1
SMASS = 0 # for NVT&NPT
TEBEG = 300
TEEND = 300
MDALGO = 2 #2for NVT 3for NPT
ALGO = Veryfast
ISYM = 0
NBLOCK = 1
IWAVPR = 11 #外扩参数详解：
分子动力学模拟开启的参数是 IBRION ，当 IBRION = 0 的时候表示开启分子动力学。
由于分子动力学计算比较消耗资源，因此这里需要设置 ALGO = Veryfast 用于加速收敛。同时可以适量减小 EDIFF 。
当 IBRION = 0 时，EDIFFG 无效。
NSW
分子动力学步数，单位：fs。
POTIM
每一步的步长，一般设置为 0.5，1，2。
TEBEG TEEND
分子动力学模拟的初始温度和终止温度，当两个参数值一样的时候，表示恒温模拟，适用于 NPT，NVT 系综。
ISYM = 0
建议跑分子动力学开启。
SMASS MDALGO NBLOCK
当 SMASS = -1 时，要把 MDALGO 注释掉，表示开启 NVT 系综的退火升温模拟，采用速度缩放算法快速升温，此时 NBLOCK 的值代表多少步长强制控制一下温度。当 SMASS >= 0 时，此时设置 MDALGO 可以设置不同的热浴器。MDALGO = 2 代表开启 Nosé-Hoover 热浴器。此时NBLOCK 的值代表多少步保存一下轨迹文件。
IWAVPR
外扩，用于加速收敛。