Figure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., 2002.

Size: px
Start display at page:

Download "Figure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., 2002."

Transcription

1 1 Traveling Fronts Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306

2 2 When mature Xenopus oocytes (frog eggs) are loaded with a Ca2+ dye such as fura-2 and stimulated by the fusion of sperm, a wave of elevated Ca2+ is produced in the large egg ( 1.2 mm in diameter). This is called a fertilization Ca2+ wave and it is important for triggering the formation of an envelope around the oocyte that prevents the fusion of other sperm (which would result in polyspermy). This process occurs in eggs of other species, from starfish to mammals. Figure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., This fertilization wave is an example of a traveling front, where the variable (Ca2+ concentration in this case) goes from a resting low level to a stimulated high level. This can be described as a bistable system, and the front is perturbing the system from the resting state to the co-existing stimulated state. The wave can be observed using a space-time plot, as shown below.

3 3 Figure 2: Space-time plot of a traveling Ca2+ wave. From Fall et al., This can be viewed in a different way by plotting the Ca2+ concentration versus time at several values of x, where x is the distance from sperm fusion, as shown in Fig. 3. Figure 3: Time dynamics for several values of the distance from sperm fusion. From Fall et al., The fertilization wave is described by the partial differential equation c = D 2c + f (c) t (1)

4 where c is the intracellular Ca 2+ concentration, D is the diffusion coefficient multiplied by the fraction of Ca 2+ that is free, and f(c) is a reaction term that describes the net Ca 2+ flux from the ER into the cytosol. The reaction term is f(c) = [k leak + k IP3 ][c ER c] k SERCA c (2) 4 where the first terms describe influx through ER membrane leak and IP 3 receptors, and the last term is efflux through SERCA pumps. This is an example of a reaction-diffusion system that gives rise to traveling front solutions for some parameter values. Rather than studying this biophysical model, we will focus on simpler reaction-diffusion systems that also produce traveling fronts. Fisher s Equation In a paper published by Fisher in 1937 the spread of an advantageous allele of a gene (call it a ) into a region in which a second allele (call it A ) was initially present was considered. The spread was assumed to be random. If p = Prob[a is present], then 1 p = Prob[A is present]. Using the standard rules of genetics (Hardy-Weinberg genetics), it can be shown that the rate at which the frequency p evolves over time is described by p t = D 2 p + αp(1 p) (3) x2

5 where D is a diffusion or mixing rate and α is the intensity of selection. This is called Fisher s Equation. We will analyze the system using D = α = 1. Then, p t = 2 p + p(1 p) (4) x2 This PDE has two constant solutions, p 1 = 0 and p 2 = 1. The stability of these solutions can be observed by plotting the reaction term (Fig. 4). Clearly, the equilibrium at 1 is stable, while the equilibrium at 0 is not. f(p) 5 1 p Figure 4: Analysis of the Fisher reaction term shows stability properties of equilibria. In terms of the biology, this means that, ultimately, the probability that an individual has the a allele is 1, and that it has the A allele is 0. This is an asymptotic result, but what about the temporal dynamics? As with the fertilization wave, the spread of the advantageous allele can occur as a traveling front. This can be viewed in the (t, x)-plane (Fig. 5). Here, p = 0 ahead of the wave front, and p = 1 behind the wave front. The

6 front is assumed to move at a constant velocity c, so the shape of the front in the (t, x)-plane is a line with slope c. This line has a vertical-intercept of z f. The value of p for any other point in this plane depends on whether it lies above or below the front line, and this in turn is determined by whether the line of slope c containing that point has a vertical-intercept z that is greater than or less than z f. That is, if a point (t, x) is on a line of slope c with vertical-intercept z, then the value of p at that point depends entirely on the value of z. Therefore, we look for a traveling front solution, which would be a function depending only on z. Here, x = ct + z (since slope of the line is c and z is the vertical-intercept), so z = x ct. (5) We call this coordinate z the traveling coordinate. The traveling x p=0 z (t,x) front 6 z f p=1 Figure 5: Traveling wave as viewed in the (t, x)-plane. The traveling coordinate z is the verticalintercept of the line with slope c through a point (t, x). t

