synchronization is a phenomenon observed in many different contexts , ranging from epidemics spreading @xcite to neuro - sciences @xcite .
the analysis of simple models pursued in the last two decades has clarified several aspects of synchronization ( for a comprehensive review on the subject see ref .
@xcite ) .    particularly interesting both from a theoretical and an applied point of view is synchronization of chaotic systems .
this is a phenomenon known since many years  @xcite with important applications , e.g. , in secure communications of digital signals , laser dynamics etc .
a simple but not trivial framework to study chaotic synchronization is represented by discrete - time dynamical systems . within this context
, the essence of the phenomenon can be captured already by considering two coupled identical maps : @xmath3 where @xmath4 is a chaotic map of the unit interval and @xmath5 the coupling constant . for coupling larger than a critical value @xmath0 a transition from a desynchronized state to a completely synchronized one along a chaotic orbit @xmath6
is observed , i.e. @xmath7 for @xmath8 .
the coupling @xmath0 can be predicted within a linear analysis as the value of @xmath5 for which the transverse lyapunov exponent ( tle ) @xmath1 , ruling the instabilities transversal to the line @xmath9 , vanishes . for model ( [ eq:2map_1 ] )
@xmath10 @xcite , and consequently @xmath11 being @xmath12 the lyapunov exponent ( le ) of the single ( uncoupled ) map .
while in low dimensional systems the basic mechanisms ruling synchronization , in the absence of strong nonlinear effects , have been well established since long time  @xcite , the studies devoted to synchronization in extended systems are still limited to few specific examples  @xcite . in the latter case ,
two main frameworks have been considered .
on one hand , the investigation of mutual synchronization between two replicas of a spatially extended system has been pursued for systems coupled either via local interaction @xcite or via spatio - temporal noise @xcite . in both cases , when the amplitude of the coupling ( noise ) overcomes a certain threshold , synchronization is observed : after a transient , the two systems follow the same spatio - temporal chaotic ( stochastic ) orbit .    on the other hand , another commonly observed situation corresponds to self synchronization of elements belonging to the same system .
this occurs for instance in large collections of coupled elements such as populations of neurons @xcite , josephson junctions @xcite or cardiac pacemaker cells @xcite . in these systems
the interaction among the elements can range from nearest neighbors to globally coupled .
therefore , a model able to encompass both these limiting cases is certainly of interest . a good candidate is represented by the following coupled map lattice with power - law interactions  @xcite @xmath13 which reduces to globally coupled maps ( gcm s ) @xcite for @xmath14 , and to standard coupled map lattices ( cml s ) with nearest neighbor coupling @xcite in the limit @xmath15 . in equation ( [ eq : cml ] ) , @xmath16 and @xmath17 are the discrete temporal and spatial indexes , @xmath18 is the lattice size ( @xmath19 ) , @xmath20 the state variable , @xmath21 $ ] measures the strength of the coupling and @xmath22 its power - law decay .
since the sum extends up to @xmath23 the model is well defined only for odd @xmath18-values , and @xmath24 is a normalization factor .
periodic boundary conditions are assumed .
the synchronization between replicas of the same spatial system has been the subject of a more systematic study than self synchronization and is now rather well understood  @xcite .
two different mechanisms of mutual synchronization have been identified according to the predominance of linear versus nonlinear effects .
nonlinearly dominated transitions are observed when the local dynamics is ruled by a discontinuous ( e.g the bernoulli map ) or `` almost discontinuous '' maps ( e.g. maps possessing very high values of the first derivatives ) . in these cases a linear analysis is no more sufficient to fully characterize the transition , because the instabilities associated with finite amplitude perturbations may desynchronize the system . in both cases
the transition to the synchronized state is typically of the second order .
however , depending on the linear or non - linear nature of the prevailing mechanism two different universality classes characterize the transition itself . for continuous maps ( e.g. logistic maps ) critical exponents
associated with the multiplicative noise ( mn ) universality class are usually found  @xcite , while for ( almost ) discontinuous maps the transition belongs to the directed percolation ( dp ) class  @xcite .
the aim of the present work is to clarify the effect of strong nonlinearities in the synchronization phenomenon , in particular in the case of self synchronization .    as a first point
, we show that strong nonlinearities may induce nontrivial effects also in low dimensional systems such as two coupled maps . in particular
, we find that even though the critical coupling @xmath0 can be always predicted within the linear framework , the transient dynamics preceding the synchronization can be strongly affected by finite scale instabilities . to properly characterize the latter we introduce a new indicator , the finite size transverse lyapunov exponent ( fstle ) , which generalizes the concept of finite size lyapunov exponent ( fsle )  @xcite to the transverse dynamics .
the fstle extends the definition of tle to finite perturbation and is therefore able to discern linearly from nonlinearly dominated transitions .
depending on the behavior of the fstle at finite amplitudes two different class of maps are singled out maps characterized by a decreasing fstle at any finite scale , and _ class ii _ maps that present a peak in the fstle for some finite amplitude value .
moreover , in the proximity of the transition the shape of the probability density function ( pdf ) , @xmath25 , of the synchronization times @xmath26 depends on the dynamics at finite scales and on the multifractal properties of the map . for class
i maps two cases can be identified . for maps , such as the symmetric tent map and the logistic map at the crisis , which are characterized by the vanishing of fluctuations of the finite time le at long times
, the pdf exhibits a fast falloff at large @xmath26 s .
in particular , for the symmetrical tent map we provide an analytical expression for @xmath25 . for maps , such as the skew tent map , exhibiting _
modulational intermittency _ @xcite ( which is related to the persistence of fluctuations of the finite time le also in the long time limit ) , @xmath25 becomes an inverse gaussian distribution @xcite originating by the diffusive motion of the ( logarithm of the ) perturbation in transverse space . for class
ii maps the pdf s display an exponential tail at long times which , differently from the previous case , is due to nonlinear effects .    as far as spatially extended systems are concerned , anteneodo et al .
@xcite have performed a linear stability analysis of model ( [ eq : cml ] ) obtaining analytically the critical coupling @xmath0 for the transition .
this prediction works perfectly for class i maps , as found for the coupled logistic maps  @xcite .
however , as shown in the present paper , linear analysis may fail in class ii maps .
in particular , for large system sizes ( @xmath27 ) an exponential divergence of the synchronization times with @xmath18 can be observed even for negative tle , as a result in the limit @xmath28 the transition takes place at a critical coupling @xmath29 .
finite scale instabilities are the key elements for observing this _ nonlinear synchronization transition_. similarly to the case of diffusively coupled map lattices it is possible to link the synchronization transition of cml with power law coupling with non - equilibrium phase transitions . indeed as shown in ref .
@xcite synchronization of replicas of cml with short range coupling belong either to the mn or to the dp universality classes . here
we shall discuss the connection of self synchronization of model ( [ eq : cml ] ) with long range spreading processes  @xcite .
the material is organized as follows .
[ sec:2 ] is devoted to synchronization of two coupled maps . in sect .
[ sec:3 ] coupled maps with power law coupling are examined . finally in sec .
[ sec:4 ] we briefly summarize the reported results .
to analyze the synchronization of two coupled maps ( [ eq:2map_1 ] ) it is useful to introduce the following variables : @xmath30 in terms of which eq .
( [ eq:2map_1 ] ) can be rewritten as @xmath31\nonumber\\ w^{t+1 } & = & \left ( \frac{1}{2}-\varepsilon \right )   \left[f(u^t+w^t ) - f(u^t - w^t ) \right ] \label{eq : diff}\,.\end{aligned}\ ] ] it is easily checked that the synchronized solution ( corresponding to @xmath32 and @xmath33 ) is an admissible solution of ( [ eq : diff ] ) for any value of the coupling @xmath5 .
the stability of such a solution can be studied by considering the linear dynamics of an infinitesimal perturbation @xmath34 of the synchronized state .
this evolves according to the linearized equation : @xmath35 where @xmath36 .
clearly the stability of the state @xmath32 is controlled by the sign of the transverse ( or conditional ) lyapunov exponent @xmath37 where @xmath12 is the lyapunov exponent of the single map @xmath4 .
therefore , by requiring @xmath38 one obtains the critical coupling @xmath0 ( [ eq : epsc ] ) above which the synchronized state is linearly stable .
notice that for @xmath39 , @xmath1 coincides with the second lyapunov exponent @xmath40 of the coupled system ( [ eq:2map_1 ] ) .
the synchronization transition has been mainly studied for continuous maps , here denoted as _ class i _ maps , in particular for the logistic map @xmath41 at the crisis @xcite and the skew tent map @xcite @xmath42 it should be noticed that inside this class of maps one has to distinguish between two situations according to the behavior of the fluctuations of the finite time lyapunov exponents in the long time limit .
the vanishing of such fluctuations is observed for the symmetric tent map and the logistic one at the crisis , while their persistence characterizes the skew tent map that represents the generic case @xcite .
recent studies however pointed out that discontinuous or `` almost '' discontinuous maps ( i.e. maps with @xmath43 in some point of the definition interval ) , here termed _ class ii _ maps , may give rise to non trivial interesting phenomena .
for instance , the synchronization transition for two replicas of class ii cmls is not driven , as usual , by infinitesimal perturbations , but by the instability ( spreading ) of finite amplitude perturbations @xcite .
this nonlinear synchronization is linked to the so - called _ stable chaos _
@xcite : a phenomenon characterized by information propagation even in the absence of chaos  @xcite .
moreover , as we shall show in the following also the self synchronization of a single cml can be strongly affected by nonlinearities .
starting from random initial conditions . from top : bernoulli map with @xmath44 , logistic map at the crisis and the skew tent map with @xmath45 .
for all maps , the coupling has been chosen in order to ensure @xmath46 .
the dashed straight lines display the @xmath47 decay . ]
representative examples of class ii maps are the generalized bernoulli map @xmath48 that is discontinuous and its continuous version @xmath49 with @xmath50 , @xmath51 , @xmath52 .
as already mentioned , the peculiarity of such maps is that they are not only unstable with respect to infinitesimal perturbations but also to finite ones ( for ( [ eq : bsc ] ) this is true for sufficiently small @xmath53 values ) .
the main differences between synchronization transitions ruled by linear mechanisms with respect to those driven by finite amplitude instabilities can be captured already by considering the evolution of @xmath54 .
the typical behavior of this quantity for @xmath55 is reported in fig .
[ fig : wt ] for the bernoulli , the logistic and the skew tent map . for the logistic map , apart from ( short time ) fluctuations due to the multiplier variations , an average linear decrease of @xmath56 with slope given by @xmath1 is observed . for the bernoulli map , despite the tle is negative and does not fluctuate in time , @xmath54 exhibits instantaneous jumps to @xmath57 values .
these jumps , which are due to finite amplitude instabilities , are more probable for large @xmath54 amplitudes .
completely analogous behaviors have been found for the map ( [ eq : bsc ] ) , where the multipliers do fluctuate .    as shown in fig.[fig : wt ] for the skew tent map with @xmath58 , @xmath59 displays an intermittent evolution during the transient preceding the synchronization , that is reflected in a diffusive motion of @xmath60 . as a consequence long transients
are observed close to the transition .
this phenomenon , known as _ modulational intermittency _
@xcite is induced by fluctuations of the finite time tle that do not average out in the long time limit and as a consequence can be explained by linear analysis .    for the sake of completeness
, we mention that long synchronization transients have been reported also for nonlinearly coupled expanding maps  @xcite . in this model resurgences of @xmath54 during the transient are due to the presence of an invariant chaotic repelling set .
we stress that for class ii maps the transient is absolutely non chaotic i.e. their transverse dynamics in tangent space is contracting at any time .      to quantitatively characterize the different _
transverse_-space dynamics ( i.e. the evolution of @xmath61 shown in fig .
[ fig : wt ] ) for maps of class i and ii , let us now introduce the finite size transverse lyapunov exponent .
the fstle , @xmath62 , generalizes the concept of transverse lyapunov exponent to finite value of the perturbation @xmath63 .    following aurell et al .
@xcite ( see also ref .
@xcite ) we have defined the fstle as follows .
we introduce a set of thresholds @xmath64 with @xmath65 , since on average a transverse expanding ( resp .
contracting ) dynamics is expected in the desynchronized ( resp .
synchronized ) regime , we choose @xmath66 ( e.g. @xmath67 ) for @xmath68 and @xmath69 ( e.g. @xmath70 ) for @xmath8 .
first , starting from a random initial conditions we wait for @xmath6 to relax onto its attractor .
then we initialize a second variable @xmath71 as @xmath72 by choosing an initial perturbation @xmath73 ( resp .
@xmath74 ) if @xmath75 ( resp .
if @xmath76 ) , and let @xmath6 and @xmath77 evolve according to eq .
( [ eq:2map_1 ] ) .
care should , of course , be taken to maintain @xmath71 within the interval of definition of the map . during the evolution
we record the time , @xmath78 , needed for @xmath54 to pass for the first time from one threshold @xmath79 to the following one @xmath80 and also the value of @xmath81 at the moment of the passage .
when the last threshold , @xmath82 , is reached the system is reinitialized with the above described procedure and this is repeated @xmath83 times . the fstle is thus defined as @xmath84 where the average @xmath85\rangle$ ] is performed over the set of @xmath83 different initial conditions . in the limit
@xmath86 the fstle converges to @xmath1 .
notice that due to its definition , the fstle can not measure at the same time expansion and contraction rates , i.e. we have limited the analysis only to consecutive contractions ( resp .
expansions ) in the synchronized ( resp .
desynchronized ) regime .
this implies that the sign of @xmath62 is always negative or positive in accordance with the investigated situation .    in fig .
[ fig : fsle ] @xmath62 versus @xmath87 is shown for different maps in the desynchronized regime . as expected , in all cases , for very small values of @xmath54
the tle is recovered .
however , at larger values of @xmath54 maps of the class ii display an increase of the growth rate , i.e. , @xmath88 for some finite @xmath87 value , while for maps of class i the fstle is monotonically decreasing with @xmath87 .
@xmath62 provides a quantitative measure of the strength of the nonlinear effects that , in maps of class ii , may in principle overwhelm the linear mechanisms as pointed out in ref .
@xcite .
vs @xmath87 in the desynchronized regime for the bernoulli map with @xmath89 ( empty circles ) , its continuous version ( empty boxes ) with @xmath90 , the logistic map ( filled circles ) and the skew tent map ( filled triangles ) .
the coupling constants have been chosen in such a way that in all systems the transverse le is equal to @xmath91 ( solid line ) . ]    in fig .
[ fig : firstarrival ] we report the behavior of @xmath62 in the synchronized regime . again , while for very small @xmath87 the ( negative ) tle is recovered , at larger @xmath87 values the bernoulli - like maps display an increase which is due to the jumps @xmath92 of the transverse perturbation .
vs @xmath87 in the synchronized regime for the bernoulli map at @xmath89 ( empty circles ) , the logistic map ( filled boxes ) .
the inset show the logistic map ( filled boxes ) and the skew tent map ( empty boxes ) .
the straight lines correspond to the transverse le , here @xmath93 for all systems . ]
a closer inspection of the small @xmath87 regime reveals that for the logistic map @xmath94 only for very small scales @xmath95 , while for the skew tent map the asymptotic @xmath94 in never exactly reached .
the latter convergence difficulties are probably due to the intermittent behavior along the transverse direction , in the proximity of the transition @xcite . in the regime
@xmath96 the finite time effects are less important and the asymptotic behavior is recovered for all the investigated maps .
the results presented in this subsection clearly indicate that finite scale instabilities can prevail on infinitesimal ones only for strongly nonlinear maps .      notwithstanding the observed differences in the fstle
, we found that the critical coupling for synchronization to occur is always given by ( [ eq : epsc ] ) for both classes .
this means that , at least for the case of two coupled maps , nonlinear effects do never modify the critical coupling .
it is natural to wonder whether other observables related to the synchronization transition could be influenced by the presence of finite amplitude instabilities . as we shall show in this section ,
this is the case for the synchronization times statistics .
let us define the synchronization time @xmath26 , as the shortest time needed to @xmath97 for decreasing below a certain threshold @xmath98 , with the further requirement that @xmath99 for a sufficiently long successive time @xmath100 .
however , if the threshold value is small enough ( e.g. @xmath101 ) @xmath26 essentially coincides with the first arrival time to the considered @xmath98 .
the quantities of interest are the pdf s , @xmath25 , and their moments ( we shall focus mainly on the first moment ) .      in fig .
[ fig : pdftent ] we show @xmath25 for the the logistic and the symmetric tent map .
the main feature is the presence of an exponential tail at short times and the faster than exponential falloff at long times .
this means that synchronization times @xmath102 are not observed . as a first result we derive , by following ref .
@xcite , an approximate analytical expression for @xmath25 , which is exact for the symmetric tent map .
let us define @xmath103 , and consider its linearized evolution , from eqs .
( [ eq : tang_diff ] ) and ( [ eq : lambdat ] ) one obtains @xmath104 formally the above equation can be applied only in the true linear regime , i.e. when @xmath105 , so that it is not appropriate in the early stages of the evolution . however , as shown in fig .
[ fig : wt ] , for maps like the logistic one , usually after a reasonably short transient , the linear regime sets in and the use of eq .
( [ eq : zet1 ] ) is justified for the successive evolution . from eq .
( [ eq : zet1 ] ) it is clear that @xmath25 is related to the pdf @xmath106 of the variable @xmath107 and to that of the local multipliers @xmath108 .
the formal solution of eq .
( [ eq : zet1 ] ) up to time @xmath109 can be written as follows : @xmath110 where @xmath111 .
for sufficiently large @xmath109 , large deviation theory tells us that the pdf of @xmath112 takes the form @xmath113 , being @xmath114 the cramer function @xcite , which is convex and has its minimum value at @xmath115 .
it is now clear the distinction between the generic case in which @xmath114 does not collapse onto a @xmath53-function , and the non - generic one in which it does , as for the symmetric tent map and the logistic one at the ulam point . in the former case
the dynamics of @xmath116 becomes a biased brownian motion , with an average drift given by @xmath1 .
for the sake of simplicity , let us start from the tent map ( eq .  ( [ eq : skt ] ) for @xmath117 ) , for which @xmath118 . in this case
@xmath26 is simply given by : @xmath119 and @xmath25 can be directly related to @xmath106 . since , in the proximity of the transition , @xmath120 assumes the form ( analytically derived in ref .
@xcite ) : @xmath121 it can be easily obtained the related distribution @xmath122 which perfectly agrees with the numerical results ( fig .
[ fig : pdftent]a ) .     as a function of the synchronization time @xmath26 : ( a ) tent map for @xmath123 , the dashed line is the analytical expression ( [ eq : p_tau ] ) ; ( b ) logistic map for @xmath124
, the solid line refers to the expression ( [ eq : p_tau_log ] ) with @xmath125 and @xmath126 . in all cases
the pdf have been obtained by averaging over @xmath127 different initial conditions and by choosing @xmath128 . ]    unfortunately the ( short time ) multiplier statistics of the logistic map is nontrivial , impeding a straightforward derivation of @xmath25 .
however , we numerically observed that @xmath129 ( @xmath130 ) with @xmath131 almost constant in a neighborhood of the transition ( namely , @xmath132 for @xmath133 ) .
this implies that @xmath134 which is in fairly good agreement with the numerically evaluated pdf ( fig .
[ fig : pdftent]b ) .
notice that @xmath135 is not simply given by @xmath136 , and a fitting procedure is needed .
we found @xmath137 , with @xmath138 and 0.65 for @xmath139 and @xmath140 , respectively .
the above reported results allow us also to predict the scaling of the average synchronization time @xmath141 with @xmath142 in the proximity of the synchronization transition .
in particular , for the symmetric tent map the following expression can be derived @xmath143 where @xcite @xmath144 being @xmath145 the euler - mascheroni constant . since @xmath146 , approximatively we have @xmath147 where @xmath148 ( being @xmath149 ) . note that ( [ eq : tmed ] ) is in perfect agreement with the numerical results for the tent map ( see fig .
[ fig : tmedtent ] ) . for the logistic map ,
a similar dependence on @xmath150 , namely @xmath151 , has been found . the interesting point in ( [ eq : tmed ] )
is the logarithmic correction to the scaling @xmath152 , which would have not been predicted by a naive guess based on ( [ eq : tau ] ) .
vs @xmath153 for the tent map ( empty boxes ) and the logistic one ( empty circles ) .
the continuous line is the prediction ( [ eq : tmed ] ) with @xmath128 , the dotted one is obtained by a best fit of the form @xmath154 , with @xmath155 and @xmath156 . the data of the logistic map have been shifted by a factor @xmath157 for plotting purposes . ]    we now briefly discuss the generic situation ( see chap .  13 of ref .
@xcite for more details ) . as mentioned above
, @xmath158 performs a biased brownian motion , so that computing @xmath25 amounts to evaluate the distribution of the first passage times to a threshold of a wiener process with drift that , as a standard result of stochastic process theory , is given by an inverse gaussian density @xcite that for long times @xmath26 becomes @xmath159 where the quadratic dependence from the tle comes from the diffusive dynamics of the perturbation .      for class ii maps
the situation is completely different .
for both the continuous bernoulli map ( fig .
[ fig : pdfbsc]a ) and the discontinuous bernoulli shift ( fig .
[ fig : pdfbsc]b ) we observe that for @xmath8 the pdf s are characterized by an exponential tail at large @xmath26 , similar to poissonian distributions .
moreover , we also observe that the pdf s , once rescaled as @xmath160 , and reported as a function of @xmath161 ( where @xmath162 ) collapse onto a common curve in the proximity of the transition .
these results are particularly striking in the case of the bernoulli map for which naively one would have expected a complete similarity with the tent map , since its multiplier @xmath163 is also constant . here
the presence of strong nonlinear effects makes the transient time statistics to be dramatically different from ( [ eq : p_tau ] ) . indeed , as seen in fig .
[ fig : wt ] , @xmath164 is subject to noticeably nonlinear amplifications induced by the ( almost ) discontinuities in the map .
therefore , the expression ( [ eq : zet1 ] ) is no more appropriate to describe the dynamics of @xmath165 , and its full nonlinear dynamics has to be taken into account .
the latter is characterized by the fact that , with a finite probability , @xmath54 can jump to @xmath92 values at any time during the transient preceding the synchronization , even if the transient is not chaotic .
this idea can be better clarified and its consequences on the @xmath25 better appreciated by considering a simple stochastic model for the dynamics of the transverse variable at @xmath8 .
this model was originally proposed in ref .
@xcite , and reads as @xmath166 where @xmath167 and @xmath168 is a shorthand notation for with probability .
the underlying idea is very simple : the transverse perturbation @xmath61 is usually contracted , but with probability proportional to its amplitude can be re - expanded to @xmath169 values , in agreement with the numerical observations .
the jumps occur when @xmath6 and @xmath77 are close but located at the opposite sides of the ( almost ) discontinuity .     as a function of @xmath170 ( where @xmath162 ) for the continuous bernoulli shift map ( [ eq : bsc ] ) with @xmath171 and @xmath172 at different distances , @xmath173 , from the critical coupling .
note the fairly good collapse , indicating a common decay as @xmath174 ( dashed line ) .
the collapse is however only approximate near the peak .
( b ) same as ( a ) for the bernoulli map with @xmath175 at different distances from the critical coupling .
the dashed line gives @xmath174 . ]    for this simple model , it is possible to derive the corresponding probability density @xmath25 , which displays the same peculiar features of the pdfs reported in figs .
[ fig : pdfbsc]a and [ fig : pdfbsc]b . in particular , the long time tail .
first , notice that the probability that an initial perturbation @xmath176 is never amplified to @xmath92 up to time @xmath16 is given by @xmath177 \to\ !
\bar q(w_0 ) \!=\!\lim_{t
\to \infty } \exp\left [ -\sum_{k=1}^\infty \frac{w_0^k { \rm e}^{-|\lambda_t| k}}{k(1- { \rm e}^{-|\lambda_t| k})}\right ] .
\label{eq : p0}\end{aligned}\ ] ] interestingly , for @xmath178 the quantity @xmath179 is always positive and strictly smaller than @xmath180 for any perturbation @xmath181 , i.e. the probability that a perturbation of amplitude @xmath182 can be amplified is finite at any time .
the pdf of the times @xmath26 needed to observe @xmath183 can be factorized as @xmath184 where @xmath185 is the probability to receive a kick at time @xmath186 , and @xmath187 is the probability of not being amplified for the following @xmath188 steps , where @xmath188 is given by @xmath189 as shown in ref.@xcite @xmath190 , and therefore at large @xmath26 s eq .
( [ eq : ptau_teo ] ) can be rewritten as @xmath191 } \label{eq : ptau_teo_2}\ ] ] that confirms the poissonian character of the pdf @xmath25 for the maps of class ii . by properly normalizing eq .
( [ eq : ptau_teo_2 ] ) , one obtains @xmath192 , that is in good agreement with the numerical results for the bernoulli map ( in particular , we considered @xmath175 and @xmath193 ) .
moreover , in the interval @xmath194 $ ] we found that the decay rate of the pdf is directly proportional to the tle , i.e @xmath195 , this should be contrasted with the quadratic dependence found in the case of class i maps , see eq .
( [ eq : inv_gaussian ] ) .
we conclude the investigation of two coupled maps by studying the nature of the transition .
for this purpose let us introduce the order parameter @xmath196 defined as ( see also ref .
@xcite ) : @xmath197 i.e. for each value of the coupling @xmath5 , one considers @xmath109 different random initial conditions and each one is iterated for a transient @xmath198 after which the time average of @xmath54 over a time lapse @xmath199 is considered and further averaged over all the initial conditions .
the value of the coupling giving the synchronization transition is then implicitly defined as @xmath200 . in principle
, @xmath2 may differ from @xmath0 defined by ( [ eq : epsc ] ) , since this expression is valid only within a linear approximation formalism .
vs @xmath5 for the bernoulli map with @xmath201 and its continuous version ( [ eq : bsc ] ) for different values of @xmath53 . for the bernoulli map one
has that according to ( [ eq : epsc ] ) @xmath202 .
the order parameter is computed according to ( [ eq : op ] ) averaging over @xmath203 different random initial conditions , discarding the first @xmath127 iteration and averaging over the following @xmath204 iterations .
@xmath196 has been estimated with an @xmath5-resolution @xmath205 .
the inset shows an enlargement of the critical region done with a coupling resolution @xmath206 , the same number of initial conditions and a longer integration transient @xmath207 iterations .
( b ) the same for the logistic map at the crisis . here
( [ eq : epsc ] ) predicts @xmath208 .
the blow up of the critical region ( in the inset ) shows that the transition is continuous .
note also that for @xmath209 we observe @xmath210.,title="fig : " ]   vs @xmath5 for the bernoulli map with @xmath201 and its continuous version ( [ eq : bsc ] ) for different values of @xmath53 .
for the bernoulli map one has that according to ( [ eq : epsc ] ) @xmath202 .
the order parameter is computed according to ( [ eq : op ] ) averaging over @xmath203 different random initial conditions , discarding the first @xmath127 iteration and averaging over the following @xmath204 iterations .
@xmath196 has been estimated with an @xmath5-resolution @xmath205 .
the inset shows an enlargement of the critical region done with a coupling resolution @xmath206 , the same number of initial conditions and a longer integration transient @xmath207 iterations .
( b ) the same for the logistic map at the crisis . here
( [ eq : epsc ] ) predicts @xmath208 .
the blow up of the critical region ( in the inset ) shows that the transition is continuous .
note also that for @xmath209 we observe @xmath210.,title="fig : " ]    in fig .
[ fig : op ] we show the behavior of the order parameter ( [ eq : op ] ) as a function of @xmath5 for both classes of systems .
two observations are in order .
first , the transition to the synchronized state always occurs at the critical coupling defined by ( [ eq : epsc ] ) ( i.e. @xmath211 ) .
second , the transition is continuous for the logistic map ( simulations show that also the tent map , which belongs to the same class , has a continuous transition ) and discontinuous for the bernoulli shift .
this confirms the results reported in ref .
however , since the transition for the continuous version of the bernoulli shift map ( [ eq : bsc ] ) is steep but clearly continuous ( see the inset of fig.[fig : op]a ) , this does not seem to be a generic property for all the maps of class ii .
the behavior of the bernoulli shift map can be easily understood by considering the map controlling the dynamics of @xmath61 , i.e. @xmath212 being @xmath213 $ ] and @xmath2141/r:1]$ ] . from the above expression
, it is evident that the attractor in the @xmath215 plane has always a finite width along the transversal direction , unless @xmath216 .
this explains the discontinuity observed at @xmath0 . in the continuous map ( [ eq : bsc ] )
the attractor in the plane @xmath215 has a `` transverse '' width that decreases continuously to zero for @xmath217 . in conclusion ,
the discontinuity in the transition is a pathology of the bernoulli map , and possibly of other discontinuous maps .
the order parameter ( [ eq : op ] ) can be seen as the time average , once discarded an initial transient , of the following quantity @xmath218 where @xmath219 indicates the average over many different initial conditions .
once @xmath220 is initialized @xmath92 , depending on whether @xmath5 is smaller or larger than @xmath0 two different asymptotic behavior are observed . in the desynchronized regime , obviously @xmath221 goes to a finite value @xmath222 , while above the synchronization transition @xmath223 with a decay law determined by the nature of the considered maps . for class
i maps , the decay is ruled by the tle : @xmath224 with @xmath167 . for class ii
maps @xmath221 shows an initial exponential decay @xmath225 ( with @xmath226 ) , due to the effect of the resurgences averaged over many different initial conditions , followed by a final linear decay @xmath227 .
the nonlinear rate @xmath228 can be simply related to the exponential tail of the pdf of the first arrival times , by assuming ( as indeed observed ) the following decay @xmath229}$ ] .
since @xmath221 is the average amplitude value of the @xmath61 at time @xmath16 , one has @xmath230 the above result tells us that @xmath231 ( if the integral does not diverge , i.e if @xmath232 ) .
this is confirmed by numerical checks .
in this section we analyze high dimensional systems , namely the power - law coupled maps defined in eq .
( [ eq : cml ] ) . in particular , for class ii maps we shall show that there exist situations where , due to the nonlinearities , the synchronization time diverges exponentially with the number of maps . as a consequence , the transition takes place at a ( nonlinear ) critical coupling @xmath2 larger than the linear value @xmath233 .
a similar phenomenon has been observed for gcm s with nonlinear coupling , where the synchronization time divergence is due to a chaotic transient  @xcite .
chaotic transients , diverging exponentially with the number of coupled elements , have been reported also for spatially extended reaction - diffusion systems  @xcite and for diluted networks of spiking neurons  @xcite . however , the emphasis of our work is on non - chaotic transients similar to the stable - chaos phenomenon  @xcite .
let us start by reviewing the linear theory developed in ref .
@xcite , which is able to account for the synchronization transition of class i maps .      when nonlinear effects are not sufficiently strong , excluding the pathological cases of chaotic transients , the critical coupling for observing the synchronization transition can be predicted by computing the lyapunov spectrum . in particular , since above the synchronization transition the transverse lyapunov exponent coincides with the second lyapunov exponent , @xmath40 , it suffices to evaluate the dependence of the latter on @xmath5 .    in order to compute the lyapunov spectrum of the model ( [ eq : cml ] ) , it is necessary to consider the tangent space evolution : @xmath234 in ref .
@xcite it has been shown that the full lyapunov spectrum can be easily obtained for the bernoulli map : @xmath235 where @xmath236 in this notation the maximal lyapunov exponent is given by @xmath237 and the second lyapunov exponent is @xmath40 .    for generic maps
it is still possible to obtain the lyapunov spectrum in an analytical way , but only in the synchronized state , by substituting @xmath238 with the maximal le of the single uncoupled map @xmath12 in ( [ eq : spectrum ] ) .    from a linear analysis point of view one
expects that the synchronization transition has to occur when @xmath239 .
therefore , the following expression for the critical line in the @xmath240-plane can be derived : @xmath241 as shown in ref .
@xcite this prediction is well verified for the logistic map , which is here used as a benchmark for our codes .
let us conclude this short review by noticing that in the limit @xmath242 synchronization can be achieved only for @xmath243  @xcite .
therefore we shall limit our analysis to this interval .
let us now study the critical line for class ii maps , exemplified by the bernoulli shift map and its continuous version .
first of all we introduce the main observables and the numerical method employed to determine it .    a meaningful order parameter for the transition is represented by time average of the following mean field quantity : @xmath244 this can be operatively defined as follows : firstly the system ( [ eq : cml ] ) is randomly initialized and iterated for a transient time @xmath245 proportional to the system size @xmath18 , then the time average of @xmath246 is computed over a time window @xmath199 , i.e. @xmath247 .
finally the state of the system is defined as synchronized if @xmath248 , being @xmath98 a sufficiently small value ( @xmath249 is usually enough ) . the coupling value corresponding to the synchronization transition
is then obtained by using a bisection method : chosen two coupling values across the transition line one corresponding to a desynchronized state ( @xmath250 ) and the other to a synchronized case ( @xmath251 ) a third value is selected as @xmath252 .
if at this new coupling value the system synchronizes ( resp .
not synchronizes ) @xmath253 is identified with the new @xmath251 ( resp .
@xmath250 ) .
the procedure is then repeated until @xmath254 ( in our simulations @xmath255 ) .
finally the critical coupling is defined as @xmath256 .
the algorithm , tested on the logistic map , was able to recover ( [ eq : epsc_cml ] ) with the required accuracy .
in general one has that , once fixed @xmath22 , @xmath257 will be a function of @xmath245 for a given @xmath18 .     at @xmath258 .
the solid lines are the analytically estimated `` linear '' transition values eq .
( [ eq : epsc_cml ] ) , while the symbols refer to the numerically obtained values .
we used @xmath259 , @xmath260 , @xmath261 and @xmath262 . ]    in refs .
@xcite it has been employed as a synchronization indicator the time average of the following mean - field quantity : @xmath263 we have verified that this indicator gives results completely analogous to @xmath264 in all the considered cases .    in ref .
@xcite the authors have reported for logistic coupled maps a very good agreement between @xmath265 given by ( [ eq : epsc_cml ] ) and @xmath257 estimated via @xmath266 . however , this is not the case for the bernoulli map , as shown in fig .  [
fig : sync.bern ] . in this case
( depending on the slope of the map @xmath267 , on @xmath22 and on the chain length @xmath18 ) strong disagreements between the linear transition line given by @xmath265 and the numerically obtained values are observed .
these disagreements are typical of class ii maps . since the nonlinear effects locally desynchronize the system even if @xmath268 , in general @xmath269 .     as a function of @xmath18 for different @xmath245 in @xmath270 for the continuous bernoulli shift map with @xmath201 and @xmath271 at @xmath272 ( a ) and @xmath273 ( b ) .
the same for the bernoulli shift with @xmath201 for @xmath272 ( c ) and for @xmath273 ( d ) .
the continuous solid lines are the theoretical values obtained in the linear analysis framework ( [ eq : epsc_cml ] ) . ]    in obtaining the data reported in fig .
[ fig : sync.bern ] , the cpu time restrictions forced us to employ a large but somehow limited transient time ( namely , @xmath259 ) .
therefore , we checked for the dependence of the results on @xmath245 for fixed @xmath18 .
in particular , we measured @xmath257 for two @xmath22-values only ( namely , @xmath272 and 0.8 ) for several chain lengths and transient times .
this analysis has been performed for the logistic map at the crisis , for the bernoulli shift map with @xmath201 and for its continuous version .
the results are reported in fig .
[ fig : transients ] . for the logistic map we found that , for any value of @xmath18 and @xmath22 , a relatively small value of @xmath274 was sufficient to observe a clear convergence of @xmath257 to the linear value @xmath265 .
this is also the case of the bernoulli map and its continuous version for @xmath272 ( though as one can see in fig .
[ fig : transients]c the discontinuous map is characterized by a slightly slower convergence than its continuous version ) . on the contrary for @xmath273 with @xmath275
we were unable to observe synchronization even for @xmath276 , while at @xmath277 the logistic map had already converged to @xmath265 . summarizing for @xmath201
we found that for @xmath278 the bernoulli shift map and its continuous version always converge toward the linear critical value , while for @xmath279 at increasing @xmath18 the critical coupling @xmath257 becomes more and more independent of the transient time @xmath245 ( see the right of fig .
[ fig : transients ] for @xmath280 ) .
however , for a smaller value of @xmath267 and for large @xmath18 ( we investigated the case @xmath281 and @xmath282 ) we observed @xmath29 for any value of @xmath283 $ ] . in particular , we stress that also in the globally coupled case ( corresponding to @xmath284 ) we have clear evidences that @xmath285 .
this is due to the fact that for @xmath286 nonlinear finite amplitude instabilities becomes more and more predominant with respect to the mechanism of linear stabilization ( see  @xcite for more details ) .
on the basis of the previous results , it is natural to conjecture that in the limit @xmath28 the non - linear transition is well defined , i.e. that the limit @xmath287 exists and is typically larger than @xmath288 .
note that in the above expression the order of the two limits is crucial , we expect that performing at fixed @xmath18 the limit @xmath289 we should always observe a convergence to @xmath0 . however , as shown in the following , the times to reach synchronization may diverge exponentially fast with @xmath18 in the region just above @xmath0 making rapidly infeasible this limit .      as for the case of two coupled maps ,
we study now the synchronization time statistics .
in particular , we consider the system above the linear transition line @xmath8 and measure the first passage times @xmath26 needed for @xmath246 ( [ indicator ] ) to decrease below a given threshold @xmath98 . in this way
we determine the corresponding pdf and the related momenta .    in the high dimensional case ,
it is fundamental to analyze the dependence of the synchronization time on the the system size @xmath18 .
however , as reported in eq .
( [ eq : epsc_cml ] ) , the critical value @xmath0 itself depends on @xmath18 . therefore to perform a meaningful comparison of systems of different sizes we considered situations characterized by the same linear behavior in the transverse space , i.e. having the same value of @xmath40
this means that for each length @xmath18 we chose the coupling strength according to @xmath290 further tests performed at different lengths but at a fixed distance @xmath291 from the critical line give essentially the same results .
vs @xmath26 for different system sizes ( see label ) at fixed @xmath292 for different maps : ( a ) logistic map at the crisis with @xmath273 ( for @xmath272 we obtained qualitatively similar results ) ; skew tent map ( [ eq : skt ] ) with @xmath45 ( b ) ; bernoulli shift map with @xmath201 for @xmath272 ( c ) and @xmath273 ( d ) . in
the latter case sizes larger than @xmath293 were not drawn for the sake of clarity of the plot .
the pdf s have been obtained in all cases by considering @xmath294 different initial conditions and with @xmath128 . ]    as shown in fig .
[ fig : pdf_log]a , for the logistic map ( but these results can be extended to all the maps belonging to class i , see fig .
[ fig : pdf_log]b ) the pdf s of the synchronization times obtained at constant @xmath40 display a very weak ( almost absent ) dependence on the system size , and are qualitatively similar to that found for two coupled maps ( compare fig .
[ fig : pdf_log]a with fig .
[ fig : pdftent ] ) .
moreover , measurements of the average synchronization times @xmath295 , done by fixing the `` distance '' @xmath296 from the critical line , exhibit a clear tendency to saturate for increasing @xmath18 ( fig .
[ fig : tau_log ] ) .
therefore , in the limit @xmath28 the synchronization time will not diverge .
, as a function of the system size @xmath18 for the logistic map at the crisis with @xmath297 estimated at a fixed distance from the linear threshold @xmath298 .
data have been averaged over @xmath204 initial conditions with @xmath128 . ]    for the maps of class ii the situation is different . here , as shown in fig .
[ fig : pdf_log]c , for values of @xmath22 sufficiently small @xmath25 is weakly dependent on the system size as found for the logistic map . on the other hand , by considering more local couplings ( e.g. @xmath273 in fig .
[ fig : pdf_log]d ) , the tail of the pdf becomes more and more pronounced as the system size increases .
results obtained by fixing the distance from the critical coupling instead of the value of the tle display qualitatively similar features .
this picture is further confirmed by examining @xmath295 as a function of @xmath18 for fixed @xmath40 ( i.e. by choosing @xmath5 according to ( [ eps_l2 ] ) ) .
as one can see in fig .
[ fig : tau_bernou ] , @xmath295 displays a dramatic dependence on @xmath18 .
in particular , at @xmath299 ( see the inset of fig .  [
fig : tau_bernou ] ) @xmath300 grows exponentially with @xmath18 , while a power - like scaling is observable for @xmath301 , at least for the chain lengths we could reach .
deeper inside the `` linear '' synchronization region , i.e. for @xmath302 , we observed a saturation of @xmath303 with the system size .
it is worth stressing here that for @xmath272 and @xmath201 , where no appreciable discrepancies between @xmath2 and @xmath0 have been observed , @xmath295 saturates for large @xmath18 .
nevertheless the corresponding pdf exhibits the usual tail at long times characteristic of maps of class ii but with weak dependence on @xmath18 ( see fig .  [
fig : pdf_log]c ) .
notice also that , in principle , the transition from the exponential dependence to the saturation can be used to estimate @xmath2 ( see ref .
@xcite ) , however for the present model a systematic study of this aspect is infeasible due to the required computational resources     as a function of @xmath18 for the bernoulli map with @xmath201 and @xmath273 for various values of the second lyapunov exponent ( see label ) . the inset shows the case of @xmath292 in linear scale .
data have been averaged over @xmath294 initial conditions and @xmath128 . ]
we thus found clear indications that a nonlinear synchronization transition can be observed for class ii maps when finite size nonlinear instabilities are sufficiently strong to overcome the effects of stabilization associated with the local linear dynamics and to the spatial coupling .
this suggests that nonlinear effects tend to decouple the single units of the chain .
this decoupling can be indeed related , by means of the following simple argument , to the observed exponential divergence with @xmath18 of the synchronization times at @xmath304 .
synchronization happens when the difference @xmath305 is contracted for ever in each site until @xmath306 decreases below threshold .
the probability that for a single site the initial difference @xmath307 will be never amplified to @xmath92 is given by @xmath308 ( see eq .
( [ eq : p0 ] ) ) , therefore by assuming that each sites is completely decoupled from its neighbors the probability of contraction for the whole chain is given by @xmath309}$ ] ( for simplicity we set @xmath310 without loss of generality ) . and the probability that the system settles onto this `` contracting '' state after @xmath311 steps is @xmath312 this quantity represents a good approximation of @xmath25 .
therefore one expects a poissonian decay for the pdf of the synchronization times with an associated average time @xmath313}$ ] exponentially diverging with the chain length , as indeed observed .
let us now characterize the synchronization transition in the spatially extended model ( [ eq : cml ] ) within the framework of non - equilibrium phase - transitions .
a similar parallel was recently established in a series of works concerning synchronization of two replicas of cml s with nearest neighbor coupling  @xcite . in these studies
it has been shown that the synchronization transition is continuous and belongs to the multiplicative noise ( resp . directed percolation )
universality class depending on the linear ( resp .
nonlinear ) nature of the prevailing mechanisms . in both cases
one observes a transition from an active phase ( characterized by @xmath314 ) to a unique absorbing state ( identified by @xmath315 ) .
we indicated with @xmath316 the order parameter , however for sake of clarity it should be said that for two replicas of a cml this corresponds to @xmath317 and is not given by expression ( [ indicator ] ) used in the present paper for characterizing self synchronization .
the rationale for using the same symbol is that the two definitions embody essentially the same information : @xmath318 is a measure of the density of non synchronized ( active ) sites .    in the proximity of the transition ( if continuous ) the scaling behavior of the saturated density @xmath319 of active sites
is given by @xmath320 while at the critical point @xmath321 the density of active sites scales as @xmath322 and the average synchronization time diverges as @xmath323 we denoted as @xmath324 the critical coupling to avoid at this stage distinctions between linearly and nonlinearly driven synchronization . since continuous non - equilibrium phase transitions are typically characterized by three independent critical exponents , once ( @xmath325 ) are known all the other scaling exponents can be derived  @xcite .     for @xmath326 ( filled circles ) @xmath327 ( empty boxes ) and @xmath328 ( crosses ) .
the solid straight line displays the power law decay @xmath329 .
data refer to the logistic map with @xmath330 and @xmath331 and are obtained after averaging over @xmath332 initial conditions . in the inset the average value of the order parameter @xmath333 as a function of @xmath334
is reported .
the straight line shows the scaling behavior @xmath335 .
data refer to the logistic map with @xmath272 and @xmath331 the average is performed over @xmath336 different initial conditions . ]
let us now analyze cml s with power - law coupling ( [ eq : cml ] ) .
we start with class i maps , for which the transition is completely characterized by the linear dynamics , i.e. @xmath337 . in particular
, we consider the logistic maps at the ulam point for three different @xmath22-values ( namely , @xmath338 and @xmath339 ) .
the first observation is that the critical properties of the model seem to be independent of @xmath22 : we measured exactly the same critical exponents for the all @xmath22 values , namely @xmath340 , @xmath341 ( see fig .
[ fig : scaling ] ) and @xmath342 ( see fig .
[ fig : tau_log ] ) .
these exponents coincide with the mean - field exponents reported in ref .
@xcite for a model of anomalous directed percolation ( dp ) .
the model differently from standard dp considers the spreading of epidemics in the case of long - range infection , namely in one spatial dimension the probability distribution for a site to be infected at distance @xmath267 decays as @xmath343 .
the critical behavior of this model can be obtained by considering a langevin equation with power - law decaying spatial coupling for the coarse grained density of infected sites . in ref .
@xcite , the authors found that the critical exponents vary continuously with @xmath344 , but below a critical value ( @xmath345 in 1d ) the exponents coincide with the corresponding mean - field results , namely @xmath346 and @xmath347 . in the limit
@xmath348 the latter exponents are identical to the ones we have found for the logistic coupled maps . in the following we shall give an argument to explain these similarities and differences .    as shown in refs .
@xcite there is a deep connection between the synchronization problem for two replicas of a diffusively coupled cml and nonequilibrium phase transitions . indeed there
the difference field @xmath349 can be mapped onto the density of infected ( active ) sites and an appropriate langevin equation describing the evolution of @xmath350 can be derived in proximity of the transition  @xcite . in our case of cmls with power - law coupling
we expect that the corresponding langevin equation for the spatio - temporal coarse grained defect density @xmath351 should contain a long - range interaction with a power law spatial coupling decaying with an exponent @xmath22
. therefore , in the proximity of the transition , it seems natural to map the dynamics of ( [ eq : cml ] ) onto the model studied in  @xcite once the exponent @xmath22 is identified with @xmath352 . however , in ( [ eq : cml ] ) the coupling is rescaled by the factor @xmath353 to avoid divergences in the limit @xmath354 . the rescaling @xmath355 amounts to consider effective interaction on large scales of the type @xmath356 independently of @xmath22 .
this may explain the fact that the exponents characterizing the transition of model ( [ eq : cml ] ) coincide with the mean - field values reported in  @xcite for @xmath357 .
a similar rescaling was performed in ref .
@xcite to show , for a chain of power - law coupled rotators , that all the equilibrium properties of the system coincide for @xmath358 , once suitably scaled .
let us also remark that since we are dealing with a mean - field case the nature of the noise entering in the langevin equation , which distinguishes directed percolation from multiplicative noise , is irrelevant @xcite .
this leads us to conclude that , for continuous synchronization transitions , we would not expect any differences in the measured exponents between linearly and nonlinearly driven transitions .     for the bernoulli map with @xmath201 , @xmath273 , @xmath331 , and different coupling strength above and below @xmath2 ( see the labels ) .
average is done over @xmath359 initial conditions , tests on shorter times and more initial conditions confirm the reliability of the statistics . ]    in examining the critical properties associated to coupled class ii maps , we confronted with hard numerical difficulties related to the fact that @xmath360 can not be measured with the required precision .
moreover , due to the presence of exponentially long transients in the proximity of the nonlinear transition , simulations of model ( [ eq : cml ] ) with class ii maps become extremely time consuming . for these reasons we were unable to find conclusive results .
however , within these limitations , our analysis ( performed up to @xmath331 and up to integration times @xmath361 ) suggests that for class ii maps ( both the bernoulli map and its continuous version ( [ eq : bsc ] ) with @xmath201 were considered ) the transition may be discontinuous as suggested from the data shown in fig .
[ fig : disco ] . in the figure the temporal evolution of the order parameter @xmath246 averaged over many initial conditions
is reported in proximity of @xmath2 and no indication of critical behaviour is observable .
this indicates that in class ii cmls : short range interactions give rise to dp - like continuous transitions , while power - law coupled maps exhibit discontinuous transitions .
we should mention that the introduction of long - range interactions in dp model can similarly modify the nature of the transition from continuous to discontinuous , as shown in a recent paper @xcite .
the authors have introduced a generalized directed percolation model , in which the activation rate of a site at the border of an inactive island of length @xmath362 is given by @xmath363 . in particular , in ref .  @xcite it has been shown that the transition is continuous for @xmath364 with exponent coinciding with dp , and discontinuous for @xmath365 .
in the present paper we have examined the influence of strong nonlinearities in the transverse synchronization transition of low and high dimensional chaotic systems . in particular , for two coupled maps we have shown that finite amplitude nonlinear instabilities can give rise to long transients , preceding the synchronization , even when the dynamics is transversally stable at any instant .
this should be contrasted with the results found for continuous maps where long transients may happen only as a result of transient chaotic or intermittent transverse dynamics . in high dimensional systems with long range ( power - law ) interactions strong nonlinearities
may invalidate the linear criterion to locate the critical coupling . in this case
the transition occurs , due to nonlinear mechanisms , at a larger coupling value .
the nonlinear transition is characterized by the emergence of transients diverging exponentially with the system size even above the linear critical coupling .
the origin of this transition is closely related to the stable chaos phenomenon .
the synchronization phenomena both in linearly and nonlinearly driven systems have been compared with models for long ranged contact processes .
we found that , for linearly driven systems , the transition is continuous and the critical exponents are given by a mean field prediction . for
nonlinearly driven systems , though the results are not conclusive , evidence of a discontinuous transition have been found .
we are grateful to w.  just , a.  politi and a.  pikovsky for useful discussions and remarks and to f.  ginelli also for a careful reading of this manuscript .
partial support from the italian firb contract n. rbne01cw3m_001 is acknowledged .
99 d. he and l. stone , proc r soc lond b biol sci .
* 270 * ( 2003 ) 1519 .
p. n. steinmetz , a. roy , p. j. fitzgerald , s. s. hsiao , k. o. johnson , and e. niebur , nature * 404 * ( 2000 ) 187 .
a. pikovsky , m. rosenblum , and j. kurths , _ synchronization a universal concept in nonlinear sciences _ , ( cambridge university press , 2001 ) l.m .
pecora and t.l .
carroll , phys .
* 64 * ( 1990 ) 821 .
kuznetsov and a. s. pikovsky , radiophys .
, * 32 * ( 1989 ) 49 .
a. s. pikovsky and p. grassberger , j. phys .
a : math . gen .
* 24 * ( 1991 ) 4587 .
w. just , physica d * 81 * ( 1995 ) 317 .
l.  baroni , r.  livi , and a.  torcini , phys .
e * 63 * ( 2001 ) 036226 . v. ahlers and a. s. pikovsky , phys .
* 88 * ( 2002 ) 254101 . c. anteneodo , s.e .
de s. pinto , a. m. batista , and r. l. viana , phys .
e , * 68 * ( 2003 ) 045202(r ) ; c. anteneodo , a.m. batista , and r. l. viana , phys . lett .
a , * 326 * ( 2004 ) 227 . m. dharmala , v.k .
jirsa , and m. ding , phys .
* 92 * ( 2004 ) 028101 .
k. wiesenfeld , p. colet and s.h .
strogatz , phys .
* 76 * ( 1996 ) 404 . c. peskin ,
_ mathematical aspects of heart physiology
_ ( courant institute of mathematical sciences , new york university , new york , 1975 ) . g. paladin and a. vulpiani , j. phys . * a 25 * ( 1994 ) 4911 .
a.  torcini and s.  lepri , phys .
e , * 55 * ( 1997 ) r3805 .
k. kaneko , physica d , * 34 * ( 1989 ) 1 .
i.  waller and r.  kapral , phys .
a * 30 * , 2047 ( 1984 ) ; k. kaneko , prog .
theor . phys .
* 72 * ( 1984 ) 980 . v. ahlers , _ scaling and synchronization in deterministic and stochastic nonlinear dynamical systems _ ,
phd thesis ( potsdam , 2001 ) f. ginelli et al .
e * 68 * ( 2003 ) 065102(r ) .
m. a. munz and r. pastor - satorras , phys .
* 90 * ( 2003 ) 204101 .
m. a. munz , in  advances in condensed matter and statistical mechanics " , eds .
e. korutcheva _
( nova science publishers , 2004 , new york ) .
h. hinrichsen , adv .
* 49 * ( 2000 ) 815 .
e. aurell , g. boffetta , a. crisanti , g. paladin and a. vulpiani , phys .
* 77 * ( 1996 ) 1262 .
g. boffetta , m. cencini , m. falcioni and a. vulpiani phys . rep .
* 356 * ( 2002 ) 367 .
m. cencini and a. torcini phys .
e * 63 * ( 2001 ) 056201 .
w.  feller , _ an introduction to probability theory and its applications _
( wiley , new york , 1974 ) ; h.c .
tuckwell , _ introduction to theoretical neurobiology - vol .
2 - nonlinear and stochastic theories _ ( cambridge university press , cambridge , 1988 ) . f. ginelli , h. hinrichsen , r. livi , d. mukamel , and a. politi , phys .
e * 71 * ( 2005 ) 026121 .
a.  politi , r.  livi , g.l .
oppo and r.  kapral , europhys .
* 22 * ( 1993 ) 571 .
a.  torcini , p.  grassberger and a.  politi , j. phys .
* a27 * ( 1995 ) 4533 .
a. politi and a. torcini , europhys .
lett . , * 28 * ( 1994 ) 545 .
g.  paladin and a.  vulpiani , phys . rep . * 156 * ( 1987 ) 147 .
h. jeffreys and b.s .
jeffreys , in _ methods of mathematical physics _ , ( cambridge , england , cambridge university press,1988 ) pp .
470 - 472 .
f. ginelli , r. livi , and a. politi , j. phys .
a : math . gen .
* 35 * ( 2002 ) 499 .
a.  lipowski and m.  droz , `` synchronization and partial synchronization of linear maps '' , cond - mat/0312067 r. wackerbauer and k. showalter , phys .
* 91 * ( 2003 ) 174103 .
a. zumdieck , m. timme , t. geisel , and f. wolf , phys .
* 93 * ( 2004 ) 244103 .
f. ginelli , r. livi , a. politi , and a. torcini phys .
e * 67 * ( 2003 ) 046217 .
h. hinrichsen and m. howard , eur .
j. b * 7 * ( 1999 ) 635 .
f. tamarit and c. anteneodo , phys .
* 84 * ( 2000 ) 208 .