<div dir="auto"><div>Open core puts the 4f into the core as spherically symmetric wavefunctions (with spin orbit). Therefore it will get you a reference state, but it will not give you the crystal field splitting.<div dir="auto"><br></div><div dir="auto">To get the splitting you have to use LDA+U (which I don't trust for 4f) or on-site hybrid (better).</div><div dir="auto"><br></div><div dir="auto">N.B., </div><div dir="auto">1) Has your calculation converged? </div><div dir="auto">2) Did you correctly alter the number of atoms -- i.e. grep :NEC *.scf?</div><div dir="auto">3) Are you running spin polarized with correct case.incup/dn (which does not seem to be described in the FAQ)? I would not trust a non spin polarized calculation much.</div><div dir="auto">4) Have you run a number of other cases first to learn? This is not trivial, I would suggest running more than TiC first, for instance NiO and some others.<br><br><div data-smartmail="gmail_signature" dir="auto">_____<br>Professor Laurence Marks<br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br><a href="http://www.numis.northwestern.edu" rel="noreferrer noreferrer" target="_blank">www.numis.northwestern.edu</a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Nov 5, 2019, 07:36 Abhijit B K <<a href="mailto:bk.abhijit1@gmail.com" rel="noreferrer noreferrer" target="_blank">bk.abhijit1@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Wien2k users,<div><br></div><div>I am a beginner for WIEN2k or any kind of DFT calculations. I am trying to follow the formalism by Pavel Novak to calculate the <b>crystal field parameters</b> of a rare earth oxide insulator. I am stuck in the first step i.e <b>open core treatment</b>. I follow all the steps mentioned in the manual as well as FAQ and I don't sift the k-mesh. What I see in the .scf file<b> </b>is <b>metallic behavior instead of insulator. </b>i.e</div><div><br></div><div><span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px">:GAP (global)   :    0.0    Ry =     0.0   eV (metal)</span>   <br></div><div><b> </b>I understand this as a product  of screening.  <br></div><div><br></div><div>Can someone suggest some ways to achieve a<b> non zero gap</b>?</div><div>I wouldn't want to use LDA+U.</div><div><br></div><div>Best wishes,</div><div><br></div><div>Abhijit  </div></div>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" rel="noreferrer noreferrer noreferrer" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=dV2WxxTROR_G-ZQ9n16BLC8F6upgKFBEM4nejgNfkMc&s=Cw4JVHx0eXJF4I_JGNIyJY0XT6CWaJ5HmT6wqVYhc8Y&e=" rel="noreferrer noreferrer noreferrer noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=dV2WxxTROR_G-ZQ9n16BLC8F6upgKFBEM4nejgNfkMc&s=Cw4JVHx0eXJF4I_JGNIyJY0XT6CWaJ5HmT6wqVYhc8Y&e=</a> <br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=dV2WxxTROR_G-ZQ9n16BLC8F6upgKFBEM4nejgNfkMc&s=d56fuMDDJYzQBZ68t9ogb3proq7ZUcapIZsLE1Gw3rU&e=" rel="noreferrer noreferrer noreferrer noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=dV2WxxTROR_G-ZQ9n16BLC8F6upgKFBEM4nejgNfkMc&s=d56fuMDDJYzQBZ68t9ogb3proq7ZUcapIZsLE1Gw3rU&e=</a> <br>
</blockquote></div>
</div></div>