Try cell_dynamics=&#39;bfgs&#39; and not &#39;damp&#39;<br>Success<br><br><div class="gmail_quote">On Tue, Jul 27, 2010 at 12:42 PM, Olga Sedelnikova <span dir="ltr">&lt;<a href="mailto:o.sedelnikova@gmail.com">o.sedelnikova@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><font color="#330099"><font size="2"><font face="comic sans ms,sans-serif">Dear All,</font></font></font></div>

<div><font color="#330099"><font size="2"><font face="comic sans ms,sans-serif">I have a problem with &#39;vc-relax&#39; option for graphene. </font></font></font><font color="#330099"><font size="2"><font face="comic sans ms,sans-serif">I have used david and cg diagonalization proceedures, different types of smearing and ion and cell dynamics but all calculations were crashed. </font></font></font></div>



<div><font color="#330099"><font size="2"><font face="comic sans ms,sans-serif">One of my input file:</font></font></font></div>
<div><font color="#330099" face="Comic Sans MS"></font> </div>
<div><font color="#330099" face="Comic Sans MS">&amp;CONTROL<br>     calculation = &#39;vc-relax&#39;,<br>     prefix = &#39;graphene&#39;,<br>     restart_mode=&#39;from_scratch&#39;<br>     tstress=.true.,<br>     tprnfor=.true.,<br>


     nstep=45,<br>      etot_conv_thr = 1.0E-4 ,<br>      forc_conv_thr = 1.0D-3 ,<br>      max_seconds = 36000 ,<br>      dt = 100,<br>   /<br>&amp;SYSTEM<br>     ibrav = 4,<br>     celldm(1)= 4.64833<br>     celldm(3)=5.447621<br>


               nat = 2, <br>               ntyp = 1, <br>               ecutwfc = 18.0,<br>               nbnd =8, <br>     occupations=&#39;smearing&#39;<br>     smearing = &#39;mv&#39;<br>     degauss = 0.005<br>       /<br>


&amp;ELECTRONS<br>     conv_thr = 1.0d-7<br>     electron_maxstep = 70<br>     diagonalization=&#39;cg&#39;<br>   /<br>&amp;IONS<br>     ion_dynamics=&#39;damp&#39;<br>/<br>&amp;CELL<br>    cell_dynamics = &#39;damp-w&#39; ,<br>


    press=0.00 ,<br>    wmass=0.0015 ,<br>/<br>  <br>ATOMIC_SPECIES<br>     C   12.01100 C.pz-vbc.UPF<br>ATOMIC_POSITIONS crystal <br>    C      0.00000000    0.00000000    0.50000000<br>    C      0.33333333    0.66666666    0.50000000<br>


K_POINTS {automatic}<br>   6 6 6 0 0 0<br></font></div>
<div><font color="#330099" face="Comic Sans MS">Inthis cause: </font></div>
<div><font color="#330099" face="Comic Sans MS">     </font></div>
<div><font color="#330099" face="Comic Sans MS">from electrons : error #         1<br>     charge is wrong<br></font></div>
<div><font color="#330099"><font size="2"><font face="comic sans ms,sans-serif">Maybe somebody has already dealed with this problem and can advice the correct options for graphene.</font></font></font></div>
<div><font color="#330099"><font size="2"><font face="comic sans ms,sans-serif">Any suggestion will be appreciated.</font></font></font></div>
<div><br clear="all">-- <br>Best wishes,<br>Olga Sedelnikova<br>Nikolaev Institute of Inorganic Chemistry of SB RAS<br></div>
<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>
<br></blockquote></div><br><br clear="all"><br>-- <br>Phillip W. Otieno Nyawere,<br>Kabarak University,<br>Dept of Physics &amp; Mathematics,<br>P.O Box Private Bag - 20157,<br>Kabarak, Kenya.<br>Tel +254728342054<br><a href="mailto:pnyawere@gmail.com" target="_blank">pnyawere@gmail.com</a>, <a href="mailto:potieno@kabarak.ac.ke" target="_blank">potieno@kabarak.ac.ke</a><br>
<br>The battle belongs to the Lord.<br>CURRENTLY IN CAPE TOWN, SOUTH AFRICA<br><br>