Dear All,<br><br>i`m working on extracting the Wannier functions of LiV2O4.<br><span style="font-family:arial,helvetica,sans-serif"><font color="#05004F">The trigonal symmetry 
splits the low-lying t_{2g} levels into one  a_{1g} and two degenerate 
e_{g}^pi levels, with the former slightly lower in energy and narrower 
in bandwidth.</font></span> The density of states looks very simple. <br>At Fermi level there are Vanadium d bands having a bandwidth of ~ 2eV.<br>Below in energy separated from the Vanadium d bands there are bands with predominant Oxygen p bands.<br>
I can extract with no problem the Vanadium d bands using an energy window sitting at the  Vanadium d bands.<br>Wannier90 bands are on top of DFT bands (i also compared mine with some known in literature).<br>As a starting point for the projection i used:<br>
<br>begin projections<br>V:dxz;dyz;dxy<br>end projections<br>Once i try to calculate Wannierfunctions for "d" and "p" orbitals i get the following error:<br><br>   A_mn = <psi_m|g_n> --> S = A.A^+ --> U = S^-1/2.A<br>
   In dis_project...  ERROR: unitarity of initial U<br> nkp=  1<br> i=  2  j=  4<br> [u_matrix_opt.transpose(u_matrix_opt)]_ij=     0.114936    -0.050292<br> Exiting.......<br> dis_project: Error in unitarity of initial U in dis_project<br>
<br>Do you have any idea for solving it?<br>Below it is my input file.<br><br>thank you in advance.<br><br>Gianluca Giovannetti<br><br>=========================<br><br>iprint = 3<br>num_bands       = 70<br>num_wann        = 36<br>
<br>num_iter        =  2000<br>num_print_cycles = 100<br><br>dis_win_max     = 9.1<br>dis_win_min     = -0.5<br><br>dis_num_iter    = 200<br>dis_mix_ratio   = 0.75<br><br>write_r2mn=.true.<br>hr_plot=.true.<br><br>kmesh_tol=0.00001<br>
<br>! SYSTEM<br><br>begin unit_cell_cart<br>bohr<br>7.789167066 0.0000 7.789167066<br>0.0000 7.789167066 7.789167066<br>7.789167066 7.789167066 0.0000<br>end unit_cell_cart<br><br>begin atoms_cart<br>Ang<br>V    0.000000    0.000000    0.000000<br>
V    2.060925    2.060925    0.000000<br>V    2.060925    0.000000    2.060925<br>V    0.000000    2.060925    2.060925<br>O     0.094555000   0.094555000  1.966370000<br>O     1.966370000   1.966370000  1.966370000<br>O     0.094555000   1.966370000  0.094555000<br>
O     1.966370000   0.094555000  0.094555000<br>O    -0.094556000   4.027295000  2.155480000<br>O     2.155480000  -1.966370000  2.155480000<br>O    -0.094556000  -1.966370000 -0.094556000<br>O     2.155480000   4.027295000 -0.094556000<br>
Li     1.030463000   5.152312000  1.030463000<br>Li     3.091388000   3.091388000  3.091388000<br>end atoms_cart<br><br>begin projections<br>V:dxz;dyz;dxy<br>O:px;py;pz<br>end projections<br><br>begin kpoint_path<br>G 0.00  0.00  0.00    X 0.00  0.50  0.50<br>
X 0.00  0.50  0.50    W 0.25  0.75  0.50<br>W 0.25  0.75  0.50    L 0.00  0.50  0.00<br>L 0.00  0.50  0.00    G 0.00  0.00  0.00<br>end kpoint_path<br>bands_plot = .true.<br><br>!To plot the WF<br>!wannier_plot = T<br>!wannier_plot_supercell = 2<br>
!wannier_plot_list = 1-12<br>!restart = plot<br><br>! KPOINTS<br><br>mp_grid : 8 8 8<br><br>begin kpoints<br>.............<br>