Ermentrout-Kopell规范模型

The Ermentrout-Kopell canonical model is better known as the "theta model" and is a simple one-dimensional model for the spiking of a neuron. It is closely related to the quadratic integrate and fire neuron. The model takes the following form:

$\displaystyle{ \label{theta} \frac{d\theta}{dt} = 1-\cos\theta + (1+\cos\theta) I(t) (1) }$

where $\displaystyle{ I(t) }$ are the inputs to the model. The variable $\displaystyle{ \theta }$ lies on the unit circle and ranges between 0 and $\displaystyle{ 2\pi\ . }$ When $\displaystyle{ \theta=\pi }$ the neuron "spikes", that is, it produces an action potential.

派生Derivation

The theta model is the normal form for the saddle-node on a limit cycle bifurcation (SNIC). (Caution! Do not confuse this with a saddle-node of limit cycles in which a pair of limit cycles collide and annihilate.) Figure 1 shows a schematic of the bifurcation as a parameter varies through the critical value of $\displaystyle{ I=0. }$ When $\displaystyle{ I\lt 0 }$ there is a pair of equilibria. One of the equilibria is a saddle point with a one-dimensional unstable manifold. The two branches of the unstable manifold form an invariant circle with the stable equilibrium point. (See saddle-node bifurcation where this is called a saddle-node homoclinic bifurcation.) In neurophysiological terms, the stable manifold of the saddle point forms a true threshold for the neuron. In Figure 1 , the stable manifold is shown in green. Any initial conditions to the left of the manifold will be attracted to the stable equilibrium (in blue), while initial data to the right of the manifold will make a large excursion around the circle before returning to the rest state.

θ模型是极限环分岔(SNIC)鞍节点的标准形式。(注意!不要将此与极限环的鞍节点混淆，在鞍节点中，一对极限环会碰撞和湮灭。)图1显示一分岔示意图，参数随临界值$\displaystyle{ I=0. }$变化，当$\displaystyle{ I\lt 0 }$时，存在一对平衡。其中一个平衡是具有一维不稳定流形unstable manifold的鞍点。不稳定流形的两个分支构成一个具有稳定平衡点的不变圆。(参见鞍节分岔saddle-node bifurcation，这被称为鞍节同宿分岔。)在神经生理学术语中，鞍点的稳定流形stable manifold为神经元形成了一个真正的阈值threshold。在图1中，稳定歧管以绿色显示。流形左边的任何初始条件都将被吸引到稳定平衡(蓝色部分)，而流形右边的初始数据将在返回到静止状态之前绕圆进行较大的偏移。

Near the transition, the local dynamics is like a saddle-node bifurcation and has the form:

$\displaystyle{ \label{sn} \frac{dx}{dt} = x^2 + I. }$

For $\displaystyle{ I\lt 0 }$ (resp $\displaystyle{ I\gt 0 }$) there are two (resp no) equilibria. In the case where $\displaystyle{ I\gt 0 }$ solutions to this differential equation "blow up" in finite time

$\displaystyle{ I\lt 0 }$(resp $\displaystyle{ I\gt 0 }$) ，存在两个(resp no)均衡。当$\displaystyle{ I\gt 0 }$的情况下，这个微分方程的解在有限时间内“爆炸”

$\displaystyle{ T_{blow} =1/2\, \left( -2\,\arctan \left( {\frac { {\it x(0)}}{\sqrt {I}} } \right) +\pi \right) {\frac {1}{\sqrt {I}}}. }$

Here $\displaystyle{ x(0) }$ is the initial condition. In particular, suppose we reset $\displaystyle{ x(t) }$ to $\displaystyle{ -\infty }$ when it blows up to $\displaystyle{ +\infty }$ Then the total transit time is

$\displaystyle{ T_{per} = \frac{\pi}{\sqrt{I}}. }$

Thus the frequency (the reciprocal of the period) goes to zero as the parameter approaches criticality from the right. This observation led Rinzel & Ermentrout to remark that this bifurcation corresponded to Hodgkin's Class I excitable membranes while the more familiar Andronov-Hopf bifurcation corresponded to Class II excitability. The latter is best exemplified by the classical Hodgkin-Huxley model for the squid axon. Neural models undergoing a SNIC bifurcation include the Connor-Stevens model for crab leg axons, the Wang-Buzsaki model for inhibitory interneurons, the Hindmarsh-Rose model, and the Morris-Lecar model under some circumstances.

The quadratic integrate and fire model is essentially equation \eqref{sn} with a finite value for the blow up and a finite reset. It is closely related to the Izhikevich neuron, which has an additional linear variable modeling the dynamics of a recovery variable.

To derive the theta model (1)from the saddle-node(2), we make a simple change of variables, $\displaystyle{ x=\tan(\theta/2) }$ from which it is simple calculus to obtain the theta model. We note that as $\displaystyle{ \theta }$ approaches $\displaystyle{ \pi }$ from the left, $\displaystyle{ x }$ goes to $\displaystyle{ +\infty\ . }$ The theta model collapses the real line to the circle. The SNIC is a global bifurcation, so that to rigorously prove the equivalence of the SNIC and the theta model requires quite a bit more work than is shown in this formal derivation. The interested reader should consult the paper by Ermentrout and Kopell (1986).

The advantage of the theta model over the quadratic integrate and fire model is that there is no reset to deal with and the resulting dynamics are smooth and stay bounded. However, as the Izhikevich neuron demonstrates, it is sometimes useful to have the freedom to reset the dynamics anywhere.

有噪声的θ模型Noisy theta models

$\displaystyle{ dx = (x^2+I(t))dt + \sigma dW }$

and make the change of variables, $\displaystyle{ x=\tan(\theta/2) \ , }$ where we are careful to account for the fact that we must use Ito Calculus. The resulting noisy theta model takes the form:

