# Synchronized dynamics of cortical neurons with time-delay feedback

- Alexandra S Landsman
^{1}Email author and - Ira B Schwartz
^{1}

**1**:2

**DOI: **10.1186/1753-4631-1-2

© Landsman and Schwartz; licensee BioMed Central Ltd. 2007

**Received: **16 March 2007

**Accepted: **05 July 2007

**Published: **05 July 2007

## Abstract

The dynamics of three mutually coupled cortical neurons with time delays in the coupling are explored numerically and analytically. The neurons are coupled in a line, with the middle neuron sending a somewhat stronger projection to the outer neurons than the feedback it receives, to model for instance the relay of a signal from primary to higher cortical areas. For a given coupling architecture, the delays introduce correlations in the time series at the time-scale of the delay. It was found that the middle neuron leads the outer ones by the delay time, while the outer neurons are synchronized with zero lag times. Synchronization is found to be highly dependent on the synaptic time constant, with faster synapses increasing both the degree of synchronization and the firing rate. Analysis shows that pre-synaptic input during the inter-spike interval stabilizes the synchronous state, even for arbitrarily weak coupling, and independent of the initial phase. The finding may be of significance to synchronization of large groups of cells in the cortex that are spatially distanced from each other.

## 1 Background

There is a significant amount of research showing that spike coincidence of neurons encodes information. Singer and colleagues showed not only that networks of neighboring neurons tend to become synchronized, but that tight synchronization occurs at the opposite hemispheres of the brain, in the presence of significant (6 to 8 milliseconds propagation delays) [1]. Zero time-lag synchronization of neural activity has been involved in such important phenomena as ability to recognize objects (by binding different attributes) [2, 3], olfactory discrimination [4], and has even been proposed as one of the neural correlates of consciousness [5]. It was found that visuomotor integration is associated with synchronization of signals recorded from the visual and parietal, and motor and parietal areas of the awake cat [6, 7].

While much of analysis focuses on network synchrony in the absence of time delays, delays are common in neural networks and there is evidence that under certain parameters [8] or coupling architectures, delays may actually contribute to synchronization. In a recent issue of Science [9], a short article describes a recent laser experiment [10] of mutually coupled lasers in a row. It is noted that when only two lasers are coupled, there is a lag in their phases equal to the amount of time it takes light to pass between them. However, when a third laser is added, the outputs of the outer lasers show zero time-lag synchronization. It has been suggested by Wolf Singer that this phenomena could shed light on synchronization of nerve signals in the brain in the presence of delays [9]. In particular, researchers indicate that this zero-lag synchronization effect in three coupled oscillators could clarify how the hemispheres of the brain synchronize [9].

With this in mind, we propose a three neuron coupling scheme, analogous to the three laser experiment described above. The synchronization phenomena of the three laser system has been recently analyzed [11], elucidating the conditions under which synchronization occurs. It was shown analytically and confirmed with numerics that increasing the delays improves synchronization of the outer lasers [11]. As will be shown in the present paper the dependence of synchronization on parameters is somewhat different in the three coupled neurons case, where it is actually the synaptic time constant, rather than delays in propagation that are important for synchronization of the outer neurons. It is perhaps not surprising that the mechanism behind synchronization in this three neuron model is different, since dynamics at each synapse are modeled by two differential equations, rather than much simpler linear coupling of the three laser scheme.

While the basic coupling architecture we choose to consider is the same as the above described laser experiment, in that we have a symmetric system of three mutually delay coupled oscillators, the projections from the middle to the outer neurons are somewhat stronger (with faster synaptic time constant) than the coupling from the outer to the middle neuron. This difference in coupling strength and synaptic time constant was inspired by the hierarchical structure of many networks in the brain, which to a first approximation is a feed-forward network, modulated by a weaker modulatory feedback [12]. For example, there is a strong forward projection from the thalamic LGN area to V1 or from V1 to MT, with weaker modulatory feedback projections that modulate the magnitude of the cell's response [3]. It has been suggested that exessively strong mutually coupled loops in the brain would promote uncontrollable oscillations, such as in epilepsy [12, 13].

The present paper investigates a simple model for synchronization of cortical cells receiving a common time-delayed input from a different area of the brain, such as the case with projections of pyramidal cells to other cortical areas or the thalamus to the cortex, and sending a weaker time-delayed feedback. The central dynamical question addressed considers when the synchronous behavior of such a neural network is stable, particularly in regard to the synaptic coupling strengths, and the synaptic time constant. We find that shorter synaptic time constant of the target cells promotes synchronization, and even increases the firing rate, for the same strength of input.

