Is there something I am missing, I have a discrepancy, but perhaps they are equivalent for reasons I can't see:
Should equal: $\frac{i}{\mu}r^{-3} p_mq_mq_j - \frac3\mu r^{-3}q_j - \frac{i}{\mu}r^{-1}p_j$
I obtained: $\frac{1}{2\mu}\left( p_mir^{-3}q_mq_j+ ir^{-3}q_mp_mq_j\right) $