计算表面态的详细过程.docVIP

外表态计算的详细过程

其实，即算外表态的方法，本质上是计算任意指定原子在能带中占有的比重。

第一步，计算。

按照一般的流程，作结构优化、静态计算、能带计算。

如果，你的结构需要考虑范德瓦耳斯作用力，那么在结构优化的计算中，选用适宜的vdw计算方法；

如果，你的计算需要考虑自旋轨道耦合，那么在静态和能带计算的时候，参加soc的参数。

在算能带的时候，根据你的需要，设定K点的轨道。

因为算的是指定原子在能带中所占的比重，建议算能带的K点一定要多一些，这样点才可以在最终画出来的图中练成线，看出图案，比方101，甚至更大，这个可以调试。

第二步，处理数据。

首先，你要有一个可以算外表态的脚本，或者由这个脚本生成的小程序。

以下是是Fortran程序编写的脚本原文，可以用来处理由VASP计算的数据：

!*********************************************************************************

!*得出每种元素对每条能带的每个K点的奉献值(大于给定的cutoff)*

!所需文件POSCAR,POTCAR,PROCAR,bnd.dat

!输出文件PRO,H//bnd_up.dat,H//bnd_down.dat

!是否考虑自旋

!根据PRO文件确定Cutoff

!*********************************************************************************

PROGRAMForPROCAR

CHARACTER(len=80)HEADER,filename

CHARACTER(len=2)H

character(len=132)line

integerspec_counter,lower,upper

parameter(n=10000)

character*4aindex(n)

character*50KP(:,:),PK(:,:)

realKK(:,:,:),BB(:,:,:),Natoms(:),B(:,:,:),Cut(:)

allocatableKK,BB,Natoms,B,Cut,KP,PK,NatomJia(:)

INTEGERA,I,J,K,Nk,Nb,Na,EleN,Nok,soc

C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

OPEN(20,FILE=PROCAR)

open(102,file=PRO)

READ(20,*)

read(20,30)Nk,Nb,Na

write(*,*)Nk,Nb,Na

30format(15x,I5,20X,I5,17x,I5)

allocate(KK(Nk,Nb,Na))

!*********************************************************************************

!*ChecktheatomlabelsandspeciesinPOTCAR*

!*********************************************************************************

open(15,file=POTCAR,status=OLD)

EleN=0

dowhile(.not.EOF(15))

read(15,(a),end=10)line

if(index(line,VRHFIN).ne.0)then

EleN=EleN+1

upper=index(line,:)-1

lower=index(line,=)+1

aindex(EleN)=line(lower:upper)

print*,aindex(EleN)

endif

enddo

close(15)

10continue

!*********************************************************************************

!*Checktheeverys