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 an SCF calculation by QE including spin-orbit coupling, I was wondering how I can extract the mentioned Bloch wave function?
I have followed the below link:
.hdf5 files. This link, presents how to collect wavefunctions.
I am not sure how to use the script together with python to obtain the desired wave function. I would appreciate your help and advice on this.