<div dir="auto"><div>A very basic question for you -- what is your surface structure? In general, just taking the bulk and cutting it is wrong (too lazy).<div dir="auto"><br></div><div dir="auto">Before what you said:</div><div dir="auto">a) Optimize the volume of the bulk cell with your functional of choice, no so or mbj.</div><div dir="auto">b) Create a representative (e.g. valence neutral) surface slab using the 2D lattice parameters in the surface plane from a).</div><div dir="auto">c) Optimise the internal positions (-min) as there will be some relaxation. Again no so or mbj. Check that you have, for instance, an insulator. You may have to recreate the surface if you did it wrong. Depending on case you will use run_lapw, runsp_c_lapw or runsp_lapw. Adding DW is optional.</div><div dir="auto"><br></div><div dir="auto">Then start with adding mbj & so on. You have to keep your fingers crossed that the structure from c) is a valid representation, since you cannot optimize positions with so & mbj.</div><div dir="auto"><br></div><div dir="auto">Depending upon care, 50% probability that you surface structure will be different from nature.</div><div dir="auto"><br><div data-smartmail="gmail_signature" dir="auto">--<br>Professor Laurence Marks (Laurie)<br>Department of Materials Science and Engineering, Northwestern University<br><a href="http://www.numis.northwestern.edu">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><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jun 20, 2023, 07:02 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" style="word-wrap:break-word"><div class="m_5195853539771538025WordSection1"><p class="MsoNormal">Dear experts, I am having a compact cluster with 100gb RAM, 36 cores and 12 TB SATA hard disk with CentOS.<u></u><u></u></p><p class="MsoNormal">I want to calculate the band structure at the surface of topological insulator. Below are the steps I have used for the detailed calculation.<u></u><u></u></p><p class="MsoNormal">After initialization, I run a pbe scf using <u></u><u></u></p><p class="MsoNormal">“run_lapw -p -ec 0.0001 -cc 0.001 -NI”<u></u><u></u></p><p class="MsoNormal">Save the result using <u></u><u></u></p><p class="MsoNormal">“save_lapw <b>pbe</b>”<u></u><u></u></p><p class="MsoNormal">In the next step <u></u><u></u></p><p class="MsoNormal">“init_so_lapw” ( took the default parameters, magnetization is 0 0 1)<u></u><u></u></p><p class="MsoNormal">“Init_nl_vdw “ (kernel type 1)</p><p class="MsoNormal">“init_mbj_lapw”</p><p class="MsoNormal">“run_lapw -p -so -i 1 -NI”</p><p class="MsoNormal">“save_lapw <b>mbj</b>”</p><p class="MsoNormal">“init_mbj_lapw”</p><p class="MsoNormal">Took the selection 2 for surface calculation, and lmbj for potential </p><p class="MsoNormal">And finally I run “run_lapw -p -lmbj -so -ec 0.0001 -cc 0.001 -i 999 -NI”</p><p class="MsoNormal">And then I save the result using save_lapw <b>lmbj</b>.</p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">I want to know whether these are the correct steps for calculating band structure at the surface or I have to used DFT-D3 for the correction of dispersion which are basically used for the surface calculation. If DFT-D3 must be used for the surface calculation to show Dirac Cone/band inversion nature, then what modification I have to do in my steps to use DFT-D3 so that I can get a clear Band structure with Dirac Cone.</p><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Thanks in Advance </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" rel="noreferrer">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer 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 noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div></div></div>