Dear Jonathan (and All),<br>
<br>
as i understand the flag xaxis helps to converge the WFs function to
the correct one but it doesn`t change the local and internal  frame
which is set from PW calculations.<br>
<br>
I had a look to the paper suggested.<br>
I get Wannier functions which have correct shape and they are the same starting from different starting guess.<br>
probably i`m in the correct minimum.<br>
<br>
Then the only solution is to rotate the unit cell in PW calculations
such to have Wannier functions of Fe states with internal reference
frame directed along the Fe-Fe bonds (then pointing along o1-----2o
bonds) from:<br>
<br>
Y<br>

|    ---------------------<br>

|    |   o1          +  |<br>

|    |                    |<br>

|    |                    |<br>

|    |   -           2o |<br>

|    ---------------------    <br>


|<br>

|------------------------------ X<br>
 <br>
to (rotating the unit cell of 45 degree):<br>
<br>
<br>
Y<br>
|    <br>
|                            +<br>
|    <br>
|                  o1               2o<br>
|          <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>Then i expect the Hlda(R)_n,m is written in the correct local frame. <br>
Is this way correct?<br>
<br>
thank you.<br>
<br>
Gianluca<br><br><br><br><div class="gmail_quote">On Sat, Feb 13, 2010 at 12:38 AM, Jonathan Yates <span dir="ltr"><<a href="mailto:jonathan.yates@materials.ox.ac.uk">jonathan.yates@materials.ox.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im"><br>
On 12 Feb 2010, at 15:46, Gianluca Giovannetti wrote:<br>
> however i want to have printed the H(R) in the file "_hr.dat" such<br>
> to have Wannier functions of Fe states with internal reference frame<br>
> 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<br>
> 3z^2 - r^2 and t_2g are xy, yz, zx then the above choice of<br>
> coordinate system gives<br>
> the correct e_g - t_2g splitting coming from the Se crystal field<br>
> (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<br>
> set x'=1,-1,0 in the file ".win".<br>
<br>
<br>
</div>The global x,y,z axes are set by your pw code. Wannier90 can use a<br>
projection of the bloch states onto a local function times spherical<br>
harmonics as an initial guesses for the MLWF. For example a p_z<br>
function. You might want the local coordinate system of this function<br>
to be different from the global x,y,z - in otherwords to have the p_z<br>
point in an arbitary direction. For this reason we allowed the<br>
possibility to specify the orientation for each projection.<br>
  So what you did to align your starting guess functions is sound.<br>
<br>
However, these functions are just starting guesses. The MLWF should be<br>
independent of these functions. No symmetry is imposed during the<br>
minimisation. Maybe the true MLWF do have the symmetry you want - in<br>
any case you can't automatically enforce that they have this symmetry.<br>
  In practice the landscape can be complex, and you might find picking<br>
a 'good' starting guess helps to get to the true minimum. (if two<br>
different initial projections lead to different sets of WF, you can at<br>
least compare the spreads).<br>
<br>
It might be worth looking at<br>
Phys. Rev. B 65, 184422 (2002)<br>
Maximally localized Wannier functions in antiferromagnetic MnO within<br>
the FLAPW formalism<br>
  In this situation the minimizing the spread was not sufficient to<br>
uniquely define the MLWF - and a further step was used to obtain MLWF<br>
with the desired symmetry.<br>
<br>
<br>
  Jonathan<br>
<br>
<br>
--<br>
Department of Materials, University of Oxford, Parks Road, Oxford, OX1<br>
3PH, UK<br>
tel: +44 (0)1865 612797                <a href="http://users.ox.ac.uk/%7Eoums0549/" target="_blank">http://users.ox.ac.uk/~oums0549/</a><br>
<br>
_______________________________________________<br>
Wannier mailing list<br>
<a href="mailto:Wannier@quantum-espresso.org">Wannier@quantum-espresso.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/wannier" target="_blank">http://www.democritos.it/mailman/listinfo/wannier</a><br>
</blockquote></div><br>