Thanks a lot. <div><span class="gmail_quote">On 4/5/07, <b class="gmail_sendername"><a href="mailto:cesards@msi.umn.edu">cesards@msi.umn.edu</a> <<a href="mailto:cesards@msi.umn.edu">cesards@msi.umn.edu</a>> wrote:

</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Javier,<br><br>After your first cell optimization, try starting again from scratch

<br>using the final configuration of the previous calculation. If the problem<br>persists, increase the cutoff energies.<br><br>Good luck.<br><br>Cesar R.S. da Silva<br><br>><br>> Message: 10<br>> Date: Wed, 4 Apr 2007 17:31:23 +0200

> From: "Javier Antonio Montoya" <<a href="mailto:j.antonio.montoya@gmail.com">j.antonio.montoya@gmail.com</a>> > To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a> > Subject: [Pw_forum] vc-relax bug?

<br>> Reply-To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>><br>> ------=_Part_26133_24163750.1175700683893<br>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>> Content-Transfer-Encoding: 7bit

<br>> Content-Disposition: inline<br>><br>> Hello dear all,<br>><br>> In recent days I have been trying to get a triclinic cell to have a<br>> diagonal<br>> stress by using vc-relax. The pw.x program reads data in the correct way,

<br>> calculates the right starting volume for the cell given as input and also<br>> seems to do right the first electronic minimization but, as soon as the<br>> Wentzcovitch Damped Cell-Dynamics Minimization starts, it gives cell

<br>> parameters and ion positions that are not consistent with the energies and<br>> stress that it prints (i.e. if you take those cells and ions to an<br>> independent scf run the final energy and stress differ from the ones in

<br>> the<br>> vc-output) the problem is that in the end, when I get a nice diagonal<br>> stress<br>> after 60 or more dynamic steps, the difference between the reported energy<br>> in the vc-output and the energy of the system that you get with the cell

> and > ion parameters shown is of almost 1 Ry and the pressures differ by more > than > 10GPa. Any idea? by the way the ph.x complain and doesn't get anything > from > such restart files.

<br>> The output files are huge for sending by e-mail because my system is 72<br>> atoms big, but I can provide them.<br>><br>> cheers,<br>><br>> JAvier.<br>><br><br><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>