<p>Dear Ayedun,</p>
<p>Try using a higher ecutwfc, like 50 for example. Also a smaller degauss value. Finally set the ecutrho to 4*ecutwfc if your pseudopotential is not US, and to 12*ecutwfc if it is.</p>
<p>It`s also a good idea to specify the number of bands, to make sure 
you got enough ones on the calculation. If unsure you can just use a 
high value and reduce it after a test calculation to compensate for the 
extra empty bands.</p>
<p>Choose a single ibrav value (if unsure, go for the experimental value
 or close to it) and play with these parameters, checking how the total 
energy changes with them. Once you get convergence on the total energy, then you'll have a 
reliable setup that will allow you to change the ibrav values and obtain
 good curves.<br>
</p>
<p>Sincerely,</p>
<p>Giovani Faccin</p>
<p>UFMS/Brazil<br>
</p>

<br><br><div class="gmail_quote">2012/4/8 funmi ayedun <span dir="ltr"><<a href="mailto:funmiayedun@yahoo.com">funmiayedun@yahoo.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font:inherit" valign="top">Dear All,<br><br>Good day all, I am a new user in quantum espresso, I am working on anatase TiO2 (with hexagonal structure) but  I have problem with convergence test  for lattice parameter for TiO2 system.  I have tried to vary celldm 1 from 8.2 to 9.0Bohr but I was unable get a convergence result. Also, I tried to vary celldm 3  from 0.2441 to 1.3441Bohr and  keeping celldm 1 constant, I was unable to get convergence as well. I have my input file below. <br>
<br>&control<br>  calculation = 'scf',<br>  restart_mode='from_scratch',<br>  prefix='TiO2'<br>  pseudo_dir='/home/funmilayo/espresso-4.2.1/pseudo',<br>  outdir='./',<br>    tprnfor = .true.<br>
  tstress = .true.<br>/<br>&system<br>  ibrav=6, celldm(1)= 8.6808, celldm(3)= 0.7450,<br>  nat=3,
 ntyp=2,occupations='smearing', smearing='marzari-vanderbilt', degauss=0.02 <br>  ecutwfc= 25.0,<br>/<br>&electrons<br> mixing_beta = 0.7<br> diagonalization='davidson'<br>conv_thr =  1.0d-8<br>
/<br>ATOMIC_SPECIES<br>Ti 47.867  Ti.blyp-sp-van_ak.UPF<br>O  15.994  O.blyp-van_ak.UPF<br>ATOMIC_POSITIONS<br>Ti 0.00000000   0.00000000   0.00000000 <br>O  0.33333333   0.33333333   0.33333333<br>O  0.80000000   0.20000000   0.50000000<br>
K_POINTS (automatic)<br>6 6 6  0  0  0<br><br>Can anyone please tell me where I have gone wrong.  <br>Thank you all.<br><br><br>Ayedun Oluwafunmilayo,<br>PhD student,<br>Department of Physics,<br>University of Agriculture,<br>
Abeokuta,<br>Nigeria.<br></td></tr></tbody></table><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>