Dynamical Systems in Neuroscience:
Dynamical Systems in Neuroscience: Dynamical Systems in Neuroscience:
424 Solutions to Exercises, Chap. 62. (a) The systemẋ = bx 2 , b ≠ 0cannot exhibit saddle-node bifurcation: It has one equilibrium for any non-zero b, or aninfinite number of equilibria when b = 0. The equilibrium x = 0 is non-hyperbolic and thenon-degeneracy condition is satisfied (a = b ≠ 0). However, the transversality conditionis not satisfied at the equilibrium x = 0. Another example is ẋ = b 2 + x 2 .(b) The systemẋ = b − x 3has a single stable equilibrium for any b. However, the point x = 0 is non-hyperbolic whenb = 0, and the transversality condition is also satisfied. The non-degeneracy condition isviolated, though.3. It is easy to check (by differentiating) that√c(b − bsn )V (t) = √ tan( √ ac(b − b sn )t)ais a solution to the system. Since tan(−π/2) = −∞ and tan(+π/2) = +∞, it takesfor the solution to go from −∞ to +∞.T =π√ac(b − bsn )4. The first system can be transformed into the second one if we use complex coordinates z = u+iv.To obtain the third system, we use polar coordinatesso thatre iϕ = z = u + iv ∈ C ,ż(c(b)+iω(b))z+(a+id)z|z| 2{ }} { { }} {ṙe iϕ + re iϕ i ˙ϕ = (c(b) + iω(b))re iϕ + (a + id)r 3 e iϕ .Next, we divide both sides of this equation by e iϕ and separate the real and imaginary partsto obtain{ṙ − c(b)r − ar 3 } + ir{ ˙ϕ − ω(b) − dr 2 } = 0 ,which we can write in the polar-coordinates form.5. (a) The equilibrium r = 0 of the systemṙ = br 3 ,˙ϕ = 1 ,has a pair of complex-conjugate eigenvalues ±i for any b, and the non-degeneracy conditionis satisfied for any b ≠ 0. However, the transversality condition is violated, andthe system does not exhibit Andronov-Hopf bifurcation (no limit cycle exists near theequilibrium).(b) The equilibrium r = 0 for b = 0ṙ = br ,˙ϕ = 1 ,has a pair of complex-conjugate eigenvalues ±i and the transversality condition is satisfied.However, the bifurcation is not of the Andronov-Hopf type because no limit cycleexists near the equilibrium for any b.
Solutions to Exercises, Chap. 6 4256. The Jacobian matrix at the equilibrium (u, v) = (0, 0) has the form( )b −1L =.1 bIt has eigenvalues b±i. Therefore, the loss of stability occurs at b = 0, and the non-hyperbolicityand transversality conditions are satisfied. Since the model can be reduced to the polarcoordinatesystem (see Ex. 4), and a ≠ 0, the non-degeneracy condition is also satisfied, andthe system undergoes an Andronov-Hopf bifurcation.7. Since(cr + ar 3 ) ′ r = c + 3ar 2 = c + 3a|c/a| =the limit cycle is stable when a < 0{c + 3|c| when a > 0,c − 3|c| when a < 0,8. The sequence of bifurcations is similar to that of the RS neuron in Fig. 8.15. The restingstate is a globally asymptotically stable equilibrium for I < 5.64. At this value a stable(spiking) limit cycle appears via a big saddle homoclinic orbit bifurcation. At I = 5.8 asmall-amplitude unstable limit cycle is born via another saddle homoclinic orbit bifurcation.This cycle shrinks to the equilibrium and makes it lose stability via subcritical Andronov-Hopfbifurcation at I = 6.5. This unstable focus becomes an unstable node when I increases, andthen it coalesces with the saddle (at I = 7.3) and disappears. Notice that there is a saddle-nodebifurcation according to the I-V relation, but it corresponds to the disappearance of an unstableequilibrium.9. The Jacobian matrix of partial derivatives has the form( )−I′L = V (V, x) −I x(V, ′ x)x ′ ,∞(V )/τ(V ) −1/τ(V )so thatandThe characteristic equationhas two solutionstr L = −{I ′ V (V, x) + 1/τ(V )}det L = {I ′ V (V, x) + I ′ x(V, x)x ′ ∞(V )}/τ(V ) = I ′ ∞(V )/τ(V ) .which might be complex-conjugate.10. Let z = re ϕi , thenλ 2 − λ tr L + det L = 0c{ }} { { √ }} {(tr L)/2 ± {(tr L)/2}2 − det Lωr ′ = ar + r 3 − r 5 ,ϕ ′ = ω .Any limit cycle is an equilibrium of the amplitude equation, i.e.,a + r 2 − r 4 = 0 .The system undergoes fold limit cycle bifurcation when the amplitude equation undergoes asaddle-node bifurcation, i.e., whena + 3r 2 − 5r 4 = 0(check the non-degeneracy and transversality conditions). The two equations have the nontrivialsolution (a, r) = (−1/4, 1/ √ 2).
- Page 384 and 385: 374 Burstingdepending on the type o
- Page 386 and 387: 376 BurstingsubcriticalAndronov-Hop
- Page 388 and 389: 378 Bursting(a)membrane potential,
- Page 390 and 391: spiking380 Burstingfoldbifurcations
- Page 392 and 393: 382 BurstingspikingsupercriticalAnd
- Page 394 and 395: 384 Burstingaction potentials cut2
- Page 396 and 397: 386 Bursting9.4.3 BistabilitySuppos
- Page 398 and 399: 388 BurstingFigure 9.49: The instan
- Page 400 and 401: esting390 Burstingspikesynchronizat
- Page 402 and 403: 392 BurstingReview of Important Con
- Page 404 and 405: 394 BurstingspikingrestingFigure 9.
- Page 406 and 407: 396 Bursting0-10membrane potential,
- Page 408 and 409: 398 Bursting0membrane potential, V
- Page 410 and 411: 400 BurstingFigure 9.61: A cycle-cy
- Page 412 and 413: 402 Bursting28. [Ph.D.] Develop an
- Page 414 and 415: 404 Synchronization (see www.izhike
- Page 416 and 417: 406 Synchronization (see www.izhike
- Page 418 and 419: 408 Solutions to Exercises, Chap. 3
- Page 420 and 421: 410 Solutions to Exercises, Chap. 3
- Page 422 and 423: 412 Solutions to Exercises, Chap. 3
- Page 424 and 425: 414 Solutions to Exercises, Chap. 4
- Page 426 and 427: 416 Solutions to Exercises, Chap. 4
- Page 428 and 429: 418 Solutions to Exercises, Chap. 4
- Page 430 and 431: 420 Solutions to Exercises, Chap. 4
- Page 432 and 433: 422 Solutions to Exercises, Chap. 5
- Page 436 and 437: 426 Solutions to Exercises, Chap. 6
- Page 438 and 439: 428 Solutions to Exercises, Chap. 8
- Page 440 and 441: 430 Solutions to Exercises, Chap. 9
- Page 442 and 443: 432 Solutions to Exercises, Chap. 9
- Page 444 and 445: 434 Solutions to Exercises, Chap. 9
- Page 446 and 447: 436 Solutions to Exercises, Chap. 9
- Page 448 and 449: 438 Solutions to Exercises, Chap. 9
- Page 450 and 451: 440 Solutions to Exercises, Chap. 9
- Page 452 and 453: 442 Referencesterneurons mediated b
- Page 454 and 455: 444 ReferencesDickson C.T., Magistr
- Page 456 and 457: 446 ReferencesGuckenheimer J. (1975
- Page 458 and 459: 448 Referencestional Journal of Bif
- Page 460 and 461: 450 ReferencesMarkram H, Toledo-Rod
- Page 462 and 463: 452 ReferencesRosenblum M.G. and Pi
- Page 464 and 465: 454 ReferencesTuckwell H.C. (1988)
- Page 466 and 467: 456 References9
- Page 468 and 469: 458 Synchronization (see www.izhike
- Page 470 and 471: 460 Synchronization (see www.izhike
- Page 472 and 473: 462 Synchronization (see www.izhike
- Page 474 and 475: 464 Synchronization (see www.izhike
- Page 476 and 477: 466 Synchronization (see www.izhike
- Page 478 and 479: 468 Synchronization (see www.izhike
- Page 480 and 481: 470 Synchronization (see www.izhike
- Page 482 and 483: 472 Synchronization (see www.izhike
Solutions to Exercises, Chap. 6 4256. The Jacobian matrix at the equilibrium (u, v) = (0, 0) has the form( )b −1L =.1 bIt has eigenvalues b±i. Therefore, the loss of stability occurs at b = 0, and the non-hyperbolicityand transversality conditions are satisfied. S<strong>in</strong>ce the model can be reduced to the polarcoord<strong>in</strong>atesystem (see Ex. 4), and a ≠ 0, the non-degeneracy condition is also satisfied, andthe system undergoes an Andronov-Hopf bifurcation.7. S<strong>in</strong>ce(cr + ar 3 ) ′ r = c + 3ar 2 = c + 3a|c/a| =the limit cycle is stable when a < 0{c + 3|c| when a > 0,c − 3|c| when a < 0,8. The sequence of bifurcations is similar to that of the RS neuron <strong>in</strong> Fig. 8.15. The rest<strong>in</strong>gstate is a globally asymptotically stable equilibrium for I < 5.64. At this value a stable(spik<strong>in</strong>g) limit cycle appears via a big saddle homocl<strong>in</strong>ic orbit bifurcation. At I = 5.8 asmall-amplitude unstable limit cycle is born via another saddle homocl<strong>in</strong>ic orbit bifurcation.This cycle shr<strong>in</strong>ks to the equilibrium and makes it lose stability via subcritical Andronov-Hopfbifurcation at I = 6.5. This unstable focus becomes an unstable node when I <strong>in</strong>creases, andthen it coalesces with the saddle (at I = 7.3) and disappears. Notice that there is a saddle-nodebifurcation accord<strong>in</strong>g to the I-V relation, but it corresponds to the disappearance of an unstableequilibrium.9. The Jacobian matrix of partial derivatives has the form( )−I′L = V (V, x) −I x(V, ′ x)x ′ ,∞(V )/τ(V ) −1/τ(V )so thatandThe characteristic equationhas two solutionstr L = −{I ′ V (V, x) + 1/τ(V )}det L = {I ′ V (V, x) + I ′ x(V, x)x ′ ∞(V )}/τ(V ) = I ′ ∞(V )/τ(V ) .which might be complex-conjugate.10. Let z = re ϕi , thenλ 2 − λ tr L + det L = 0c{ }} { { √ }} {(tr L)/2 ± {(tr L)/2}2 − det Lωr ′ = ar + r 3 − r 5 ,ϕ ′ = ω .Any limit cycle is an equilibrium of the amplitude equation, i.e.,a + r 2 − r 4 = 0 .The system undergoes fold limit cycle bifurcation when the amplitude equation undergoes asaddle-node bifurcation, i.e., whena + 3r 2 − 5r 4 = 0(check the non-degeneracy and transversality conditions). The two equations have the nontrivialsolution (a, r) = (−1/4, 1/ √ 2).