VASP+BoltzTrap热电计算程序安装

1.准备的文件：

（1）.VASP5.2.11程序包；

（2）.BoltzTrap_vasp 程序包（来源 Wenqing Zhang课题组）；

 

2.VASP编译安装

（1）.修改VASP的main.F

在 CALL DUMP_ALLOCATE(IO%IU6)前添加以下代码，目的是产生SYSMETRY文件：

----------------------------------------------------------

!added by John. Yang for TE calculations

     OPEN(UNIT=38,FILE="SYMMETRY",STATUS="REPLACE")

     WRITE(38,"(I12)")NROTK

     DO NEDOS=1,NROTK   ! I DON'T WANT TO DEFINE A NEW INTEGER VARIABLE, SO HERE I USE NEDOS

        WRITE(38,"(3F10.5)") IGRPOP(1,1,NEDOS)*1.0,IGRPOP(1,2,NEDOS)*1.0,IGRPOP(1,3,NEDOS)*1.0

        WRITE(38,"(3F10.5)") IGRPOP(2,1,NEDOS)*1.0,IGRPOP(2,2,NEDOS)*1.0,IGRPOP(2,3,NEDOS)*1.0

        WRITE(38,"(3F10.5)") IGRPOP(3,1,NEDOS)*1.0,IGRPOP(3,2,NEDOS)*1.0,IGRPOP(3,3,NEDOS)*1.0

        WRITE(38,*)""

     ENDDO

     CLOSE(38)

-----------------------------------------------------------------------------

·        对vasp5.2.11版本，以上修改经过测试，可编译成功。

·        对vasp4.6版本，在"CALL TIMING"前添加，测试可编译成功。

 

（2）编译VASP,具体设置参照之前的文章：vasp5.2安装(fftw3.3+Gotoblas)

 

3.BoltzTrap_vasp程序包编译：

（1）参照组里的编译环境，对Makefile进行如下修改：主要改动BLAS,LAPACK,LIBS这几个参数，以下仅显示makefile的部分内容；

------------------------------------------------------------

SHELL = /bin/sh

FC =ifort

FOPT = -FR -mp -C -g

FOPT = -FR -mp

LINKER  =       $(FC)

LFLAGS  =      

FGEN =

BLAS= /usr/local/lib/libgoto2.so

LAPACK= ../vasp5.2.11/vasp.5.lib/lapack_double.o

LIBS  = -L../vasp5.2.11/vasp.5.lib -ldmy \

    ../vasp5.2.11/vasp.5.lib/linpack_double.o $(LAPACK) \

    $(BLAS)

LDFLAGS = -L/export/mathlib/mkl/8.1.1/lib/em64t -Vaxlib -pthread

DESTDIR = .

EXECNAME = BoltzTrap_vasp

------------------------------------------------------------

（2）cp BoltzTrap_vasp ~/bin/

 

4,计算，

（1）VASP自洽计算（待考证是自洽计算还是DOS计算）

（2）建立文件：文件夹名.intrans

文件夹名.intrans(以下参数设置待考证，可以参考UserGuide.pdf)

--------------------------------------------------------------

VASP                                                                        

0 0 0 0.0                          

0.233 0.0001 0.4  52

CALC                      # CALC (calculate expansion coeff), NOCALC read from file            

5                         # lpfac, number of latt-points per k                                  

BOLTZ                     # run mode (only BOLTZ is supported)                                  

0.15                       # (efcut) energy range of chemical potential                        

900. 900.                 # Tmax, temperature grid                                            

0.             # energyrange of bands given individual DOS output sig_xxx and dos_xxx (xxx is band number)

------------------------------------

（3）执行：/home/zjma/setup/BoltzTrap_vasp/x_trans BoltzTrap_vasp

（4）结果在 文件夹名.trace及.condtens 文件中。

转自：http://blog.sciencenet.cn/home.php?mod=space&uid=90975&do=blog&id=772892