The Polymer should know how to calculate all internal polymer energies, and the field should know how to calculate all field energies. The Polymer should automagically log all contributors to field energies.
In order to make this fast, we should have the polymer return a JIT compiled function that takes numpy array representations of the field and polymers and computes the energy of those states. We should then just use that energy function in mc_sim directly.