# Noise and Synchronization Analysis of the Cold-Receptor Neuronal Network Model

Authors: Ying Du; Rubin Wang; Jingyi Qu
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

This paper analyzes the dynamics of the cold receptor neural network model. First, it examines noise effects on neuronal stimulus in the model. FromISI plots, it is shown that there are considerable differences between purely deterministic simulations and noisy ones. The ISI-distance is used to measure the noise effects on spike trains quantitatively. It is found that spike trains observed in neural models can be more strongly affected by noise for different temperatures in some aspects; meanwhile, spike train has greater variability with the noise intensity increasing. The synchronization of neuronal network with different connectivity patterns is also studied. It is shown that chaotic and high period patterns are more difficult to get complete synchronization than the situation in single spike and low period patterns. The neuronal network will exhibit various patterns of firing synchronization by varying some key parameters such as the coupling strength. Different types of firing synchronization are diagnosed by a correlation coefficient and the ISI-distance method. The simulations show that the synchronization status of neurons is related to the network connectivity patterns.


---

## Body

## 1. Introduction

It is expected that the processed neural information can be encoded in the structure of inter-spike-interval (ISI) series; that is, the neural firing activities can be represented by the patterns of neural spike trains [1–4].Since the introduction of the concept of stochastic resonance [5], a positive role of noise has become an important issue in physics and related neural systems. Nonlinear systems can change their dynamics with the addition of noise [1, 6]. Most neural models consist of differential equations that can reproduce observed neuronal behavior. When these models are stimulated by noise, a variety of phenomena including stochastic resonance and coherence resonance will be observed [7–10]. Besides the deterministic state that triggers the membrane potential of neurons, firing spike train is also affected by noise [11]. For example, for subthreshold oscillations near the threshold of spike generation, it is dependent on noise whether a spike is triggered or not. Spike activity with skipping can be fully understood only when noise is considered. Transitions of spike train patterns are also affected by noise [12]. With the addition of noise, it was observed there were a variety of spike train patterns in temperature sensitive skin receptors for information encoding [13, 14]. The corresponding data analyses suggest that noise can extend the dynamical behavior of the system [15, 16]. Similar results can also be observed in experiment [17–20], which indicate that the nervous system may take advantage of the benefits of noises, including internal thermal fluctuations and external noisy inputs.Except the contribution of a single neuron on information transmitting, synchronization of a set of interacting units has been intensively studied in the natural world [21, 22]. The synchronization of neuronal signals has been proposed as one of the mechanisms to encode the signal [23, 24]. It is suggested that theoretical studies of synchronized behaviors in neuronal system play an important role in information processing of olfactory systems. Synchronous behavior of two electrically coupled neurons was studied by Postnova et al. [25]. Asynchronous and other synchronous states were investigated too. The results showed that the coupling strength could have important role on the synchronous states. Synchronization of fast-spiking neurons connected by electrical synapses was investigated by Nomura and his team [26].Many different measures have been introduced to characterize effects of noise and the synchronization transaction of neurons [26–29]. There exist different connectivity patterns including ring-like neuronal network and grid-like neuronal network. The former only considers the coupling of nearest neighbor neurons, while the latter includes all the nearest neighbor connected couplings. Synchronization studies between two coupled neurons and multiple neural networks are studied [30–34].Except the traditional methods, in this paper,ISI distance, a new method recently introduced by Kreuz et al. [35, 36], is used to characterize effects of noise in a Hodgkin-Huxley type model of temperature encoding and synchronization. This is a simple approach that extracts information from the interspike intervals by evaluating the ratio of the instantaneous firing rates, which is complementary to the spike-based approaches. With the addition of noise, the neural model can generate a variety of new different patterns compared to deterministic situation; synchronization status of neural network with different connection types is also investigated, and the simulations show that the synchronization status is much related to the network connectivity patterns.The paper is organized as follows: the model of neural system is introduced in Section2. Noise effects on neural system are characterized in Section 3, including bifurcation plots andISI distance about the noise effects. Section 4 presents the synchronization results of two-coupled neurons and neuronal networks with different connection patterns, showing the different states of synchronization by varying control parameters, the coupling strength. Finally, a brief conclusion is given in Section 5.


