Misplaced Pages

In mathematics and mathematical biology, the Mackey–Glass equations, named after Michael Mackey and Leon Glass, refer to a family of delay differential equations whose behaviour manages to mimic both healthy and pathological behaviour in certain biological contexts, controlled by the equation's parameters. Originally, they were used to model the variation in the relative quantity of mature cells in the blood. The equations are defined as:

${\displaystyle {\frac {dP(t)}{dt}}={\frac {\beta _{0}\theta ^{n}}{\theta ^{n}+P(t-\tau )^{n}}}-\gamma P(t)}$

(Eq. 1)

and

${\displaystyle {\frac {dP(t)}{dt}}={\frac {\beta _{0}\theta ^{n}P(t-\tau )}{\theta ^{n}+P(t-\tau )^{n}}}-\gamma P(t)}$

(Eq. 2)

where ${\displaystyle P(t)}$ represents the density of cells over time, and ${\displaystyle \beta _{0},\theta ,n,\tau ,\gamma }$ are parameters of the equations.

Equation (2), in particular, is notable in dynamical systems since it can result in chaotic attractors with various dimensions.

Introduction

There exist an enormous number of physiological systems that involve or rely on the periodic behaviour of certain subcomponents of the system. For example, many homeostatic processes rely on negative feedback to control the concentration of substances in the blood; breathing, for instance, is promoted by the detection, by the brain, of high CO₂ concentration in the blood. One way to model such systems mathematically is with the following simple ordinary differential equation:

