I did a test in phonon.f90. The code is sth like this:<br><br>  USE ions_base, only : nat,amass<br>  complex(dp), allocatable :: work(:,:)<br>  allocate(work(3*nat,3*nat))<br>  <br>  DO ipert = 1,3*nat<br>     DO jpert = 1,3*nat

<br>        DO i=1,3*nat<br>        work(ipert,jpert) = work(ipert,jpert) + conjg(dyn(i,ipert))*amass(1)*dyn(i,jpert)<br>        ENDDO<br>     ENDDO<br>  ENDDO<br>  write(*,*) 'dynmat:'<br>  write(*,"(6(f12.8))") work(:,:)

<br>  deallocate(work)<br>------------------------------------------<br>The output for example07 is:<br>  1.00000000  0.00000000  0.00000000  0.00000000  0.00000000  0.00000000<br>  0.00000000  0.00000000  1.00000000  0.00000000

  0.00000000  0.00000000<br>  0.00000000  0.00000000  0.00000000  0.00000000  1.00000000  0.00000000<br> The relation dyn^{\diag} * amass * dyn = I sustains indeed.<br> Actually we're talking about the same thing. The inner product of eigenvector U^{\diag} M U gives scalar 1.

<br> While dyn^{\diag} * amass * dyn gives unit matrix.<br><br><div><span class="gmail_quote">On 12/1/06, <b class="gmail_sendername">degironc</b> <<a href="mailto:degironc@sissa.it">degironc@sissa.it</a>> wrote:</span>

<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Jin Zhang wrote:<br><br>> 2. dyn(iat, jat): The dynamical matrix which obeys -- dyn^{\diag} *

> amass * dyn = 1 I don't think dymanical matric should obey this relation... EIGENVECTORS,  U,  of the generalized problem Dyn U = M U omega^2  will obey U^\dagger M U = 1 or something like that...

<br>while<br>EIGENVECTORS, V=\sqrtM U of problem 1/sqrtM Dyn 1/sqrtM V = omega^2 V<br>will obey V^\dagger V = 1<br><br>stefano<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">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br>Dep. of Physics, Peking University

<br>Beijing, P.R.China<br>Phone: 86-10-62-753-190