## 2 Basic Model and Numerics

*V*and

*R*are voltage and recovery variables, respectively, and

*τ*

_{ R }= 5.6 ms. The

*dR*/

*dt*equation is written as a sum of a linear term, for the normal

*N a*

^{+}and

*K*

^{+}currents, and a quadratic terms in

*V*to approximate the transient potassium current contributions [15]. The above model has been optimized to provide an accurate quantitative fit to the shape of a regular spiking neuron potentials obtained from human neocortical neurons [16]. Figure 2 shows the limit cycle of a cortical neuron in Eq. (1) when the applied current,

*I*, is above the bifurcation value, resulting in a saddle-node bifurcation. Due to a quadratic term in the recovery variable, the spike rate of a cortical neuron can be arbitrarily low, for low currents, and increases as

*I*is increased. The above neurons can be coupled by adding an additional term to

*dV*/

*dt*in Eq. (1) proportional to

*g*(

*V*-

*E*

_{ syn }), where

*g*is the synaptic conductance variable,

*E*

_{ syn }= 0 for excitatory synapses and -0.92 for inhibitory. The synaptic conductance,

*g*, is obtained from the following equations, commonly used for synaptic coupling [14],

*H*

_{ step }(

*x*) = 1, if

*x*> 0 and zero if

*x*< 0,

*V*

_{ pre }is the voltage of the presynaptic neuron and

*τ*

_{ syn }is the synaptic conductance time constant. In numerical simulation, Ω = -0.20 mV was chosen [14]. The reason for using two synaptic equations is that, depending on

*τ*

_{ syn }, the conductance will peak after

*V*

_{ pre }, continuing to depolarize the membrane after the end of the presynaptic spike (see Figure 3). This type of response is consistent with physiological data. For brief stimulus spike at

*t*= 0, the two synaptic equations produce a response proportional to [13]

Figure 3 shows *g* for different values of *τ*_{
syn
}. In each case, the area under the curve is the same (equal to one), with *τ*_{
syn
}determining the width of the post-synaptic spike.

where *i* is the index of each neuron, with outer neurons receiving the same current, *I*_{1,3} ≡ *I*_{
s
}, and the same synaptic strength, *δ*_{1} = *δ*_{3} ≡ *δ*. The current input to the middle neuron is higher than to the outer ones: *I*_{2} > *I*_{
s
}, leading to the lower uncoupled spiking rate for the outer cells. This was done so that the higher activity of the inner cell drives the outer ones, as might be the case in a typical hierarchical network [3]. The synaptic input to the middle neuron is weaker, *δ*_{2} <*δ*. The synaptic time-constant in Eq. (2) is the same for the outer neurons, *τ*_{syn1 }= *τ*_{syn3 }≡ *τ*_{
syn
}, but longer for the middle neuron, to model the affect of a slowly-varying modulatory feedback. The conductances, *g*_{
i
}, are obtained from Eq. (2), with presynaptic voltage to the middle cell given by: *V*_{pre2 }= *V*_{1} (*t* - *τ*_{
d
}) + *V*_{3} (*t* - *τ*_{
d
}), and to the outer cells by *V*_{pre1,3 }= *V*_{2} (*t* - *τ*_{
d
}). The delay in the propagation of the signal is given by *τ*_{
d
}.

### 2.1 The effect of delays on correlation and synchronization

*τ*

_{ d }, does not seem to effect the degree of synchronization, it has a substantial effect on correlations and phase relations between neurons. Figure 4 shows correlations in spiking output between the middle and outer neuron, and between outer neurons for

*τ*

_{ d }= 10. This delay time is similar, for instance, to the propagation delay between the two hemispheres of the brain [1]. The

*x*-axis indicates the time-shift at which the correlations function was computed. The outer neurons are synchronized, since

*C*

_{13}= 1 at

*t*= 0. It can be seen that the greatest correlations between the inner and outer neurons occur when

*t*is shifted by the delay time,

*τ*

_{ d }. The delay also creates spikes in correlation at intervals of twice the delay time, as can be seen in Fig. 4. Since, in this case, the outer neurons are synchronized, the spikes indicate that there are self-correlations in the time series of a single neuron at intervals of 2

*τ*