$\displaystyle{ d\theta = (1-\cos\theta + [1+\cos\theta](I(t)-\frac{\sigma^2}{2}\sin\theta))dt + \sigma(1+\cos\theta)dW. }$

Note that the sine term in the equation comes from the Ito change of variables. For small noise, that is, $\displaystyle{ \sigma \ll 1 \ , }$ this term can be neglected and one gets the equation analyzed in Gutkin and Ermentrout.

θ模型的相位重置曲线The phase resetting curve for the theta model

In the oscillatory regime, the phase resetting curve (PRC) can be computed. Izhikevich computed the PRC for finite size stimuli by adding an instantaneous pulse of size $\displaystyle{ a }$ to the quadratic version of the model. From this, he obtained a map from the old phase to the new phase:

$\displaystyle{ \theta \to 2 \arctan (\tan \frac{\theta}{2} + a) }$

Note that the PRC is set in phase coordinates rather than in time coordinates.

The adjoint or infinitesimal PRC is very easy to compute using the quadratic version of the model (Ermentrout, 1996) . For any scalar oscillator model, $\displaystyle{ du/dt=f(u) \ , }$ the adjoint is $\displaystyle{ u_a(t)=1/du/dt \ . }$ Since the "periodic" solution to the quadratic model is

$\displaystyle{ u(t) = -\sqrt{I}\cot(\sqrt{I}t) }$

the PRC is

PRC是

$\displaystyle{ PRC(t) =\frac{1}{du/dt} = \frac{1}{2\sqrt{I}}(1-\cos(2\sqrt{I}t)). }$

This is non-negative and has been suggested as the signature of neurons undergoing a SNIC bifurcation.

与其他模型的联系Relation to Other Models

The canonical model described here is closely related to other phase models arising in applications. For example, the classical description of forced oscillations in a damped pendulum is given by the differential equation

$\displaystyle{ \label{pendulum} \mu\ddot\theta+f \dot\theta+\ddot H(t)\cos\theta+\ddot V(t)\sin\theta=\omega }$

where $\displaystyle{ \theta }$ is the angle between the down direction and the radius through the center of mass, $\displaystyle{ \mu }$ is the mass, $\displaystyle{ f }$ is the coefficient of friction (damping), $\displaystyle{ \ddot H }$ and $\displaystyle{ \ddot V }$ are the horizontal and vertical accelerations of the support point, and $\displaystyle{ \omega(t) }$ is the torque applied to the support point (see Chester (1975)). This model has been applied to describe mechanical systems (eg., pendulums), micro-electromechanical systems (Hoppensteadt-Izhikevich (2001)), rotating electrical machinery (Stoker (1951)), power systems (Salam (1984)), electronic circuits, such as phase-locked loops (Viterbi (1966)) and parametric amplifiers (Horowitz-Hill (1980)), quantum mechanical devices (Feynman (1963)), and neurons (see VCON).

The model（1） is equivalent to（3） when $\displaystyle{ \mu\to 0 }$ and $\displaystyle{ H(t)\ , }$ $\displaystyle{ V(t)\ , }$ and $\displaystyle{ \omega }$ are chosen appropriately.

$\displaystyle{ \mu\to 0 }$$\displaystyle{ H(t)\ , }$$\displaystyle{ V(t)\ , }$$\displaystyle{ \omega }$选择适当时，模型(1)与(3)等价。

Although \eqref{theta} may have been referred to as being the theta-equation, this causes confusion when working with theta rhythms in the brain, and so is not preferred. Hoppensteadt and Izhikevich (1997) suggested to call it the Ermentrout-Kopell canonical model.

参考文献References

• W. Chester, The forced oscillations of a simple pendulum, J. Inst. Maths. Appl. 15 (1975) 298-306.
• G.B. Ermentrout and N. Kopell, Parabolic bursting in an excitable system coupled with a slow oscillation, SIAM-J.-Appl.-Math, 46 (1986), 233-253.
• B. Ermentrout, Type I membranes, phase resetting curves, and synchrony. Neural-Comput. 8, (1996) 979-1001
• F.C. Hoppensteadt, E.M. Izhikevich, Synchronization of MEMS resonators and mechanical neurocomputing, IEEE Trans. Circs. Systems, 48 (2001) 133-138.
• F. C. Hoppendsteadt and E.M. Izhikevich (1997) Weakly Connected Neural Networks. Springer-Verlag, NY.
• J. Stoker, Nonlinear Vibrations, Interscience, 1951.
• F.M. Salam, J.E. Marsden, P.P. Varaiya, Arnold diffusion in the swing equations of a power system, IEEE Trans. Circuits and Systems, 31 (1984) 673-688.
• A.J. Viterbi, Principles of Coherent Communication, McGraw-Hill, New York, 1966.
• P. Horowitz, W. Hill, The Art of Electronics, Camb. U. Press, 1989.
• R.P. Feynman, R. B. Leighton, M. Sands, The Feynman Lectures on Physics, Addison-Wesley, Menlo Park, 1965.

Internal references

• John Guckenheimer (2007) Bifurcation. Scholarpedia, 2(6):1517.
• Valentino Braitenberg (2007) Brain. Scholarpedia, 2(11):2918.
• Eugene M. Izhikevich (2007) Equilibrium. Scholarpedia, 2(10):2014.
• James Murdock (2006) Normal forms. Scholarpedia, 1(10):1902.
• Jeff Moehlis, Kresimir Josic, Eric T. Shea-Brown (2006) Periodic orbit. Scholarpedia, 1(7):1358.
• Philip Holmes and Eric T. Shea-Brown (2006) Stability. Scholarpedia, 1(10):1838.