1. Introduction
Neurons are the basic structure and function of the nervous system; they are the foundation of the brain. There are about
neurons in the brain, which can perceive changes in the environment and complete the reception and transmission of information. And these neurons have some crucial properties, such as excitability, conductivity, and plasticity. Therefore, it is necessary to study the dynamical behavior of neurons. Many researchers tend to use mathematical models with parameters to describe neurons and use these models to study the complex nonlinear behavior of neurons [
1]. With the further consideration and improvement of the neuron model by researchers, the Hodgkin–Huxley neuron, Fitzhugh–Nagumo neuron, Morris–Lecar neuron, Chay neuron, Hindmarsh–Rose neuron, Izhikevich neuron, Hopfield neural network, and TrueNorth chip models [
2,
3,
4,
5,
6,
7,
8,
9,
10] were proposed successively, and these models have become major tools for computational neuroscience research.
In recent years, many works have concentrated on the discharge activity and rich nonlinear dynamics of single neurons and neuronal networks [
11,
12,
13]. These studies provided an analytical mechanism for neurons to participate in information coding and energy metabolism during brain operation. Lakshmanan et al. [
14] analyzed the stability, bifurcation, and chaos of the Hindmarsh–Rose neuron model with time delay, it was found that the neuron model can show different discharge behaviors by setting corresponding bifurcation parameters, and the enactment of synchronization criteria can make the system achieve global asymptotic stability. In [
15], the authors took into account the method in which two FHN neuron circuits are coupled with a Josephson junction and realized the synchronization stability of the whole neural circuit by adjusting the parameters of the Josephson junction in the coupling channel. And the investigation results revealed that the synchronization between the neural circuits is determined by the coupling method and also depends on the physical properties of the nonlinear oscillator in the coupling channel. In fact, the activation of magnetic field coupling is the result of the continuous release and propagation of intracellular and extracellular ions. Therefore, magnetic field coupling may play a crucial role in modulating collective behavior among neurons.
Importantly, the fluctuation in ion concentration inside and outside the cell membrane would cause changes in electromagnetic field distribution inside and outside the cell, resulting in electromagnetic induction. The influence of the electromagnetic effect should be considered in many investigations of neurons [
16,
17]. A memristor is a kind of circuit element with memory characteristics which is defined by magnetic flux and charge. A memristor can be used to simulate the synapses of biological neurons, and memristor synapses can be used to characterize this electromagnetic effect [
18]. Xu et al. [
19] established a bi-neuron Rulkov network with a memristive synapse and explored the multi-stability and phase synchronization of a neural network; it was also observed that the dynamical behaviors of the neural network depended on the setting of the synaptic coupling strength and initial value, and the neural network could produce rich discharge patterns under the influence of different conditions. Kafraj et al. [
20] discussed the dynamical behaviors of neurons under electromagnetic induction and noise by virtue of the memristive Izhikevich neuron model. Their studies revealed that the neuron model showed complex discharge modes under different parameter conditions. Mondal et al. [
21] established a single Izhikevich neuron under periodic signal stimulation; the chaotic resonance of the model could produce corresponding dynamic responses in different chaotic states by relying on Lyapunov exponent diagrams and bifurcation diagrams. Fang et al. [
22] studied the discharge patterns in the three-dimensional memristive Izhikevich neuron model and constructed an MIZH neuronal network. It was found that it could show rich collective dynamical behaviors, and MIZH neuronal networks may remember and retrieve information efficiently. These investigations revealed that neuronal networks can exhibit several types of complex phenomena, including symmetry-breaking and period-doubling scenarios. It is well known that the nervous system has a large number of neurons with different biological structures and functions [
23]. Therefore, it is necessary to develop more different memristive neuron and neural network models based on the different biological neuronal systems as well as different mathematical models.
As we know, energy is consumed so that neurons can save normal and continuous electric activity during the metabolic process of the neuronal system. Thus, it is interesting to detect the energy transmission and release dependencies on the electric activity mode in these neuronal models. In other words, the supply and consumption of energy is a key step in the metabolism and life activities of organisms. There also exists energy absorption and release in neurons and neural networks; the work of the brain is supported by energy [
24,
25,
26,
27]. In [
28], the authors used Helmholtz’s theorem to derive the Hamilton energy function of the dynamics system, which provided a theoretical basis for the calculation of energy. Lu et al. [
29] estimated the dynamical behaviors and Hamilton energy of the improved memristive Hindmarsh–Rose neuron model under the effect of an external stimulation current and magnetic field. In their numerical simulations, the results showed that the transition of neuronal discharge patterns is closely related to energy changes, and the electrical activity patterns stimulated by periodic signals are more complex than those stimulated by mixed signals. Yang et al. [
30] investigated the effect of periodic current and high–low-frequency electromagnetic radiation on neuronal discharge activity and energy; the model was built by a single Izhikevich neuron model under electromagnetic induction. The results revealed that the injection of stimulus signals can bring energy which is able to make the discharge mode transform to bursting. Besides the above-mentioned studies, there are many studies on the rich dynamic behavior of neuron models. However, the Hamilton energy of coupled neurons has rarely been reported in previous investigations; the study of the interaction between coupled neurons is also a key step in exploring the mysteries of the brain.
Inspired by the above contributions, in the present paper, we will investigate the complex dynamical behavior and Hamilton energy of the coupled neuron model. The five-dimensional neuronal system with memristor synapses is presented in
Section 2. The global exponential synchronization of the neuron system is proved in
Section 3. Moreover, the Hamilton energy functions of the improved Izhikevich neuron system and the error system are derived by Helmholtz’s theorem. In
Section 4, numerical simulations are discussed, which describes the effects of the external stimulus current and parameters on the complex dynamics and energy in the improved system. Finally,
Section 5 concludes the paper.
2. Model Description
To further understand how the brain works, we need to combine experimental studies of human nervous systems with the numerical simulation of large-scale brain models. A neuron model with a simple structure called the Izhikevich model was proposed in 2003, which has rich physiological meanings and unfolds all the known neural behaviors of the cortical neuron [
31]. This model is suitable for large-scale simulations. Depending on four parameters, the model reproduces the spiking and bursting behavior of known types of cortical neurons.
The expression of the Izhikevich model is most consistent with the behavior of cortical neurons. The mathematical expression of this model can be described as follows:
with the post-spike resetting
where
is time, the variables
and
are the membrane potential and recovery variable in the Izhikevich neuron model (1), and
is the applied current. The control parameter
represents the rate of change in recovery variable
and
represents the sensitivity of recovery variable
to the fluctuation of membrane potential
near the threshold.
represent the reset value of membrane potential
and recovery variable
after neuron discharges, which accounts for the activation of
ionic currents and the inactivation of
ionic currents, and it provides negative feedback to
. After the spike reaches its apex (+30 mV), the membrane voltage and the recovery variable are reset. Synaptic currents or injected dc currents are delivered via variable
I. The
part was obtained by fitting the spike initiation dynamics of a cortical neuron so that the membrane potential
has an
scale and the time has a scale. The resting potential in the model is between −70 and −60 mV depending on the value of
. Depending on the history of the membrane potential prior to the spike, the threshold potential can be as low as −55
or as high as −40
. By adjusting the value of these four parameters, the model can simulate spiking or bursting behaviors of neurons.
In the following, considering the existence of electromagnetic induction, a five-dimension neuron system is proposed by coupling two Izhikevich neurons with memristor synapses. The improved Izhikevich model can be written as follows:
with the post-spike resetting
where
and
indicate membrane potentials and recovery variables in the coupled Izhikevich neuron model and
describes the magnetic flux.
and
are the induced currents generated by electromagnetic induction in bidirectionally coupled neurons;
indicates memory conductance controlled by a magnetic flux in the memristor [
32];
denotes the intensity of modulation of
by induced current;
is constant conductance;
is the feedback rate of magnetic flux;
is the effect induced by
on magnetic flux;
means the gain of the electromagnetic induction;
is the magnetic leakage;
means the flux feedback coefficient modulated by the neuron itself.
Analysis of Synchronous Stability
In this section, the global exponential synchronization of the corresponding model is proved by the Lyapunov stability theorem. In order to achieve the global exponential stability for the error system, the drive-response synchronization approach will be employed.
The drive system of the improved Izhikevich model is
Then, the response system is described as follows:
In Equation (4),
is the controllers which are the nonlinear functions designed. Let
,
,
,
,
, so the error system of the drive system and the response system can be calculated. The synchronization manifold is as follows:
Theorem 1. The given initial value of the drive system (3) is , and the initial value of the response system (4) is taken as .
The controllers are designed as follows:
Here,
,
, and the zero solution of the error system (5) satisfies the following inequality:
; here,
is a constant which depends on
, and then, the zero solution of the error system (5) is globally exponentially stable, and the drive system (3) and the response system (4) can be realized as the global exponential synchronization.
Proof. According to the description of the controllers
, Equation (5) can be rewritten in the following format:
We construct the Lyapunov function as follows:
Suppose that
is the smallest eigenvalue of the matrix
,
is the largest eigenvalue of the matrix
, and
; thus,
where
Assume that
is the largest eigenvalue of matrix
, since
; then,
.
Here, , .
Therefore, the zero solution of the error system (5) is globally exponentially stable, and the drive system (3) and the response system (4) achieve global exponential synchronization. □
3. Hamilton Energy of the Izhikevich Model
A neuron under an electromagnetic environment is accompanied by the conversion and migration of energy. And the support of energy in the brain is essential for the discharge activity of a neuron. Here, the Hamilton energy functions are determined, which are associated with the coupled Izhikevich neuron model and the corresponding error system. Based on Helmholtz’s theorem [
33], the Hamilton energy function of the neuron model could be calculated. Helmholtz’s theorem decomposes an arbitrary electromagnetic field
into the superposition of a gradient field
and a vortex field
, and its dynamical equation can be expressed by Equation (12).
where
represents the conservative component and
represents the dissipative component, and
is a skew-symmetric matrix and
is a symmetric matrix.
indicates the gradient matrix of the Hamilton energy function
.
The Hamilton energy function must be satisfied as
Therefore, system (2) is rewritten as
with
By combining Equation (13) with Equation (14), one can obtain
The general solution of Equation (15) is expressed as
The derivative of the Hamilton energy function versus time is denoted by
Substituting the values of system (2), we find that
According to the above derivation, satisfies the existence condition of the Hamilton energy function in Equation (13), and Equation (16) is the Hamilton energy function requested.
Next, we derive the Hamilton energy function
of the error system (6) to study the flow of energy in the maintenance of synchronization. According to the above theorem, we rewrite system (6) as follows:
with
Then, we obtain
the general solution of Equation (21) is expressed as
the derivative of the Hamilton energy function versus time is denoted by
Substituting the values of system (6), we can obtain
Therefore, Equation (22) is the Hamilton energy function of the error system.
4. Numerical Simulations
In this section, to better understand the dynamical behavior of the improved Izhikevich neuron system, numerical simulations are performed under different values of the parameters. Here, we use MATLAB software to plot all the figures. The integration is performed using the fourth-order Runge–Kutta algorithm with a time step of .
The parameters are taken as
,
,
,
,
, and
. In addition, the values of the parameters related to the excitatory neuron and the inhibitory neuron are selected as
,
,
, and
and
,
,
, and
, respectively.
Figure 1 shows the phase portrait of the improved Izhikevich neuron system (2); system (2) exhibits dynamics from stability to chaos. Chaos occurs when the value of membrane potential is larger than 30 mV, which is due to the auxiliary after-spike resetting of the Izhikevich neuron model.
The variations in discharge patterns in the excitatory and inhibitory neurons are explored by setting three groups’ initial values in the improved Izhikevich neuron system (2). The time series of membrane potential with different initial conditions are plotted in
Figure 2. It is shown that the discharge patterns exhibit chattering peaks if the coupled neurons are all excited, and the internal chattering frequency and the time interval of chattering are different within the three groups of initial values in
Figure 2a. And the discharge patterns of the neuron appear as inhibition spiking if the coupled neurons are all inhibited. It is shown that the frequencies of spiking peaks become higher in comparison with the excitatory neuron in
Figure 2b. As a consequence, changing the initial conditions can produce different discharge modes in the different types of neurons in system (2).
As shown in
Figure 3, the time series of the membrane potential are plotted by applying an external stimulus current at amplitude fixation or frequency fixation within the coupled neurons in system (2). At amplitude fixation, the neurons exist in a refractory period which is affected by the frequency of the stimulation current. The coupled excitatory neuronal system presents the chattering pattern, and the chattering frequency is the inequality in
Figure 3a. When the frequency is fixed, the refractory period of neurons disappears gradually. The time series of the membrane potential transform from an inhibition bursting to a tonic bursting action sequence, the time interval between chattering peaks decreases, and the discharge mode becomes more regular, as shown in
Figure 3b. It is found that the frequency of the external stimulation current has little effect on inhibitory neurons, as shown in
Figure 3c. As depicted in
Figure 3d, the discharge mode transforms from spiking to inhibition bursting, and the frequency of the discharge spike increases. Thereby, it is concluded that the discharge mode of the neuron system changes with the amplitude and frequency of the external stimulus current.
In order to further investigate the overall influence of amplitude
and frequency
on the discharge activity and energy of the excitatory neuron in system (2), the bifurcation diagram (ISI) and the largest Lyapunov exponent diagram (LLE) are plotted in
Figure 4. The effect of the frequency and amplitude of the external stimulus current on the dynamics behavior of neurons is consistent in terms of excitability and inhibitability. From
Figure 4a,b, it can be observed that the neuron model shifts from chaos to stability, and its corresponding largest Lyapunov exponent turns from greater than zero to equal to zero. That is, the coupled neuron system transforms from chaotic to periodic as the amplitude increases. And as can be seen from
Figure 4b,d, chaotic and periodic behaviors coexist in the coupled neuron system as the frequency increases. Therefore, the amplitude of external current stimulation can promote periodic discharge activity of the neuron model.
In
Figure 5, the amplitude
and frequency
are chosen as variables to characterize the transformation of the energy in system (2). The evolution of energy is consistent with the dynamical behavior of system (2). It is shown that the Hamilton energy with respect to the dynamical behavior decreases as the amplitude
increases, which corresponds exactly to the amplitude control of the system in
Figure 5a. And the energy shows a smooth change as frequency
increases in
Figure 5b. The two-parameter figures of the Hamilton energy of system (2) are plotted in
Figure 5c,d; the energy varies with a periodic state when the amplitude
and frequency
are varied simultaneously. These results demonstrate that the energy of system (2) depends on the transformation of amplitude
or frequency
, and the energy requirement for the chaotic discharge activities is more than the other one for the periodic discharge activities in the neuronal system.
On the other hand, it is noticed that the modulation intensity of the induced current
and the gain in the electromagnetic induction
play important roles in the expression of Hamilton energy function. The following
Figure 6a–d show the complex dynamical behavior of the coupled neurons with the changing parameters
and
. When the coupled neurons are all excited, the coupled neurons present a mixture of periodic and chaotic states. Nevertheless, when the coupled neurons are all inhibited, system (2) transforms from periodic to chaotic as
increases, and it maintains multi-periodicity all the time as
increases, as shown in
Figure 6c,d. Moreover, the evolutions of the Hamilton energy function related to
and
are depicted in
Figure 7. We find that the Hamilton energy of system (2) increases as
changes, whereas parameter
does not cause the energy to change. Compared to
Figure 6 and
Figure 7, the results imply that the gain in the electromagnetic induction can modulate the energy released by the neuron model.
5. Conclusions
In summary, the dynamical behavior of the coupled Izhikevich neuron model using the memristor synapse has been investigated. Based on the Lyapunov stability theorem, the condition of the global exponential synchronization of the improved Izhikevich neuron model has been derived. Meanwhile, we analyzed the Hamilton energy functions associated with the improved Izhikevich neuron model and its error system according to Helmholtz’s theorem.
Furthermore, the influence of the parameters of the external stimulus currents and the electromagnetic induction has been investigated. Numerical computations indicated that the amplitude of external current stimulation can promote periodic discharge activity of the neuron model. In addition, the amount of energy for the chaotic state could be more than that for the periodic state during the discharge activity in the neuronal system. The dynamical behavior of the neuron model transforms from periodic to chaotic as the modulation intensity of the induced current increases, while the gain in the electromagnetic induction can maintain the periodic discharge pattern and modulate the energy released by the neuron model.
Our approach can be extended to investigate the energy aspects of different oscillatory regimes of nonidentical oscillators. The theorem results will be useful to study the energy aspects of other coupled complex systems. Possible extensions to neural networks can provide finer insight into the energy modulation mechanisms of various biological systems. Moreover, these results may help us better understand the relationship between electrical activity and energy and further understand the energy characteristics of complex systems under different dynamical behaviors.