<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>

Dear Xijun,<div><br></div><div>Thanks very much for your advice. I will do the necessary changes and let you know how things go.</div><div><br></div><div>Regards</div><div><br></div><div>Elie<br><br><div><div id="ecxSkyDrivePlaceholder"></div>&gt; Date: Wed, 7 Sep 2011 15:58:37 -0400<br>&gt; From: xijunw@gmail.com<br>&gt; To: pw_forum@pwscf.org<br>&gt; Subject: Re: [Pw_forum] problems with vc-relax<br>&gt; <br>&gt; On Wed, Sep 7, 2011 at 3:28 PM, Elie Moujaes &lt;elie.moujaes@hotmail.co.uk&gt; wrote:<br>&gt; &gt; Dear all,<br>&gt; &gt;<br>&gt; <br>&gt; Dear Elie,<br>&gt; <br>&gt; &gt;<br>&gt; &gt; I am relaxing a 60 atom supercell (graphene grain boundary) to make the<br>&gt; &gt; total force on atoms mainly 0.001 Ry/au. However the vc-relax stopped after<br>&gt; &gt; 4 days of execution and without reaching the desired force threshold. The<br>&gt; &gt; following output was obtained :<br>&gt; &gt;<br>&gt; &gt; Please can anyone advice me on this? Shall I start with the new<br>&gt; &gt; configuration and start relaxing again?<br>&gt; &gt;<br>&gt; <br>&gt; I'm also doing some cell relaxation. As I understand, cell<br>&gt; optimization is not very easy, and it is common that it does not<br>&gt; converge after a very long run.<br>&gt; <br>&gt; Surely you should restart the calculation to achieve your convergence<br>&gt; criteria. You can just change the "restart_mod = 'from_scratch' " to<br>&gt; "restart_mode = 'restart' ", and submit the job again to continue from<br>&gt; where you stopped. Please be advised that pw.x will read cell<br>&gt; parameters / atomic positions, initial density, pseudopotentials from<br>&gt; the folder of GBphonon.save, so make sure you copy the folder to<br>&gt; computing-nodes.<br>&gt; <br>&gt; However, I highly suggest you do something to accelerate the<br>&gt; convergence. For example, if you do not expect the cell angle to<br>&gt; change (i.e. keep angles=90), you may add the following keywords to<br>&gt; the &amp;CELL section.<br>&gt; <br>&gt;    cell_dofree = xyz<br>&gt; <br>&gt; This will greatly reduce the uncertainty of the geometry optimization.<br>&gt; and reduce energy/force oscillation.<br>&gt; <br>&gt; Another great idea is to manually choose x, y and z, then do a<br>&gt; relaxation of atomic positions only. and plot the total energy changes<br>&gt; with x, y and z. However, be noted that different volume means<br>&gt; different numbers of PW basis function, and energy comparison in this<br>&gt; way could be a little confusing. Anyway, if volume change is not<br>&gt; large, still this would help a lot.<br>&gt; <br>&gt; Since I'm also a learner, any more comments or correction to my<br>&gt; comments are welcome. Thanks<br>&gt; <br>&gt; Xijun<br>&gt; <br>&gt; &gt;<br>&gt; &gt; Regards<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; Elie Moujaes<br>&gt; &gt;<br>&gt; &gt; University of Nott<br>&gt; &gt;<br>&gt; &gt; University Park<br>&gt; &gt;<br>&gt; &gt; NGT 3RD<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; _______________________________________________<br>&gt; &gt; Pw_forum mailing list<br>&gt; &gt; Pw_forum@pwscf.org<br>&gt; &gt; http://www.democritos.it/mailman/listinfo/pw_forum<br>&gt; &gt;<br>&gt; &gt;<br>&gt; <br>&gt; <br>&gt; <br>&gt; -- <br>&gt; Department of Chemistry and Biochemistry, Concordia University<br>&gt; 7141 Sherbrooke Ouest, Montreal, QC H4B 1R6<br>&gt; Tel: 514-848-2424-#5835 (Lab) TA office: SP175.23<br>&gt; _______________________________________________<br>&gt; Pw_forum mailing list<br>&gt; Pw_forum@pwscf.org<br>&gt; http://www.democritos.it/mailman/listinfo/pw_forum<br></div></div>                                               </div></body>
</html>