Post-inhibitory rebound (PIR) spike induced by the negative stimulation, which plays important roles and presents counterintuitive nonlinear phenomenon in the nervous system, is mainly related to the Hopf bifurcation and hyperpolarization-active caution ($ I_h $) current. In the present paper, the emerging condition for the PIR spike is extended to the bifurcation of the big homoclinic (BHom) orbit in a model without $ I_h $ current. The threshold curve for a spike evoked from a mono-stable or coexisting steady state surrounds the steady state from left, to below, and to right, because the BHom orbit is big enough to surround the steady state. The right part of the threshold curve coincides with the stable manifold of the saddle and acts the threshold for the spike induced by the positive stimulation, resembling that of the saddle-node bifurcation on an invariant cycle, and the left part acts the threshold for the PIR spike, resembling that of the Hopf bifurcation. The bifurcation curve and a codimension-2 bifurcation point related to the BHom orbit are acquired in the two-parameter plane. The results present a comprehensive viewpoint to the dynamics near the BHom orbit bifurcation, which presents a novel threshold curve and extends the conditions for the PIR spike.
Citation: Xianjun Wang, Huaguang Gu, Bo Lu. Big homoclinic orbit bifurcation underlying post-inhibitory rebound spike and a novel threshold curve of a neuron[J]. Electronic Research Archive, 2021, 29(5): 2987-3015. doi: 10.3934/era.2021023
Post-inhibitory rebound (PIR) spike induced by the negative stimulation, which plays important roles and presents counterintuitive nonlinear phenomenon in the nervous system, is mainly related to the Hopf bifurcation and hyperpolarization-active caution ($ I_h $) current. In the present paper, the emerging condition for the PIR spike is extended to the bifurcation of the big homoclinic (BHom) orbit in a model without $ I_h $ current. The threshold curve for a spike evoked from a mono-stable or coexisting steady state surrounds the steady state from left, to below, and to right, because the BHom orbit is big enough to surround the steady state. The right part of the threshold curve coincides with the stable manifold of the saddle and acts the threshold for the spike induced by the positive stimulation, resembling that of the saddle-node bifurcation on an invariant cycle, and the left part acts the threshold for the PIR spike, resembling that of the Hopf bifurcation. The bifurcation curve and a codimension-2 bifurcation point related to the BHom orbit are acquired in the two-parameter plane. The results present a comprehensive viewpoint to the dynamics near the BHom orbit bifurcation, which presents a novel threshold curve and extends the conditions for the PIR spike.
[1] | Local control of postinhibitory rebound spiking in CA1 pyramidal neuron dendrites. J. Neurosci. (2010) 30: 6434-6442. |
[2] | Dynamics and bifurcations in a silicon neuron. IEEE Trans. Biomed. Circuits Syst. (2010) 4: 320-328. |
[3] | Postinhibitory rebound during locomotor-like activity in neonatal rat motoneurons in vitro. J. Neurophysiol. (1998) 79: 342-351. |
[4] | P. Channell, G. Cymbalyuk and A. Shilnikov, Origin of bursting through homoclinic spike adding in a neuron model, Phys. Rev. Lett., 98 (2007), 134101. doi: 10.1103/PhysRevLett.98.134101 |
[5] | Dynamics of neurons in the pre-Bötzinger complex under magnetic flow effect. Nonlinear Dyn. (2018) 94: 1961-1971. |
[6] | Bursting and two-parameter bifurcation in the Chay neuronal model. Discrete Contin. Dyn. Syst. Ser. B (2011) 16: 445-456. |
[7] | B. Ermentrout, Simulating, Analyzing, and Animating Dynamical Systems. A Guide to XPPAUT for Researchers and Students, SIAM, Philadelphia, 2002. doi: 10.1137/1.9780898718195 |
[8] | Improving control effects of absence seizures using single-pulse alternately resetting stimulation (SARS) of corticothalamic circuit. Appl. Math. Mech-Engl. (2020) 41: 1287-1302. |
[9] | Mathematical models of threshold phenomena in the nerve membrane. Bull. Math. Biophys. (1955) 17: 257-278. |
[10] | Synchronization and rhythmic processes in physiology. Nature (2001) 410: 277-284. |
[11] | Slow and persistent postinhibitory rebound acts as an intrinsic short-term memory mechanism. J. Neurosci. (2010) 30: 4687-4692. |
[12] | L. Guan, B. Jia and H. Gu, A novel threshold across which negative stimulation evokes action potential near a saddle-node bifurcation in a neuronal model with $I_h$ current, Internat. J. Bifur. Chaos Appl. Sci. Engrg., 29 (2019), 1950198, 26 pp. doi: 10.1142/S0218127419501980 |
[13] | Homoclinic orbits of the FitzHugh-Nagumo equation: The singular-limit. Discrete Contin. Dyn. Syst. Ser. S (2009) 2: 851-872. |
[14] | Global Hopf bifurcation analysis of a six-dimensional FitzHugh-Nagumo neural network with delay by a synchronized scheme. Discrete Contin. Dyn. Syst. Ser. B (2011) 16: 457-474. |
[15] | The local electric changes associated with repetitive action in a non-medullated axon. J. Physiol. (1948) 107: 165-181. |
[16] | Neural excitability, spiking and bursting. Internat. J. Bifur. Chaos Appl. Sci. Engrg. (2000) 10: 1171-1266. |
[17] | Which model to use for cortical spiking neurons?. IEEE T. Neural Network (2004) 15: 1063-1070. |
[18] | (2007) Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting.MIT press. |
[19] | Dynamics analysis of the hippocampal neuronal model subjected to cholinergic action related with alzheimer's disease. Cogn. Neurodyn. (2020) 14: 483-500. |
[20] | T. Malashchenko, A. Shilnikov and G. Cymbalyuk, Bistability of bursting and silence regimes in a model of a leech heart interneuron, Phys. Rev. E, 84 (2011), 041910. doi: 10.1103/PhysRevE.84.041910 |
[21] | Dynamic transitions of the Fitzhugh-Nagumo equations on a finite domain. Discrete Contin. Dyn. Syst. Ser. B (2018) 23: 3935-3947. |
[22] | J. Mitry, M. McCarthy, N. Kopell and M. Wechselberger, Excitable neurons, firing threshold manifolds and canards, J. Math. Neurosci., 3 (2013), Art. 12, 32 pp. doi: 10.1186/2190-8567-3-12 |
[23] | J. Rinzel and G. Ermentrout, Analysis of neural excitability and oscillations, in Methods in Neuronal Modeling (eds. C. Koch and I. Segev), MIT Press, Cambridge, (1998), Chap. 7,251–291. |
[24] | The potassium A-current, low firing rates and rebound excitation in Hodgkin-Huxley models. Bull. Math. Biol. (1995) 57: 899-929. |
[25] | Codimension-two bursting analysis in the delayed neural system with external stimulations. Nonlinear Dyn. (2012) 67: 309-328. |
[26] | Extrinsic modulation and motor pattern generation in a feeding network: A cellular study. J. Neurosci. (2001) 21: 1767-1778. |
[27] | Synchronization transitions induced by partial time delay in an excitatory-inhibitory coupled neuronal network. Nonlinear Dyn. (2017) 89: 2509-2520. |
[28] | A. Tonnelier, Threshold curve for the excitability of bidimensional spiking neurons, Phys. Rev. E, 90(2014), 022701. doi: 10.1103/PhysRevE.90.022701 |
[29] | F. Q. Wu and H. Gu, Bifurcations of negative responses to positive feedback current mediated by memristor in neuron model with bursting patterns, Internat. J. Bifur. Chaos Appl. Sci. Engrg., 30 (2020), 2030009, 22 pp. doi: 10.1142/S0218127420300098 |
[30] | F. Wu, H. Gu and Y. Li, Inhibitory electromagnetic induction current induced enhancement instead of reduction of neural bursting activities, Commun. Nonlinear Sci. Numer. Simul., 79 (2019), 104924, 15 pp. doi: 10.1016/j.cnsns.2019.104924 |
[31] | Stochastic dynamics of conduction failure of action potential along nerve fiber with hopf bifurcation. Sci. China Technol. Sci. (2019) 62: 1502-1511. |
[32] | A feasible neuron for estimating the magnetic field effect. Nonlinear Dyn. (2020) 102: 1849-1867. |
[33] | Bifurcation analysis of bursting solutions of two Hindmarsh-Rose neurons with joint electrical and synaptic coupling. Discrete Contin. Dyn. Syst. Ser. B (2011) 16: 637-651. |
[34] | Z. Zhao and H. Gu, Transitions between classes of neuronal excitability and bifurcations induced by autapse, Sci. Rep., 7 (2017), 6760. doi: 10.1038/s41598-017-07051-9 |
[35] | Z. Zhao, L. Li and H. Gu, Dynamical mechanism of hyperpolarization-activated non-specific cation current induced resonance and spike-timing precision in a neuronal model, Front. Cell. Neurosci., 12 (2018), 62. doi: 10.3389/fncel.2018.00062 |