<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:large" class="gmail_default">Dear Dr. Tran,</div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default">I tested two different approaches for mBJ+SO  mentioned at [1]</div><div style="font-size:large" class="gmail_default">In the present case the ground state energy differed by an amount  0.00062938Ry and the band gap by a number 0.001eV which is as per the difference in the :ENE.</div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default">I know this difference is negligible and can be ignored but could you please make any comment on it, why this difference occurred? <br></div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default">I did the calculation in Wien2k_18.2, i5 machine with composer_xe_2015.0.090 ifort and CC compiler.</div><div style="font-size:large" class="gmail_default">The system is binary and no RLO.<br></div><div style="font-size:large" class="gmail_default"></div><div style="font-size:large" class="gmail_default"><span style="color:rgb(0,0,255)">pbe then mbj with -SO</span><br>TOTAL ENERGY IN Ry =       -43189.10920210<br>GAP =  0.006765 Ry =     0.092 eV<br><br><span style="color:rgb(0,0,255)">pbe followed by  mbj and then followed by scf with -SO</span><br>TOTAL ENERGY IN Ry =       -43189.10857272<br>GAP = 0.006866 Ry =     0.093 eV</div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default"><b>mBJ parameters:</b></div><div style="font-size:large" class="gmail_default"> -1.200000000000000E-002<br>   1.02300000000000<br>  0.500000000000000</div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default"><b>case.inso</b></div><div style="font-size:large" class="gmail_default">WFFIL<br>4  0  0                 llmax,ipr,kpot<br>-10  1.9                Emin, Emax<br>    0 0 1                           h,k,l (direction of magnetization)<br> 0                       number of atoms with RLO<br>1 2      number of atoms without SO, atomnumbers<br></div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default"><b>scf information</b><br></div><div style="font-size:large" class="gmail_default">rkmax=7, PBE, mesh size 10 10 10 (shifted, however no change with shifted or non-shifted k-mesh).</div><div style="font-size:large" class="gmail_default">-ec 0.0001 -cc 0.001<br></div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default">[1] <a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg18427.html">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg18427.html</a></div><div style="font-size:large" class="gmail_default"><br></div><div style="font-size:large" class="gmail_default">regards</div><div style="font-size:large" class="gmail_default">Bhamu<br></div></div></div></div></div></div></div></div>