<P>Dear all:<BR>I'm using QE-5.0.1 to calculate phonons of a hypothetical crystal.<BR>It is known to all that an imaginary phonon frequency at non-Gamma point means<BR>the crystal structure is not stable. My question is how to adjust the crystal<BR>structure to eliminate the imaginary phonon frequency? My method is as follows. </P>
<P>(1) calculate the dynamical matrix at the q-point which has imaginary phonon frequency;<BR>(2) diagonalize the dynamical matrix by dynmat.x;<BR>(3) plot the imaginary phonon modes by Xcrysden to see the movements of atoms;<BR>(4) adjust the inner-coordinations of atoms towards the direction of motion,<BR>corresponding to the imaginary phonon modes. Relax this new structure and recalculate<BR>phonons to see whether there is some improvement or not.</P>
<P>Is the procedures correct? But after relaxation, atoms do not form a stable structure and<BR>go back to the initial positions, which has imaginary phonon frequency. Any suggestions?<BR>Any help will be appreciated.</P>
<P>Best</P>
<P>Miao Gao</P>
<P>---------------------------------------<BR>Ph.D. student of Department of Physics, <BR>Renmin University of China.<BR>---------------------------------------<BR><BR></P><SPAN></SPAN><br><br>