<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
<font style="font-size: 10pt;" size="2" face="Verdana">Dear </font><font style="font-size: 10pt;" size="2" face="Verdana">Dear Wien2k users and dear Prof. Blaha,</font><font style="font-size: 10pt;" size="2" face="Verdana"><br></font><pre><font style="font-size: 10pt;" size="2" face="Verdana">I have tried to rotation x,y by 45 degree around the z axis for a transition metal ion in an </font><font style="font-size: 10pt;" size="2" face="Verdana">octahedron</font><font style="font-size: 10pt;" size="2" face="Verdana"> according<br>to your direction in previous discussion:<br>1. NSYM=1 with </font><font style="" face="Verdana"><font style="font-size: 10pt;" size="2">identit</font></font><font style="" face="Verdana">y in the </font><font style="font-size: 10pt;" size="2" face="Verdana">first operation;<br>2. x kgen with full k-mesh;<br>3. x lapw1 -up/dn;<br>4. change Local Rotation Martrix;<br>5. x lapw2 -qtl -up/dn.<br><br></font><font style="font-size: 10pt;" size="2" face="V!
 erdana"><font style="font-size: 10pt;">It is okey </font></font><font style="font-size: 10pt;" size="2" face="Verdana">for every </font><font style="font-size: 10pt;" size="2" face="Verdana">none spin polarized calculation or MMTOT=0, but </font><font style="font-size: 10pt;" size="2" face="Verdana"><font style="font-size: 10pt;">problem happened</font></font><font style="font-size: 10pt;" size="2"><font style="" face="Verdana"> </font></font><font style="font-size: 10pt;" size="2" face="Verdana"><font style="font-size: 10pt;">as we execute the<br>sp calculation and even </font></font><font style="font-size: 10pt;" size="2" face="Verdana">more problem for the none-zero MMTOT:<br>1. It can </font><font style="font-size: 10pt;" size="2" face="Verdana"><font style="font-size: 10pt;">only</font></font><font style="font-size: 10pt;" size="2" face="Verdana"> </font><font style="font-size: 10pt;" size="2" face="Verdana">be succeed without parallel calculation, such as we could use!
  x lapw1 -up/dn, x lapw2 -qtl -up/dn,<br>but not </font><font style="f
ont-size: 10pt;" size="2" face="Verdana"><font style="font-size: 10pt;">x lapw1 -up/dn -p, x lapw2 -qtl -up/dn</font></font><font style="font-size: 10pt;" size="2" face="Verdana"> -p;</font><font style="font-size: 10pt;" size="2" face="Verdana"><br></font><font style="font-size: 10pt;" size="2" face="Verdana">2. As we check the QTL results after rotation in sp calculation, we found that the quantities is about 2 times of<br>the "right" value;</font><font style="font-size: 10pt;" size="2" face="Verdana"><br></font><font style="font-size: 10pt;" size="2" face="Verdana">3. We think that the local rotation can not change the charge density of the compound, so the x,y rotated by<br>45 degree will show</font><font style="font-size: 10pt;" size="2" face="Verdana"><font style=""> </font></font><font style="font-size: 10pt;" size="2" face="Verdana">an EXCHANGE</font><font style="font-size: 10pt;" size="2" face="Verdana"> in QTL between orbital dXY and dX2Y2, dZ2 should keep the SAME !
 value<br>as before, and the charge density will projected to the new orbital direction of dXZ, dYZ. The results</font><font style="font-size: 10pt;" size="2"><font style="" face="Verdana"> of</font></font><font style="font-size: 10pt;" size="2" face="Verdana"> MMTOT=0<br>coincide with what we analysis, but not for none-zero MMTOT.</font><font style="font-size: 10pt;" size="2"><font style="" face="Verdana"><br></font></font><font style="font-size: 10pt;" size="2" face="Verdana">I tried the same rotation in TiO2 for different MMTOT,<br>it shows the same problem. How could we explain this calculation?</font><font style="font-size: 10pt;" size="2" face="Verdana"> Is there any error in our rotation?<br><br>Below is the QTL results for Ti atom in TiO2 before and after rotation in the two condition of MMTOT=0 and MMTOT=2,<br>the old local rotation matrix is <br>LOCAL ROT MATRIX:  0.7071068 0.7071068 0.0000000<br>                            -0.7071068 0.7071068 0.0000000<br>       !
                       0.0000000 0.0000000 1.0000000<br>and the new one
 is <br></font><font style=""><font style="font-size: 10pt;" size="2" face="Verdana">LOCAL ROT MATRIX:  1.0000000 0.0000000 0.0000000<br>                             0.0000000 1.0000000 0.0000000<br>                             0.0000000 0.0000000 1.0000000</font></font><br><br><font style="font-size: 10pt;" size="2" face="Verdana">S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ<br>MMTOT=0<br></font><font style=""><font style="font-size: 10pt;" size="2" face="Verdana">Before</font></font><br><font style="font-size: 10pt;" size="2" face="Verdana">QTL001: 1.0406 2.9865 0.6504 0.0232 0.9942 0.9941 0.9980 0.1061 0.1543 0.0958 0.2028 0.0915<br>QTL001: 1.0406 2.9865 0.6504 0.0232 0.9942 0.9941 0.9980 0.1061 0.1543 0.0958 0.2028 0.0915<br></font><font style=""><font style="font-size: 10pt;" size="2" face="Verdana">After</font></font><br><font style="font-size: 10pt;" size="2" face="Verdana">QTL001: 2.0813 5.9727 1.3008 0.0465 1.9883 1.9883 1.9961 0.2122 0.1916 0.3085 0.2943 0.2943<br>Q!
 TL001: 1.0407 2.9864 0.6504 0.0233 0.9942 0.9942 0.9981 0.1061 0.0958 0.1543 0.1472 0.1472<br><br>MMTOT=2<br>Before<br>QTL001: 1.0394 2.9814 0.9011 0.0232 0.9924 0.9924 0.9966 0.1547 0.1968 0.2326 0.1751 0.1418<br>QTL001: 1.0380 2.9787 0.5225 0.0209 0.9909 0.9915 0.9964 0.0889 0.1270 0.0650 0.1733 0.0685<br>After<br>QTL001: 2.0768 5.9613 1.0590 0.0459 1.9842 1.9842 1.9928 0.1834 0.1348 0.2586 0.2411 0.2411<br>QTL001: 1.0383 2.9812 0.5307 0.0234 0.9923 0.9923 0.9967 0.0900 0.0683 0.1274 0.1225 0.1225<br><br></font><font style="font-size: 10pt;" size="2" face="Verdana">Thank you in advance,</font><font style="font-size: 10pt;" size="2" face="Verdana"><br></font><font style="font-size: 10pt;" size="2" face="Verdana">Xueli Wang</font><br></pre>                                               </body>
</html>