Synchronization in time-varying random networks with vanishing connectivity

A sufficiently connected topology linking the constituent units of a complex system is usually seen as a prerequisite for the emergence of collective phenomena such as synchronization. We present a random network of heterogeneous phase oscillators in which the links mediating the interactions are constantly rearranged with a characteristic timescale and, possibly, an extremely low instantaneous connectivity. We show that with strong coupling and sufficiently fast rewiring the network reaches partial synchronization even in the vanishing connectivity limit. In particular, we provide an approximate analytical argument, based on the comparison between the different characteristic timescales of our system in the low connectivity regime, which is able to predict the transition to synchronization threshold with satisfactory precision beyond the formal fast rewiring limit. We interpret our results as a qualitative mechanism for emergence of consensus in social communities. In particular, our result suggest that groups of individuals are capable of aligning their opinions under extremely sparse exchanges of views, which is reminiscent of fast communications that take place in the modern social media. Our results may also be relevant to characterize the onset of collective behavior in engineered systems of mobile units with limited wireless capabilities.


I. INTRODUCTION
The emergence of collective phenomena in complex systems is related to the interplay between interaction topology and local dynamics [1][2][3][4] .Stationary connections can lead to coherent dynamical patterns, typically studied in the framework of network theory, with the local dynamics taking places on individual nodes and interactions modelled as links.In the context of complex networks, conditions of minimal connectivity are know for enabling the emergence of collective dynamics 1,5 .A prominent example is synchronization in networks of oscillators [6][7][8][9] , where the connectivity thresholds for a wide range of different network topologies have been determined in great detail 10,11 .However, many complex systems, and particularly social and engineered ones, may not maintain a constant connectivity, but rather yield a topology defined by a time-dependent connectivity matrix A t ij .Examples range from animal groups 12,13 and time-dependent plasticity in neural networks 14,15 to robot swarms 16 , human social networks 17 and communication networks of moving units 18 .
Synchronization in time-varying networks received considerable attention in the control and nonlinear dynamics literature [19][20][21][22][23] .Yet, these efforts almost exclu-sively concentrated on systems composed of homogeneous units, largely relying on standard linear stability analysis.Introducing quenched disorder, that is, considering systems composed by many heterogeneous oscillators is however more challenging, especially for a finite number of units, as the stability of the partially synchronized state cannot be typically treated by simple linear stability analysis 24 .In this situation, averaging theorems 25 may not be trivially applicable, so that different approaches may be needed.
A first step towards the study of time-varying networks of heterogeneous oscillators is provided by Ref. 26 , which focused on two populations of oscillators switching between two fixed topologies at a given frequency.Interestingly, analysis of this "blinking" network revealed that high-frequency switching may induce synchronization, even when the two individual topologies can only sustain an incoherent phase.
While these findings provide a first hint that results for time-varying networks of homogeneous units can be extended to the heterogeneous case, here we wish to take a step further and study a time-varying network of heterogeneous oscillators, where individual nodes interact randomly (and possibly quite seldomly) in both time and oscillator space In particular, here we ask under which conditions macroscopic synchronization may emerge in Erdös-Rény networks with random rewiring and arbitrarily small instantaneous connectivity.We thus consider N heterogeneous agents interacting randomly, with a bidirectional and typically sparse connectivity matrix in a regime of strong coupling.We will see that our system is characterized by three different timescales: interacting agents quickly converge towards a common state on a short local syncronization timescale τ LS , while each agent may randomly rewire all his connections with a typical rewiring timescale T .When two connected agents are separated, their internal states diverge, with yet another local de-sinchronization timescale τ LD which depends on their heterogeneity and it is typically larger than τ LS .One can interpret this setup as a crude model of social interactions, where individuals interact in time with different subsets of their common social network.When interacting, and despite their intrinsic differences, they tend to quickly converge towards a common opinion, but when separated, their differences take over again and their opinion diverges.
Model parameters allow to control the separation between these characteristic times, which enables a detailed study of the emergence of synchronization in relation to the interplay between different timescales.In the following, we show that -provided the links are rewired frequently enough -dynamics can permanently achieve a partially synchronized state, even when the instantaneous connectivity is far smaller than what is needed to synchronize stationary networks.In particular, via numerical simulations and approximate analytical arguments of a concrete model, we show that for a sufficiently strong coupling and a sufficiently fast rewiring, our system reaches and maintains a macroscopic (partially) synchronized state, even in the limit of vanishing connectivity: it is the high frequency blinking of links that prevents the system from relaxing into an incoherent state as would happen with stationary topologies.This paper is organized as follows: in Section II we define a precise model for our time-varying network and sketch its synchronization phase diagram through direct numerical simulations.In Section III we focus on the low connectivity regime.Analysing the characteristic timescales of the system and invoking an averaging theorem in the limit T → 0, we provide an approximate expression for the synchronization threshold which compares favorably with numerical estimates.In Section IV we first discuss higher connectivities, where the instantaneous network topology is characterized by a system spanning giant component, and then argue that the transition to synchronization belongs to the standard Kuramoto class in the entire phase diagram.Conclusions are finally drawn in Section V.