${\displaystyle y'(t)=k-cy(t)}$

where ${\displaystyle k}$ is the rate at which a "substance" is produced, and ${\displaystyle c}$ controls how the current level of the substance discourages the continuation of its production. The solutions of this equation can be found via an integrating factor, and have the form:

${\displaystyle y(t)={\frac {k}{c}}+f(y_{0})e^{-ct}}$

where ${\displaystyle y_{0}}$ is any initial condition for the initial value problem.

However, the above model assumes that variations in the substance concentration is detected immediately, which often not the case in physiological systems. In order to ease this problem, Mackey, M.C. & Glass, L. (1977) proposed changing the production rate to a function ${\displaystyle k(y(t-\tau ))}$ of the concentration at an earlier point ${\displaystyle t-\tau }$ in time, in hope that this would better reflect the fact that there is a significant delay before the bone marrow produces and releases mature cells in the blood, after detecting low cell concentration in the blood. By taking the production rate ${\displaystyle k}$ as being:

${\displaystyle {\frac {\beta _{0}\theta ^{n}}{\theta ^{n}+P(t-\tau )^{n}}}~~{\text{ or }}~~{\frac {\beta _{0}\theta ^{n}P(t-\tau )}{\theta ^{n}+P(t-\tau )^{n}}}}$

we obtain Equations (1) and (2), respectively. The values used by Mackey, M.C. & Glass, L. (1977) were ${\displaystyle \gamma =0.1}$, ${\displaystyle \beta _{0}=0.2}$ and ${\displaystyle n=10}$, with initial condition ${\displaystyle P(0)=0.1}$. The value of ${\displaystyle \theta }$ is not relevant for the purpose of analyzing the dynamics of Equation (2), since the change of variable ${\displaystyle P(t)=\theta \cdot Q(t)}$ reduces the equation to:

${\displaystyle Q'(t)={\frac {\beta _{0}Q(t-\tau )}{1+Q(t-\tau )^{n}}}-\gamma Q(t).}$

This is why, in this context, plots often place ${\displaystyle Q(t)=P(t)/\theta }$ in the ${\displaystyle y}$-axis.

Dynamical behaviour

It is of interest to study the behaviour of the equation solutions when ${\displaystyle \tau }$ is varied, since it represents the time taken by the physiological system to react to the concentration variation of a substance. An increase in this delay can be caused by a pathology, which in turn can result in chaotic solutions for the Mackey–Glass equations, especially Equation (2). When ${\displaystyle \tau =6}$, we obtain a very regular periodic solution, which can be seen as characterizing "healthy" behaviour; on the other hand, when ${\displaystyle \tau =20}$ the solution gets much more erratic.

The Mackey–Glass attractor can be visualized by plotting the pairs ${\displaystyle (P(t),P(t-\tau ))}$. This is somewhat justified because delay differential equations can (sometimes) be reduced to a system of ordinary differential equations, and also because they are approximately infinite dimensional maps.

See also

• Circadian rhythm
• Circadian oscillator
• Neural oscillation

References

1. ^ Mackey, M.C.; Glass, L. (1977). "Oscillation and chaos in physiological control systems". Science. 197 (4300): 287–9. Bibcode:1977Sci...197..287M. doi:10.1126/science.267326. PMID 267326.
2. ^ "Mackey-Glass equation". Wolfram Demonstrations Project. Retrieved 10 August 2020.
3. ^ Kantz, H.; Schreiber, T. (2004). Nonlinear time series analysis. Vol. 7. Cambridge University Press.
4. Glass, L. (2001). "Synchronization and rhythmic processes in physiology". Nature. 410 (6825): 277–84. Bibcode:2001Natur.410..277G. doi:10.1038/35065745. PMID 11258383. S2CID 4379463.
5. Specht, H.; Fruhmann, G. (1972). "Incidence of periodic breathing in 2000 subjects without pulmonary or neurological disease". Bulletin de physio-pathologie respiratoire. 8 (5): 1075–1083. PMID 4657862.
6. Rubin, R.; Strayer, D.S.; Rubin, E. (2008). Rubin's pathology: clinicopathologic foundations of medicine. Lippincott Williams & Wilkins.
7. Junges, L.; Gallas, J.A. (2012). "Intricate routes to chaos in the Mackey–Glass delayed feedback system". Physics Letters A. 376 (30–31): 2109–2116. Bibcode:2012PhLA..376.2109J. doi:10.1016/j.physleta.2012.05.022.

v t e Chaos theory
Concepts	Core Attractor Bifurcation Fractal Limit set Lyapunov exponent Orbit Periodic point Phase space Anosov diffeomorphism Arnold tongue axiom A dynamical system Bifurcation diagram Box-counting dimension Correlation dimension Conservative system Ergodicity False nearest neighbors Hausdorff dimension Invariant measure Lyapunov stability Measure-preserving dynamical system Mixing Poincaré section Recurrence plot SRB measure Stable manifold Topological conjugacy Theorems Ergodic theorem Liouville's theorem Krylov–Bogolyubov theorem Poincaré–Bendixson theorem Poincaré recurrence theorem Stable manifold theorem Takens's theorem
Theoretical branches	Bifurcation theory Control of chaos Dynamical system Ergodic theory Quantum chaos Stability theory Synchronization of chaos
Chaotic maps (list)	Discrete Arnold's cat map Baker's map Complex quadratic map Coupled map lattice Duffing map Dyadic transformation Dynamical billiards outer Exponential map Gauss map Gingerbreadman map Hénon map Horseshoe map Ikeda map Interval exchange map Irrational rotation Kaplan–Yorke map Langton's ant Logistic map Standard map Tent map Tinkerbell map Zaslavskii map Continuous Double scroll attractor Duffing equation Lorenz system Lotka–Volterra equations Mackey–Glass equations Rabinovich–Fabrikant equations Rössler attractor Three-body problem Van der Pol oscillator
Physical systems	Chua's circuit Convection Double pendulum Elastic pendulum FPUT problem Hénon–Heiles system Kicked rotator Multiscroll attractor Population dynamics Swinging Atwood's machine Tilt-A-Whirl Weather
Chaos theorists	Michael Berry Rufus Bowen Mary Cartwright Chen Guanrong Leon O. Chua Mitchell Feigenbaum Peter Grassberger Celso Grebogi Martin Gutzwiller Brosl Hasslacher Michel Hénon Svetlana Jitomirskaya Bryna Kra Edward Norton Lorenz Aleksandr Lyapunov Benoît Mandelbrot Hee Oh Edward Ott Henri Poincaré Itamar Procaccia Mary Rees Otto Rössler David Ruelle Caroline Series Yakov Sinai Oleksandr Mykolayovych Sharkovsky Nina Snaith Floris Takens Audrey Terras Mary Tsingou Marcelo Viana Amie Wilkinson James A. Yorke Lai-Sang Young
Related articles	Butterfly effect Complexity Edge of chaos Predictability Santa Fe Institute

• Chaotic maps