<div dir="auto">Dear Prof Blaha,<div dir="auto"><br></div><div dir="auto">I have performed non magnetic spin polarized calculations using runsp_c_lapw -orb. I have not shifted the k-mesh for FS calculations. I couldn't use the command x kgen -so for FS calculations as I got the following error. </div><div dir="auto">forrtl: severe (24): end of file during read, unit 20, file/case.ksym</div><div dir="auto">error command: /home/kbmaiti/wien2k/kgen kgen.def failed.</div><div dir="auto">So I used x kgen 30000 kpoints no shift</div><div dir="auto">x lapw1 - up</div><div dir="auto">x lapwso </div><div dir="auto">cp case.outputso case.output1</div><div dir="auto">x lapw2 - so - fermi </div><div dir="auto">I don't think there is any problem with the FS calculations as it is matching with the bandstructure calculations. The only problem is the lowering of the symmetry with SOC from 48 to 16 number of symmetry operations though the structure is still cubic for magnetization direction along 0 0 1 axis. Due to this lowering of symmetry the band structure plot along G-X and X-G is different. But experimentally if one measures the bandstructure or the FS of this material it should be the same because experimentally the structure is cubic. So my question is whether non magnetic spin polarized calculations with U and SOC is at all reliable to match or analyze the experimental results. The only reason for doing non magnetic spin polarized calculations was to take into account the strong correlation of R 4f states. I am confused. Please clarify. </div><div dir="auto"><br></div><div dir="auto">Sincerely,</div><div dir="auto">Anup Pradhan Sakhya</div><div dir="auto"><br></div><div dir="auto"><br><div dir="auto"><br></div><div dir="auto"><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Nov 5, 2018, 1:23 PM Anup Shakya <<a href="mailto:npshakya31@gmail.com" target="_blank" rel="noreferrer">npshakya31@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 All,<div>I have done non-magnetic GGA+SOC+U calculations with U and J for a material containing rare earth element. Without SOC and U the Fermi surface is 4 fold symmetric but when I performed GGA+SOC+U calculations the Fermi surface is not 4 fold symmetric. Before applying SOC the position of the atoms were as follows:</div><div>Atom X Mult =1  0.0      0.0       0.0</div><div>Atom Y Mult 6  0.207   0.5      0.5</div><div>After SOC the a,b,c values and the symmetry i.e, Pm3m is the same as before but now the Y positions has split into two non-equivalent atoms as shown below.</div><div>Atom Y 2 Mult 4 0.207   0.5  0.5</div><div>Atom Y 3 Mult 2 0.5      0.5   0.207.</div><div>The position and the multiplicity of X atom is unchanged. I think probably because of SOC the equivalent atoms have changed into non-equivalent atoms and led to the anisotropy in the band structure along the G-X and X-G direction and also 4 fold symmetry is not there in the Fermi surface calculations. But if I want to compare the obtained band structure or the Fermi surface with the ARPES results then it doesn't make sense because the Fermi surface should be 4 fold symmetric. I do not know what to do. If you need any more information then please let me know. Please provide some suggestions as to preserve the 4 fold symmetry along with GGA+SOC+U calculations. Any suggestion would be very much valuable for me. </div><div><br></div><div>Sincerely,</div><div>Anup Pradhan Sakhya,</div><div>TIFR, India</div><div><br></div><div> <br clear="all"><div><br></div></div></div>
</blockquote></div>