7 7 wave solution (if it exists) has the form p(x, t) = u(z) = u(x ct) (6) where c is the unknown front speed. Our goal is to pick the true front speed, and then find the traveling front solution u(z). Now So the traveling front would satisfy p = du dz t dz dt = cdu dz (7) 2 p x = d2 u 2 dz 2 (8) p(1 p) = u(1 u). (9) c du dz = d2 u + u(1 u). (10) dz2 This can be converted into two first-order ODEs using v = du dz, so du dz = v (11) dv dz = cv u(1 u). (12) ( ( 0 1 The system Eqs. 11, 12 has two steady states, u 1 = and u 0) 2 =. 0) For stability we linearize the system and look at the eigenvalues of the Jacobian matrix, J = ( ) 0 1 2u 1 c. (13) Evaluated at u 1 the eigenvalues are λ = ( c ± c 2 4)/2 and since c > 0 they have negative real parts. They are real when c 2 and

8 8 complex otherwise. Hence, u 1 = { stable node when c 2 stable spiral when c < 2. (14) At u 2 the eigenvalues are λ = ( c± c 2 + 4)/2. There will be one positive real eigenvalue and one negative real eigenvalue, so u 2 is a saddle point. Hence, u 1 has a 2-dimensional stable manifold and u 2 has 1-dimensional stable and unstable manifolds. v u 2 u Figure 6: The Fisher model with c < 2. A heteroclinic orbit (green) connects the saddle point at u 2 to the stable spiral at u 1. A couple of additional trajectories are shown (black). This is non-physical since it would require some negative values of u. In the case of c < 2 (Fig. 6) the trajectories will have some negative values of u as they spiral in to u 1. This includes the heteroclinic orbit that connects one branch of the u 2 unstable manifold to the stable spiral at u 1. Since u is a probability, this is non-physical. In the case of c 2 (Fig. 7) the trajectories enter u 1 without spiraling, so u can remain positive, and thus provide physical solutions.

9 9 v u 2 u Figure 7: The Fisher model with c 2. A heteroclinic orbit (green) connects the saddle point at u 2 to the stable node at u 1. This can be done without negative u values. The traveling wave solution starts with the system in a state in which the population is dominated by one allele of the gene and is eventually dominated by the invading allele. In terms of the system of 2 ODEs, this corresponds to the heteroclinic orbit in Fig. 7. In this case, u(z) 1 as z and u(z) 0 as z. A sketch of the heteroclinic orbit as a function of z is shown in Fig. 8. This wave has the shape of a moving u 1 z Figure 8: The Fisher model with c 2. The heteroclinic trajectory is shown as a function of the traveling coordinate z.

10 front. For populations on the left side of the diagram the invading allele a is dominant (u is near 1); for those on the right the original allele A is dominant (u is near 0). The analysis above was performed in the frame of reference of traveling coordinate z = x ct. In this frame, the wave is stationary. In the frame of reference of a stationary observer there is a traveling front moving from left to right (since c > 0, Fig. 9). One can also view this by looking at time courses at different spatial locations (Fig. 10). u 1 10 t t t x Figure 9: Traveling front moving from left to right in the stationary frame of reference (t 0 < t 1 < t 2 ). u 1 x 0 x 1 Figure 10: Time courses showing the traveling pulse at two different locations (x 0 < x 1 ). t

11 Eventually u 1 for any location x and thus the invading allele a dominates the original (and not as fit) allele A in the population. In this model, traveling front solutions exist for any front velocity c 2. That is, there is a continuum of traveling fronts, each with a different front speed. Which one actually occurs depends upon the initial conditions. Spreading Microorganisms In a typical laboratory experiment on yeast, one fills the bottom of a petri dish with a solidified gel called agar, which permits the diffusion of small molecules like glucose. Then, a small number of yeast cells are put on top of the agar. The glucose below provides the necessary nutrients for the yeast. The yeast reproduce and migrate in directions where glucose is available. A 1-dimensional model of yeast spread can be constructed by keeping track of the density of yeast and concentration of glucose: n(x, t) = density of yeast g(x, t) = glucose concentration. As in the last example, we will consider an infinite domain, so x (, ) with no boundaries. Then the yeast PDE is n t = D 2 n n x + Kn(g g 1) (15) 2 11

12 12 where D n 2 n x 2 reflects random diffusion-like migration of yeast and Kn(g g 1 ) reflects yeast proliferation. Cells are assumed to increase proportionally to the amount of glucose in excess of g 1. The glucose PDE is g t = D 2 g g x ρkn(g g 1). (16) 2 The D g 2 g x 2 term represents glucose diffusion. The ρkn(g g 1 ) term represents glucose depletion by yeast at a rate of ρ units of glucose per new cell made. To simplify the analysis we omit the yeast diffusion term, which should be much slower than glucose diffusion. That is, the spread of yeast will be due entirely to the diffusion of glucose. Also, we define ĝ(x, t) g(x, t) g 1. Then, n t ĝ t = Knĝ (17) 2 ĝ = D g ρknĝ. (18) x2 We now seek traveling front solutions, using N(z) = n(x, t) (19) G(z) = ĝ(x, t) (20) where z = x ct. Inserting N and G into the differential equations, we

13 13 obtain the following ordinary differential equations: c dn dz c dg dz = kng (21) = D g d 2 G ρkng. (22) dz2 This could be converted into a set of three first-order ODEs, but then phase plane analysis would not be possible. Instead, multiply the first equation by ρ and add it to the second to obtain Integrate this from to z, cρ dn dz cdg dz = D d 2 G g. (23) dz 2 or cρ z dn ds ds c z dg z ds ds = D g d 2 G ds (24) ds2 cρn(z)+cρn( ) cg(z)+cg( ) = D g dg dz D g dg ( ) (25) dz What does z mean? It means that ct x, or that the location x is far behind the wave front. At this location, one expects the yeast to have reached some carrying capacity N o and the glucose to be at a level at which yeast exhibit neither net proliferation nor death. That is, G = 0 since then g = g 1. Also, one would not expect much change in G as z is

14 14 increased a small amount. Altogether, so Eq. 25 becomes N( ) = N o (26) G( ) = 0 (27) dg ( ) = 0 dz (28) cρn(z) + cρn o cg(z) = D g dg dz. (29) We thus have two first-order ODEs, dg dz = c G cρ N + cρ N o D g D g D g (30) dn dz = k NG. c (31) Any traveling front solutions must be heteroclinic orbits of these equations. The traveling front equations have two equilibria, (0, N o ) and (ρn o, 0). The first has yeast at its carrying capacity and glucose g = g 1 (so that G = 0). The second has no yeast and glucose at g = ρn o + g 1 (so G = ρn o ). The first equilibrium is a saddle point while the second is a stable node. The heteroclinic orbit is drawn in green in Fig. 11, connecting the saddle point to the stable node. It is sketched in Fig. 12 as a function of the traveling coordinate. So this example again shows how a traveling front can come about, in this case for the population of yeast. Since the traveling

15 15 N No ρno G Figure 11: Phase portrait for the density of yeast N and glucose G. The heteroclinic orbit is drawn in green. coordinate is z = x ct and c > 0, the front is moving to the right when viewed as N vs. x at different points in time. Hence, at any location x, the glucose level will eventually be g = g 1 and the yeast number will be at the carrying capacity N = N o. N G Figure 12: Traveling front of cell number (N) and glucose concentration (G) as a function of the traveling coordinate z. z In this example, the only constraint on the front speed c is that it be positive so that N moves to its stable carrying capacity as t. In both examples of traveling fronts the system was described in terms

16 of planar ODEs. In the second example, we were lucky; the reduction to two dimensions was the result of a reduction of order. Other PDEs and systems of PDEs can have traveling front solutions, but the analysis is harder since the system of ODEs that must be analyzed has more than two dimensions. 16

17 17 References Mathematical Models in Biology, L. Edelstein-Keshet, SIAM Classics in Applied Mathematics 46, 2005.

Lecture 5: Travelling Waves

Lecture 5: Travelling Waves Computational Biology Group (CoBI), D-BSSE, ETHZ Lecture 5: Travelling Waves Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015 26. Oktober 2016 2 / 68 Contents 1 Introduction to Travelling Waves

More information

= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :

= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F : 1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change

More information

Travelling waves. Chapter 8. 1 Introduction

Travelling waves. Chapter 8. 1 Introduction Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part

More information

Introduction LECTURE 1

Introduction LECTURE 1 LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in

More information

1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics

1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics 1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics Theorem. Consider the equation u t = Du xx + f(u) with f(0) = f(1) = 0, f(u) > 0 on 0 < u < 1, f (0)

More information

Math 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( )

Math 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( ) #7. ( pts) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, λ 5 λ 7 t t ce The general solution is then : 5 7 c c c x( 0) c c 9 9 c+ c c t 5t 7 e + e A sketch of

More information

Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells

Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University

More information

Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics

Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:

More information

Classification of Phase Portraits at Equilibria for u (t) = f( u(t))

Classification of Phase Portraits at Equilibria for u (t) = f( u(t)) Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the

More information

Answers to Problem Set Number 02 for MIT (Spring 2008)

Answers to Problem Set Number 02 for MIT (Spring 2008) Answers to Problem Set Number 02 for 18.311 MIT (Spring 2008) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139). March 10, 2008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics,

