Two-loop QCD renormalization and anomalous dimension of the scalar diquark operator