neurons are known to exhibit a plethora of dynamical behaviors , represented by the generation of action potential patterns .
one of such patterns is bursting , defined by the repeated firing of action potentials followed by quiescent periods .
hence the dynamics of bursting neurons has two timescales : a fast scale related to spiking and a slow scale of bursting itself .
these timescales are related to different biophysical mechanisms occuring at the level of neuron membrane : there are fast ionic currents ( chiefly @xmath0 and @xmath1 ) responsible for spiking activity and slower @xmath2 currents that modulate this activity .
most neurons exhibit bursting behavior if conveniently stimulated .
for example , in the neocortical layer @xmath3 pyramidal neurons , when stimulated with dc current pulses , fire an initial burst of spikes followed by shorter bursts @xcite . in layers @xmath4 , @xmath5 , and @xmath6 chattering neurons fire high - frequency bursts of @xmath7 spikes with a short interburst period @xcite .
cortical interneurons have been found to exhibit bursting as a response to dc pulses @xcite .
pyramidal neurons in the @xmath8 region of hippocampus produce high - frequency bursts after current injection @xcite .
thalamocortical neurons and reticular thalamic nucleus inhibitory neurons exhibit bursting as well @xcite .
purkinje cells in cerebellum can burst when their synaptic input is blocked @xcite .
bursting is also an important feature of sensory systems , because bursts can increase the reliability of synaptic transmission @xcite . in some systems ,
bursts improve the signal - to - noise ratio of sensory responses and might be involved in the detection of specific stimulus features @xcite .
due to both synaptic coupling and common inputs among neurons there are many types of synchronization , which can be generally regarded as the presence of a consistent temporal relationship between their activity patterns @xcite . a strong form of the latter relationship is complete synchronization , where neurons spikes at the same time , i.e. a precise temporal coincidence of events .
a weaker relationship is bursting synchronization , in which only the beginning of bursting is required to occur at the same time , even though the repeated spiking may not occur synchronously .
there has been observed bursting synchronization in cell cultures of cortical neurons , where uncorrelated firing appeared withing the first three days and transformed progressively into synchronized bursting within a week @xcite .
large - scale bursting synchonization in the @xmath9 range has been found in the thalamus during slow - wave sleep , partially originated in the thalamus and gated by modulatory input from the brainstem @xcite .
various areas of the basal ganglia have been found to exhibit bursting synchronization related to parkinson s disease and resting tremor @xcite .
there exists sound neurophysiological evidence that hypokinetic motor symptoms of parkinson s disease such as slowness and rigidity of voluntary movements are closely related to synchronized bursting in the @xmath10 range @xcite .
the connection between bursting synchronization and pathological conditions like parkinson s disease , essential tremor and epilepsy has led to the proposal of many control strategies aiming to suppress or mitigate bursting synchronization @xcite .
one of such strategies is deep - brain stimulation ( dbs ) , which consists of the application of an external high - frequency ( @xmath11 ) electrical signal by depth electrodes implanted in target areas of the brain like the thalamic ventralis intermedius nucleus or the subthalamic nucleus @xcite .
the effect of dbs would be similar to that produced by tissue lesioning and has proved to be effective in suppression of the activity of the pacemaker - like cluster of synchronously firing neurons , and achieving a suppression of the peripheral tremor @xcite .
there is strong clinical evidence that dbs is a highly effective technique for treatment of patients with parkinson s disease @xcite .    in spite of these results ,
dbs is yet far from being completely understood .
many results in this field have been obtained from empirical observations made during stereotaxic neurosurgery , but further progress can be obtained with proper mathematical modelling of dbs @xcite .
the effects of dbs in networks of bursting neurons have been investigated when dbs is implemented through an harmonic external current @xcite and a delayed feedback signal @xcite .    on modelling the response of a neuronal network to an external perturbation like dbs it is of paramount importance to keep the model simple enough such that large - scale simulations ( using a large number of neurons ) can be performed in a reasonable computer time .
in such reductionist point of view a minimal model could be one in which we can assign a geometrical phase to the bursting activity .
the bursting neuron is thus regarded as a phase oscillator undergoing spontaneous oscillations with a given frequency .
thus bursting synchronization becomes a special case of phase synchronization , a phenomenon well understood for coupled oscillators with and without external excitation @xcite .
a simple model for the dynamics of nonlinearly coupled phase oscillators is the kuramoto model , which in its original version considers a global ( all - to - all ) coupling @xcite .
it can be generalized by considering an arbitrary coupling architecture ( generalized kuramoto model ) @xcite .
the particular interest in such models is that many analytical and numerical results are known for them , specially the global case for which a mean - field theory exists for the transition between a non - synchronized to a ( phase-)synchronized behavior @xcite . for generalized kuramoto models it is possible to derive analytical expressions for the critical value of the coupling strength for which the abovementioned transition occurs @xcite .
hence such a body of knowledge can be applied to networks of bursting neurons , helping to design strategies of synchronization control and/or suppression like dbs .
the main goal of this paper is to show , using analytical and numerical arguments , that a system of coupled bursting neurons described by rulkov s model can be reduced to a generalized kuramoto model .
this reduction is valid as long as phase synchronization is concerned , since for frequency synchronization the behaviors can be quite different , though .
we consider , in particular , some widely used connection topologies , like random ( erds - renyi ) , small - world , and scale - free networks .
we show that the analytical results for the critical coupling strength to synchronized behavior , originally derived for the generalized kuramoto model , can be used to describe the synchronization transition also for networks of bursting neurons .    as a matter of fact ,
since bursting activity presents two timescales it can be also approached from the point of view of a relaxation oscillator @xcite . in our work , however , we describe bursting using a single phase . this simplification is justified since phase synchronization of bursting is chiefly related to the slow timescale . in other words , the fast spikes can be nonsynchronized even though the slow dynamics is synchronized .
this paper is organized as follows : in section ii we describe the model we used to describe bursting neurons .
section iii considers networks of coupled bursting neurons using different connection topologies .
section iv reviews some results on the generalized kuramoto model , and section v includes the comparisons we made between kuramoto model and the network of bursting neurons .
our conclusions are left to the final section .
the choose of a suitable model describing the dynamics of biological neurons is dictated by some requirements .
first the model must take into account the kind of dynamics one wishes to describe @xcite .
for example , if all one needs is to describe a spiking neuron , for the sake of neural coding simulations for example , a simple leaky integrate - and - fire ( lif ) model would be enough @xcite .
however if one needs to describe the interplay between different ionic currents flowing through the neuron membrane , the hodgkin - huxley ( hh ) model would be a natural choice @xcite .
on the other hand , the hh model would require far more computational power than the lif since the former involves four complicated first - order differential equations whereas the latter just one simple equation .    with bursting neurons this criterion also holds .
given that bursting results from the interplay between fast and slow ionic currents , hodking - huxley - type models would need at least one more equation to describe slow @xmath12 modulation @xcite .
a model of thermally sensitive neurons exhibiting bursting has been proposed by huber and braun @xcite , which describes spike train patterns experimentally observed in facial cold receptors and hypothalamic neurons of the rat @xcite , electro - receptors organs of freshwater catfish @xcite , and caudal photo - receptor of the crayfish @xcite .
however , the huber - braun model has @xmath3 differential equations for each neuron , and computational limitations impose restrictions to its use for large networks @xcite .    if numerical simulations do not need to take into account the effect of system parameters and only the phenomenological aspects of bursting are relevant , then a good choice is the two - dimensional mapping equations proposed by rulkov @xcite @xmath13}^2 } + y(n ) , \\
\label{rulkovy } y(n+1 ) & = & y(n ) - \sigma x(n ) - \beta,\end{aligned}\ ] ] where @xmath14 is the fast and @xmath15 is the slow dynamical variable , whose values are taken at discrete time @xmath16 , with @xmath17 and @xmath18 .     time evolution of the ( a ) fast and ( b ) slow variables in the rulkov map ( [ rulkovx])-([rulkovy ] ) for @xmath19 , @xmath20.,scaledwidth=60.0% ]    the parameter @xmath21 affects directly the spiking timescale , its values being chosen in the interval @xmath22 $ ] so as to produce chaotic behavior for the evolution of the fast variable @xmath23 , characterized by an irregular sequence of spikes .
the parameters @xmath24 and @xmath25 , on their hand , describe the slow timescale represented by the bursts , and take on small values ( namely @xmath26 ) so as to model the action of an external dc bias current and the synaptic inputs on a given isolated neuron @xcite .
the rulkov model was derived using dynamical rather than biophysical hypotheses .
we choose the parameter @xmath21 so as to yield chaotic behavior for the characteristic spiking of the fast variable @xmath27 [ fig .
[ rulkovfig](a ) ] .
the bursting timescale , on the other hand , comes about the influence of the slow variable @xmath28 , which provide a modulation of the spiking activity due to a saddle - node bifurcation [ fig .
[ rulkovfig](b ) ] . on comparing the dynamics of the rulkov map with similar results of the huber - braun model one is led to an approximate correspondence between variables of the rulkov map and variables with biophysical significance : the discrete time @xmath29 in the map corresponds to @xmath30 of the continuous time ; and @xmath14 stands for the membrane potential , in such a way that each unit of @xmath14 in fig .
[ rulkovfig](a ) corresponds to _ circa _ @xmath31 .
the rulkov model has been used in several numerical investigation of coupled neuron models @xcite    due to the chaoticity of the @xmath14-dynamics of the rulkov map , the duration of each burst suffers a slight variability .
the beginning of each burst can be chosen rather arbitrarily , but it turns out that a useful choice is to consider the local maxima of the variable @xmath15 : let @xmath32 be the time at which the @xmath33th burst begins [ see fig .
[ rulkovfig](b ) ] .
hence the duration of this burst is @xmath34 .
we can define a geometric bursting phase by considering a variable @xmath35 that increases of @xmath36 after a bursting event .
a linear interpolation gives @xcite [ phase ] ( n ) = 2 + 2 , ( n_n n_+1 ) .
we can also define a bursting ( angular ) frequency , which gives the time rate of the phase evolution : @xmath37
in the rulkov model ( [ rulkovx])-([rulkovy ] ) the variable @xmath27 plays the role of the membrane potential at discrete time @xmath16 , where @xmath17 .
hence the difference @xmath38 can be interpreted as the time derivative of the potential .
if the membrane capacitance is scaled to the unity , it amounts to the transmembrane current .
hence the effect of coupling is to inject a synaptic current in the equation for the @xmath14-variable ( [ rulkovx ] ) .
let us denote by @xmath39 the membrane potential of the @xmath40th neuron ( @xmath41 ) .
the equations for a network of coupled neurons are @xmath42}^2 } + y_i(n ) + \varepsilon \sum_{j=1}^{n } a_{ij }   x_j(n ) , \\
\label{y } y_i(n+1 ) & = & y_i(n ) - \sigma x_i(n ) - \beta,\end{aligned}\ ] ] in which @xmath43 is the adjacency matrix , which defines the type of synapse which connects the @xmath40th and @xmath44th neurons ( respectively post - synaptic and pre - synaptic ) .    if the spatial distance between them is taken into account , we can model an electrical synapse by a band - diagonal adjacency matrix which includes only near - neighbors of a given neuron .
chemical synapses are thus represented here by off - band - diagonal elements , since they include the effect of distant neurons .
this is a simplified model , though , since it does not take into account the number of open channels in the post - synaptic neuron . for studies of synaptic - dependent phenomena , like plasticity , memory storage , learning , pattern recognition , etc .
this simplified form of the coupling term may not be sufficient for numerical simulations .
moreover , we consider all interactions as bidirectional , in such a way that the adjacency matrix is symmetric ( @xmath45 ) , what is more likely to be the case in electrical than chemical synapses .
however these simplifications are acceptable in a model as long as we are interested only in the effect of the network topology on the neuron dynamics .
an example is the investigation of the dbs effects related to tremor associated with bursting synchronization in the thalamus @xcite .
the only parameter in our model that is to characterize the intensity of the synaptic connections is the coupling strength @xmath46 .
it takes on positive values for excitatory synapses and negative values for inhibitory ones . since @xmath46 is the same for all neurons in our model we can not consider here the case in which part of the synapses are excitatory ( @xmath47 ) and part inhibitory . a modification such as this ,
although simple to implement in principle , would make the model more difficult to compare with minimal models of phase dynamics .
the only restriction on the values of @xmath46 is that the coupling term itself can not be too large so as to drive the neuron dynamics off a bursting state .
this parameter must be adjusted by trial - and - error .    since in neuronal assemblies
the neurons are likely to exhibit some diversity , we choose randomly the values of @xmath21 within the interval @xmath22 $ ] according to a specified probability distribution function ( pdf ) @xmath48 @xcite . in spite of this we keep the other parameters with the same values ( namely @xmath49 ) since this do not change our results in an appreciable way .
we can compute the bursting phase @xmath50 and the corresponding angular frequencies @xmath51 using ( [ phase ] ) and ( [ frequency ] ) , respectively , for the coupled neurons in the same way as we did for isolated ones .
when the @xmath21 parameter is chosen in the interval @xmath52 $ ] there is a linear relation between @xmath21 and the mean burst frequency .
for this reason when we choose the @xmath21 parameter for each neuron according with a given pdf @xmath48 we generate a pdf for the bursting frequency of uncoupled neurons denoted by @xmath53 with the same shape as @xmath54 . in this work
we consider two types of distributions : ( i ) a waterbag ( or uniform ) distribution : [ uniform ] g_w ( ) = \ {    ccc & & ab + 0 &    . , where @xmath55 and @xmath56 ; and ( ii ) a truncated cauchy distribution : [ cauchy ] g_t ( ) = ^-1 1 + ( ) ^2 ^-1 , where @xmath55 , @xmath56 , @xmath57 is the position of the distribution peak , and @xmath58 is the half - width at half - maximum of the pdf @xcite
.     time evolution of ( a ) mean field and ( c , e ) fast variables of two selected neurons in a network of uncoupled rulkov neurons .
( b , d , f ) stand for coupled neurons.,scaledwidth=70.0% ]    we shall defer the discussion of the nature and properties of the adjacency matrix to the following section . for the moment
let us assume that a convenient form of @xmath43 has been given to the system of coupled rulkov neurons ( [ x])-([y ] ) .
one distinctive effect of synaptic coupling is bursting synchronization : two or more neurons begin their bursting activity at the same times ( up to a given tolerance ) , regardless of whether or not the ensuing sequence of spikes coincides .
a example is provided by fig
. [ bs ] , where we consider two randomly selected neurons in a network : in the uncoupled case the neurons burst at different and uncorrelated times [ fig .
[ bs](c ) and ( e ) ] , what makes the mean field [ meanfield ] x(n ) = _
i=1^n x_i(n ) , to have small - amplitude fluctuations [ fig .
[ bs](a ) ] . for large enough values of @xmath46 the coupled system of neurons present bursting synchronization as illustrated by the two selected neurons [ fig .
[ bs](d ) and ( f ) ] .
the mean field in this case undergoes large - amplitude oscillations with the same frequency as of the bursting neurons [ fig .
[ bs](b ) ] .
another quantitative characterization is provided by the order parameter magnitude and its time average given , respectively , by @xcite [ order ] r(n ) = _
i=1^n e^i_i(n ) , = _ n=1^n r(n ) , where @xmath59 is chosen so as to yield stationary values of @xmath60 . if the neurons burst exactly at the same time their phases coincide and hence the normalized sum in ( [ order ] ) gives @xmath61 . on the other hand ,
if the neuron bursting is so uncorrelated that the times at which they burst are practically random the summation gives nearly zero and @xmath62 for such an extreme case .
in finite networks there is likely to be chance correlations , hence we consider @xmath63 as a threshold for the transition from a non - synchronized to a partially synchronized state .
intermediate cases ( @xmath64 ) thus represent partially synchronized states .
neurons are connected by axons and dendrites , so that we can regard those neurons as embedded in a three - dimensional lattice . however , due to the high connectivity of the neurons it is necessary to use complex networks to describe the properties of neuronal assemblies .
these complex networks can be viewed in two basic levels : a microscopic , neuroanatomic level , and a macroscopic , functional level .
studies in the former level are limited to those few examples in which there is available data on the neuronal connectivity , as the worm _ c. elegans
_ @xcite .    in the second ( macroscopic )
level of description of neural networks , the use of non - invasive techniques as electroencephalography , magnetoencephalography and functional magnetic resonance imaging provides anatomical and functional connectivity patterns between different brain areas @xcite .
this information provides a way to study the brain cortex , considering the latter as being divided into anatomic and functional areas , linked by axonal fibers .
scannell and coworkers have investigated the anatomical connectivity matrix of the visual cortex for the macaque monkey and the cat @xcite .
the basic elements of a complex network are nodes and links .
a link connects two or more nodes , and these connections can be unidirectional or bidirectional . in the language of neuronal networks a node can be either a single neuron or a cortical area , depending on the level of characterization .
if we consider a microscopic description the links are synapses ( electrical or chemical ) , whereas they are axonal fibers in networks of cortical areas or simply they stand for the functional relationship between the corresponding regions in the cortex .
a complex network can be characterized by various topological and metrical properties @xcite . for our purposes
we will need only two of them : the average path length @xmath65 and the average clustering coefficient @xmath66 .
the path length between two nodes @xmath67 and @xmath68 is the minimum number of links that must be traversed through the network to travel from @xmath67 to @xmath68 .
the average path length @xmath65 is the path length averaged over all pairs of nodes in the network .
the average clustering coefficient @xmath66 of a network is defined as the average fractions of pairs of neighbors of a given node which happen also to be neighbors of each other @xcite .
for example , if a node @xmath67 is connected to nodes @xmath68 and @xmath69 , and if @xmath68 and @xmath69 are themselves connected , then @xmath67 , @xmath68 and @xmath69 form a triangular cluster .
the value of @xmath66 turns to be the fraction of triangular clusters that exist in the network with respect to the total possible number of such clusters .    from the graph - theoretical perspective
a complex network can be described by an adjacency matrix @xmath43 , whose elements are equal to the unity if the nodes are connected and zero otherwise . in studies of cortical area networks , both in the anatomic and functional levels , each link is sometimes given a specific weight , which can be even time - varying in investigations of plasticity .
since we discard self - interactions , the diagonal elements of the adjacency matrix are zero .      in a globally coupled network
all nodes are connected with all other nodes . since in this case
the contribution of the coupling would increase with the number of nodes , the coupling term is usually divided by the number of nodes . in such a network each node can be regarded as being coupled to the mean field of other nodes .
the adjacency matrix of global networks is constant : @xmath70 for all @xmath71 .
there have been numerical investigations of such networks with respect to strategies of dbs @xcite
. the coupling term in ( [ x ] ) will be divided by the number of neurons , hence we shall redefine it as [ xi ] =     ( color online ) variation of the order parameter magnitude with the coupling strength for a globally coupled ( a ) , and a randomly coupled ( b ) network of @xmath72 rulkov neurons and probability @xmath73 .
black and red lines stand , respectively , for a waterbag and a truncated cauchy distribution of values of the @xmath21 parameter in the @xmath74 $ ] .
the green curves represent a polynomial fit given by eq .
( [ rfit ] ) .
the insets zoom the behavior near the transition to bursting synchronization.,scaledwidth=100.0% ]    the variation of the order parameter magnitude of the coupled network with the coupling strength above is depicted in fig .
[ g](a ) , where we plot @xmath75 as a function of @xmath76 when the values taken on by the parameter @xmath77 are drawn from a waterbag pdf ( black line ) and a truncated cauchy distribution ( red line ) . since we used randomly chosen initial conditions for the neurons
@xmath78 we considered averages over @xmath79 different realizations of the initial pattern .    in both cases the behavior is qualitatively similar : for small coupling values there is no bursting synchronization and , after a critical coupling value @xmath80 we begin to observe partial synchronization and , further on , complete phase synchronization .
the approach to the latter differs with respect to the pdf @xmath54 , being slightly faster for a truncated cauchy distribution than for a waterbag pdf [ see the inset of fig . [ g](a ) ] .
the `` physical '' distance between two neurons does not play any role in globally coupled networks , since the contribution is counted in the same way for all nodes .
this deficiency has been circumvented by considering a network in which the coupling strength decreases with the lattice distance as a power - law , with exponent @xmath81 @xcite . if @xmath82 the lattice distance does not matter and we recover the globally coupled case .
as the exponent @xmath81 tends to infinity we approach a locally coupled network , where only the nearest - neighbor nodes have to be taken into account .
unlike global networks , random network present a typically small density of links .
random ( erds - renyi ) networks are obtained by from @xmath83 initially uncoupled nodes and building @xmath84 links between randomly chosen pairs with a uniform probability @xmath85 ( avoiding self - interactions ) @xcite .
since the total possible number of links is @xmath86 one has [ prober ] p = .
the average path length of er networks is typically very small since it scales as the logarithm of the network size ( @xmath87 ) .
moreover the clustering coefficient of er networks is likewise small since @xmath88 .
since the contribution of the coupling term is typically small we do not need here to rescale it as we did before .
we show , in fig .
[ g](b ) , the variation of the order parameter magnitude with the coupling strength @xmath46 , showing a qualitatively similar picture but with a considerably lower critical coupling strength @xmath89 for the transition from a non - synchronized to a partially synchronized bursting .
if the probability @xmath85 is too small there may not occur bursting synchronization at all , and we require a minimum number of links to observe such phenomenon , as shown in fig .
[ g](b ) , where @xmath73 and we have @xmath90 links .      random networks of the erds - renyi type represent one limit case of a spectrum of networks of which the other end comprises regular networks , which are lattices of @xmath83 nodes for which each node has links to its @xmath91 nearest and next - to - the - nearest neighbors , hence there are local connections only .
the average path length of one - dimensional regular lattices is @xmath92 and the clustering coefficient is given by @xmath93/[4(z-1)]$ ] @xcite .    in between those limiting cases we have small - world networks , for which we typically have @xmath94 and @xmath95 .
it is possible to obtain them from regular lattices essentially by two procedures : ( i ) watts - strogatz ( ws ) and ( ii ) newman - watts ( nw ) ones .
ws networks are obtained from regular lattices by going through each of its links and , with some specified probability @xmath85 , randomly rewiring that link , moving one of its ends to a new node randomly chosen from the rest of the lattice @xcite .
the coordination number of the network is still @xmath91 on average , but the number of links per node may be greater or smaller than @xmath91 due to the possible existence of nonlocal shortcuts , in addition to the local links .
one possible disadvantage of this construction is that , if @xmath83 is small enough , the rewiring process may create clusters disconnected to the rest of the network , for which nodes the path length would be obviously infinite .
the latter problem can be circumvented in nw networks , that are constructed similarly to the ws ones , but the nonlocal shortcuts are added to the lattice with probability @xmath85 , instead of being rewired @xcite .    a small - world network of either ws or nw types is essentially described by the probability of nonlocal shortcuts @xmath85 . by computing both @xmath65 and @xmath66 we can get a range of @xmath85 for which the small - world property is fulfilled by the network .
the small - world property requires the network to exhibit a relatively small path length while retaining an appreciable degree of clustering , i.e. the following conditions @xmath94 and @xmath95 may hold for many different types of networks of the ws or nw types @xcite .
the adjacency matrix for such networks is symmetric and band diagonal ( with zero diagonal elements ) and the non - diagonal parts are sparse , most of their elements being equal to zero .
( color online ) ( a ) dependence of the ratios @xmath96 ( black line ) and @xmath97 ( red line ) with the probability of shortcuts in a small - world network of the newman - watts type with @xmath72 and @xmath98 ; ( b ) dependence of the merit figure @xmath24 with @xmath85 ; ( c ) variation of the order parameter magnitude with the coupling strength for a small - world network of rulkov neurons with probability @xmath99 .
black and red lines stand , respectively , for a waterbag and a truncated cauchy distribution of values of the @xmath21 parameter in the @xmath22 $ ] .
the green curve represents a polynomial fit given by eq .
( [ rfit ] ) .
the insets zoom the behavior near the transition to bursting synchronization.,scaledwidth=90.0% ]    if the probability of shortcuts is zero we have the average path length and cluster coefficient of a regular networks , denoted as @xmath100 and @xmath101 , respectively , for a network of @xmath102 neurons with @xmath103 local connections . in fig .
[ s](a ) we plot the ratios @xmath96 and @xmath97 as a function of the probability of randomly chosen shortcuts @xmath85 in the nw procedure .
the small - world property holds as long as the ratio @xmath96 is large and @xmath97 is small , what yields an interval @xmath104 $ ] . a further confirmation of this interval is provided by the merit figure @xmath24 , defined by [ sigma ] = , where [ ratios ] = , = , in such a way that the larger is @xmath24 the better the small - world property holds for the network .
we plot the merit figure as a function of @xmath85 in fig .
[ s](b ) , which indeed assumes larger values in the interval @xmath104 $ ] .
the variation of the order parameter magnitude for a small - world network with @xmath99 is depicted in fig .
[ s](c ) as a function of @xmath46 , showing a transition to partial bursting synchronization for @xmath105 .      in scale - free networks
the connectivity @xmath106 ( the number of links per node @xmath40 ) satisfies a power - law pdf [ sf ] p(k ) ~k^- , ( > 1 ) , in such a way that highly connected neurons are connected , on the average , with highly connected ones , a property also found in many social and computer networks @xcite .
functional magnetic resonance imaging experiments have suggested that some brain activities can be assigned to scale - free networks , with a scaling exponent @xmath107 between @xmath108 and @xmath109 , with a mean connectivity @xmath110 @xcite .
in fact , this scale - free property is consistent with the fact that the brain network increases its size by the addition of new neurons , and the latter attach preferentially to already well - connected neurons .    in this paper
we use the barabsi - albert ( ba ) coupling prescription through a sequence of steps , starting from an initial random ( erds - rnyi ) network of size @xmath111 nodes and @xmath112 random connections @xcite .
every step we add a new node to the network which makes two connections : the first is determined by a uniform probability of connection among the vertices in the network and the second link is chosen such that the probability of connection decays according with the degree of connectivity of each node .
hence this second link is more likely to be attached to the most connected node in the network than to the less connected one .
when the network reaches the desired size @xmath83 we stop adding new nodes .
( color online ) ( a ) probability distribution function for the number of connections per site in a network obtained from ba procedure and @xmath72 nodes ; ( b ) variation of the order parameter magnitude with the coupling strength for a small - world network of rulkov neurons with probability @xmath99 .
black and red lines stand , respectively , for a waterbag and a truncated cauchy distribution of values of the @xmath21 parameter in the @xmath74 $ ] .
the green curve represents a polynomial fit given by eq .
( [ rfit ] ) .
the insets zoom the behavior near the transition to bursting synchronization.,scaledwidth=100.0% ]    we applied the ba procedure until the final network has @xmath72 nodes .
we then obtained a numerical approximation for the pdf of connections per node @xmath113 , depicted in fig .
[ law](a ) .
the solid red line is a least squares fit giving a power - law dependence of the form ( [ sf ] ) with an exponent @xmath114 , confirming that the network is scale - free indeed .
the variation of the order parameter magnitude with @xmath46 is plotted in fig .
[ law](b ) , illustrating the transition to bursting synchronization with @xmath115 .
the behavior of the average order parameter magnitude with @xmath46 after the transition to bursting synchronization at @xmath116 can be fitted , when the frequency distribution is a truncated cauchy pdf , by the following expression [ rfit ] ^s , where the exponents are different according to the type of network [ table [ fitab ] ] .
.[fitab ] critical coupling strength and fitting parameters for some types of networks with @xmath102 nodes and a truncated cauchy distribution . [ cols="^,^,^,^,^",options="header " , ]     after having chosen adequate values of the parameters appearing in the pdf of both rulkov and kuramoto models , we can compare the results for them with respect to the transition to bursting synchronization . in figs .
[ compfig](a ) to ( d ) we compute the order parameter magnitude for the rulkov ( black line : waterbag , red line : cauchy ) and kuramoto ( green line : waterbag , blue line : cauchy ) as a function of the coupling strength for the network topologies considered in this paper .
the evolution of @xmath60 is nearly the same for both models , with small differences .
for example the value of @xmath60 in the gkm reaches its maximum slightly before ( i.e. with smaller coupling strengths ) the rulkov network does .
this precedence of gkm in comparison with rulkov is observed for all connection topologies we have considered in this work .
( color online ) order parameter magnitude vs. coupling strength for a network of rulkov neurons ( black : waterbag , red : cauchy ) and kuramoto oscillators ( green : waterbag , blue : cauchy ) with @xmath72 nodes , for ( a ) global coupling , ( b ) erds - rnyi ; ( c ) small - world , and ( d ) scale - free.,scaledwidth=100.0% ]
in this work we propose minimal models for the description of bursting neurons using some of the most used connection topologies for neuronal networks which can be used in numerical experiments .
one of the simplest models that describe bursting is a two - dimensional discrete - time mapping proposed by rulkov , and we considered networks of coupled rulkov neurons using global , random ( erds - renyi ) , small - world ( newman - watts ) , and scale - free ( barabasi - albert ) topologies .
we can define a geometrical phase for coupled rulkov neurons , such that bursting synchronization is equivalent to some form of phase synchronization . using the time rate of the bursting phase we can define a bursting frequency and thus investigate frequency synchronization in such networks .    using a complex order parameter
we quantify the state of phase synchronization of the system . as the coupling strength
is increased , there is a transition from a non - synchronized to a partially phase - synchronized state , which evolves to a completely phase - synchronized one for larger values of the coupling strength .
we computed the corresponding critical values of the coupling strength for the connection topologies considered in this work .
the existence of a phase reduction for the coupled rulkov neuron system suggests the usefulness of a generalized kuramoto model which , in the present context , can be considered a minimal model for bursting neurons .
we presented in this paper some quantitative results that justify this claim . in order to compare both models , however , it is necessary to adjust the corresponding probability distribution functions for the natural frequencies .
we did so by using theoretical estimates for the critical coupling strength for the generalized kuramoto model in different connection topologies . after choosing the pdf parameters in a suitable way
the behavior of the order parameter is very similar for both the rulkov and kuramoto models .
in this appendix we show that , in the scenario near a phase synchronized regime , one can perform a phase reduction from a network of coupled oscillators to a generalized kuramoto model , what explains why the properties of a coupled rulkov neuronal network are similar to the kuramoto model .
we stress that this similarity does not necessarily apply to the situation near a frequency synchronized regime , since the dynamical requirements are different .
our presentation follows closely that in ref . @xcite and the general treatment given in ref .
@xcite .
let a general @xmath117-dimensional flow be given by [ flow ] = * f*(*x * ) , where @xmath118 is a @xmath117-dimensional vector in the system phase space and @xmath119 a vector field .
we assume that there is a stable period-@xmath120 orbit [ orbit ] * x*_0(t ) = * x*_0(t+t ) .
the `` slow '' dynamics along this periodic orbit can be described by a phase @xmath121 which time evolution is given by [ phasedef ] ( * x * ) = _ * x * ( * x * ) = 1 .
now let us consider a network of coupled oscillators with a slight mismatch in their parameters described by [ coupled ] _
i = * f*(*x*_i ) + * f*_i(*x*_i ) + _ j a_ij * v*(*x*_i,*x*_j ) , ( i=1,2,
n ) , where @xmath122 is different for each oscillator and stands for the vector field part containing slightly mismatched parameters , @xmath123 is the coupling strength , @xmath124 the adjacency matrix , and @xmath125 is a coupling function .    from ( [ phasedef ] ) , the slow phase @xmath126 of the coupled oscillators is implicitely defined by the function [ z ] * z*(_i ) = _ * x * ( * x*_0(_i ) ) , where @xmath127 is a period-@xmath120 stable orbit .
the time evolution of the phase is then governed by = 1 + _
j a_ij * z*(_i ) * v*(*x*_i,*x*_j ) + * z*(_i)*f*_i(*x*_i ) , for @xmath128 .    on introducing an auxiliary phase @xmath129 and using a time average over a period @xmath120 , by keeping @xmath130 fixed we can write , in a first - order approximation , an equation governing the time evolution of @xmath130 : [ psi ] _
i = w_i + _ j a_ij ( _ i,_j ) , where [ w ] w_i & = & _ 0^t * z*(t+_i)*f*_i(*x*_i ) dt , + [ gama ] ( _ i,_j ) & = & _ 0^t * z*(t+_i)*v*(*x*_0(t+_i),*x*_0(t+_j ) ) dt , play the roles of frequencies and coupling functions , respectively , for the auxiliary phases @xmath130 .    if we are close to a global phase synchronized state , for which @xmath131 for any pairs of oscillators @xmath132 , we have an approximate form for the coupling function @xmath133 .
introducing the time - dependent variable @xmath134 and supposing that @xmath135 there results [ gama2 ] ( _ i,_j ) = _ 0^t+_j * z*(+_i-_j ) * v*(*x*_0(+_i-_j),*x*_0(+_j ) ) dt .    by expanding the coupling function @xmath136 in a power series and assuming that @xmath137 is nearly constant over the periodic orbit @xmath127 , we get [ gama3 ] ( _ i,_j ) a + b ( _ i-_j ) , where [ a ] a & = & _ 0^t * z * * v*(*x*_0(),*x*_0 ( ) ) d , + [ b ] b & = & _ 0^t * z * ( _ * x**v*(*x*,*y * ) ) _ * x*=*y*=*x*_0 ( ) d , notice that this approximation holds whenever [ cond ] _ * x**v*(*x*,*y * ) 0 for any @xmath118 and @xmath138 belonging to the periodic orbit @xmath127 .    substituting ( [ gama3 ] ) into ( [ psi ] ) yields [ psi1 ] _ i _ i + s_i a + b _ j a_ij ( _ i-_j ) , where @xmath139 is the intensity of the @xmath40th node of the network . as long as we deal with networks
for which the intensities @xmath140 present only a small variation over the network , we can take @xmath140 as practically constant , i.e. , independent of @xmath40 .
moreover , since we are by hypothesis near a global phase synchronized state , the phase difference @xmath141 is small enough to justify the replacement @xmath142 , in such a way that the equation governing the time evolution of the auxiliary phases ( near a phase - synchronized situation ) is a generalized kuramoto model [ kuramoto1 ] _
i w_i + _ j a_ij ( _ i-_j ) , where @xmath143 and @xmath144 .
this work was made possible through financial support from the brazilian research agencies , capes and cnpq .
fabiano ferrari thanks to ruedi stoop and the institute of neuroinformatics ( uzh , ethz ) for having received him as a visiting student during 2014 through the program `` science without boundaries '' .
we thank dr .
j. c. p. coninck for fruitful discussions and helpful suggestions .
00 b. w. connors and m. j. gutnick , _ intrinsic firing patterns of diverse neocortical neurons _ , trends in neuroscience * 13 * , 99 - 104 ( 1990 ) d. a. blank and r. stoop , _ collective bursting in populations of intrinsically nonbursting neurons _ , z. naturforsch . a , * 54 * , 617 - 627 ( 1999 ) . c. m. gray and d. a. mccormick , _ chattering cells : superficial pyramidal neurons contributing to the generation of synchronous oscillations in the visual cortex _ , science * 274 * , 109 - 113 ( 1996 ) .
r. stoop , d. blank , a. kern , j.j .
vyver , m. christen , s. lecchini , c. wagner , _ collective bursting in layer iv synchronization by small thalamic inputs and recurrent connections _
, cognitive brain research , * 13 * , 293 - 304 ( 2002 ) .
h. markram , m. toledo - rodriguez , y. wang , a. gupta , g. silberberg , and c. wu , _ interneurons of the neocortical inhibitory system _ , nature reviews neuroscience * 5 * , 793 - 807 ( 2004 ) .
h. su , g. alroy , e. d. kirson , and y. yaari , _ extracellular calcium modulates persistent sodium current - dependent burst - firing in hippocampal pyramidal neurons _ , journal of neuroscience , * 21 * , 4173 - 4182 ( 2001 ) .
e. j. ramcharan , j. w. gnadt , and s. murray sherman , _ burst and tonic firing in thalamic cells of unanesthetized , behaving monkeys _ , visual neuroscience * 17 * , 55 - 62 ( 2000 ) .
m. d. womack and k. khodakhah , _ active contribution of dendrites to the tonic and trimodal patterns of activity in cerebellar purkinje neurons _ , j. neurosci . * 22 * , 10603 - 10612 ( 2002 ) .
r. krahe and f. gabbiani , _ burst firing in sensory systems _ , nature reviews neuroscience * 5 * , 13 - 23 ( 2004 ) .
w. metzner , c. koch , r. wessel , and f. gabbiani , _ feature extraction by burst - like spike patterns in multiple sensory maps _ , j. neurosci . * 18 * , 22832300 ( 1998 ) .
r. c. elson , a. i. selverston , r. huerta , n. f. rulkov , m. i. rabinovich and h. d. i. abarbanel . _ synchronous behavior of two coupled biological neurons .
* 81 * , 56925695 ( 1998 ) .
p. varona , j. j. torres , h. d. i. abarbanel , m. i. rabinovich , and r. c. elson . _
dynamics of two electrically coupled chaotic neurons : experimental observations and model analysis . _ biological cybernetics * 84 * , 91101 ( 2001 ) .
v. makarenko and r. llins .
_ experimentally determined chaotic phase synchronization in a neuronal system .
sci . * 95 * , 1574715752 ( 1998 ) .
m. i. rabinovich , p. varona , j. j. torres , r. huerta and h. d. i. abarbanel .
_ slow dynamics and regularization phenomena in ensembles of chaotic neurons_. physica a , * 263 * , 405414 ( 1999 ) .
h. kamioka , e. maeda , y. jimbo , h. p. c. robinson , and a. kawana , _ spontaneous periodic synchronized bursting during formation of mature patterns of connections in cortical cultures _ , neuroscience letters , * 206 * , 109112 ( 1996 ) . m. steriade , d. a. mccormick , and t. j. sejnowski , _ thalamocortical oscillations in the sleeping and aroused brain _ , science ,
* 262 * , 679685 ( 1993 ) .
m. bevan , p. magill , d. terman , j. bolam , and c. wilson , _ move to the rhythm : oscillations in the subthalamic nucleus - external globus pallidus network _ , trends neurosci . *
25 * , 525 - 531 ( 2002 ) .
p. j. uhlhaas and w. singer , _ neural synchrony in brain disorders : relevance for cognitive dysfunctions and pathophysiology _ , neuron * 52 * , 155 - 168 ( 2006 ) .
p. brown , _ abnormal oscillatory synchronisation in the motor system leads to impaired movement _ , curr opin neurobiol , * 17 * , 656 - 664 ( 2007 ) .
w. d. hutchison , j.o .
dostrovsky , j.r .
walters and _ et al .
_ , _ neuronal oscillations in the basal ganglia and movement disorders : evidence from whole animal and human recordings _ , j. neurosci . * 24 * , 9240 - 9243 ( 2004 ) . c. park , r. m. worth , and l. l. rubchinsky , _ fine temporal structure of beta oscillations synchronization in subthalamic nucleus in parkinson s disease _ , j. neurophysiol . * 103 * , 2707 - 2716 ( 2010 ) .
c. hammond , c. , h. bergman , and p. brown , _ pathological synchronization in parkinson s disease : networks , models and treatments _ , trends neurosci . *
30 * , 357 - 364 ( 2007 ) .
a. l. benabid , p. pollak , c. gervason , d. hoffmann , d. m. gao , m. hommel , j. e. perret , and j. de rougemont , _ long - term suppression of tremor by chronic stimulation of the ventral intermediate thalamic nucleus _
, lancet * 337 * , 403 - 406 ( 1991 ) .
s. blond , d. caparros - lefevre , f. parker , r. assaker , h. petit , j .-
d . guieu , and j .- l .
christiaens , _ control of tremor and involuntary movement disorders by chronic stereotaxic stimulation of the ventral intermediate thalamic nucleus _
, journal of neurosurgery * 77 * , 62 - 68 ( 1992 ) .
m. c. rodriguez - oroz , j. a. obeso , a. e. lang , j .- l . houeto _
_ , _ bilateral deep brain stimulation in parkinson s disease : a multicentre study with @xmath6 years follow - up _ , brain * 128 * , 2240 - 2249 ( 2005 ) .
a. albanese and l. romito , _ deep brain stimulation for parkinson s disease : where do we stand ?
_ , frontiers in neurology , * 2 * , 33 ( 2011 ) .
p. a. tass , _ a model of desynchronizing deep brain stimulation with a demand - controlled coordinated reset of neural subpopulations _ , biol .
cybern . * 89 * , 81 - 88 ( 2003 ) j .- p .
pfister , p.a .
tass , _ stdp in oscillatory recurrent networks : theoretical conditions for desynchronization and applications to deep - brain stimulation _ , frontiers in computational neuroscience , * 4 * , 22 ( 2010 ) . c. hauptmann , o. popovych , p.a .
tass , _ effectively desynchronizing deep brain stimulation based on a coordinated delayed feedback stimulation via several sites : a computational study _ , biological cybernetics * 93 * , 463 - 470 ( 2005 ) . c. a. s. batista , a. m. batista , j. c. a. de pontes , r. l. viana and s. r. lopes .
_ chaotic phase synchronization in scale - free networks of bursting neurons _ , phys .
e * 76 * , 016218 ( 2007 ) .
c. a. s. batista , s. r. lopes , r. l. viana , and a. m. batista , _ delayed feedback control of bursting synchronization in a scale - free neuronal network _ , neural networks * 23 * , 114 - 124 ( 2010 ) .
a. s. pikovsky , m. g. rosenblum , g. osipov , and j. kurths .
_ phase synchronization of chaotic oscillators by external driving_. physica d * 104 * , 219 - 238 ( 1997 ) .
y. kuramoto , _ chemical oscillations , waves , and turbulence _
( dover publications , new york , 1984 ) .
j. a. acebrn , l. l. bonilla , c. j. p. vicente , f. ritort , and r. spigler , _ the kuramoto model : a simple paradigm for synchronization phenomena _ , reviews of modern physics * 77 * , 137 - 185 ( 2005 ) s. strogatz , _ from kuramoto to crawford : exploring the onset of synchronization in populations of coupled oscillators _ , physica d * 143 * , 120 ( 2000 ) .
g. restrepo , e. ott , b. r. hunt , _ onset of synchronization in large networks of coupled oscillators _ , phys .
e , * 71 * , 036151 ( 2005 ) .
d. l. wang . _
the time dimension for scene analysis_. ieee transcation on neural networks , * 16 * , 1401 ( 2005 ) .
b. ibarz , j. m. casado , and m. a. f. sanjun .
_ map - based models in neuronal dynamics_. phys
* 501 * , 1 ( 2011 ) . c. koch and i. segev ,
_ methods in neuronal modeling : from ions to networks _
( the mit press , cambridge , mass .
1999 ) a. l. hodgkin and a. f. huxley , _ a quantitative description of membrane current and its application to conduction and excitation in nerve _
, j. physiology * 117 * , 500 - 544 ( 1952 ) .
r. e. plant and m. kim .
_ mathematical descriptions of a bursting pacemaker neuron by a modification of hodgkin - huxley equations
. _ biophys .
j. * 16 * , 227 - 244 ( 1976 ) .
p. r. shorten and d. j. n. wall .
_ a hodgkin
huxley model exhibiting bursting oscillations _ , bulletin of mathematical biology , * 62 * , 695 - 715 ( 2000 )
. h. a. braun , m. t. huber , m. dewald , k. schfer , and k. voigt._computer simulations of neuronal signal transduction : the role of nonlinear dynamics and noise_. int .
j. bifurcat .
chaos * 8 * , 881 ( 1998 ) .
w. braun , b. eckhardt , h. a. braun , and m. huber .
_ phase - space structure of a thermoreceptor_. phys .
e * 62 * , 6352 - 6360 ( 2000 ) .
h. a. braun , m. t. huber , n. anthes , k. voigt , a. neiman , x. pei , and f. moss .
_ noise - induced impulse pattern modifications at different dynamical period - one situations in a computer model of temperature encoding_. biosystems , * 62 * , 99 - 112 ( 2001 ) h. a. braun , m. dewald , k. schfer , k. voigt , x. pei , k. dolan , and f. moss .
_ low - dimensional dynamics in sensory biology 2 : facial cold receptors of the rat .
_ j. comput
. neurosci .
* 7 * , 17 - 32 ( 1999 ) .
k. schfer , h. a. braun , r. c. peters , and f. bretschneider .
_ periodic firing pattern in afferent afferent discharges from electroreceptor organs of catfish_. pflgers arch .
j. physiol .
* 429 * , 378 - 385 ( 1995 ) .
u. feudel , a. neiman , x. pei , w. wojtenek , h. braun , m. huber , and f. moss .
_ homoclinic bifurcation in a hodgkin - huxley model of thermally sensitive neurons_. chaos * 10 * , 231 - 239 ( 2000 ) .
n. f. rulkov .
_ regularization of synchronized chaotic bursts_. phys .
* 86 * , 183 - 186 ( 2001 ) .
d. s. andres , f. gomez , f. s. ferrari , d. cerquetti , m. merello , r. l. viana , and r. stoop , _ diffusive coupling modulates neuronal coding in the parkinsonian globus pallidus _ , submitted .
n. f. rulkov .
_ modeling of spiking - bursting neural behavior using two - dimensional map_. phys .
e * 65 * , 041922 ( 2002 ) .
s. nadarajah and s. kotz , _ a truncated cauchy distribution _ , internat .
techn . * 37 * , 605 - 608 ( 2006 ) .
n. f. rulkov , i. timofeev , and m. bazhenov .
_ oscillations in large - scale cortical networks : map - based model_. j. comp .
* 17 * , 203 - 223 ( 2004 ) . c. a. s. batista , a. m. batista , j. c. a.
de pontes , s. r. lopes , and r. l. viana .
_ bursting synchronization in scale - free networks . _ chaos , solit . &
fractals , * 41 * , 22202225 ( 2009 ) .
h. yu , j. wang , b. deng , x. wei , _ et al .
_ _ chaotic phase synchronization in small - world networks of bursting neurons
. _ chaos , * 21 * , 013127 ( 2011 ) .
r. l. viana , a. m. batista , c. a. s. batista , j. c. a. de pontes , f. a. dos s. silva , and s. r. lopes .
_ bursting synchronization in networks with long - range coupling mediated by a diffusing chemical substance . _ comm .
simulat . , * 17 * , 2924 - 2942 ( 2012 ) .
m. v. ivanchenko , g. v. osipov , v. d. shalfeev , j. kurths . _ phase synchronization in ensembles of bursting oscillators_. phys
lett . * 93 * , 134101 ( 2004 ) . m. g. rosenblum and a. s. pikovsky .
_ delayed feedback control of collective synchrony : an approach to suppression of pathological brain rhythms .
e * 70 * , 041904 ( 2004 ) .
a. arenas , a. fernndez , and s. gmez , _ a complex network approach to the determination of functional groups in the neural system of caenorhabditis elegans .
notes comput .
* 515 * , 918 ( 2008 ) . c. c. hilgetag and m. kaiser , organization and function of complex cortical networks ,
_ in _ `` lectures in supercomputational neuroscience ( dynamics in complex brain networks ) '' , eds .
p. b. graben , c. zhou , m. thiel , and j. kurths ( springer , berlin - heidelberg - new york , 2008 ) .
g. zamora - lpez , c. zhou , and j. kurths .
_ exploring brain function from anatomical connectivity_. frontiers in neuroscience * 5 * , 83 ( 2011 ) .
j. w. scannell and m. p. young .
_ the connectional organization of neural systems in the cat cerebral - cortex_. curr . biol . * 3 * , 191 - 200 ( 1993 )
. j. w. scannell , c. blakemore , and m. p. young .
_ analysis of connectivity in the cat cerebral - cortex_. j. neurosci .
* 15 * , 1463 ( 1995 ) .
r. albert and a .-
_ statistical mechanics of complex networks_. rev .
phys . * 74 * , 47 ( 2002 ) .
d. j. watts , _ small worlds _ ( princeton university press , princeton , 2000 ) .
p. erds and a. rnyi .
_ on the evolution of random graphs_. bulletin of the international statistical institute , * 38 * , 343 - 347 ( 1959 ) .
m. e. j. newman .
_ models of the small world_. j. stat .
* 101 * , 819 - 841 ( 2000 ) .
d. j. watts and s. h. strogatz . _
collective dynamics of small - world networks_. nature , * 393 * , 440 - 442 ( 1998 ) .
m. e. j. newman and d. j. watts .
_ renormalization group analysis of the small - world network model_. phys .
lett . a * 263 * , 341 - 346 ( 1999 ) .
j. c. a. de pontes , r. l. viana , s. r. lopes , c. a. s. batista , and a. m. batista , _ bursting synchronization in non - locally coupled maps _ , physica a , * 387 * , 4417 - 4428 ( 2008 ) .
d. r. chialvo , _ critical brain networks _ , physica a * 340 * , 756 - 765 ( 2004 ) .
o. sporns , d. r. chialvo , m. kaiser , and c. c. hilgetag , _ organization , development and function of complex brain networks _ , trends cogn .
* 8 * , 418 - 425 ( 2004 ) . v. m. eguiluz , d. r. chialvo , g. a. cecchi , m. baliki , and a. v. apkarian , _ scale - free brain functional network _ , phys .
lett . * 94 * , 018102 ( 2005 ) .
a. l. barabasi and r. albert , _ emergence of scaling in random networks .
_ science * 286 * , 509512 ( 1998 ) .
m. breakspear , s. heitmann , and a. daffertshofer , _ generative models of cortical oscillations : neurobiological implications of the kuramoto model . _
frontiers in human neuroscience , * 4 * , 114 ( 2010 ) .
l. timms and l. q. english .
_ synchronization in phase - coupled kuramoto oscillator networks with axonal delay and synaptic plasticity _ , physical review e , * 89 * , 032906 ( 2014 ) .
a. arenas , a. daz - guilera , j. kurths , y. moreno , and c. zhou , _ synchronization in complex networks_. phys . rep . * 469 * , 93153 ( 2008 ) .
m. jalili , _ failure tolerance of spike phase synchronization in coupled neural networks .
_ chaos * 21 * , 033126 ( 2011 ) . c. a. s. batista , e. l. lameu , a. m. batista , s. r. lopes , t. pereira , g. zamora - lpez , j. kurths , and r. l. viana .
_ phase synchronization of bursting neurons in clustered small - world networks .
e * 86 * , 016211 ( 2012 ) .
a. nordenfelt , j. used , and miguel a.f .
sanjun , _ bursting frequency versus phase synchronization in time - delayed neuron networks . _
e * 87 * , 052903 ( 2013 ) .
t. de l. prado , s. r. lopes , c. a. s. batista , j. kurths , and r. l. viana , _ synchronization of bursting hodgkin - huxley - type neurons in clustered networks _ , phys .
e * 90 * , 032818 ( 2014 ) .