<HTML><BODY style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><BR><DIV><DIV>On Aug 20, 2006, at 3:15 AM, lan haiping wrote:</DIV><BR class="Apple-interchange-newline"><BLOCKQUOTE type="cite"><DIV>Dear  Prof Nicola,</DIV></BLOCKQUOTE><DIV><BR class="khtml-block-placeholder"></DIV>just a note about the "western" usage of names. (I would be glad to hear an explanation on how names are used in Chinese!). "Nicola" is the "given" name, i.e. it is the name that parents give to the newborn baby. It is also called "first" name because in written it is the first name that must be used. The other name "Marzari" is the family name which is transmitted from one generation to the next in a patrilinear way (a baby gets the fathers' family name by default). It is a form of confidentiality to address somebody by his/her first name. It is perfect here because: i) it is common usage on the web and; ii) it is common practice among scientists (particularly in the US). However, addressing somebody by his/her title is _very_ formal. Not very appropriate in this forum, but above all kind of funny if accompanied by the much more confidential usage of the first name. Resuming: "Dear Nicola" is confidential and "Dear Prof. Marzari" is formal - both ar OK, the first more appropriate to this forum. "Dear Prof Nicola" is just funny! ... I was not joking saying that I am curious about _your_ usage of "given", "family", or whatever names. If you have a  few spare minutes, just let me know ...</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>Science, now!</DIV><DIV><BR><BLOCKQUOTE type="cite"> <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 ?<BR></DIV></BLOCKQUOTE><DIV><BR class="khtml-block-placeholder"></DIV>The code must have some default threshold for the atomic forces (and, maybe, stress too if you also relax the unit cell) such that  when forces are smaller, it stops. The threshold is usually more than adequate for most practical purposes, _but to learn_ What I mean is that it is a good idea to practice and to see _with one's own eyes_ how a given property (phonon frequencies, in the present case) depend on external parameter and/or convergence criteria. As it is possible to vary by hands the values of the convergence thresholds, it is a good idea to try a same calculation with different values of the relevant thresholds and to see how the frequencies depend on them.  I have a comment, though. As Nicola said in a previous post, it is unlikely (though not impossible) that phonons go soft because the structure is poorely optimized. It is worth stopping a while and understand carefully why it is both unlikely _and_ possible. It is possible because squared phonon frequencies are related to the Hessian (matrix of second derivatives) of the energy with respect to the nuclear coordinates. It is perfectly possible that, out of equilibrium, one or more of the eigenvalues of this matrix are negative. I even believe (but would like not to comment on this statement) that the ratio of the configuration space corresponding to all positive eigenvalues to the total configuration space decreases exponentially with the number of atoms. HOWEVER, phonon frequencies only depend on second derivatives of the energy, and not directly on first derivatives. So, if you are not too far from the minimum, not too good a convergence (i.e. the existence of nonvanishing forces) may certainly affect the accuracy of phonon calculations (it would not if the crystal were perfectly harmonic), but not the stability of the system (i.e. the sign of any of the squared frequencies).</DIV><DIV><BR><BLOCKQUOTE type="cite"><DIV>    2, I knew from most refs that ecutwfs  30 Ry for ultrasoft-PP is enough in pw calculations , Is it right ? <BR></DIV></BLOCKQUOTE><DIV><BR class="khtml-block-placeholder"></DIV>This statement is meaningless. 30 Ry may be too large or too small a cutoff, depending on the PP that you are using. It is true that many US PP calculations are done with cutoffs in this range, but whether or not this value is appropriate in your case ONLY YOU CAN TELL. i.e. you have to perform careful convergence tests in order to be sure that your PW basis set is adequate to describe the properties that you want with the accuracy that you need (the same holds, of course, for many other devil's details of the calculation: k-points, gaussian smearing, etc.)</DIV><DIV><BR><BLOCKQUOTE type="cite"><DIV>But i donot have much idear about ecutwfs in ph calculations .  <BR></DIV></BLOCKQUOTE><DIV><BR class="khtml-block-placeholder"></DIV>By and large, (squared) phonon frequencies are second derivatives of the crystal energy for a given kinetic-energy cutoff. If a cutoff is OK for geometry optimization, it  is also so for phonon calculation (see below)</DIV><DIV> <BR><BLOCKQUOTE type="cite"><DIV>Would you please give me some refs of more detail settings about PH calculations ? </DIV></BLOCKQUOTE><DIV><BR class="khtml-block-placeholder"></DIV>S. Baroni, A. Dal Corso, S. de Gironcoli, and P. Giannozzi, Rev. Mod. Phys. 73, 515 (2001).</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>Enjoy!</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>Stefano B.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>PS: for us, it is considered polite to sign with our first name, or with the first name followed by the family name. It would be unpolite (except for formal business correspondence) to prepend the title to the signature - "Stefano" or "Stefano B.": confidential, OK; "Stefano Baroni" a bit less confidential, also OK in most cases; "Prof. Stefano Baroni": very formal and kind of assuming and not-so-polite in many cases - OK in a few</DIV><DIV><BR><BLOCKQUOTE type="cite"> <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 <BR>than GGA calculations (the exchange-correlation terms, calculated in<BR>real space, introduce very high frequencies in reciprocal space).<BR><BR>                       nicola<BR><BR><BR>>> Dear All,<BR>>> I also have calculated gamma phonon frequencies of anatase TiO2 . <BR>>> I run ph.x caculation after relaxing the structure, but  9 of 18 frequencies<BR>>> are minus.<BR>>><BR>>> Would you please give me some hints and help?<BR>>><BR>>> 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></BLOCKQUOTE></DIV><BR><DIV> <SPAN class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><SPAN class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><SPAN class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">---</SPAN></SPAN></SPAN></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">Stefano Baroni - SISSA</SPAN></SPAN></SPAN></FONT><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">  </SPAN></SPAN></SPAN></FONT><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">&</SPAN></SPAN></SPAN></FONT><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">  </SPAN></SPAN></SPAN></FONT><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">DEMOCRITOS National Simulation Center - Trieste</SPAN></SPAN></SPAN></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">[+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)</SPAN></SPAN></SPAN></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; ; font-size: 12px; "><BR style="; font-size: 12px; "></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">Please, if possible, don't</SPAN></SPAN></SPAN></FONT><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">  </SPAN></SPAN></SPAN></FONT><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">send me MS Word or PowerPoint attachments</SPAN></SPAN></SPAN></FONT></DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">Why? See:</SPAN></SPAN></SPAN></FONT><FONT class="Apple-style-span" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;"><SPAN class="Apple-style-span" style="font-size: 12px; "><SPAN class="Apple-style-span" style="font-size: 12px; ">  </SPAN></SPAN></SPAN></FONT><A href="http://www.gnu.org/philosophy/no-word-attachments.html"><FONT class="Apple-style-span" color="#0023E9" size="3"><SPAN class="Apple-style-span" style="font-size: 12px;; color: rgb(0, 35, 233); -khtml-text-decorations-in-effect: underline; "><SPAN class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; "><SPAN class="Apple-style-span" style="color: rgb(0, 35, 233); font-size: 12px; -khtml-text-decorations-in-effect: underline; ">http://www.gnu.org/philosophy/no-word-attachments.html</SPAN></SPAN></SPAN></FONT></A></DIV><DIV><BR class="khtml-block-placeholder"></DIV><BR class="Apple-interchange-newline"></SPAN></SPAN></SPAN> </DIV><BR></BODY></HTML>