本次计算借助二维材料2dmatpedia数据库给出输入文件,并结合自己的工具进行能带计算。
第一步,下载所需要计算体系的输入文件
理论上需要首先进行结构优化,比如这里的relax1,relax2。
其中结构优化的作用是重新得到POSCAR,结构优化形成的文件写在CONTCAR中
执行命令
cp CONTCAR POSCAR
然后开始静态计算,这里要特别注意,这里二维材料的静态计算中static考虑所谓一个二维的范德瓦斯力,这样需要额外的文件,否则计算量特别大,所以这里使用novdwstatic计算,由于直接借助了别人的结论,所以我们可以直接进行静态自洽计算,不需要进行结构优化。
最后,进行能带计算,能带计算使用vaspkit生成对称的KPOINTS文件,数据库中没有提供此文件,拷贝自洽计算的WAVECAR和CHGCAR文件到能带计算中进行计算。
计算中暂时使用PBE赝势能。
最后拷贝Vasprun.xml 和KPOINS文件到Python目录下,使用pymatgen 工具进行能带分析。