  The increasing availability of GPUs for scientific computing has prompted
interest in accelerating quantum chemical calculations through their use. The
complexity of integral kernels for high angular momentum basis functions
however often limits the utility of GPU implementations with large basis sets
or for metal containing systems. In this work, we report implementation of $f$
function support in the GPU-accelerated TeraChem software package through the
development of efficient kernels for the evaluation of Hamiltonian integrals.
The resulting code is very efficient, and we demonstrate this with density
functional theory calculations on increasingly large alkanes as well as coupled
cluster singles and doubles calculations on water clusters. Preliminary
investigations into Ni(I) catalysis and the photochemistry of MnH(CH$_3$) are
also carried out. Overall, our GPU-accelerated software appears to be
well-suited for fast simulation of large transition metal containing systems.
