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.

252 Excitability(a) <strong>in</strong>tegrator(b) resonator?new rest<strong>in</strong>gstatethresholdspikethresholdsetnew rest<strong>in</strong>gstatePSPspikeABold rest<strong>in</strong>gstate, I=0old rest<strong>in</strong>gstate, I=0Figure 7.28: Integrators have well-def<strong>in</strong>ed rheobase current while resonators may not.spond<strong>in</strong>g to partial-amplitude action potentials, such as those <strong>in</strong> Fig. 7.26. No s<strong>in</strong>glecurve separates small potentials from action potentials, so there is no well-def<strong>in</strong>edthreshold manifold.FitzHugh (1955) noticed that the threshold set can be quite th<strong>in</strong> <strong>in</strong> some models,<strong>in</strong>clud<strong>in</strong>g the Hodgk<strong>in</strong>-Huxley model. In particular, the difference between the trajectoriescorrespond<strong>in</strong>g to small potentials and action potentials can be as small as 0.0001mV, which is smaller than the noisy fluctuations of the membrane potential. Thus,to observe an <strong>in</strong>termediate-amplitude spike <strong>in</strong> such models, one needs to simulate themodels with accuracy beyond the limits of uncerta<strong>in</strong>ty which appear when the physical<strong>in</strong>terpretation of the model is considered. As a result, for any practical purposesuch models exhibit all-or-none behavior with the threshold set look<strong>in</strong>g like a thresholdmanifold. FitzHugh referred to this as be<strong>in</strong>g a quasithreshold phenomenon.Quasi-thresholds are related to the special canard trajectory depicted <strong>in</strong> Fig. 7.27d.The trajectory follows the unstable branch of the cubic nullcl<strong>in</strong>e all the way to the rightknee po<strong>in</strong>t P. The flow near the trajectory is highly unstable; any small perturbationpushes the state of the system to the left or to the right, result<strong>in</strong>g <strong>in</strong> a “subthreshold”or “superthreshold” response. The solutions depicted <strong>in</strong> Fig. 7.26, right, try to followsuch a trajectory. An easy way to compute the trajectory <strong>in</strong> two-dimensional relaxationoscillators is to start with the po<strong>in</strong>t P and <strong>in</strong>tegrate the system backwards (t → −∞).We discuss canard (French duck) solutions <strong>in</strong> detail <strong>in</strong> Sect. 6.3.4.7.2.6 RheobaseNeuronal rheobase, i.e., the m<strong>in</strong>imal amplitude of a current of <strong>in</strong>f<strong>in</strong>ite duration thatmakes the neuron fire, measures the “current threshold” of the neuron. Integrators havea well-def<strong>in</strong>ed rheobase while resonators may not. To see this, consider an <strong>in</strong>tegratorneuron <strong>in</strong> Fig. 7.28a receiv<strong>in</strong>g a current step that changes <strong>in</strong>stantly its phase portrait.In particular, the current moves the equilibrium from the old location correspond<strong>in</strong>g

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

Saved successfully!

Ooh no, something went wrong!