<div dir="auto"><div>40 iterations is not much for a problem like this. I have no idea why your calculation went strange at the end of the first run, and I doubt that anyone can guess.</div><div dir="auto"><br></div><div dir="auto">First, I assume that you bulk material is an insulator. Is your slab also an insulator -- look at the top of case.scf2up? It is easy to get the number of atoms wrong. I suggest checking the stoichiometry.</div><div dir="auto"><br></div><div dir="auto">Second, did you check the bond valence sums? Depending upon what the atoms are you can get these from "grep Bond *outputnn" -- they are produced by x nn.</div><div dir="auto"><br></div><div dir="auto">Third, why are you using the bulk lattice parameters? Particularly for surfaces you need to use the DFT ones, especially in the surface plane.</div><div dir="auto"><br></div><div dir="auto">Fourth, why are you not relaxing? Not only is this right, big movements will tell you that your slab is wrong. Note: you cannot relax with -so.</div><div dir="auto"><br></div><div dir="auto">I would not expect a TMDC to be hard, but I don't know details. Of course mistakes in the model make things hard.<br><br><div data-smartmail="gmail_signature" dir="auto">--<br>Professor Laurence Marks<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><br><div class="gmail_quote" dir="auto"><div dir="ltr" class="gmail_attr">On Mon, Feb 20, 2023, 05:09 pluto via Wien <<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Prof. Marks, dear All,<br>
<br>
I am calculating a large slab of one of the 1T TMDC materials, 30 unique <br>
atoms (this makes it 20 tri-layers, i.e. 20 vdW layers).<br>
<br>
The slab of the same material of half size, 15 unique atoms (10 <br>
tri-layers) converged well.<br>
<br>
I am using TEMP 0.002 (first I tried with TETRA, but restated with TEMP <br>
after some iterations). It is a spin-polarized calculation with SOC. <br>
Inversion is present. Experimental lattice constant, slab made with "x <br>
supercell", no relaxation. Slab looks good in xcrysden ("x supercell" <br>
sometimes puts surface atoms wrongly, but I corrected it.). I used 30 <br>
Bohr vacuum between layers. I use 36 k-points (klist pasted below). <br>
Other parameters are WIEN2k default.<br>
<br><br><br>
</blockquote></div></div></div>