<div dir="ltr"><div>Hi users,</div><div><br></div><div>I would like to compute < \psi_n | U^\dagger H U | \psi_m>, where H is the KS Hamiltonian, \psi_n, \psi_m are the KS orbitals, and U is some unitary operator, by modifying the source code (not done separately through a separate Python script, for example), and wanted to ask for some advice. I think I can get the \psi_n, \psi_m in reciprocal space, but am not sure how to get the bra-ket. <br></div><div><br></div><div>Best regards,</div><div>Andrew<br></div></div>