<div>                Dear Prof. P. Blaha<br><br>I apologies for my basic knowledge, I am new to this field and I accept that I have many things to learn. I try my best to learn all these correctly.<br><br>I checked the DOS in my scf-sp and hubbard-U added case.<br><br> In runsp_lapw, it is metallic and conduction bands are dominated by f-electrons only, and it doesn't show any f-electron in valance band.  the QTL001and QTL002 are;<br><br>:PCS001: PARTIAL CHARGES SPHERE =  1 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ<br>:QTL001: 0.9965 2.0077 5.1500 0.5380 0.7525 0.7120 0.7040 1.0287 1.1474 1.0191 1.0006 1.0300<br><br>:PCS002: PARTIAL CHARGES SPHERE =  2 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ<br>:QTL002: 0.9876 2.0052 5.1575 0.5790 0.7022 0.7209 0.7028 1.0530 1.0304 1.0279 1.0188 1.0203<br><br><br>in runsp_lapw -orb case:<br>Still metal. But there are f-bands in both valance and conduction bands now<br><br>:PCS001: PARTIAL CHARGES SPHERE =  1 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ<br>:QTL001: 0.9077 2.13283 5.1590 0.3702 0.7082 0.7452 0.7073 1.1299 1.0491 1.0198 1.0193 1.0390<br><br>:PCS002: PARTIAL CHARGES SPHERE =  2 S,P,D,F,PX,PY,PZ,D-Z2,D-X2Y2,D-XY,D-XZ,D-YZ<br>:QTL002: 0.9154 2.1210 5.1754 0.5125 0.7100 0.7008 0.7116 1.1582 1.0402 1.0320 1.0218 1.0239<br><br><br><br>In both runsp_lapw and orb cases, I couldn't see d-orbital contribution in DOS.<br><br>Thank you<br><br>Prasad<br><br><br><br><br><br><br>            </div>            <div class="yahoo_quoted" style="margin:10px 0px 0px 0.8ex;border-left:1px solid #ccc;padding-left:1ex;">                        <div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">                                <div>                    On Thursday, August 1, 2019, 12:32:52 a.m. CST, Peter Blaha <pblaha@theochem.tuwien.ac.at> wrote:                </div>                <div><br></div>                <div><br></div>                <div><div dir="ltr">I guess you should first "understand" the results of a calculation, <br clear="none">before you start to experiment with more complicated tasks.<br clear="none"><br clear="none">After runsp_lapw  you saved the calculation.<br clear="none">What was the result of this calculation ??? Was it an insulator ??<br clear="none">Did you plot a DOS (also partial DOS including U-f states ?<br clear="none">Where are the 5f states ?<br clear="none">What was the occupation of the U-5f states (:QTL0XX line in the saved <br clear="none">scf file).<br clear="none"><br clear="none">If even a GGA calculation gives "zero" (or very small) 5f occupation, a <br clear="none">GGA+U calculation can ONLY give you even smaller occupation, since a <br clear="none">filling of less than 0.5 e/orbital will shift these states further up <br clear="none">with GGA+U.<br clear="none">How much did the gap change with U ??<br clear="none">How are the occupations ??<br clear="none"><br clear="none">I don't think that at least for your U6+ there is anything to play.<br clear="none"><br clear="none">For the U+5 it could be different, but again, first one would need more <br clear="none">info in order to be able to guide you.<br clear="none"><div class="yqt8684475281" id="yqtfd69260"><br clear="none">Am 31.07.2019 um 18:27 schrieb prasad jayasena:<br clear="none">> Dear experts<br clear="none">> <br clear="none">> I am trying a test case to learn the correct procedure to check local minima.<br clear="none">> Dear experts<br clear="none">> <br clear="none">> I am sending this followup questions related to my previous post on " electron occupancy in dmatup/dn" in the mailing list.<br clear="none">> <br clear="none">> <br clear="none">> In the firt step I finished a "runsp_lapw" and saved. I checked case.dmatup/dn and they were empty.<br clear="none">> Then "runsp_lapw -orb -dm" is employed.<br clear="none">> <br clear="none">> the dmatup/dn both have values similar to the one in my previous email. I believe this work flow is corresponding to a random localization.<br clear="none">> <br clear="none">> I appreciate your support on following questions:<br clear="none">> <br clear="none">> (1)<br clear="none">> If I want to fix initial density matrix, is this the correct place to do it? Is the only file I have to edit case.dmatup/dn?<br clear="none">> <br clear="none">> Or,  do I have to do that without trying a random PBE+U as I have done? But in that case since the case.dmatup/dn from PBE is empty, I wonder how to create that file correctly.<br clear="none">> <br clear="none">> (2)<br clear="none">> If I need to include Spin orbital coupling, how I should proceed after this?<br clear="none">> <br clear="none">> (3)<br clear="none">> In my case I expected to have d-orbital electrons. But I do not see them in the case.dmatup/dn. Are there any special flow to find d-matrix in the output?<br clear="none">> <br clear="none">> Thank you in advance<br clear="none">> <br clear="none">> Prasad</div><br clear="none">> _______________________________________________<br clear="none">> Wien mailing list<br clear="none">> <a shape="rect" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br clear="none">> <a shape="rect" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br clear="none">> SEARCH the MAILING-LIST at:  <a shape="rect" href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br clear="none">> <br clear="none"><br clear="none">-- <br clear="none">--------------------------------------------------------------------------<br clear="none">Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br clear="none">Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br clear="none">Email: <a shape="rect" ymailto="mailto:blaha@theochem.tuwien.ac.at" href="mailto:blaha@theochem.tuwien.ac.at">blaha@theochem.tuwien.ac.at</a>    WIEN2k: <a shape="rect" href="http://www.wien2k.at" target="_blank">http://www.wien2k.at</a><br clear="none">WWW: <br clear="none"><a shape="rect" href="http://www.imc.tuwien.ac.at/tc_blaha------------------------------------------------------------------------- " target="_blank">http://www.imc.tuwien.ac.at/tc_blaha------------------------------------------------------------------------- </a><br clear="none"><br clear="none">_______________________________________________<br clear="none">Wien mailing list<br clear="none"><a shape="rect" ymailto="mailto:Wien@zeus.theochem.tuwien.ac.at" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br clear="none"><a shape="rect" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br clear="none">SEARCH the MAILING-LIST at:  <a shape="rect" href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><div class="yqt8684475281" id="yqtfd49814"><br clear="none"></div></div></div>            </div>                </div>