II. KURAMOTO MODEL ON TIME-VARYING NETWORKS II.A. Model definition
We first introduce our model.Let us consider a network of N Kuramoto oscillators, where the state of the i-th node is represented by a phase variable ϕ i ∈ [0, 2π].Each oscillator is characterized by a quenched natural frequency ω i , drawn from a zero-mean Gaussian distribution with standard deviation σ.Oscillators interact with each others according to a time-varying adjacency matrix A t ij , with m t i = j A t ij being the instantaneous degree of node i.For simplicity, we chose the adjacency matrix to be symmetric and with binary values A t ij = 0, 1, leaving other cases for future studies.Hence, the dynamics of the oscillators obey the following equation: where epsilon quantifies the strength of the coupling 29 .Obviously, when no edges at all insist on node i we have simply φi = ω i .The dynamics of A t ij (T ) is determined as follows.At each moment, the adjacency matrix corresponds to a random, or Erdös-Rény (ER) network, defined by the vertex number N and the linking probability p 1 .The random rewiring of edges is then modelled as a Poissonian process, with each individual node rewiring synchronously all its incident edges with probability rate 1/T , with T being the typical rewiring time.In the instantaneous rewiring process of vertex i, all edges incident on i are first deleted; all the potential links of vertex i are then considered, and new edges i − j are created with probability p.It is well known that the topological properties of ER networks are essentially determined by the mean degree connectivity m = (N − 1)p, so that in the following we find convenient to define q = p N ≈ m and adopt q ≈ m as the relevant connectivity parameter.
Note finally that Eq. ( 1) is invariant under the following rescaling: (with α ∈ R + ), provided also the rewiring time is rescaled accordingly, T = αT .Due to this invariance, it is easy to show that the dynamics is actually controlled by the two dimensionless quantities T /σ and ε/σ and by the connectivity q.
One can interpret this setup as a crude model for several natural/social phenomena.Consider for example social interactions, where individuals interact in time with different subsets of their common social network, according to a certain frequency of personal encounters/interactions.When interacting, despite intrinsic differences of their opinions (i.e.different quenched natural frequencies), individuals tend to quickly con- verge towards a common opinion.However, when separated, their differences take over again, and their opinions slowly diverge.The spectrum of natural frequencies ω i can thus represent the range of "unperturbed" opinions of a population, while in Eq. ( 1) interactions with other persons (nodes) leads to the effective frequencies ω ef f i (t) = φi , representing the actual opinion of the agents.
It is well known that the Kuramoto model with stationary network connectivity, either globally connected 24 or with other sufficiently connected topologies 10 , displays a synchronized solution for large enough couplings .In this synchronized state, a macroscopic fraction of oscillators share a common effective frequency, reaching macroscopic consensus in our point of view.The degree of synchronization can be evaluated through the standard instantaneous Kuramoto order parameter which is finite for synchronized states and tends to zero as 1/ √ N in the absence of macroscopic synchronization.In the following, we will typically consider its average over time and disorder (i.e.different natural frequencies realization), ∆ = R t ,ω , and make use of its different finite size scaling behavior to better estimate the transition between (partial) synchronization and disorder.

