make polymerization return the modification concentration
Factor out species_t into a common header
move exceptions to their own file
reserve count vectors before resizing

This way, we don't use too much memory for them.