1 / 38

BME 6938 Neurodynamics

BME 6938 Neurodynamics. Instructor: Dr Sachin S Talathi. Recap . One Dimensional Flow: Bifurcations, Normal form, Stability One dimensional neuron model: Inward activation Inward Inactivation Outward Activation Outward Inactivation XPPAUTO for neuron modeling Auto

Download Presentation

BME 6938 Neurodynamics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. BME 6938Neurodynamics Instructor: Dr Sachin S Talathi

  2. Recap • One Dimensional Flow: Bifurcations, Normal form, Stability • One dimensional neuron model: • Inward activation • Inward Inactivation • Outward Activation • Outward Inactivation • XPPAUTO for neuron modeling • Auto • Bifurcation diagram, bistability, hysterisis

  3. Two Dimensional Dynamical Systems • Also called planar systems • f and g describe the evolution of two-dimensional state variables: x(t), y(t) • f and g are vector fields that describe how the trajectory will evolve

  4. Two dimensional vector fields:

  5. Nullclines • Nullclines are curves along which the vector fields are completely horizontal or completely vertical • The set of points where the vector field changes its horizontal direction defines the x-nullcline • x-nullcline partitions the phase space into two parts where either x increases or x decreases • x-nullcline is defined through • Similarly y-nullcline is defined through

  6. Example • Consider the following two dimensional ODE • Draw the nullclines for the dynamical system. Infer the stability of the fixed point from the phase plot by looking at the evolution of the vector field

  7. XPPAUTO to see Nullclines/Direction Fields of a Neuron Model V-Nullcline n-Nullcline

  8. Setting up XppAuto for TwoD-Neuron Model (Ex3.ode) • The ODE file: Copy the ode file from course website (Ex3.ode) • Neglect the comment in the ode file on type-1 and type-II dynamics for now. We will look into it later in the class. • For time being select the parameters for type-1 dynamics • In the nUmerics menu select Dt=0.1; Total Time=100; Bounds=1000 • Go to main menu and press (W)(W) to set X-axis:[0 100] and y-axis [-80 120] • Click (H)(C) to open a new plot window. Click (v)(2d) to set the axes: x-axis -80:100 y-axis:-.5:1

  9. We can get oscillations in Two-D (Limit cycles and ) • Limit cycles: Isolated periodic orbits in the phase space • Center: Continuum of periodic orbits in the phase space • Limit cycles are inherently a nonlinear phenomenon • Usually it is difficult to guess the existence of limit cycles by looking at the set of ODEs • We will learn more about limit cycles and their role in neuronal dynamics Simple example • Can be observed in linear systems • The amplitude of oscillations typically depend on initial conditions. Simple example:

  10. Relaxation oscillators: Intuitive way of thinking about the origin of spikes by neurons • Relaxation oscillators are special because of huge separation of time scales • Neurons also exhibit fast and slow time scales • Neurons are not relaxation oscillators, because mu is not small enough; • There is finite rise time for membrane potential • Time permitting we will consider relaxation oscillators in details when we study bursting neuron models

  11. Recap • Two Dimensional Flow • Direction Field; Nullclines • Mentioned limit cycles in passing (closed trajectories in phase space) • XPPAUTO to draw nullclines and direction fields.

  12. Phase portrait for two dimensional system: General Features A B C D Fixed points A, B and C satisfy D represents a closed trajectory; solution for which

  13. Two Dimensional Phase PortraitAdditonal Properties • Trajectories in phase portrait do not cross each other • (Consequence of existence and uniqueness property of solution for initial value problem in dynamical systems) • Fate of a trajectory enclosed in a close in a closed orbit • Poincare-Bendixon Theorem: If a trajectory is enclosed by a closed orbit and there is no fixed point inside the trajectory, then the trajectory will eventually approach the closed orbit • Also tells us that a closed orbit encloses a fixed-point (Another theorem, btw)

  14. Two Dimensional Linear systems A is constant matrix With initial conditions General solution is

  15. Simple example Uncoupled two-d system where Solution:

  16. Phase portrait-Bifurcation Trajectory approaches the stable fixed point in direction tangent to the slower axis

  17. General Case: Some Linear Algebra Diagonalize A: where =2x2 diagonal matrix OR

  18. What choice for P will diagonalize A? • Select P to be the eigen-vector of A that satisfies the following linear equation Eigen value of A Eigen vector of A Find eigen vectors and eigen values by solving the characteristics equation

  19. 3 possibilities OR

  20. Choice of for three cases considered We have: Case I.

  21. Example Eigen values are: Eigen vectors are: For a general invertible 2x2 matrix The solution to Easy to show

  22. Complex conjugate eigen-values • Case II • Claim: There exists P such that Let z be a complex eigenvector of A such that Let

  23. Example Eigen values are: Eigen vectors corresponding to is The solution to & Verify and show

  24. Degenerate eigenvalues • Case III. • Trivial case: for any choice of e1 and e2 • Nontrivial case: There exists matrix P=[v1 v2] such that Note:If we are willing to work with complex eigen vectors then it is possible to dia gonalize any matrix A with distinct eigen values

  25. Example Eigen values are: Lets choose Note: & Verify and show

  26. Intuitive way to think of degenerate case • The trajectories are almost trying to spiral into the fixed point (if stable) but do not quite make it • It represents a border line case between a node and spiral

  27. Classification Scheme for fixed points • We saw how to get explicit solution to any 2-dimensional linear ODE. • However if we are interested in global properties of the trajectories: explicit solutions are unnecessary • The eigen values of A give us all the information about the stability properties of the fixed points of the system and we can devise classification scheme based on the eigenvalues for the system

  28. Classification Scheme for fixed points • Case I: saddle node • Case II: • (a) either stable or unstable node • (b) either stable or unstable spiral • (c) border line between nodes and spirals • Case III: one of the eigenvalues is zero; no isolated fixed point but a series of fixed points (centers)

  29. Summary of the classification scheme centers Degenerate nodes

  30. Dynamics around fixed point based on previous analysis

  31. Linear Stability Analysis for Nonlinear Two-Dimensional System Evaluated at the fixed point (x*, y*)

  32. Few comments • The Linearized system does represent the dynamics of the nonlinear system locally around the fixed point correctly (Stable manifold theorem or the Hartman-Grobman Theorem) when the real part of eigenvalues are non-zero. • Fixed point in these case are referred to as the hyperbolic fixed points. The contribution from nonlinear higher order terms is negligible locally around hyperbolic fixed points. • Nonhyperbolic fixed points are those for which the real part of eigenvalue is zero. They are more sensitive to higher order nonlinear terms eg: Center’s Bifurcation points etc..

  33. Invariant Manifolds • The eigenvectors of A corresponding to the cases with non-degenerate eigenvalues considered earlier represent invariant manifolds for the dynamical system. Eg. Lets say the phase space is 2 dimensional made up of dynamical variables x and y. If initial condition is on x-axis and the flow as the system evolves remains on x-axis then x-axis is the invariant manifold of the dynamical system. In other words, orbits that start on the manifold remain in it.

  34. Examples of Invariant Manifold • Invariant manifolds of saddle (1-Dimensional manifolds) • Spirals (II-Dimensional manifolds)

  35. Stable and Unstable Manifolds of Saddle Unstable manifold: v1 Stable manifold: v2 Stable manifold of saddle is also referred to as the seperatrix; since it separates the phase plane into different regions of long term behavior

  36. Revisit a simple example We saw the phase space for this system earlier in our class. Lets Revisit and now Identify the stable and unstable manifolds of the saddle node (This time using XPPAUTO)

  37. Important of stable manifold of saddle in Neurodynamics Note: Threshold is not a single voltage value. It is a curve in Phase space; defined by the Stable manifold Of saddle

  38. Homoclinic and Heteroclinic Trajectories A trajectory is homoclinic if it originates from and terminates at the same equilibrium point A trajectory is heteroclinic if it originates at one equilibrium and terminates at a different equilibrium point

More Related