<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">Others will probably give you suggestions about converging mBJ. Some deeper comments.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">The most common reason that calculations behave badly is user error. Sometimes this is doing the initialization wrong, often it is creating an inappropriate model. Just because one can use x supercell or related codes to create atomic positions does <i>not</i> make them sensible.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">Your slab has a composition Bi12 Te19, i.e. it is Te rich. You therefore have a reduced, n-type semiconductor with a small gap (about 0.1eV) which will behave badly. If you look at your BVS you will see that atom Te7 is severely underco-ordinated (unstable).</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">I doubt that your slab will ever converge to anything which should be published. Surfaces are not simple, it is unfortunately too easy to set something up and get it wrong. (There are many papers in the literature of dubious merit.)</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 21, 2023 at 10:47 PM Burhan Ahmed <<a href="mailto:burhan.ahmed@aus.ac.in">burhan.ahmed@aus.ac.in</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg2255841863679317883"><div lang="EN-US" style="overflow-wrap: break-word;"><div class="m_2255841863679317883WordSection1"><p class="MsoNormal">Dear experts, I am doing an scf calculation taking lmbj potential. My system is a slab of 6ql with a vacuum of 40 ang along c-axix. Whenever I try to run the scf calculation including SOC, after 999 iteration the scf is still not converged. I have analyzed the scf file, it shows the fluctuating nature. At first the selected Rmt was 2.5 then I reduces it to 2.34 and again convergence failed after 999 cycle. I have included -hdlo and -lvns 8 switch in init_lapw as it shows heavy atom. I am attaching the case.struct file. Hoping for any suggestion/solution. </p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Regards<u></u><u></u></p><p class="MsoNormal">Burhan Ahmed<u></u><u></u></p><p class="MsoNormal"><b>Research Scholar, AUS <u></u><u></u></b></p><p class="MsoNormal"><u></u> <u></u></p></div></div>_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</div></blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr">Professor Laurence Marks (Laurie)<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a><br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Györgyi</div></div>