redefine mass action kinetics
2 views (last 30 days)
I am looking for a general way to redefine the mass-action kinetics. The way they are implemented in Matlab is that for a reaction of type 2A->B with reaction constant k, the rate of the reaction is given by k*A^2.
I would like to change the reaction rate to k*A*(A-1).
I know how to do it for a specific type of reactions using the sbioaddtolibrary function. But this solution doesn't take care of a reaction of type 2A+2B->C.
Thank you for your help.
Arthur Goldsipe on 16 Dec 2017
It's not currently possible to create a custom kinetic law that depends on the reaction's stoichiometry. You will need to manually write each reaction rate.
Based on the form of your reaction rate, I suspect you're trying to model a system containing a small numbers of reacting molecules. Have you considered using one of SimBiology's stochastic solvers? That would likely be a better model of your system than using deterministic solver with the reaction rates you describe.