<font color='black' size='2' face='arial'><font size="2" color="black" face="arial">

<div> Dear Yuyang, <br>
<br>
If I may ask (I did not follow the discussion, only now I saw your e-mail, so I don't know what was discussed previously), what are the values of the "conv_thr" and "tr2_ph" parameters used in your case? I never did phonon calculations for your specific case, but using USPP for Ag/Cu/Au/Pt, the phonon frequencies are already converged at wfc and charge cutoffs of 40/480, as long as one performs a very accurate scf and phonon calculation. For silver, for example, I had to use a "conv_thr" as low as 10-12, and tr2_ph about 10-19.<br>
<br>
Fabio Negreiros Ribeiro<br>
Post-doc at CNR<br>
Pisa - Italy</div>

<div> <br>

</div>



<div style="font-family:arial,helvetica;font-size:10pt;color:black">-----Original Message-----<br>

From: YY &lt;xrhino@gmail.com&gt;<br>

To: Paolo Giannozzi &lt;giannozz@democritos.it&gt;<br>

Cc: pw_forum &lt;pw_forum@pwscf.org&gt;<br>

Sent: Wed, Apr 18, 2012 5:26 pm<br>

Subject: Re: [Pw_forum] ecutrho in phonon calculations<br>

<br>









<div id="AOLMsgPart_0_26c81fa7-6029-40b9-b55c-5fc81da07bad" style="margin: 0px;font-family: Tahoma, Verdana, Arial, Sans-Serif;font-size: 12px;color: #000;background-color: #fff;">

<pre style="font-size: 9pt;"><tt>Hi Paolo,

Thanks for your reply.  For the phonon calculation, I use a PBE-US PP
(Fe.pbe-sp-van_ak.UPF, and Se.pbe-van.UPF).
According to your suggestion, I increase the ecutwfc, and do several
tests.  Unfortunately, it seems there are still convergence problem:

ecutwfc  ecutrho  f1(cm-1) ... f12(cm-1)
100        1200      106.5    ...  329.5
100        2400      127.2    ...  339.8
200        2400      127.8    ...  341.2
200        4800        32.9    ...  319.8

Could there be any other reasons to cause this convergence problem?

Yuyang@RPI

On Mon, Apr 16, 2012 at 11:52 AM, Paolo Giannozzi
&lt;<a __removedlink__527898086__href="mailto:giannozz@democritos.it">giannozz@democritos.it</a>&gt; wrote:
&gt; On Thu, 2012-04-12 at 16:07 -0400, YY wrote:
&gt;
&gt;&gt; But I find there is no simple convergence trend while increasing
&gt;&gt; ecutrho. &nbsp;When I use 35 Ry ecutwfc, and fix all the other parameters,
&gt;&gt; the results are like this:
&gt;&gt;
&gt;&gt; &nbsp; &nbsp;ecutrho &nbsp; &nbsp; &nbsp; f1(cm-1) &nbsp; &nbsp; &nbsp;.... &nbsp;f12(cm-1)
&gt;&gt; &nbsp; &nbsp;420 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;40.475 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;324.062
&gt;&gt; &nbsp; &nbsp;600 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 123.746 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 350.727
&gt;&gt; &nbsp; &nbsp;900 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; -251.876 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;306.721
&gt;
&gt; unfortunately in at least some [ and I am afraid, in many ] cases,
&gt; a rather high cutoff is needed anyway to get decent phonons. If
&gt; this happens, the advantage of using USPP becomes dubious, given
&gt; their significant computational overhead wrt Norm-Conserving PPs.
&gt; Which USPP did you use ?
&gt;
&gt; P.
&gt;
&gt; --
&gt; Paolo Giannozzi, IOM-Democritos and University of Udine, Italy
&gt;
&gt;
_______________________________________________
Pw_forum mailing list
<a __removedlink__527898086__href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a __removedlink__527898086__href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>

</tt></pre>
</div>

 <!-- end of AOLMsgPart_0_26c81fa7-6029-40b9-b55c-5fc81da07bad -->



</div>

</font></font>