这里写自定义目录标题
算了11层黑磷结构的gamma点的价带,用getwave程序导出了可以被xcrysden、vesta等打开的xsf类型的文件,xsf跟bxsf一样,其实就是一个三维密度图,有grid。在这个过程中出现了两个错误
1. getwave的时候报错
报错信息:
number of generated plane wave coeffecients isn’t equal that from WAVECAR!
有一个方法就是在运行getwave的时候加入-g参数,把OUTCAR里面的NGX、NGY、NGZ作为输入网格输入进来。
getwave -k 1 -b 55 -g 24 24 2500 > G_VBM.xsf
由于z轴就1个K点,所以网格比较密,这一步花的时间就比较长。然后生成了个xsf后缀的文件。
2. 用vesta打开xsf的时候报错
把已经生成的xsf文件导入到vesta里面,却发现闪退了,一开始以为是文件太大了,有100多M的原因,但是后来尝试着打开了一个100多M的CHGCAR文件,VESTA是能打开的,所以就肯定是xsf文件本身的问题。
用写字板打开xsf文件,发现是由于数值太大,而预设的两个变量值之间的空间太小,导致两个数碰到一起,变成了一个数。
如图,黄色部分是我改好的了,原本里面是红色部分的样子。
在xsf后缀的文件中,网格太密可能会导致出现这种现象的地方一般在晶格矢量部分、原子坐标部分、BEGIN_DATAGRID_3D后面的三个网格数部分,要特别注意检查。
我手动的在这些地方加了一个空格,然后VESTA就能读取啦。完美解决。