Feynman Rules for massive vector boson interactions I am stuck at the beginning of a problem where I am given an interaction term that modifies the regular QED Lagrangian. It involves the interaction between a photon field and a massive vector boson: $\mathcal{L}_{int_1} = \frac{1}{4}g_1G_{\mu \nu}F^{\mu \nu}$. Here, $F^{\mu \nu}$ is the electromagnetic field tensor. Similarly, $G_{\mu \nu} = \partial_{\mu}B_{\nu}-\partial_{\nu}B_{\mu}$ where $B_{\mu}$ is the massive vector field.
I am trying to derive the Feynman rule for this vertex. I am not sure what to do. Usually, I just bring down a 4-momentum when I have a field with derivatives in the Lagrangian, but in this case, it looks like two of the 4-momenta would end up forming a dot product together, while the other two would not, and I would be left with a term without any indices and two terms with indices, which doesn't make sense. Does anyone know what the Feynman rule for this vertex is and why it is so?
EDIT: I think the Feynman rule might be $i\mathcal{M} = -i g_1 (g_{\mu \nu}-\frac{k_\mu a_\nu}{k \cdot q})$. I obtained this by an analogous computation for finding the photon propagator. Is doing this allowed?
The main problem that I am trying to solve involves the diagram $f\bar{f} \to \phi\phi^*$. From comments, I realize that I have two options: 1) I can diagonalize the mass matrix and use the new mass as the mass in the propagator between the two end vertices or 2) I can consider the infinite series of $A$ propagator > $B$ propagator > ... > $A$ propagator > $B$ propagator and obtain the same answer. The problem I am having with diagonalizing the mass matrix is that in previous problems where I have diagonalized a mass matrix, my interaction term did not contain derivative terms. So my attempt at a solution (which I am trying now) is to just bring down a 4-momentum and go from there. But then wouldn't that make my mass momentum dependent? With the second method, I think I would need to still know the feynman rule for the interaction vertex, which circles back to my original problem.
 A: Your first interaction term is bilinear in 2 gauge fields. Terms of this form indicate that you need to diagonalize your mass matrix. So you could work out the Feynman rules for your first interaction, but your question is sort of a moot point since you won't have to do perturbation theory for this term, it will get assimilated into the propagators in the diagonal basis. 
EDIT: 
Just so you get your orginial question answered, the Feynman rule for your first interaction is obtained by Fourier transforming the Action:
$S =  \frac{g_1}{4}\int d^4 x F^{\mu} G_{\mu \nu}   
= \frac{g_1}{4}\int d^4 x \left( \partial_\mu F_\nu -\partial_\nu F_\mu \right)  \left(  \partial^\mu G^\nu - \partial^\nu G^\mu  \right)
$
$=\frac{g_1}{2}\int d^4 x \left( F_\nu \partial_\mu \partial^\nu G_\mu - F_\nu \Box G_\nu    \right) = \int d^4 p \tilde{F}_{\nu} (-p) \left[ \frac{g_1}{2} \left( p^2 \eta^{\nu \mu} - p^\nu p^\mu  \right)   \right] \tilde{G}_\mu (p)
$
where the Feynman verte rule is in the square brackets now. 
Now, to rephrase my original point, the Lagrangian
$\mathcal{L} = \frac{1}{4} F^{\mu \nu} F_{\mu \nu} + \frac{1}{4} G^{\mu \nu} G_{\mu \nu} +m^2 G^\mu G_\mu+ g_1 \frac{1}{4} F^{\mu \nu} G_{\mu \nu}$
is trivial, in that it is free. You can see this by doing perturbation theory, or you can make things easy on yourself and diagonalize the quadratic portion of the action.
Let me know if you want further clarification.
