<div dir="ltr"><div>Hi, </div>Thanks for all and sorry for the late follow up. <div>Yet I have not been successful in doing this after adopting the suggestions. </div><div><br></div><div>Take the SrVO3 as the example. </div><div><br></div><div>I first did "setrmt -r 3". The screen output is attached. </div><div>It created the SrVO3.struct_setrmt with RMT = 1.84 / 2.47 / 1.66 for V / Sr / O. </div><div><br></div><div>Then I did "init_lapw -prec 2". The screen output is also attached. </div><div>It resets the RMT of Sr from 2.47 to 2.35. </div><div>The chunk of lstart-related output is like this: </div><div><br></div><div>next is lstart<br> 3 Atoms found:  with labels V    Sr   O<br>generate atomic configuration for atom 1 : V<br>generate atomic configuration for atom 2 : Sr<br>generate atomic configuration for atom 3 : O<br>  SELECT XCPOT:<br>  recommended: PBE    [(13) GGA of Perdew-Burke-Ernzerhof 96]<br>               LDA    [( 5)]<br>               WC     [(11)  GGA of Wu-Cohen 2006]<br>               PBESOL [(19) GGA of Perdew etal. 2008]<br>  SELECT ENERGY to separate core and valence states:<br>  recommended: -6.0 Ry (check how much core charge leaks out of MT-sphere)<br>  ALTERNATIVELY: specify charge localization (between 0.97 and 1.0) to select core state<br>STOP LSTART ENDS</div><div><br></div><div>Here I can't find the "charge leakage" information in the screen output of lstart. </div><div>It suggests "(check how much core charge leaks out of MT-sphere)", but where should I look for it? </div><div><br></div><div>Then, doing "run_lapw" still results in the same error. </div><div>STOP  LAPW0 END<br>STOP  LAPW1 END<br>STOP  LAPW2 END<br>STOP  CORE  END<br>STOP ERROR: NEC01 charge leakage too large<br></div><div><br></div><div>Please suggest me to do more about this! </div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Nov 12, 2024 at 11:43 AM Peter Blaha <<a href="mailto:peter.blaha@tuwien.ac.at">peter.blaha@tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi,<br>
<br>
You have to learn a few wien2k "rules".<br>
<br>
i) How did you create the struct files ?  Always set the sphere sizes <br>
according to the recommendations (see below).<br>
<br>
ii) It is VERY important that you use     setrmt  to set proper RMT <br>
radii. If you want to change lattice parameters or optimize positions, <br>
you may want to reduce the RMT  using   setrmt -r 3    or similar.<br>
Remember: in order to compare total energies, RMTs (and RKMAX) should be <br>
identical (or "equivalent").<br>
<br>
iii) During initialization observe the output of lstart, which lists all <br>
orbitals of all atoms, their energies and in particular the "charge <br>
leakage". If it is too large, it will create a -lcore file for you, or <br>
with -prec 2 will rerun lstart with a smaller -ecut automatically to <br>
avoid any core leakage.<br>
<br>
iv) The init_lapw   has now these   -prec 0-3/0-3n (for semiconductors) <br>
switches and it is highly recommended to use them (and not set emin, <br>
rkmax, k-points manually - except you "know what you are doing").<br>
For such small cases I would always use -prec 2 (2n), larger one can <br>
start with -prec 1. When you do large force minimizations -prec 0 might <br>
be a good choice (at the beginning) until the forces are small.<br>
<br>
-----------------<br>
The error occurred in MIXER, thus you should not only provide the "stop <br>
Error" statement, but check    *.error  files and *scfm, *outputm (m for <br>
mixer) and *dayfile.<br>
You may find more information which may give you a hint what was wrong.<br>
<br>
<br>
<br>
Am 12.11.2024 um 04:25 schrieb Long Zhang via Wien:<br>
> Dear Wien users,<br>
> <br>
> I am running WIEN2k 23.2 with the attached two simple cases.<br>
> The initialization runs OK in both cases.<br>
> init_lapw -vxc 13 -ecut -6.0 -rkmax 7 -numk 500<br>
> <br>
> Then I got the following error when executing run_lapw in both cases.<br>
> <br>
> STOP  LAPW0 END<br>
> STOP  LAPW1 END<br>
> STOP  LAPW2 END<br>
> STOP  CORE  END<br>
> STOP ERROR: NEC01 charge leakage too large<br>
> <br>
>  >   stop error<br>
> <br>
> I searched the mail archive and there's a similar question in 2023<br>
> (<a href="https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/" rel="noreferrer" target="_blank">https://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/</a> <br>
> msg22819.html <<a href="https://www.mail-archive.com/" rel="noreferrer" target="_blank">https://www.mail-archive.com/</a> <br>
> <a href="http://wien@zeus.theochem.tuwien.ac.at/msg22819.html" rel="noreferrer" target="_blank">wien@zeus.theochem.tuwien.ac.at/msg22819.html</a>>)<br>
> where it's suggested to re-run the initialization with only "-prec 1" <br>
> with no other options.<br>
> I did it, for both cases, but it's not helping.<br>
> What else can I do to fix the too large charge leakage?<br>
> <br>
> Best regards,<br>
> Long<br>
> <br>
> <br>
> _______________________________________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
<br>
-- <br>
-----------------------------------------------------------------------<br>
Peter Blaha,  Inst. f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-158801165300<br>
Email: <a href="mailto:peter.blaha@tuwien.ac.at" target="_blank">peter.blaha@tuwien.ac.at</a><br>
WWW:   <a href="http://www.imc.tuwien.ac.at" rel="noreferrer" target="_blank">http://www.imc.tuwien.ac.at</a>      WIEN2k: <a href="http://www.wien2k.at" rel="noreferrer" target="_blank">http://www.wien2k.at</a><br>
-------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>