_{ d }. This is the round-trip or feed-back time, since its the minimum time that it would take for a signal to travel from one of the neurons, affect the target, and get back to that same neuron. It follows that time delays lead to self-correlations in the spike trains that are not observed when delays are absent, and which may lead to more regular patterns in the time-series data.

For reasons explained in the introduction, *δ* > *δ*_{2} and *τ*_{
syn
}<*τ*_{syn2 }were used to model the affect of a stronger forward and a weaker modulatory feedback. For this type of coupling, the middle neuron leads the outer ones by the delay time, *τ*_{
d
}. This type of phase-locking behavior has been observed in other time-delay systems, such as lasers [17, 18]. Thus the time-shift between two correlated spike trains should be directly related to the delays in transmission, with the input leading the output by the delay time.

### 2.2 The effect of the coupling strength on firing rate and synchronization

*I*

_{ s }= 0.22 and

*I*

_{2}= 0.5 in simulations) leads to a big difference in firing rate. Figure 5 shows a typical voltage trace for the uncoupled,

*δ*= 0, and coupled neurons,

*δ*= 4, where coupling is sufficiently strong to cause phase-locking. Phase-locking or 1 : 1 frequency locking occurs for

*δ*> 3.4. This type of behavior has been observed for mutually coupled neurons in the absence of delays [14]. For lesser values of the coupling strength, different frequency locked behaviors are observed, with frequency ratio between 1 and

*f*

_{ i }/

*f*

_{ o }, where

*f*

_{ i }and

*f*

_{ o }are the frequencies of inner and outer neurons, respectively, in the absence of synaptic coupling. Figure 6 shows 1 : 4 frequency locking that occurs for weak synaptic coupling:

*δ*= 0.5,

*δ*

_{2}= 0.3. When the coupling is too weak, the outer neurons become desynchronized, as shown in Figure 7. The bifurcation value of

*δ*, however depends on the synaptic time constant,

*τ*

_{ syn }, of the outer-neurons. Thus Figs. 6 and 8, which have a very short time constant of

*τ*

_{ syn }= 0.03, show synchronization at a lower coupling strength of

*δ*= 0.5, compared to a minimum

*δ*= 1.03 needed for synchronization of neurons in Fig. 7, where

*τ*

_{ syn }= 0.5 (a more realistic value for fast synapses). This sensitive dependence of synchronization on the coupling strength and synaptic time constant of the outer-neurons is explored analytically in Section III.

### 2.3 The effect of synaptic time-constant on synchronization and firing rate

*τ*

_{ syn }, is decreased. This increase in firing rate may be surprising, since the time constant only controls the width of the conductance spike and not the area under the curve, as shown in Fig. 3. Thus the contribution of a pre-synaptic spike to a change in post-synaptic voltage during an inter-spike interval is largely independent of the synaptic time constant. This can be seen by integrating the

*g*(

*V*

_{ i }-

*E*

_{ syn }) term in Eq. (4) over the interval of conductance change. Figure 10 shows a fluctuation in conductance for a system given by Eqs. (2) and (4), with

*τ*

_{ syn }= 1 and

*τ*

_{syn2 }= 3. A lower firing rate that occurs for slower synapses may be the result of the decay of any increase in voltage during the inter-spike interval back to the limit cycle trajectory (see the bottom of Fig. 8). A pre-synaptic spike from the inner neuron can trigger a spike from the outer one, when it is delivered toward the end of the inter-spike interval. Thus a more narrow jump in conductance and the resultant jump in voltage,

*V*

_{1,3}, may mean that there is less decay before the critical threshold is reached, thereby increasing the likelihood of a spike.

*τ*is greater synchronization. This is shown in Fig. 11, where the outer neurons become progressively synchronized as

*τ*

_{ syn }is decreased from 0.5 to 0.2. It can be seen that for relatively weak coupling of

*δ*= 1,

*δ*

_{2}= 0, and

*τ*

_{ syn }= 0.2, the outer neurons are completely synchronized, after the transients die out. The next section analyzes the affect of a presynaptic spike on a fast synapse (small

*τ*

_{ syn }) in synchronizing two nearby trajectories of the outer neurons.

## 3 Analysis of synchronization for weakly coupled, fast synapses

*V*

_{1},

*R*

