Thank you Dr. Blaha<div>This work well.</div><div>All the best for you.<br><br><div id="ymail_android_signature"><a id="ymail_android_signature_link" href="https://mail.onelink.me/107872968?pid=nativeplacement&c=Global_Acquisition_YMktg_315_Internal_EmailSignature&af_sub1=Acquisition&af_sub2=Global_YMktg&af_sub3=&af_sub4=100000604&af_sub5=EmailSignature__Static_">Sent from Yahoo Mail on Android</a></div> <br> <blockquote style="margin: 0 0 20px 0;"> <div style="font-family:Roboto, sans-serif; color:#6D00F6;"> <div>On Wed, Nov 15, 2023 at 0:32, Peter Blaha</div><div><peter.blaha@tuwien.ac.at> wrote:</div> </div> <div style="padding: 10px 0 0 20px; margin: 10px 0 0 0; border-left: 1px solid #6D00F6;"> Well, at the end it is exactly as I said:<br clear="none"><br clear="none">Your manual RMT settings are very bad. One of the rules is, that <br clear="none">smallest and largest RMTs must not be too different.<br clear="none">With your spheres you get "effective" R(Si)Kmax of more than 12 and this <br clear="none">gives numerical linear dependency.<br clear="none"><br clear="none">With your struct file: RMT(N)=1.2  (Si)=2.1<br clear="none">init<br clear="none">set RKmax=7<br clear="none">run<br clear="none">----> gives the semicore error you describe due to some ghostbands.<br clear="none"><br clear="none">init<br clear="none">run<br clear="none">save  rkm6<br clear="none">set RKmax=7<br clear="none">run<br clear="none">---> runs through and converges.<br clear="none"><br clear="none">However, increasing RKmax to 8 gives ghost bands again.<br clear="none"><br clear="none">---------------------------------<br clear="none"><br clear="none">setrmt<br clear="none">cp case.struct_setrmt case.struct  ( similar RMTs for Si and N, around 1.6)<br clear="none">init -ecut -8   (to avoid core leakage)<br clear="none">set rkmax=7<br clear="none">run<br clear="none">---> converges without problems. One can increase Rkmax further to 8 and <br clear="none">even 9.<br clear="none"><br clear="none">----------------------------------<br clear="none"><br clear="none">PS: These ghostbands are located in the interstital, give no qtl-b <br clear="none">errors. Once such a state is taken into the density, you get these <br clear="none">"select"-errors.<br clear="none"><br clear="none">Am 14.11.2023 um 19:19 schrieb hajar.nejatipoor--- via Wien:<br clear="none">> Dr. Blaha<br clear="none">> sometimes, *semicore* error appears in iteration3, sometime in 6, and <br clear="none">> ... (with changing rmts).<br clear="none">> <br clear="none">> I tried with the struct attached here and the default init_lapw. After <br clear="none">> finalized initialization, I changed RKm=7 (and Emax=3) in case.in1c <br clear="none">> attached here, and run dstart. This time, semicore error was appeared in <br clear="none">> the *first* iteration!<br clear="none">> I have been confused why?<br clear="none">> (I see that default init_lapw for above setting (RKm=7 and Emax=3)) <br clear="none">> contains 24 kpoints. and error is appeared.<br clear="none">> (In another scf with above struct (RKm=7 and Emax=1.5), default number <br clear="none">> of kpoints was 7! ). and scf ended successfully.<br clear="none">> <br clear="none">> May this error be dependent on the number of k-point and the number of <br clear="none">> cores in .machines file?<br clear="none">> .machines file in two calculations contained 7 cores.<br clear="none">> On Tuesday, November 14, 2023 at 03:08:50 PM GMT+3:30, Peter Blaha <br clear="none">> <<a shape="rect" ymailto="mailto:peter.blaha@tuwien.ac.at" href="mailto:peter.blaha@tuwien.ac.at">peter.blaha@tuwien.ac.at</a>> wrote:<br clear="none">> <br clear="none">> <br clear="none">> Again, your message gets too big. You must delete the older content.<br clear="none">> -------------------------------------------------------<br clear="none">> <br clear="none">> grep :DIS in case.scf:<br clear="none">> <br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0122755 for atom    3 spin 1)<br clear="none">> 0.0083335<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0117894 for atom    3 spin 1)<br clear="none">> 0.0077543<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0405700 for atom    1 spin 1)<br clear="none">> 0.0200036<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.2010006 for atom    1 spin 1)<br clear="none">> 0.0741310<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0164221 for atom    1 spin 1)<br clear="none">> 0.0107305<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.1052329 for atom    1 spin 1)<br clear="none">> 0.0370176<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0075476 for atom    1 spin 1)<br clear="none">> 0.0021153<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0848258 for atom    1 spin 1)<br clear="none">> 0.0300654<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0018758 for atom    1 spin 1)<br clear="none">> 0.0007564<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0008796 for atom    3 spin 1)<br clear="none">> 0.0006306<br clear="none">> :DIS  :  CHARGE DISTANCE      (  0.0013281 for atom    3 spin 1)<br clear="none">> 0.0008331<br clear="none">> <br clear="none">> after iteration 3, the semicore error is appeared for rkm=7.<br clear="none">> ------------------------------------------------------------------<br clear="none">> <br clear="none">> Nobody knows what you were sending. Is this from the RKM=6 calculation ?<br clear="none">> <br clear="none">> You have done a scf with RKmax=6.  This should be saved.<br clear="none">> <br clear="none">> Then you have an empty scffile, and you should increase RKmax to 7 and<br clear="none">> run the scf.<br clear="none">> If the error occurs after iteration 3, we expect to see exactly 3 lines.<br clear="none">> <br clear="none">> ???<br clear="none">> Did you ever save the rkm6 results ?<br clear="none">> <br clear="none">> restore them in a new directory, increase rkmax and then run_lapw.<br clear="none">> When it crashes, show us :dis.<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<br clear="none">> Email: <a shape="rect" ymailto="mailto:peter.blaha@tuwien.ac.at" href="mailto:peter.blaha@tuwien.ac.at">peter.blaha@tuwien.ac.at</a> <mailto:peter.blaha@tuwien.ac.at>    <br clear="none">> WIEN2k: <a shape="rect" href="http://www.wien2k.at" target="_blank">http://www.wien2k.at</a> <<a shape="rect" href="http://www.wien2k.at" target="_blank">http://www.wien2k.at</a>><br clear="none">> WWW: <a shape="rect" href="http://www.imc.tuwien.ac.at" target="_blank">http://www.imc.tuwien.ac.at</a> <<a shape="rect" href="http://www.imc.tuwien.ac.at" target="_blank">http://www.imc.tuwien.ac.at</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> <mailto:Wien@zeus.theochem.tuwien.ac.at><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">> <<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: <br clear="none">> <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">> <<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">> 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="yqt4672033803" id="yqtfd01928"><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<br clear="none">Email: <a shape="rect" ymailto="mailto:peter.blaha@tuwien.ac.at" href="mailto:peter.blaha@tuwien.ac.at">peter.blaha@tuwien.ac.at</a>    WIEN2k: <a shape="rect" href="http://www.wien2k.at" target="_blank">http://www.wien2k.at</a><br clear="none">WWW:   <a shape="rect" href="http://www.imc.tuwien.ac.at" target="_blank">http://www.imc.tuwien.ac.at</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><br clear="none"></div> </div> </blockquote></div>