Dear Cesar,<br>
<br>
Thank you for your reply. According to your suggestion, I change "dt"
and conv_thr to 80 and 1.D-10, the convergence is not improved,
although the calculation is running.  The  calculated 
total  force  and  stress  is  given here as a
attachment because it is too long. As you can see from this file, the
final cell parameters is changed  largely, which also affects the
k-points meshing during my calculation. Another point added here is
that my systems only contain some light elements like H, Li. How to
further optimize these parameters?<br>
<br>
Regards,<br>
C.H.<br><br><div><span class="gmail_quote">On 10/23/07, <b class="gmail_sendername">Cesar R.S. da Silva</b> <<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;">
Dear Chao,<br><br>Your dt setting should work OK. wmass could be reduced to somthing like<br>0.002, provided atomic masses are all set to the average mass (amu), but<br>your value is not that bad. The only thing I think could produce
<br>convergence problems is your conv_thr. Try conv_thr = 1e-9. Variable cell<br>requires precise stress calculation, and a more stringent convergence<br>criteria for the SCF cycle is required.<br><br>Cesar.<br><br>
</blockquote></div><br>-- <br>============================<br>C.H. Hu<br>Postdoctoral fellow<br>Chimie et Physico-Chimie appliquees<br>Institut Francais du Petrole (IFP)<br>Rueil-Malmaison, France<br>E_mail: <a href="mailto:chaohao.mse@gmail.com">
chaohao.mse@gmail.com</a><br>============================