Dear All,
I calculate band structure and hopping matrix of monolayer Graphene by using wannier90-2.0.1. There are two atoms in one unit cell, the following is input information :
INCAR :
IBRION = -1
ECUT = 500
EDIFF = 1e-5
IALGO = 38
PREC = Accurate
ISMEAR = -5
SIGMA = 0.05
LREAL = F
#LWAVE = T
NBANDS = 16
##
ALGO = None
NELM = 0
LCHARG = F
LWAVE = F
LWANNIER90=.TRUE.
wannier90.win :
num_wann=8
num_bands=16
# for GW uncomment
#exclude_bands 16-64
Begin Projections
C:s,px,py,pz
End Projections
# parameters for disentanglement
dis_froz_max=5
dis_froz_min=-5
dis_num_iter=1500
dis_conv_tol = 1.0d-10
dis_conv_window = 20
# parameters for wannierisation
num_iter = 2000
conv_tol = 1.0d-10
conv_window = 20
guiding_centres=true
# output hamiltonian matrix
#write_hr = true
hr_plot = true
# Bandstructure plot
#restart = plot
bands_plot = true
begin kpoint_path
G 0.0000000 0.0000000 0.0000000 K 0.3333333 0.6666667 0.0000000
K 0.3333333 0.6666667 0.0000000 M 0.5000000 0.5000000 0.5000000
M 0.5000000 0.5000000 0.5000000 G 0.0000000 0.0000000 0.0000000
end kpoint_path
bands_num_points 40
bands_plot_format gnuplot xmgrace
begin unit_cell_cart
2.4667090 0.0000000 0.0000000
1.2333545 2.1362327 0.0000000
0.0000000 0.0000000 15.1023000
end unit_cell_cart
begin atoms_cart
C 0.6166773 0.3560388 0.0000000
C 3.0833862 1.7801939 0.0000000
end atoms_cart
mp_grid = 6 6 1
begin kpoints
0.000000000000 0.000000000000 0.000000000000
0.166666666667 0.000000000000 0.000000000000
0.333333333333 0.000000000000 0.000000000000
0.500000000000 0.000000000000 0.000000000000
0.333333333333 0.166666666667 0.000000000000
0.500000000000 0.166666666667 0.000000000000
-0.333333333333 0.333333333333 0.000000000000
-0.166666666667 -0.000000000000 0.000000000000
-0.000000000000 -0.166666666667 0.000000000000
0.000000000000 0.166666666667 0.000000000000
-0.166666666667 -0.166666666667 0.000000000000
0.166666666667 0.166666666667 0.000000000000
-0.333333333333 0.000000000000 0.000000000000
0.000000000000 -0.333333333333 0.000000000000
0.000000000000 0.333333333333 0.000000000000
-0.333333333333 -0.333333333333 0.000000000000
0.333333333333 0.333333333333 0.000000000000
-0.000000000000 -0.500000000000 0.000000000000
-0.500000000000 -0.500000000000 0.000000000000
-0.333333333333 -0.166666666667 0.000000000000
-0.166666666667 -0.333333333333 0.000000000000
0.166666666667 0.333333333333 0.000000000000
0.166666666667 -0.166666666667 0.000000000000
-0.166666666667 0.166666666667 0.000000000000
-0.500000000000 -0.166666666667 0.000000000000
-0.166666666667 -0.500000000000 0.000000000000
0.166666666667 0.500000000000 0.000000000000
0.166666666667 -0.333333333333 0.000000000000
-0.166666666667 0.333333333333 0.000000000000
-0.333333333333 -0.500000000000 0.000000000000
0.333333333333 0.500000000000 0.000000000000
0.333333333333 -0.166666666667 0.000000000000
-0.333333333333 0.166666666667 0.000000000000
0.500000000000 0.333333333333 0.000000000000
-0.500000000000 -0.333333333333 0.000000000000
0.333333333333 -0.333333333333 0.000000000000
end kpoints
part of output wannier_hr.dat :
0 0 0 1 1 #5.975871 -0.000000
0 0 0 2 1 -1.089392 0.557442
0 0 0 3 1 -0.139528 -0.085332
0 0 0 4 1 0.155557 -0.074815
0 0 0 5 1 0.048334 0.012936
0 0 0 6 1 -0.040701 0.012192
0 0 0 7 1 -0.000254 -0.005978
0 0 0 8 1 0.140810 -0.003572
0 0 0 1 2 -1.089392 -0.557442
0 0 0 2 2 #-3.873225 -0.000000
0 0 0 3 2 3.091729 1.850427
0 0 0 4 2 2.671730 0.889670
0 0 0 5 2 -0.057204 -0.033336
0 0 0 6 2 0.183550 -0.012852
0 0 0 7 2 0.063698 -0.025852
0 0 0 8 2 0.199665 0.137492
0 0 0 1 3 -0.139528 0.085332
0 0 0 2 3 3.091729 -1.850427
0 0 0 3 3 # -9.031542 0.000000
0 0 0 4 3 -2.032979 0.353178
0 0 0 5 3 -0.003251 0.006119
0 0 0 6 3 0.238048 -0.035893
0 0 0 7 3 -0.085293 0.054229
0 0 0 8 3 0.377565 0.008021
0 0 0 1 4 0.155557 0.074815
0 0 0 2 4 2.671730 -0.889670
0 0 0 3 4 -2.032979 -0.353178
0 0 0 4 4 #-10.109918 -0.000000
0 0 0 5 4 -0.036041 0.016305
0 0 0 6 4 -0.777184 0.094481
0 0 0 7 4 0.054059 -0.024728
0 0 0 8 4 -1.182728 -0.226013
0 0 0 1 5 0.048334 -0.012936
0 0 0 2 5 -0.057204 0.033336
0 0 0 3 5 -0.003251 -0.006119
0 0 0 4 5 -0.036041 -0.016305
0 0 0 5 5 #5.601534 -0.000000
0 0 0 6 5 0.229079 -0.103372
0 0 0 7 5 -0.141783 0.103966
0 0 0 8 5 0.697744 -0.056069
0 0 0 1 6 -0.040701 -0.012192
0 0 0 2 6 0.183550 0.012852
0 0 0 3 6 0.238048 0.035893
0 0 0 4 6 -0.777184 -0.094481
0 0 0 5 6 0.229079 0.103372
0 0 0 6 6 -#0.667422 -0.000000
0 0 0 7 6 -0.247113 0.097277
0 0 0 8 6 3.423367 1.060382
0 0 0 1 7 -0.000254 0.005978
0 0 0 2 7 0.063698 0.025852
0 0 0 3 7 -0.085293 -0.054229
0 0 0 4 7 0.054059 0.024728
0 0 0 5 7 -0.141783 -0.103966
0 0 0 6 7 -0.247113 -0.097277
0 0 0 7 7 #-0.698032 -0.000000
0 0 0 8 7 0.083201 0.011267
0 0 0 1 8 0.140810 0.003572
0 0 0 2 8 0.199665 -0.137492
0 0 0 3 8 0.377565 -0.008021
0 0 0 4 8 -1.182728 0.226013
0 0 0 5 8 0.697744 0.056069
0 0 0 6 8 3.423367 -1.060382
0 0 0 7 8 0.083201 -0.011267
0 0 0 8 8 #-9.878974 0.000000
Some numbers should be equal : h11=h55, h22=h66, h33=h77, h44=h88, but why my results not show? any wrong input parameters? (Mark # in the *_hr.dat)