I am currently trying to figure out how to compute band structures for my system, using the hybrid functional HSE06. I'm doing this on Quantum Espresso. As I understand, there are a handful of ways to do this:
1) Generate maximally localized wannier functions (MLWFs) with the HSE functional turned on
2) The 'fake scf' procedure: This is similar to what is listed in Vaspwiki(here).
3) Generate scf data for a coarse 'q' mesh (of the Fock operator) and then interpolate these bands using softwares like BoltztraP2.
I am not familiar with the Wannier module on Quantum ESPRESSO, hence option (1) is sort of last resort for me, since Wannier takes some time to learn.
The 'fake scf' procedure doesn't seem to work. The k-point mesh for my original calculation (my system is a bilayer TMD - Transition metal dichalcogenide) is 12x9x1. I tried using a commensurate 'q' grid and tried the 'fake scf' procedure, but I couldn't get it to work.
I would be grateful if anyone could help me figure out how to compute these bandstructures - either through one of the methods I mentioned, or any easy method as you see fit. I would like to add that I have access to considerable computational resources, so doing expensive calculations should not be much of an issue.
Also, I would like to say that I am open to methods that use other softwares - including VASP or CASTEP. I have both licenses (CASTEP through Materials Studio). But things can get a little tricky here because VASP for example, uses a different implementation of the Hubbard 'U' and it would corrupt my results if I port my QE input to VASP to calculate the band-structures.