# Bloch wave function extraction for the specific k-points and specific band from QE calculation [closed]

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:

What is the best program to manipulate numerical DFT wavefunctions to calculate custom matrix elements?

to get .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.