Equilibrium analysis and phase synchronization of two coupled HR neurons with gap junction

Research on cascading high-dimensional isomorphic chaotic maps

In order to overcome the security weakness of the discrete chaotic sequence caused by small Lyapunov exponent and keyspace, a general chaotic construction method by cascading multiple high-dimensional isomorphic maps is presented in this paper. Compared with the original map, the parameter space of the resulting chaotic map is enlarged many times. Moreover, the cascaded system has larger chaotic domain and bigger Lyapunov exponents with proper parameters. In order to evaluate the effectiveness of the presented method, the generalized 3-D Hénon map is utilized as an example to analyze the dynamical behaviors under various cascade modes. Diverse maps are obtained by cascading 3-D Hénon maps with different parameters or different permutations. It is worth noting that some new dynamical behaviors, such as coexisting attractors and hyperchaotic attractors are also discovered in cascaded systems. Finally, an application of image encryption is delivered to demonstrate the excellent performance of the obtained chaotic sequences.

Cluster burst synchronization in a scale-free network of inhibitory bursting neurons

We consider a scale-free network of inhibitory Hindmarsh-Rose (HR) bursting neurons, and make a computational study on coupling-induced cluster burst synchronization by varying the average coupling strength J 0 . For sufficiently small J 0 , non-cluster desynchronized states exist. However, when passing a critical point J c ∗ ( ≃ 0.16 ) , the whole population is segregated into 3 clusters via a constructive role of synaptic inhibition to stimulate dynamical clustering between individual burstings, and thus 3-cluster desynchronized states appear. As J 0 is further increased and passes a lower threshold J l ∗ ( ≃ 0.78 ) , a transition to 3-cluster burst synchronization occurs due to another constructive role of synaptic inhibition to favor population synchronization. In this case, HR neurons in each cluster make burstings every 3rd cycle of the instantaneous burst rate R w ( t ) of the whole population, and exhibit burst synchronization. However, as J 0 passes an intermediate threshold J m ∗ ( ≃ 5.2 ) , HR neurons fire burstings intermittently at a 4th cycle of R w ( t ) via burst skipping rather than at its 3rd cycle, and hence they begin to make intermittent hoppings between the 3 clusters. Due to such intermittent intercluster hoppings via burst skippings, the 3 clusters become broken up (i.e., the 3 clusters are integrated into a single one). However, in spite of such break-up (i.e., disappearance) of the 3-cluster states, (non-cluster) burst synchronization persists in the whole population, which is well visualized in the raster plot of burst onset times where bursting stripes (composed of burst onset times and indicating burst synchronization) appear successively. With further increase in J 0 , intercluster hoppings are intensified, and bursting stripes also become dispersed more and more due to a destructive role of synaptic inhibition to spoil the burst synchronization. Eventually, when passing a higher threshold J h ∗ ( ≃ 17.8 ) a transition to desynchronization occurs via complete overlap between the bursting stripes. Finally, we also investigate the effects of stochastic noise on both 3-cluster burst synchronization and intercluster hoppings.

Nonlinear optimal control for the synchronization of biological neurons under time-delays

The article proposes a nonlinear optimal control method for synchronization of neurons that exhibit nonlinear dynamics and are subject to time-delays. The model of the Hindmarsh-Rose (HR) neurons is used as a case study. The dynamic model of the coupled HR neurons undergoes approximate linearization around a temporary operating point which is recomputed at each iteration of the control method. The linearization procedure relies on Taylor series expansion of the model and on computation of the associated Jacobian matrices. For the approximately linearized model of the coupled HR neurons an H-infinity controller is designed. For the selection of the controller's feedback gain an algebraic Riccati equation is repetitively solved at each time-step of the control algorithm. The stability properties of the control loop are proven through Lyapunov analysis. First, it is shown that the H-infinity tracking performance criterion is satisfied. Moreover, it is proven that the control loop is globally asymptotically stable.

Burst synchronization in a scale-free neuronal network with inhibitory spike-timing-dependent plasticity

