12.07.2015 Views

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

236 Excitabilitysaddle-node bifurcationsubcritical Andronov-Hopf bifurcation<strong>in</strong>hibitorypulseattractiondoma<strong>in</strong>n-nullcl<strong>in</strong>espik<strong>in</strong>g limit cycleexcitatorypulsenV-nullcl<strong>in</strong>eVV-nullcl<strong>in</strong>eexc.pulsespik<strong>in</strong>g limit cycleunstablen-nullcl<strong>in</strong>e<strong>in</strong>hibitorypulseattractiondoma<strong>in</strong>of spik<strong>in</strong>glimit cycleFigure 7.11: Co-existence of stable equilibrium and spik<strong>in</strong>g limit cycle attractor <strong>in</strong> theI Na,p +I K -model. Left: The rest state is about to disappear via saddle-node bifurcation.Right: The rest state is about to lose stability via subcritical Andronov-Hopfbifurcation. Right (left) arrows denote the location and the direction of an excitatory(<strong>in</strong>hibitory) pulse that switches spik<strong>in</strong>g behavior to rest<strong>in</strong>g.Fig. 7.10d and e, where we consider the model’s phase portraits. Notice the co-existenceof the rest<strong>in</strong>g state and a limit cycle attractor. The rest<strong>in</strong>g state loses stability via subcriticalAndronov-Hopf bifurcation at I = 5.25, so the emerg<strong>in</strong>g spik<strong>in</strong>g has non-zerofrequency at I ≈ 5.25. However, <strong>in</strong>ject<strong>in</strong>g steps of current results <strong>in</strong> transitions to thelimit cycle even before the rest<strong>in</strong>g state loses its stability. The limit cycle <strong>in</strong> the modelappears via saddle homocl<strong>in</strong>ic orbit bifurcation at I ≈ 3.8866, its period is quite largeresult<strong>in</strong>g <strong>in</strong> the Class 1 response to steps of current. The F-I curves for homocl<strong>in</strong>icbifurcations have logarithmic scal<strong>in</strong>g, so small frequency oscillations are difficult tocatch numerically let alone experimentally.The surpris<strong>in</strong>g discrepancy <strong>in</strong> Fig. 7.10a occurs because the rest<strong>in</strong>g state of theI Na,p +I K -model is near the Bogdanov-Takens bifurcation, i.e., the model is near atransition from resonator to <strong>in</strong>tegrator. Such a bifurcation was recorded, though <strong>in</strong>directly,<strong>in</strong> some neocortical pyramidal neurons, as we will show later <strong>in</strong> this chapter and<strong>in</strong> Chap. 8. Another surpris<strong>in</strong>g example of Andronov-Hopf bifurcation with Class 1excitability is presented <strong>in</strong> Ex. 6. To avoid such surprises, we adopt the ramp def<strong>in</strong>itionof excitability throughout the book.7.1.6 BistabilityWhen transition from rest<strong>in</strong>g to spik<strong>in</strong>g states occurs via saddle-node (off <strong>in</strong>variantcircle) or subcritical Andronov-Hopf bifurcation, there is a co-existence of a stableequilibrium and a stable limit cycle attractor just before the bifurcation, as we illustrate<strong>in</strong> Fig. 7.11. We refer to such systems as bistable. They have a remarkable

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!