Mixed Mode Oscillations as a Mechanism for Pseudo-Plateau Bursting

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

PURE MATHEMATICS A-LEVEL PAPER 1

1985 AP Calculus BC: Section I

ln x = n e = 20 (nearest integer)

Mixed mode oscillations as a mechanism for pseudo-plateau bursting

Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells

A Mathematical Study of Electrical Bursting in Pituitary Cells

ECE594I Notes set 6: Thermal Noise

LECTURE 13 Filling the bands. Occupancy of Available Energy Levels

Bifurcation Theory. , a stationary point, depends on the value of α. At certain values

Washington State University

MONTGOMERY COLLEGE Department of Mathematics Rockville Campus. 6x dx a. b. cos 2x dx ( ) 7. arctan x dx e. cos 2x dx. 2 cos3x dx

2008 AP Calculus BC Multiple Choice Exam

H2 Mathematics Arithmetic & Geometric Series ( )

(Reference: sections in Silberberg 5 th ed.)

EE 232 Lightwave Devices Lecture 3: Basic Semiconductor Physics and Optical Processes. Optical Properties of Semiconductors

Restricted Factorial And A Remark On The Reduced Residue Classes

Time : 1 hr. Test Paper 08 Date 04/01/15 Batch - R Marks : 120

Probability & Statistics,

Option 3. b) xe dx = and therefore the series is convergent. 12 a) Divergent b) Convergent Proof 15 For. p = 1 1so the series diverges.

CDS 101: Lecture 5.1 Reachability and State Space Feedback

Chapter 2 Infinite Series Page 1 of 11. Chapter 2 : Infinite Series

CDS 101: Lecture 5.1 Reachability and State Space Feedback

Discrete Fourier Transform (DFT)

Session : Plasmas in Equilibrium

S- AND P-POLARIZED REFLECTIVITIES OF EXPLOSIVELY DRIVEN STRONGLY NON-IDEAL XENON PLASMA

Hydrogen Atom and One Electron Ions

NET/JRF, GATE, IIT JAM, JEST, TIFR

Ordinary Differential Equations

1973 AP Calculus AB: Section I

They must have different numbers of electrons orbiting their nuclei. They must have the same number of neutrons in their nuclei.

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

Digital Signal Processing, Fall 2006

A Prey-Predator Model with an Alternative Food for the Predator, Harvesting of Both the Species and with A Gestation Period for Interaction

Folding of Hyperbolic Manifolds

Hadamard Exponential Hankel Matrix, Its Eigenvalues and Some Norms

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Ideal crystal : Regulary ordered point masses connected via harmonic springs

CHAPTER CHAPTER. Discrete Dynamical Systems. 9.1 Iterative Equations. First-Order Iterative Equations. (b)

07 - SEQUENCES AND SERIES Page 1 ( Answers at he end of all questions ) b, z = n

Chapter 3 Fourier Series Representation of Periodic Signals

z 1+ 3 z = Π n =1 z f() z = n e - z = ( 1-z) e z e n z z 1- n = ( 1-z/2) 1+ 2n z e 2n e n -1 ( 1-z )/2 e 2n-1 1-2n -1 1 () z

Available online at Energy Procedia 4 (2011) Energy Procedia 00 (2010) GHGT-10

Introduction to Quantum Information Processing. Overview. A classical randomised algorithm. q 3,3 00 0,0. p 0,0. Lecture 10.

Blackbody Radiation. All bodies at a temperature T emit and absorb thermal electromagnetic radiation. How is blackbody radiation absorbed and emitted?

ASSERTION AND REASON

Chapter At each point (x, y) on the curve, y satisfies the condition

Chapter 6 Folding. Folding

15/03/1439. Lectures on Signals & systems Engineering

t i Extreme value statistics Problems of extrapolating to values we have no data about unusually large or small ~100 years (data) ~500 years (design)

ELECTRONIC APPENDIX TO: ELASTIC-PLASTIC CONTACT OF A ROUGH SURFACE WITH WEIERSTRASS PROFILE. Yan-Fei Gao and A. F. Bower

y = 2xe x + x 2 e x at (0, 3). solution: Since y is implicitly related to x we have to use implicit differentiation: 3 6y = 0 y = 1 2 x ln(b) ln(b)

Normal Form for Systems with Linear Part N 3(n)

Mock Exam 2 Section A

Review Exercises. 1. Evaluate using the definition of the definite integral as a Riemann Sum. Does the answer represent an area? 2

DFT: Discrete Fourier Transform

Physics 2D Lecture Slides Lecture 14: Feb 3 rd 2004

Electromagnetics Research Group A THEORETICAL MODEL OF A LOSSY DIELECTRIC SLAB FOR THE CHARACTERIZATION OF RADAR SYSTEM PERFORMANCE SPECIFICATIONS

