# Comparison and Regulation of Neuronal Synchronization for Various STDP Rules

Authors: Yanhua Ruan; Gang Zhao
Journal: Neural Plasticity
 (2009)

---

## Abstract

We discuss effects of various experimentally supported STDP learning rules on frequency synchronization of two unidirectional coupled neurons systematically. First, we show that synchronization windows for all STDP rules cannot be enhanced compared to constant connection under the same model. Then, we explore the influence of learning parameters on synchronization window and find optimal parameters that lead to the widest window. Our findings indicate that synchronization strongly depends on the specific shape and the parameters of the STDP update rules. Thus, we give some explanations by analyzing the synchronization mechanisms for various STDP rules finally.


---

## Body

## 1. Introduction

Synchronous activity is a basic characteristic in the brain. It exists in many regions of the brain, such as CA1 of the hippocampus [1], visual cortex [2]  and cortical areas correlating with conscious perception [3]. It is known that synchronization is very important for information processing, such as predicting sensory input [4], and information codes [5]. Moreover, synchronous activity plays a crucial role in epileptic activity [6, 7], modulation of neurons about attention [8], memory and learning [9, 10], and cognitive functions [11].Since the discovery of long term potentiation (LTP) and LTD (long term depression) [12–14], it has been debated how synaptic modifications are correlated to neuron activities. Spike timing-dependent plasticity (STDP) is a form of synaptic modification discovered relatively recently, which depends on the relative timing of pre- and post-synaptic action potentials at a millisecond time scale [15, 16]. Many experiments have proved the existence of STDP, such as in neocortical slices [17], hippocampus slice [18], hippocampal cell cultures [19], and tadpole rectum in vivo [20]. In addition, STDP provides powerful mechanisms  for models of temporal pattern recognition [21], temporal sequence learning [22, 23], a continuous-time associative memory [24], coincidence detection [16, 25], navigation [26, 27] and direction selectivity [28].The interaction among neurons relies much on synaptic modification in which STDP is the only one that greatly expands the capability of Hebbian learning to address temporally sensitive computational tasks. STDP in synchronization has attracted wide interests. For example, the result of learning-induced synchronization of a neural network at various developing stages using STDP rule is consistent with recent experimental observations [29]. Furthermore, the comparison of synchronization between discontinuous anti-STDP(dc-aSTDP, see Section 2) and constant connection has been investigated [30]. Following it, the continuous STDP(c-STDP, see Section 2) has also been studied [31] by the same authors. They suggest that a functional role of STDP might be enhancing synchronization. Motivated by their work, we systemically discuss the roles of four types of STDP rules(c-STDP, dc-STDP, dc-aSTDP and in-STDP, see Section 2)  in frequency synchronization in the present paper, employing the same model [31] with only values of some parameters different, such as Aplus, tsyn, Vslope, gmax (see Section 3).We find, however, not all STDP rules facilitate synchronization. It encourages us to trace the reason. We then consider if the learning curves, which characterize the STDP rule, have certain effects on synchronization. Results indicate that synchronization strongly depends on the specific shape and the parameters of the STDP rule. However, the optimal synchronization ranges for dc-STDP and in-STDP, got from regulating learning parameters, are not wider than those for the corresponding strongest constant connection respectively. As a result, when we seek the reason, we discover that the synchronization mechanisms of above four STDP rules can be classified into two categories: (i) c-STDP and dc-aSTDP rules; (ii) dc-STDP and in-STDP rules. The synchronization mechanisms of the two categories are different. For c-STDP and dc-aSTDP rules, two neurons’ synchronization either relies on the balancing out potentiation and depression during one cycle consistent with the perspective of Nowotny et al. or relies on  the maximal synaptic conductance. However, for dc-STDP and in-STDP rules, the synchronization windows are completely provided by the respective maximal synaptic conductance. As regards this finding, we offer an intuitive explanation finally.


## 2. Models and Method

