How to plot hybrid functional band structure

Hello all,

To preface, I am using VASP and the HSE06 functional to get the electronic band structure of a bulk structure. In pymatgen it is simple to plot a band structure by reading in the vasprun.xml file using BSVasprun(), getting the band structure with get_band_structure, and putting the structure on a BSPlotter(). Provided that the KPOINTS file is given (in line-mode) with commented labels, everything goes well for non-hybrid calculations.

My issue is that I want to plot a hybrid functional band structure, which requires the use of 0-weight (fake) SC k-point procedure, which involves explicitly stating k-points from the IBZKPT file that VASP produces and appending high-symmetry points at the end of the file. My KPOINTS file is shown below.

Automatically generated mesh
22
Reciprocal lattice
0.00000000000000 0.00000000000000 0.00000000000000 1
0.25000000000000 0.00000000000000 0.00000000000000 6
0.50000000000000 0.00000000000000 0.00000000000000 3
0.25000000000000 0.25000000000000 0.00000000000000 12
0.50000000000000 0.25000000000000 0.00000000000000 12
0.50000000000000 0.50000000000000 0.00000000000000 3
0.25000000000000 0.25000000000000 0.25000000000000 8
0.50000000000000 0.25000000000000 0.25000000000000 12
0.50000000000000 0.50000000000000 0.25000000000000 6
0.50000000000000 0.50000000000000 0.50000000000000 1
0.00000000000000 0.00000000000000 0.00000000000000 0 \Gamma
0.00000000000000 0.50000000000000 0.00000000000000 0 X
0.00000000000000 0.50000000000000 0.00000000000000 0 X
0.50000000000000 0.50000000000000 0.00000000000000 0 M
0.50000000000000 0.50000000000000 0.00000000000000 0 M
0.00000000000000 0.00000000000000 0.00000000000000 0 \Gamma
0.00000000000000 0.00000000000000 0.00000000000000 0 \Gamma
0.50000000000000 0.50000000000000 0.50000000000000 0 R
0.50000000000000 0.50000000000000 0.50000000000000 0 R
0.50000000000000 0.50000000000000 0.00000000000000 0 M
0.00000000000000 0.50000000000000 0.00000000000000 0 X
0.50000000000000 0.50000000000000 0.50000000000000 0 R

Now, I have seen that there is a parameter for get_band_structure(): force_hybrid_mode=True. I have tried using this parameter but pymatgen produces a graph that does not show the band structure but some squished plot instead. Below is the plot I get.

How can I fix this and get a plot that shows the band structure for a hybrid calculation using the 0-weight kpoint method in VASP?

I managed to fix this issue.

The trick is that the plot will squeeze together if not enough unlabeled points are used in between each labeled one. So it has to be spaced like so for each symmetry point

0.00000000 0.00000000 0.00000000 0 \Gamma
0.00000000 0.03333333 0.00000000 0
0.00000000 0.06666667 0.00000000 0
0.00000000 0.10000000 0.00000000 0
0.00000000 0.13333333 0.00000000 0
0.00000000 0.16666667 0.00000000 0
0.00000000 0.20000000 0.00000000 0
0.00000000 0.23333333 0.00000000 0
0.00000000 0.26666667 0.00000000 0
0.00000000 0.30000000 0.00000000 0
0.00000000 0.33333333 0.00000000 0
0.00000000 0.36666667 0.00000000 0
0.00000000 0.40000000 0.00000000 0
0.00000000 0.43333333 0.00000000 0
0.00000000 0.46666667 0.00000000 0
0.00000000 0.50000000 0.00000000 0 X

Alternatively, removing all of the labels also works, albeit the wave vector axis will not include the symmetry point labels.