<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi,<div class=""><br class=""></div><div class="">first of all a little note: unless my e-mail server is going crazy, I think it is the sixth e-mail on the same topic sent in two or three days. This is a forum, where people voluntarily share their (little or big) knowledge and experience</div><div class="">with others. As such, once a question is made, it would be kind (this is my opinion, of course!) to wait for someone that, after reading the message, finds some time to spent to provide an answer. Posting multiple messages,</div><div class="">usually delays the answer!!! </div><div class=""><br class=""></div><div class="">Moreover, I think it is quite inefficient to post messages with no subject, since this does not allow to immediately identify which reply concerns a given thread.</div><div class=""><br class=""></div><div class="">This being said, this is what I could notice as far as your question is concerned:</div><div class=""><br class=""></div><div class="">1) your original question, that is, “is -92.306…… Ry” the energy of an Na atom?, is not really well posed, since the total energy depends on cutoff, pseudo potential, and so on, so just the value</div><div class="">of the total energy might be meaningless, much better to share input file (as you did later) to understand how such an energy is calculated</div><div class=""><br class=""></div><div class="">2) if you need such energy to calculate the adsorption energy, you should decide if your reference is an hypothetic gas of Na atoms or bulk Na, in the latter case you should calculate the total</div><div class="">energy per Na atom of bulk Na, instead</div><div class=""><br class=""></div><div class="">3) you are using an ultra soft pseudo potential, that means that the default value of ecutrho=4*ecutwfc might be insufficient or, in some cases, even provide truly wrong results. As such, you should check</div><div class="">the convergence also with respect to ecutrho, in any run where ultra soft pseudos are involved. Consider that even the cutoff is maybe too small. In the pseudo file you’re using, you should find two lines</div><div class="">like:</div><div class=""> Suggested minimum cutoff for wavefunctions: 66. Ry<br class=""> Suggested minimum cutoff for charge density: 323. Ry</div><div class="">So 35 Ry, in comparison, looks really small!</div><div class=""><br class=""></div><div class="">4) as already suggested, If the total energy of an isolated Na atom is exactly what you need, you should also look for a spin polarised ground state, to check if including the spin degree of freedom leads</div><div class="">to a lower energy ground state</div><div class=""><br class=""></div><div class="">5) again, If the total energy of an isolated Na atom is exactly what you need, your unit cell is likely to be too small, only 4 A separate two neighbour periodic replicas of the Na atom, an this is not enough</div><div class="">to ensure no interaction between the periodic replicas</div><div class=""><br class=""></div><div class="">6) again, If the total energy of an isolated Na atom is exactly what you need, you don’t have to sample the Brillouin zone with an 8x8x1 grid, because if the atom is truly isolated (see previous point),</div><div class="">the energy levels are flat, so just gamma point is enough</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Giovanni</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 15 Jun 2017, at 17:03, nipesh dulal <<a href="mailto:nipeshdulal25@gmail.com" class="">nipeshdulal25@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div style="color:rgb(49,49,49);word-spacing:1px" dir="auto" class=""><div class=""><div class="">Hello Everyone <br class=""><br class=""></div>I have calculated the ground state energy of Sodium Atom. The final result is -92.3063689424Ry. I have use below input file. Is it correct?<br class=""><br class=""> &control<br class=""> calculation='relax',<br class=""> verbosity='high'<br class=""> prefix='sodiumrelax'<br class=""> outdir='/home/quantumespresso/Downloads/graphene/adsorption'<br class=""> pseudo_dir = '/home/quantumespresso/Downloads/graphene/pseudopotentials',<br class=""> tprnfor=.true.<br class=""> forc_conv_thr=1.0d-3<br class=""> tstress=.true.<br class=""> etot_conv_thr=1.0d-4 <br class=""> <br class=""> <br class=""> /<br class=""> &SYSTEM<br class=""> ibrav=3<br class=""> celldm(1)=8.906<br class=""> nat=1<br class=""> ntyp=1<br class=""> ecutwfc=35<br class=""> occupations='smearing'<br class=""> smearing='mv'<br class=""> nbnd =12<br class=""> degauss=0.002<br class=""> /<br class=""> &ELECTRONS<br class=""> electron_maxstep = 100<br class=""> mixing_beta=0.6<br class=""> conv_thr = 1.0D-8<br class=""> /<br class=""> &ions<br class=""> /<br class=""> ATOMIC_SPECIES<br class=""> Na 23.0 Na.pbe-spn-rrkjus_psl.0.2.UPF<br class="">ATOMIC_POSITIONS (angstrom)<br class=""> Na 0.000000000 0.000000000 0.0000000<br class=""> <br class=""> K_POINTS {automatic}<br class=""> 8 8 1 0 0 0<br class=""><br class=""><br class=""></div>With Best Regards<br class=""></div><span style="color:rgb(49,49,49);word-spacing:1px;background-color:rgb(255,255,255)" class=""><div dir="auto" class="">NIpesh Dulal </div><div dir="auto" class="">Kathmandu,Nepal</div><div dir="auto" class="">Trihbhuwan university </div></span>
_______________________________________________<br class="">Pw_forum mailing list<br class=""><a href="mailto:Pw_forum@pwscf.org" class="">Pw_forum@pwscf.org</a><br class="">http://pwscf.org/mailman/listinfo/pw_forum</div></blockquote></div><br class=""><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">-- <br class=""><br class="">Giovanni Cantele, PhD<br class="">CNR-SPIN<br class="">c/o Dipartimento di Fisica<br class="">Universita' di Napoli "Federico II"<br class="">Complesso Universitario M. S. Angelo - Ed. 6<br class="">Via Cintia, I-80126, Napoli, Italy<br class="">e-mail: <a href="mailto:giovanni.cantele@spin.cnr.it" class="">giovanni.cantele@spin.cnr.it</a><br class="">Phone: +39 081 676910<br class="">Skype contact: giocan74<br class=""><br class="">ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" class="">http://www.researcherid.com/rid/A-1951-2009</a><br class="">Web page: <a href="http://people.fisica.unina.it/~cantele" class="">http://people.fisica.unina.it/~cantele</a><br class=""></div>
</div>
<br class=""></div></body></html>