## 2. The Neural Model

A simplified Hodgkin-Huxley model [16] is used to simulate the spike train activity of cold receptor response to temperature changes in the skin. This model consists of two sets of ionic conductance, each related to de- and repolarizing currents:
(1)
C
M
d
V
d
t
=
f
(
V
)
≡
-
g
l
(
V
-
V
l
)
-
α
(
I
d
+
I
r
)
-
β
(
I
s
d
+
I
s
r
)
,
(2)
I
i
=
ρ
g
i
a
i
(
V
-
V
i
)
,
(
i
=
d
,
r
,
s
d
,
s
r
)
,
(3)
a
i
∞
=
1
(
1
+
exp
⁡
(
-
s
i
(
V
-
V
0
i
)
)
)
,
(4)
d
a
i
d
t
=
ϕ
(
a
i
∞
-
a
i
)
τ
i
,
(5)
a
d
=
a
d
∞
,
where V is the membrane potential, C
M is the membrane capacitance, and the leakage currents are given by g
l
(
V
-
V
l
). I
d and I
r are the fast depolarizing current and the repolarizing one for spike generation, and two slow currents for subthreshold oscillations are given by I
s
d and I
s
r. α and β are parameters used to alter the spike currents, especially for a
s
r,
(6)
d
a
s
r
d
t
=
ϕ
(
-
η
I
s
d
-
k
a
s
r
)
τ
s
r
.
V
i are the equilibrium potentials, g
i are the conductance, and V
0
i and s
i are the half-activation potentials and the slopes of sigmoidal steady-state activation curves. The constant η is the coupling contrast, while k is the relaxation factor. Temperature dependencies are introduced by the parameters ρ and ϕ, respectively,
(7)
ρ
=
1
.
3
(
T
-
T
0
/
10
)
,
ϕ
=
3
.
0
(
T
-
T
0
)
/
10
,
where T is the temperature and T
0 is the reference temperature.


## 3. Noise Effects on the Neural System