_{1}} and {

*V*

_{3},

*R*

_{3}} are interchanged. It follows that the synchronized regime, where the symmetric variables are exactly equal: {

*V*

_{1}=

*V*

_{3}=

*V*

_{0},

*R*

_{1}=

*R*

_{3}=

*R*

_{0}}, is a solution [19]. For an uncoupled system, this solution would not be a stable one, since any perturbation along the limit cycle would result in a phase-difference. A spiking input,

*V*

_{2}, from the middle neuron stabilized the synchronous state when the synaptic time-constant,

*τ*

_{ syn }is sufficiently short. To study the affect of a single pre-synaptic spike on two identical neurons, with nearby trajectories along the limit cycle, introduce new variables: $\tilde{V}$ =

*V*

_{1}-

*V*

_{3}and $\tilde{R}$ =

*R*

_{1}-

*R*

_{3}. These new variables correspond to a perturbation transverse to the synchronized state: {

*V*

_{0},

*R*

_{0}}. Using Eqs. (1), (3) and (4), the linearized dynamics are:

*n*(

*V*

_{0}) = 26 (

*V*

_{0}+ 0.95),

*h*(

*V*

_{0},

*R*

_{0}) = (-5.03 + 101.4${V}_{0}^{2}$ - 32.45

*V*

_{0}+ 26

*R*

_{0}), and

*m*(

*V*

_{0}) = -(3.8 + 6.6

*V*

_{0}).

*t*

_{ s }denotes the time since the arrival of the last spike from the center to outer neurons. Equation (5) is valid for a short synaptic time constant, when only the effect of the last presynaptic spike is significant (see Fig. 3). During an inter-spike interval of the outer cells,

*n*,

*h*and

*m*are all positive. This can be easily confirmed by using Fig. 2, where -0.8 <

*V*< -0.6 and .18 <

*R*< 0.4 during the inter-spike interval, and calculating the lowest possible values of

*n*,

*h*and

*m*for a given range of

*V*

_{0}and

*R*

_{0}. Since {$\tilde{V}$, $\tilde{R}$} denote the difference between the two nearby trajectories along the limit cycle, there is a relationship between the two variables given by

where *l* is the negative of the slope of the limit cycle at {*V*_{0}(*t*), *R*_{0}(*t*)}. From Fig. 2, *l*(*V*_{0}) ≈ 2 during the inter-spike interval. The above equation is valid as long as two nearby trajectories remain on the limit cycle.

*dt*over to the right-hand side of the equation and integrating. For a small synaptic time constant,

*τ*

_{ syn }≤ 1, the width of a duration of a post-synaptic spike, Δ

*t*, is short compared to the time-scale of neuronal dynamics of outer cells during the inter-spike interval. Here, Δ

*t*measures the duration of a conductance spike, (

*t*

_{ s }/${\tau}_{syn}^{2}$)

*exp*(-

*t*

_{ s }/

*τ*

_{ syn }), which is quite narrow for sufficiently fast synapses (see Fig. 3). Thus most of the change in $\tilde{V}$ during the narrow post-synaptic spike is due to the spike itself. This can be seen by considering a voltage difference of two nearby trajectories at time

*t*+ Δ

*t*in the

**absence**of synaptic input: ${\tilde{V}}_{u}(t+\Delta t)\approx \tilde{V}(t)\mathrm{exp}\phantom{\rule{0.5em}{0ex}}(\Delta t(nl-h))$, with

*n*and

*h*defined after Eq. (6) and

*l*given in Eq. (7). This expression was obtained by substituting Eq. (7) into Eq. (5) and integrating over Δ

*t*while assuming that

*V*

_{0}stays almost constant over Δ

*t*. Since Δ

*t*is small, we have |Δ

*t*(

*nl*-

*h*)| ≪ {1,

*δ*}. So that only the conductance term in Eq. (5) makes a significant contribution to $\tilde{V}$ during the duration, Δ

*t*, of a post-synaptic spike. Dividing Eq. (5) by $\tilde{V}$ and integrating, we get

