处理计算本征值.pdfVIP

处理VASP计算得到的本征值文件

在VASP计算能带结构时，输出的本征值文件为EIGENVAL，它的格式如下：例子

2211

0.3484800E+030.1200000E-080.1200000E-080.4840000E-090.5000000E-15

1.00000000000000D-004

CAR

Al

62615

0.0000000E+000.0000000E+000.0000000E+000.3846154E-01

1-10.1916

2-5.0683

3-4.2171

4-4.2171

5-3.1342

6-0.4730

70.4006

80.6271

91.0307

101.0307

111.1635

121.3672

131.3672

142.0885

152.1586

。。。。。。。。。。。。。

第一行，前三个整数无意义，第四个整数，如果是2,表示是自旋极化的计算，如果是1,表示非

自旋极化的计算。

第2至5行的数据含义不大明确，可以不管它。

第6行的数据表示：第一个数表示体系总的价电子数目，第二个数表示的计算能带时总的k点数

目，第三个数表示的是计算能带时计算了多少条能带。

第8行的前三个数是k点的坐标，第四个数是相应k点的权重。

第9行给出的是该k点对应的本征值的序号(即第几条能带)，及相应的本征值。

下面是根据pwscf提供的band_plot.f90文件略改写后针对VASP的EIGENVAL文件，把它转换

为origin或xmgrace软件能直接画图的数据。

下面是band_plotvasp.f90源代码：

!forVASP

!

programprog

real,allocatable::e(:,:),eup(:,:),edn(:,:)

real,allocatable::k(:,:)

real,dimension(3)::k0,a

real,dimension(6)::xxxx

character(len=32)::xx,yy

write(6,*)fermilevel(eV)

read(5,*)ef

open(10,file=EIGENVAL,status=old)

open(11,file=bnd.dat)

read(10,*)iii,iii,iii,ispin

read(10,*)(xxxx(i),i=1,5)

read(10,*)xxxx(6)

read(10,*)xx

read(10,*)yy

read(10,*)nn,nk,nbands

allocate(e(nk,nbands))

allocate(k(nk,3))

if(ispin.eq.2)then

doi=1,nk

read(10,*)

read(10,*)(k(i,j),j=1,3),wtk

doj=1,nbands

read(10,*)jj,eup(i,j),edn(i,j)

enddo

write(13,9030)(eup(i,j),j=1,nbands)

write(14,9030)(edn(i,j),j=1,nbands)

enddo

write(6,*)(k(i,j),j=1,3)

read(10,*)(e(i,n),n=1,nbands)

write(13,9030)(e(i,n),n=1,nbands)

else

doi=1,nk

read(10,*)

read(10,*)(k(i,j),j=1,3),wtk

doj=1,nbands

read(