In case of ideal CRNs, the fluxes \(j^{\pm \rho }({z})\)  are expressed in terms of mass-action kinetics [1]}, [2]}, [3]}:
\(j^{\pm \rho }({z}) = {k}^{\text{id}}_{\pm \rho }{z}^{\nu _{\pm \rho }}\,,\) 
