“…Recently, there has been significant research effort afforded to porting electronic structure software to the GPU (Gordon et al, 2020). In the case of large-scale calculations, much work has gone into the development of massively parallel GPU implementations of methods based on plane wave (Maintz et al, 2011;Wang et al, 2011;Jia et al, 2019), real space (Andrade and Aspuru-Guzik, 2013;Hakala et al, 2013), finite element (Das et al, 2019;Motamarri et al, 2020), and various other discretizations (Genovese et al, 2009;van Schoot and Visscher, 2016;Yoshikawa et al, 2019;Huhn et al, 2020) of the Kohn-Sham equations. In this work, we consider the Gaussian basis set discretization of the Kohn-Sham equations (Pople et al, 1992), which poses a number of challenges for GPU implementations.…”