I would like to use QE to obtain u_K(r), which is the Bloch wave function of spin-up maximum valance band at K point. If I define the K point in the input file and do a scf calculation by QE including spin-orbit coupling, I was wondering how I can extract the mentioned Bloch wave function?
I have...