In order to study the stochastic disturbances on the model, noise is implemented to the membrane potential Equation (1) in each time step Δ
t according to the Box-Mueller algorithm [37]:
(8)
V
t
+
Δ
t
=
V
t
+
(
f
(
V
)
C
M
)
Δ
t
+
g
w
,
g
w
=
(
-
4
d
ln
⁡
(
a
)
Δ
t
)
1
/
2
cos
⁡
(
2
π
b
)
,
where a, b are random numbers between 0 and 1, d is the noise intensity, and the noise g
w is determined by its Gaussian statistical features.In the deterministic situation without noise, the neural model described above exhibits different types of spiking patterns when the parameter of temperatureT is varying. At low temperature, spike pattern is period single firing, and then for slightly higher temperature the firing pattern changes into chaos firing. When T is increasing to 1
5
°
C, spike train patterns become bursts, and the number of spikes per burst will be reduced if the temperature is further increased (up to T
=
2
0
°
C), for high temperatures again regular single-spike activity occurs (T
=
3
5
°
C), but with higher spiking frequency. These varying processes of spike train pattern can be seen from theISI plot shown in Figure 1(a), in which successive interspike intervals are plotted versus time as a function of temperature scaling. Starting form high temperatures, the behavior of the spike firing patterns changes from subthreshold oscillations to high-frequency regular single spike. At lower temperatures, there are several abrupt transitions from period one to period four with the number of spikes per burst increases. At temperature around 1
0
°
C, increasing the number of bifurcations occurs followed by chaotic patterns. At still lower temperature, regular single-spike activity with lower frequency appears again.ISI plots with respect to variation of the temperature T in the deterministic situation and the noise situation with d
=
0.5.
(a)
(b)In the noisy simulation, the neural model exhibits temperature dependencies that well follow the main characteristics of the experimental data. When especiallyd is taken as 0.5, some fluctuations of spike timing are visible as shown in Figure 1(b). At T
=
3
5
°
C, instead of regular single-spike generation, two spikes per oscillation cycle appear. At mid temperatures, the abrupt transitions of deterministic simulation are completely smoothed out by noise, which can be seen in the interval plot as shown in Figure 1(b). The qualitative differences occur in the lower and higher temperature situations, in which skipping in the noisy simulations appears and the transitions to chaotic interval sequences in deterministic simulations disappear.ISI distance is an useful method to qualitative the effects of noise on the neuron system [35]. This approach is based on the time intervals between successive spikes as the basic element for measuring the difference of the two spike trains qualitatively. Due to the self-adaptation of this method, it can naturally describe the characterization of neuronal spike trains. TheISI distance between two spike trains x
(
t
) and y
(
t
) is introduced briefly by the following procedure. First, define the interspike interval of the spike train x
(
t
) at time t by
(9)
x
i
s
i
(
t
)
=
min
⁡
(
t
i
x
∣
t
i
x
>
t
)
-
max
⁡
(
t
i
∣
t
i
x
<
t
)
,
x
i
s
i
(
t
)
=
x
i
s
i
(
t
)
=
i
=
1
,
…
,
M
,
t
1
x
<
t
<
t
M
x
,
where t
i
x is the time of the ith spike of the spike train x
(
t
). In the same manner, define the function of ith spike of y
(
t
), and then the normalized ratio I
(
t
) is taken as follows:
(10)
I
(
t
)
=
x
i
s
i
(
t
)
-
y
i
s
i
(
t
)
max
⁡
(
x
i
s
i
(
t
)
,
y
i
s
i
(
t
)
)
.Finally, theI
S
I distance D
I is computed through integrating over time as follows:
(11)
D
I
=
∫
0
t
M
d
t
|
I
(
t
)
|
.Here, this measure is used to qualitatively characterize the difference between deterministic and noisy spike trains of the cold receptor. Figure2 shows the change ofISI distances as a continuous function of temperature. As shown in the picture, at each fixed temperature, theISI distance increases with the noise intensity becoming larger. For example, at T
=
1
0
°
C, when d is taken as 0.005, 0.05, and 0.5, theISI distance is 0.212, 0.457, and 0.603, respectively, which means that the higher noise intensity causes stronger fluctuations of the times in the spike train.Figure 2
The change ofISI distance as a function of temperature T at different noise densities. Solid line: d
=
0.5; dashed line: d
=
0.05; dot-dashed line: d
=
0.5.From the study of noise effect in the neural system, we can obtain that, with the addition of noise, the neural model can generate a variety of different spike train patterns and allow continuous transitions from bursting to subthreshold oscillations. Thus, the system with noise may be more close to actual dynamics underlying the experimental data.


## 4. Synchronization Study in Neural Network System

In the process of information transmitting and code, synchronization of firing plays an important role. Hence in this section, synchronization of two coupled neurons and network is studied.Dynamics of multiple coupled neurons in this paper are controlled as shown by the following differential equations:(12)
k
C
M
d
V
i
d
t
=
f
(
V
)
C
M
d
V
i
d
t
h
≡
-
g
l
(
V
i
-
V
l
)
-
α
(
I
d
+
I
r
)
-
β
(
I
s
d
+
I
s
r
)
I
couple
(
i
)
i
k
+
I
couple
(
i
)
,
I
couple
(
i
)
=
∑
j
g
c
(
V
i
-
V
j
)
,
where g
c is coupling strength and V
i and V
j are voltages between adjacent neuron i and j. For two bidirectional coupled neurons (Figure 3(a)), i
,
j
=
1,2; for the ring-like network with N neurons (Figure 3(b)), I
couple
(
i
) is the sum of the adjacent neurons including previous and next ones:
(13)
I
couple
(
i
)
=
g
c
(
V
i
-
V
i
-
1
)
+
g
c
(
V
i
-
V
i
+
1
)
i
=
1,2
,
…
,
N
.Structure of network connectivity patterns. (a) Two bidirectional coupled neurons. (b) A coupled ring-like neuronal network. (c) A coupled grid-like neuronal network.
(a)
(b)
(c)For the grid-like network (Figure3(c)), the I
couple
(
i
,
j
) of a neuron at position (
i
,
j
) is the sum of the input currents from the nearest neighbor neurons:
(14)
I
couple
(
i
,
j
)
=
∑
g
c
(
V
i
,
j
-
V
i
+
n
,
j
+
m
)
i
,
j
=
1,2
,
…
,
N
.


