<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><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:SimSun;
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:SimSun;
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";
color:black;}
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;}
pre
{mso-style-priority:99;
mso-style-link:"HTML \9884\8BBE\683C\5F0F Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";
color:black;}
span.EmailStyle17
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.HTMLChar
{mso-style-name:"HTML \9884\8BBE\683C\5F0F Char";
mso-style-priority:99;
mso-style-link:"HTML \9884\8BBE\683C\5F0F";
font-family:"Courier New";
color:black;}
span.EmailStyle20
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@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]--></head><body bgcolor=white lang=ZH-CN link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>Input file (O) almost the same as O2<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>/<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> &SYSTEM<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> ibrav = 12,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> celldm(1) = 16.0348541363689,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> celldm(2) = 1.52752767,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> celldm(3) = 2.665651575,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> celldm(4) = -0.3273333438252,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> nat = 1,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> ntyp = 1,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> ecutwfc = 45,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> ecutrho = 400,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> occupations = 'smearing' ,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> degauss = 0.01,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> smearing = 'methfessel-paxton' ,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> nspin = 2,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> starting_magnetization(1) = 1,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>/<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> &ELECTRONS<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> electron_maxstep = 200,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> conv_thr = 1.D-6,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> mixing_beta = 0.4,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> mixing_mode = 'local-TF' ,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> /<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> &IONS<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> ion_dynamics = 'bfgs' ,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> /<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>ATOMIC_SPECIES<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> O 15.9994 O.pbe-rrkjus.UPF<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>ATOMIC_POSITIONS crystal<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>O 0.217969593 0.389882242 0.401960243<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>K_POINTS automatic<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> 5 5 2 0 0 0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>nspin and starting_magnetization(1) have been set<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'>thank you for your reply!<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'> Best wishes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='color:#1F497D'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal align=left style='text-align:left'><b><span lang=EN-US style='font-size:11.0pt;color:windowtext'>From:</span></b><span lang=EN-US style='font-size:11.0pt;color:windowtext'> pw_forum-bounces@pwscf.org [mailto:pw_forum-bounces@pwscf.org] <b>On Behalf Of </b>Lorenzo Paulatto<br><b>Sent:</b> Monday, April 28, 2014 4:24 PM<br><b>To:</b> PWSCF Forum<br><b>Subject:</b> Re: [Pw_forum] binding energy of molecular oxygen<o:p></o:p></span></p></div></div><p class=MsoNormal align=left style='text-align:left'><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal align=left style='margin-bottom:12.0pt;text-align:left'><span lang=EN-US>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</span><span lang=EN-US style='font-size:12.0pt'><o:p></o:p></span></p><div><p class=MsoNormal><span lang=EN-US>On 04/28/2014 10:06 AM, Junxiang Chen wrote:<o:p></o:p></span></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><p class=MsoNormal><span lang=EN-US>Dear all,<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 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 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 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><p class=MsoNormal align=left style='text-align:left'><span lang=EN-US style='font-size:12.0pt;font-family:"Times New Roman","serif"'><br><br><br><o:p></o:p></span></p><pre><span lang=EN-US>_______________________________________________<o:p></o:p></span></pre><pre><span lang=EN-US>Pw_forum mailing list<o:p></o:p></span></pre><pre><span lang=EN-US><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><o:p></o:p></span></pre><pre><span lang=EN-US><a href="http://pwscf.org/mailman/listinfo/pw_forum">http://pwscf.org/mailman/listinfo/pw_forum</a><o:p></o:p></span></pre></blockquote><p class=MsoNormal align=left style='text-align:left'><span lang=EN-US style='font-size:12.0pt;font-family:"Times New Roman","serif"'><br><br><o:p></o:p></span></p><pre><span lang=EN-US>-- <o:p></o:p></span></pre><pre><span lang=EN-US>Dr. Lorenzo Paulatto<o:p></o:p></span></pre><pre><span lang=EN-US>IdR @ IMPMC -- CNRS & Université Paris 6<o:p></o:p></span></pre><pre><span lang=EN-US>+33 (0)1 44 275 084 / skype: paulatz<o:p></o:p></span></pre><pre><span lang=EN-US><a href="http://www-int.impmc.upmc.fr/~paulatto/">http://www-int.impmc.upmc.fr/~paulatto/</a><o:p></o:p></span></pre><pre><span lang=EN-US>23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05<o:p></o:p></span></pre></div></body></html>