Dear Prof. Marks,<br><br>I&#39;ve sent the whole output file to your private e-mail as you suggested.<br>Here I just put the result of grep command:<br><br>CVS version:<br>!ššš total energyššššššššššššš =ššš -242.65271405 Ry<br>
šššš number of bfgs stepsššš =šš 0<br>!ššš total energyššššššššššššš =ššš -242.66331605 Ry<br>šššš number of bfgs stepsššš =šš 1<br>!ššš total energyššššššššššššš =ššš -242.66341971 Ry<br>šššš number of bfgs stepsššš =šš 2<br>
!ššš total energyššššššššššššš =ššš -242.66342267 Ry<br>šššš bfgs converged inšš 4 scf cycles andšš 3 bfgs steps<br>!ššš total energyššššššššššššš =ššš -242.66284562 Ry<br><br>v.4.2.1<br>!ššš total energyššššššššššššš =ššš -242.65271405 Ry<br>
šššš number of bfgs stepsššš =šš 0<br>!ššš total energyššššššššššššš =ššš -242.64460983 Ry<br>šššš number of bfgs stepsššš =šš 1<br>!ššš total energyššššššššššššš =ššš -242.66331847 Ry<br>šššš number of bfgs stepsššš =šš 1<br>
!ššš total energyššššššššššššš =ššš -242.66341565 Ry<br>šššš number of bfgs stepsššš =šš 2<br>!ššš total energyššššššššššššš =ššš -242.66342170 Ry<br>šššš number of bfgs stepsššš =šš 3<br>!ššš total energyššššššššššššš =ššš -242.66342336 Ry<br>
šššš bfgs converged inšš 6 scf cycles andšš 4 bfgs steps<br><br>Best regards,<br>Maxim.<br><br><br><div class="gmail_quote">2011/4/18 Laurence Marks <span dir="ltr">&lt;<a href="mailto:L-marks@northwestern.edu">L-marks@northwestern.edu</a>&gt;</span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">The latest BFGS has a curvature trap which can change the behavior<br>
particularly if you have too few k-points. Can you please do<br>
<br>
grep -e &quot;! š štot&quot; -e bfgs WHATEVER.out<br>
<br>
(where WHATEVER is relevant) and send it to me, perhaps at my private<br>
email (or the full .out file). I am not an expert with QE, but might<br>
be able to say something.<br>
<br>
2011/4/18 νΑΛΣΙΝ πΟΠΟΧ &lt;<a href="mailto:max.n.popov@gmail.com">max.n.popov@gmail.com</a>&gt;:<br>
<div><div></div><div class="h5">&gt; Dear QE developers,<br>
&gt;<br>
&gt; comparing vc-relax calculations between BFGS(4.2.1) and the latest BFGS(from<br>
&gt; CVS) I facedš with a following problem:<br>
&gt; the newest QE calculates one extra point after achieving convergence<br>
&gt; criteria (etot_conv_thr,force_conv_thr,press_conv_thr),<br>
&gt; which leads to slightly higher total energy.<br>
&gt;<br>
&gt; Here are the two latest energies from outputs:<br>
&gt; šššššš v.4.2.1šššššššš š š š š š ššš CVS(18.04.2011)<br>
&gt;<br>
&gt; -242.66342170 Ryšššššš -242.66342267 Ry<br>
&gt; -242.66342336 Ryšššššš -242.66284562 Ry<br>
&gt;<br>
&gt; One can see that there is a good agreement between the last energy by<br>
&gt; v.4.2.1 and next to the last by CVS-version, whereas<br>
&gt; the last energy by CVS is higher by ~ 0.5 mRy.<br>
&gt;<br>
&gt; I tried relax (cell shape fixed) also - it works fine, no extra point<br>
&gt; calculation.<br>
&gt;<br>
&gt; --<br>
&gt; Best regards, Max Popov<br>
&gt; Ph.D. student<br>
&gt; Materials center Leoben (MCL), Leoben, Austria.<br>
&gt;<br>
</div></div>&gt; _______________________________________________<br>
&gt; Pw_forum mailing list<br>
&gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
&gt;<br>
&gt;<br>
<br>
<br>
<br>
--<br>
Laurence Marks<br>
Department of Materials Science and Engineering<br>
MSE Rm 2036 Cook Hall<br>
2220 N Campus Drive<br>
Northwestern University<br>
Evanston, IL 60208, USA<br>
Tel: (847) 491-3996 Fax: (847) 491-7820<br>
email: L-marks at northwestern dot edu<br>
Web: <a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a><br>
Chair, Commission on Electron Crystallography of IUCR<br>
<a href="http://www.numis.northwestern.edu/" target="_blank">www.numis.northwestern.edu/</a><br>
Research is to see what everybody else has seen, and to think what<br>
nobody else has thought<br>
Albert Szent-Gyorgi<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" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div>