From Eq. (8), a single presynaptic spike from a center neuron acts to decrease the voltage difference, $\tilde{V}$, of the outer neurons by a factor of *exp*(-*δ*). This presynaptic spike also decreases $\tilde{R}$ by decreasing the positive contribution from the *m*/*l* > 0 term in Eq. (9). Thus the immediate effect of a presynaptic input is to decrease the perturbation from synchronized state of the outer neurons, pushing their trajectories closer in phase-space.

Eqs. (8) and (9) show that a spiking input from the middle neuron has a stabilizing affect on the synchronized state when the synaptic time constant is short. Since the difference in trajectories, $\tilde{V}$, $\tilde{R}$ is taken along a limit cycle, the maximum Lyapunov exponent in the absence of synaptic coupling would be zero (corresponding to the displacement along a trajectory in phase-space), and the transverse exponents must be negative since the trajectory collapses onto a limit cycle. It follows that, for sufficiently small *τ*_{
syn
}, a common synaptic input acting during the inter-spike interval should eventually synchronize the neurons, even for weak synaptic coupling.

Figure 8 shows the synchronization of the outer neurons for a very short synaptic time constant, *τ*_{
syn
}= 0.03 and weak coupling, *δ* = 0.5. After the transients die out, the outer neurons become synchronized, thereby falling on a straight line in the *V*_{1} vs *V*_{3} plot. As can be seen in Figure 6, for weak coupling, there is a big difference in firing rate between the outer and the middle neuron, due to differences in injected current. It follows that the synchronization is not due to phase-locking between the middle and the outer cells. The effect of the spiking input on the limit cycle trajectory can be seen at the bottom of Figure 8. There is an integer ratio between the inner and the outer frequencies, whereby the onset of a spike in the outer neuron is triggered by spiking input from the middle one, delivered toward the end of the inter-spike interval. The phenomena is similar to the subharmonic resonance where the limit cycle responds at a subharmonic of the stimulus frequency [14].

Eqs. (8) and (9) show a sensitive dependence of synchronization on the coupling strength from the middle to the outer neurons. This is confirmed by numerics. Figure 7 shows *V*_{1} vs *V*_{3} for two slightly different values of *δ*, *δ* = 1.02 and *δ* = 1.03. A slight change in *δ* leads to an onset of synchronization between the outer neurons.

## 4 Conclusion

Synchronization of nearby cells is often the result of receiving common input, such as when a pyramidal cell sends projections to a targeted area in the cortex [3]. While pyramidal cells tend to target specific areas, matrix projection cells from the thalamus reach in a diffuse manner into adjacent cortical areas helping to synchronize the activity of large populations of cells [3, 20]. The three neuron scheme investigated in this paper is a simple model for studying these types of hierarchical networks, since it incorporates this phenomena of synchronization of certain areas of the brain (outer neurons) due to common input from a different area (middle neuron), and since synchronization of the two outer neurons would indicate synchronization of many neurons, if coupled to the middle neuron in the same way as the two outer neurons in a three neuron model.

This model of three mutually coupled cortical neurons with delays was studied using analysis and numerical simulation. The outer neurons were stimulated with smaller current and had a much lower firing frequency in the uncoupled case, with their frequency significantly increasing depending on the strength of synaptic coupling with the middle neuron. At higher values of the synaptic coupling constant, typical phase-locked behavior and 1 : 1 frequency locking was found between the middle and the outer neurons, with different frequency locking ratios as the synaptic strength was lowered. It was found that delays affected the time-series data by introducing correlations at the time-scale of the delay. While the spiking behavior in the synchronized case was fairly regular, this effect would be interesting to explore for a more complicated, chaotic spike train that can be achieved by incorporating slow adaptation currents into the neuron model. In the case of phase or frequency locking, the middle neuron leads the outer by the delay time, *τ*_{
d
}.

While synchronization of outer neurons was sensitive to the synaptic strength, the synaptic time constant of outer neurons, *τ*_{
syn
}was also highly significant. It was found that shorter synaptic constant substantially improves correlations, leading to zero-lag synchronization of end neurons even when the coupling strength is very weak. A short synaptic constant was also able to significantly increase firing rate, to the point of inducing 1 : 1 frequency locking with the middle neuron, at a much weaker mutual coupling than would otherwise occur for the input currents used. Analysis of dynamics for fast synapses showed that fast synaptic input during the inter-spike interval stabilized the synchronization manifold, even for arbitrarily weak coupling, and independent of the phase relationship between the inner and outer cells. This indicates that even a very weak synaptic input can synchronize cells, as long as the synaptic time constant is suficiently short. The finding may have significance in synchronizing large groups of cells in the cortex via weak synaptic input from other areas, such as the thalamus, or other areas in the cortex proper.

## Declarations

### Acknowledgements

This work was supported by a grant from the Office of Naval Research. ASL is currently a post doctoral fellow with the National Research Council.

## Authors’ Affiliations

## References

- Engel AK, Konig P, Kreiter AK, Singer W: Interhemispheric synchronization of oscillatory neuronal responses in cat visual-cortex. Science. 1991, 252: 1177-1179. 10.1126/science.252.5009.1177.View ArticleADSGoogle Scholar
- Gerstner W, Kistler W: Spiking Neuron Models. 2002, Cambridge: Cambridge University PressView ArticleMATHGoogle Scholar
- Koch C: The Quest for Consciousness. 2004, Englewood, Colorado: Roberts and Company PublishersGoogle Scholar
- Stopfer M, Bhagavan S, Smith BH, Laurent G: Impaired odour discrimination on desynchronization of odour-encoding neural assemblies. NATURE. 1997, 390 (6655): 70-74. 10.1038/36335.View ArticleADSGoogle Scholar
- Crick F, Koch C: Why neuroscience may be able to explain consciousness. Scientific American. 1995, 273 (6): 84-85.Google Scholar
- Roelfsema PR, Engel AK, Konig P, Singer W: Visuomotor integration is associated with zero time-lag synchronization among cortical areas. Nature. 1997, 385: 157-161. 10.1038/385157a0.View ArticleADSGoogle Scholar
- Konig P, Engel AK, Singer W: Relation between oscillatory activity and long-range synchronization in cat visual-cortex. Proceedings Of The National Academy Of Sciences Of The United States Of America. 1995, 92: 290-294. 10.1073/pnas.92.1.290.View ArticleADSGoogle Scholar
- Dhamala M, Jirsa VK, Ding MZ: Enhancement of neural synchrony by time delay. Physical Review Letters. 2004, 92: 074104-10.1103/PhysRevLett.92.074104.View ArticleADSGoogle Scholar
- Cho A: Nonlinear dynamics – bizarrely, adding delay to delay produces synchronization. Science. 2006, 314: 37-37. 10.1126/science.314.5796.37.View ArticleMathSciNetGoogle Scholar
- Fischer I, Vicente R, Buldu JM, M P, Mirasso CR, Torrent MC, Garcia-Ojalvo J: Zero-lag long-range synchronization via dynamical relaying. 2006, [http://arxiv.org/abs/nlin/0612053]Google Scholar
- Landsmand AS, Schwartz IB: Complete chaotic synchronziation in mutually coupled time-delay systems. Physical Review E. 1997, 75: 026201-10.1103/PhysRevE.75.026201.View ArticleADSGoogle Scholar
- Crick F, Koch C: Constraints on cortical and thalamic projections: The no-strong-loops hypothesis. Nature. 1998, 391 (6664): 245-250. 10.1038/34584.View ArticleADSGoogle Scholar
- Koch C: Biophysics of Computation. 1999, New York: Oxford University PressGoogle Scholar
- Wilson H: Spikes decisions and actions. 1999, New York: Oxford University PressMATHGoogle Scholar
- Rose RM, Hindmarsh JL: The assembly of ionic currents in a thalamic neuron 1. The three dimensional model. Proceedings of the Royal Society of London Series B-Biological Sciences. 1989, 237 (1288): 267-View ArticleADSGoogle Scholar
- Foehring RC, Lorenzon NM, Herron P, Wilson CJ: Correlation of physiologically and morphologically identified neuronal types in human association cortex in vitro. Journal of Neurophysiology. 1991, 66 (6): 1825-1837.Google Scholar
- Heil T, Fischer I, Elsasser W, Mulet J, Mirasso CR: Chaos synchronization and spontaneous symmetry-breaking in symmetrically delay-coupled semiconductor lasers. Physical Review Letters. 2001, 86 (5): 795-798. 10.1103/PhysRevLett.86.795.View ArticleADSGoogle Scholar
- White JK, Matus M, Moloney JV: Achronal generalized synchronization in mutually coupled semiconductor lasers. Physical Review E. 2002, 65 (3 Pt 2A): 036229-10.1103/PhysRevE.65.036229. art. no. 036229View ArticleADSGoogle Scholar
- Pikovsky A, Kurths JM: Synchronization: A universal concept in nonlinear science. 2001, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Shepherd GM: The Synaptic Organization of the Brain. 2004, New York: Oxford University PressView ArticleGoogle Scholar
- Traub RD, Whittington MA, Stanford IM, Jefferys JGR: A mechanism for generation of long-range synchronous fast oscillations in the cortex. Nature. 1996, 383: 621-624. 10.1038/383621a0.View ArticleADSGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.