<div>Dear  Prof Nicola,</div>

<div> </div>

<div>Thanks for your long suggestion.</div>

<div> </div>

<div>I may not follow some of your words, would you please clarify  ?</div>

<div> 1,  How to relax the ions more ?    </div>

<div> 2, I knew from most refs that ecutwfs  30 Ry for ultrasoft-PP is enough in pw calculations , Is it right ? But i donot have much idear about ecutwfs in ph calculations .  Would you please give me some refs of more detail settings about PH calculations ?  

</div>

<div> </div>

<div> </div>

<div>I am trying to find more problems related to this wrong result following your suggestion.</div>

<div>Thank you again, it really helps me much.</div>

<div> </div>

<div>p.s.  It is nice for me listening  your lectures in Beijing , i really learn much from your 3 lectures . </div>

<div> </div>

<div> </div>

<div>Regards,</div>

<div> </div>

<div>Hai-Ping<br><br> </div>

<div><span class="gmail_quote">On 8/20/06, <b class="gmail_sendername">Nicola Marzari</b> <<a href="mailto:marzari@mit.edu">marzari@mit.edu</a>> wrote:</span>

<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><br><br><br>Dear Haiping,<br><br><br>there are two possibilities that lead to negative phonon frequencies at

<br>gamma.<br><br>The first one is less likely:<br><br>1) you have not carefully relaxed the *internal* degrees of freedom -<br>i.e. the ionic coodinates. Then, the phonon code discovers that moving<br>the atoms along a certain pattern (an "eigenvector" of the dynamical

<br>matrix at gamma) the energy goes down. Perfect - just relax the ions<br>more. Relaxation of the unit cell has nothing to do with this<br><br>The second one is more likely:<br><br>2) the calculation is not converged with respect

<br>to the plane wave cutoff needed, and the number of k-points needed,<br>for the small or zero frequencies at gamma. Obtaining an accurate,<br>faithful description of low energy modes requires cutoffs and especially<br>

duals much higher than anything you are used to (to give you an example,<br>for an ultrasoft C for diamond, you need 100 Ry in the wavefunction and<br>2800 Ry in the charge density to get a convergence, at gamma,<br>of the zero modes down to the last cm-1). Does it mean you need

<br>to use 100/2800 Ry for diamond ? No ! By enforcing the translation sum<br>rules that the interatomic force constants must satisfy, you can monitor<br>(after enforcing said rules) the convergence of the low frequency modes,

<br>and discover you can use much more reasonable cutoffs.<br><br>Still, even more than a regular calculations, phonons require a<br>through control on the parameters that define your calculation.<br><br>Note that, as a warm-up, it is easier to converge LDA calculations

than GGA calculations (the exchange-correlation terms, calculated in real space, introduce very high frequencies in reciprocal space).                        nicola >> Dear All, >> I also have calculated gamma phonon frequencies of anatase TiO2 .

>> I run ph.x caculation after relaxing the structure, but  9 of 18 frequencies >> are minus. >> >> Would you please give me some hints and help? >> >> input settings are below, thank you in advance !

<br>>><br>>> *phonons of TiO2 at Gamma<br>>>  &inputph<br>>>   tr2_ph=1.0d-14,<br>>>   recover=.true.,<br>>>   epsil=.true.,<br>>>   prefix='tio2',<br>>>   fildyn='TiO2-gamma.dyn

',<br>>>   outdir='/home/haiping/tmp/'<br>>>  /<br>>> 0.0 0.0 0.0*<br><br><br>---------------------------------------------------------------------<br>Prof Nicola Marzari   Department of Materials Science and Engineering

<br>13-5066   MIT   77 Massachusetts Avenue   Cambridge MA 02139-4307 USA<br>tel 617.4522758 fax 2586534 <a href="mailto:marzari@mit.edu">marzari@mit.edu</a> <a href="http://quasiamore.mit.edu">http://quasiamore.mit.edu</a>

<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>