II.B. Direct numerical simulations
In this work, numerical simulation are performed using a standard 4th order Runge-Kutta integrator of step dt.After each Runge-Kutta time-step, each vertex may undergo a rewiring event (as defined above) with Poissonian probability For small values of q the system is strongly dependent on the value of the rewiring time T and the phase diagram shows a clear transition from partial synchronization to disorder as T is increased beyond a critical value Tc(q).At larger q values, the transition approaches a vertical asymptote, roughly located at q = q = 1.66(6) (dashed black line).For q > q the dynamics achieves partial synchronisation regardless of the value of T .(b) Zoomed view of panel (a) in the range q ∈ [0.5, 1.7].
We use a time-step of at most dt = 10 −2 .When investigating fast network dynamics however, we are forced to adopt time-steps smaller than the network rearrangement timescale T , that is dt ≈ T /10.In order to illustrate the behavior of our time-varying network dynamics in a strong coupling regime, = 8, σ = 1, we begin presenting numerical simulations of the dynamics (1) for a network of N = 10 4 elements and a mean connectivity q = 0.8.
As it is shown in Fig. 1, no synchronization emerges when the rewiring is sufficiently slow (T ≈ 6.3 in this example).As the rewiring time is lowered past a synchronization threshold, we observe macroscopic synchronization with an increasing order parameter R(t).This shows that sufficiently fast network rewiring can overcome the effects of low network connectivity, inducing partial synchronization on the network.
We next want to characterize with more details the parameter space (q, T ).We do that by repeating the above computation for a grid (lattice) of different values of q and T .For each of these values we calculate ∆ = R t ,ω by averaging over 10 random realizations of the quenched natural frequencies ω i and different random initial phases.Time averages are performed over the stationary part of R(t), after a proper initial transient has been discarded.The results are shown via colorplot in Fig. 2a (the lighter the color, the larger the value of ∆).We first observe that for a sufficiently large connectivity, q > q, the system always reaches macroscopic synchronization, regardless of the rewiring time T .Analysis of the averaged order parameter ∆ in the large T limit, as reported in more details in Section IV.A, suggests q = 1.66 (6).Here, we just wish to point out that q is clearly larger than q = 1, the threshold for the emergence of a giant connected component in ER graphs 30 .In this regime, synchronization is indeed to be expected in the strong coupling limit, due to sufficient interactions among the oscillators.
For smaller values of q, on the other hand, where no large components characterize the instantaneous network topologies, sufficiently fast rewiring is needed to achieve synchronization, at least for q > 0.5.A transition line T c (q) separating partial synchrony from incoherence (i.e. the violet border between the dark and the bright zone) in the plane (q, T ) can be roughly identified from this colour plot.Indeed, a closer look at the phase diagram, as reported in Fig. 2b, suggests that the transition line T c (q) separating partial synchrony from incoherence in the plane (q, T ) is initially characterized by a linear behaviour.For larger connectivity values, on the other hand, T c (q) grows faster than linear, finally diverging as a vertical asymptote is approached at q = q.Note however that the transition line is characterized by a non zero intercept at q = q 0 ≈ 0.5 with the T = 0 axis.Thus, for smaller connectivity values (q 0.5), no synchronization is possible for the coupling /σ = 8, no matter how fast is the rewiring.
In the following section we will proceed to better characterize the transition to synchrony in the low average connectivity region q < 1 by means of approximate analytical arguments and detailed numerical simulations.

