<html><head></head><body><div style="color:#000; background-color:#fff; font-family:標楷體, dfkai-sb;font-size:16px"><div id="yui_3_16_0_ym19_1_1509705151887_2898"><span id="yui_3_16_0_ym19_1_1509705151887_3005">Dear Sir Peter/Stefaan/Gerhard/Gavin,</span></div><div id="yui_3_16_0_ym19_1_1509705151887_3085"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3278"><span id="yui_3_16_0_ym19_1_1509705151887_3005">Thanks you all for your helpful suggestions.</span></div><div id="yui_3_16_0_ym19_1_1509705151887_3214"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3215"><span id="yui_3_16_0_ym19_1_1509705151887_3005">I could manage to calculate Si atomization energy for both sp and without -sp case.</span></div><div id="yui_3_16_0_ym19_1_1509705151887_3277"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div dir="ltr" id="yui_3_16_0_ym19_1_1509705151887_3348"><span id="yui_3_16_0_ym19_1_1509705151887_3005">Si_atom_nsp Ene=-1157.44977178 Ry</span></div><div dir="ltr" id="yui_3_16_0_ym19_1_1509705151887_3355"><br></div><div dir="ltr" id="yui_3_16_0_ym19_1_1509705151887_3356"><span id="yui_3_16_0_ym19_1_1509705151887_3357">Si_atom_sp Ene= -1157.44961225 Ry</span></div><div dir="ltr" id="yui_3_16_0_ym19_1_1509705151887_3351"><span id="yui_3_16_0_ym19_1_1509705151887_3005"></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3262"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3483"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div dir="ltr" id="yui_3_16_0_ym19_1_1509705151887_3276"><span id="yui_3_16_0_ym19_1_1509705151887_3005">and the Cohesive Energy I obtained is: 8.037453341 eV per f.u.</span></div><div dir="ltr" id="yui_3_16_0_ym19_1_1509705151887_3648"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3083"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3841"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3649"><span id="yui_3_16_0_ym19_1_1509705151887_3005">But I am not able to run O2 mulecule.</span></div><div id="yui_3_16_0_ym19_1_1509705151887_3703"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3704"><span id="yui_3_16_0_ym19_1_1509705151887_3005">Below os my O2 structure</span></div><div id="yui_3_16_0_ym19_1_1509705151887_3908"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3705"><span id="yui_3_16_0_ym19_1_1509705151887_3005"><br></span></div><div id="yui_3_16_0_ym19_1_1509705151887_3008"><span></span></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: 標楷體, dfkai-sb; font-size: 16px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, Sans-Serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> On Friday, 3 November 2017 2:23 PM, Stefaan Cottenier <Stefaan.Cottenier@UGent.be> wrote:<br></font></div> <br><br> <div class="y_msg_container"><div dir="ltr">Hello Gerhard,<br clear="none"><br clear="none">I'm reasoning by analogy now, not by deep understanding, but indeed, if you would apply Hund's rule to the relativistic p-shell, it looks to me it would be as you describe. The p-1/2 would then play the same role as an s-orbital for Hund's rules.<br clear="none"><br clear="none">Anyway, this is not relevant for wien2k free atom calculations because (1) Si and Ge are too light for the fully relativistic quantum numbers to apply, and (2) the calculations are done at the scalar-relativistic level.<br clear="none"><br clear="none">Best,<br clear="none">Stefaan<br clear="none"><br clear="none"><br clear="none">> -----Oorspronkelijk bericht-----<br clear="none">> Van: Wien [mailto:<a shape="rect" ymailto="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>] Namens<br clear="none">> Fecher, Gerhard<br clear="none">> Verzonden: vrijdag 3 november 2017 9:29<br clear="none">> Aan: A Mailing list for WIEN2k users <<a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>><br clear="none">> Onderwerp: Re: [Wien] open shell case [Si and Ge]<br clear="none">> <br clear="none">> Dear Stefaan<br clear="none">> What would happen if you treat all electrons relativistic, I would expect that<br clear="none">> the two p valence electrons occupy the p1/2,+1/2 and p1/2,-1/2 states<br clear="none">> leaving the 4 p3/2 states unoccupied and there would be no spinpolarisation<br clear="none">> <br clear="none">> I did not check, maybe I am wrong.<br clear="none">> <br clear="none">> Ciao<br clear="none">> Gerhard<br clear="none">> <br clear="none">> DEEP THOUGHT in D. Adams; Hitchhikers Guide to the Galaxy:<br clear="none">> "I think the problem, to be quite honest with you, is that you have never<br clear="none">> actually known what the question is."<br clear="none">> <br clear="none">> ====================================<br clear="none">> Dr. Gerhard H. Fecher<br clear="none">> Institut of Inorganic and Analytical Chemistry Johannes Gutenberg -<br clear="none">> University<br clear="none">> 55099 Mainz<br clear="none">> and<br clear="none">> Max Planck Institute for Chemical Physics of Solids<br clear="none">> 01187 Dresden<br clear="none">> ________________________________________<br clear="none">> Von: Wien [<a shape="rect" ymailto="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>] im Auftrag von<br clear="none">> Stefaan Cottenier [<a shape="rect" ymailto="mailto:Stefaan.Cottenier@UGent.be" href="mailto:Stefaan.Cottenier@UGent.be">Stefaan.Cottenier@UGent.be</a>]<br clear="none">> Gesendet: Freitag, 3. November 2017 06:52<br clear="none">> An: A Mailing list for WIEN2k users<br clear="none">> Betreff: Re: [Wien] open shell case [Si and Ge]<br clear="none">> <br clear="none">> Following Hund's rules, the atomic ground state for Si and Ge is spin-<br clear="none">> polarized ('open shell'):<br clear="none">> <br clear="none">> <a shape="rect" href="https://www.webelements.com/silicon/atoms.html" target="_blank">https://www.webelements.com/silicon/atoms.html</a><br clear="none">> <br clear="none">> Hence, yes, you need -sp to find a meaningful value for the free atom total<br clear="none">> energy.<br clear="none">> <br clear="none">> Stefaan<br clear="none">> <br clear="none">> <br clear="none">> Van: Wien [mailto:<a shape="rect" ymailto="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>] Namens<br clear="none">> Gavin Abo<br clear="none">> Verzonden: vrijdag 3 november 2017 5:42<br clear="none">> Aan: <a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a><br clear="none">> Onderwerp: Re: [Wien] open shell case [Si and Ge]<br clear="none">> <br clear="none">> Feel free to correct me if I'm wrong, but I think Si and Ge have an even<br clear="none">> number of electrons (or paired electrons).<br clear="none">> <br clear="none">> For Si, 4 electrons in the 3s^2 3p^2. For Ge, 4 electrons in the 4s^2 4p^2. [1]<br clear="none">> <br clear="none">> This making them closed shell [2].<br clear="none">> <br clear="none">> In Prof. Blaha's example for free atoms [3], spin-polarized (runsp_lapw) is<br clear="none">> used for open shell (or non-closed shell [4]) Li (1 electron in 2s^1) and B (3<br clear="none">> electron in 2s^2 2p^1) and non-spin polarized is used for closed shell Be (2<br clear="none">> electron in 2s^2).<br clear="none">> <br clear="none">> I don't know why it is fluctuating, diverging? Is the box you used containing<br clear="none">> the free atom large enough? I remember that 12x12x12 angstrom usually<br clear="none">> might be large enough [5].<br clear="none">> <br clear="none">> [1]<br clear="none">> <a shape="rect" href="https://en.wikipedia.org/wiki/Electron_configurations_of_the_elements_" target="_blank">https://en.wikipedia.org/wiki/Electron_configurations_of_the_elements_</a>(da<br clear="none">> ta_page)<br clear="none">> [2] <a shape="rect" href="https://www.mail-/" target="_blank">https://www.mail-</a><br clear="none">> archive.com/<a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>/msg11320.html<br clear="none">> [3] <a shape="rect" href="https://www.mail-/" target="_blank">https://www.mail-</a><br clear="none">> archive.com/<a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>/msg16691.html<br clear="none">> [4] <a shape="rect" href="http://susi.theochem.tuwien.ac.at/reg_user/faq/cohesive.html" target="_blank">http://susi.theochem.tuwien.ac.at/reg_user/faq/cohesive.html</a><br clear="none">> [5] <a shape="rect" href="https://www.mail-/" target="_blank">https://www.mail-</a><br clear="none">> archive.com/<a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>/msg12815.html<br clear="none">> On 11/2/2017 7:06 PM, chin Sabsu wrote:<br clear="none">> Dear Peter Sir,<br clear="none">> <br clear="none">> Do Si [S^23P^4] and Ge [4S^24P^2] are also an open shell case?<br clear="none">> <br clear="none">> Because my Si atomization energy is fluctuating without -sp switch after 55<br clear="none">> scf cycles also.<br clear="none">> <br clear="none">> <br clear="none">> <br clear="none">> Sincerely<br clear="none">> <br clear="none">> Chin<br clear="none">> _______________________________________________<br clear="none">> Wien mailing list<br clear="none">> <a shape="rect" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br clear="none">> <a shape="rect" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br clear="none">> SEARCH the MAILING-LIST at: <a shape="rect" href="http://www.mail-/" target="_blank">http://www.mail-</a><br clear="none">> archive.com/<a shape="rect" ymailto="mailto:wien@zeus.theochem.tuwien.ac.at" href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>/index.html<div class="yqt2549219375" id="yqtfd70941"><br clear="none">_______________________________________________<br clear="none">Wien mailing list<br clear="none"><a shape="rect" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br clear="none"><a shape="rect" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br clear="none">SEARCH the MAILING-LIST at: <a shape="rect" href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br clear="none"></div></div><br><br></div> </div> </div> </div></div></body></html>