12.07.2015 Views

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

Dynamical Systems in Neuroscience:

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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Simple Models 29920 mV100 ms49 pArecovery variablespik<strong>in</strong>g limitcycleBABA-60 mVmembrane potentialFigure 8.16: Stutter<strong>in</strong>g behavior of an RS neuron (data were provided by Dr. Klaus M.Stiefel, P28-36 adult mouse, coronal slices, 300µm, layer II/III pyramid, visual cortex).attraction doma<strong>in</strong> of the focus (shaded region <strong>in</strong> the figure) is bounded by the stablemanifold of the saddle. As I <strong>in</strong>creases, the stable manifold makes a loop and becomesa big homocl<strong>in</strong>ic orbit giv<strong>in</strong>g birth to a spik<strong>in</strong>g limit cycle attractor. When I = 125,stable rest<strong>in</strong>g and spik<strong>in</strong>g states co-exist, which plays an important role <strong>in</strong> expla<strong>in</strong><strong>in</strong>gthe paradoxical stutter<strong>in</strong>g behavior of some neocortical neurons discussed later. AsI <strong>in</strong>creases, the saddle quantity, i.e., the sum of eigenvalues of the saddle, becomespositive. When the stable manifold makes another, smaller loop, it gives birth toan unstable limit cycle, which then shr<strong>in</strong>ks to the rest<strong>in</strong>g equilibrium and results <strong>in</strong>subcritical Andnronov-Hopf bifurcation.What is the excitability class of the RS model neuron <strong>in</strong> Fig. 8.15? If a slow ramp ofcurrent is <strong>in</strong>jected, the rest<strong>in</strong>g state of the neuron becomes a stable focus and then losesstability via subcritical Andronov-Hopf bifurcation. Hence the neuron is a resonatorexhibit<strong>in</strong>g Class 2 excitability. Now suppose steps of dc-current of amplitude I = 125pA or less are <strong>in</strong>jected. The trajectory starts at the <strong>in</strong>itial po<strong>in</strong>t (v, u) = (−60, 0),which is the rest<strong>in</strong>g state when I = 0, and then approaches the spik<strong>in</strong>g limit cycle.Because the limit cycle was born via a homocl<strong>in</strong>ic bifurcation to the saddle, it has alarge period and hence the neuron is Class 1 excitable. Thus, depend<strong>in</strong>g on the natureof stimulation, i.e., ramps vs. pulses, we can observe small or large spik<strong>in</strong>g frequencies,at least <strong>in</strong> pr<strong>in</strong>ciple.In practice, it is quite difficult to catch homocl<strong>in</strong>ic orbits to saddles because they aresensitive to noise. Injection of a constant current just below the neuron’s rheobase <strong>in</strong>Fig. 8.15 would result <strong>in</strong> random transitions between the rest<strong>in</strong>g state and a periodicspik<strong>in</strong>g state. Indeed, the two attractors co-exist and are near each other, so weakmembrane noise can push the trajectory <strong>in</strong> and out of the shaded region result<strong>in</strong>g <strong>in</strong>a stutter<strong>in</strong>g spik<strong>in</strong>g, illustrated <strong>in</strong> Fig. 8.16, m<strong>in</strong>gled with subthreshold oscillations.Such a behavior is also exhibited by FS <strong>in</strong>terneurons studied later <strong>in</strong> this section.

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

Saved successfully!

Ooh no, something went wrong!