Sometimes it is also difficult to converge the total energy. So usually it is better to check the convergence of the quantity you are interested in.<br /><br />But anyway it is really strange that total energy doesn&#39;t converge at such abnormally high cut off values.<br /><br />Did you test the pseudopotentials individually?? For eg. did you try calculating the bond lengths of say a nitrogen molecule & see how the results are??<br /><br />With regards,<br /><br />Prasenjit.<br /><br />On Nov 6, 2008 5:45pm, Timothy Mason &lt;thmason24@yahoo.com&gt; wrote:<br />&gt; Dear Lorenzo,<br />&gt; <br />&gt; here are results for the lower energy range.  Still doesn&#39;t seem to be converging.<br />&gt; <br />&gt; ecut     energy              energy differnce<br />&gt; 30        -48.19258627    <br />&gt; 35        -48.19503183       0.00244556<br />&gt; 40        -48.197306           0.00227417<br />&gt; 45        -48.20217820       0.0048722<br />&gt; <br />&gt; thanks<br />&gt; <br />&gt; --- On Thu, 11/6/08, Lorenzo Paulatto paulatto@sissa.it&gt; wrote:<br />&gt; From: Lorenzo Paulatto paulatto@sissa.it&gt;<br />&gt; Subject: Re: [Pw_forum] Total energy does not coverge with ecut<br />&gt; To: thmason24@yahoo.com, &quot;PWSCF Forum&quot; pw_forum@pwscf.org&gt;<br />&gt; Date: Thursday, November 6, 2008, 10:08 AM<br />&gt; <br />&gt; Timothy Mason ha scritto:<br />&gt; &gt;<br />&gt; &gt;     Gentlemen,<br />&gt; &gt;<br />&gt; &gt; My structures do not see to be converging with increased ecutwfc.  <br />&gt; &gt; I&#39;m setting ecutrho to ecutwfc * 12 and using US PP.    below is a<br />&gt; &gt; table of ecut values vs energies vs difference from previous steps<br />&gt; &gt; followed by my input file followed by part of the output file.    It<br />&gt; &gt; appears to begin to converge and 120 ecut but then jumps at 130 and<br />&gt; &gt; then more at 140.<br />&gt; &gt;<br />&gt; Dear Timothy,<br />&gt; these are ridiculous cutoff energy for ultrasoft pseudopotentials - they<br />&gt; are very high for norm-conserving too! I suspect the presence of a<br />&gt; high-energy ghost-state in one of the pseudopotential; these kind<br />&gt;  of<br />&gt; states aren&#39;t usually harmful but can suddenly pop-up under some<br />&gt; specifics conditions like high pressure or very high cutoff.<br />&gt; <br />&gt; What is happening when you try some reasonable cutoff in the 20-50 Ry range?<br />&gt; <br />&gt; regards<br />&gt; <br />&gt; -- <br />&gt; Lorenzo Paulatto<br />&gt; SISSA  &  DEMOCRITOS (Trieste)<br />&gt; phone: +39 040 3787 511<br />&gt; skype: paulatz<br />&gt; www:   http://people.sissa.it/~paulatto/<br />&gt; <br />&gt; Due to budget cuts approved by Italian Minister of Finance and Italian<br />&gt; Minister &quot;for&quot; Education, University and Research this may be the<br />&gt; last<br />&gt; email you receive from me.<br />&gt; <br />&gt; Please send any future correspondence, possibly via pigeon, to<br />&gt; mr. Lorenzo Paulatto<br />&gt; under Roiano Railway Bridge - Pillar 2<br />&gt; 34136, Trieste<br />&gt; <br />&gt; Donations in food, coins, booze and processor time are warmly welcome.<br />&gt; <br />&gt; <br />&gt; <br />&gt;