The effective masses of the holes were computed using a finite-difference (FD) derivative of the band energies around the VBM. A 5-point FD stencil was used. The step size between the stencils was chosen as 1/30 of the distance between Γ and the high-symmetry points 𝐾 or 𝑀 to adapt to the changing size of the Brillouin zone. The validity of this choice of the grid was confirmed by converging the effective masses with respect to the step size Δ between the stencils. This is shown in Figure .