Hello,<br><br>I am having trouble understanding the output from ph.x in the matdyn file.  I am looking at Barium Titanate in the cubic phase, but with broken symmetry and then relaxed.  <br>I'm trying to compare the linear response frequencies from ph.x to frozen phonon frequencies.  From what I understand, the diagonalized dynamical matrix should give displacements of each atom.  So in the output:<br>
<br>omega( 4) =       3.763989 [THz] =     125.553988 [cm-1]<br> (  0.097362  0.000000 -0.116603  0.000000 -0.288197  0.000000 ) <br> ( -0.066349  0.000000  0.077204  0.000000  0.200853  0.000000 ) <br> ( -0.169282  0.000000  0.204883  0.000000  0.421687  0.000000 ) <br>
 ( -0.159112  0.000000  0.218980  0.000000  0.474827  0.000000 ) <br> ( -0.184271  0.000000  0.191581  0.000000  0.474679  0.000000 ) <br><br>each column gives the real and imaginary components of the displacements in cartesian coordinates.  I could not find what units these values have, first of all.  I assume they are unitless and are a fraction of the lattice parameter.  <br>
Also, each column should add up to zero to preserve symmetry, should it not?  For higher frequencies in this calculation the columns do indeed add up to zero, but for lower frequencies (and negative frequencies) this is not the case.  Does this indicate a problem with my pw calculation or something else? <br>
<br>Is there a resource for good documentation on the matdyn file I have overlooked? <br><br>Thanks for helping,<br><br>Dennis Jackson<br>Oregon State University<br><a href="mailto:denny.jackson@gmail.com">denny.jackson@gmail.com</a><br>