III. SYNCHRONIZATION FOR LOW AND VANISHING CONNECTIVITY III.A. Characteristic time scales and the onset of synchronization
We next seek to understand the physical mechanism leading to synchronisation in the low connectivity region via switching.For this we need to grasp the three characteristic timescales governing information flow and the dynamics of our system.The first time scale is the local synchronisation time τ LS , related to the synchronisation of a connected pair of oscillators.The second is the the local desynchronisation time τ LD , related to the typical desynchronisation time as the link between two synchronized oscillators is severed.The third one, finally, is the effective rewiring time τ ER , describing the typical time needed for an oscillator to establish a new link after a rewiring event.
We focus on the limit in which τ LS is much smaller than both τ LD and τ ER .In this regime, oscillators couples quickly synchronize when connected by a link, starting to loose their relative synchrony when their mutual link is deleted in a rewiring event.In practice, oscillators tend to loose the information gained when linked with the characteristic timescale τ LD .Two possibilities are then in order for low connectivity.Either a new link is forged by one of these two oscillators with a third node in a time shorter than τ LD , propagating the information it carries from its previous local synchronization to a new node, or no link at all is established before this information is completely lost.We argue that global synchronization will take place when, on average, the information gained by local synchronization events does not get lost but is rather able to propagate through the entire network.This will happen when τ LD τ ER .On the other hand, when τ ER τ LD , no information can propagate through the network, and macroscopic synchronization cannot take place.The transition from the desynchronized to the synchronized regimes will thus take place when Note that a similar argument, based on the characteristic timescales of information transfer, has been previously successfully applied to estimate the transition line separating disordered from collective motion in the well known Vicsek model for flocking [31][32][33] .
We now proceed to estimate the three timescales introduced above.First consider the local synchronization scale τ LS , that is, the time needed by two oscillators i and j sharing a non-directed link to synchronize their effective frequencies.In the low connectivity approximation one can assume for a couple of oscillators m t i = m t j = 1 , i.e. that they are only connected one to each other.Hence, from Eq. ( 1), one immediately gets for their mutual phase difference δϕ = ϕ i − ϕ j the dynamics where δω = ω i − ω j is the difference between their natural frequencies.In the strong coupling regime we are interested into, ε σ and one readily sees that the phase difference converges exponentially fast towards the asymptotic solution δϕ = δω/(2ε) while the two oscillators effective frequencies synchronize with a time scale τ LS ≈ (2ε) −1 .In the following we first assume τ LS T , that is, once a link is established oscillators typically synchronize before being rewired.Once the link is removed in a rewiring event, nodes can be left without any link, so that the phase of previously connected and synchronized oscillators will start to drift away one from each other due to their natural frequencies difference δω, loosing any information regarding their previous mutual synchronization when their phase difference approaches π/2.This allows one to define the typical local desynchronization timescale τ LD such that with being the average natural frequency difference.For Gaussian distributed natural frequency one of course has which finally yields Before proceeding further, one comment is in order about our estimate of the typical local desynchronization timescale.We have computed it as the time required by a typical pair of oscillators to desynchronize.This is of course different from the average of individual couples desynchronization times π/(2δω) , which is dominated by oscillators couples with almost degenerate natural frequencies, δω ≈ 0. These latters, however, characterized by a very large local desinchronization time, are far from being representative of the typical behavior of random oscillators couples.
We finally estimate the effective rewiring timescale τ ER .The probability for an oscillator to be linked to at least one other oscillator is equal to: where P not link is the probability of not having any link at all, that is Substituting back into Eq.( 10) we get with lowest order corrections of order q 2 and q/N .We thus evaluate the effective rewiring time in the low connectivity limit as Summing up, the synchronization condition (5) yields a linear relation between the rewiring time T and the connectivity q, yielding the synchronization line We now compare our predictions with numerical simulations.We determine the synchronization threshold by finite size analysis, comparing the averaged order parameter ∆(N ) for system sizes N 1 = 1000 and N 2 = 2000.In the presence of macroscopic synchronization one expects ∆(N 1 ) ≈ ∆(N 2 ), while in the disordered phase we have An example of our procedure is given in Fig. 3 for q = 0.8, where we have estimated T c = 0.53 (6).
Numerical estimates of the synchronization threshold are reported in Fig. 4 for σ = 1 and different values of the coupling constant ε.They confirm the linear relation between T c and q in the low connectivity regime, predicting the actual slope s = π/(2 √ 2 σ) within numerical accuracy (see inset).However, it is clear that for finite values of the coupling, it is always possible to find sufficiently small values of q such that synchronization cannot be achieved, no matter how small is T .Said differently, the critical line T c has a non-zero intercept q 0 ( ) with the T = 0 axis.Interestingly, the value of q 0 ( ) of the intercept decreases towards zero as ε increases, suggesting that Eq. ( 14) can be fully recovered as → ∞.
This is equivalent to the strong coupling limit under which we have derived Eq. ( 14): By taking the limit → ∞ first, in fact, we assure that the condition T is verified for any non-zero rewiring time T .On the other hand, numerical simulations with a finite coupling constant show that one can always find a sufficiently low connectivity q such that T c (q) τ LS = (2 ) −1 and our approximation breaks down.
In the next section, we will attempt to better understand this regime and the behavior of the intercept q 0 ( ) by means of averaging considerations.

