A massively multicore parallelization of the Kohn-Sham energy gradients