How can we expand the operator $$O_{12} = \frac{1}{r^3}[({\bf \sigma_{1}\cdot{\bf r}})({\bf \sigma_{2}\cdot{\bf r}}) - r^2{\bf \sigma_{1}}\cdot {\bf \sigma_{2}}]$$ where ${\bf \sigma} = (\sigma_{x}, \sigma_{y}, \sigma_{z})$ and ${\bf r}$ is an ordinary 3-vector. Subscript 1, 2 refers to the first and second spin.
I have tried using the relation $(\sigma . {\bf a})(\sigma . {\bf b})$ but that only involves one of the spins.
Can someone help me?