vasp 非自洽计算 包括能带计算吗
怎么画面心立方晶格的倒格子呢?是由原胞基矢画呢,还是怎么画。求详细解答 面心立方的倒格子是体心立方的形式。倒格子是相对于原胞基失的,一般情况下是不能直观地画出来的,而应该计算。先写出原胞基失三个向量,然后按公式(固体书上必然有)求出。
能带结构的程序 计算材料的能带结构即色散曲线E(k),步骤为(并以计算fcc结构Al的能带结构为例进行说明):根据特殊k点的走向,选取特殊k点及特殊k点间的分割点数,准备好产生k点的输入文件syml6。特殊k点的个数20 20 20 10 20。特殊k点间的分割点数X 0.5 0.0 0.5。特殊k点的坐标,相对于倒格子矢量G 0.0 0.0 0.0L 0.5 0.5 0.5W 0.5 0.25 0.75K 0.375 0.375 0.75G 0.0 0.0 0.0。下面三行,前三列是正格子基矢,后三列是倒格子基矢0.000000000 1.987500000 1.987500000-0.251572327 0.251572327 0.2515723271.987500000 0.000000000 1.987500000 0.251572327-0.251572327 0.2515723271.987500000 1.987500000 0.000000000 0.251572327 0.251572327-0.25157232720.0 15.0。在画能带结构时,每个特殊k点所对应的竖线的能量范围7.068339。费米能级用程序gk.x产生k点,得到KPOINTS文件。注释:程序gk.x是由gk.f文件编译后得到的目标文件,其输入文件为syml,输出文件为KPOINTS,inp.kpt。紧接着利用前面计算得到的自洽电荷密度作一次非自洽的计算。采用命令解压保存的电荷密度文件chg.tgz:tar xzvf chg.tgz另外设置ISTART=1,ICHARG=11,并增加NBANDS的值,ISMEAR采用默认值SYSTEM=。