### 4.1. Synchronization of Two Coupled Neurons

First, synchronization status of two identical coupled neurons is considered. TheISI distance method introduced above is used to estimate the degree of synchrony between two spike trains. As introduced in Section 3, the quantity ofISI ratio I
1,2
(
t
) will be zero if two spike trains are completely alike, or it will approach −1 and 1, respectively, if two spike trains differ greatly. So the more synchronous the coupled neurons are, the less theISI distance D
I is, and D
I is equal to 0 when the complete synchronization state of the coupled neurons is achieved.For two identical coupled neurons as Figure3(a), synchronization of the two-coupled neurons both at regular single spiking is investigated and the numerical result is shown in Figure 4 by varying the coupling strengths. In the figure, the change of theISI distance ((10) and (11)) as a function of the increasing coupling strength is shown. It increases to its maximal value when g
c
=
0.013
ms
/
c
m
2, before that, the I
S
I distance is increasing due to chaotic behavior of the neuron. Then it decreases gradually to 0 in the range of 0.013
<
g
c
<
0.03
ms
/
c
m
2. Finally, it remains as 0 when g
c is larger. TheISI distance here shows that the synchronous degree of two neurons first decreases and then increases to complete synchronization when g
c is large enough. Similar results can be found in other identical coupled neurons with different firing patterns (bursting, chaos firing, etc.).Figure 4
Plot of theISI distance of two coupled neurons as a function of coupling strength g
c.Another method to characterize the synchrony degree between two spike trains is the correlation coefficient(
CC
), and the CC is calculated as follows:
(15)
CC
=
∑
i
-
1
n
|
(
V
i
1
-
〈
V
i
1
〉
)
|
|
(
V
i
2
-
〈
V
i
2
〉
)
|
∑
i
-
1
n
(
V
i
1
-
〈
V
i
1
〉
)
2
(
V
i
2
-
〈
V
i
2
〉
)
2
,
where V
i
1 (or V
i
2) is the sampling of the membrane potential V
1
(
t
) (or V
2
(
t
)). 〈
·
〉 calculates the average over the number of the sampling. So the larger the correlation coefficient is, the more synchronous the two coupled neurons are, and the complete synchronization will be achieved when correlation coefficient is equal to 1.As illustrated in Figure5, the correlation coefficient value (15) increases gradually from 0 to 1 when g
c reaches 0.032 and maintains it thereafter.Figure 5
Plot of correlation coefficient of two coupled neurons as a function of coupling strengthg
c.Except synchronization of two identical neurons, several different synchronized states of distinct neurons are also investigated. Similar results are obtained that all neurons will get synchronized wheng
c is larger enough as shown in Figure 6, in which we put all correlation coefficients with respect to the coupling strength g
c together. But it can be observed from the local enlargement of the figure in the right corner that there are still some differences. In Figure 6, different states correspond to different line styles. It can be seen that CC of neurons (
V
1
(
T
=
3
°
C
) versus V
2
(
T
=
5
°
C
)
) and (
V
1
(
T
=
5
°
C
) versus V
2
(
T
=
7
°
C
)
) gets synchronized rapidly after being coupled with g
c, while these two states are all synchronized in single-spike pattern, which means that this state is more easy to be synchronized. CC of neurons (
V
1
(
T
=
3
2
°
C
) versus V
2
(
T
=
3
5
°
C
)
) is also increased rapidly, but slower than the above two cases because of their high frequencies. The most slowly synchronized neurons are (
V
1
(
T
=
5
°
C
) versus V
2
(
T
=
1
0
°
C
)
), in which two neurons are all in chaos states.Figure 6
Correlation coefficient with respect to the coupling strengthg
c of all the states.


