While trying to use the gammapy backend, I stomped over some error related to the geometries of the arrays when evaluating the model.
It appears that Synchrotron.evaluate_tau_ssa() is flattening the emissivity array over energies, causing this error:
ValueError: operands could not be broadcast together with shapes (38,12) (456,)
Here 38 is the number of energy bins