<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><div>Dear Prof. Peter Blaha,</div><div><br></div><div>I think I found another issue with kgen -1 mode especially concerned with I4/mmm symmetry.</div><div><br></div><div>As previously mentioned, the '-1' mode of kgen generates uneven k divisions, as you rightly pointed out. However, I've observed potential issues with the weight of the uneven k mesh, particularly in the context of I4/mmm symmetry (and possibly other symmetries as well). This discrepancy appears to lead to inaccuracies in determining the Fermi energy.</div><div><br></div><div>Here is how I benchmarked:</div><div><br></div><div>1. run_lapw -p  till converge for CaFe2As2 and 10 10 10 shift kmesh. save_lapw -d scf</div><div>2. x kgen generate 9 9 9 non shifted mesh | x lapw1 -p; x lapw2 -p -fermi | get Ef=0.5076760931</div><div>3. x kgen generate 9 9 13 non shifted mesh | x lapw1 -p; x lapw2 -p -fermi | get Ef=0.5028504280</div><div>2. x kgen -fbz generate full 9 9 9 non shifted mesh | x lapw1 -p; x lapw2 -p -fermi | get Ef=0.5076760047</div><div>3. x kgen -fbz generate full 9 9 13 non shifted mesh | x lapw1 -p; x lapw2 -p -fermi | get Ef=0.5077417577</div><div><br></div><div>You can observe the significant difference in the Ef obtained from the 9 9 13 kmesh in the irreducible Brillouin zone (IBZ). However, the Ef calculated from the full Brillouin zone (FBZ) using the same mesh appears normal. Consequently, I suspect there may be an issue with the weight of k-points in the 9 9 13 mesh generated by your modified kgen. On the other hand, equal division kmesh 9 9 9 gives consistent Ef for both ibz and fbz.</div><div><br></div><div>What is more, for 9 9 13 IBZ kmesh, when viewing fermi surface. xcrysden 1.6.2 gives error. xcrysden 1.5.60 gives completely wrong fermi surface. </div><div><br></div><div>Please help, thank you very much!</div><div><br></div><div>best regards</div><div><br></div><div><hr align="left" style="margin: 0 0 10px 0;border: 0;border-bottom:1px solid #E4E5E6;height:0;line-height:0;font-size:0;padding: 20px 0 0 0;width: 50px;"><div style="font-size:14px;font-family:Verdana;color:#000;"><p></p></div></div><div> </div><div style="position: relative;"><div><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"                                                                                    <peter.blaha@tuwien.ac.at>;</div><div><b>Date:</b> Mon, Mar 11, 2024 08:33 PM</div><div><b>To:</b> "wien"<wien@zeus.theochem.tuwien.ac.at>;<wbr></div><div></div><div><b>Subject:</b> Re: [Wien] Inconsistency in kgen</div></div><div><br></div>Ups.<br>Here it comes.<br><br>Am 11.03.2024 um 13:26 schrieb balabi via Wien:<br>> Dear Prof. Peter Blaha<br>> <br>> Thank you so much for your reply! But I can not find your attachment of <br>> bravai.f.gz<br>> <br>> best regards<br>> <br>> ------------------------------------------------------------------------<br>> <br>> <br>> <br>> ------------------ Original ------------------<br>> *From:* "A Mailing list for WIEN2k users" <peter.blaha@tuwien.ac.at>;<br>> *Date:* Mon, Mar 11, 2024 04:03 PM<br>> *To:* "wien"<wien@zeus.theochem.tuwien.ac.at>;<br>> *Subject:* Re: [Wien] Inconsistency in kgen<br>> <br>> Hi,<br>> <br>> Thank you very much for your report. I can confirm the problem.<br>> <br>> Both, for bct and bco lattices (body-centered tetragonal or <br>> orthorhombic) kgen enforced in default modes equal divisions of the <br>> reciprocal lattice vectors (as it should be for bcc). This was not a <br>> good choice and the selection made by option "-1" (mesh density in <br>> bohr^-1) is correct.<br>> <br>> I attach a modified   bravai.f.gz  file, which should be copied and <br>> unziped in SRC_kgen, then recompiled (make; cp kgen ..).<br>> <br>> PS: The prevous setting was not a problem if your k-mesh is converged <br>> (besides the larger computational effort), but may lead to extra <br>> inaccuracy for non-converged meshes.<br>> <br>> Peter Blaha<br>> <br>> Am 10.03.2024 um 14:48 schrieb °ÍÀ­±È via Wien:<br>>> Dear wien2k developers and users,<br>>><br>>> I am using wien2k 23.2 and working with CaFe2As2 structure which has <br>>> I4/mmm symmetry. I am trying to generate klist using kgen. The kgen <br>>> has several mode:<br>>><br>>> the 1st mode is to specify k-mesh density<br>>><br>>> /NUMBER OF K-POINTS IN WHOLE CELL: (0 for 3 divisions of K, -1 for <br>>> delta-K)/<br>>> /-1/<br>>> / length of reciprocal lattice vectors (bohr^-1):  0.898   0.898   1.203/<br>>> /  Specify density of k-mesh in bohr^-1:/<br>>> /0.2/<br>>> /  Shift of k-mesh allowed. Do you want to shift: (0=no, 1=shift)/<br>>> /0/<br>>> /          17  k-points generated, ndiv=           5      5           7/<br>>> / delta-K (bohr^-1):     0.1796    0.1796    0.1719/<br>>> /KGEN ENDS/<br>>> /0.004u 0.016s 1:03.31 0.0%      0+0k 0+88io 0pf+0w/<br>>><br>>> As you can see, the ndiv=5 5 7<br>>><br>>> The 2nd mode is to specify number of k points<br>>><br>>> /NUMBER OF K-POINTS IN WHOLE CELL: (0 for 3 divisions of K, -1 for <br>>> delta-K)/<br>>> /1000/<br>>> / length of reciprocal lattice vectors (bohr^-1):  0.898   0.898   1.203/<br>>> /  Shift of k-mesh allowed. Do you want to shift: (0=no, 1=shift)/<br>>> /0/<br>>> /         102  k-points generated, ndiv=          10     10          10/<br>>> / delta-K (bohr^-1):     0.0898    0.0898    0.1203/<br>>> /KGEN ENDS/<br>>> /0.026u 0.003s 0:09.41 0.2%      0+0k 0+344io 0pf+0w/<br>>><br>>> as you can see unlike -1 mode, the ndiv=10 10 10 which is even.<br>>><br>>> The 3rd mode is to specify ndiv explicitly, and here comes the problem<br>>><br>>> /NUMBER OF K-POINTS IN WHOLE CELL: (0 for 3 divisions of K, -1 for <br>>> delta-K)/<br>>> /0/<br>>> / length of reciprocal lattice vectors (bohr^-1):  0.898   0.898   1.203/<br>>> /  Specify 3 mesh-divisions (n1,n2,n3):/<br>>> /5,5,7/<br>>> / Lattice symmetry requires equal mesh in x and z direction/<br>>> /  Specify 3 mesh-divisions (n1,n2,n3):/<br>>><br>>> If I set ndiv as 5 5 7 as -1 mode. It just does not work. kgen insists <br>>> equal mesh on x and z. So I must input even ndiv like 5 5 5.<br>>><br>>> So the current issue is why the ndiv obtained with the -1 mode in kgen <br>>> is not uniform for the CaFe2As2 system, whereas ndiv obtained with the <br>>> k-point mode is uniform, and using the 0 mode, it's impossible to set <br>>> non-uniform ndiv at all£¿<br>>><br>>> finally, I am attaching the struct file of CaFe2As2 as below<br>>><br>>> /blebleble /<br>>> /B   LATTICE,NONEQUIV.ATOMS:  3 139 I4/mmm /<br>>> /             RELA /<br>>> /  7.383336  7.383336 21.922849 90.000000 90.000000 90.000000/<br>>> /ATOM  -1: X=0.00000000 Y=0.00000000 Z=0.00000000/<br>>> /          MULT= 1          ISPLIT=-2/<br>>> /Ca1        NPT=  781  R0=0.00005000 RMT= 2.50000     Z: 20.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>>> /ATOM  -2: X=0.50000000 Y=0.00000000 Z=0.25000000/<br>>> /          MULT= 2          ISPLIT=-2/<br>>> /      -2: X=0.00000000 Y=0.50000000 Z=0.25000000/<br>>> /Fe1        NPT=  781  R0=0.00005000 RMT= 2.28        Z: 26.0 /<br>>> /LOCAL ROT MATRIX:    0.7071068-0.7071068 0.0000000/<br>>> /                     0.7071068 0.7071068 0.0000000/<br>>> /                     0.0000000 0.0000000 1.0000000/<br>>> /ATOM  -3: X=0.00000000 Y=0.00000000 Z=0.63490625/<br>>> /          MULT= 2          ISPLIT=-2/<br>>> /      -3: X=0.00000000 Y=0.00000000 Z=0.36509375/<br>>> /As1        NPT=  781  R0=0.00005000 RMT= 2.17        Z: 33.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>>> /  16      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>>> / 0-1 0 0.00000000/<br>>> / 1 0 0 0.00000000/<br>>> / 0 0 1 0.00000000/<br>>> /       3/<br>>> / 0 1 0 0.00000000/<br>>> /-1 0 0 0.00000000/<br>>> / 0 0 1 0.00000000/<br>>> /       4/<br>>> /-1 0 0 0.00000000/<br>>> / 0 1 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /       5/<br>>> / 1 0 0 0.00000000/<br>>> / 0-1 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /       6/<br>>> / 0 1 0 0.00000000/<br>>> / 1 0 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /       7/<br>>> / 0-1 0 0.00000000/<br>>> /-1 0 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /       8/<br>>> /-1 0 0 0.00000000/<br>>> / 0-1 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /       9/<br>>> / 1 0 0 0.00000000/<br>>> / 0 1 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /      10/<br>>> / 0 1 0 0.00000000/<br>>> /-1 0 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /      11/<br>>> / 0-1 0 0.00000000/<br>>> / 1 0 0 0.00000000/<br>>> / 0 0-1 0.00000000/<br>>> /      12/<br>>> / 1 0 0 0.00000000/<br>>> / 0-1 0 0.00000000/<br>>> / 0 0 1 0.00000000/<br>>> /      13/<br>>> /-1 0 0 0.00000000/<br>>> / 0 1 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>>><br>>><br>>> best regards<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>> Peter Blaha,  Inst. f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>> Phone: +43-158801165300<br>> Email:peter.blaha@tuwien.ac.at           <br>> WWW:http://www.imc.tuwien.ac.at       WIEN2k:http://www.wien2k.at<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>Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>Phone: +43-1-58801-165300<br>Email: peter.blaha@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><br></div>