<div dir="ltr"><br><div class="gmail_quote"><div dir="ltr"><div>Dears all:</div><div> Recently I have calculated density matrix for bcc-iron using WIEN2k. In case.dmatup for atom 1 and l=2, there should be a 5 by 5 matrix,</div><div><br></div><div>with m labelling the column and row. As I know, there should not be nondiagonal element for when we integrate the density matrix within the atomic sphere of atom 1 as the spherical harmonics are othorgonal sets,i.e.,integral corresponding to different m's should be zero. However, I get nonzero matrix elements in the position(-2,2) and (2,-2) from the case.dmatup. Maybe something is wrong within my understanding of the Bloch eigenfunction?</div><div><br></div><div> The eigenfunction can be partitioned to two parts : one is in the intersitial region one is in the atoms spheres. I just consider that in the atomic spheres, which can be a linear combination of spherical harmonics. So is there something wrong?</div><div><br></div><div> Best regards!</div><div><br></div><div> Feng Tang</div><div><br></div><div> From Nanjing University, China</div><div><br></div></div>
</div><br></div>