<div dir="ltr"><div>Dear Wien2K community,</div><div>             I am currently using the latest version of the WIEN2k package (<b>Version: WIEN2k_23.2</b>). At certain point of my work, I am in need of momentum matrix elements ( <b><<i>nk</i>|vx|<i>mk</i>></b> ), where <i>nk</i> and <i>m</i>k are the indices of Bloch states. Here, <i><b>n</b></i> & <i><b>m</b></i> denotes the band-indices and <i><b>k</b></i> stands for the particular <i><b>k</b></i>-point of the Brillouin zone. I found from here (<a href="https://triqs.github.io/dft_tools/latest/faqs/faqs.html#x-optic-does-not-write-a-case-pmat-file" target="_blank">https://triqs.github.io/dft_tools/latest/faqs/faqs.html#x-optic-does-not-write-a-case-pmat-file</a>) that one can obtain the momentum matrix elements using the <b><i>x optic</i></b> module. For this, I need to modify the <b><i>case.inop</i></b> file as mentioned in the above link. The momentum matrix elements will be contained in the <b><i>case.pmat</i></b> file.</div><div><br></div><div>I would like to highlight that <b>TaAs doesn't have inversion symmetry but possess considerable strength of SOC</b>. Also, here I am performing non-spin polarized calculations. Regarding this, the userguide reads: "<i><b>In cases of non-spinpolarized spin-orbit calculations WITHOUT inversion symmetry one must do some tricks and “mimick” a spinpolarized calculation</b></i>". In the present case, to generate the TaAs.pmat file, I have followed the following steps.</div><div><ul><li><b><i>Took a scf-converged folder (with SOC)</i></b></li><li><b><i>x kgen -fbz   (generated 8 points, just for checking purpose)</i></b></li><li><b><i>run_lapw -so -s lapw1 -e lcore</i></b></li><li><b><i>cp TaAs.vsp TaAs.vspup</i></b></li><li><b><i>cp TaAs.vsp TaAs.vspdn</i></b></li><li><b><i>cp TaAs.vectorso TaAs.vectorsoup</i></b></li><li><b><i>cp TaAs.energyso TaAs.energysoup</i></b></li><li><b><i>x lapw2 -fermi -so -up</i></b></li><li><b><i>x optic -so -up</i></b></li></ul></div><div><br></div><div>My <i style="font-weight:bold">TaAs.inop </i>file look-like:</div><div><br></div><div><b><i>500 1       number of k-points, first k-point<br>0.572903715 1.057994381 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</i></b><br><br></div><div>This produces the pmat file as (some of the initial lines of the <b><i>TaAs.pmatup</i></b> file) :</div><div>       1      79     102<br>  0.000000000000  0.000000000000  0.000000000000         1<br> (7.532302704505345E-002,8.167223407989921E-013)<br> (-2.452818645937255E-002,-5.617391534568084E-012)<br> (-2.736308582684184E-009,7.097544424488960E-006)<br> (-3.998615698872405E-002,2.777014533732855E-002)<br> (1.302109113438884E-002,8.527865972100164E-002)<br> (-1.026926403886681E-010,1.200070342578796E-010)<br> (-1.791111955708441E-002,5.236942743794264E-002)<br> (-3.899495437248428E-002,2.316018225850822E-002)<br> (3.763225250509643E-002,-7.797121255922022E-003)<br> (7.310146318006450E-002,-1.059709078143053E-002)<br> (-3.652737960637320E-002,-7.161172405686120E-002)<br> (1.068047033704979E-002,-2.846154200524592E-002)<br> (7.610044364078097E-002,6.031752549762000E-002)<br> (-3.292267660532026E-003,7.254003661389334E-002)<br> (-3.375310650422014E-010,-6.926246219766473E-009)<br> (-4.035365458096263E-002,3.022214437391330E-002)<br> (6.674399181287456E-002,-4.736686924751896E-002)<br> (7.125860558467264E-010,-4.391983650614438E-009)<br> (-1.849323572944582E-003,9.216247127675071E-005)<br> (9.836142143746737E-004,1.454250471303621E-003)<br></div><div><br></div><div>Now, I have performed the same steps on another computer. Here, I want to specially highlight that <b><i>I</i></b> <i style="font-weight:bold">Took the same scf-converged folder (with SOC) </i>and <i style="font-weight:bold">performed the same steps as mentioned above (with the same TaAs.inop file and the TaAs.klist file). Only I have done these calculations on a different PC. </i>To my surprise, now I am getting different values in the <b><i>TaAs.pmatup</i></b> file, which is shown below. </div><div>       1      79     102<br>  0.000000000000  0.000000000000  0.000000000000         1<br> (-1.214159087495217E-002,-5.600747865703096E-014)<br> (8.045331706420837E-004,-5.000950391645387E-012)<br> (1.836020562752319E-010,7.097544236721407E-006)<br> (8.866257188263187E-002,5.929833225268443E-003)<br> (-5.875000758651897E-003,8.949005485342118E-002)<br> (-9.539115301314537E-011,1.978797799446000E-010)<br> (4.155953914982438E-002,-3.846832311555534E-002)<br> (3.909488619796690E-002,5.101528227938528E-002)<br> (1.510824184241179E-003,-3.533078749450909E-002)<br> (5.739778691473076E-002,-4.492185969678203E-002)<br> (4.437764328880048E-002,4.918477575252587E-002)<br> (1.739430397746836E-003,3.387537549328652E-002)<br> (-4.779075848660323E-002,6.299577197014268E-002)<br> (6.329168974258584E-002,3.712633656957084E-002)<br> (3.923592277339212E-010,1.128917508969145E-009)<br> (3.378630332312889E-002,-6.765602685801606E-002)<br> (-6.810286018561787E-002,-4.414831302200305E-002)<br> (3.059394194239015E-010,5.748344979517646E-010)<br> (-1.368362945645098E-003,-9.656884465217919E-004)<br> (9.797714419467955E-004,-1.155838456028902E-003)<br></div><div><br></div><div>To me, this is not making any sense at all. Thus, there is a high chance that I am missing something or doing something wrong.</div><div><br></div><div><div>Kindly help me with this. Also, please correct me if I am misinterpreting something at some point.</div><div><br></div><b>Thanks and Regards</b><br><br><b>Yours Sincerely</b><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></div><div> <br></div></div>