I'm trying to calculate the neutrino magnetic moment in the theory with this additional term in the Lagrangian: $\frac{a}{M^2}(\bar{\nu}\sigma_{\mu\nu}\nu)(\bar{e}\sigma^{\mu\nu}e)$, where $\sigma^{\mu\nu}=\frac{i}{2}[\gamma^{\mu},\gamma^{\nu}]$. This is the interaction of neutrino with electron (for simplicity I'm considering only electrons and electron neutrinos) As far as I understand, I should calculate a vertex $\Gamma^{\mu}$, corresponding to this diagram.
Then I should split it up to obtain a term proportional to $\frac{i\sigma^{\mu\nu}q_{\nu}}{2m}$ ($q = \bar{p}-p$), which would be magnetic moment.
But I don't quite understand how to write down an integral for this diagram. I'm thinking of something like this ($p$ for neutrino, $\bar{p}$ for antinutrino): $$\int\frac{d^4k}{(2\pi)^4}\frac{\gamma_{\mu}}{\gamma (p+k)-m+i\epsilon}\frac{\sigma^{\mu\nu}}{\gamma (\bar{p}+k)-m+i\epsilon}$$ Howewer I'm not sure about upper and lower indices, since in usual cases there is a boson propogator proportional to $g_{\mu\nu}$, which provides the correct contraction, but not in this case. On the other hand two indices must be contracted because $\Gamma^{\mu}$ has one upper index.
P.S. I know, that the integral is divergent, but nevertheless I'd like to know how to write it down.