### 4.2. Synchronization of Neuronal Network

First, synchronization effects in neural networks coupled in ring-like manner are investigated. The network is connected by 25 coupled neurons as ring-like manner (Figure3(b)). AnISI-distance method for multiple spike trains is used to give the quantitative results for synchronized spike trains [35]. First, theISI for each spike train is assigned and by computing over all absoluteISI ratios |
I
m
,
n
(
t
)
|, the instantaneous average A
(
t
) is obtained as follows:
(16)
A
(
t
)
=
1
N
(
N
-
1
)
/
2
∑
n
=
1
N
∑
m
=
n
+
1
N
|
I
m
,
n
(
t
)
|
,
D
I
a
=
1
T
∫
t
=
0
T
A
(
t
)
d
t
.
Average over time, D
I
a, is achieved. It is obvious that the more synchronous the multiple spike trains are, the smaller theISI distance D
I
a is; when complete synchronization state of the neuronal network is achieved, D
I
a is equal to 0. Then the effect of two parameters on synchronization of ring-like network is investigated. The contour graph of theISI distance in the (
g
c
,
T
) parameter plane is shown in Figure 7. The plot shows that when coupling strength is bigger than 0.02
ms
/
c
m
2, no matter what the value of T is, the neurons are in a complete synchronous state. Also the region thatISI distance almost equals 0 is much larger than that of only changing one parameter g
c shown as Figure 4. It indicates that the realization of synchronization in the ring-like neuron network is easier than that of two coupled neurons with changing only one parameter. Another obvious conclusion is that, no matter what the temperature T is, when increasing the coupling strength, the neurons in the ring-like network can finally reach synchronization.Figure 7
The contour plot of theISI distance in the (
g
c
,
T
) plane for ring-like network.Next, to understand the relationship between synchronization and the network connection patterns, network connection type is changed to grid network. The number of neurons in the grid-like network is still 25. The coupling current of a neuron is the sum of the currents for the nearest neighbor neurons (M
=
N
=
3). The corresponding contour graph of theISI distance in the (
g
c
,
T
) parameter plane is shown in Figure 8. As illustrated in the figure, the map differs from that of the ring-like neuronal network. It can be observed that the region thatISI distance almost equals 0 is larger than that of the ring-like neuronal network; when coupling strength is bigger than 0.012
ms
/
c
m
2, all neurons are in a complete synchronous state. It indicates that synchronization state is related to the network connection type; that is, the synchronization in coupled grid-like neuronal network is easier than ring-like connected neurons. Also we can see from the simulation that, no matter what the T is (the uncoupled neurons are in all kinds of firing patterns including regular single spike, chaotic, or bursting status), when the coupling strength is bigger enough, the grid-like neuronal network can eventually reach the complete synchronization, which is similar to ring-like network.Figure 8
The contour plot of the averageISI distance in the (
g
c
,
T
) plane for grid network.Compared to Figure1, it can be observed from Figures 7 and 8 that g
c of chaotic and high period pattern are larger than the ones in spike and low period patterns, which means that chaotic and high period patterns are more difficult to get complete synchronization than the situation in single spike and low period patterns. It can be also seen from these two figures that the tendency of I
S
I with temperature T
(
g
c
<
0.02 in Figure 7 and g
c
<
0.015 in Figure 8) is consistent with the situation of single neuron, that is, increasing first and then decreasing. This means that it is possible to establish the relation between network dynamics and single-neuron dynamics.


## 4.1. Synchronization of Two Coupled Neurons