We are concerned about burst synchronization (BS), related to neural information processes in health and disease, in the Barabási-Albert scale-free network (SFN) composed of inhibitory bursting Hindmarsh-Rose neurons. This inhibitory neuronal population has adaptive dynamic synaptic strengths governed by the inhibitory spike-timing-dependent plasticity (iSTDP). In previous works without considering iSTDP, BS was found to appear in a range of noise intensities for fixed synaptic inhibition strengths. In contrast, in our present work, we take into consideration iSTDP and investigate its effect on BS by varying the noise intensity. Our new main result is to find occurrence of a Matthew effect in inhibitory synaptic plasticity: good BS gets better via LTD, while bad BS get worse via LTP. This kind of Matthew effect in inhibitory synaptic plasticity is in contrast to that in excitatory synaptic plasticity where good (bad) synchronization gets better (worse) via LTP (LTD). We note that, due to inhibition, the roles of LTD and LTP in inhibitory synaptic plasticity are reversed in comparison with those in excitatory synaptic plasticity. Moreover, emergences of LTD and LTP of synaptic inhibition strengths are intensively investigated via a microscopic method based on the distributions of time delays between the pre- and the post-synaptic burst onset times. Finally, in the presence of iSTDP we investigate the effects of network architecture on BS by varying the symmetric attachment degree l ∗ and the asymmetry parameter Δ l in the SFN.

Effect of spike-timing-dependent plasticity on stochastic burst synchronization in a scale-free neuronal network

We consider an excitatory population of subthreshold Izhikevich neurons which cannot fire spontaneously without noise. As the coupling strength passes a threshold, individual neurons exhibit noise-induced burstings. This neuronal population has adaptive dynamic synaptic strengths governed by the spike-timing-dependent plasticity (STDP). However, STDP was not considered in previous works on stochastic burst synchronization (SBS) between noise-induced burstings of sub-threshold neurons. Here, we study the effect of additive STDP on SBS by varying the noise intensity D in the Barabási-Albert scale-free network (SFN). One of our main findings is a Matthew effect in synaptic plasticity which occurs due to a positive feedback process. Good burst synchronization (with higher bursting measure) gets better via long-term potentiation (LTP) of synaptic strengths, while bad burst synchronization (with lower bursting measure) gets worse via long-term depression (LTD). Consequently, a step-like rapid transition to SBS occurs by changing D, in contrast to a relatively smooth transition in the absence of STDP. We also investigate the effects of network architecture on SBS by varying the symmetric attachment degree [Formula: see text] and the asymmetry parameter [Formula: see text] in the SFN, and Matthew effects are also found to occur by varying [Formula: see text] and [Formula: see text]. Furthermore, emergences of LTP and LTD of synaptic strengths are investigated in details via our own microscopic methods based on both the distributions of time delays between the burst onset times of the pre- and the post-synaptic neurons and the pair-correlations between the pre- and the post-synaptic instantaneous individual burst rates (IIBRs). Finally, a multiplicative STDP case (depending on states) with soft bounds is also investigated in comparison with the additive STDP case (independent of states) with hard bounds. Due to the soft bounds, a Matthew effect with some quantitative differences is also found to occur for the case of multiplicative STDP.

Dynamical Analysis of the Hindmarsh-Rose Neuron With Time Delays

This brief is mainly concerned with a series of dynamical analyses of the Hindmarsh-Rose (HR) neuron with state-dependent time delays. The dynamical analyses focus on stability, Hopf bifurcation, as well as chaos and chaos control. Through the stability and bifurcation analysis, we determine that increasing the external current causes the excitable HR neuron to exhibit periodic or chaotic bursting/spiking behaviors and emit subcritical Hopf bifurcation. Furthermore, by choosing a fixed external current and varying the time delay, the stability of the HR neuron is affected. We analyze the chaotic behaviors of the HR neuron under a fixed external current through time series, bifurcation diagram, Lyapunov exponents, and Lyapunov dimension. We also analyze the synchronization of the chaotic time-delayed HR neuron through nonlinear control. Based on an appropriate Lyapunov-Krasovskii functional with triple integral terms, a nonlinear feedback control scheme is designed to achieve synchronization between the uncontrolled and controlled models. The proposed synchronization criteria are derived in terms of linear matrix inequalities to achieve the global asymptotical stability of the considered error model under the designed control scheme. Finally, numerical simulations pertaining to stability, Hopf bifurcation, periodic, chaotic, and synchronized models are provided to demonstrate the effectiveness of the derived theoretical results.

