VASP
首先，做分子动力学的模型要原子足够多，可以建立超胞来处理比如小编可以在MS中做一个超胞，然后导出坐标,步骤如下：
打开MS，建立初始模型，这里小编以TiO2为例
点击 build-symmetry-supercell，因为小编的晶胞有12个原子，所以做了个3x3x3的超胞，也就是9x12=108个原子
然后把里面的坐标导出：Moduless-CASTEP-Calculation-files
注意这一步不要选择高对称性：如下图，选择no
然后在工作目录里找到.cell文件，里面就是坐标了，可以复制出来作为你的输入文件POSCAR来使用
接下来做INCAR文件：设置如下：ISYM = 0 PREC = Normal ENCUT = 400 ALGO = Fast EDIFF = 1E-4 EDIFFG = -0.02 LREAL = Auto ISIF = 0 NELM = 200 NELMIN = 5 NSW = 5000 IBRION = 0 TEBEG = 300 TEEND = 300 SMASS = 0 POTIM = 1.0 ISMEAR = 0 SIGMA = 0.2 NWRITE = 2 LWAVE = .FALSE. LCHARG = .FALSE.
其中第一行是体系名字
设置好之后，建议KPOINTS选取1 1 1 的gamma设置，来加快计算速度
TEBEG = 300 TEEND = 300这个参数是设置初始温度和终止温度的，一般一样，如果这两个参数相差大的话，计算特别慢哦，比如你想从200到300来计算，这样的计算量特别大哦
希望对你们有帮助，谢谢！