III.B. Average network for very fast rewiring
We now consider the limit of extremely fast rewiring, where T and τ ER ≈ T /q are much smaller than the local synchronization and desynchronization times.In this regime, one expects the instantaneous order parameter to be approximately constant over a timescales τ av min(τ LS , τ LD ), so that Following the argument of Ref. 26 , we may invoke a well known result from Ott and Antonsen 34 to argue that the low dimensional dynamics of the Kuramoto order parameter is essentially controlled by the time-averaged interaction matrix This result, stating that for T → 0 the dynamics of Eq. ( 1) is the same as the one of the time-average network with stationary connectivity, can be essentially seen as a form of the averaging theorem 25 .While the latter typically involves periodic systems, a recent extension to non-periodic systems has been discussed, for instance, in Ref. 35 .
In the limit T → 0 the average in Eq. ( 17) is computed over arbitrarly many rewiring events and we have where is the probability that node i has an active link with node j and exactly k − 1 other links.Note that the binomial factor accounts for all the different configurations in which the k − 1 active link can be chosen out of N − 2 potential ones after the one between i and j has been activated.By recalling that p = q/N , using Eqs.( 18)-( 19) one can find where .
In the limit N 1 we have and therefore, to lowest order in 1/N , so that the average network is characterized by a globally connected topology.Therefore, under our conjecture, in large networks the fast rewiring (T → 0) dynamics (1) can be replaced by the averaged one that is, a globally coupled Kuramoto model with coupling constant which exhibits macroscopic synchronization for J > J c , with the critical point J c depending on the natural frequency distribution.For T → 0, therefore, macroscopic synchronization can only be achieved provided > J c and for connectivities q > q 0 with We conclude that, according to Eq. ( 27), in the strong coupling limit /σ → ∞ synchronization can be achieved for arbitrarily small connectivity q.
Furthermore, we can interpret q 0 as the intercept of the transition line T c (q 0 ) with the T = 0 axis.In particular, for a Gaussian distribution of natural frequencies with unit standard deviation we have J c = 8/π 24 , which allows us to to compare Eq. ( 27) with the intercept values obtained by extrapolating the best linear fit for the transition lines of Fig. 4. Direct comparison (see Fig. 5) shows excellent agreement in the coupling range ∈ [8, 32] we have probed.We can now correct Eq. ( 14) by adding a constant term such that T c (q 0 ) = 0, thus obtaining with for Gaussian distributed natural frequencies.This is exactly the linear formula we plotted in Fig. 4 for unit variance (σ = 1), showing good comparison with the numerical transition values T c (q, ) in the small T regime.As q (and thus T c ) grows larger, however, deviations from the linear behavior are clearly visible.Indeed, as q is increased, the critical line T c (q) bends upwards to meet the vertical asymptote at q.In this regime, contributions from nodes with more than one link at the time becomes relevant, and the simple arguments leading to the linear relation (28) are expected to break down.

IV. BEHAVIOR AT FINITE CONNECTIVITY
For completeness, in this section we briefly discuss the synchronization transition at finite connectivity q.