First, synchronization status of two identical coupled neurons is considered. TheISI distance method introduced above is used to estimate the degree of synchrony between two spike trains. As introduced in Section 3, the quantity ofISI ratio I
1,2
(
t
) will be zero if two spike trains are completely alike, or it will approach −1 and 1, respectively, if two spike trains differ greatly. So the more synchronous the coupled neurons are, the less theISI distance D
I is, and D
I is equal to 0 when the complete synchronization state of the coupled neurons is achieved.For two identical coupled neurons as Figure3(a), synchronization of the two-coupled neurons both at regular single spiking is investigated and the numerical result is shown in Figure 4 by varying the coupling strengths. In the figure, the change of theISI distance ((10) and (11)) as a function of the increasing coupling strength is shown. It increases to its maximal value when g
c
=
0.013
ms
/
c
m
2, before that, the I
S
I distance is increasing due to chaotic behavior of the neuron. Then it decreases gradually to 0 in the range of 0.013
<
g
c
<
0.03
ms
/
c
m
2. Finally, it remains as 0 when g
c is larger. TheISI distance here shows that the synchronous degree of two neurons first decreases and then increases to complete synchronization when g
c is large enough. Similar results can be found in other identical coupled neurons with different firing patterns (bursting, chaos firing, etc.).Figure 4
Plot of theISI distance of two coupled neurons as a function of coupling strength g
c.Another method to characterize the synchrony degree between two spike trains is the correlation coefficient(
CC
), and the CC is calculated as follows:
(15)
CC
=
∑
i
-
1
n
|
(
V
i
1
-
〈
V
i
1
〉
)
|
|
(
V
i
2
-
〈
V
i
2
〉
)
|
∑
i
-
1
n
(
V
i
1
-
〈
V
i
1
〉
)
2
(
V
i
2
-
〈
V
i
2
〉
)
2
,
where V
i
1 (or V
i
2) is the sampling of the membrane potential V
1
(
t
) (or V
2
(
t
)). 〈
·
〉 calculates the average over the number of the sampling. So the larger the correlation coefficient is, the more synchronous the two coupled neurons are, and the complete synchronization will be achieved when correlation coefficient is equal to 1.As illustrated in Figure5, the correlation coefficient value (15) increases gradually from 0 to 1 when g
c reaches 0.032 and maintains it thereafter.Figure 5
Plot of correlation coefficient of two coupled neurons as a function of coupling strengthg
c.Except synchronization of two identical neurons, several different synchronized states of distinct neurons are also investigated. Similar results are obtained that all neurons will get synchronized wheng
c is larger enough as shown in Figure 6, in which we put all correlation coefficients with respect to the coupling strength g
c together. But it can be observed from the local enlargement of the figure in the right corner that there are still some differences. In Figure 6, different states correspond to different line styles. It can be seen that CC of neurons (
V
1
(
T
=
3
°
C
) versus V
2
(
T
=
5
°
C
)
) and (
V
1
(
T
=
5
°
C
) versus V
2
(
T
=
7
°
C
)
) gets synchronized rapidly after being coupled with g
c, while these two states are all synchronized in single-spike pattern, which means that this state is more easy to be synchronized. CC of neurons (
V
1
(
T
=
3
2
°
C
) versus V
2
(
T
=
3
5
°
C
)
) is also increased rapidly, but slower than the above two cases because of their high frequencies. The most slowly synchronized neurons are (
V
1
(
T
=
5
°
C
) versus V
2
(
T
=
1
0
°
C
)
), in which two neurons are all in chaos states.Figure 6
Correlation coefficient with respect to the coupling strengthg
c of all the states.


## 4.2. Synchronization of Neuronal Network

