xiaoqiugood的个人博客分享 http://blog.sciencenet.cn/u/xiaoqiugood

博文

报错信息收集:wannier函数处理能带结构等

已有 18202 次阅读 2014-6-17 09:12 |系统分类:科研笔记

关注:

1) 从报错当中寻找wannier函数成功处理能带结构的信息

 

1. 如果wannier90.win中注释掉exclude,则num_bands的设置(num_bands=30)将被程序忽略,直接采用NBANDS的指定数量(NBANDS=48)

 

 

 

 

 

num_wann=24

num_bands=30


# for GW uncomment
#exclude_bands 31-48

Begin Projections
Sc:l=0;l=1;l=2
H: l=0
End Projections


dis_froz_max=9
dis_num_iter=1000
guiding_centres=true


begin unit_cell_cart

 

运行wannier90.x wannier90

产生.chk与.wout文件,下面的三个文件由早前的VASP计算产生

du -sh wannier90.*
14M     wannier90.amn
540K    wannier90.eig
160M    wannier90.mmn

 

1.  qsub提交时等号后面不能留有空

#PBS -N vasp
#PBS -l nodes=1:ppn=16
#PBS -l walltime=7200:00:00

 

 

2.  cat  wannier.win

 

num_wann =    8

num_bands=    8

# Bandstructure plot
restart         =  plot
bands_plot      =  true


begin kpoint_path
G 0.000  0.000  0.000   A 0.000  0.000  0.500
A 0.000  0.000  0.500   H -0.333  0.667  0.500
H -0.333  0.667  0.500  K -0.333  0.667  0.000
K -0.333  0.667  0.000  G 0.000  0.000  0.000
G 0.000  0.000  0.000   M 0.000  0.500  0.000
M 0.000  0.500  0.000   L 0.000  0.500  0.500
L 0.000  0.500  0.500   H -0.333  0.667  0.500
end kpoint_path


bands_num_points 40
bands_plot_format gnuplot xmgrace

use_bloch_phases = T

 

 

wannier90.wout

 

 

Ignoring <num_bands> in input file
Ignoring <mp_grid> in input file
Exiting.......
 Error: Cannot use bloch phases for disentanglement [解开纠结]

3. 采用如下的设置成功得到band.dat

 

 

# #######################step 4.2:  Analysis of hse band, by using  LWANNIER90_RUN = .TRUE.  #################################

cat > INCAR << EOF

System  = SrVO3-hse-band

ISMEAR = -5
EMIN = -20 ; EMAX = 20 ; NEDOS = 1000  # usefull energy range for density of states

EDIFF = 1E-8                           # high precision for groundstate calculation

#KPAR = 3
#LHFCALC = .TRUE.  ; HFSCREEN = 0.2  ;
NBANDS = 64
#PRECFOCK = Fast   ; NELM = 1
#ALGO = Eigenval

LWAVE = .FALSE.                        # do not write the wave functions

#LORBIT = 11                           # for DOS

LWANNIER90_RUN = .TRUE.     # for band calculation

EOF

cat > wannier90.win << EOF

#num_wann =    8

#num_bands=    8   [由程序自动指出该参数]

#use_bloch_phases = T

# Bandstructure plot
restart         =  plot
bands_plot      =  true


begin kpoint_path
G 0.000  0.000  0.000   A 0.000  0.000  0.500
A 0.000  0.000  0.500   H -0.333  0.667  0.500
H -0.333  0.667  0.500  K -0.333  0.667  0.000
K -0.333  0.667  0.000  G 0.000  0.000  0.000
G 0.000  0.000  0.000   M 0.000  0.500  0.000
M 0.000  0.500  0.000   L 0.000  0.500  0.500
L 0.000  0.500  0.500   H -0.333  0.667  0.500
end kpoint_path


bands_num_points 40
bands_plot_format gnuplot xmgrace
EOF
# start calculation

$PARA -n $NP -machinefile .nodelists.$$ $EXEC  > STDOUT
 
cp OUTCAR OUTCAR-4-2
cp STDOUT STDOUT-4-2

 

 

运行完毕后,wannier90.win具有如下形式:

#num_wann =    8

#num_bands=    8

#use_bloch_phases = T

# Bandstructure plot
restart         =  plot
bands_plot      =  true


begin kpoint_path
G 0.000  0.000  0.000   A 0.000  0.000  0.500
A 0.000  0.000  0.500   H -0.333  0.667  0.500
H -0.333  0.667  0.500  K -0.333  0.667  0.000
K -0.333  0.667  0.000  G 0.000  0.000  0.000
G 0.000  0.000  0.000   M 0.000  0.500  0.000
M 0.000  0.500  0.000   L 0.000  0.500  0.500
L 0.000  0.500  0.500   H -0.333  0.667  0.500
end kpoint_path


bands_num_points 40
bands_plot_format gnuplot xmgrace
num_wann =    64  ! set to NBANDS by VASP

use_bloch_phases = T

begin unit_cell_cart

.....(下略)

 



https://blog.sciencenet.cn/blog-567091-804012.html

上一篇:GW计算及输出文件:VASP手册说明文档解读
下一篇:自旋轨道耦合计算:VASP手册实例分析
收藏 IP: 61.157.130.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-15 16:22

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部