  Electric dipole moments are sensitive probes of new phases in the Higgs
Yukawa couplings. We calculate the complete two-loop QCD anomalous dimension
matrix for the mixing of CP-odd scalar and tensor operators and apply our
results for a phenomenological study of CP violation in the bottom and charm
Yukawa couplings. We find large shifts of the induced Wilson coefficients at
next-to-leading-logarithmic order. Using the experimental bound on the electric
dipole moment of the neutron, we update the constraints on CP-violating phases
in the bottom and charm quark Yukawas.
