<div dir="ltr"><div>Hello, <br> I am calculating the band structure of bulk Bi2Te3 using GW method. I have recompiled the vasp code with vasp2wannier included and the band structure was calculated using wannier interpolation. The topology of the band structure with GW is almost same as the one without GW, but the Fermi level is not correct. That is, the Fermi level is now ~0.25 eV below the original VBM site. Could anyone give me some suggestion?  Thank you very muck!<br>
 The following are the wannier90.win file and the INCAR file in the last step I use:    <br> INCAR: <br><br> ISPIN = 2 <br> LSORBIT = .TRUE. <br> ISMEAR = 0 <br> SIGMA = 0.05 <br><br> LORBIT = 11 <br> NPAR = 16 <br> EDIFF = 1E-8 <br>
 NBANDS  = 288 <br> EMIN = -20; EMAX = 20 <br> NEDOS=2000 <br><br> ALGO = GW0 ; LSPECTRAL = .TRUE. ; NOMEGA = 288 <br> PRECFOCK = Fast <br> NKRED = 3 <br> LWANNIER90=.TRUE. <br><br><br> wannier90.win: <br><br> num_wann=30 <br>
 num_bands=60 <br><br> num_iter  = 200 <br><br> # for GW uncomment <br> exclude_bands 1-10,71-288 <br><br> #kmesh_tol = 0.0000001 <br> #use_bloch_phases = T <br><br> Begin Projections <br> Bi:px;py;pz <br> Te:px;py;pz <br>
 End Projections <br><br> dis_win_max   =  10.0 <br> dis_froz_max  =  7 <br> dis_mix_ratio = 1.0 <br><br> #dis_froz_max=20 <br> dis_num_iter= 1000 <br><br> spinors = .true. <br><br> #guiding_centres=true <br><br> # Bandstructure plot <br>
 restart         =  plot <br> bands_plot      =  true <br> begin kpoint_path <br> G 0.00000  0.00000 0.0000 Z 0.50000  0.50000 0.5000 <br> Z 0.50000  0.50000 0.5000 F 0.50000  0.50000 0.0000 <br> F 0.50000  0.50000 0.0000 G 0.00000  0.00000 0.0000 <br>
 G 0.00000  0.00000 0.0000 L 0.00000  0.50000 0.0000 <br> end kpoint_path <br> bands_num_points 40 <br> bands_plot_format gnuplot xmgrace <br><br> begin unit_cell_cart <br>      2.1909999999999998   -1.2649999999999999   10.1579999999999995 <br>
      0.0000000000000000    2.5299999999999998   10.1579999999999995 <br>     -2.1909999999999998   -1.2649999999999999   10.1579999999999995 <br> end unit_cell_cart <br><br> begin atoms_frac <br> Bi  0.4002480527838191  0.4002475805055647  0.4002480685003403 <br>
 Bi  0.5997519485504189  0.5997524155539223  0.5997519357215851 <br> Te -0.0000000014683620  0.0000000026871744 -0.0000000028952512 <br> Te  0.2089282289616552  0.2089253259473623  0.2089276850399828 <br>  Te  0.7910717711724687  0.7910746753059763  0.7910723136333432 <br>
 end atoms_frac <br><br> mp_grid =     6     6     6 <br><br> begin kpoints <br>       0.000000000000      0.000000000000      0.000000000000 <br>       0.166666666667      0.000000000000      0.000000000000 <br>       0.333333333333      0.000000000000      0.000000000000 <br>
 ...... <br> end kpoints </div><div> </div><div>Best!</div><div>Hongyan</div></div>