Dear QE Developers and Users,<div><br></div><div>I am developing code to interpolate band structures, which necessitates the use of orthogonal wavefunctions. When working with PAW or ultrasoft pseudopotentials, the wavefunctions generated are nonorthogonal. I understand that orthogonalization can theoretically be achieved using a transformation operator T, such that psi_orthogonal = T psi_nonorthogonal. However, I am uncertain how to implement or obtain this transformation operator T within the QE code.</div><div><br></div><div>Could anyone provide insights or methods for computing T in the QE framework? Either real space or reciprocal space would be acceptable.</div><div><br></div><div>Thank you for your assistance and valuable time.</div><div><br></div><div>Best regards,</div><div>Kai Wu</div>