<div class="gmail_quote">
<div> </div>
<div>Dear Yin:</div>
<div> </div>
<div> </div>
<div>Thank you for your answers. Sorry, I didn't make myself clear.</div>
<div> </div>
<div>For A1:</div>
<div> </div>
<div>QHA is what I want. Thank you.</div>
<div> </div>
<div>For A2:</div>
<div>The soft modes are not from numerical error.  It is high temperature phase.</div>
<div> </div>
<div>For A3:</div>
<div>It is not what I want.  The discontinuity happens when two phonon bands cross. The phonon eigenvalues in the output file xxx.freq are always arranged ascendingly. When ploting with gnuplot, the cross two bands turn to two ugly separate lines. Another represention, how do you plot just a single phonon band among the phonon dispersion ?</div>

<div> </div>
<div>Zhou  Bo</div>
<div>2009.8.9</div>
<div> </div>
<div> </div></div>