Entrainment and Chaos in the Hodgkin-Huxley Oscillator
|
|
- Ginger Stewart
- 6 years ago
- Views:
Transcription
1 Entrainment and Chaos in the Hodgkin-Huxley Oscillator Kevin K. Lin klin Courant Institute, New York University Mostly Biomath p.1/42
2 Overview (1) Goal: Show that the Hodgkin-Huxley neuron model, driven by a periodic impulse train, can exhibit entrainment, transient chaos, and fully chaotic behavior. Mostly Biomath p.2/42
3 Overview (2) Why? 1. Suggested by general, rigorous perturbation theory of kicked oscillators (Qiudong Wang & Lai-Sang Young). 2. Hodgkin-Huxley is a paradigm for excitable biological systems where pulse-like inputs and outputs are natural. Mostly Biomath p.3/42
4 Overview (3) Results: 1. Entrainment and chaos are readily observable in the pulse-driven Hodgkin-Huxley system. 2. The pulse-driven Hodgkin-Huxley system prefers entrainment. 3. Strong expansion is caused by invariant structures. Mostly Biomath p.4/42
5 Outline Classical Hodgkin-Huxley neuron model Kicked nonlinear oscillators & Wang-Young theory Pulse-driven Hodgkin-Huxley neuron model Mostly Biomath p.5/42
6 Squid giant axon Mostly Biomath p.6/42
7 Schematic (rest state) References: Abbott and Dayan, Theoretical Neuroscience, MIT Press 2001 Cronin, Mathematical Aspects of Hodgkin-Huxley Neural Theory, Cambridge University Press 1987 Mostly Biomath p.7/42
8 Hodgkin-Huxley equations v = C 1 [ I ḡ K n 4 (v v K ) ḡ Na m 3 h(v v Na ) ḡ leak (v v leak ) ] ṁ = αm(v)(1 m) βm(v)m ṅ = αn(v)(1 n) βn(v)n ḣ = α h (v)(1 h) β h (v)h αm(v) = Ψ ( ) v+25 10, β m(v) = 4 exp (v/18), αn(v) = 0.1Ψ ( ) v+10 10, β n(v) = exp (v/80), 1 α h (v) = 7 exp (v/20), β h (v) = 1+exp( v ), v Ψ(v) = exp(v) 1. Mostly Biomath p.8/42
9 Equivalent circuit [ ] v = C 1 I ḡ K n 4 (v v K ) ḡ Na m 3 h(v v Na ) ḡ leak (v v leak ) Mostly Biomath p.9/42
10 Parameters In this study: All parameters take on original squid values except the injected current I This guarantees stable oscillations Mostly Biomath p.10/42
11 Parameters (cont d) v Na = 115 mv, ḡ Na = 120 mω 1 /cm 2, v K = +12 mv, ḡ K = 36 mω 1 /cm 2, v leak = mv, ḡ leak = 0.3 mω 1 /cm 2, C = 1 µf/cm 2, I = Mostly Biomath p.11/42
12 Parameters (cont d) Stable fixed point Unstable fixed point V Unstable cycle -80 Limit cycle I Mostly Biomath p.12/42
13 Dynamics without kicks v -4 n t v -2 Mostly Biomath p.13/42
14 Outline Classical Hodgkin-Huxley neuron model Kicked nonlinear oscillators & Wang-Young theory Pulse-driven Hodgkin-Huxley neuron model Mostly Biomath p.14/42
15 Kicked oscillators A stable, nonlinear oscillator is a flow with a limit cycle γ (period=t γ ) and basin of attraction U. A kick instantaneously changes the system s state: Mostly Biomath p.15/42
16 Examples of kicked oscillators Circadian rhythm, phase reset experiments (Winfree). Possible approach to disrupting synchronous firing of neuron (Tass). Mostly Biomath p.16/42
17 Simple Example 1.0 y x ṙ = (µ αr 2 )r+ 1 2 sin(3θ) δ(t nt) n Z θ = ω + βr 2 Mostly Biomath p.17/42
18 Effect of Kick-and-Flow on Phase Space t = 0 t = 0 (after kick) t = t = 2 t = 2 (after kick) t = 4 Mostly Biomath p.18/42
19 Discrete time map Define F T : R 4 R 4 by F T (x) = Φ T (K(x)), where K(x) represents kicks Φ T (x) = flow map T = period of kicks. Continuous time Discrete time Entrainment F T has sinks Chaos F T chaotic Mostly Biomath p.19/42
20 Reduction to 1-D Wang and Young start with F T and 1. Reduces from the map F T on R n to a circle map f T : lim F T+nT γ (x) n induces a map f T on γ S 1. We refer to f T as the singular limit or the phase resetting curve. 2. Analyze f T and infer properties of F T. Mostly Biomath p.20/42
21 Wang-Young Conditions If 1. Kicks do not send limit cycle to bad places, i.e. K(γ) does not go outside the basin of γ 2. Kicks are in the right directions (e.g. not along W ss (x)) to take advantage of shear Mostly Biomath p.21/42
22 Wang-Young Consequences Then for different kick amplitude A & kick period T the discrete-time system F T can have 1. Rotation-like behavior (small A) 2. Sinks and sources (for all A large enough) 3. Transient chaos / horseshoes (for large A & T) 4. Strange attractor & chaos (for large A & T 1) Mostly Biomath p.22/42
23 Wang-Young Theory (cont d) Notes: 1. The conditions are satisfied for the simple example. 2. For Hodgkin-Huxley there is not too much choice if we want to stay close to physical interpretation of model. Mostly Biomath p.23/42
24 Lyapunov exponents The largest Lyapunov exponent λ of F T is useful for distinguishing different scenarios numerically: Rotations λ = 0 Sinks λ < 0 Chaos λ > 0 Mostly Biomath p.24/42
25 Outline Classical Hodgkin-Huxley neuron model Kicked nonlinear oscillators & Wang-Young theory Pulse-driven Hodgkin-Huxley neuron model Mostly Biomath p.25/42
26 Hodgkin-Huxley equations v = C 1 [ I ḡ K n 4 (v v K ) ḡ Na m 3 h(v v Na ) ḡ leak (v v leak ) ] +A n Z δ (t nt) ṁ = αm(v)(1 m) βm(v)m ṅ = αn(v)(1 n) βn(v)n ḣ = α h (v)(1 h) β h (v)h Prior work: Winfree, Best on null space and degree-change. Mostly Biomath p.26/42
27 Entrainment A = 10, T = v t v t Mostly Biomath p.27/42
28 Entrainment (cont d) -2 A = 10, T = 81.0 v t v1(t)-v2(t) t Mostly Biomath p.28/42
29 Entrainment (cont d) Time-T map: F T = Φ T K 2 v1(n)-v2(n) Multiple of T (n) 2 Mostly Biomath p.29/42
30 Chaos v A = 10, T = t 2.0 Log10(dist) t Mostly Biomath p.30/42
31 λ (F T ) versus T Largest Lyapunov exponent of F_T T / T_gamma Mostly Biomath p.31/42
32 λ (F T ) versus A ROTATIONS 0.6 SINKS CHAOS UNKNOWN Drive amplitude Chaos: Prob( ˆλ > 3ǫ) Rotations: Prob( ˆλ < ǫ/3) Sinks: Prob( ˆλ < 3ǫ) Unknown: everything else Mostly Biomath p.32/42
33 Phase resetting curves ( f T ) A=5, T=101.5 A=10, T= A=20, T= Mostly Biomath p.33/42
34 Plateau and fixed points The first return map R ft to the interval [4, 10] (enclosing the plateau), for A = 10 and T = Mostly Biomath p.34/42
35 Plateau and fixed points Drive amplitude A Probability of sink near plateau 5 41% 10 58% 20 68% 30 76% Mostly Biomath p.35/42
36 Zooming into the kink Mostly Biomath p.36/42
37 Why the kink? h 0.34 n v m But Hodgkin-Huxley lives in R 4... Mostly Biomath p.37/42
38 Why the kink? (cont d) Approaching critical A crit : A=13.58 A= A= A= A=13.59 A= Mostly Biomath p.38/42
39 Why the plateau? Graph of f T for A = 10, around plateau f_t phase phase Black: log 10 f Blue: log 10 (Ess(K(γ(θ)), γ(θ))) Mostly Biomath p.39/42
40 Finding horseshoes Horseshoes can produce transient chaos: A = 10, T = Mostly Biomath p.40/42
41 Summary Can observe entrainment and chaos in the pulse-driven Hodgkin-Huxley neuron model. The pulse-driven Hodgkin-Huxley model prefers entrainment. This can be explained. Complex phase response can arise from kicks going near invariant structures. Mostly Biomath p.41/42
42 References 1. Eric N. Best, Null space in the Hodgkin-Huxley equations, Biophys. J. 27 (1979) 2. Kevin K. Lin, Entrainment and chaos in the Hodgkin-Huxley oscillator, in preparation 3. Qiudong Wang and Lai-Sang Young, Strange attractors in periodically-kicked limit cycles and Hopf bifurcations, Comm. Math. Phys. 240 (2003) 4. Arthur Winfree, The Geometry of Biological Time, 2nd Edition, Springer-Verlag (2000) Acknowledgements I am grateful to Lai-Sang Young for her help with this work, and to Eric Brown, Adi Rangan, Alex Barnett, and Toufic Suidan for critical comments. Many thanks to Charlie Peskin for the invitation. This work is supported by the National Science Foundation. Mostly Biomath p.42/42
MATH 3104: THE HODGKIN-HUXLEY EQUATIONS
MATH 3104: THE HODGKIN-HUXLEY EQUATIONS Parallel conductance model A/Prof Geoffrey Goodhill, Semester 1, 2009 So far we have modelled neuronal membranes by just one resistance (conductance) variable. We
More informationElectrophysiology of the neuron
School of Mathematical Sciences G4TNS Theoretical Neuroscience Electrophysiology of the neuron Electrophysiology is the study of ionic currents and electrical activity in cells and tissues. The work of
More informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationDynamical 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 informationSynchronization and Phase Oscillators
1 Synchronization and Phase Oscillators Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 Synchronization
More informationDynamical systems in neuroscience. Pacific Northwest Computational Neuroscience Connection October 1-2, 2010
Dynamical systems in neuroscience Pacific Northwest Computational Neuroscience Connection October 1-2, 2010 What do I mean by a dynamical system? Set of state variables Law that governs evolution of state
More informationPhase Response Curves, Delays and Synchronization in Matlab
Phase Response Curves, Delays and Synchronization in Matlab W. Govaerts and B. Sautois Department of Applied Mathematics and Computer Science, Ghent University, Krijgslaan 281-S9, B-9000 Ghent, Belgium
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best Koch, Chapters 6, 8, 9 See also Hodgkin and Huxley, J. Physiol. 117:500-544 (1952. (the source Clay, J. Neurophysiol. 80:903-913
More informationNeural Modeling and Computational Neuroscience. Claudio Gallicchio
Neural Modeling and Computational Neuroscience Claudio Gallicchio 1 Neuroscience modeling 2 Introduction to basic aspects of brain computation Introduction to neurophysiology Neural modeling: Elements
More informationEntrainment Alex Bowie April 7, 2004
Entrainment Alex Bowie April 7, 2004 Abstract The driven Van der Pol oscillator displays entrainment, quasiperiodicity, and chaos. The characteristics of these different modes are discussed as well as
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
More informationVoltage-clamp and Hodgkin-Huxley models
Voltage-clamp and Hodgkin-Huxley models Read: Hille, Chapters 2-5 (best) Koch, Chapters 6, 8, 9 See also Clay, J. Neurophysiol. 80:903-913 (1998) (for a recent version of the HH squid axon model) Rothman
More informationStrange Nonchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley Neuron
Journal of the Korean Physical Society, Vol. 57, o. 1, July 2010, pp. 23 29 Strange onchaotic Spiking in the Quasiperiodically-forced Hodgkin-Huxley euron Woochang Lim and Sang-Yoon Kim Department of Physics,
More informationBursting and Chaotic Activities in the Nonlinear Dynamics of FitzHugh-Rinzel Neuron Model
Bursting and Chaotic Activities in the Nonlinear Dynamics of FitzHugh-Rinzel Neuron Model Abhishek Yadav *#, Anurag Kumar Swami *, Ajay Srivastava * * Department of Electrical Engineering, College of Technology,
More informationSingle neuron models. L. Pezard Aix-Marseille University
Single neuron models L. Pezard Aix-Marseille University Biophysics Biological neuron Biophysics Ionic currents Passive properties Active properties Typology of models Compartmental models Differential
More informationChaos in the Hodgkin Huxley Model
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 1, No. 1, pp. 105 114 c 2002 Society for Industrial and Applied Mathematics Chaos in the Hodgkin Huxley Model John Guckenheimer and Ricardo A. Oliva Abstract. The
More informationPhase 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 informationDynamical phase transitions in periodically driven model neurons
Dynamical phase transitions in periodically driven model neurons Jan R. Engelbrecht 1 and Renato Mirollo 2 1 Department of Physics, Boston College, Chestnut Hill, Massachusetts 02467, USA 2 Department
More informationCOMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL
More informationDynamical modelling of systems of coupled oscillators
Dynamical modelling of systems of coupled oscillators Mathematical Neuroscience Network Training Workshop Edinburgh Peter Ashwin University of Exeter 22nd March 2009 Peter Ashwin (University of Exeter)
More informationMATH 723 Mathematical Neuroscience Spring 2008 Instructor: Georgi Medvedev
MATH 723 Mathematical Neuroscience Spring 28 Instructor: Georgi Medvedev 2 Lecture 2. Approximate systems. 2. Reduction of the HH model to a 2D system. The original HH system consists of 4 differential
More informationModeling Action Potentials in Cell Processes
Modeling Action Potentials in Cell Processes Chelsi Pinkett, Jackie Chism, Kenneth Anderson, Paul Klockenkemper, Christopher Smith, Quarail Hale Tennessee State University Action Potential Models Chelsi
More informationLearning Cycle Linear Hybrid Automata for Excitable Cells
Learning Cycle Linear Hybrid Automata for Excitable Cells Sayan Mitra Joint work with Radu Grosu, Pei Ye, Emilia Entcheva, I V Ramakrishnan, and Scott Smolka HSCC 2007 Pisa, Italy Excitable Cells Outline
More informationShear-Induced Chaos. Kevin K. Lin and Lai-Sang Young Courant Institute of Mathematical Sciences New York University.
Shear-Induced Chaos Kevin K. Lin and Lai-Sang Young Courant Institute of Mathematical Sciences New York University February, Abstract Guided by a geometric understanding developed in earlier works of Wang
More informationThe Spike Response Model: A Framework to Predict Neuronal Spike Trains
The Spike Response Model: A Framework to Predict Neuronal Spike Trains Renaud Jolivet, Timothy J. Lewis 2, and Wulfram Gerstner Laboratory of Computational Neuroscience, Swiss Federal Institute of Technology
More informationNeuroscience applications: isochrons and isostables. Alexandre Mauroy (joint work with I. Mezic)
Neuroscience applications: isochrons and isostables Alexandre Mauroy (joint work with I. Mezic) Outline Isochrons and phase reduction of neurons Koopman operator and isochrons Isostables of excitable systems
More informationLinearization of F-I Curves by Adaptation
LETTER Communicated by Laurence Abbott Linearization of F-I Curves by Adaptation Bard Ermentrout Department of Mathematics, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A. We show that negative
More informationPseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics
J Comput Neurosci (1) 8:7 66 DOI 1.17/s187-9 Pseudo-s and predictability of Hodgkin-Huxley neuronal network dynamics Yi Sun Douglas Zhou Aaditya V. Rangan David Cai Received: 16 May 9 / Revised: 3 November
More informationConductance-Based Integrate-and-Fire Models
NOTE Communicated by Michael Hines Conductance-Based Integrate-and-Fire Models Alain Destexhe Department of Physiology, Laval University School of Medicine, Québec, G1K 7P4, Canada A conductance-based
More informationPhysics 235 Chapter 4. Chapter 4 Non-Linear Oscillations and Chaos
Chapter 4 Non-Linear Oscillations and Chaos Non-Linear Differential Equations Up to now we have considered differential equations with terms that are proportional to the acceleration, the velocity, and
More information1 Introduction and neurophysiology
Dynamics of Continuous, Discrete and Impulsive Systems Series B: Algorithms and Applications 16 (2009) 535-549 Copyright c 2009 Watam Press http://www.watam.org ASYMPTOTIC DYNAMICS OF THE SLOW-FAST HINDMARSH-ROSE
More informationOn Parameter Estimation for Neuron Models
On Parameter Estimation for Neuron Models Abhijit Biswas Department of Mathematics Temple University November 30th, 2017 Abhijit Biswas (Temple University) On Parameter Estimation for Neuron Models November
More informationShear-Induced Chaos. Kevin K. Lin and Lai-Sang Young Courant Institute of Mathematical Sciences New York University. May 22, 2007
Shear-Induced Chaos arxiv:75.394v [math.ds] 3 May 7 Kevin K. Lin and Lai-Sang Young Courant Institute of Mathematical Sciences New York University May, 7 Abstract Guided by a geometric understanding developed
More informationWHAT IS A CHAOTIC ATTRACTOR?
WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties
More informationThreshold, excitability and isochrones in the Bonhoeffer van der Pol system
CHAOS VOLUME 9, NUMBER 4 DECEMBER 1999 Threshold, excitability and isochrones in the Bonhoeffer van der Pol system A. Rabinovitch and I. Rogachevskii Physics Department and Department of Mechanical Engineering,
More informationStochastic Oscillator Death in Globally Coupled Neural Systems
Journal of the Korean Physical Society, Vol. 52, No. 6, June 2008, pp. 19131917 Stochastic Oscillator Death in Globally Coupled Neural Systems Woochang Lim and Sang-Yoon Kim y Department of Physics, Kangwon
More informationNonlinear Dynamics: Synchronisation
Nonlinear Dynamics: Synchronisation Bristol Centre for Complexity Sciences Ian Ross BRIDGE, School of Geographical Sciences, University of Bristol October 19, 2007 1 / 16 I: Introduction 2 / 16 I: Fireflies
More informationChaotic Vibrations. An Introduction for Applied Scientists and Engineers
Chaotic Vibrations An Introduction for Applied Scientists and Engineers FRANCIS C. MOON Theoretical and Applied Mechanics Cornell University Ithaca, New York A WILEY-INTERSCIENCE PUBLICATION JOHN WILEY
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More informationDynamics of Periodically Perturbed Homoclinic Solutions
Dynamics of Periodically Perturbed Homoclinic Solutions Qiudong Wang University of Arizona Part I Structure of Homoclinic Tangles Unperturbed equation Let (x,y) R 2 be the phase variables. dx dt = αx+f(x,y),
More informationProblem 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 informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationFrom neuronal oscillations to complexity
1/39 The Fourth International Workshop on Advanced Computation for Engineering Applications (ACEA 2008) MACIS 2 Al-Balqa Applied University, Salt, Jordan Corson Nathalie, Aziz Alaoui M.A. University of
More informationOn the Dynamics of Delayed Neural Feedback Loops. Sebastian Brandt Department of Physics, Washington University in St. Louis
On the Dynamics of Delayed Neural Feedback Loops Sebastian Brandt Department of Physics, Washington University in St. Louis Overview of Dissertation Chapter 2: S. F. Brandt, A. Pelster, and R. Wessel,
More informationIntegral Equation and PDE Neuronal Models:
Integral Equation and PDE Neuronal Models: 1952-2015 W. C. Troy Integral Equation and PDE Neuronal Models: 1952-2015 p.1/51 Part 1. 1D Single Neuron Models. (i) Hodgkin-Huxley (1952), Fitzhugh-Nagumo (1961)
More informationNonlinear systems, chaos and control in Engineering
Nonlinear systems, chaos and control in Engineering Module 1 block 3 One-dimensional nonlinear systems Cristina Masoller Cristina.masoller@upc.edu http://www.fisica.edu.uy/~cris/ Schedule Flows on the
More informationTopics in Neurophysics
Topics in Neurophysics Alex Loebel, Martin Stemmler and Anderas Herz Exercise 2 Solution (1) The Hodgkin Huxley Model The goal of this exercise is to simulate the action potential according to the model
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationComputational 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 informationSingle-Cell and Mean Field Neural Models
1 Single-Cell and Mean Field Neural Models Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 The neuron
More informationChimera states in networks of biological neurons and coupled damped pendulums
in neural models in networks of pendulum-like elements in networks of biological neurons and coupled damped pendulums J. Hizanidis 1, V. Kanas 2, A. Bezerianos 3, and T. Bountis 4 1 National Center for
More informationSingle-Compartment Neural Models
Single-Compartment Neural Models BENG/BGGN 260 Neurodynamics University of California, San Diego Week 2 BENG/BGGN 260 Neurodynamics (UCSD) Single-Compartment Neural Models Week 2 1 / 18 Reading Materials
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary Figure S1. Pulses >3mJ reduce membrane resistance in HEK cells. Reversal potentials in a representative cell for IR-induced currents with laser pulses of 0.74 to
More informationCanonical Neural Models 1
Canonical Neural Models 1 Frank Hoppensteadt 1 and Eugene zhikevich 2 ntroduction Mathematical modeling is a powerful tool in studying fundamental principles of information processing in the brain. Unfortunately,
More informationSynchrony in Stochastic Pulse-coupled Neuronal Network Models
Synchrony in Stochastic Pulse-coupled Neuronal Network Models Katie Newhall Gregor Kovačič and Peter Kramer Aaditya Rangan and David Cai 2 Rensselaer Polytechnic Institute, Troy, New York 2 Courant Institute,
More informationCHAPTER 1. Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting
CHAPTER 1 Bifurcations in the Fast Dynamics of Neurons: Implications for Bursting John Guckenheimer Mathematics Department, Cornell University Ithaca, NY 14853 Joseph H. Tien Center for Applied Mathematics,
More informationInvariant manifolds of the Bonhoeffer-van der Pol oscillator
Invariant manifolds of the Bonhoeffer-van der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
More informationarxiv:nlin/ v1 [nlin.cd] 4 Oct 2005
Synchronization of Coupled Chaotic Dynamics on Networks R. E. Amritkar and Sarika Jalan Physical Research Laboratory, Navrangpura, Ahmedabad 380 009, India. arxiv:nlin/0510008v1 [nlin.cd] 4 Oct 2005 Abstract
More informationNeural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback
Neural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback Gautam C Sethia and Abhijit Sen Institute for Plasma Research, Bhat, Gandhinagar 382 428, INDIA Motivation Neural Excitability
More informationNeuronal Dynamics: Computational Neuroscience of Single Neurons
Week 4 part 5: Nonlinear Integrate-and-Fire Model 4.1 From Hodgkin-Huxley to 2D Neuronal Dynamics: Computational Neuroscience of Single Neurons Week 4 Recing detail: Two-dimensional neuron models Wulfram
More informationSaturation of Information Exchange in Locally Connected Pulse-Coupled Oscillators
Saturation of Information Exchange in Locally Connected Pulse-Coupled Oscillators Will Wagstaff School of Computer Science, Georgia Institute of Technology, Atlanta, Georgia 30332, USA (Dated: 13 December
More informationBME 5742 Biosystems Modeling and Control
BME 5742 Biosystems Modeling and Control Hodgkin-Huxley Model for Nerve Cell Action Potential Part 1 Dr. Zvi Roth (FAU) 1 References Hoppensteadt-Peskin Ch. 3 for all the mathematics. Cooper s The Cell
More informationxt+1 = 1 ax 2 t + y t y t+1 = bx t (1)
Exercise 2.2: Hénon map In Numerical study of quadratic area-preserving mappings (Commun. Math. Phys. 50, 69-77, 1976), the French astronomer Michel Hénon proposed the following map as a model of the Poincaré
More informationDynamical Systems in Neuroscience: The Geometry of Excitability and Bursting
Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting Eugene M. Izhikevich The MIT Press Cambridge, Massachusetts London, England Contents Preface xv 1 Introduction 1 1.1 Neurons
More informationNumerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a Mathematical Model of Snail Neurons
International Journal of Theoretical and Mathematical Physics 2015, 5(5): 145-150 DOI: 10.5923/j.ijtmp.20150505.08 Numerical Simulation of Bistability between Regular Bursting and Chaotic Spiking in a
More informationA MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS
International Journal of Bifurcation and Chaos, Vol. 18, No. 5 (2008) 1567 1577 c World Scientific Publishing Company A MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS ZERAOULIA ELHADJ Department
More informationHandling of Chaos in Two Dimensional Discrete Maps
Handling of Chaos in Two Dimensional Discrete Maps Anil Kumar Jain Assistant Professor, Department of Mathematics Barama College, Barama, Assam, Pincode-781346, India (Email: jainanil965@gmail.com) Abstract:
More informationarxiv: v1 [nlin.cd] 20 Jul 2010
Invariant manifolds of the Bonhoeffer-van der Pol oscillator arxiv:1007.3375v1 [nlin.cd] 20 Jul 2010 R. Benítez 1, V. J. Bolós 2 1 Departamento de Matemáticas, Centro Universitario de Plasencia, Universidad
More informationExperimental and Simulated Chaotic RLD Circuit Analysis with the Use of Lorenz Maps
Experimental and Simulated Chaotic RLD Circuit Analysis with the Use of Lorenz Maps N.A. Gerodimos, P.A. Daltzis, M.P. Hanias, H.E. Nistazakis, and G.S. Tombras Abstract In this work, the method Ed. Lorenz
More informationDESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS
Letters International Journal of Bifurcation and Chaos, Vol. 8, No. 8 (1998) 1733 1738 c World Scientific Publishing Company DESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS I. P.
More informationOn the Response of Neurons to Sinusoidal Current Stimuli: Phase Response Curves and Phase-Locking
On the Response of Neurons to Sinusoidal Current Stimuli: Phase Response Curves and Phase-Locking Michael J. Schaus and Jeff Moehlis Abstract A powerful technique for analyzing mathematical models for
More informationExample Chaotic Maps (that you can analyze)
Example Chaotic Maps (that you can analyze) Reading for this lecture: NDAC, Sections.5-.7. Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January
More informationAN ELECTRIC circuit containing a switch controlled by
878 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 46, NO. 7, JULY 1999 Bifurcation of Switched Nonlinear Dynamical Systems Takuji Kousaka, Member, IEEE, Tetsushi
More informationMath 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction
Math 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction Junping Shi College of William and Mary November 8, 2018 Neuron Neurons Neurons are cells in the brain and other subsystems
More informationIs the Hénon map chaotic
Is the Hénon map chaotic Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology, Poland, galias@agh.edu.pl International Workshop on Complex Networks and Applications
More informationThe selection of mixed-mode oscillations in a Hodgkin-Huxley model with multiple timescales
CHAOS 18, 015105 2008 The selection of mixed-mode oscillations in a Hodgkin-Huxley model with multiple timescales Jonathan Rubin Department of Mathematics, University of Pittsburgh, Pittsburgh, Pennsylvania
More informationEffect of various periodic forces on Duffing oscillator
PRAMANA c Indian Academy of Sciences Vol. 67, No. 2 journal of August 2006 physics pp. 351 356 Effect of various periodic forces on Duffing oscillator V RAVICHANDRAN 1, V CHINNATHAMBI 1, and S RAJASEKAR
More information3 Action Potentials - Brutal Approximations
Physics 172/278 - David Kleinfeld - Fall 2004; Revised Winter 2015 3 Action Potentials - Brutal Approximations The Hodgkin-Huxley equations for the behavior of the action potential in squid, and similar
More informationTitle. Author(s)Fujii, Hiroshi; Tsuda, Ichiro. CitationNeurocomputing, 58-60: Issue Date Doc URL. Type.
Title Neocortical gap junction-coupled interneuron systems exhibiting transient synchrony Author(s)Fujii, Hiroshi; Tsuda, Ichiro CitationNeurocomputing, 58-60: 151-157 Issue Date 2004-06 Doc URL http://hdl.handle.net/2115/8488
More informationTopological analysis of chaos in neural spike train bursts
CHAOS VOLUME 9, NUMBER 3 SEPTEMBER 1999 Topological analysis of chaos in neural spike train bursts R. Gilmore Department of Physics and Atmospheric Science, Drexel University, Philadelphia, Pennsylvania
More informationChaos synchronization of complex Rössler system
Appl. Math. Inf. Sci. 7, No. 4, 1415-1420 (2013) 1415 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/070420 Chaos synchronization of complex Rössler
More informationCHALMERS, 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: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationChapter 3. Gumowski-Mira Map. 3.1 Introduction
Chapter 3 Gumowski-Mira Map 3.1 Introduction Non linear recurrence relations model many real world systems and help in analysing their possible asymptotic behaviour as the parameters are varied [17]. Here
More informationBursting Oscillations of Neurons and Synchronization
Bursting Oscillations of Neurons and Synchronization Milan Stork Applied Electronics and Telecommunications, Faculty of Electrical Engineering/RICE University of West Bohemia, CZ Univerzitni 8, 3064 Plzen
More informationRecent new examples of hidden attractors
Eur. Phys. J. Special Topics 224, 1469 1476 (2015) EDP Sciences, Springer-Verlag 2015 DOI: 10.1140/epjst/e2015-02472-1 THE EUROPEAN PHYSICAL JOURNAL SPECIAL TOPICS Review Recent new examples of hidden
More informationSTUDENT PAPER. Santiago Santana University of Illinois, Urbana-Champaign Blue Waters Education Program 736 S. Lombard Oak Park IL, 60304
STUDENT PAPER Differences between Stochastic and Deterministic Modeling in Real World Systems using the Action Potential of Nerves. Santiago Santana University of Illinois, Urbana-Champaign Blue Waters
More informationChaotifying 2-D piecewise linear maps via a piecewise linear controller function
Chaotifying 2-D piecewise linear maps via a piecewise linear controller function Zeraoulia Elhadj 1,J.C.Sprott 2 1 Department of Mathematics, University of Tébéssa, (12000), Algeria. E-mail: zeraoulia@mail.univ-tebessa.dz
More informationStabilization of Hyperbolic Chaos by the Pyragas Method
Journal of Mathematics and System Science 4 (014) 755-76 D DAVID PUBLISHING Stabilization of Hyperbolic Chaos by the Pyragas Method Sergey Belyakin, Arsen Dzanoev, Sergey Kuznetsov Physics Faculty, Moscow
More informationObservable Chaos in Switch-Controlled Chua s Circuit
Observable Chaos in Switch-Controlled Chua s Circuit Ali Oksasoglu, Dongsheng Ma, and Qiudong Wang, This research is partially supported by a grant from NSF. A. Oksasoglu is with Honeywell, Inc., Tucson,
More informationPhase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable
Printed from the Mathematica Help Browser 1 1 of 10 Phase Locking A neuron phase-locks to a periodic input it spikes at a fixed delay [Izhikevich07]. The PRC's amplitude determines which frequencies a
More informationAll-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model
All-or-None Principle and Weakness of Hodgkin-Huxley Mathematical Model S. A. Sadegh Zadeh, C. Kambhampati International Science Index, Mathematical and Computational Sciences waset.org/publication/10008281
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationA Highly Chaotic Attractor for a Dual-Channel Single-Attractor, Private Communication System
A Highly Chaotic Attractor for a Dual-Channel Single-Attractor, Private Communication System Banlue Srisuchinwong and Buncha Munmuangsaen Sirindhorn International Institute of Technology, Thammasat University
More informationLecture 10 : Neuronal Dynamics. Eileen Nugent
Lecture 10 : Neuronal Dynamics Eileen Nugent Origin of the Cells Resting Membrane Potential: Nernst Equation, Donnan Equilbrium Action Potentials in the Nervous System Equivalent Electrical Circuits and
More informationThis script will produce a series of pulses of amplitude 40 na, duration 1ms, recurring every 50 ms.
9.16 Problem Set #4 In the final problem set you will combine the pieces of knowledge gained in the previous assignments to build a full-blown model of a plastic synapse. You will investigate the effects
More informationNonlinear 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