Abstract:
Oscillations in cytoplasmic calcium concentration are a crucial control mechanism in almost every cell type. Two important classes of oscillation are of particular interest: solitary and periodic waves. Both types of waves are commonly observed in physical experiments and found in mathematical models of calcium dynamics and other excitable systems. In this thesis, we try to understand these two classes of wave solutions. We first investigate wave solutions of the canonical excitable model, the FitzHugh-Nagumo (FHN) equations. We analyze the FHN equations using geometric singular perturbation theory and numerical integration, and find some new codimension-two organizing centres of the overall dynamics. Many analytical results about the FHN model in its classical form have already been established. We devise a transformation to change the form of the FHN equations we study into the classical form to make use of the results. This enables us to show how basic features of the bifurcation structure of the FHN equations arise from the singular limit. We then study waves of a representative calcium model. We analyze the dynamics of the calcium model in the singular limit, and show how homoclinic and Hopf bifurcations of the full system arise as perturbations of singular homoclinic and Hopf bifurcations. We compare the wave solutions in the FHN model and the calcium model, and show that the dynamics of the two models differ in some respects (most importantly, in the way in which diffusion enters the equations). We conclude that the FHN model should not uniformly be used as a prototypical model for calcium dynamics. Motivated by phenomena seen in the FHN and calcium models, we then investigate reduction techniques for excitable systems, including the quasi-steady state approximation and geometric singular perturbation theory, and show that criticality of Hopf bifurcations may be changed when applying these reduction methods to slow-fast biophysical systems. This suggests that great care should be taken when using reduction techniques such as these, to ensure that spurious conclusions about the dynamics of a model are not drawn from the dynamics of a reduced version of the model. Finally, we describe the class of numerical algorithms used to compute features of the detailed bifurcation sets for the FHN and calcium models, and show how these were used to locate a non-structurally stable heteroclinic connection between periodic orbits in a calcium model; this is the first time such a global bifurcation has been computed.