Self-interaction mass formula that relates all leptons and quarks to the electron

Global Chaos Synchronization of the Hyperchaotic Qi Systems by Sliding Mode Control

1973 AP Calculus BC: Section I

Chapter 1. Chapter 10. Chapter 2. Chapter 11. Chapter 3. Chapter 12. Chapter 4. Chapter 13. Chapter 5. Chapter 14. Chapter 6. Chapter 7.

ONLINE SUPPLEMENT Optimal Markdown Pricing and Inventory Allocation for Retail Chains with Inventory Dependent Demand

1997 AP Calculus AB: Section I, Part A

Periodic Structures. Filter Design by the Image Parameter Method

Chapter 4 - The Fourier Series

to the SCHRODINGER EQUATION The case of an electron propagating in a crystal lattice

ELG3150 Assignment 3

Construction of Mimetic Numerical Methods

Lectures 9 IIR Systems: First Order System

Discrete Fourier Transform. Nuno Vasconcelos UCSD

Calculus II (MAC )

MILLIKAN OIL DROP EXPERIMENT

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 12

Elements of Statistical Thermodynamics

AS 5850 Finite Element Analysis

Land concolidation planning

Australian Journal of Basic and Applied Sciences, 4(9): , 2010 ISSN

The Interplay between l-max, l-min, p-max and p-min Stable Distributions

Solution evolutionary method of compressible boundary layers stability problems

On the approximation of the constant of Napier

NEW VERSION OF SZEGED INDEX AND ITS COMPUTATION FOR SOME NANOTUBES

Thomas Whitham Sixth Form

Journal of Modern Applied Statistical Methods

WBJEEM MATHEMATICS. Q.No. μ β γ δ 56 C A C B

2011 HSC Mathematics Extension 1 Solutions

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

Triple Play: From De Morgan to Stirling To Euler to Maclaurin to Stirling

Empirical Study in Finite Correlation Coefficient in Two Phase Estimation

Chapter 11.00C Physical Problem for Fast Fourier Transform Civil Engineering

INC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor

Lecture contents. Bloch theorem k-vector Brillouin zone Almost free-electron model Bands Effective mass Holes. NNSE 508 EM Lecture #9

SOLUTIONS TO CHAPTER 2 PROBLEMS

