Dear All,<br><br>i have performed LDA+U calculations and i want to know the occupancies of d states.<br>can i get it from occupation matrix printed in the output file?<br>For example PW gives me:<br><br>atom  4   Tr[ns(na)]=   7.4651394<br>
atom  4  spin  1<br>eigenvalues:  0.4830637 0.5075106 0.5192560 0.6802564 0.6865322<br> eigenvectors<br> 1   0.0000000  0.0000000  1.0000000  0.0000000  0.0000000<br> 2   0.0000000  0.0000000  0.0000000  0.0000000  1.0000000<br>
 3  -0.7160583  0.0000000  0.0000000  0.6980405  0.0000000<br> 4   0.0000000 -1.0000000  0.0000000  0.0000000  0.0000000<br> 5  -0.6980405  0.0000000  0.0000000 -0.7160583  0.0000000<br> occupations<br> 0.601  0.000  0.000  0.084  0.000<br>
 0.000  0.680  0.000  0.000  0.000<br> 0.000  0.000  0.483  0.000  0.000<br> 0.084  0.000  0.000  0.605  0.000<br> 0.000  0.000  0.000  0.000  0.508<br>atom  4  spin  2<br>eigenvalues:  0.8843330 0.9084008 0.9152566 0.9213736 0.9591566<br>
 eigenvectors<br> 1  -0.0108812  0.0000000  0.0000000 -0.9999408  0.0000000<br> 2   0.0000000  0.0000000  0.0000000  0.0000000  1.0000000<br> 3   0.0000000  1.0000000  0.0000000  0.0000000  0.0000000<br> 4   0.0000000  0.0000000 -1.0000000  0.0000000  0.0000000<br>
 5  -0.9999408  0.0000000  0.0000000  0.0108812  0.0000000<br> occupations<br> 0.959  0.000  0.000 -0.001  0.000<br> 0.000  0.915  0.000  0.000  0.000<br> 0.000  0.000  0.921  0.000  0.000<br>-0.001  0.000  0.000  0.884  0.000<br>
 0.000  0.000  0.000  0.000  0.908<br><br>can i say:<br><br>n_down_z2=0.959,<br>n_down_xz=0.915<br>n_down_yz=0.921<br>n_down_x2-y2=0.884<br>n_down_xy=0.908                 ?<br><br><br>thank you.<br><br>Gianluca<br><br><br>
<br>