In the DFT community, it is common practice to use regular k-point grids (Monkhorst-Pack, MP) for Brillouin zone integration. Recently Wisesa et. al. [1] and Morgan et. al. [2] demonstrated that generalized regular (GR) grids offer advantages over traditional MP grids. GR grids have not been widely adopted because one must search through a large number of candidate grids. This work describes an algorithm that can quickly search over GR grids for those that have the most uniform distribution of points and the best symmetry reduction. The grids are ~60% more efficient, on average, than MP grids and can now be generated on the fly in seconds.

1. P. Wisesa, K. A. McGill, and T. Mueller, Phys. Rev. B 93, 155109 (2016).
2 W. S. Morgan, J. J. Jorgensen, B. C. Hess, and G. L. W. Hart, Computational Materials Science 153, 424 (2018).