Dear All,<br><br>i`m a beginner using W90.<br><br>i`m working with FeSe.<br><br>i have the following problem to understand how the flag xaxis works in the file ".win".<br><br>i have a unit cell in the PW simulation with two Fe (o) sites as here (+ and - are Se sites):<br>
<br>Y<br>|    ---------------------<br>|    |   o1          +  |<br>|    |                    |<br>|    |                    |<br>|    |   -           2o |<br>|    ---------------------    <br>
|<br>|------------------------------ X<br><br>X=1,0,0<br>Y=0,1,0<br>(Z=0,0,1 is pointing out of the plane)<br><br>however i want to have printed the H(R) in the file "_hr.dat" such to have Wannier functions of Fe states with internal reference frame directed along the Fe-Fe bonds (then pointing along o1-----2o bonds).<br>
Formally, if we want to maintain that e_g orbitals are x^2 - y^2 and 3z^2 - r^2 and t_2g are xy, yz, zx then the above choice of
coordinate system gives <br>the correct e_g - t_2g splitting coming from
the Se crystal field (e_g go lower in energy and t_2g higher).<br><br>As my hamiltonian must be written in this set, my idea would be to set x'=1,-1,0 in the file ".win".<br><br>The output and the WFs look ok.<br>
<br>here the input file:<br><br>-----------------------------------------------------------------------------------------------------------------------------------------------------------------------<br>num_wann        = 10<br>
num_iter        = 5000<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-16,27-50<br><br>!bands_plot_project : 1,2,3,4,5,6,7,8,9,10<br><br>write_r2mn=.true.<br>
hr_plot=.true.<br><br>kmesh_tol=0.00001<br><br>begin atoms_cart<br>ang<br>Fe   2.82697500  0.94232500  2.74305000<br>Fe   0.94232500  2.82697500  2.74305000<br>Se   0.94232500  0.94232500  4.24459557<br>Se   2.82697500  2.82697500  1.24150443<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<br>
!restart = plot<br><br>guiding_centres = T<br><br>begin projections<br>Fe:dxy;dxz;dyz;dx2-y2;dz2:x=1,-1,0<br>end projections<br><br><br>begin unit_cell_cart<br>bohr<br>7.12294416914163 0.0 0.0<br>0.0 7.12294416914163 0.0<br>
0.0 0.0 10.367225746512<br>end unit_cell_cart<br><br>mp_grid : 12 12 8<br><br>begin kpoints<br>(.....)<br>end kpoints<br><br>-----------------------------------------------------------------------------------------------------------------------------------------------------------------------<br>
<br>is this correct?<br>do i get what i want in this way?<br><br>thank you.<br><br>Gianluca <br><br><br><br> <br>