
Dear Prof. P. Blaha<br><div><div><br></div><div>Thank you so much for your detailed answer. I tested CuOH according to your suggestion and the bands are fine.</div><div><br></div><div>best regards<br></div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ Original ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>From:</b>                                                                                                                        "A Mailing list for WIEN2k users"                                                                                    <pblaha@theochem.tuwien.ac.at>;</div><div><b>Date:</b> Sat, Jan 23, 2021 04:13 PM</div><div><b>To:</b> "wien"<wien@zeus.theochem.tuwien.ac.at>;<wbr></div><div></div><div><b>Subject:</b> Re: [Wien] [EXTERNAL] Bands are often not smooth when zoom in</div></div><div><br></div>If your basis set is not converged, it is clear that the eigenvalues <br>will change (go down in energy) when the basis set is increased.<br><br>The cutoff is defined by RKMAX, i.e. we take all plane waves such that <br>|(k+K)| < KMAX.  In a bandstructure plot k varies slightly, but it can <br>then happen that at a certain k, the maximal K is too large and gets <br>reduced. Then a "jump" in the bandstructure will appear.<br>In this description, "k" is the vector of the k-point, "K" is a full <br>reciprocal lattice vector.<br><br>PS: Yes, RKMAX=7 is of course too small for a 3d (Cu) element, while <br>RKMAX=7 is much too large for H. It is a compromise.<br><br>On our website www.wien2k.at at reg.user and faq there are some <br>suggestions for a "minimal" RKMAX for each element. You look into your <br>structure and identify the atom which has the smallest RMT. This atom <br>determines RKMAX according to the list given on this site. However, <br>these values are kind of "minimal" RKMAX values (quite good for eg. <br>optimization of atomic positions (-min), but for highest precision they <br>should be increased by 1-3.<br><br>Just a simple example:<br>You do Cu and need RKMAX=9 for full convergence.<br>You do CuOH with RMT=1.8; 1.2 and 0.6 for the 3 elements.<br>RKMAX=3  is counting for H; gives "effective RKMAX=6 for O and 9 for Cu.<br>Hence RKMAX=3 is fully converged for CuOH (when using RMTs according to <br>the suggestion of setrmt).<br><br><br>Am 23.01.2021 um 04:49 schrieb 것윗궐:<br>> Dear professor Jianxin Zhu,<br>> <br>> Thank you so much for your quick reply.<br>> <br>> Enlarging RKmax solves the problem! I found I have to set RKmax=9 which <br>> is pretty large compared to default 7 to completely eliminate all <br>> sawtooth artifacts in bands.<br>> <br>> But I still want to know why a smaller RKmax would cause eigenvalues to <br>> be discontinuous (in essence hamiltonian is discontinuous) at several k <br>> points.<br>> <br>> best regards<br>> <br>> <br>> <br>> ------------------ Original ------------------<br>> *From:* "A Mailing list for WIEN2k users" <jxzhu@lanl.gov>;<br>> *Date:* Fri, Jan 22, 2021 11:57 PM<br>> *To:* "A Mailing list for WIEN2k users"<wien@zeus.theochem.tuwien.ac.at>;<br>> *Subject:* Re: [Wien] [EXTERNAL] Bands are often not smooth when zoom in<br>> <br>> What rkmax value are you using? Usually a larger rkmax helps make the <br>> bands smooth.<br>> <br>> Jianxin<br>> <br>> <br>> Sent with BlackBerry Work<br>> (www.blackberry.com)<br>> <br>> <br>> *From: *Wien <wien-bounces@zeus.theochem.tuwien.ac.at <br>> <mailto:wien-bounces@zeus.theochem.tuwien.ac.at>> on behalf of: 것윗궐 <br>> <balabi@qq.com <mailto:balabi@qq.com>><br>> *Date: *Friday, Jan 22, 2021, 8:15 AM<br>> *To: *wien <wien@zeus.theochem.tuwien.ac.at <br>> <mailto:wien@zeus.theochem.tuwien.ac.at>><br>> *Subject: *[EXTERNAL] [Wien] Bands are often not smooth when zoom in<br>> <br>> Dear wien2k developers,<br>> <br>> I found bands obtained from wien2k are often not smooth and have many <br>> small sawtooth bumps especially when you zoom in. I can reproduce this <br>> kind of band artifacts in version 19.2 and 17.1, on differnet computers <br>> and link to different version of MKL. I'll take simple copper as an <br>> example to show how I obtain the bands:<br>> <br>> the structure file I use is as:<br>> <br>> blebleble<br>> F   LATTICE,NONEQUIV.ATOMS:  1 225 Fm-3m<br>>               RELA<br>>    6.872726  6.872726  6.872726 90.000000 90.000000 90.000000<br>> ATOM   1: X=0.00000000 Y=0.00000000 Z=0.00000000<br>>            MULT= 1          ISPLIT= 2<br>> Cu1        NPT=  781  R0=0.00005000 RMT=    2.0000   Z: 29.0<br>> LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000<br>>                       0.0000000 1.0000000 0.0000000<br>>                       0.0000000 0.0000000 1.0000000<br>>    48      NUMBER OF SYMMETRY OPERATIONS<br>> 1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>> 0 0 1 0.00000000<br>>         1<br>> -1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>> 0 0 1 0.00000000<br>>         2<br>> -1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>> 0 0-1 0.00000000<br>>         3<br>> 1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>> 0 0-1 0.00000000<br>>         4<br>> 0 0 1 0.00000000<br>> 1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>>         5<br>> 0 0 1 0.00000000<br>> -1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>>         6<br>> 0 0-1 0.00000000<br>> -1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>>         7<br>> 0 0-1 0.00000000<br>> 1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>>         8<br>> 0 1 0 0.00000000<br>> 0 0 1 0.00000000<br>> 1 0 0 0.00000000<br>>         9<br>> 0-1 0 0.00000000<br>> 0 0 1 0.00000000<br>> -1 0 0 0.00000000<br>>        10<br>> 0 1 0 0.00000000<br>> 0 0-1 0.00000000<br>> -1 0 0 0.00000000<br>>        11<br>> 0-1 0 0.00000000<br>> 0 0-1 0.00000000<br>> 1 0 0 0.00000000<br>>        12<br>> 0 1 0 0.00000000<br>> 1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>>        13<br>> 0-1 0 0.00000000<br>> -1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>>        14<br>> 0 1 0 0.00000000<br>> -1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>>        15<br>> 0-1 0 0.00000000<br>> 1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>>        16<br>> 1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>> 0-1 0 0.00000000<br>>        17<br>> -1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>> 0 1 0 0.00000000<br>>        18<br>> -1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>> 0-1 0 0.00000000<br>>        19<br>> 1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>> 0 1 0 0.00000000<br>>        20<br>> 0 0 1 0.00000000<br>> 0 1 0 0.00000000<br>> -1 0 0 0.00000000<br>>        21<br>> 0 0 1 0.00000000<br>> 0-1 0 0.00000000<br>> 1 0 0 0.00000000<br>>        22<br>> 0 0-1 0.00000000<br>> 0 1 0 0.00000000<br>> 1 0 0 0.00000000<br>>        23<br>> 0 0-1 0.00000000<br>> 0-1 0 0.00000000<br>> -1 0 0 0.00000000<br>>        24<br>> -1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>> 0 0-1 0.00000000<br>>        25<br>> 1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>> 0 0-1 0.00000000<br>>        26<br>> 1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>> 0 0 1 0.00000000<br>>        27<br>> -1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>> 0 0 1 0.00000000<br>>        28<br>> 0 0-1 0.00000000<br>> -1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>>        29<br>> 0 0-1 0.00000000<br>> 1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>>        30<br>> 0 0 1 0.00000000<br>> 1 0 0 0.00000000<br>> 0-1 0 0.00000000<br>>        31<br>> 0 0 1 0.00000000<br>> -1 0 0 0.00000000<br>> 0 1 0 0.00000000<br>>        32<br>> 0-1 0 0.00000000<br>> 0 0-1 0.00000000<br>> -1 0 0 0.00000000<br>>        33<br>> 0 1 0 0.00000000<br>> 0 0-1 0.00000000<br>> 1 0 0 0.00000000<br>>        34<br>> 0-1 0 0.00000000<br>> 0 0 1 0.00000000<br>> 1 0 0 0.00000000<br>>        35<br>> 0 1 0 0.00000000<br>> 0 0 1 0.00000000<br>> -1 0 0 0.00000000<br>>        36<br>> 0-1 0 0.00000000<br>> -1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>>        37<br>> 0 1 0 0.00000000<br>> 1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>>        38<br>> 0-1 0 0.00000000<br>> 1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>>        39<br>> 0 1 0 0.00000000<br>> -1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>>        40<br>> -1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>> 0 1 0 0.00000000<br>>        41<br>> 1 0 0 0.00000000<br>> 0 0-1 0.00000000<br>> 0-1 0 0.00000000<br>>        42<br>> 1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>> 0 1 0 0.00000000<br>>        43<br>> -1 0 0 0.00000000<br>> 0 0 1 0.00000000<br>> 0-1 0 0.00000000<br>>        44<br>> 0 0-1 0.00000000<br>> 0-1 0 0.00000000<br>> 1 0 0 0.00000000<br>>        45<br>> 0 0-1 0.00000000<br>> 0 1 0 0.00000000<br>> -1 0 0 0.00000000<br>>        46<br>> 0 0 1 0.00000000<br>> 0-1 0 0.00000000<br>> -1 0 0 0.00000000<br>>        47<br>> 0 0 1 0.00000000<br>> 0 1 0 0.00000000<br>> 1 0 0 0.00000000<br>>        48<br>> <br>> then I do step by step<br>> 1. init_lapw -b<br>> 2. run_lapw -in1new 3   # I use -in1new 3 to avoid the QTL-B value <br>> warning message in the last few cycles of case.scf file<br>> 3. using xcrysden to generate a k path along GAMMA-W-X-L with 400 k points.<br>>      since the klist_band file is too long. I share it on <br>> https://pastebin.com/ZaYuFMZE<br>> 4. x lapw1 -band<br>> 5. x lapw2 -band -qtl<br>> <br>> finally, I use the data in case.qtl to plot bands<br>> in [-10eV,10eV] range, it looks pretty normal as shared in below link<br>> https://pasteboard.co/JKPC1tB.png<br>> in [-2eV,-1eV] zoom in range, the look of bands becomes abnormal with <br>> full of sawtooth bump as shared in link<br>> https://pasteboard.co/JKPCB7Y.png<br>> <br>> So I am wondering if this is a bug? Or what I have done wrong to have <br>> this kind of bands(it appears in all my other bands calculation)? I <br>> inspected the k path and found it is pretty fine. So I can not think of <br>> a reason why the eigenvalue solutions would be discontinous at two <br>> nearby k points.<br>> <br>> best regards<br>> <br>> <br>> <br>> <br>> <br>> <br>> <br>> _______________________________________________<br>> Wien mailing list<br>> Wien@zeus.theochem.tuwien.ac.at<br>> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<br>> SEARCH the MAILING-LIST at:  http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html<br>> <br><br>-- <br>--------------------------------------------------------------------------<br>Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>Email: blaha@theochem.tuwien.ac.at    WIEN2k: http://www.wien2k.at<br>WWW:   http://www.imc.tuwien.ac.at<br>-------------------------------------------------------------------------<br>_______________________________________________<br>Wien mailing list<br>Wien@zeus.theochem.tuwien.ac.at<br>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<br>SEARCH the MAILING-LIST at:  http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html<br></div>