fig4a_dft_GM.dat: Data for the DFT bands on the G-M path in Fig.4(a)
column 1: k-path, column 2: energy (eV)
(k-path G: 0.00000000, M: 0.13751346)

fig4a_dft_MK.dat: Data for the DFT bands on the M-K path in Fig.4(a)
column 1: k-path, column 2: energy (eV)
(k-path M: 0.00000000, K: 0.07939343)

fig4a_dft_KG.dat: Data for the DFT bands on the K-G path in Fig.4(a)
column 1: k-path, column 2: energy (eV)
(k-path K: 0.00000000, G: 0.15878687)

fig4a_wannier.dat: Data for the TB model in Fig.4(a)
column 1: k-path, column 2: energy (eV)
(Note: In the plot, the energy is shifted by 0.0985472 eV, and the k-path is divided by 2pi.) 

fig4b.dat: Data for Fig.4(b)
column 1: k-path, columns 2-33: energies of bands 1-32 (eV), columns 34-65: weights of Gd d orb. on bands 1-32 
(k-path G: 0.00000000, M: 0.36666667, K: 0.57666667, G: 1.00000000)

wannier90_hr_single_GdGaI.dat: Hopping parameters from wannier90
(Ga-1 p orb.: 1-3, Ga-2 p orb.: 4-6, Gd-1 d orb.: 7-11, Gd-2 d orb.: 12-16)