IV.A. Behaviour for large connectivity
We have already seen that, in order to synchronize for arbitrarily large rewiring times, the connectivity q should be larger than a threshold q > 1, so that the typical emergences of connected components of macroscopic (black), σ 2 = 1/4 (red), σ 2 = 1/2 (green), σ 2 = 1 (blue) and σ 2 = 2 (brown).Bottom: Ratios ∆(N 1 )/∆(N 2 ) as a function of q.The two horizontal dotted lines mark the ratios √ 2 (disordered phase) and 1 (synchronized phase).The colours coding for the variance is the same as in the top panel.The vertical dashed line q = 1 marks the emergence of giant connected components.Data has been averaged over Ω = 10 different realizations.size, taking place for q > 1, is not sufficient for the onset of synchronization.In particular, we have seen that for /σ = 8 we have q = 1.66 (6).We now show numerical evidence that in the large coupling limit, /σ → ∞, we have q = 1 + , that is the onset of synchronization do coincide with the emergence of giant connected components in the graph topology.
Next we show that the phase transition thresholds depend on the distribution of frequencies ω when the coupling ε is fixed, which demonstrates that the critical point is not dependent on the topology (and the percolation threshold for the giant connected component), but only on the dynamics.
In the following, we analyze numerically the synchronization transition at in the region q ≈ 1 through the finite size analysis of the averaged parameter ∆ at large rewiring times T as the ratio /σ is progressively increased.Fixing , we increase σ between 1/ √ 8 and √ 2. In order to evaluate ∆ at large enough rewiring times, in agreement with the scaling relation (2) we fix T such that the dimensionless parameter T σ = 20π 36 , and compare the order parameter at two different system sizes N 1 < N 2 .As already remarked in Sec.III.A, we can distinguish the synchronized from the disordered phase by the ratio ∆(N 1 )/∆(N 2 ).This is, for instance, how in Sec.II.B we have estimated q = 1.66 (6) for /σ = 8 from the data of Fig. 6 (blue symbols).More in general, numerical simulations, reported in Fig. 6, clearly indicates that, as σ is lowered and the strong coupling regime is approached, the synchronization threshold approaches the onset of giant connected components, i.e. q → 1.These results indicate that for finite couplings, in the regime 1 < q < q, giant connected components may be unable to synchronize when large enough rewiring times T are considered.This effect is indeed due to the interaction between the giant component topology and the quenched disorder.For q 1 the giant component should be characterized by a large number of bridges (i.e.links whose deletion would split the giant component in two disconnected parts).When these bridges insist on nodes characterized by extreme natural frequencies (i.e.lying in the tail of the distribution P (ω)) which do escape partial synchronization, they act as effective obstacles to information spreading, splitting the topologically connected giant component into different synchronized subcomponents which are, however, not mutually synchronized.This mechanism, which clearly prevent macroscopic synchronization to emerge in the slow switching regime, is however going to become less and less important as the connectivity q is increased and the number of bridges in the giant connected component is reduced, allowing for a more efficient information flow, eventually leading to global synchronization as q > q.
Fast switching, on the other hand, allows information to travel through the network by rearranging the giant cluster quickly enough, preventing instantaneous bridges from acting as effective roadblocks.Therefore, for 1 < q < q, a transition to macroscopic synchronization is eventually observed as the rewiring time is decreased.A precise analytical estimate of T c (q) in this regime, however, is beyond the scope of this work.

IV.B. Critical behavior
We have finally verified that, as expected, the phase transition to synchronization belongs to the usual Kuramoto model class.In Fig. 7 we report numerical simulations for both (relatively) slow and fast rewiring times T , showing that the average order parameter follows the usual Kuramoto model scaling, ∆ ∼ q − q c (T ) for q > q c 24 , with q c (T ) being the (T dependent) critical connectivity parameter.Numerical results suggest this to be true for any finite q c , as expected given that synchronization seems to be essentially guided by the properties of the globally connected time-averaged connectivity matrix.

