than the typical application where e.g. aromatic rings would push the ratio upward. This example does employ d-functions in the basis set for oxygen and the associated integral code kernels are thus here brought into action. The reason why this deserves to be mentioned is that these kernels are much more extended (number of source code lines) than those involving only lower angular momenta and thereby more difficult to launch efficiently on GPUs. This example does, however, not employ any diffuse basis functions that would hamper screening and push the true scaling of the implementation to be observed for larger system sizes. We argue, however, that DFT linear response properties for large-scale 3D systems often are reasonably well described with a basis set such as def2-SVP such that this example becomes relevant for real case studies. The observed scaling with respect to system size (or number of basis functions) is N , which is a testament of the quadratic scaling of the J matrix mixed with the incipient linear scaling of the K matrix.