<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
How did you compute EO? You need to use spin polarization and to
break symmetry in order to get the lowest possible energy (i.e.
Hund's rule).<br>
<br>
kind regards<br>
<br>
<div class="moz-cite-prefix">On 04/28/2014 10:06 AM, Junxiang Chen
wrote:<br>
</div>
<blockquote cite="mid:007901cf62b8$c3afda20$4b0f8e60$@gmail.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:\5B8B\4F53;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@\5B8B\4F53";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
text-align:justify;
text-justify:inter-ideograph;
font-size:10.5pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
/* Page Definitions */
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Dear all,</span><span
lang="EN-US"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The same problem as <a
moz-do-not-send="true"
href="http://qe-forge.org/pipermail/pw_forum/2009-November/089525.html">http://qe-forge.org/pipermail/pw_forum/2009-November/089525.html</a>
, the binding energy of detached O2 is 6.57eV, which is far
from experiment value ( 5.12eV ), or 5.7eV in this page.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Input file ( O2 ):<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&SYSTEM<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ibrav = 12,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
celldm(1) = 16.0348541363689,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
celldm(2) = 1.52752767,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
celldm(3) = 2.665651575,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
celldm(4) = -0.3273333438252,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
nat = 2,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ntyp = 1,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ecutwfc = 45,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ecutrho = 400,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
occupations = 'smearing' ,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
degauss = 0.01,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
smearing = 'methfessel-paxton' ,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
nspin = 2,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
starting_magnetization(1) = 1,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&ELECTRONS<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
electron_maxstep = 200,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
conv_thr = 1.D-6,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
mixing_beta = 0.4,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
mixing_mode = 'local-TF' ,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">&IONS<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">
ion_dynamics = 'bfgs' ,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">/<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ATOMIC_SPECIES<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> O 15.9994
O.pbe-rrkjus.UPF<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">ATOMIC_POSITIONS crystal<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O 0.217969593
0.389882242 0.401960243<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">O 0.217971037
0.389881524 0.456565745<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">K_POINTS automatic<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"> 5 5 2 0 0 0<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The binding energy is
calculated by BE=EO2-2*EO* and BE=6.57eV <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Could anyone tells me
what’s wrong? <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">-------------------------------------------------------<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Junxiang Chen <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Ph.D. candidate<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Department of chemistry
and molecular<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Wuhan University <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Wuhan, China<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Email: <a
moz-do-not-send="true" href="mailto:cjxxjc729@gmail.com">cjxxjc729@gmail.com</a><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">-------------------------------------------------------<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Junxiang Chen <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Ph.D. candidate<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Department of chemistry
and molecular<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Wuhan University <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Wuhan, China<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Email:
<a class="moz-txt-link-abbreviated" href="mailto:cjxxjc729@gmail.com">cjxxjc729@gmail.com</a><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">QQ: 95508097<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://pwscf.org/mailman/listinfo/pw_forum">http://pwscf.org/mailman/listinfo/pw_forum</a></pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Dr. Lorenzo Paulatto
IdR @ IMPMC -- CNRS & Université Paris 6
+33 (0)1 44 275 084 / skype: paulatz
<a class="moz-txt-link-freetext" href="http://www-int.impmc.upmc.fr/~paulatto/">http://www-int.impmc.upmc.fr/~paulatto/</a>
23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05</pre>
</body>
</html>