Consider the invariant matrix element
$$\mathcal M=2ige\frac{\epsilon^{\mu \nu \rho \sigma}k_{1 \mu}\epsilon_\nu(k_1,s_1)\bar{u}(p_2,r_2)\gamma_\sigma u(p_1,r_1)q_\rho}{q^2}, \quad q=k_2-k_1$$ corresponding to the following diagram for the process $e^-+\gamma\to e^-+\pi^0$:
My goal is to compute the quantity $|\mathcal M|^2$. However, I'm not sure what is the best way to go about this calculation, as I haven't done a lot of algebra involving these objects.
A trick that I've seen before is to rewrite $\mathcal M$ after noticing that since $\nu$ is the index of a photon polarization, it must be spatial, and therefore
\begin{align*} \mathcal M=\epsilon^{0ijk}k_{10}\epsilon_i(k_1)\bar u(p_2)\gamma_ju(p_1)q_k+\epsilon^{ij0k}k_{1i}\epsilon_j(k_1)\bar u(p_2)\gamma_0 u(p_1)q_k+\epsilon^{ijk0}k_{1i}\epsilon_j(k_1)\bar u(p_2)\gamma_k u(p_1)q_0 \\ =\epsilon^{ijk}\left(-|\mathbf k_1|\epsilon_i(k_1)\bar u(p_2)\gamma_ju(p_1)q_k+-k_{1i}\epsilon_j(k_1)\bar u(p_2)\gamma_0 u(p_1)q_k+k_{1i}\epsilon_j(k_1)\bar u(p_2)\gamma_k u(p_1)q_0\right) \end{align*}
but this doesn't seem to lead to something much simpler. Any suggestions?