<div dir="ltr"><div dir="ltr"><span style="color:rgb(0,0,0)">Hello,</span></div><div dir="ltr"><span style="color:rgb(0,0,0)"><br></span></div><div dir="ltr"><span style="color:rgb(0,0,0)">I just finished some calculation with PBE/PBE+SO/mBJ+SO.</span></div><div><span style="color:rgb(0,0,0)">I know that with mBJ we can not do force optimization.</span></div><div><span style="color:rgb(0,0,0)">I know how to do force optimization with -SO</span></div><div><br></div><div><font color="#000000">I have done PBE+mBJ calculation </font></div><div dir="ltr"><span style="color:rgb(0,0,0)">save_lapw -d case</span></div><div><span style="color:rgb(0,0,0)">initso_lapw and </span><span style="color:rgb(0,0,0)">fixed position </span>for<span style="color:rgb(0,0,0)"> some atoms and then </span></div><div><font color="#000000">started min -j "run_lapw -ec xxx -cc xxx -fc x -i x"</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">My first question:</font></div><div><font color="#000000">Is the above process is okay, can we do that? </font></div><div><font color="#000000">From above method, my band gap is in the reasonable agreement with Experimental band gap (still optimization is going on).</font></div><div><font color="#000000"><br></font></div><div><span style="color:rgb(0,0,0)">As we know the standard process for the mBJ is:</span><br></div><div><font color="#000000">1.init_mbj_lapw</font></div><div><font color="#000000">2. run_lapw -i 1 -NI</font></div><div><font color="#000000">3. save_lapw **</font></div><div><font color="#000000">4. init_mbj_lapw</font></div><div><font color="#000000">5. 0/1/2/3</font></div><div><font color="#000000">6. run_lapw .......</font></div><div><br></div><div><font color="#000000">So my second question: </font></div><div><font color="#000000">If we do a general mBJ and then do -SO calculation after save_lapw in mBJ directory then -so will be assumed from the step 6 only. What would be the effect of step 2 "</font><span style="color:rgb(0,0,0)">run_lapw -</span>i<span style="color:rgb(0,0,0)"> 1 -NI" (that have been done in the case of mBJ and the -SO assumed from step 6).</span></div><div><span style="color:rgb(0,0,0)">From above two question, I am not able to find a proper way to follow (mBJ)+(SO) or (mBJ+SO)</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)">We can do in two ways if I understood it correctly,  first run an mBJ and then do scf with -SO is different then doing mBJ+SO in a single scf step. From one we can relax the structure while from two we can not.</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)">Hope I have framed this draft in detail and you will be able to understand what I want to know.</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div>Sincerely</div><div>Fatima</div></div>