<div dir="ltr"><div dir="ltr"><div><div><div><div>Dear Dr. Blaha,<br></div>For some reason, my message below did not get posted.<br></div><br>Would you please kindly check it to see whether my email address is in your database?<br></div>Thanks!<br><br></div>Guoping<br>-------------------------------------------<br><div>Dear Prof. Blaha and Wien2k users,<br><br>Greetings!<br><br>I have a short question on the density of states decomposition for<br>CrO2 with spin polarized calculation. It has rutile TiO2 structure,<br>but with a small distortion along the apex oxygen atom.<br><br>My case.qtlup file has the following heading<br><br> LATTICE CONST.=  8.3560  8.3560  5.5146   FERMI ENERGY=   0.50182<br> 1048 < NMAT < 1090   SPIN=2   NAT=  2      SO 0<br> JATOM  1  MULT= 2  ISPLIT= 8  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3  <br> JATOM  2  MULT= 4  ISPLIT= 8  tot,0,1,PX,PY,PZ,2,DZ2,DX2Y2,DXY,DXZ,DYZ,3  <br><br>According to the UG, "The m-decomposed DOS (e.g. pz , py , px ) is<br>given with respect to the local coordinate system", but when I plot<br>DOS of DZ2 for Cr (atom 1) and compare it with Fig. 4 of PRL Vol. 80,<br>4305 (1998), Dz2 looks more like their Dxy. Then I check the rotation<br>matrix of Cr atom (which is the same as in the Appendix of UG) to see<br>whether the decomposition is for the rotated atom, but the rotation is<br>45 degree rotation with respect to the $c$ axis, so there is no way<br>Dz2 could become Dxy.<br><br>Therefore, I wonder whether WIEN2k uses a different local coordinate<br>from the above literature where the z axis is from the Cr at the body<br>center to the apex oxygen.<br><br>I have searched the mail list, but could not find an easy answer. I<br>also tried to simply rotate d orbitals using 5x5 rotation matrix, but<br>I have to know which axis orbitals DZ2,DX2Y2,DXY,DXZ,DYZ use.<br><br>Thanks a lot in advance!<br><br>Best regards,<br><br>Guoping<br><br></div></div></div>