( ) ( ) ( ) 2011 HSC Mathematics Solutions ( 6) ( ) ( ) ( ) π π. αβ = = 2. α β αβ. Question 1. (iii) 1 1 β + (a) (4 sig. fig.

امتحانات الشهادة الثانوية العامة فرع: العلوم العامة

Solid State Device Fundamentals

DTFT Properties. Example - Determine the DTFT Y ( e ) of n. Let. We can therefore write. From Table 3.1, the DTFT of x[n] is given by 1

For more important questions visit :

Numbering Systems Basic Building Blocks Scaling and Round-off Noise. Number Representation. Floating vs. Fixed point. DSP Design.

Transcription:

Mixd Mod Oscillatios as a Mchaism for Psudo-Platau Burstig Richard Brtram Dpartmt of Mathmatics Florida Stat Uivrsity Tallahass, FL

Collaborators ad Support Thodor Vo Marti Wchslbrgr Joël Tabak Uivrsity of Sydy Florida Stat Uivrsity Full articl: Vo t al., J. Comput. Nurosci., 2010 Supportd by NIH grat DA 43200 ad NSF grat DMS 0917664

Ctral Aim Udrstad th dyamics of a ovl form of burstig i a modl for xcitabl docri clls

Psudo-Platau Burstig Oft Occurs i Edocri Clls Prforatd patch rcordig from a GH4 pituitary cll li Simulatio from a mathmatical modl (Toporikova t al., 2008)

Thr ar Svral Modls That Produc Psudo-Platau Burstig Lactotrophs: Tabak t al., J. Comput. Nurosci., 2007 Toporikova t al., Nural Computatio, 2008 Somatotrophs: Tsava-Ataasova t al., J. Nurophysiol., 2007 Corticotrophs: LBau t al., J. Thortical Biology, 1998 Short t al., J. Thortical Biology, 2000 Sigl β-clls: Zhag t al., Biophysical J., 2003

Typical Psudo-Platau Fast- Subsystm Bifurcatio Structur V (mv) [Ca] (μm) Ulik squar-wav burstig (platau burstig), th top brach of th z-curv is stabilizd; thr is o stabl spikig brach.

O Form of Psudo-Platau Burstig has a Diffrt Structur V (mv) [Ca] (μm) What is th mchaism of this burstig oscillatio?

Th Mathmatical Modl Calcium coctratio has b rmovd; ot cssary for th burstig. C dv dt = g Ca m ( V )( V V ) + g ( V V ) Ca K K + g A a ( V ) ( V K V ) + g L ( V K V ) τ d dt = ( V ) d τ = dt ( V ) Th ad variabls chag o a slowr tim scal tha V. Thr ar 2 slow variabls ad 1 fast variabl.

Burstig Occurs Ovr a Rag of g K Valus burstig spikig V (mv) 10 20 50 80 PD 1 HP 1 PD 2 SN 1 HM SN 2 0 5 10 15 20 g K (S) Black = statioary Rd = priodic (spikig)

Burstig Occurs i a Rgio of th g K - g A Paramtr Spac g A (S) 25 20 15 10 5 0 dp/hyp dp SN 1 hyp HP 1 /PD 1 PD 2 burstig spikig 3 4 5 6 7 g k (S) Paramtrs ar th maximum coductacs corrspodig to th two slow variabls.

Mchaism of Burstig: Go to th Sigular Limit ( C 0) Paramtrs st to produc spikig. Th spikig solutio bcoms a rlaxatio oscillatio o th critical maifold; this is th quasi-quilibrium surfac for th V variabl.

Th Rducd ad Dsigularizd Systms ) ( ),, ( L A K Ca I I I I V f + + + RHS of V-ODE: { } 0 ),, ( : ),, ( 3 = R V f V S Critical maifold: 0 ),, ( dt d V f dt d = Dyamics o S: Dsigularizd systm: ),, ( V F f f d dv + = τ τ τ V f d d = τ τ t V f 1 τ with Rducd systm: f f dt dv V f + = τ τ = V dt d τ ) ( ), ( V = with =0 o folds Goal: Driv quatios for th flow o th critical maifold.

Foldd Sigularity A foldd sigularity of th rducd systm is a quilibrium of th dsigularizd systm that occurs o a fold curv, ad satisfis f ( V,, ) = 0 o th critical maifold F( V,, ) = 0 V tim drivativ is 0 i dsigularizd systm f V = 0 o a fold curv of S

Foldd Nod Sigularity A foldd od sigularity (FN) is a foldd sigularity with gativ ral igvalus. For small valus of C (larg, but ot ifiit, tim scal sparatio) th slow maifold is twistd i th ighborhood of th FN. From Dsrochs t al., Chaos, v. 18, 2008

Th Sigular Ful Th sigular ful of a foldd od is dlimitd by th fold curv (L + ) ad th strog sigular caard (SC). This is th trajctory that is tagt to th igdirctio of th strog igvalu of th FN. SC FN ful L + I th sigular limit, trajctoris trig th ful mov through th FN i fiit tim ad follow th middl sht of S for som tim. For small, but o-zro C, corrspodig trajctoris xhibit small oscillatios du to th twistd slow maifold.

Rlaxatio Oscillatios (Spikig) do ot Etr th Sigular Ful g A =0.2 S ful

Mixd-Mod Oscillatios (Burstig) Etr th Sigular Ful ful g A =4 S For small C, th trajctory oscillats oc it jumps up. Th small oscillatios i combiatio with th larg jumps is a mixd-mod oscillatio. Th small oscillatios ar th spiks of th psudo-platau burst.

μ = 0 Caard Thory Provids th Burstig Bordrs μ = 0 FN bcoms a foldd saddl, oly o caard For C>0 this is a Hopf bifurcatio μ is igvalu ratio λ μ = 1 λ 2 ( 0,1) =1 μ Dgrat FN. To th right it bcoms a foldd focus, with o caards. δ is th distac of th sigular orbit from th SC, withi th sigular ful. δ>0 i th MMO rgio.

Spiks Emrg as C is Icrasd Du to jump up Du to th FN

Maximum Numbr of Small Oscillatios For C sufficitly small, is Giv by a Formula s max = μ + 1 2μ whr [ ] is th gratst itgr fuctio. For drivatio s Wchslbrgr, SIAM J. Dy. Syst., 2005. This is issitiv to chags i g A. Blu: s max =1 Gr: s max =2 Yllow: s max =5 Rd: s max =12

Actual Numbr of Small Oscillatios Dtrmid by Whr th Sigular Trajctory Etrs th Ful δ For C small ad S max =3: SC S=1 S=2 scodary caards FN L + S=3 Numbr of small oscillatios icrass with δ, th distac from th SC.

Numrical Simulatio Agrs with th Caard Thory Numbr of spiks pr burst, with C=2 pf. s max δ Small, dark blu = spikig Small, light blu = burst with 2 spiks Larg, rd = burst with 53 spiks

Thak You! (Postdoc positio availabl)