I was wondering if anyone could help guide me in finding the Feynman rules for a real pseudoscalar field ($\phi$) interacting with the electromagnetic field $(F^{\mu\nu})$.
The (effective) interaction Lagrangian is: $$ \phi \epsilon^{\alpha\mu\beta\nu} F_{\mu\alpha} F_{\nu\beta} $$
In my textbooks, the interaction with a complex scalar field is taken into account by imposing gauge symmetry and using the covariant derivative. Can I use a similar method for real pseudoscalar fields?
I'm hoping to evaluate the amplitude for an incoming pseudoscalar decaying to 2 photons. Much like the decay of a neutral pion to 2 photons, but using an effective Lagrangian instead of evaluating the triangle diagram.
I'm still fairly new to quantum field theory so apologies if this is ill posed or you need more information.