of a lumped model of neocortex to study epileptiform activity Sid Visser Hil Meijer Stephan van Gils March 21, 2012
What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized by increased probability of recurring seizures Seizures The activity of braincells hypersynchronizes Observed as large oscillations
What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized by increased probability of recurring seizures Seizures The activity of braincells hypersynchronizes Observed as large oscillations
EEG during a seizure
Goal of this research Goals Study several models of neural activity in neocortex Identify correspondences between these models Extrapolate results from simpler models to more complex models Results so far A large, detailed model describing activity of individual neurons in both normal and epileptiform states A simple, lumped model that has shown to have some similarity [Visser et al., 2010]
Motivation of model Simulation on Youtube
Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
Simplifications Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Make system symmetric: µ 1 = µ 2 := µ F 1 = F 2 := F G 1 = G 2 := G Choose F(x) and G(x) as rescalings of: S = (tanh(x a) tanh( a)) cosh 2 ( a)
Time rescaling Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Time is rescaled to non-dimensionalize system: ẋ 1 (t) = x 1 (t) α 1 S(β 1 x 1 (t τ 1 )) + α 2 S(β 2 x 2 (t τ 2 )) ẋ 2 (t) = x 2 (t) α 1 S(β 1 x 2 (t τ 1 )) + α 2 S(β 2 x 1 (t τ 2 )) Introduce vector notation: ẋ(t) = f(x t ), with x t C([ h, 0], R 2 ) and h = max(τ 1, τ 2 ). Main parameters of interest: representing connection strengths between populations.
Time rescaling Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Time is rescaled to non-dimensionalize system: ẋ 1 (t) = x 1 (t) α 1 S(β 1 x 1 (t τ 1 )) + α 2 S(β 2 x 2 (t τ 2 )) ẋ 2 (t) = x 2 (t) α 1 S(β 1 x 2 (t τ 1 )) + α 2 S(β 2 x 1 (t τ 2 )) Introduce vector notation: ẋ(t) = f(x t ), with x t C([ h, 0], R 2 ) and h = max(τ 1, τ 2 ). Main parameters of interest: representing connection strengths between populations.
Fixed points Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Corollary The origin is always a fixed point of the DDE. Theorem The DDE only has only symmetric fixed points, i.e. f(x) = 0 = x = (x, x ).
Linearization Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient The DDE is linearized about an equilibrium: in which u 1 (t) = u 1 (t) k 1 u 1 (t τ 1 ) + k 2 u 2 (t τ 2 ), u 2 (t) = u 2 (t) k 1 u 2 (t τ 1 ) + k 2 u 1 (t τ 2 ), k 1 := α 1 β 1 S (β 1 x ) k 2 := α 2 β 2 S (β 2 x ). Note that k 1 and k 2 depend on the equilibrium at which we linearize.
Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Characteristic equation and eigenvalues Consider solutions of the form u(t) = e λt c, c R 2. Non-trivial solutions of the system correspond with (λ)c = 0, for: [ λ + 1 + k1 e (λ) = λτ 1 k 2 e λτ ] 2 k 2 e λτ 2 λ + 1 + k 1 e λτ 1 Non-trivial solutions exist when det (λ) = 0, or: (λ + 1 + k 1 e λτ 1 + k 2 e λτ 2 ) (λ + 1 + k }{{} 1 e λτ 1 k 2 e λτ 2 ) = 0. }{{} := +(λ) := (λ) This characteristic function has an infinite (but countable) number of roots.
Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Characteristic equation and eigenvalues Consider solutions of the form u(t) = e λt c, c R 2. Non-trivial solutions of the system correspond with (λ)c = 0, for: [ λ + 1 + k1 e (λ) = λτ 1 k 2 e λτ ] 2 k 2 e λτ 2 λ + 1 + k 1 e λτ 1 Non-trivial solutions exist when det (λ) = 0, or: (λ + 1 + k 1 e λτ 1 + k 2 e λτ 2 ) (λ + 1 + k }{{} 1 e λτ 1 k 2 e λτ 2 ) = 0. }{{} := +(λ) := (λ) This characteristic function has an infinite (but countable) number of roots.
Symmetries of solutions Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Roots of correspond to symmetric solutions, whereas roots of + relate to asymetric solutions. Proof. Let Z 2 act on C 2 so that 1 Z 2 acts as ξ(x, y) : (x, y) (y, x), then: { (λ)v = 0 (λ) = 0 ξv = v { (λ)v = 0 + (λ) = 0 ξv = v
Symmetries of solutions Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Roots of correspond to symmetric solutions, whereas roots of + relate to asymetric solutions. Proof. Let Z 2 act on C 2 so that 1 Z 2 acts as ξ(x, y) : (x, y) (y, x), then: { (λ)v = 0 (λ) = 0 ξv = v { (λ)v = 0 + (λ) = 0 ξv = v
Roots of (λ) Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem For λ C 0 e λh and λ > C, with C C 0, (λ) > 0. Hence, (λ) has a zero free right half-plane [Bellman and Cooke, 1963].
Minimal Stability Region Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem For k 1 + k 2 < 1 the system has no eigenvalues with positive real part.
Bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Note: As in ODEs, an equilibrium of a DDE loses its stablity when either a single real eigenvalue passes through zero (fold/transcritical), or a pair of complex eigenvalues passes through the imaginary axis (Hopf).
Fold/transcritical bifurcation Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Equilibria undergo a fold or transcritical bifurcation on the lines 1 + k 1 ± k 2 = 0 in (k 1, k 2 )-space.
Hopf bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Two curves, h + (ω) and h (ω), exist along which the linearized system has a pair of complex eigenvalues λ = ±iω. Proof. For now, we consider only + (iω) = 0; (iω) is similar. iω + 1 + k 1 e iωτ 1 + k 2 e iωτ 2 = 0 [ ] [ ] [ ] cos(ωτ1 ) cos(ωτ 2 ) k1 1 = sin(ωτ 1 ) sin(ωτ 2 ) ω k 2
Hopf bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Proof cntd. The matrix is invertible if det = sin(ω(τ 2 τ 1 )) 0, yielding: [ ] [ ] [ ] k1 1 sin(ωτ2 ) cos(ωτ = h + (ω) := 2 ) 1. sin(ω(τ 2 τ 1 )) sin(ωτ 1 ) cos(ωτ 1 ) ω k 2
Bifurcations in (k 1, k 2 )-space Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Very sensitive, complex branch structure Intersections correspond with codim-2 bifurcations. k 2 k 1 (τ 1 = 11.6, τ 2 = 20.3)
Returning branches Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem For a given ω 0 > 0, (λ) has no roots λ = ±iω, ω > ω 0 inside the square k 1 + k 2 < 1 + ω0 2. Proof. First, assume iω is a root of + (a similar argument holds for ): 0 = 1 + iω + k 1 e iωτ 1 + k 2 e iωτ 2, 1 + iω k 1 k 2. Then, this root lies outside the designated square for ω > ω 0.
Revised stability region Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Conjecture For the considered parameters τ 1 and τ 2 the full stability region is bounded, connected and contains the origin.
The first Lyapunov coefficient Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Stability region partly bounded by curve of Hopfs The first Lyapunov coefficient is determined along the boundary Normal given by [Diekmann et al, 1995]: c 1 = 1 2 qt D 3 f(0)(φ, φ, φ) + q T D 2 f(0)(e 0 (0) 1 D 2 f(0)(φ, φ), φ) + 1 2 qt D 2 f(0)(e 2iω (2iω) 1 D 2 f(0)(φ, φ), φ).
The first Lyapunov coefficient Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Stability region partly bounded by curve of Hopfs The first Lyapunov coefficient is determined along the boundary Normal given by [Diekmann et al, 1995]: c 1 = 1 2 qt D 3 f(0)(φ, φ, φ) + q T D 2 f(0)(e 0 (0) 1 D 2 f(0)(φ, φ), φ) + 1 2 qt D 2 f(0)(e 2iω (2iω) 1 D 2 f(0)(φ, φ), φ).
Codim 2 bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Generalized Hopf Lypunov coefficient passes through zero, Generalized Hopf bifurcation at boundary
Summary Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Equilibria Only symmetric equilibria Stability region identified Bifurcations of stability region Fold and transcritical bifurcations Hopf-bifurcations, both sub- and supercritical Zero-Hopf and Hopf-Hopf
Summary Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Equilibria Only symmetric equilibria Stability region identified Bifurcations of stability region Fold and transcritical bifurcations Hopf-bifurcations, both sub- and supercritical Zero-Hopf and Hopf-Hopf
Numerical bifurcation analysis Overview One parameter Two parameters One parameter Identify stable and unstable manifolds of non-trivial equilibrium and periodic solutions for varying α 2. Two parameters Continue boundaries of stable solutions in α 1 and α 2. Software DDE-BIFTOOL [Engelborghs et al., 2002] PDDE-CONT/Knut [Roose and Szalai, 2007]
Numerical bifurcation analysis Overview One parameter Two parameters One parameter Identify stable and unstable manifolds of non-trivial equilibrium and periodic solutions for varying α 2. Two parameters Continue boundaries of stable solutions in α 1 and α 2. Software DDE-BIFTOOL [Engelborghs et al., 2002] PDDE-CONT/Knut [Roose and Szalai, 2007]
Overview One parameter Two parameters Bifurcations in α 2 0.2 A 2 B 0.1 1.8 0 1.6 0.1 0 3 C 2 100 200 1.4 0 100 200 3 D 2 1 1 0 0 1 0 100 200 1 0 100 200
Overview One parameter Two parameters Bifurcations in α 2 3 Bifurcation diagram II 2.5 I Maximal value 2 III 1.5 IV 1 0.5 PD1 0 H1 B1 0.5 0.4 0.5 0.6 0.7 0.8 0.9 1 α 2 Detail I 2.482 2.6 Detail II LPC4
Overview One parameter Two parameters Bifurcations in α 2 l k h Maximal value j g c i d f e α 2 a b
Fixed points in (k 1, k 2 )-space Overview One parameter Two parameters 1.6 1.2 F1 B1 k 2 0.8 H2 H3 H4 H1 Non-trivial Trivial 0.4 0 0 0.1 0.2 0.3 0.4 k 1
Overview One parameter Two parameters Bifurcations of stable solutions in α 1 and α 2 1 0.4 0.3 GH1 0.2 0.1 0 ZH2 ZH1 HH1 0.4 0.5 0.6 0.7 0.8 2
Overview One parameter Two parameters Bifurcations of stable solutions in α 1 and α 2 1 0.4 CP1 0.3 GH1 0.2 FF1 0.1 0 ZH2 ZH1 HH1 0.4 0.5 0.6 0.7 0.8 2
Overview One parameter Two parameters Bifurcations of stable solutions in α 1 and α 2 1 0.4 CP1 0.3 GH1 0.2 FF1 CP3 0.1 R1_1 0 ZH2 ZH1 HH1 0.4 0.5 0.6 0.7 0.8 2
Complex bifurcation structure Overview One parameter Two parameters α 1 R2_1 FF2 R2_2 FF3 CP4 CP2 α 2
Comparison Introduction The bifurcation analysis is compared with a survey on the detailed model:
Introduction Comparison The bifurcation analysis is compared with a survey on the detailed model:
Summary and conclusions Summary We studied a simple, lumped model for neural activity Parameter regions of steady states and periodic solutions are identified The behavior of both models varies similarly for changes of parameters Regions of multistability are of interest for studying epilepsy
Future work Introduction Details Study the model s codim-2 bifurcations Proofs rather than conjectures Next steps Expand model (e.g. break symmetry) Parameter estimation