<div dir="auto"><div>I am pretty certain that you have done this wrong, double counting the +U. Unless something has changed, the default is to combine the +U with SOC via x lapwso -orb, no "-orb" in lapw1. Check how this was done for the scf iteration (in :log).<br><br><div data-smartmail="gmail_signature">_____<br>Professor Laurence Marks<br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br><a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Jul 18, 2018, 1:21 AM Anup Shakya <<a href="mailto:npshakya31@gmail.com">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>
<div dir="ltr">
<div class="gmail_quote">
<div dir="ltr"><span style="color:rgb(255,0,0)"><span style="color:rgb(0,0,0)">Dear All,</span></span>
<div class="gmail_quote">
<div dir="ltr">
<div><span style="color:rgb(255,0,0)"><span style="color:rgb(0,0,0)"><br>
</span></span></div>
<div><span style="color:rgb(255,0,0)"><span style="color:rgb(0,0,0)">I have performed non magnetic GGA+U+SOC calculations with U =6ev. In the DOS plots in the region between 0 and -1 eV there are three very sharp DOS peaks in both spin up and spin down regions
 because of the contribution of Rare earth 4f state. But when I compare the DOS plots with the bandstructure plots then they are not matching. The three sharp peaks in the DOS region should be visible as three localized bands in the same energy range but I
 am getting these three localized bands above the Fermi level in between 0 to 1.3 eV. For the band structure I have used the Fermi energy in case.insp from case .scf file using the command
<span>grep :FER *.scf. For a non magnetic GGA+U+SO  calculation do we have to subtract the Fermi energy in DOS also from case.qtlup/dn file after using x lapw2 -qtl -up -orb -so -up/dn or is there some other problem with my calculation. The command that I have
 used for obtaining band structure are as follows: </span><br>
</span></span></div>
<span>
<div><span style="color:rgb(255,0,0)"><span style="color:rgb(0,0,0)"><br>
</span></span></div>
<div><span style="color:rgb(0,0,0)">x lapw1 -band -orb -up -p</span></div>
<div><span style="color:rgb(0,0,0)">x lapw1 -band -orb -dn -p</span></div>
<div><span style="color:rgb(0,0,0)">x lapwso -orb -up -p</span></div>
</span>
<div><span style="color:rgb(0,0,0)">copied the EF from file in case.insp file and then did</span></div>
<span>
<div><span style="color:rgb(0,0,0)">x spaghetti -up -orb -so -p</span></div>
<div><span style="color:rgb(0,0,0)"><br>
</span></div>
</span>
<div><span style="color:rgb(0,0,0)">Please provide some suggestions.<br>
</span></div>
<div><span style="color:rgb(0,0,0)"></span><br clear="all">
</div>
<div class="m_-9002515692931328656m_-294828331030561199m_1524528041548694693m_3785448746472339875gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr" style="font-size:12.8px"><font color="#333333"><font face="tahoma, sans-serif">Anup Pradhan Sakhya (Ph.D.)</font></font></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>

</blockquote></div></div></div>