<div dir="ltr">Dear Wien2K community,<br><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  I am currently using the latest version of the WIEN2k package (<b>Version: WIEN2k_23.2</b>). I am calculating the <b>momentum matrix elements</b> using the OPTIC module as implemented in the WIEN2k package. From userguide I found that the momentum matrix elements will be written in <b><i>case.mommat2</i></b> file. </div><div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  I am performing SOC included calculations on TaAs using the <b><i>PBESol</i></b> as an exchange correlation functional. My <b>TaAs.inop</b> file looks like:</div><div><br></div><div>99999 1 Â  Â  Â  number of k-points, first k-point<br>-0.5 0.5 9999 Emin, Emax for matrix elements, NBvalMAX<br>2 Â  Â  Â  Â  Â  Â  number of choices (columns in *outmat): 2: hex or tetrag. case<br>1 Â  Â  Â  Â  Â  Â  Re xx<br>3 Â  Â  Â  Â  Â  Â  Re zz<br>ON Â  Â  Â  Â  Â  Â ON/OFF Â  writes MME to unit 4<br>1<br></div><div><br></div><div>I have observed that <b>once the scf converged ground state is obtained</b>, there is <b>no change in the momentum matrix elements values </b>(in TaAs.mommat2 file) <b>if we repeat the calculations</b> of TaAs.mommat2 file. <b>However, if I rerun the scf calculations again from starting and then obtain the TaAs.mommat2 file, the data in the file gets changed (considerable changes are observed).</b> Note: All the parameters during these two calculations are kept the same. I am showing below the output file (TaAs.mommat2 file) obtained during the two calculations.</div><div><br></div><div><b>First calculation:</b></div><div> </div><div>  Â Complex momentum matrix elements: (Re_x,Im_x),(Re_y,Im_y),(Re_z,Im_z)<br><br>  Â KP: Â  Â  1 NEMIN NEMAX : Â  Â 65 Â  70 dE: Â -0.5 Â  0.5 K: Â  Â  Â  Â  1<br><br>  65 Â  65 -6.791572E-14 Â 4.326127E-16 -1.908323E-12 -6.827043E-12 -1.230687E-12 Â 4.691060E-07 Â  Â 0.00000000<br>  65 Â  66 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00000000<br>  65 Â  67 -1.122490E-04 Â 1.272331E-04 Â 7.983179E-05 -3.217857E-05 Â 4.305635E-11 -9.848311E-12 Â  Â 0.17874679<br>  65 Â  68 Â 1.660215E-03 Â 5.409568E-04 Â 5.499541E-04 -1.663694E-03 Â 2.075423E-12 Â 3.826721E-12 Â  Â 0.17874679<br>  65 Â  69 Â 3.936429E-04 Â 1.240500E-03 Â 1.266557E-03 -3.814258E-04 Â 3.255843E-13 Â 4.740868E-12 Â  Â 0.51386663<br>  65 Â  70 -1.237442E-04 -3.543114E-04 -2.453771E-04 Â 1.574218E-04 Â 4.178337E-12 Â 6.320158E-12 Â  Â 0.51386663<br>  66 Â  66 -6.791572E-14 Â 4.326127E-16 -1.908323E-12 -6.827043E-12 -1.230687E-12 Â 4.691060E-07 Â  Â 0.00000000<br>  66 Â  67 Â 1.660215E-03 Â 5.409568E-04 Â 5.499541E-04 -1.663694E-03 Â 2.075423E-12 Â 3.826721E-12 Â  Â 0.17874679<br>  66 Â  68 -1.122490E-04 Â 1.272331E-04 Â 7.983179E-05 -3.217857E-05 Â 4.305635E-11 -9.848311E-12 Â  Â 0.17874679<br>  66 Â  69 -1.237442E-04 -3.543114E-04 -2.453771E-04 Â 1.574218E-04 Â 4.178337E-12 Â 6.320158E-12 Â  Â 0.51386663<br>  66 Â  70 Â 3.936429E-04 Â 1.240500E-03 Â 1.266557E-03 -3.814258E-04 Â 3.255843E-13 Â 4.740868E-12 Â  Â 0.51386663<br>  67 Â  67 -3.062120E-11 -1.361805E-15 -1.900231E-11 -1.083939E-11 -1.835530E-10 -1.298187E-05 Â  Â 0.00000000<br>  67 Â  68 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00000000<br>  67 Â  69 Â 8.836885E-03 -5.639363E-03 Â 5.375814E-03 Â 8.636869E-03 Â 4.880596E-02 Â 1.287101E-01 Â  Â 0.33511984<br>  67 Â  70 -1.671701E-03 Â 1.109599E-03 -2.253235E-03 -2.312311E-03 Â 1.325112E-01 Â 4.650952E-01 Â  Â 0.33511984<br>  68 Â  68 -3.062120E-11 -1.361805E-15 -1.900231E-11 -1.083939E-11 -1.835530E-10 -1.298187E-05 Â  Â 0.00000000<br>  68 Â  69 -1.671701E-03 Â 1.109599E-03 -2.253235E-03 -2.312311E-03 Â 1.325112E-01 Â 4.650952E-01 Â  Â 0.33511984<br>  68 Â  70 Â 8.836885E-03 -5.639363E-03 Â 5.375814E-03 Â 8.636869E-03 Â 4.880596E-02 Â 1.287101E-01 Â  Â 0.33511984<br>  69 Â  69 -7.367419E-12 Â 1.678033E-15 -1.891869E-11 Â 1.300217E-12 -4.739680E-11 Â 1.787176E-05 Â  Â 0.00000000<br>  69 Â  70 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00000000<br>  70 Â  70 -7.367419E-12 Â 1.678033E-15 -1.891869E-11 Â 1.300217E-12 -4.739680E-11 Â 1.787176E-05 Â  Â 0.00000000<br><br>  Â KP: Â  Â  2 NEMIN NEMAX : Â  Â 65 Â  72 dE: Â -0.5 Â  0.5 K: Â  Â  Â  Â  2<br><br>  65 Â  65 Â 9.561577E-02 Â 6.302006E-07 Â 9.561578E-02 -6.302075E-07 -3.861626E-10 Â 7.639556E-07 Â  Â 0.00000000<br>  65 Â  66 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00030716<br>  65 Â  67 -1.500039E-03 Â 1.544209E-03 -1.500158E-03 -1.544114E-03 Â 4.595806E-09 Â 5.467578E-04 Â  Â 0.16416515<br>  65 Â  68 -1.734816E-03 -2.038025E-03 -1.734780E-03 Â 2.038026E-03 Â 2.598629E-08 Â 1.625489E-04 Â  Â 0.17013061<br></div><div><br></div><div><b>Second calculation:</b></div><div><br></div><div>  Â Complex momentum matrix elements: (Re_x,Im_x),(Re_y,Im_y),(Re_z,Im_z)<br><br>  Â KP: Â  Â  1 NEMIN NEMAX : Â  Â 65 Â  70 dE: Â -0.5 Â  0.5 K: Â  Â  Â  Â  1<br><br>  65 Â  65 -7.845272E-12 Â 5.460652E-17 Â 6.845644E-13 -6.827471E-12 -3.414512E-12 Â 4.691060E-07 Â  Â 0.00000000<br>  65 Â  66 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00000000<br>  65 Â  67 -3.691215E-04 Â 1.629442E-03 Â 1.510884E-03 Â 3.903687E-04 -2.048646E-11 -2.391975E-12 Â  Â 0.17874679<br>  65 Â  68 -1.935785E-04 -5.696257E-04 Â 8.347352E-04 -1.460673E-04 Â 3.717129E-11 -1.078993E-11 Â  Â 0.17874679<br>  65 Â  69 Â 1.272132E-04 -1.157650E-03 -4.132527E-04 -3.733865E-04 -1.266540E-12 Â 4.856729E-13 Â  Â 0.51386663<br>  65 Â  70 Â 3.871718E-04 -4.810268E-04 -1.183538E-03 -1.633025E-04 Â 1.166324E-12 -1.468572E-12 Â  Â 0.51386663<br>  66 Â  66 -7.845272E-12 Â 5.460652E-17 Â 6.845644E-13 -6.827471E-12 -3.414512E-12 Â 4.691060E-07 Â  Â 0.00000000<br>  66 Â  67 -1.935785E-04 -5.696257E-04 Â 8.347352E-04 -1.460673E-04 Â 3.717129E-11 -1.078993E-11 Â  Â 0.17874679<br>  66 Â  68 -3.691215E-04 Â 1.629442E-03 Â 1.510884E-03 Â 3.903687E-04 -2.048646E-11 -2.391975E-12 Â  Â 0.17874679<br>  66 Â  69 Â 3.871718E-04 -4.810268E-04 -1.183538E-03 -1.633025E-04 Â 1.166324E-12 -1.468572E-12 Â  Â 0.51386663<br>  66 Â  70 Â 1.272132E-04 -1.157650E-03 -4.132527E-04 -3.733865E-04 -1.266540E-12 Â 4.856729E-13 Â  Â 0.51386663<br>  67 Â  67 -6.209961E-11 -5.437752E-15 Â 7.061251E-11 -1.084517E-11 -1.171486E-10 -1.298187E-05 Â  Â 0.00000000<br>  67 Â  68 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00000000<br>  67 Â  69 -5.944409E-03 -1.682884E-03 -1.372490E-03 Â 6.934560E-03 -3.516587E-01 -1.020410E-01 Â  Â 0.33511984<br>  67 Â  70 Â 6.567929E-03 Â 1.343172E-03 Â 1.661758E-03 -5.511497E-03 -3.165542E-01 -1.362602E-01 Â  Â 0.33511984<br>  68 Â  68 -6.209961E-11 -5.437752E-15 Â 7.061251E-11 -1.084517E-11 -1.171486E-10 -1.298187E-05 Â  Â 0.00000000<br>  68 Â  69 Â 6.567929E-03 Â 1.343172E-03 Â 1.661758E-03 -5.511497E-03 -3.165542E-01 -1.362602E-01 Â  Â 0.33511984<br>  68 Â  70 -5.944409E-03 -1.682884E-03 -1.372490E-03 Â 6.934560E-03 -3.516587E-01 -1.020410E-01 Â  Â 0.33511984<br>  69 Â  69 Â 2.024387E-11 Â 5.052165E-15 Â 2.031190E-11 Â 1.307417E-12 -2.125244E-11 Â 1.787176E-05 Â  Â 0.00000000<br>  69 Â  70 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00000000<br>  70 Â  70 Â 2.024387E-11 Â 5.052165E-15 Â 2.031190E-11 Â 1.307417E-12 -2.125244E-11 Â 1.787176E-05 Â  Â 0.00000000<br><br>  Â KP: Â  Â  2 NEMIN NEMAX : Â  Â 65 Â  72 dE: Â -0.5 Â  0.5 K: Â  Â  Â  Â  2<br><br>  65 Â  65 Â 9.561577E-02 Â 6.302003E-07 Â 9.561577E-02 -6.302073E-07 -4.163895E-10 Â 7.639556E-07 Â  Â 0.00000000<br>  65 Â  66 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â 0.000000E+00 Â  Â 0.00030716<br>  65 Â  67 Â 1.500038E-03 -1.544209E-03 Â 1.500158E-03 Â 1.544114E-03 -5.230746E-09 -5.467574E-04 Â  Â 0.16416515<br>  65 Â  68 -1.734810E-03 -2.038025E-03 -1.734786E-03 Â 2.038025E-03 Â 2.622505E-08 Â 1.625494E-04 Â  Â 0.17013061<br></div><div><br></div><div>I am observing that almost all the values are getting changed. So, if I am using these values to calculate some transport, I am getting different result each time (<b>I have verified it</b>). Please guide me regarding what is going on here. <b>Also, how can I get the correct values of momentum matrix so that I can get reliable results upon using them in some other calculation.</b></div><div><br></div><div><b>Thanks and Regards</b><br><br><b>Yours Sincerely</b><font color="#888888"><font color="#888888"><font color="#888888"><br>Vivek Pandey<br>Research Scholar<br>Indian Institute of Technology, Mandi<br></font><div>Himachal Pradesh, India </div></font></font></div></div>