<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt">Hi 2 All<br>I want the chemical potential of oxygen molecule , but the binding energy that I get is too far from&nbsp; experiment and theoretical results. Could anyone help me How can I improve my calculation.<br>&amp;CONTROL<br>&nbsp; calculation&nbsp; = "relax",<br>&nbsp; restart_mode='from_scratch'<br>&nbsp; prefix&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "o2",<br>&nbsp; pseudo_dir&nbsp;&nbsp; = "./",<br>&nbsp; outdir&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "/p5/batch/kazempou/",<br>&nbsp; verbosity='high'<br>/<br>&amp;SYSTEM<br>&nbsp; ibrav&nbsp;&nbsp;&nbsp;&nbsp; = 1,<br>&nbsp; celldm(1)=20<br>&nbsp; nat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2,<br>&nbsp; nosym=.true.<br>&nbsp; ntyp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1,<br>&nbsp; ecutwfc&nbsp;&nbsp; = 45.D0,<br>&nbsp; ecutrho&nbsp;&nbsp; =
 400.D0,<br>&nbsp;occupations='smearing', smearing='mp', degauss=<br>0.001,nspin=2,starting_magnetization(1)= 1.0,<br><br>/<br>&amp;ELECTRONS<br>&nbsp; conv_thr&nbsp;&nbsp;&nbsp; = 1.D-5,<br>&nbsp; mixing_beta = 0.7D0,<br>/<br>&amp;IONS<br>/<br>ATOMIC_SPECIES<br>O&nbsp; 1.00 o.optgga1.fhi.UPF<br>ATOMIC_POSITIONS crystal<br>O 0.5&nbsp;&nbsp; 0.5&nbsp; 0.5 0 0 0<br>O&nbsp; 0.5&nbsp; 0.5&nbsp; 0.0<br>K_POINTS {Gamma}<br><br>thanks a lot<br><div>&nbsp;</div>Ali Kazempour<br><br><br><br><div><br>Fritz-Haber-Institut              fax   : ++49-30-8413 4701<br>der Max-Planck-Gesellschaft<br>Faradayweg 4-6                    e-mail: kazempou@fhi-berlin.mpg.de<br>D-14 195 Berlin-Dahlem / German<div><br></div></div>
<!-- cg14.c4.mail.gq1.yahoo.com compressed/chunked Mon Nov  9 05:15:07 PST 2009 -->
</div><br>

      </body></html>