eval_rate_law = reaction.rate_laws[0].expression._parsed_expression.eval({
                        wc_lang.Species: mito_species_counts,
                        wc_lang.Compartment: {
                            mitochondria.id: mitochondria.init_volume.mean * \
                                mitochondria.init_density.value}