Dear All,<br><br>i`m calculating Wannier functions for a molecular crystal (sites are molecule with C and H ions within each molecule).<br><br>i want to plot Wannier functions of the 4 bands at Fermi level.<br>the system is metallic.<br>
these four bands come from HOMO and LUMO molecular states of the two inequivalent molecules in the unit cell.<br><br>here it is my input file:<br><br>--------------------------------<br><br>num_wann        = 4<br>num_iter        = 2000<br>
<br>dis_num_iter    =  500<br><br>iprint = 3<br><br>!here we exclude all bands except the O2p bands<br>exclude_bands : 1-126,131-180<br><br>!bands_plot_project : 1,2,3,4<br><br>write_r2mn=.true.<br>hr_plot=.true.<br><br>kmesh_tol=0.00001<br>
<br>begin atoms_cart<br>ang<br>[......]<br>end atoms_cart<br><br>begin kpoint_path<br>G 0.00000000 0.00000000 0.00000000    M  0.50000000 0.00000000 0.00000000<br>M 0.50000000 0.00000000 0.00000000    X  0.50000000 0.50000000 0.00000000<br>
X 0.50000000 0.50000000 0.00000000    G  0.50000000 0.00000000 0.50000000<br>G 0.50000000 0.00000000 0.50000000    Z  0.00000000 0.00000000 0.00000000<br>Z 0.00000000 0.00000000 0.00000000    R  0.50000000 0.50000000 0.00000000<br>
R 0.50000000 0.50000000 0.00000000    A  0.50000000 0.00000000 0.50000000<br>A 0.50000000 0.00000000 0.50000000    Z  0.50000000 0.50000000 0.50000000<br>Z 0.50000000 0.50000000 0.50000000    C  0.50000000 0.00000000 0.00000000<br>
C 0.50000000 0.00000000 0.00000000    F  0.50000000 0.50000000 0.00000000<br>end kpoint_path<br>bands_plot =T<br><br>!fermi_surface_plot = .true.<br><br>!To plot the WF<br>!wannier_plot = T<br>!wannier_plot_supercell = 2<br>
!wannier_plot_list = 1 2 3 4<br>!restart = plot<br><br>guiding_centres = T<br><br>begin projections<br>random<br>end projections<br><br><br>begin unit_cell_cart<br>bohr<br>16.4538 0.000 0.000<br>0.000 0.000 -11.172<br>-1.18447 24.481 0.000<br>
end unit_cell_cart<br><br>mp_grid : 8 8 8<br><br>begin kpoints<br>[...]<br>end kpoints<br><br><br>--------------------------------<br><br>However i get Wannier eigenvalues which are not exactly the same as the one i have calculated Quantum espresso (the discrepancy is not big but finite at some k-points).<br>
<br>My feeling is that the setting in the win file:<br><br>begin projections<br>
random<br>
end projections<br><br>is not really good as starting point.<br><br>Indeed the CONV is 1.0E-03.<br><br>HOMO and LUMO are localized on the molecules.<br><br>What is the best strategy in this case?<br><br>thank you in advance.<br>
<br>Gianluca<br><br><br>