- 1School of Science, Xuchang University, Xuchang, China
- 2School of Mathematics and Statistics, North China University of Water Resources and Electric Power, Zhengzhou, China
- 3School of Mathematics, Northwest University, Xi’an, China
- 4School of Mathematics and Statistics, Northwestern Polytechnical University, Xi’an, China
Neuron modelling helps to understand the brain behavior through the interaction between neurons, but its mechanism remains unclear. In this paper, the spatiotemporal patterns is investigated in a general networked Hindmarsh-Rose (HR) model. The stability of the network-organized system without delay is analyzed to show the effect of the network on Turing instability through the Hurwitz criterion, and the conditions of Turing instability are obtained. Once the analysis of the zero-delayed system is completed, the critical value of the delay is derived to illustrate the profound impact of the given network on the collected behaviors. It is found that the difference between the collected current and the outgoing current plays a crucial role in neuronal activity, which can be used to explain the generation mechanism of the short-term memory. Finally, the numerical simulation is presented to verify the proposed theoretical results.
Introduction
Neuron modelling plays a vital role in understanding the brain behavior with its dynamic mechanisms Sun et al. (2018); Wang et al. (2021a),Wang et al. (2021b), Wang et al. (2022); Perc (2006). The HR model is the extended Fitzhugh–Nagumo model, which could exhibit various firing activities (periodic, chaotic, bursting firings, etc.) of the membrane potential Hindmarsh and Rose (1982), Hindmarsh and Rose (1984); Gu and Pan (2015). The dynamic bifurcation of the HR model was calculated and analyzed to illustrate the effects of different parameters on the dynamic behavior Wu et al. (2016). And the HR model with electromagnetic induction by analyzing the stability of equilibrium points was studied to show that electromagnetic induction could regularize chaotic regimes by the number of spikes Goulefack et al. (2022). The HR model with the slow intracellular exchange of calcium ions was investigated to present the effects of the coupling strength and the forcing current on the network behavior Rajagopal et al. (2019). Then the role of delay in the bifurcation behavior was considered in the fractional-order HR model with delay, and the condition for the existence of Hopf bifurcation was present Shi et al. (2020). Also, noise-induced resonances were obtained to describe the response of the HR model to noisy signals and intrinsic oscillations Baltanas and Casado (2002). A networked-organized HR model with delay was developed to explain the spatial relations between neurons and model the process of neural migration Lepek and Fronczak (2018). Although the phenomenon of spatiotemporal patterns was often considered to show the firing activities in a general networked HR model with the delay Umesh and Ambika (2021); Shi et al. (2022); Santos et al. (2017), the interaction mechanism of network nodes remains unclear.
Turing instability is a kind of collective spatial behavior, and the pattern formation could describe the interaction between species in a reaction-diffusion system and explain the biological mechanism (the formation of zebra and zebrafish skin, the nonlocal synaptic interactions) Turing (1990); Kolodina et al. (2021); Kondo et al. (2021). A nonlocal evolution equation was developed to explain color patterns on a guppy fish skin through the bifurcation theory Cygan (2021). The effect of landscape heterogeneity on pattern formation is studied to find the mechanisms of diffusion-driven instabilities in the predator-prey interactions Zaker et al. (2021). And the firing behavior of neurons can be represented by the pattern formation in the network-organized FitzHugh-Nagumo model Zheng and Shen (2020). The elimination of spiral waves is investigated to prevent some brain disorders by the transition between pattern formations in an HR model Eteme et al. (2019). Pattern formation is a well-studied phenomenon in neural field models, which describes the collected groups of neurons and presents the generation of spatial distribution in a dynamical system Kriener et al. (2013). Sometimes, the collective behavior is easier to keep the network in activity Erichsen and Brunnet (2009). Also, short-term memory is a kind of the collective behavior of neurons Zheng et al. (2020). Meanwhile, the network topology is predominant in the critical features of Turing systems Diego et al. (2018). Therefore, the general networked HR model should be considered to show the dynamical mechanism of neuronal activity.
It is well known that the neurons link each other through a network, and the brain activity is the collective behavior of the neurons rather than a single neuron. Pattern formation and bifurcation Tian et al. (2021); Yang (2022); Ma et al. (2021) is a crucial tool to elaborate on the dynamic and biological mechanism of the collective behavior of the neurons. In this paper, an HR model with a random network is considered to show the spatiotemporal patterns of collective behaviors. The effect of the network on pattern formation is presented through the Hurwitz criterion, and the conditions of Turing instability are derived. Then, Hopf bifurcation illustrates the profound impact of network and delay on the collected behavior. It is found that the difference in network topology plays a vital role in neuronal activity. Finally, the numerical simulation explains the generation mechanism of the short-term memory.
Model Description
In this paper, we consider the following HR model on networks,
where xi, yi and zi (i = 1, … , n) denote the membrane potential, recovery variable and bursting variable at node (neuron) i, respectively. I is the external stimulus, A is the adjacent matrix of the interaction between nodes, f (xj, xi), g (yj, yi) is the diffusive function. The equilibrium point of system (1) can be derived from the following Equation 2
The linear term of system (1) at (x0, y0, z0) can be expressed as
where Aij is the symmetric adjacent matrix to show the interaction of nodes on network. xj − qxi, yj − qyi is the linear part of f (xj, xi), g (yj, yi), respectively. q can be treated as the difference between the collected current and the outgoing current.
A general solution of system (3) is
Substituting system (4) into system (3), one has
where
Then the characteristic equation is
where
When τ = 0, the characteristic equation is,
According to Hurwitz criterion, the sufficient and necessary conditions for stable system (1) are,
The converse-negative proposition is the condition of Turing instability, namely, Turing instability induced by network when one of Hi(i = 1, 2, 3) holds.
More precisely, for H1, Turing instability occurs when
where
According to the properties of quadratic equations of one variable, the minimum value of p1 (Λk) is
For H2, Turing instability occurs when
where, the analysis process refers to p1 (Λk), and
For H3, Turing instability occurs when
where
Suppose r1, r2 (r1 ≤ r2) is the roots of
where p3 (r1), p3 (r2) are local maximum and minimum respectively. According to the properties of cubic equations of one variable, Turing instability occurs when r1 < 0 and p3 (r1) > 0; Turing instability occurs when r1 > 0 and p3 (0) > 0.
Finally, we consider the stability of system (1), namely,
Suppose that the pure imaginary root λ = j ω (j represents the imaginary unit) exists and we substitute it into the above characteristic equation, we have
Separating the real and imaginary parts
Solving cos (ωτ), sin (ωτ) to get
Due to cos (ωτ)2 + sin (ωτ)2 = 1, one has
where
If a positive root exists xi (i = 1, 2, 3, 4) at least in system (9), one has
where
The transversality condition
where
If
Turing instability occurs when τ0 > τ > τc. τ0 > τ make system (1) without network stable, and network induce Turing instability, namely, τ0 > τ > τc.
Through the above analysis, we can draw the following sufficient conditions for Turing instability.
Theorem 1. Turing instability occurs in the network-organized system when a Hi (i = 1, 2, 3) holds; Turing instability occurs in the network-organized system with delay when τ ∈ (τ0, τc).
Proof 1. The proof process can refer to the above analysis.
Results and Discussion
In this section, these parameters a = 1, b = 3, c = 1, d = 5, r = 0.01, s = 4, I = 1, xr = −1.6 Hindmarsh and Rose (1984) and d1 = 0.1, d2 = 0.3 are set. The adjacent matrix A is generated by random network Zheng and Shen (2020); Erdos and Renyi (1959) and the link probability p and nodes n = 100 is the initial value. For node i and node j, if the random number
From Figure 1, system (1) is stable when d1 = 0, d2 = 0, τ = 0, which is the precondition of Turing instability induced by network and delay. Namely, Hi > 0 when Λi = 0 (Figure 1A,B,C) and the current values below the activity threshold of the single neuron. Also, Λi can not lead to Turing instability when Λi < 0, which means it is difficult to keep the network in activity when the collected current is the same as the outgoing current. Let’s take H3 (H1 and H2 are relatively simple, and H3 is more representative.) as an example to illustrate the stability of system (1) (Figure 2–6). Although system (1) without delay is periodic behavior (Figure 2B) under the initial stimulus, it ultimately tends to the rest state (Figure 2A) when no stimulus is added. Short-term memory attributes to a fixed point attractor Wang (2001), which could persistent neuronal activity when the remembered stimulus is removed Goldman (2009). The above results bring into correspondence with our results (Figure 2). When the collected current is larger than the outgoing currentq = 0.96, H3 < 0 holds (Figure 3A) and Turing instability occurs (Figure 3B). In general, most of the neurons will be in the rest state. Meanwhile, a few persistent neuronal activities because of the constant external stimuli (Figure 3B). When the link probability p increases and the difference q becomes larger, H3 < 0 holds (Figure 4A) and the corresponding pattern formation (Figure 4B) shows the periodic neural activity. Namely, constant external stimulation from other neurons is necessary to keep the neural activity. For the short-term memory, noise could induce the switch of different memories Zheng et al. (2020), and the constant external stimulation from other neurons can also leads to the switch of different memories. Only the link probability could not causes instability (Figure 5A), but it could keep the neural activity longer (Figure 5B). When the difference between the collected current and the outgoing current exceeds the threshold value of the neural activity (Figure 6A), Turing instability occurs (Figure 6B). Namely, the short-term memory generation requires enough stimulation, which is in line with the actual situation. Also, these above results can be verified by the bifurcation (Figure 7):system (1) is always stable when q = 0 (Figure 7A) and q could induce Turing instability (Figure 7B).
FIGURE 1. The distribution of Hi(i = 1, 2, 3) when a = 1, b = 3, c = 1, d = 5, r = 0.01, s = 4, I = 1, xr = −1.6, d1 = 0.1, d2 = 0.3, τ = 0. (A) Turing instability occurs when H1 < 0, or else stable. (B) Turing instability occurs when H2 < 0, or else stable. (C) Turing instability occurs when H3 < 0, or else stable. (D) The bifurcation about τ when d1 = d2 = 0.
FIGURE 2. The stability of system (1) when p = 0.05, τ = 0. (A) The distribution of eigenvalues Λi (labeled by *) when q = 0. (B) The corresponding pattern formation is stable in system (1) when no eigenvalue Λi make H3 < 0 hold.
FIGURE 3. The stability of system (1) when p = 0.05, q = 0.96, τ = 0. (A) The distribution of eigenvalues Λi (labeled by *). (B) Turing instability occurs in system (1) when some eigenvalue Λi make H3 < 0 hold.
FIGURE 4. The stability of system (1) when p = 0.5, q = 0.96, τ = 0. (A) The distribution of eigenvalues Λi (labeled by *). (B) Turing instability occurs in system (1) when a eigenvalue Λi make H3 < 0 hold.
FIGURE 5. The stability of system (1) when p = 0.1, q = 0, τ = 0. (A) The distribution of eigenvalues Λi (labeled by *). (B) The corresponding pattern formation is stable in system (1) when no eigenvalue Λi make H3 < 0 hold.
FIGURE 6. The stability of system (1) when p = 0.1, q = 0.82, τ = 0. (A) The distribution of eigenvalues Λi (labeled by *). (B) Turing instability occurs in system (1) when some eigenvalue Λi make H3 < 0 hold.
FIGURE 7. The bifurcation about p, q when τ = 0. (A) The inexistence bifurcation about p when q = 0. (B) The bifurcation about q when p = 0.1.
Delay plays a vital role in the rise or fall of neural activity. From Figure 8, the dynamic behavior (Figure 8A) is different from system (1) without delay (Figure 2), but the link probability nearly has no effect on the stability of system (1) (Figure 8B), because the collected current always equal to the outgoing current when q = 0. When q ≠ 0, the periodic behavior of the neuron may occur (Figure 8C,D). Meanwhile, the periodic behavior of the neuron under different p (Figure 8C,D). It is found that the link probability p could make system (1) network synchronization and the synchronicity becomes stronger with the increase of p, but the neuronal activity is low (Figure 9A,B,C, Figure 10A). Of course, q could increases neuronal activity (Figure 9D), which can also be represented in the bifurcation (Figure 10B). The phenomenon of spatiotemporal patterns Lepek and Fronczak (2018); Umesh and Ambika (2021); Shi et al. (2022); Santos et al. (2017) was often considered in the network-organized system when q = 1. The difference q < 1 between the collected current and the outgoing current plays a vital role in the neuronal activity through the above analysis, which may further show the firing mechanism in a general networked HR model with delay. Meanwhile, short-term memory results from external stimuli Zheng et al. (2020), but not all stimuli result in short-term memory. The difference q < 1 may theoretically explain why some stimuli can’t lead to short-term memory.
FIGURE 8. The pattern formation of system (1) when τ = 65. (A) The pattern formation when p = 0.05, q = 0. (B) The pattern formation when p = 0.1, q = 0. (C) The pattern formation when p = 0.1, q = 0.96. (D) The pattern formation when p = 0.15, q = 0.96.
FIGURE 9. The pattern formation of system (1) when τ = 68. (A) The pattern formation when p = 0, q = 0. (B) The pattern formation when p = 0.01, q = 0. (C) The pattern formation when p = 0.1, q = 0. (D) The pattern formation when p = 0.1, q = 0.96.
FIGURE 10. The bifurcation about p, q when τ = 65. (A) The inexistence bifurcation about p when q = 0. (B) The bifurcation about q when p = 0.1.
Conclusion
In this paper, spatiotemporal patterns are investigated to illustrate the collected behavior of neurons and the generation mechanism of short-term memory. We obtain the algebraic expressions for Turing instability to occur in any HR network setup. Then, we derive the critical value of Hopf bifurcation to present a profound impact of both network and delay on the Turing instability. Also, we find that the collected current and outgoing current play a vital role in neuronal activity, especially in the generation mechanism of the short-term memory. Meanwhile, the collected behavior may fire when the input to the neuron is below a certain threshold, and the output reaches a stationary regime. Finally, we try to explain the generation mechanism of the short-term memory through the theoretical results and numerical simulation.
Data Availability Statement
The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author (Some related code can be found: https://github.com/zhengqianqian35).
Author Contributions
YX and JS conceived of the study, designed the study. The remaining authors contributed to refining the ideas, carrying out additional analyses and finalizing this paper. All authors contributed to the writing and revisions.
Funding
This work is supported by National Natural Science Foundation of China (12002297), Basic research Project of Universities in Henan Province (21zx009), Program for Science& Technology Innovation Talents in Universities of Henan Province (22HASTIT018), The Scientific Research Innovation Team of Xuchang University (2022CXTD002), Key scientific research projects of Henan Institutions of Higher learning in 2021 (21B130004), Outstanding Young Backbone Teacher of Xuchang University (2022).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Baltanás J. P., Casado J. M. (2002). Noise-induced Resonances in the Hindmarsh-Rose Neuronal Model. Phys. Rev. E Stat. Nonlin Soft Matter Phys. 65, 041915. doi:10.1103/PhysRevE.65.041915
Cygan S. (2021). Pattern Formation in Nonlocal Kondo Model. Math. Meth Appl. Sci. 44, 10804–10822. doi:10.1002/mma.7448
Diego X., Marcon L., Muller P., Sharpe J. (2018). Key Features of Turing Systems Are Determined Purely by Network Topology. Phys. Rev. X 8, 021071. doi:10.1103/physrevx.8.021071
Erichsen R., Brunnet L. G. (2009). Multistability in Networks of Hindmarsh-Rose Neurons. Phys. Rev. E Stat. Nonlin Soft Matter Phys. 78, 061917. doi:10.1103/PhysRevE.78.061917
Etémé A. S., Tabi C. B., Mohamadou A., Kofané T. C. (2019). Elimination of Spiral Waves in a Two-Dimensional Hindmarsh-Rose Neural Network under Long-Range Interaction Effect and Frequency Excitation. Phys. A Stat. Mech. its Appl. 533, 122037. doi:10.1016/j.physa.2019.122037
Goldman M. S. (2009). Memory without Feedback in a Neural Network. Neuron 61, 621–634. doi:10.1016/j.neuron.2008.12.012
Goulefack L. M., Chamgoue A. C., Anteneodo C., Yamapi R. (2022). Stability Analysis of the Hindmarsh-Rose Neuron under Electromagnetic Induction. Nonlinear Dyn. 108, 2627–2642. doi:10.1007/s11071-022-07331-8
Gu H., Pan B. (2015). A Four-Dimensional Neuronal Model to Describe the Complex Nonlinear Dynamics Observed in the Firing Patterns of a Sciatic Nerve Chronic Constriction Injury Model. Nonlinear Dyn. 81, 2107–2126. doi:10.1007/s11071-015-2129-7
Hindmarsh J. L., Rose R. M. (1984). A Model of Neuronal Bursting Using Three Coupled First Order Differential Equations. Proc. R. Soc. Lond B Biol. Sci. 221, 87–102. doi:10.1098/rspb.1984.0024
Hindmarsh J. L., Rose R. M. (1982). A Model of the Nerve Impulse Using Two First-Order Differential Equations. Nature 296, 162–164. doi:10.1038/296162a0
Kolodina K., Wyller J., Oleynik A., Sørensen M. P. (2021). Pattern Formation in a 2-population Homogenized Neuronal Network Model. J. Math. Neurosci. 11, 9–38. doi:10.1186/s13408-021-00107-1
Kondo S., Watanabe M., Miyazawa S. (2021). Studies of Turing Pattern Formation in Zebrafish Skin. Philosophical Trans. R. Soc. A 379, 20200274. doi:10.1098/rsta.2020.0274
Kriener B., Helias M., Rotter S., Diesmann M., Einevoll G. T. (2013). How Pattern Formation in Ring Networks of Excitatory and Inhibitory Spiking Neurons Depends on the Input Current Regime. Front. Comput. Neurosci. 7, 187. doi:10.3389/fncom.2013.00187
Lepek M., Fronczak P. (2018). Spatial Evolution of Hindmarsh-Rose Neural Network with Time Delays. Nonlinear Dyn. 92, 751–761. doi:10.1007/s11071-018-4088-2
Ma J., Xu Y., Li Y., Tian R., Ma S., Kurths J. (2021). Quantifying the Parameter Dependent Basin of the Unsafe Regime of Asymmetric Lévy-Noise-Induced Critical Transitions. Appl. Math. Mech.-Engl. Ed. 42, 65–84. doi:10.1007/s10483-021-2672-8
Rajagopal K., Khalaf A. J. M., Parastesh F., Moroz I., Karthikeyan A., Jafari S. (2019). Dynamical Behavior and Network Analysis of an Extended Hindmarsh-Rose Neuron Model. Nonlinear Dyn. 98, 477–487. doi:10.1007/s11071-019-05205-0
Santos M., Szezech J., Borges F., Iarosz K., Caldas I., Batista A., et al. (2017). Chimera-like States in a Neuronal Network Model of the Cat Brain. Chaos, Solit. Fractals 101, 751–761. doi:10.1016/j.chaos.2017.05.028
Shi M., Yu Y., Xu Q. (2020). Hopf Bifurcation of the Fractional-Order Hindmarsh-Rose Neuron Model with Time-Delay. Rocky Mt. J. Math. 50, 2213–2222. doi:10.1216/rmj.2020.50.2213
Shi X., Wang Z., Wang Z., Zhuang L. (2022). Spatiotemporal Pattern in a Neural Network with Non-smooth Memristor. era 30, 715–731. doi:10.3934/era.2022038
Sun X., Perc M., Kurths J., Lu Q. (2018). Fast Regular Firings Induced by Intra- and Inter-time Delays in Two Clustered Neuronal Networks. Chaos 28, 106310. doi:10.1063/1.5037142
Tian R., Zhao Z., Xu Y. (2021). Variable Scale-Convex-Peak Method for Weak Signal Detection. Sci. China Technol. Sci. 64, 331–340. doi:10.1007/s11431-019-1530-4
Turing A. (1990). The Chemical Basis of Morphogenesis. Bull. Math. Biol. 52, 153–197. doi:10.1016/s0092-8240(05)80008-4
Umesh K. V., Ambika G. (2021). Emergent Dynamics and Spatio Temporal Patterns on Multiplex Neuronal Networks. Front. Comput. Neurosci. 15, 774969. doi:10.3389/fncom.2021.774969
Wang X.-J. (2001). Synaptic Reverberation Underlying Mnemonic Persistent Activity. Trends Neurosci. 24, 455–463. doi:10.1016/s0166-2236(00)01868-3
Wang Z., Li Y., Xu Y., Kapitaniak T., Kurths J. (2022). Coherence-resonance Chimeras in Coupled Hr Neurons with Alpha-Stable Levy Noise. J. Stat. Mech. Theory Exp. 2022, 1–15. doi:10.1088/1742-5468/ac6254
Wang Z., Xu Y., Li Y., Kapitaniak T., Kurths J. (2021a). Chimera States in Coupled Hindmarsh-Rose Neurons with α-stable Noise. Chaos, Solit. Fractals 148, 110976. doi:10.1016/j.chaos.2021.110976
Wang Z., Xu Y., Li Y., Kurths J. (2021b). The Probability Density Function of Interspike Intervals in an Fhn Model with α-stable Noise. Eur. Phys. J. Plus 136, 1–14. doi:10.1140/epjp/s13360-021-01245-x
Wu K., Luo T., Lu H., Wang Y. (2016). Bifurcation Study of Neuron Firing Activity of the Modified Hindmarsh-Rose Model. Neural Comput. Applic 27, 739–747. doi:10.1007/s00521-015-1892-1
Yang W. (2022). Bifurcation and Dynamics in Double-Delayed Chua Circuits with Periodic Perturbation. Chin. Phys. B 31, 020201. doi:10.1088/1674-1056/ac1e0b
Zaker N., Cobbold C. A., Lutscher F. (2021). The Effect of Landscape Fragmentation on Turing-Pattern Formation. Math. Biosci. Eng. 19, 2506–2537. doi:10.3934/mbe.2022116
Zheng Q., Shen J. (2020). Turing Instability Induced by Random Network in Fitzhugh-Nagumo Model. Appl. Math. Comput. 381, 125304. doi:10.1016/j.amc.2020.125304
Keywords: HR model, pattern formation, network, matrix, turing instability, delay
Citation: Zheng Q, Shen J, Zhang R, Guan L and Xu Y (2022) Spatiotemporal Patterns in a General Networked Hindmarsh-Rose Model. Front. Physiol. 13:936982. doi: 10.3389/fphys.2022.936982
Received: 05 May 2022; Accepted: 03 June 2022;
Published: 28 June 2022.
Edited by:
Kelly Cristiane Iarosz, University of São Paulo, BrazilReviewed by:
José Trobia, Universidade Estadual de Ponta Grossa, BrazilArthur Valencio, University of Campinas, Brazil
Copyright © 2022 Zheng, Shen, Zhang, Guan and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Xu, aHN1eDNAbndwdS5lZHUuY24=; Jianwei Shen, eGNqd3NoZW5AZ21haWwuY29t