<div dir="ltr"><div dir="ltr"><div>Dear all,</div><div><br></div><div>I have some doubts about plotting the Fermi surface of a GGA+U+SO calculation in Xcrysden. After a self-consistence run with runsp_lapw -p -orb -so, I followed the steps in [1] to calculate the inputs for rendering the Fermi surface, <br></div><div>x lapw1 -up -p</div><div>x lapw1 -dn -p</div><div>x lapwso -orb -up -p</div><div>x lapw2 -so -fermi -up -p</div><div>x lapw2 -so -fermi -dn -p</div><div><br></div><div>Then I concatenated the case.outputso_* files to a single case.output1up file.  In Xcrysden by clicking render the Fermi surface for the spin-up sector, it generates the Fermi surface for half of the number of the bands that cross the Fermi energy. If I want to get the Fermi surface for another half of the bands, is it correct to just simply concatenate the case.outputso* files to case.output1dn and render the Fermi surface for the spin-dn sector in Xcrysden? or some special cares have to be taken?<br></div></div><div dir="ltr"><br><div>Regards,</div><div>Wing Chi<br></div><div><br></div><div>[1] <a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg17837.html">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg17837.html</a></div><div><br></div></div></div>