<DIV>
<DIV>In order (to try) to explain both the written expression in the usersguide, "The EFG tensor is given in both, the ``local-rotation-matrix'' coordinate system, and then diagonalized. The resulting eigenvectors of this rotation are given by columns.", which is in brief taking into account the columns and the programmed expression in the efg.f given by Peter, </DIV>
<DIV>"anorm=sqrt(eivec(1,1)**2+eivec(1,2)**2+eivec(1,3)**2)<BR>winkel=acos(eivec(1,1)/anorm) * 180.d0 / acos(-1.d0)<BR>WRITE(6,1355) ((eivec(i,j),j=1,3),i=1,3),jatom,winkel"</DIV>
<DIV>which is in brief taking into account the rows instead of columns are correct, let us back to the active and passive transformations.</DIV>
<DIV>We know that in an active transformation the coordinate system is fixed and this is the vector that rotates applying the rotation matrix.</DIV>
<DIV>However, in a passive transformation the vector is fixed and this is the coordinate system that rotates.</DIV>
<DIV>Now one can interpret the passive transformation as an active transformation taking transpose of the matrix into account which is identical to exchanging the columns and rows.<BR></DIV>
<DIV> </DIV>
<DIV>So it looks that nothing is wrong with both the usresguied and the program efg.f.</DIV>
<DIV> </DIV>
<DIV>Obviously taking x = (1,0,0) and x' = (1.0000 0.2760 -0.3709) and then calculating winkel = arccos(x . x'/|x|*|x'|)* 180.d0 / acos(-1.d0) gives the expected 24.8 result. </DIV>
<DIV> </DIV>
<DIV>I will be thankful to anyone telling us if something is wrong here.</DIV>
<DIV>Your,</DIV>
<DIV>Saeid Jalali.</DIV>
<DIV><BR><B><I>Peter Blaha <pblaha@zeus.theochem.tuwien.ac.at></I></B> wrote:</DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">I'm sorry, at the moment I cannot clarify what is ok and what is wrong:<BR><BR>The UG says it is a column vector; however, in the code (efg.f) the<BR>angle is calculated from a "row":<BR><BR>anorm=sqrt(eivec(1,1)**2+eivec(1,2)**2+eivec(1,3)**2)<BR>winkel=acos(eivec(1,1)/anorm) * 180.d0 / acos(-1.d0)<BR>WRITE(6,1355) ((eivec(i,j),j=1,3),i=1,3),jatom,winkel<BR><BR>It could well be that "winkel" was programmed for a case where the EFG was<BR>fixed along z and only x,y could rotate (in this case column or row does<BR>not matter).<BR><BR>Please check it! (Take part of efg.f ; initialize mat with "hand" selected <BR>values and diagonalize with eigen3 and check the eigenvectors/angles.<BR><BR>I'd be gratefull if you could confirm which is right and which is wrong.<BR><BR>Regards<BR>> <BR>> Dear Wien2k users,<BR>> <BR>> I'm trying to calculate EFG and have some questions !
<BR>>
concerning information about EFG in the case.scf_mini file. <BR>> A typical entry looks like below:<BR>> <BR>> MAIN DIRECTIONS OF THE EFG 1.0000 0.2760 -0.3709 <BR>> 0.4072 -0.5385 1.0000 <BR>> -0.0567 1.0000 0.6408 <BR>> :ANG014: ANGLE WITH OLD X-AXIS = 24.8 <BR>> <BR>> :ETA014: ASYMM. ETA = 0.56833<BR>> <BR>> <BR>> I found in the usersguide that, this gives the eigenvectors of <BR>> EFG (in colums). If I'm not wrong this the axes of the "EFG cordinate <BR>> system".<BR>> Thats why if the first column will be the x' axis (x axis in the <BR>> EFG system) the angle with the old x axis could be calculated as <BR>> follows:<BR>> alpha = arccos(x . x'/|x|*|x'|). However when I take <BR>> x = (1,0,0) and x' = (1.000, 0.4072, -0.0567)<BR>> I get alpha = 22.3 instead of 24.8. <BR>> <BR>> I've checked this for several results and I have never get<BR>> exact angle shown in *.scf_mini. <BR>> <BR>> I will be appre!
ciated if
somebody could explain me how the angle <BR>> between the axes is calculated, or if there is an error in <BR>> my impretation.<BR>> <BR>> Greetings<BR>> Daniel Kmiec<BR>> <BR>> -- <BR>> ______________________________________________________________________<BR>> Mag. Daniel Kmiec http://www.ap.univie.ac.at/users/kmiec/<BR>> <BR>> _______________________________________________<BR>> Wien mailing list<BR>> Wien@zeus.theochem.tuwien.ac.at<BR>> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<BR>> <BR><BR><BR>P.Blaha<BR>--------------------------------------------------------------------------<BR>Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<BR>Phone: +43-1-58801-15671 FAX: +43-1-58801-15698<BR>Email: blaha@theochem.tuwien.ac.at WWW: http://info.tuwien.ac.at/theochem/<BR>--------------------------------------------------------------------------<BR><BR>_______________________________________________<BR>Wien !
mailing
list<BR>Wien@zeus.theochem.tuwien.ac.at<BR>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<BR></BLOCKQUOTE></DIV><p>__________________________________________________<br>Do You Yahoo!?<br>Tired of spam? Yahoo! Mail has the best spam protection around <br>http://mail.yahoo.com