Hi Gianluca,<br><br>I think your interpretation is correct. In fact the occupation matrix is almost diagonal...<br><br>Matteo<br><br><div class="gmail_quote">On Mon, Jul 26, 2010 at 10:34 AM, Gianluca Giovannetti <span dir="ltr"><<a href="mailto:gianluca.giovannetti@gmail.com">gianluca.giovannetti@gmail.com</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;">Dear All,<br><br>i have made self-consistent for my system (FeSe).<br><br>Printing out the density matrix in LDA+U within U=0 at Fe sites i believe i can get the occupancies of d atomic orbitals.<br>
<br>For example in the output i get (non spin polarized calculations):<br>
-----------------------------------------------------------------------------------------------------------<br>atom  4  spin  1<br>eigenvalues:  0.7294836 0.7306191 0.7326025 0.7629065 0.8283887<br> eigenvectors<br> 1   0.0000000  0.0000000 -1.0000000  0.0000000  0.0000000<br>

 2   0.0000000 -1.0000000  0.0000000  0.0000000  0.0000000<br> 3  -0.0087135  0.0000000  0.0000000  0.9999620  0.0000000<br> 4   0.0000000  0.0000000  0.0000000  0.0000000 -1.0000000<br> 5  -0.9999620  0.0000000  0.0000000 -0.0087135  0.0000000<br>

 occupations<br> 0.828  0.000  0.000  0.001  0.000<br> 0.000  0.731  0.000  0.000  0.000<br> 0.000  0.000  0.729  0.000  0.000<br> 0.001  0.000  0.000  0.733  0.000<br> 0.000  0.000  0.000  0.000  0.763<br>-----------------------------------------------------------------------------------------------------------<br>

<br>looking at the eigenvalues i can say for the occupancies n of d orbital states:<br><br>   z2-r2         xz            yz             x2-y2        xy<br>n 0.8283887 0.7306191 0.7294836 0.7326025 0.7629065<br><br>is this correct?<br>

<br>thank you in advance.<br><font color="#888888"><br>Gianluca Giovannetti<br>
</font><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br>