<div dir="ltr">Dear Sonu Kumar<br>Thank you for your guidelines. I will test them.<br>Best regards,Payam Norouzzadeh<br><br><div class="gmail_quote">On Thu, Mar 8, 2012 at 7:11 AM, Sonu Kumar <span dir="ltr"><<a href="mailto:1009ukumar@gmail.com">1009ukumar@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote"><div class="im">On Thu, Mar 8, 2012 at 4:11 AM, Payam Norouzzadeh <span dir="ltr"><<a href="mailto:payam.norouzzadeh@gmail.com" target="_blank">payam.norouzzadeh@gmail.com</a>></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 dir="ltr">Dear QE users<br><br>My system has 54 atoms with a BCC structure. Its primitive unit cell includes 27 atoms. <br>I am trying to find its lattice constant by using vc-relax.<br>After running it several times and for more than 3 weeks still the convergence thresholds are not satisfied. <br>



The forces don't like to get smaller than a specific value (around 0.02).<br>Is something wrong with my input file? Please let me know your comments.<br><br>Best regards,Payam Norouzzadeh<br><br><br>&CONTROL<br>  title = Eight,<br>



  calculation  = "vc-relax",<br>  restart_mode = 'from_scratch' ,<br>  prefix       = "Eight",<br>  pseudo_dir   = "/...",<br>  wfcdir = '...' ,<br>  outdir = '/...',<br>



  nstep =  100  ,<br>  etot_conv_thr = 1.0E-9  ,<br>  forc_conv_thr = 5.0D-5 ,<br>  max_seconds = 172800 ,<br>  dt = 100 ,<br>/<br></div></blockquote></div><div><br>force and energy threshold seems to very high, especially <br>


the latter. reduce them.<br> </div><div class="im"><blockquote class="gmail_quote" style="border-left:1px solid rgb(204,204,204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex"><div dir="ltr">&SYSTEM<br>  ibrav = 3, <br>  celldm(1)=19.622,<br>


  nat  = 27,<br>  ntyp = 3,<br>  ecutwfc   = 70.D0,<br>
  ecutrho   = 700.D0,<br>  occupations='smearing', <br>  smearing='gaussian',<br>  degauss=0.02,<br>/<br>&ELECTRONS<br>diagonalization = 'cg',<br>  conv_thr    = 1.D-9,<br>  mixing_beta = 0.5D0,<br>



/<br></div></blockquote></div><div><br>use diagonalization= david<br> </div><div class="im"><blockquote class="gmail_quote" style="border-left:1px solid rgb(204,204,204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex"><div dir="ltr">
&IONS<br> ion_dynamics = 'damp' ,<br>

/<br></div></blockquote></div><div><br>use default bfgs algorithm for dynamics <br> 
<br></div><div class="im"><blockquote class="gmail_quote" style="border-left:1px solid rgb(204,204,204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex"><div dir="ltr">&CELL<br>cell_dynamics = 'damp-w',<br>/<br></div>
</blockquote>

</div><div><br>similarly, use default bfgs. <br></div><blockquote class="gmail_quote" style="border-left:1px solid rgb(204,204,204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<br></blockquote></div><br>GOOD LUCK !!<br><br>SKumar<br>IITD<br>
</blockquote></div><br></div>