The use of delay differential equations in chemical kinetics

Journal of Physical Chemistry

Volumn 100, Issue 20, 1996, Pages 8323-8330

  Roussel, Marc R ^a,b  

^a MCGILL UNIVERSITY   (Canada)

^b UNIVERSITY OF LETHBRIDGE   (Canada)

Author keywords

[No Author keywords available]

Indexed keywords

CATALYSIS; DIFFERENTIAL EQUATIONS; ENZYMES; MATHEMATICAL MODELS;

DELAYED EFFECTS; LAW OF MASS ACTION; OREGONATOR MODELS;

REACTION KINETICS;

EID: 0030134446     PISSN: 00223654     EISSN: None     Source Type: Journal    
DOI: 10.1021/jp9600672     Document Type: Article

References (50)

1. Schell, M.; Ross, J. J. Chem. Phys. 1986, 85, 6489-6503.
2. Epstein, I. R. J. Chem. Phys. 1990, 92, 1702-1712.
3. Epstein, I. R.; Luo, Y. J. Chem. Phys. 1991, 95, 244-254.
4. Cheng, Z.; Hamilton, I. P.; Teitelbaum, H. J. Phys. Chem. 1991, 95, 4929-4931.
5. Chevalier, T.; Freund, A.; Ross, J. J. Chem. Phys. 1991, 95, 308-316.
6. Bar-Eli, K.; Noyes, R. M. J. Phys. Chem. 1992, 96, 7664-7670.
7. Wei, J. J. Chem. Phys. 1962, 36, 1578-1584.
8. Horn, F.; Jackson, R. Arch. Rational Mech. Anal. 1972, 47, 81-116.
9. Perelson, A. S.; Wallwork, D. J. Chem. Phys. 1977, 66, 4390-4394.
10. Schuster, S.; Schuster, R. Z. Phys. Chem. 1990, 271, 337-345.
11. Curtiss, C. F.; Hirschfelder, J. O. Proc. Natl. Acad. Sci. U.S.A. 1952, 38, 235-243.
12. Heinrich, R.; Rapoport, S. M.; Rapoport, T. A. Prog. Biophys. Mol. Biol. 1977, 32, 1-82.
13. van Kampen, N. G. Phys. Rep. 1985, 124, 69-160.
14. Okuda, M. Prog. Theor. Phys. 1981, 66, 90-100.
15. Fraser, S. J. J. Chem. Phys. 1988, 88, 4732-4738.
16. Roberts, A. J. SIAM J. Math. Anal. 1989, 20, 1447-1458.
17. Saaty, T. L. Modern Nonlinear Eqtiations; Dover: New York, 1981; Chapter 5.
18. Rather than constrain the class of equations at the outset, we can always modify the equations by a judicious, if somewhat ad hoc, application of Heaviside functions, as was done in ref 2. A more elegant general solution is sought here.
19. max may be infinite in the case of distributed delays.
20. Brown, A. J. J. Chem. Soc. 1902, 81, 373-388.
21. Henri, V. C. R. Hebd. Seances Acad. Sci. 1902, 135, 916-919.
22. Okuda, M. Prog. Theor. Phys. 1982, 68, 1827-1840.
23. Nguyen, A. H.; Fraser, S. J. J. Chem. Phys. 1989, 91, 186-193.
24. Roussel, M. R.; Fraser, S. J. J. Phys. Chem. 1991, 95, 8762-8770.
25. Roussel, M. R. Ph.D. Thesis, University of Toronto, 1994.
26. Quellet, L.; Laidler, K. J. Can. J. Chem. 1956, 34, 146-150.
27. Darvey, I. G. J. Theor. Biol. 1968, 19, 215-231.
28. Ryde-Pettersson, U. Eur. J. Biochem. 1989, 186, 145-148.
29. Roussel, M. R.; Fraser, S. J. J. Chem. Phys. 1990, 93, 1072-1081.
30. Roussel, M. R.; Fraser, S. J. J. Chem. Phys. 1991, 94, 7106-7113.
31. A very simple strategy was adopted for numerical integration: The integrator was a straightforward constant-step-size implicit trapezoidal code. Values obtained by integration were stored in a circular buffer. The routine which returns the function to be integrated extracts the history-dependent , terms from the buffer by linear interpolation. To verify the plausibility of the trajectories obtained, the results of the integration were compared to analytic results whenever possible.
32. Fraser, S. J.; Roussel, M. R. Can. J. Chem. 1994, 72, 800-812.
33. These are functions with a random component rather than truly random functions. These functions were constructed in the following manner: Since a DDE integrator uses stored points and interpolation to compute the history-dependent terms, it is necessary only to store a set of points to initialize the integrator. Random initial values for the concentrations are chosen in some positive interval. These values are then incremented by random amounts uniformly selected in a small interval at each mesh point. Negative values, when they arise, are replaced by small positive values in order to maintain the nonnegativity constraint imposed on these models. Finally, the stored values are scaled to satisfy any additional normalization conditions.
34. This procedure is similar to numerical methods for constructing the slow manifold of enzymic mechanisms; see Roussel, M. R. M.Sc. Thesis, University of Toronto, 1990. It is also closely related to seminumerical methods for constructing integral manifolds of PDE systems. The boundary of the manifold is constructed (analogous to finding an initial function which generates the slow manifold), and the manifold is then found by integration from the boundary; see Constantin, P.; Foias, C.; Nicolaenko, B.; Temam, R. Integral Manifolds and Inertial Manifolds for Dissipative Partial Differential Equations; Springer-Verlag: New York, 1989.
35. Wischert, W.; Wunderlin, A.; Pelster, A.; Olivier, M.; Groslambert, J. Phys. Rev. E 1994, 49, 203-219.
36. Field, R. J.; Noyes, R. M. J. Chem. Phys. 1974, 60, 1877-1884.
37. Field, R. J.; Försterling, H.-D. J. Phys. Chem. 1986, 90, 5400-5407.
38. Luo, Y.; Epstein, I. R. Adv. Chem. Phys. 1990, 79, 269-299.
39. This scaling was chosen to make the steady-state values of the concentration variables of unit magnitude.
40. Because of the extreme stiffness of the Oregonator and delayed Oregonator equations, a slightly different numerical integration strategy had to be adopted than for the other models studied here. A third-order Gear method was used with a conservative step-size adjustment policy: The step size was decreased whenever either the solution threatened to become negative or the functional iteration routine at the core of the integrator failed to converge reasonably quickly. After a suitably large number of steps, the step size was increased again. Otherwise, the method was exactly as described in note 31 above.
41. Busenberg, S. N.; Travis, C. C. J. Math. Anal. Appl. 1982, 89, 46-66.
42. Cooke, K. L.; Grossman, Z. J. Math. Anal. Appl. 1982, 86, 592-627.
43. Macdonald, N. Biological Delay Systems: Linear Stability Theory; Cambridge: Cambridge, 1989; Chapter 10.
44. The universal acceptance of this idea is relatively modern. In the early days of enzymology, it was believed by many that enzymes catalyzed only exothermic reactions. For a discussion, see Fruton, J. S. Molecules and Life; Wiley: New York, 1972; p 83.
45. Noyes, R. M. J. Phys. Chem. 1977, 81, 2315-2316.
46. Hearon, J. Z. Ann. N.Y. Acad. Sci. 1963, 108, 36-68.
47. Shear, D. J. Theor. Biol. 1967, 16, 212-228.
48. Higgins, J. J. Theor. Biol. 1968, 21, 293-304.
49. Shear, D. B. J. Chem. Phys. 1968, 48, 4144-4147.
50. Bélair, J. J. Dyn. Diff. Eq. 1993, 5, 607-623.