More information

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This

More information

8.1 Bifurcations of Equilibria

8.1 Bifurcations of Equilibria 1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

More information

Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

More information

Dynamical Systems in Neuroscience: Elementary Bifurcations

Dynamical Systems in Neuroscience: Elementary Bifurcations Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

+ i. cos(t) + 2 sin(t) + c 2.

+ i. cos(t) + 2 sin(t) + c 2. MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,

More information

Problem set 7 Math 207A, Fall 2011 Solutions

Problem set 7 Math 207A, Fall 2011 Solutions Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

1 Assignment 1: Nonlinear dynamics (due September

1 Assignment 1: Nonlinear dynamics (due September Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.

More information

The Rela(onship Between Two Fast/ Slow Analysis Techniques for Burs(ng Oscilla(ons

The Rela(onship Between Two Fast/ Slow Analysis Techniques for Burs(ng Oscilla(ons The Rela(onship Between Two Fast/ Slow Analysis Techniques for Burs(ng Oscilla(ons Richard Bertram Department of Mathema(cs Florida State University Tallahassee, Florida Collaborators and Support Joël

More information

Local Phase Portrait of Nonlinear Systems Near Equilibria

Local Phase Portrait of Nonlinear Systems Near Equilibria Local Phase Portrait of Nonlinear Sstems Near Equilibria [1] Consider 1 = 6 1 1 3 1, = 3 1. ( ) (a) Find all equilibrium solutions of the sstem ( ). (b) For each equilibrium point, give the linear approimating

More information

Mathematical Models of Biological Systems

Mathematical Models of Biological Systems Mathematical Models of Biological Systems CH924 November 25, 2012 Course outline Syllabus: Week 6: Stability of first-order autonomous differential equation, genetic switch/clock, introduction to enzyme

More information

Modeling Microbial Populations in the Chemostat

Modeling Microbial Populations in the Chemostat Modeling Microbial Populations in the Chemostat Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Modeling Microbial Populations in the Chemostat MBI, June 3, 204 / 34 Outline Why

More information

F n = F n 1 + F n 2. F(z) = n= z z 2. (A.3)

F n = F n 1 + F n 2. F(z) = n= z z 2. (A.3) Appendix A MATTERS OF TECHNIQUE A.1 Transform Methods Laplace transforms for continuum systems Generating function for discrete systems This method is demonstrated for the Fibonacci sequence F n = F n

More information

Nonlinear Dynamics of Neural Firing

Nonlinear Dynamics of Neural Firing Nonlinear Dynamics of Neural Firing BENG/BGGN 260 Neurodynamics University of California, San Diego Week 3 BENG/BGGN 260 Neurodynamics (UCSD) Nonlinear Dynamics of Neural Firing Week 3 1 / 16 Reading Materials

More information

Continuous Threshold Policy Harvesting in Predator-Prey Models

Continuous Threshold Policy Harvesting in Predator-Prey Models Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July

More information

Math 5490 November 12, 2014

Math 5490 November 12, 2014 Math 5490 November 12, 2014 Topics in Applied Mathematics: Introduction to the Mathematics of Climate Mondays and Wednesdays 2:30 3:45 http://www.math.umn.edu/~mcgehee/teaching/math5490-2014-2fall/ Streaming

More information

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12 Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

More information

ODE, part 2. Dynamical systems, differential equations

ODE, part 2. Dynamical systems, differential equations ODE, part 2 Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Dynamical systems, differential equations Consider a system of n first order equations du dt = f(u, t),

More information

MATH 215/255 Solutions to Additional Practice Problems April dy dt

MATH 215/255 Solutions to Additional Practice Problems April dy dt . For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the

More information

Magnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation

Magnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation Center for Turbulence Research Annual Research Briefs 006 363 Magnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation By S. Fedotov AND S. Abarzhi 1. Motivation

More information

FAST AND HETEROCLINIC SOLUTIONS FOR A SECOND ORDER ODE

FAST AND HETEROCLINIC SOLUTIONS FOR A SECOND ORDER ODE 5-Oujda International Conference on Nonlinear Analysis. Electronic Journal of Differential Equations, Conference 14, 6, pp. 119 14. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu

More information

A review of stability and dynamical behaviors of differential equations:

A review of stability and dynamical behaviors of differential equations: A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary

More information

Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model

Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model Letter Forma, 15, 281 289, 2000 Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model Yasumasa NISHIURA 1 * and Daishin UEYAMA 2 1 Laboratory of Nonlinear Studies and Computations,

More information

1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model

1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model 1 Hodgkin-Huxley Theory of Nerve Membranes: The FitzHugh-Nagumo model Alan Hodgkin and Andrew Huxley developed the first quantitative model of the propagation of an electrical signal (the action potential)

More information

Math 312 Lecture Notes Linearization

Math 312 Lecture Notes Linearization Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior

More information

Qualitative Analysis of Tumor-Immune ODE System

Qualitative Analysis of Tumor-Immune ODE System of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System

More information

THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL

THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd

More information

1 The pendulum equation

1 The pendulum equation Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating

More information

Stability of Dynamical systems

Stability of Dynamical systems Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u

More information

Problem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form

Problem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form Problem set 6 Math 207A, Fall 2011 s 1 A two-dimensional gradient sstem has the form x t = W (x,, x t = W (x, where W (x, is a given function (a If W is a quadratic function W (x, = 1 2 ax2 + bx + 1 2

More information

Dispersion relations, stability and linearization

Dispersion relations, stability and linearization Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential

More information

Mixed mode oscillations as a mechanism for pseudo-plateau bursting

Mixed mode oscillations as a mechanism for pseudo-plateau bursting J Comput Neurosci (2010) 28:443 458 DOI 10.1007/s10827-010-0226-7 Mixed mode oscillations as a mechanism for pseudo-plateau bursting Theodore Vo Richard Bertram Joel Tabak Martin Wechselberger Received:

More information

Outline. Learning Objectives. References. Lecture 2: Second-order Systems

Outline. Learning Objectives. References. Lecture 2: Second-order Systems Outline Lecture 2: Second-order Systems! Techniques based on linear systems analysis! Phase-plane analysis! Example: Neanderthal / Early man competition! Hartman-Grobman theorem -- validity of linearizations!

More information

Dynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs

Dynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs Dynamics of a opulation Model Controlling the Spread of lague in rairie Dogs Catalin Georgescu The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

Phase Oscillators. and at r, Hence, the limit cycle at r = r is stable if and only if Λ (r ) < 0.

Phase Oscillators. and at r, Hence, the limit cycle at r = r is stable if and only if Λ (r ) < 0. 1 Phase Oscillators Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 2 Phase Oscillators Oscillations

More information

Lecture 15: Biological Waves

Lecture 15: Biological Waves Lecture 15: Biological Waves Jonathan A. Sherratt Contents 1 Wave Fronts I: Modelling Epidermal Wound Healing 2 1.1 Epidermal Wound Healing....................... 2 1.2 A Mathematical Model.........................

More information

MATH3203 Lecture 1 Mathematical Modelling and ODEs

MATH3203 Lecture 1 Mathematical Modelling and ODEs MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1

More information

Math 5490 November 5, 2014

Math 5490 November 5, 2014 Math 549 November 5, 214 Topics in Applied Mathematics: Introduction to the Mathematics of Climate Mondays and Wednesdays 2:3 3:45 http://www.math.umn.edu/~mcgehee/teaching/math549-214-2fall/ Streaming

More information

E209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions

E209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions E09A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions Michael Vitus Stanford University Winter 007 Problem : Planar phase portraits. Part a Figure : Problem a This phase portrait is correct.

More information

ASYMPTOTIC AND BIFURCATION ANALYSIS OF WAVE-PINNING IN A REACTION-DIFFUSION MODEL FOR CELL POLARIZATION

ASYMPTOTIC AND BIFURCATION ANALYSIS OF WAVE-PINNING IN A REACTION-DIFFUSION MODEL FOR CELL POLARIZATION ASYMPTOTIC AND BIFURCATION ANALYSIS OF WAVE-PINNING IN A REACTION-DIFFUSION MODEL FOR CELL POLARIZATION YOICHIRO MORI, ALEXANDRA JILKINE, AND LEAH EDELSTEIN-KESHET Abstract. We describe and analyze a bistable

More information

Gerardo Zavala. Math 388. Predator-Prey Models

Gerardo Zavala. Math 388. Predator-Prey Models Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked

More information

MathQuest: Differential Equations

MathQuest: Differential Equations MathQuest: Differential Equations Geometry of Systems 1. The differential equation d Y dt = A Y has two straight line solutions corresponding to [ ] [ ] 1 1 eigenvectors v 1 = and v 2 2 = that are shown

More information

Existence, Uniqueness Solution of a Modified. Predator-Prey Model

Existence, Uniqueness Solution of a Modified. Predator-Prey Model Nonlinear Analysis and Differential Equations, Vol. 4, 6, no. 4, 669-677 HIKARI Ltd, www.m-hikari.com https://doi.org/.988/nade.6.6974 Existence, Uniqueness Solution of a Modified Predator-Prey Model M.

More information

4.3 Intracellular calcium dynamics

4.3 Intracellular calcium dynamics coupled equations: dv dv 2 dv 3 dv 4 = i + I e A + g,2 (V 2 V ) = i 2 + g 2,3 (V 3 V 2 )+g 2, (V V 2 )+g 2,4 (V 4 V 2 ) = i 3 + g 3,2 (V 2 V 3 ) = i 4 + g 4,2 (V 2 V 4 ) This can be written in matrix/vector

More information

TMA 4195 Mathematical Modelling December 12, 2007 Solution with additional comments

TMA 4195 Mathematical Modelling December 12, 2007 Solution with additional comments Problem TMA 495 Mathematical Modelling December, 007 Solution with additional comments The air resistance, F, of a car depends on its length, L, cross sectional area, A, its speed relative to the air,

More information

Introduction and the Hodgkin-Huxley Model

Introduction and the Hodgkin-Huxley Model 1 Introduction and the Hodgkin-Huxley Model Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Reference:

More information

Introduction to Physiology IV - Calcium Dynamics

Introduction to Physiology IV - Calcium Dynamics Introduction to Physiology IV - Calcium Dynamics J. P. Keener Mathematics Department Introduction to Physiology IV - Calcium Dynamics p.1/26 Introduction Previous lectures emphasized the role of sodium

More information

Math 216 First Midterm 19 October, 2017

Math 216 First Midterm 19 October, 2017 Math 6 First Midterm 9 October, 7 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that

More information

Lotka Volterra Predator-Prey Model with a Predating Scavenger

Lotka Volterra Predator-Prey Model with a Predating Scavenger Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics

More information

MA 777: Topics in Mathematical Biology

MA 777: Topics in Mathematical Biology MA 777: Topics in Mathematical Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.math.uky.edu/~dmu228/ma777/ Spring 2018 David Murrugarra (University of Kentucky) Lecture

More information

ENGI Duffing s Equation Page 4.65

ENGI Duffing s Equation Page 4.65 ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing

More information

APPPHYS217 Tuesday 25 May 2010

APPPHYS217 Tuesday 25 May 2010 APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 14, 2019, at 08 30 12 30 Johanneberg Kristian

More information

A. Correct! Genetically a female is XX, and has 22 pairs of autosomes.

A. Correct! Genetically a female is XX, and has 22 pairs of autosomes. MCAT Biology - Problem Drill 08: Meiosis and Genetic Variability Question No. 1 of 10 1. A human female has pairs of autosomes and her sex chromosomes are. Question #01 (A) 22, XX. (B) 23, X. (C) 23, XX.

More information

Effects of transport memory and nonlinear damping in a generalized Fisher s equation

Effects of transport memory and nonlinear damping in a generalized Fisher s equation PHYSICAL REVIEW E, VOLUME 64, 066615 Effects of transport memory and nonlinear damping in a generalized Fisher s equation G. Abramson 1, * A. R. Bishop, and V. M. Kenkre 1 1 Center for Advanced Studies

More information

Nonlinear Control Lecture 2:Phase Plane Analysis

Nonlinear Control Lecture 2:Phase Plane Analysis Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53

More information

Problem Set Number 02, j/2.036j MIT (Fall 2018)

Problem Set Number 02, j/2.036j MIT (Fall 2018) Problem Set Number 0, 18.385j/.036j MIT (Fall 018) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139) September 6, 018 Due October 4, 018. Turn it in (by 3PM) at the Math. Problem Set

More information

LIMIT CYCLE OSCILLATORS

LIMIT CYCLE OSCILLATORS MCB 137 EXCITABLE & OSCILLATORY SYSTEMS WINTER 2008 LIMIT CYCLE OSCILLATORS The Fitzhugh-Nagumo Equations The best example of an excitable phenomenon is the firing of a nerve: according to the Hodgkin

More information

Differential Equation Types. Moritz Diehl

Differential Equation Types. Moritz Diehl Differential Equation Types Moritz Diehl Overview Ordinary Differential Equations (ODE) Differential Algebraic Equations (DAE) Partial Differential Equations (PDE) Delay Differential Equations (DDE) Ordinary

More information

Invasion slowing and pinning due to spatially heterogeneous selection

Invasion slowing and pinning due to spatially heterogeneous selection Invasion slowing and pinning due to spatially heterogeneous selection Georgetown University May 2016 Cane toad Evolving: leg length Purple loosestrife Evolving: flowering time, vegetative size at first

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for COMPUTATIONAL BIOLOGY A. COURSE CODES: FFR 110, FIM740GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for COMPUTATIONAL BIOLOGY A. COURSE CODES: FFR 110, FIM740GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for COMPUTATIONAL BIOLOGY A COURSE CODES: FFR 110, FIM740GU, PhD Time: Place: Teachers: Allowed material: Not allowed: June 8, 2018, at 08 30 12 30 Johanneberg Kristian

More information

Computational Neuroscience. Session 4-2

Computational Neuroscience. Session 4-2 Computational Neuroscience. Session 4-2 Dr. Marco A Roque Sol 06/21/2018 Two-Dimensional Two-Dimensional System In this section we will introduce methods of phase plane analysis of two-dimensional systems.

More information

Objective: Building a Bioreactor with intracellular reactions

Objective: Building a Bioreactor with intracellular reactions www.optience.com BioCell Model Objective: Building a Bioreactor with intracellular reactions This example shows a more advanced bioreactor model where both intracellular and extracellular concentrations

More information

Modeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool

Modeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht

More information

Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations

Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function

More information

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian

More information

PDE: The Method of Characteristics Page 1

PDE: The Method of Characteristics Page 1 PDE: The Method of Characteristics Page y u x (x, y) + x u y (x, y) =, () u(, y) = cos y 2. Solution The partial differential equation given can be rewritten as follows: u(x, y) y, x =, (2) where = / x,

More information

Exam in TMA4195 Mathematical Modeling Solutions

Exam in TMA4195 Mathematical Modeling Solutions Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Exam in TMA495 Mathematical Modeling 6..07 Solutions Problem a Here x, y are two populations varying with time

More information

A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term

A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science

More information

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems. Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics

More information

Introduction to First Order Equations Sections

Introduction to First Order Equations Sections A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Introduction to First Order Equations Sections 2.1-2.3 Dr. John Ehrke Department of Mathematics Fall 2012 Course Goals The

More information

Qualitative Analysis of Tumor-Immune ODE System

Qualitative Analysis of Tumor-Immune ODE System of Tumor-Immune ODE System LG de Pillis and AE Radunskaya August 15, 2002 This work was supported in part by a grant from the WM Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1 Simplified System of

More information

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium

More information

Multiscale Analysis of Many Particle Systems with Dynamical Control

Multiscale Analysis of Many Particle Systems with Dynamical Control Michael Herrmann Multiscale Analysis of Many Particle Systems with Dynamical Control joint work with Barbara Niethammer and Juan J.L. Velázquez Kinetic description of multiscale phenomena Archimedes Center

More information

Final 09/14/2017. Notes and electronic aids are not allowed. You must be seated in your assigned row for your exam to be valid.

Final 09/14/2017. Notes and electronic aids are not allowed. You must be seated in your assigned row for your exam to be valid. Final 09/4/207 Name: Problems -5 are each worth 8 points. Problem 6 is a bonus for up to 4 points. So a full score is 40 points and the max score is 44 points. The exam has 6 pages; make sure you have

More information

Lecture 18: Bistable Fronts PHYS 221A, Spring 2017

Lecture 18: Bistable Fronts PHYS 221A, Spring 2017 Lecture 18: Bistable Fronts PHYS 221A, Spring 2017 Lectures: P. H. Diamond Notes: Xiang Fan June 15, 2017 1 Introduction In the previous lectures, we learned about Turing Patterns. Turing Instability is

More information

Math 20D: Form B Final Exam Dec.11 (3:00pm-5:50pm), Show all of your work. No credit will be given for unsupported answers.

Math 20D: Form B Final Exam Dec.11 (3:00pm-5:50pm), Show all of your work. No credit will be given for unsupported answers. Turn off and put away your cell phone. No electronic devices during the exam. No books or other assistance during the exam. Show all of your work. No credit will be given for unsupported answers. Write

More information

Lecture 7. Heat Conduction. BENG 221: Mathematical Methods in Bioengineering. References

Lecture 7. Heat Conduction. BENG 221: Mathematical Methods in Bioengineering. References BENG 221: Mathematical Methods in Bioengineering Lecture 7 Heat Conduction References Haberman APDE, Ch. 1. http://en.wikipedia.org/wiki/heat_equation Lecture 5 BENG 221 M. Intaglietta Heat conduction

More information

Def. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)

Def. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable) Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information

Microevolution Changing Allele Frequencies

Microevolution Changing Allele Frequencies Microevolution Changing Allele Frequencies Evolution Evolution is defined as a change in the inherited characteristics of biological populations over successive generations. Microevolution involves the

More information

Local Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)

Local Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid) Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two

More information

Action Potentials and Synaptic Transmission Physics 171/271

Action Potentials and Synaptic Transmission Physics 171/271 Action Potentials and Synaptic Transmission Physics 171/271 Flavio Fröhlich (flavio@salk.edu) September 27, 2006 In this section, we consider two important aspects concerning the communication between

More information

2D-Volterra-Lotka Modeling For 2 Species

2D-Volterra-Lotka Modeling For 2 Species Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose

More information

Combustion fronts in porous media with two layers

Combustion fronts in porous media with two layers Combustion fronts in porous media with two layers layer 1 layer 2 Steve Schecter North Carolina State University (USA) Jesus Carlos da Mota Universidade Federal de Goiás (Brazil) 1 Overview Subject: Propagation

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

McGill University April 20, Advanced Calculus for Engineers

McGill University April 20, Advanced Calculus for Engineers McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student

More information

MS: Nonlinear Wave Propagation in Singular Perturbed Systems

MS: Nonlinear Wave Propagation in Singular Perturbed Systems MS: Nonlinear Wave Propagation in Singular Perturbed Systems P. van Heijster: Existence & stability of 2D localized structures in a 3-component model. Y. Nishiura: Rotational motion of traveling spots

More information