<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>