1、结构优化需要四个文件 incar 、po衡痕贤伎scar、 potcar、 kpoints。首先通过MS来建晶体模型，从而获得poscar里面的焊剑杂锭参数。如图所示，打开MS，选择file,new,选择3D模型。然后选择build、crystal，创建晶体。选择空间群和设置晶格常数。再选择添加原子和设置原子的坐标。完成以后再菜单选择导出，导出格式为cif
2、将上一步导出的cif文件导入vesta软件，在菜单栏选择export date，导出格式为vasp，确认后会提示选择坐标，根据需要选择分数或者笛卡尔坐标，推荐大家选择分数坐标。
3、打开导出的vasp文件，里面的内容即为poscr的内容，建立一个poscar，将其复制进入即可。
4、poscar建立好以后，建立incar，incar的内容如图，一次优化和固定体积不变优化的不同就在于足毂忍珩incar中ISIF的选择不同，首先说一次性优化，设ISIF=3。incar设置好后设置Kpoints，内容如图，k点的大小可根据需要自己选择。potcar为雁势文件，可在vasp中找到相对应的原子进行复制即可。四个文本文件建好之后提交到vasp进行计算，就可以一次性对结构进行优化。优化完成以后，可将contcar中的文件导入vesta，看优化后的结构是否符合自己的要求。
5、接下来说通过BM方程拟合进行优化，除了前面的蚵播噍疸四个文件以外，还需要一个脚本文件，脚本文件的内容如图，通过设置一系列的值，保葆扒豳峤持体积不变，即在INCAR中设置ISIF=4,通过脚本计算一系列的体积和能量的数据，然后通过BM方程拟合找到最低能量所对应的体积，即为该晶体最稳定的结构。
6、得到E-V数据后，将数据放入original,用自定义函数拟合的方法进行拟合。首咸犴孜稍先建立bm方程，选择tools，fitting function 泠贾高框organizer，输入BM方程，并设置好参数，进行编译，编译成功后保存。
7、bm方程建立好以后，做出散点图，然后按图所示，选择BM方程进行拟合，拟合之前要设置参数初始值，拟合完成后就可得到最优结构的体积。
