The CP-odd Rayleigh operators \( \mathcal {Q}^{(7)}_{13}, {\cal Q}^{(7)}_{14}\) 
in Eq. (REF ) do not mix into the pseudo-scalar operators in
Eq. (REF ), contrary to their QCD counterparts, due to the
absence of a QED anomaly. The RG evolution of the operators in
Eqs. (REF )-(REF ) is also implemented in the
DirectDM code [1]}.
