<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">Some comments:</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">a) 99.999% probability that you messed something up when you created the supercell, or you have FOR in case.in2(c) in one case and not the other.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">b) Putting 4f into the core is at best dubious. What valence are you using, for instance Ce4+, Ce3+ or Ce0+ ? You are deciding what the electronic structure is. While you might be right,....</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">c) -eece, as I find +U to be dubious for 4f.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Sep 3, 2021 at 10:40 AM Artem Tarasov <<a href="mailto:artem.tarasov@spbu.ru">artem.tarasov@spbu.ru</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">Dear wien2k users,<br>
<br>
Recently, I have noticed a strange problem related to the forces <br>
inside the supercell for compounds containing 4f-atoms with 4f-states <br>
placed on the core level ("open core approximation" for Tb, Dy, Ho, Ce <br>
etc.). If we relax the slab of a cell with a period of 1x1 and a <br>
symmetry group of P4mm, and then make a cell with a doubled lateral <br>
period of 2x2 from the relaxed structure, then after converging the <br>
density for the 2x2 cell the forces on all atoms in the supercell grow <br>
significantly compared to the case of the 1x1 cell. Up to a jump from <br>
0.05 to 68 mRy/bohr on individual atoms. I thought that the problem <br>
may be related to the shift of the k-grid, because when I move to the <br>
supercell, I change the k-grid from 12x12 to 6x6, and I had supposed <br>
that such a k-grid may not be enough for an accurate calculation. <br>
However, even if I don't use the k-grid shift in both cases, this <br>
problem persists. My calculation is a slab of CeB6 with symmetry P4mm <br>
with Ce4f states given in case. inc as<br>
<br>
16 0.70  0  NUMBER OF ORBITALS (EXCLUDING SPIN), SHIFT, IPRINT<br>
…<br>
4, 3,0.43               ( N,KAPPA,OCCUP)<br>
4,-4,0.57               ( N,KAPPA,OCCUP)<br>
<br>
Thank you in advance for your help.<br>
<br>
Best regards,<br>
Artem Tarasov<br>
Department of Solid State Electronics<br>
Saint Petersburg University<br>
_______________________________________________<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="https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!FheN273SY3aaOingm6LKFziHLi0fvV1CEUYdiE7uTRRjc0mhHXMgh_Z5FU3hj1THOiEnew$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!FheN273SY3aaOingm6LKFziHLi0fvV1CEUYdiE7uTRRjc0mhHXMgh_Z5FU3hj1THOiEnew$</a> <br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!FheN273SY3aaOingm6LKFziHLi0fvV1CEUYdiE7uTRRjc0mhHXMgh_Z5FU3hj1RitPehgw$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!FheN273SY3aaOingm6LKFziHLi0fvV1CEUYdiE7uTRRjc0mhHXMgh_Z5FU3hj1RitPehgw$</a> <br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Professor Laurence Marks<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><div>"Research is to see what everybody else has seen, and to think what nobody else has thought" Albert <span style="font-family:Arial,Helvetica,sans-serif;font-size:12.8px">Szent-</span><span style="font-family:Arial,Helvetica,sans-serif;font-size:small;color:rgb(34,34,34)">Györgyi</span></div></div></div>