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><br>Gianluca Giovannetti<br>