<div dir="ltr"><div>Hello Everyone,</div><div><br></div><div>The common interband matrix element <m| dH/dkx \pm i dH\dky |n> seems to be calculable using DFT software. However, I am trying to calculate the tangential matrix element <m| dH/dt | n> \equiv <m| dH/dkx dkx/dt + dH/dky dky/dt |n> for some parameter t that parameterizes a closed loop in k-space: (kx(t),ky(t)) = (kx0+cos(t), ky0+sin(t)), where the loop center is (kx0,ky0). With my basic understanding of QE, I thought of a way to do this using a monolayer MoS2 example: <br></div><div>I first sample k-points along a closed loop. This gives me wavefunctions in the .dat format. I was hoping to somehow extract wavefunction data to smoothen the gauge using the maximally localized Wannier function method, and manually calculate the matrix element I want using a central difference method with MATLAB or Python. <br></div><div>However, I cannot seem to read the data from the .dat files. I saw several solutions in the archives that suggest using different existing routines that involve the evec variable. However, I am not familiar with Fortran90 and so have trouble figuring out how exactly to use such routines and custom code to get eigenvector data in a readable format. Does anyone have any advice or resources? I am using QE v6.6<br></div><div>Alternatively, are there other more-direct ways to go about calculating the quantity I want? <br></div><div>Thank you very much for your time!</div><div><br></div><div>Tharindu W. Fernando</div><div>PhD Student, University of Washington<br></div><div><br></div></div>