Dear Ben Thanks for your right answer. The square matrix in oversqr.dat is not the direct overlap matrix in formula C^TST=Identity. After reading your reply, I do a transpose (without conjugate) of matrix obtained from oversqr.dat and then obtain the right answer. Thanks a lot for your help!:-D Fenglin