Route searching based on neural networks and heuristic reinforcement learning

In this paper, an improved and much stronger RNH-QL method based on RBF network and heuristic Q-learning was put forward for route searching in a larger state space. Firstly, it solves the problem of inefficiency of reinforcement learning if a given problem's state space is increased and there is a lack of prior information on the environment. Secondly, RBF network as weight updating rule, reward shaping can give an additional feedback to the agent in some intermediate states, which will help to guide the agent towards the goal state in a more controlled fashion. Meanwhile, with the process of Q-learning, it is accessible to the underlying dynamic knowledge, instead of the need of background knowledge of an upper level RBF network. Thirdly, it improves the learning efficiency by incorporating the greedy exploitation strategy to train the neural network, which has been testified by the experimental results.

Hopf Bifurcations in Directed Acyclic Networks of Linearly Coupled Hindmarsh-Rose Systems

This paper addresses the existence of Hopf bifurcations in a directed acyclic network of neurons, each of them being modeled by a Hindmarsh-Rose (HR) neuronal model. The bifurcation parameter is the small parameter corresponding to the ratio of time scales between the fast and the slow dynamics. We first prove that, under certain hypotheses, the single uncoupled neuron can undergo a Hopf bifurcation. Hopf bifurcation occurrences in a directed acyclic network of HR neurons are then discussed. Numerical simulations are carried out to observe these bifurcations and to illustrate the theoretical results.

Frequency-domain order parameters for the burst and spike synchronization transitions of bursting neurons

We are interested in characterization of synchronization transitions of bursting neurons in the frequency domain. Instantaneous population firing rate (IPFR) [Formula: see text], which is directly obtained from the raster plot of neural spikes, is often used as a realistic collective quantity describing population activities in both the computational and the experimental neuroscience. For the case of spiking neurons, a realistic time-domain order parameter, based on [Formula: see text], was introduced in our recent work to characterize the spike synchronization transition. Unlike the case of spiking neurons, the IPFR [Formula: see text] of bursting neurons exhibits population behaviors with both the slow bursting and the fast spiking timescales. For our aim, we decompose the IPFR [Formula: see text] into the instantaneous population bursting rate [Formula: see text] (describing the bursting behavior) and the instantaneous population spike rate [Formula: see text] (describing the spiking behavior) via frequency filtering, and extend the realistic order parameter to the case of bursting neurons. Thus, we develop the frequency-domain bursting and spiking order parameters which are just the bursting and spiking "coherence factors" [Formula: see text] and [Formula: see text] of the bursting and spiking peaks in the power spectral densities of [Formula: see text] and [Formula: see text] (i.e., "signal to noise" ratio of the spectral peak height and its relative width). Through calculation of [Formula: see text] and [Formula: see text], we obtain the bursting and spiking thresholds beyond which the burst and spike synchronizations break up, respectively. Consequently, it is shown in explicit examples that the frequency-domain bursting and spiking order parameters may be usefully used for characterization of the bursting and the spiking transitions, respectively.

Noise-induced burst and spike synchronizations in an inhibitory small-world network of subthreshold bursting neurons

We are interested in noise-induced firings of subthreshold neurons which may be used for encoding environmental stimuli. Noise-induced population synchronization was previously studied only for the case of global coupling, unlike the case of subthreshold spiking neurons. Hence, we investigate the effect of complex network architecture on noise-induced synchronization in an inhibitory population of subthreshold bursting Hindmarsh-Rose neurons. For modeling complex synaptic connectivity, we consider the Watts-Strogatz small-world network which interpolates between regular lattice and random network via rewiring, and investigate the effect of small-world connectivity on emergence of noise-induced population synchronization. Thus, noise-induced burst synchronization (synchrony on the slow bursting time scale) and spike synchronization (synchrony on the fast spike time scale) are found to appear in a synchronized region of the [Formula: see text]-[Formula: see text] plane ([Formula: see text]: synaptic inhibition strength and [Formula: see text]: noise intensity). As the rewiring probability [Formula: see text] is decreased from 1 (random network) to 0 (regular lattice), the region of spike synchronization shrinks rapidly in the [Formula: see text]-[Formula: see text] plane, while the region of the burst synchronization decreases slowly. We separate the slow bursting and the fast spiking time scales via frequency filtering, and characterize the noise-induced burst and spike synchronizations by employing realistic order parameters and statistical-mechanical measures introduced in our recent work. Thus, the bursting and spiking thresholds for the burst and spike synchronization transitions are determined in terms of the bursting and spiking order parameters, respectively. Furthermore, we also measure the degrees of burst and spike synchronizations in terms of the statistical-mechanical bursting and spiking measures, respectively.

