Abstract
Finding the slow manifold for two-variable ordinary differential equation (ODE) models of chemical reactions with a single equilibrium is generally simple. In such planar ODEs the slow manifold is the unique trajectory corresponding to the slow relaxation of the system as it moves towards the equilibrium point. One method of finding the slow manifold is to use direct iteration of a functional equation; another method is to obtain a series solution of the trajectory differential equation of the system. In some cases these two methods agree order-by-order in the singular perturbation parameter controlling the fast relaxation of the intermediate (complex). However, de la Llave has found a model ODE where the series method always diverges. Bimolecular association is an example of a chemical reaction where the series method for finding the slow manifold diverges but the iterative method converges. In this mechanism a complex is formed which can then undergo unimolecular decay, i.e., [reaction: see text]. The kinetics of this reaction are investigated and its properties compared with two other two-step mechanisms where series expansion and iteration methods are equivalent: the Michaelis-Menten mechanism for enzyme kinetics, and the Lindemann-Christiansen mechanism of unimolecular decay in gas kinetics.
Collapse