Dear All,<br><br>i`m performing PW+W90 spin polarized calculations.<br><br>I first do a scf run:<br>---------------------------------------------------------------------------------------<br>FeSe<br> &control<br>    calculation     =  'scf'<br>

    restart_mode    =  'from_scratch'<br>    prefix          =  'fese'<br>    tprnfor         =  .true.<br>    pseudo_dir      =  './'<br>    outdir          =  './'<br>    iprint          =   2<br>

/<br> &system<br>    ibrav           =  8<br>    celldm(1)       =  10.0733642480265<br>    celldm(2)       =  1.00000000000000<br>    celldm(3)       =  1.0291721306383<br>    nat             =  8<br>    ntyp            =  5<br>

    ecutwfc         =  25.0<br>    ecutrho         =  300.0<br>    nbnd            =  80<br>    occupations     =  'smearing'<br>    smearing        =  'methfessel-paxton'<br>    degauss         =  0.01<br>

    nspin = 2<br>    starting_magnetization(1) =  1.0<br>    starting_magnetization(2) =  1.0<br>    starting_magnetization(3) = -1.0<br>    starting_magnetization(4) = -1.0<br>    starting_magnetization(5) =  0.0<br>    lda_plus_u      =.true.<br>

    Hubbard_U(1)    =  1.d-8<br>    Hubbard_U(2)    =  1.d-8<br>    Hubbard_U(3)    =  1.d-8<br>    Hubbard_U(4)    =  1.d-8<br>    Hubbard_alpha(1) = 1.d-8<br>    Hubbard_alpha(2) = 1.d-8<br>    Hubbard_alpha(3) = 1.d-8<br>

    Hubbard_alpha(4) = 1.d-8<br>/<br> &electrons<br>    conv_thr        =   1.0d-8<br>    diagonalization =  'cg'<br>/<br>ATOMIC_SPECIES<br> Fe1 55.845   Fe.pbe-sp-van.UPF<br> Fe2 55.845   Fe.pbe-sp-van.UPF<br>

 Fe3 55.845   Fe.pbe-sp-van.UPF<br> Fe4 55.845   Fe.pbe-sp-van.UPF<br> Se 78.96    Se.pbe-van.UPF<br>ATOMIC_POSITIONS { angstrom }<br>Fe1   0.00000000  2.66529759  2.74305000<br>Fe2   0.00000000  0.00000000  2.74305000<br>

Fe3   2.66529759  2.66529759  2.74305000<br>Fe4   2.66529759  0.00000000  2.74305000<br>Se   1.33264880  3.99794639  4.02515157<br>Se   3.99794639  3.99794639  1.46094843<br>Se   3.99794639  1.33264880  4.02515157<br>Se   1.33264880  1.33264880  1.46094843<br>

K_POINTS {automatic}<br>12 12 8 0 0 0<br>---------------------------------------------------------------------------------------<br><br>then a nscf run:<br>---------------------------------------------------------------------------------------<br>

FeSe<br> &control<br>    calculation     =  'nscf'<br>!    restart_mode    =  'from_scratch'<br>    prefix          =  'fese'<br>    tprnfor         =  .true.<br>    pseudo_dir      =  './'<br>

    outdir          =  './'<br>    iprint          =   2<br>    wf_collect = .true.<br>/<br> &system<br>    ibrav           =  8<br>    celldm(1)       =  10.0733642480265<br>    celldm(2)       =  1.00000000000000<br>

    celldm(3)       =  1.0291721306383<br>    nat             =  8<br>    ntyp            =  5<br>    ecutwfc         =  25.0<br>    ecutrho         =  300.0<br>    nbnd            =  80<br>    occupations     =  'smearing'<br>

    smearing        =  'methfessel-paxton'<br>    degauss         =  0.01<br>    nspin = 2<br>    starting_magnetization(1) =  1.0<br>    starting_magnetization(2) =  1.0<br>    starting_magnetization(3) = -1.0<br>

    starting_magnetization(4) = -1.0<br>    starting_magnetization(5) =  0.0<br>    lda_plus_u      =.true.<br>    Hubbard_U(1)    =  1.d-8<br>    Hubbard_U(2)    =  1.d-8<br>    Hubbard_U(3)    =  1.d-8<br>    Hubbard_U(4)    =  1.d-8<br>

    Hubbard_alpha(1) = 1.d-8<br>    Hubbard_alpha(2) = 1.d-8<br>    Hubbard_alpha(3) = 1.d-8<br>    Hubbard_alpha(4) = 1.d-8<br>    nosym=.true.<br>    noinv=.true.<br>/<br> &electrons<br>    conv_thr        =   1.0d-8<br>

    diagonalization =  'cg'<br>/<br>ATOMIC_SPECIES<br> Fe1 55.845   Fe.pbe-sp-van.UPF<br> Fe2 55.845   Fe.pbe-sp-van.UPF<br> Fe3 55.845   Fe.pbe-sp-van.UPF<br> Fe4 55.845   Fe.pbe-sp-van.UPF<br> Se 78.96    Se.pbe-van.UPF<br>

ATOMIC_POSITIONS { angstrom }<br>Fe1   0.00000000  2.66529759  2.74305000<br>Fe2   0.00000000  0.00000000  2.74305000<br>Fe3   2.66529759  2.66529759  2.74305000<br>Fe4   2.66529759  0.00000000  2.74305000<br>Se   1.33264880  3.99794639  4.02515157<br>

Se   3.99794639  3.99794639  1.46094843<br>Se   3.99794639  1.33264880  4.02515157<br>Se   1.33264880  1.33264880  1.46094843<br>K_POINTS crystal<br>1152<br>[...here i removed the list....]<br>---------------------------------------------------------------------------------------<br>

<br>then i start with Wannier90.<br>Let look only at the spin channel "up".<br><br>I run at first:<br><br>/work/aha/giovanne/ESPRESSO/espresso-4.0.5/W90/wannier90.x -pp  FeSe_up<br><br>the FeSe_up input file looks like:<br>

---------------------------------------------------------------------------------------<br>num_wann        = 20<br>num_iter        = 5000<br><br>dis_num_iter    =  500<br><br>iprint = 3<br><br>exclude_bands : 1-32,53-80<br>

<br>write_r2mn=.true.<br>hr_plot=.true.<br><br>kmesh_tol=0.00001<br><br>spin = up<br><br>begin atoms_cart<br>ang<br>Fe   0.00000000  2.66529759  2.74305000<br>Fe   0.00000000  0.00000000  2.74305000<br>Fe   2.66529759  2.66529759  2.74305000<br>

Fe   2.66529759  0.00000000  2.74305000<br>Se   1.33264880  3.99794639  4.02515157<br>Se   3.99794639  3.99794639  1.46094843<br>Se   3.99794639  1.33264880  4.02515157<br>Se   1.33264880  1.33264880  1.46094843<br>end atoms_cart<br>

<br>begin kpoint_path<br>G 0.0000 0.0000 0.000 X 0.5000 0.0000 0.0000<br>X 0.5000 0.0000 0.000 M 0.5000 0.5000 0.000<br>M 0.5000 0.5000 0.000 Y 0.0000 0.5000 0.0000<br>Y 0.0000 0.5000 0.000 G 0.0000 0.0000 0.0000<br>G 0.0000 0.0000 0.000 Z 0.0000 0.0000 0.5000<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 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20<br>

!restart = plot<br><br>guiding_centres = T<br><br>begin projections<br>Fe:dxy;dxz;dyz;dx2-y2;dz2<br>end projections<br><br><br>begin unit_cell_cart<br>bohr<br>10.0733642480265 0.0000 0.0000<br>0.0000 10.0733642480265 0.0000<br>

0.0000 0.0000   10.3672257458371<br>end unit_cell_cart<br><br>mp_grid : 12 12 8<br><br>begin kpoints<br>[....]<br>end kpoints<br>---------------------------------------------------------------------------------------<br>
then:<br>
<br>mpirun -np 8 /work/aha/giovanne/ESPRESSO/espresso-4.0.5/bin/pw2wannier90.x < FeSe_up.pw2wan<br><br>with FeSe_up.pw2wan:<br>---------------------------------------------------------------------------------------<br>

&inputpp<br>   outdir = './'<br>   prefix = 'fese'<br>   seedname = 'FeSe_up'<br>   spin_component = 'up'<br>   write_mmn = .true.<br>   write_amn = .true.<br>   write_unk = .true.<br>
/<br>
---------------------------------------------------------------------------------------<br>once i re-run W90 as:<br><br>/work/aha/giovanne/ESPRESSO/espresso-4.0.5/W90/wannier90.x FeSe_up<br><br>i get the following error:<br>

<br>---------------------------------------------------------------------------------------<br>Found a mismatch in FeSe_up.eig<br>Wanted band  : 1 found band  : 21<br>Wanted kpoint: 2 found kpoint: 1<br><br>A common cause of this error is using the wrong<br>

number of bands. Check your input files.<br>If your pseudopotentials have shallow core states remember<br>to account for these electrons.<br><br> Exiting.......<br> param_read: mismatch in FeSe_up.eig<br>---------------------------------------------------------------------------------------<br>

<br>i think i should have 20 eigenvalues in FeSe_up.eig but they are 40.<br><br>Could you please help me to understand it?<br><br>thank you.<br><font color="#888888"><br>Gianluca</font>