Neuronal spike initiation modulated by extracellular electric fields

Based on a reduced two-compartment model, the dynamical and biophysical mechanism underlying the spike initiation of the neuron to extracellular electric fields is investigated in this paper. With stability and phase plane analysis, we first investigate in detail the dynamical properties of neuronal spike initiation induced by geometric parameter and internal coupling conductance. The geometric parameter is the ratio between soma area and total membrane area, which describes the proportion of area occupied by somatic chamber. It is found that varying it could qualitatively alter the bifurcation structures of equilibrium as well as neuronal phase portraits, which remain unchanged when varying internal coupling conductance. By analyzing the activating properties of somatic membrane currents at subthreshold potentials, we explore the relevant biophysical basis of spike initiation dynamics induced by these two parameters. It is observed that increasing geometric parameter could greatly decrease the intensity of the internal current flowing from soma to dendrite, which switches spike initiation dynamics from Hopf bifurcation to SNIC bifurcation; increasing internal coupling conductance could lead to the increase of this outward internal current, whereas the increasing range is so small that it could not qualitatively alter the spike initiation dynamics. These results highlight that neuronal geometric parameter is a crucial factor in determining the spike initiation dynamics to electric fields. The finding is useful to interpret the functional significance of neuronal biophysical properties in their encoding dynamics, which could contribute to uncovering how neuron encodes electric field signals.

Synchronization and stochastic resonance of the small-world neural network based on the CPG

According to biological knowledge, the central nervous system controls the central pattern generator (CPG) to drive the locomotion. The brain is a complex system consisting of different functions and different interconnections. The topological properties of the brain display features of small-world network. The synchronization and stochastic resonance have important roles in neural information transmission and processing. In order to study the synchronization and stochastic resonance of the brain based on the CPG, we establish the model which shows the relationship between the small-world neural network (SWNN) and the CPG. We analyze the synchronization of the SWNN when the amplitude and frequency of the CPG are changed and the effects on the CPG when the SWNN's parameters are changed. And we also study the stochastic resonance on the SWNN. The main findings include: (1) When the CPG is added into the SWNN, there exists parameters space of the CPG and the SWNN, which can make the synchronization of the SWNN optimum. (2) There exists an optimal noise level at which the resonance factor Q gets its peak value. And the correlation between the pacemaker frequency and the dynamical response of the network is resonantly dependent on the noise intensity. The results could have important implications for biological processes which are about interaction between the neural network and the CPG.

Dynamic behavior analysis of fractional-order Hindmarsh-Rose neuronal model

Previous experimental work has shown that the firing rate of multiple time-scales of adaptation for single rat neocortical pyramidal neurons is consistent with fractional-order differentiation, and the fractional-order neuronal models depict the firing rate of neurons more verifiably than other models do. For this reason, the dynamic characteristics of the fractional-order Hindmarsh-Rose (HR) neuronal model were here investigated. The results showed several obvious differences in dynamic characteristic between the fractional-order HR neuronal model and an integer-ordered model. First, the fractional-order HR neuronal model displayed different firing modes (chaotic firing and periodic firing) as the fractional order changed when other parameters remained the same as in the integer-order model. However, only one firing mode is displayed in integer-order models with the same parameters. The fractional order is the key to determining the firing mode. Second, the Hopf bifurcation point of this fractional-order model, from the resting state to periodic firing, was found to be larger than that of the integer-order model. Third, for the state of periodically firing of fractional-order and integer-order HR neuron model, the firing frequency of the fractional-order neuronal model was greater than that of the integer-order model, and when the fractional order of the model decreased, the firing frequency increased.