V. CONCLUSIONS
We discussed a time-varying network of heterogeneous Kuramoto phase oscillators characterized by links being randomly switched on and off with a Poissonian probability distribution.The network dynamics exhibits three well defined time scales associated respectively to local synchronization, local desynchronization and effective rewiring, whose separation is controlled by model parameters.Numerical simulations and analytical arguments show that this system is able to achieve statistically stable macroscopic synchronization even for arbi- trarily small net connectivity (i.e., for a sparse and infrequent coupling among the oscillators), provided sufficiently fast switching and strong couplings are considered.In the formal fast switching limit, T → 0, we have argued that the synchronization dynamics is fully captured by the time-averaged connectivity matrix, suggesting that results from the averaging theorem can be applied to our Kuramoto setup.At finite T , on the other hand, our analytical arguments, based on the comparison between the different timescales at play, are indeed able to predict with a satisfactory precision the synchronization transition line in the small connectivity regime.This switching-induced synchronization maintains the same qualitative characteristics of its static counterpart, such as the Kuramoto order parameter scaling ∆ ∼ √ q − q c at q q c 11 .For larger connectivity values, beyond the onset of giant connected components (q > 1), we have finally shown that the interaction between instantaneous connectivity topology and quenched disorder may prevent the onset of synchronization for sufficiently slow rewiring times and large but finite couplings.
Our findings are primarily intended as a theoretical contribution to the field of synchronization in time-varying complex networks, and in particular nonequilibrium synchronization models with alternative mechanisms giving rise to synchronization.However, we can still envisage several lines of potential applications for our results.For instance, one can think of engineered systems of heterogeneous (and possibly mobile) units 37 -a simple paradigma for the "Internet of Things" 38 where maintaining constant connectivity could be costly, yet the system is still required to exhibit synchronization or other collective properties.Our model could help develop alternatives to constant interactions, able to gener-ate the same collective dynamics albeit a sparse and seldom connectivity.Also, as mentioned earlier, this setup could be seen as a crude model for social interactions.As such, our model could be used to qualitatively model the emergence of consensus in a community where different individuals are only interacting with a few of their contacts at any time.
Finally, these results open several avenues of future work.Rather than rewiring links at random, one can consider a network where links are rewired preferentially to nodes with a similar instantaneous dynamical state, thus favoring interactions with "like-minded" individuals.This set-up could be used, for instance, to investigate qualitatively the "echo-chambers" phenomenon in social media, which as been recently suggested to be a possible source of an increased polarization in political opinions.To this regard, one can also wish consider different distribution of quenched frequencies, such as uniformly distributed ones, which do not favour middle natural frequencies (i.e.opinions) as the Gaussian one.

1 RFigure 1
Figure1Lin-log plot of the order parameter R as function of time for a network with N = 10 4 and q = 0.8, with ε = 8 and σ = 1.The three curves correspond to three different values of switching time T : green T = 6.28, red T = 0.63, and black T = 0.31.

Figure 2 (
Figure 2 (a) Stationary values of the order parameter ∆ (colorcoded according to the right vertical bar) as function of the rewiring frequency (T ) and network connectivity (q).Simulations have been performed for a network of N = 10 4 Kuramoto oscillators with ε = 8 and σ = 1.Values have been averaged over Ω = 10 different realizations.For small values of q the system is strongly dependent on the value of the rewiring time T and the phase diagram shows a clear transition from partial synchronization to disorder as T is increased beyond a critical value Tc(q).At larger q values, the transition approaches a vertical asymptote, roughly located at q = q = 1.66(6) (dashed black line).For q > q the dynamics achieves partial synchronisation regardless of the value of T .(b) Zoomed view of panel (a) in the range q ∈ [0.5, 1.7].

6 ∆Figure 3
Figure 3 Finite size determination of the transition point Tc for q = 0.8, = 8 and σ = 1.The average order parameter ∆(N ) is evaluated for two different system sizes, respectively N 1 = 1000 (black dots) and N 2 = 2000 (red dots).Tc is estimated as the midpoint between the largest value of T such that the values of ∆(N 1 ) and ∆(N 2 ) overlap, and the smallest value of T such that the scaling ∆(N 1 )/∆(N 2 ) ≈ √ 2 is satisfied.To facilitate the comparison, the black dashed line marks the value ∆(N 1 )/ √ 2. Vertical dashed lines mark the estimated transition point (red) and its confidence interval (green).Error bars report the standard error for the average computed over Ω = 20 independent realizations.

Figure 4
Figure 4 Critical rewiring time Tc as function of q for q < 1 for σ = 1 and different values of the coupling constant (increasing along the cyan arrow).Respectively, from left to right: ε = 32 (blue circles ), ε = 16 (red circles) and ε = 8 (black circles).Error bars give the estimated upper and lower boundaries for Tc(q) as discussed in the main text.The dashed straight lines (same color coding) mark the linear prediction of Eq. 28 (see Section III.B).(Inset): The slope s of each curve, evaluated by linear regression of the main panel data, is compared with the theoretical estimate s = π/(2 √ 2) (see Eq. (14)).Data has been averaged over Ω = 20 different realizations and error bars measure one standard error.

Figure 5
Figure 5 Numerically estimated transition line intercepts q 0 as a function of the coupling constant (black dots) are compared with the analytical prediction given by Eq. (27) (dashed red line).Data as in Fig.4.Error bars represents error in the linear extrapolation process (see main text).