Comparison and Regulation of Neuronal Synchronization for Various STDP Rules

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.


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].
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 A plus , t syn , V slope , g max (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.

Models and Method
We consider two HH neurons with unidirectional activitydependent 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], where i = 1, 2. Each of the activation and inactivation variables y i (t) = {n i (t), m i (t), h i (t)}, i = 1, 2 satisfies first-order kinetics, The parameters in these equations are given in [31], I stim is a constant input current forcing each neuron to spike with a constant, I stim -dependent period, labeled as T 1 and T 2 . The postsynaptic neuron would show another firing period T 2 1 , when it is driven by the synaptic current, which is dependent on the postsynaptic potential V 2 (t), the reversal potential V rev , the activation variable S(t) and its maximal conductance g(t), where , The time-dependent synaptic coupling strength g(t) nS is Therefore g(t) always have values between 0 nS and g max . 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. g raw is modified by STDP rules that are introduced in the next paragraph. The initial value of g raw is 20 nS. The parameters of the model are g Na = 360 μS, E Na = 50 mV, g K = 70 μS, The time-dependent synaptic coupling strength g(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 neocortexlayer 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].
Δg raw is a function of Δt = t postspike − t prespike , 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 c-STDP: dc-STDP: dc-aSTDP: in-STDP: Synchronization of pre-and post-synaptic neurons occurs when |T 1 − T 2 1 | is limited in an acceptable range. We set the criteria of synchronization as |T 1 − T 1 2 | < 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 V 2 , S could result in different outputs, that is, the post-synaptic neuron sometimes 4 Neural Plasticity 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 T 2 . The standard deviation of |T 1 − T 2 1 |, indicating how precisely the neurons are synchronized, represents the quality of synchronization. Range of T 2 , in which post-synaptic neuron is successfully entrained by the pre-synaptic neuron, that is, |T 1 − T 2 1 | < 1.5 in all 40 simulations, is defined as the synchronization window.

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 neuron T 1 , we evaluate the coupled period of postsynaptic neuron T 2 1 when it is driven by the pre-synaptic neuron. The values of learning parameters used in c-STDP, dc-STDP, dc-aSTDP are A plus = 9 nS, A sub = 6 nS, t plus = 100 milliseconds, t sub = 200 milliseconds, additionallyτ 0 = 30 milliseconds in c-STDP, and in in-STDP are A plus = A sub = 8 nS, t plus = 100 milliseconds, t sub = 200 milliseconds. Our model and most values of parameters are from the model of Nowotny for c-STDP, except A plus , t syn , V slope , g max are different [31]. Especially, T 1 is fixed at 171 milliseconds in our simulations while T 2 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 Figure 2 as an example. We scan T 2 from 150 milliseconds to 320 milliseconds. The upper panel shows the number of synchronization times in 40 simulations, for each T 2 . It is clear that, in certain range of T 2 , simulations from different initial values may have different results. Only when T 2 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 T 2 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 Figure 2(b) presents the quality of synchronization against the ratios of uncoupled periods. In Figure 2(c), we fix T 2 = 265 milliseconds which falls into the probabilistic synchronization window. There are three states of the coupled period of post-synaptic neuron T 2 1 when we carry out 40 stimulation times: (i) keeping the initial period (squares); (ii) oscillating (circles); (iii) synchronizing with the pre-synaptic neuron (dots).
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 Figure 3, 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. 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 larger T 2 , but loses a portion of smaller T 2 . 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 T 1 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 Figure 3. 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.

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-takeall 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 A plus , A sub , t plus , t sub . With three of them fixed and only one parameter changing, we could explore its influence on the width of the synchronization window. For example, A plus 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 A plus is saturated around 20 nS. Further increasing A plus brings no more effect. Other parameters are fixed as: A sub = 6 nS, t plus = 100 milliseconds, t sub = 200 milliseconds, T 2 = 233 milliseconds. In these conditions, we present the effect of A plus on the location of the synchronization window with the dc-STDP rule. We carry out simulations 40 times, each from different initial values. Figure 4(a) shows the value of ARP (average change of relative period ratio) = T 2 − T 2 1 /(T 2 − T 1 ) for different A plus . 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 T 2 = 233.
We find that the absolute synchronization range is from 10 nS to 20 nS in Figure 4(b). According to the definition of synchronization window of T 2 , we can similarly define 10 nS to 20 nS as the synchronization window of A plus , with dc-STDP rule and other fixed parameters. The boundary of this synchronization window indicates where synchronyasynchrony transition happens when changing A plus .
From the results of Figure 4, the reason why synchronization window of constant synapse is wider than that of STDP rules (Figure 3 According to this global picture, the intersection of those synchronization ranges of A plus , which is from 10 nS to 11 nS, identifies the range of A plus 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 as A plus , we study the effects of other three parameters on synchronization. Figure 5 presents the situation of A sub varying from 1 nS to 20 nS with the step 1 nS (circles, Figure 5(a)), t plus (points, Figure 5(b)) and t sub (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 T 2 , 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 Figure 6 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   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.

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 strength g = 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. Figures 7(a) and 7(b) show the average Δt = t postspike − t prespike , 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 T 2 , 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 T 2 , 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.
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 few T 2 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 Δt 1 and Δt 2 , where Δt 1 − Δt 2 = T 1 = T 2 1 and Δg 1 − Δg 2 = 0 ( Figure 8) at this state. If post-synaptic neuron fires faster, Δt 1 becomes smaller. Synaptic strength will be depressed, due to Δg 1 − Δg 2 < 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 Δg 1 − Δg 2 > 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.
For in-STDP rule, we can easily find that 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 Δg 1 − Δg 2 > 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 larger T 2 to achieve the same period with pre-synaptic neuron for various STDP rules. Thus, with other three learning parameters fixed, larger A plus values, which are corresponding to the stronger stable synaptic strength for c-STDP rule,can cause larger T 2 synchronization. Similarly, smaller A sub values will give rise to larger T 2 synchronization for c-STDP ( Figure 5(a)). Moreover, from Figure 5(b), moderate t plus values can also make larger T 2 synchronized to T 1 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 t plus be variable and let other parameters keep constant. Δg raw is viewed as the function of variable t plus . By calculating the derivative of Δg raw , we can find that Δg raw is a first increasing and then decreasing function when t plus increases gradually. Thus, medial values of t plus can result in stronger synaptic strength. By the same reasoning for t sub ,we can conclude that smaller or larger t sub can make larger T 2 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 larger T 2 . However, smaller T 2 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 T 2 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.

Pre-synaptic neuron
Post-synaptic neuron One is the maximum synaptic strength which can make larger T 2 synchronize. The other is the balance of depression and potentiation which can make smaller and moderate T 2 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.

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][49][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 smaller T 2 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 of T 2 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, V slope , t syn , and g max . Furthermore, T 1 is fixed at 171 milliseconds in our simulation while T 2 is fixed at 300 milliseconds in theirs. When T 2 is fixed, the range of T 1 values leading to synchrony is limited from 0 to T 2 no matter how the strength of constant connection changes. However, when T 1 is fixed, the range of T 2 values leading to synchrony can vary from T 1 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 T 2 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.