First, synchronization effects in neural networks coupled in ring-like manner are investigated. The network is connected by 25 coupled neurons as ring-like manner (Figure3(b)). AnISI-distance method for multiple spike trains is used to give the quantitative results for synchronized spike trains [35]. First, theISI for each spike train is assigned and by computing over all absoluteISI ratios |
I
m
,
n
(
t
)
|, the instantaneous average A
(
t
) is obtained as follows:
(16)
A
(
t
)
=
1
N
(
N
-
1
)
/
2
∑
n
=
1
N
∑
m
=
n
+
1
N
|
I
m
,
n
(
t
)
|
,
D
I
a
=
1
T
∫
t
=
0
T
A
(
t
)
d
t
.
Average over time, D
I
a, is achieved. It is obvious that the more synchronous the multiple spike trains are, the smaller theISI distance D
I
a is; when complete synchronization state of the neuronal network is achieved, D
I
a is equal to 0. Then the effect of two parameters on synchronization of ring-like network is investigated. The contour graph of theISI distance in the (
g
c
,
T
) parameter plane is shown in Figure 7. The plot shows that when coupling strength is bigger than 0.02
ms
/
c
m
2, no matter what the value of T is, the neurons are in a complete synchronous state. Also the region thatISI distance almost equals 0 is much larger than that of only changing one parameter g
c shown as Figure 4. It indicates that the realization of synchronization in the ring-like neuron network is easier than that of two coupled neurons with changing only one parameter. Another obvious conclusion is that, no matter what the temperature T is, when increasing the coupling strength, the neurons in the ring-like network can finally reach synchronization.Figure 7
The contour plot of theISI distance in the (
g
c
,
T
) plane for ring-like network.Next, to understand the relationship between synchronization and the network connection patterns, network connection type is changed to grid network. The number of neurons in the grid-like network is still 25. The coupling current of a neuron is the sum of the currents for the nearest neighbor neurons (M
=
N
=
3). The corresponding contour graph of theISI distance in the (
g
c
,
T
) parameter plane is shown in Figure 8. As illustrated in the figure, the map differs from that of the ring-like neuronal network. It can be observed that the region thatISI distance almost equals 0 is larger than that of the ring-like neuronal network; when coupling strength is bigger than 0.012
ms
/
c
m
2, all neurons are in a complete synchronous state. It indicates that synchronization state is related to the network connection type; that is, the synchronization in coupled grid-like neuronal network is easier than ring-like connected neurons. Also we can see from the simulation that, no matter what the T is (the uncoupled neurons are in all kinds of firing patterns including regular single spike, chaotic, or bursting status), when the coupling strength is bigger enough, the grid-like neuronal network can eventually reach the complete synchronization, which is similar to ring-like network.Figure 8
The contour plot of the averageISI distance in the (
g
c
,
T
) plane for grid network.Compared to Figure1, it can be observed from Figures 7 and 8 that g
c of chaotic and high period pattern are larger than the ones in spike and low period patterns, which means that chaotic and high period patterns are more difficult to get complete synchronization than the situation in single spike and low period patterns. It can be also seen from these two figures that the tendency of I
S
I with temperature T
(
g
c
<
0.02 in Figure 7 and g
c
<
0.015 in Figure 8) is consistent with the situation of single neuron, that is, increasing first and then decreasing. This means that it is possible to establish the relation between network dynamics and single-neuron dynamics.


## 5. Conclusion

First, the noise effects on the neuronal system are studied. With the addition of noise, the neural model can generate a variety of new different patterns compared to deterministic situation. For a qualitative computing of noise effect, theI
S
I distance is introduced to characterize how noise affects the variability of the neural spike train. It is shown that the I
S
I distance will increase with the noise intensity becoming larger, which means that the higher noise intensity can cause stronger fluctuations in the spike train. Meanwhile, noise has less effect on spike trains for higher temperature within the bursting regime. These stochastic phenomena mean that signal encoding can be different with the change of noise level.Then synchronization status of neural network with different connection types is investigated. The simulations show that the synchronization status is much related to the network connectivity patterns. By tuning both temperature and coupling strength, there exists a larger synchronous region for grid-like neuronal network than that of ring-like network. It is also shown that chaotic and high period patterns are more difficult to get complete synchronization than the situation in single spike and low period pattern.The method and results in this study may be instructive for the study of neural encoding of spike trains under noisy and network environment. Network study results show that it may be possible to establish the relation between network dynamics and single-neuron dynamics.



---

*Source: 173894-2014-06-15.xml*