We consider two HH neurons with unidirectional activity-dependent excitatory or inhibitory synaptic coupling. Although such a configuration is too simple to find applications in brain information processing, it serves as a staring point for many model researches. The neurons are modeled with standard Na, K, and “leak” currents [32],(1)CdVi(t)dt=-gNa·mi(t)3·hi(t)·(Vi(t)-ENa)-gK·ni(t)4·(Vi(t)-EK)-gL·(Vi(t)-EL)-Isyn(t)+Istim,
where i=1,2.Each of the activation and inactivation variablesyi(t)={ni(t),mi(t),hi(t)}, i=1,2 satisfies first-order kinetics,(2)dyi(t)dt=αy[Vi(t)][1-yi(t)]-βy[Vi(t)]yi(t),i=1,2.The parameters in these equations are given in [31],(3)αn=0.032(-50-V)exp((-50-V)/5)-1,βn=0.5exp((-55-V)40),αm=0.32(-52-V)exp((-50-V)/4)-1,βm=0.28(25+V)exp((25+V)/5)-1,αh=0.128exp(-48-V18),βh=4exp((-25-V)/5)+1.Istim is a constant input current forcing each neuron to spike with a constant, Istim-dependent period, labeled as T1 and T2. The postsynaptic neuron would show another firing period T21, when it is driven by the synaptic current, which is dependent on the postsynaptic potential V2(t), the reversal potential Vrev, the activation variable S(t) and its  maximal conductance g(t),(4)Isyn(t)=g(t)S(t)(V2(t)-Vrev),
where(5)dS(t)dt=S∞(V1(t))-S(t)tsyn·1-(S∞(V1(t))),S∞(V)={tanh(V-VthVslope),forV>Vth,0,otherwise.The time-dependent synaptic coupling strengthg(t)nS is(6)g(t)=gmax2(tanh(graw-gmidgslope)+1).Thereforeg(t) always have values between 0 nS and gmax. The bound imposed on g(t) is artificially set to avoid unrealistically high synaptic conductance and negative conductance. In order to obtain biologically plausible synaptic conductance, several methods have been employed to limit the synaptic strength in literature, such as a negative total integral [33], artificial bounds [34], and self-limitation [31]. Unless otherwise stated, we employ the self-limitation method which is characterized by a function “tanh” in our simulation.graw is modified by STDP rules that are introduced in the next paragraph. The initial value of graw is 20 nS. The parameters of the model are(7)C=30μF,gL=1μS,EL=-64mV,gNa=360μS,ENa=50mV,gK=70μS,EK=-95mV,Vth=-20mV,tsyn=25ms,Vslope=15mV,gmax=25nS,gmid=12gmax,gslope=gmid,Vrev=20mV.The time-dependent synaptic coupling strengthg(t) is determined by the spike-timing of pre- and postsynaptic spikes. We consider four types of activity-dependent couplings that have been found in experiments: (1) an excitatory synapse with continuous STDP (c-STDP). There are two forms of c-STDP from two different experiments. One (Figure 1(a)) is from the recording of the neocortex-layer 5 Xenopus tectum hippocampus [31, 35], and the other is from the neocortex-layer 4 spiny stellates [36]. The latter form will not be considered here, because it introduces persistent decrease to synaptic strength that would result in none synchronization if two neurons have different inherent periods. (2) an excitatory synapse with discontinuous STDP (dc-STDP, Figure 1(b)) [17, 36]; (3) an excitatory synapse with discontinuous anti-STDP [36, 37] (dc-aSTDP, Figure 1(c)); (4) an inhibitory synapse with STDP (in-STDP, Figure 1(d)) [18, 19, 36].Different types of STDP curves are presented. (a) (c-STDP ) (b) discontinuous STDP(dc-STDP); (c) discontinuous anti-STDP(dc-aSTDP); (d) inhibitory STDP(in-STDP).
(a)(b)(c)(d)Δgraw is a function of Δt=tpostspike-tprespike, time difference between the times of postsynaptic and presynaptic spikes. The learning rules corresponding to Figures 1(a), 1(b), 1(c), 1(d)) are provided as follows
(8)c-STDP:Δgraw={AplusΔt-τ0tplus·e-Δt/tplus,Δt>τ0,AsubΔt-τ0tsub·eΔt/tsub,Δt≤τ0,dc-STDP:Δgraw={Aplus·e-Δt/tplus,Δt>0,-Asub·eΔt/tsub,Δt≤0,dc-aSTDP:Δgraw={-Aplus·e-Δt/tplus,Δt>0,Asub·eΔt/tsub,Δt≤0,in-STDP:Δgraw={Aplus·(e-Δt/tplus-0.5),Δt>0,Asub·(eΔt/tsub-0.5),Δt≤0,Aplus=Asub.Synchronization of pre- and post-synaptic neurons occurs when|T1-〈T21〉| is limited in an acceptable range. We set the criteria of synchronization as |T1-〈T21〉|<1.5 milliseconds. Although there is some arbitrariness in setting the criteria of synchronization, there is no qualitative change in our results if the criteria change in two folds. Each simulation runs 20 seconds, average is taken in the final 4000 milliseconds. We have observed that simulations from different initial values of V2, S could result in different outputs, that is, the post-synaptic neuron sometimes synchronizes with the pre-synaptic one, sometimes keep its original period, or sometimes fires with an oscillating period (see Figure 2(c)). We therefore carry out 40 times of simulations, from randomly selected initial values, for every T2. The standard deviation of |T1-〈T21〉|, indicating how precisely the neurons are synchronized, represents the quality of synchronization. Range of T2, in which post-synaptic neuron is successfully entrained by the pre-synaptic neuron, that is, |T1-〈T21〉|<1.5 in all 40 simulations, is defined as the synchronization window.Synchronization results for dc-STDP rule carrying out 40 times. Figure2(a) shows the number of synchronization times (ST) of fixed T2 varying from 150 milliseconds to 320 milliseconds. Synchronization window is from 194 to 221. The probabilistic synchronization window is from 222 to 289. Figure 2(b) presents the quality of synchronization against the ratios of uncoupled periods. In Figure 2(c), we fix T2=265 milliseconds which falls into the probabilistic synchronization window. There are three states of the coupled period of post-synaptic neuron T21 when we carry out 40 stimulation times: (i) keeping the initial period (squares); (ii) oscillating (circles); (iii) synchronizing with the pre-synaptic neuron (dots).
(a)(b)(c)


## 3. Results

### 3.1. Comparison of Synchronization Windows of Different Types of STDP

We investigate the width of synchronization window of various STDP curves, with the same set of parameters. The period of the pre-synaptic neuron is chosen to be 171 milliseconds, which falls into the range of theta waves. Several reasons make us choose such a long period. First, it has fairly wide synchronization windows which allow comparisons in a relative precise manner and can provide clearer information about synchronization windows of various STDP learning rules. Second, the slow theta waves always involve many neurons that fire synchronously [38, 39]. Also, theta waves have many interesting implications. For example, theta waves are normally absent in healthy awake adults, but appear during the state of meditation [40].  During emotional arousal and various types of rhythmic activities during sleep, neurons in the amygdala produce theta activity [41, 42]. And it is known that coherent theta activity (4–8 Hz) in amygdala-hippocampal circuits is deeply involved in fear memory [43].With fixed period of the pre-synaptic neuronT1, we evaluate the coupled period of postsynaptic neuron T21 when it is driven by the pre-synaptic neuron. The values of learning parameters used in c-STDP, dc-STDP, dc-aSTDP are Aplus=9 nS, Asub=6 nS, tplus=100 milliseconds, tsub=200 milliseconds, additionallyτ0=30 milliseconds in c-STDP, and in in-STDP are  Aplus=Asub=8 nS, tplus=100 milliseconds, tsub=200 milliseconds. Our model and most values of parameters are from the model of Nowotny for c-STDP, except Aplus, tsyn, Vslope, gmax are different [31]. Especially, T1 is fixed at 171 milliseconds in our simulations while T2 is set to constant value 300 milliseconds of Nowotny’s work.The window of synchronization (upper panel) and quality (middle panel) of dc-STDP are presented in Figure2 as an example. We scan T2 from 150 milliseconds to 320 milliseconds. The upper panel shows the number of synchronization times in 40 simulations, for each T2. It is clear that, in certain range of T2, simulations from different initial values may have different results. Only when T2 falls into the segment from 194 to 221, the post-synaptic neuron can synchronize with the pre-synaptic neuron from any initial value. It is easily found that there are some T2 corresponding to the number of synchronization times between 1 and 39. In this situation, we present the three possible states of post-synaptic neuron’s firing in Figure 2 (lower panel)—keeping the initial period (squares), oscillating (circles), and synchronizing with the pre-synaptic neuron (dots). Obviously, these states are independent of the synchronization criteria we set.For the purpose of discussing the function role of STDP rules in synchronization, the synchronization windows of various type of STDP are plotted in Figure3, in which the case of constant synaptic conductance is also included as a comparative tool. The same parameter values used in simulations ensure a fair comparison. In our simulation studies, the synapse strength is between 0 nS and 25 nS. We choose the maximal synaptic strength and the middle synaptic strength of STDP synapse as the synaptic strength of constant synapse in this study. Interestingly, because the synchronization windows for STDP rules are narrower than synchronization window for constant synapse g=25 nS in Figure 3, these results, opposite to previous reports, indicate that all STDP rules do not enhance synchronization comparing with the constant synapse under the chosen parameters.Figure 3
The synchronization window of different types of learning curves. The top five lines (points) are excitatory synapses and the bottom three lines (circles) are inhibitory synapses. See legends in the figure.Several other points are worthy of detailed describing. First, increasing the excitatory constant synaptic connections from 12.5 nS to 25 nS, leads to a wider synchronization window. However, the wider window could not totally contain the smaller one. It extends toward largerT2, but loses a portion of smaller T2. The case of increasing inhibitory constant synaptic connections is alike. Second, the widest ranges of synchronization window are achieved by the excitatory constant connection g=25 nS and c-STDP rules. However, the lower boundary of the synchronization window of c-STDP is much nearer to T1 than that of excitatory constant synapse. Thirdly, although c-STDP and dc-STDP are fitted from the same set of experimental data, c-STDP has a much wider synchronization window than dc-STDP.We conclude that all STDPs do not give rise to enhanced synchronization and the window of dc-STDP is surprisingly narrow under the chosen parameters in Figure3. Therefore, the questions about what bring about these results inspire us to study further. At the same time, Figure 2 shows a large part of probabilistic synchronization, whose range is a subset of the synchronization window of constant synaptic connection with g=25 nS. We are interested in if the probabilistic synchronization could be enhanced into absolute synchronization by modulation of learning curves. These are the theme of the next section.


### 3.2. The Effect of Learning Parameters on Synchronization

In order to establish the functional role of STDP clearly, we consider if the learning parameters for each STDP rule have important effect on synchronization. In addition, synchrony-asynchrony transition plays important role in the brain. An increase in the degree of synchrony of a uniform input can cause transitions between memorized activity patterns in the order presented during learning. However, if synchronous input is at a low level, transitions cannot occur [44]. The synchrony-asynchrony transition have also been implemented in controlling winner-take-all competition [45], the next recalled time of associative memory [46] and the fine structure of cell assemblies [47].In this section, we will discuss the flexibility of the synchronization window, by exploring regulation of width of the synchronization window,whose boundary indicates the synchrony-asynchrony transition. We take the modulation of learning parameters as the method to regulate the synchronization window. There are four parameters that determine a learning rule Aplus, Asub, tplus, tsub. With three of them fixed and only one parameter changing, we could explore its influence on the width of the synchronization window. For example, Aplus increases from an adequately small value 1 nS to 20 nS with a step of 1 nS. We have also scanned values that are beyond 20 nS, but find that the effect of increasing Aplus is saturated around 20 nS. Further increasing Aplus brings no more effect. Other parameters are fixed as: Asub = 6 nS, tplus = 100 milliseconds, tsub = 200 milliseconds, T2=233 milliseconds. In these conditions, we present the effect of Aplus on the location of the synchronization window with the dc-STDP rule. We carry out simulations 40 times, each from different initial values.Figure4(a) shows the value of ARP (average change of relative period ratio) = 〈T2-〈T21〉〉/(T2-T1) for different Aplus. Some points have value 0 or 1, which means post-synaptic neuron keeping initial period or achieving synchronization with the pre-synaptic one, respectively. Some points have values other than 0 or 1. Figure 4(b) gives an explanation that these points correspond to probabilistic synchronizations with fixed T2=233.It shows that the regulation ofAplus to synchronization for the dc-STDP carrying out 40 times of stimulations. Asub = 6 nS, tplus = 100 milliseconds, tsub = 200 milliseconds, T2 = 233 milliseconds. Top: the value of ARP = 〈T2-〈T21〉〉/(T2-T1) is 0 or 1 which, respectively, means post-synaptic neuron keeping initial period or achieving synchronization. Bottom: the number of synchronization times (ST) in 40 stimulations againstAplus.
(a)(b)We find that the absolute synchronization range is from 10 nS to 20 nS in Figure4(b). According to the definition of synchronization window of T2, we can similarly define 10 nS to 20 nS as the synchronization window of Aplus, with dc-STDP rule and other fixed parameters. The boundary of this synchronization window indicates where synchrony-asynchrony transition happens when changing Aplus.From the results of Figure4, the reason why synchronization window of constant synapse is wider than that of STDP rules (Figure 3) may be explained by learning parameters. To figure out a global picture of the effect of Aplus on synchronization for c-STDP rule, we then determine the synchronization window of Aplus with different T2. Aplus increases from 1 nS to 20 nS with a step of 1 nS while other three parameters keep initial values: Asub=6 nS, tplus=100 milliseconds, tsub = 200 milliseconds. We choose some typical values of T2 to character the global picture. The points in Figure 5(a) show the synchronization range of Aplus with T1 divided by the chosen values of T2. The lower boundaries, as well as those upper boundaries that are other than 20 nS, indicate the position of synchrony-asynchrony transitions. For example, when T2 is 177 milliseconds equivalent to T1/T2 = 0.966, Aplus outside of the points range from 8 nS to 11 nS cannot lead to synchronization between the two neurons. In addition, for those points marked on the horizontal axis, synchronization could not be established no matter what values Aplus take.The range of learning parameters leading to synchronization for all 40 stimulations for c-STDP is presented. The points are the range ofAplus and tplus which can make neuron synchronization. The circles are the range of Asub and tsub. We choose some values of T2 to investigate the effect of Aplus, Asub, tplus and tsub on global synchronization. Aplus and Asub vary from 1 nS to 20 nS. Let tplus and tsub vary from 10 milliseconds to 400 milliseconds. There are some values of T1/T2 marked on the horizontal axis that cannot be entrained to achieve synchronization in Figures 5(a) and 5(b).
(a)(b)According to this global picture, the intersection of those synchronization ranges ofAplus, which is from 10 nS to 11 nS, identifies the range of Aplus that would lead to the widest synchronization window which is from 177 to 289. This optimal synchronization window for c-STDP is wider than constant synapse g = 25 nS. Thus, the learning parameters strongly influence the role of STDP on synchronization.Taking the same method asAplus, we study the effects of other three parameters on synchronization. Figure 5 presents the situation of Asub varying from 1 nS to 20 nS with the step 1 nS (circles, Figure 5(a)), tplus (points, Figure 5(b)) and tsub (circles, Figure 5(b)) both varying from 10 milliseconds to 400 milliseconds with the step 10 milliseconds. When changing one parameter to explore how the range of synchrony evolves with T2, other three parameters keep their initial values as in Figure 3. Similarly, the global modulation picture of four learning parameters for other learning rules can be got. We only give the situation for c-STDP rule in Figure 5.The optimal synchronization windows for various learning rules are presented in Figure6 comparing with previous synchronization windows obtained in Figure 3. The parameters used for optimal synchronization windows are presented in Table 1. The parameters are derived according to regulating one parameter while other three parameters keep initial values. Under the optimal parameters, the synchronization windows for STDP rules are not narrower than constant synapse g = 25 nS. Consequently, the reported important role of STDP in synchronization should be dependent on learning parameters.Table 1
The parameters for optimal synchronization windows.
STDP ruleAplus (nS)Asub (nS)tplus (ms)tsub (ms)c-STDP96100270dc-STDP9610050dc-aSTDP96100350in-STDP88100200Figure 6
We present the optimal synchronization window for various learning rules compared with Figure3. The top and the bottom lines of this figure are the synchronization windows of excitatory and inhibitory constant synapse, respectively. In addition, there are four pair lines in the middle panel for STDP rules. Each pair includes optimal synchronization window (top line) and previous synchronization window in Figure 3 (bottom line). The left range is optimal synchronization range and the range in round bracket is previous synchronization range.For the first three STDP rules, the optimal synchronization windows are got by regulatingtsub while other three fixed learning parameters keep initial values. For the last STDP rule, the initial values are the optimal parameters.


### 3.3. The Synchronization Mechanism

It is important to understand the properties which neural synchronization depends on. We take into account this problem from two aspects.On the one hand, synchronization correlates with the chosen STDP rules. For excitatory synapse, we find that the optimal synchronization windows for c-STDP rule and dc-aSTDP rule are almost equal, and are wider than constant synapses. However, comparing with these two rules, the optimal synchronization window for dc-STDP rule is much narrower. For inhibitory STDP, the synchronization window is the same as that for constant synapse under connection strengthg=25 nS.Accordingly, what makes the optimal synchronization windows for various STDP rules different deserves an explanation. The stationary synaptic conductance is a necessary condition for stationary synchronized state [31]. We find that the mechanisms of synchronization caused by c-STDP (or dc-aSTDP) and dc-STDP (or in-STDP) are different in our model.Figures7(a) and 7(b) show the average Δt=tpostspike-tprespike, and synaptic strength after an episode of coupling time for c-STDP rule. Parameters in Figure 7 are same as those in Figure 3. There are two types of behavior for Δt when synchronization occurs (Figure 7(a)). In a section of constant Δt, the synaptic strength does not achieve the maximal value. Apparently, in this situation, postsynaptic neuron achieves synchronization with the pre-synaptic neuron depending on the balance between potentiation and depression of synaptic conductance. In the rest part of synchronization window, the synaptic strength achieves the maximal value (Figure 7(b)). It indicates that for larger T2, postsynaptic neuron achieves synchronization depending on the effect of maximal synaptic conductance. For dc-aSTDP rule, the synchronization mechanisms are similar to c-STDP rule. When the post-synaptic neuron synchronizes with the pre-synaptic neuron under small T2, the change of synaptic potentiation and depression cancel each other. However, for the small portion of synchronization window at the right side, synaptic conductance gets the maximum at the stationary synchronized state.Four pictures are obtained from the same values of parameters and model of Figure3. The front four pictures take the function of tanh to limit synaptic strength for c-STDP and dc-STDP rules. Figures 7(a) and 7(c) show the average spike time interval of postneuron’s and preneuron’s spike time for c-STDP and dc-STDP over some time after a period time of coupling, respectively. Figures 7(b) and 7(d) present the average synaptic strength for c-STDP and dc-STDP rule, respectively. Each subplot has two dash lines what indicate the boundary of synchronization window.
(a)(b)(c)(d)For dc-STDP and in-STDP rules, the synchronization mechanisms may be different with the above two STDP rules. Because the potentiation and depression of synaptic conductance cancel each other,Δt must be a fixed value for the selected STDP rule. However, for the dc-STDP rule (Figures 7(c) and 7(d)), Δt keeps varying which means that the potentiation and depression of synaptic conductance do not cancel each other at the synchronization state. But, it is easily found that the synaptic conductance is at the stationary maximum for dc-STDP rule. Thus, postsynaptic neuron achieves synchronization completely depending on the effect of the maximal synaptic conductance for dc-STDP rule. The state of in-STDP rule is similar to dc-STDP rule.As a result, neural synchronization mechanism can be different for various STDP rules. For the fewT2 at the right side of the synchronization window, the synaptic conductance achieves the maximum with c-STDP and dc-aSTDP rules. This result is obvious. Because the frequency mismatch is larger, the synapse needs to be stronger to entrain the post-synaptic neuron. But for most part of synchronization window at the left, c-STDP and dc-aSTDP rules rely on the balance of potentiation and depression. Instead of balancing out potentiation and depression during one cycle, dc-STDP and in-STDP rules depend on synaptic strength achieving its maximum. It is important to understand why c-STDP and dc-aSTDP rely on the balance of potentiation and depression while dc-STDP and in-STDP do not.Nowotny et al. have introduced the mechanisms behind the enhancement of neural synchronization by c-STDP rule which rely on the balance of potentiation and depression. The synapse strength remains stable regardless of postsynaptic neuron firing later or earlier attributed to the specific shape of c-STDP curve. The situation of c-STDP is similar to dc-aSTDP. We adopt the similar analysis method [31] here for the dc-aSTDP and dc-STDP. The time lags are recorded as Δt1 and Δt2, where Δt1-Δt2=T1=T21 and Δg1-Δg2=0 (Figure 8) at this state. If post-synaptic neuron fires faster, Δt1 becomes smaller. Synaptic strength will be depressed, due to Δg1-Δg2<0 for dc-aSTDP rule, so that the post-synaptic neuron is less excited and goes back into the synchronized state (Figure 8(b)). The other direction can be analyzed in the same way for dc-aSTDP. But for dc-STDP, when post-synaptic neuron fires faster, synaptic strength will be increased due to Δg1-Δg2>0 in this case. The post-synaptic neuron is more excited and cannot go back into the synchronized state (Figure 8(c)). The opposite direction is the same case for dc-STDP and cannot go back into the synchronized state. Therefore, the synchronization mechanisms between these two rules are different.Different synchronization mechanisms for different STDP rules. (a) Shows the situation ofΔt1-Δt2=T1=T2, where Δg1-Δg2=0. The solid lines in (b) and (c) are same as (a). (b) and (c) present the change of synaptic strength of dc-aSTDP and dc-STDP rules, respectively.
(a)(b)(c)For in-STDP rule, we can easily find thatΔg1-Δg2 is always positive, where Δg1-Δg2=Aplus*(exp(-x/tplus)-0.5)-Asub*(exp(x-171/tsub)-0.5) and values of parameters are the same as Figure 3. It means that the potentiation and depression of synaptic conductance during one period cannot achieve balance. The synaptic strength must achieve the maximum resulting from Δg1-Δg2>0.On the other hand, learning parameters also play important role in neural synchronization. We try to explain the role of a learning parameter by considering how it influences the synapse conductance, which is a major factor for synchronizing neurons with a given mismatch of intrinsic frequencies.It is obvious that if synaptic conductance becomes stronger, it can make largerT2  to achieve the same period with pre-synaptic neuron for various STDP rules. Thus, with other three learning parameters fixed, larger Aplus values, which are corresponding to the stronger stable synaptic strength for c-STDP rule,can cause larger T2 synchronization. Similarly, smaller Asub values will give rise to larger T2 synchronization for c-STDP (Figure 5(a)). Moreover, from Figure 5(b), moderate tplus values can also make larger T2 synchronized to T1 because these values bring about stronger synaptic strength. We can prove this perspective by simple calculus reasoning. Based on the expression of c-STDP rule in Section 2, let tplus be variable and let other parameters keep constant. Δgraw is viewed as the function of variable tplus. By calculating the derivative of Δgraw, we can find that Δgraw is a first increasing and then decreasing function when tplus increases gradually. Thus, medial values of tplus can result in stronger synaptic strength. By the same reasoning for tsub,we can conclude that smaller or larger tsub can make larger T2 synchronization for c-STDP (Figure 5(b)).The effects of learning parameters on synchronization about other learning rules are similar to the c-STDP rule for largerT2. However, smaller T2 values leading to synchrony are only related to c-STDP and dc-aSTDP, because they rely on the balance of depression and potentiation which could lead to an appropriate low stable synaptic strength. Therefore, only proper learning parameters got by regulating the effect of learning parameters on synchronization are required for smaller T2 achieved synchronization.Finally, we conclude why the widest synchronization windows for some STDP rules are different. From the above statements, we find that two aspects affect synchronization. One is the maximum synaptic strength which can make largerT2 synchronize. The other is the balance of depression and potentiation which can make smaller and moderate T2 synchronize. For c-STDP and dc-aSTDP, they can achieve the widest synchronization windows through modulating the two aspects. However, dc- STDP and in- STDP rules, due to their specific shape, can only make use of the first one. This implies that the widest synchronization windows for dc-STDP and in-STDP rules cannot exceed the synchronization windows for the maximum constant connection strength under the same model respectively. Therefore, the optimal synchronization windows for c-STDP and dc-aSTDP are wider than those for dc-STDP and in-STDP.


## 3.1. Comparison of Synchronization Windows of Different Types of STDP

We investigate the width of synchronization window of various STDP curves, with the same set of parameters. The period of the pre-synaptic neuron is chosen to be 171 milliseconds, which falls into the range of theta waves. Several reasons make us choose such a long period. First, it has fairly wide synchronization windows which allow comparisons in a relative precise manner and can provide clearer information about synchronization windows of various STDP learning rules. Second, the slow theta waves always involve many neurons that fire synchronously [38, 39]. Also, theta waves have many interesting implications. For example, theta waves are normally absent in healthy awake adults, but appear during the state of meditation [40].  During emotional arousal and various types of rhythmic activities during sleep, neurons in the amygdala produce theta activity [41, 42]. And it is known that coherent theta activity (4–8 Hz) in amygdala-hippocampal circuits is deeply involved in fear memory [43].With fixed period of the pre-synaptic neuronT1, we evaluate the coupled period of postsynaptic neuron T21 when it is driven by the pre-synaptic neuron. The values of learning parameters used in c-STDP, dc-STDP, dc-aSTDP are Aplus=9 nS, Asub=6 nS, tplus=100 milliseconds, tsub=200 milliseconds, additionallyτ0=30 milliseconds in c-STDP, and in in-STDP are  Aplus=Asub=8 nS, tplus=100 milliseconds, tsub=200 milliseconds. Our model and most values of parameters are from the model of Nowotny for c-STDP, except Aplus, tsyn, Vslope, gmax are different [31]. Especially, T1 is fixed at 171 milliseconds in our simulations while T2 is set to constant value 300 milliseconds of Nowotny’s work.The window of synchronization (upper panel) and quality (middle panel) of dc-STDP are presented in Figure2 as an example. We scan T2 from 150 milliseconds to 320 milliseconds. The upper panel shows the number of synchronization times in 40 simulations, for each T2. It is clear that, in certain range of T2, simulations from different initial values may have different results. Only when T2 falls into the segment from 194 to 221, the post-synaptic neuron can synchronize with the pre-synaptic neuron from any initial value. It is easily found that there are some T2 corresponding to the number of synchronization times between 1 and 39. In this situation, we present the three possible states of post-synaptic neuron’s firing in Figure 2 (lower panel)—keeping the initial period (squares), oscillating (circles), and synchronizing with the pre-synaptic neuron (dots). Obviously, these states are independent of the synchronization criteria we set.For the purpose of discussing the function role of STDP rules in synchronization, the synchronization windows of various type of STDP are plotted in Figure3, in which the case of constant synaptic conductance is also included as a comparative tool. The same parameter values used in simulations ensure a fair comparison. In our simulation studies, the synapse strength is between 0 nS and 25 nS. We choose the maximal synaptic strength and the middle synaptic strength of STDP synapse as the synaptic strength of constant synapse in this study. Interestingly, because the synchronization windows for STDP rules are narrower than synchronization window for constant synapse g=25 nS in Figure 3, these results, opposite to previous reports, indicate that all STDP rules do not enhance synchronization comparing with the constant synapse under the chosen parameters.Figure 3
The synchronization window of different types of learning curves. The top five lines (points) are excitatory synapses and the bottom three lines (circles) are inhibitory synapses. See legends in the figure.Several other points are worthy of detailed describing. First, increasing the excitatory constant synaptic connections from 12.5 nS to 25 nS, leads to a wider synchronization window. However, the wider window could not totally contain the smaller one. It extends toward largerT2, but loses a portion of smaller T2. The case of increasing inhibitory constant synaptic connections is alike. Second, the widest ranges of synchronization window are achieved by the excitatory constant connection g=25 nS and c-STDP rules. However, the lower boundary of the synchronization window of c-STDP is much nearer to T1 than that of excitatory constant synapse. Thirdly, although c-STDP and dc-STDP are fitted from the same set of experimental data, c-STDP has a much wider synchronization window than dc-STDP.We conclude that all STDPs do not give rise to enhanced synchronization and the window of dc-STDP is surprisingly narrow under the chosen parameters in Figure3. Therefore, the questions about what bring about these results inspire us to study further. At the same time, Figure 2 shows a large part of probabilistic synchronization, whose range is a subset of the synchronization window of constant synaptic connection with g=25 nS. We are interested in if the probabilistic synchronization could be enhanced into absolute synchronization by modulation of learning curves. These are the theme of the next section.


## 3.2. The Effect of Learning Parameters on Synchronization

In order to establish the functional role of STDP clearly, we consider if the learning parameters for each STDP rule have important effect on synchronization. In addition, synchrony-asynchrony transition plays important role in the brain. An increase in the degree of synchrony of a uniform input can cause transitions between memorized activity patterns in the order presented during learning. However, if synchronous input is at a low level, transitions cannot occur [44]. The synchrony-asynchrony transition have also been implemented in controlling winner-take-all competition [45], the next recalled time of associative memory [46] and the fine structure of cell assemblies [47].In this section, we will discuss the flexibility of the synchronization window, by exploring regulation of width of the synchronization window,whose boundary indicates the synchrony-asynchrony transition. We take the modulation of learning parameters as the method to regulate the synchronization window. There are four parameters that determine a learning rule Aplus, Asub, tplus, tsub. With three of them fixed and only one parameter changing, we could explore its influence on the width of the synchronization window. For example, Aplus increases from an adequately small value 1 nS to 20 nS with a step of 1 nS. We have also scanned values that are beyond 20 nS, but find that the effect of increasing Aplus is saturated around 20 nS. Further increasing Aplus brings no more effect. Other parameters are fixed as: Asub = 6 nS, tplus = 100 milliseconds, tsub = 200 milliseconds, T2=233 milliseconds. In these conditions, we present the effect of Aplus on the location of the synchronization window with the dc-STDP rule. We carry out simulations 40 times, each from different initial values.Figure4(a) shows the value of ARP (average change of relative period ratio) = 〈T2-〈T21〉〉/(T2-T1) for different Aplus. Some points have value 0 or 1, which means post-synaptic neuron keeping initial period or achieving synchronization with the pre-synaptic one, respectively. Some points have values other than 0 or 1. Figure 4(b) gives an explanation that these points correspond to probabilistic synchronizations with fixed T2=233.It shows that the regulation ofAplus to synchronization for the dc-STDP carrying out 40 times of stimulations. Asub = 6 nS, tplus = 100 milliseconds, tsub = 200 milliseconds, T2 = 233 milliseconds. Top: the value of ARP = 〈T2-〈T21〉〉/(T2-T1) is 0 or 1 which, respectively, means post-synaptic neuron keeping initial period or achieving synchronization. Bottom: the number of synchronization times (ST) in 40 stimulations againstAplus.
(a)(b)We find that the absolute synchronization range is from 10 nS to 20 nS in Figure4(b). According to the definition of synchronization window of T2, we can similarly define 10 nS to 20 nS as the synchronization window of Aplus, with dc-STDP rule and other fixed parameters. The boundary of this synchronization window indicates where synchrony-asynchrony transition happens when changing Aplus.From the results of Figure4, the reason why synchronization window of constant synapse is wider than that of STDP rules (Figure 3) may be explained by learning parameters. To figure out a global picture of the effect of Aplus on synchronization for c-STDP rule, we then determine the synchronization window of Aplus with different T2. Aplus increases from 1 nS to 20 nS with a step of 1 nS while other three parameters keep initial values: Asub=6 nS, tplus=100 milliseconds, tsub = 200 milliseconds. We choose some typical values of T2 to character the global picture. The points in Figure 5(a) show the synchronization range of Aplus with T1 divided by the chosen values of T2. The lower boundaries, as well as those upper boundaries that are other than 20 nS, indicate the position of synchrony-asynchrony transitions. For example, when T2 is 177 milliseconds equivalent to T1/T2 = 0.966, Aplus outside of the points range from 8 nS to 11 nS cannot lead to synchronization between the two neurons. In addition, for those points marked on the horizontal axis, synchronization could not be established no matter what values Aplus take.The range of learning parameters leading to synchronization for all 40 stimulations for c-STDP is presented. The points are the range ofAplus and tplus which can make neuron synchronization. The circles are the range of Asub and tsub. We choose some values of T2 to investigate the effect of Aplus, Asub, tplus and tsub on global synchronization. Aplus and Asub vary from 1 nS to 20 nS. Let tplus and tsub vary from 10 milliseconds to 400 milliseconds. There are some values of T1/T2 marked on the horizontal axis that cannot be entrained to achieve synchronization in Figures 5(a) and 5(b).
(a)(b)According to this global picture, the intersection of those synchronization ranges ofAplus, which is from 10 nS to 11 nS, identifies the range of Aplus that would lead to the widest synchronization window which is from 177 to 289. This optimal synchronization window for c-STDP is wider than constant synapse g = 25 nS. Thus, the learning parameters strongly influence the role of STDP on synchronization.Taking the same method asAplus, we study the effects of other three parameters on synchronization. Figure 5 presents the situation of Asub varying from 1 nS to 20 nS with the step 1 nS (circles, Figure 5(a)), tplus (points, Figure 5(b)) and tsub (circles, Figure 5(b)) both varying from 10 milliseconds to 400 milliseconds with the step 10 milliseconds. When changing one parameter to explore how the range of synchrony evolves with T2, other three parameters keep their initial values as in Figure 3. Similarly, the global modulation picture of four learning parameters for other learning rules can be got. We only give the situation for c-STDP rule in Figure 5.The optimal synchronization windows for various learning rules are presented in Figure6 comparing with previous synchronization windows obtained in Figure 3. The parameters used for optimal synchronization windows are presented in Table 1. The parameters are derived according to regulating one parameter while other three parameters keep initial values. Under the optimal parameters, the synchronization windows for STDP rules are not narrower than constant synapse g = 25 nS. Consequently, the reported important role of STDP in synchronization should be dependent on learning parameters.Table 1
The parameters for optimal synchronization windows.
STDP ruleAplus (nS)Asub (nS)tplus (ms)tsub (ms)c-STDP96100270dc-STDP9610050dc-aSTDP96100350in-STDP88100200Figure 6
We present the optimal synchronization window for various learning rules compared with Figure3. The top and the bottom lines of this figure are the synchronization windows of excitatory and inhibitory constant synapse, respectively. In addition, there are four pair lines in the middle panel for STDP rules. Each pair includes optimal synchronization window (top line) and previous synchronization window in Figure 3 (bottom line). The left range is optimal synchronization range and the range in round bracket is previous synchronization range.For the first three STDP rules, the optimal synchronization windows are got by regulatingtsub while other three fixed learning parameters keep initial values. For the last STDP rule, the initial values are the optimal parameters.


## 3.3. The Synchronization Mechanism

It is important to understand the properties which neural synchronization depends on. We take into account this problem from two aspects.On the one hand, synchronization correlates with the chosen STDP rules. For excitatory synapse, we find that the optimal synchronization windows for c-STDP rule and dc-aSTDP rule are almost equal, and are wider than constant synapses. However, comparing with these two rules, the optimal synchronization window for dc-STDP rule is much narrower. For inhibitory STDP, the synchronization window is the same as that for constant synapse under connection strengthg=25 nS.Accordingly, what makes the optimal synchronization windows for various STDP rules different deserves an explanation. The stationary synaptic conductance is a necessary condition for stationary synchronized state [31]. We find that the mechanisms of synchronization caused by c-STDP (or dc-aSTDP) and dc-STDP (or in-STDP) are different in our model.Figures7(a) and 7(b) show the average Δt=tpostspike-tprespike, and synaptic strength after an episode of coupling time for c-STDP rule. Parameters in Figure 7 are same as those in Figure 3. There are two types of behavior for Δt when synchronization occurs (Figure 7(a)). In a section of constant Δt, the synaptic strength does not achieve the maximal value. Apparently, in this situation, postsynaptic neuron achieves synchronization with the pre-synaptic neuron depending on the balance between potentiation and depression of synaptic conductance. In the rest part of synchronization window, the synaptic strength achieves the maximal value (Figure 7(b)). It indicates that for larger T2, postsynaptic neuron achieves synchronization depending on the effect of maximal synaptic conductance. For dc-aSTDP rule, the synchronization mechanisms are similar to c-STDP rule. When the post-synaptic neuron synchronizes with the pre-synaptic neuron under small T2, the change of synaptic potentiation and depression cancel each other. However, for the small portion of synchronization window at the right side, synaptic conductance gets the maximum at the stationary synchronized state.Four pictures are obtained from the same values of parameters and model of Figure3. The front four pictures take the function of tanh to limit synaptic strength for c-STDP and dc-STDP rules. Figures 7(a) and 7(c) show the average spike time interval of postneuron’s and preneuron’s spike time for c-STDP and dc-STDP over some time after a period time of coupling, respectively. Figures 7(b) and 7(d) present the average synaptic strength for c-STDP and dc-STDP rule, respectively. Each subplot has two dash lines what indicate the boundary of synchronization window.
(a)(b)(c)(d)For dc-STDP and in-STDP rules, the synchronization mechanisms may be different with the above two STDP rules. Because the potentiation and depression of synaptic conductance cancel each other,Δt must be a fixed value for the selected STDP rule. However, for the dc-STDP rule (Figures 7(c) and 7(d)), Δt keeps varying which means that the potentiation and depression of synaptic conductance do not cancel each other at the synchronization state. But, it is easily found that the synaptic conductance is at the stationary maximum for dc-STDP rule. Thus, postsynaptic neuron achieves synchronization completely depending on the effect of the maximal synaptic conductance for dc-STDP rule. The state of in-STDP rule is similar to dc-STDP rule.As a result, neural synchronization mechanism can be different for various STDP rules. For the fewT2 at the right side of the synchronization window, the synaptic conductance achieves the maximum with c-STDP and dc-aSTDP rules. This result is obvious. Because the frequency mismatch is larger, the synapse needs to be stronger to entrain the post-synaptic neuron. But for most part of synchronization window at the left, c-STDP and dc-aSTDP rules rely on the balance of potentiation and depression. Instead of balancing out potentiation and depression during one cycle, dc-STDP and in-STDP rules depend on synaptic strength achieving its maximum. It is important to understand why c-STDP and dc-aSTDP rely on the balance of potentiation and depression while dc-STDP and in-STDP do not.Nowotny et al. have introduced the mechanisms behind the enhancement of neural synchronization by c-STDP rule which rely on the balance of potentiation and depression. The synapse strength remains stable regardless of postsynaptic neuron firing later or earlier attributed to the specific shape of c-STDP curve. The situation of c-STDP is similar to dc-aSTDP. We adopt the similar analysis method [31] here for the dc-aSTDP and dc-STDP. The time lags are recorded as Δt1 and Δt2, where Δt1-Δt2=T1=T21 and Δg1-Δg2=0 (Figure 8) at this state. If post-synaptic neuron fires faster, Δt1 becomes smaller. Synaptic strength will be depressed, due to Δg1-Δg2<0 for dc-aSTDP rule, so that the post-synaptic neuron is less excited and goes back into the synchronized state (Figure 8(b)). The other direction can be analyzed in the same way for dc-aSTDP. But for dc-STDP, when post-synaptic neuron fires faster, synaptic strength will be increased due to Δg1-Δg2>0 in this case. The post-synaptic neuron is more excited and cannot go back into the synchronized state (Figure 8(c)). The opposite direction is the same case for dc-STDP and cannot go back into the synchronized state. Therefore, the synchronization mechanisms between these two rules are different.Different synchronization mechanisms for different STDP rules. (a) Shows the situation ofΔt1-Δt2=T1=T2, where Δg1-Δg2=0. The solid lines in (b) and (c) are same as (a). (b) and (c) present the change of synaptic strength of dc-aSTDP and dc-STDP rules, respectively.
(a)(b)(c)For in-STDP rule, we can easily find thatΔg1-Δg2 is always positive, where Δg1-Δg2=Aplus*(exp(-x/tplus)-0.5)-Asub*(exp(x-171/tsub)-0.5) and values of parameters are the same as Figure 3. It means that the potentiation and depression of synaptic conductance during one period cannot achieve balance. The synaptic strength must achieve the maximum resulting from Δg1-Δg2>0.On the other hand, learning parameters also play important role in neural synchronization. We try to explain the role of a learning parameter by considering how it influences the synapse conductance, which is a major factor for synchronizing neurons with a given mismatch of intrinsic frequencies.It is obvious that if synaptic conductance becomes stronger, it can make largerT2  to achieve the same period with pre-synaptic neuron for various STDP rules. Thus, with other three learning parameters fixed, larger Aplus values, which are corresponding to the stronger stable synaptic strength for c-STDP rule,can cause larger T2 synchronization. Similarly, smaller Asub values will give rise to larger T2 synchronization for c-STDP (Figure 5(a)). Moreover, from Figure 5(b), moderate tplus values can also make larger T2 synchronized to T1 because these values bring about stronger synaptic strength. We can prove this perspective by simple calculus reasoning. Based on the expression of c-STDP rule in Section 2, let tplus be variable and let other parameters keep constant. Δgraw is viewed as the function of variable tplus. By calculating the derivative of Δgraw, we can find that Δgraw is a first increasing and then decreasing function when tplus increases gradually. Thus, medial values of tplus can result in stronger synaptic strength. By the same reasoning for tsub,we can conclude that smaller or larger tsub can make larger T2 synchronization for c-STDP (Figure 5(b)).The effects of learning parameters on synchronization about other learning rules are similar to the c-STDP rule for largerT2. However, smaller T2 values leading to synchrony are only related to c-STDP and dc-aSTDP, because they rely on the balance of depression and potentiation which could lead to an appropriate low stable synaptic strength. Therefore, only proper learning parameters got by regulating the effect of learning parameters on synchronization are required for smaller T2 achieved synchronization.Finally, we conclude why the widest synchronization windows for some STDP rules are different. From the above statements, we find that two aspects affect synchronization. One is the maximum synaptic strength which can make largerT2 synchronize. The other is the balance of depression and potentiation which can make smaller and moderate T2 synchronize. For c-STDP and dc-aSTDP, they can achieve the widest synchronization windows through modulating the two aspects. However, dc- STDP and in- STDP rules, due to their specific shape, can only make use of the first one. This implies that the widest synchronization windows for dc-STDP and in-STDP rules cannot exceed the synchronization windows for the maximum constant connection strength under the same model respectively. Therefore, the optimal synchronization windows for c-STDP and dc-aSTDP are wider than those for dc-STDP and in-STDP.


## 4. Conclusion

STDP plays important functional role in neural synchronization. The mechanism of STDP in neuronal synchronization is still not completely clear. Inspired by previous experiments and theoretical researches, we study the important aspects of STDP-induced synchronization in this paper, such as the role of various STDP rule in synchronization, the widest synchronization window through regulating learning parameters, and synchronization mechanism.In order to explore the functional role of STDP in synchronization, we compare synchronization windows of different types of STDP rules with that of constant synapse under the same model parameters. For the given parameters, not all synchronization windows are enhanced by STDP rules.Synchronized responses have a stronger influence on cells at subsequent processing stages than nonsynchronized responses [48–50]. And the enhanced precise synchronization is important in improving a rapid and reliable transmission of information about sensory changes [2, 51]. Recent researches have reported various methods to enhance synchronization, such as, selective attention [52] and time delay [53]. Here, we present the effect of modulation of learning parameters on synchronization and optimal synchronization window which are not narrower than constant synapse. The optimal synchronization windows by c-STDP and dc-aSTDP rules are much wider than constant synapse. It indicates that the function role of STDP rule in synchronization depends on the learning parameters.The synchronization mechanism is also described here. Different shapes of STDP rule can cause different optimal synchronization windows. The optimal synchronization windows of c-STDP and dc-aSTDP are wider than that of dc-STDP for excitatory synapse. For c-STDP rule and dc-aSTDP rule, a stationary synchronized state completely depends on the balance between potentiation and depression or the maximal synaptic conductance. However, for dc-STDP and in-STDP rule, the stable synchronized state depends on the maximal synaptic conductance under the self-limitation of synaptic strength. If we change the type of bound of synaptic strength for dc-STDP rule from self-limitation to artificial bounds, we find that the synchronization mechanism does not change. In a word, on one hand, the synchronization range of dc-aSTDP and c-STDP can achieve the optimal synchronization window of dc-STDP, depending on the maximal synaptic conductance. On the other hand, dc-aSTDP and c-STDP can extend the synchronization windows to include smallerT2 by the balance between potentiation and depression.The firing pattern of neurons is regular in this paper. Many neurons in brain areas present regular firing. For example, neurons in cat area 17 can be grouped in 4 different electrophysiological cell classes, including regular spiking [54]. And spontaneous, regular action potentials were observed both with cell-attached patch recordings as well as with whole cell current-clamp recordings for cholinergic neurons in the parabigeminal nucleus of the rat midbrain [55]. Neuronal synchronization properties with regular firing neurons have been studied. For example, whether pyramidal neurons in different cortical layers exhibit similar tendencies to synchronize is studied [56]. Based on this point and the functional role of STDP in synchronization, we explore the synchronization windows of various STDP rules from the view of neurons’ regular firing.STDP-mediated synchronization is a remarkably robust phenomenon against strong noise [30, 31]. Although our simulation is not under the noise environment, our results may represent some predictions for STDP-mediated synchronization in noisy environment. In addition, from Figure 2(a), we can clearly see that some synchronization number is between 1 and 39. We estimate that it may be related to phase.Our results present that the range ofT2 values leading to synchrony increases strongly if the constant synaptic connection is increased from 12.5 nS to 25 nS. Nowotny et al. find that the extent of synchrony does not change considerably by doubling the synaptic conductance [31]. Their result is not conflict with our result. There are three parameters that are different, Vslope, tsyn, and gmax. Furthermore, T1 is fixed at 171 milliseconds in our simulation while T2 is fixed at 300 milliseconds in theirs. When T2 is fixed, the range of T1 values leading to synchrony is limited from 0 to T2 no matter how the strength of constant connection changes. However, when T1 is fixed, the range of T2 values leading to synchrony can vary from T1 to very large value due to the increase of the constant connection. At the same time, if we adopt the same parameters with Nowotny’s paper, the similar result can be got. Furthermore, when constant synaptic connection is 0 nS, it is clear that two neurons with different initial periods cannot synchronize. This situation means that the length of synchronization window is 0. Along with the increase of strength of constant connection, some T2 must cause the synchronization. Thus, it is easily found that synchronization window must become wider by increasing the constant connection to some degree.We mainly discuss synchronization for different STDP rules in this paper. The question of how the time windows of various STDP rules are biophysically regulated remains relatively unexplored. There are some experiments using neuromodulators to study the time window for STDP [57]. We are interested in building molecular kinetic equations for STDP to explain our results.It has been proposed that conscious perception depends on the transient synchronization of widely distributed neural assemblies [58]. And long-distance synchronization plays a role in triggering the cognitive processes associated with conscious awareness [59]. The changed learning parameters by neuromodulators may influence the cognitive processes. In addition, some diseases and the function of brain are related with synchronization mentioned above, especially in theta (4–8 Hz) rhythm synchronization during fear memory retrieval [43] which is consistent with what we considered here. Therefore, our work may advance understanding of synchronization to some extent. And we expect that our simulation results will provide some help for related diseases treatment.



---

*Source: 704075-2009-07-16.xml*