We develop a variational lower bound to the free energy for stochastic reaction models with monomolecular reactions which can be used for approximate inference. This bound is based on simpler partit