Arxiv:1612.06644V2 [Physics.Soc-Ph]

arXiv:1612.06644v2 [physics.soc-ph] 7 Nov 2017 rbe fcnrlig(oil ewr yaisi by is dynamics network the study (social) conveniently controlling to of way problem one VM, the specifically, nodes hub nodes. procedure driver this as net- avoided with whole are the Surprisingly, minimum control the structurally [11]. identify to work to nodes share driver nodes) not of end do set which and links start of same set maximum the match- the maximum i.e., the set, setting, ing as problem (such heterogeneous this arose From for concepts interesting strategy [14]. control proposed [13]. simple was in a networks discussed result, of been a applicability has As the networks [12], real-world to problem which eigenvalue Specifi- an reformu- as [11–14]. was lated controllability theory network of control problem ad- the of commonly cally, is concepts opinion) by consensus dressed certain a (like state the [7–10]. takes opinion (node) discrete individual structure, a an network of of a role state on health models the epidemic where are related toy Closely VM instructive the [3–6]. science net- between to social and bridge social physics in a on models represents opinions is and (VM) individual dy- works model of the voter describing behavior model the namical idealized others, systems investigated Among widely complex a and 2]. network [1, of theory application popular a rmteprpcieo pno yaisad more and, dynamics opinion of perspective the From desired a to system given a driving of problem The is networks social on dynamics opinion of study The osa nttt o lmt matRsac,PO o 01 60 Box P.O. Research, Impact Climate for Institute Potsdam osa nttt o lmt matRsac,PO o 01 60 Box P.O. Research, Impact Climate for Institute Potsdam osa nttt o lmt matRsac,PO o 01 60 Box P.O. Research, Impact Climate for Institute Potsdam tchl eiineCnr,SokomUiest,Kräf University, Stockholm Centre, Resilience Stockholm eltyEet nOiinDnmc nteAatv oe Mo Voter Adaptive the in Dynamics Opinion on Effects Zealotry ubro oncin ihohrnds respectively. adjusting dominan nodes, complete other nodes: a with remaining avoid connections the to of for number order d strategies In average possible the zealots. two and the probability of opinion that rewiring zealots’ with the link of c the spreading of on nume the set dependent the that low-dimensional verify imply a We results yielding Our opinion. model the their of of spre proximation spreading increa to efficient an sufficient exhibit an must is for zealots zealots point, simu a of transition Numerical density the of Beyond low case connectedness. a the and threshold opinion to density investigate tation and zealotry initial nodes of their the con on of study only state the the system, extend with the evolving we throughout Here, opinion o fixed network. studies their Past spread networks. to social time-evolving on processes h dpievtrmdlhsbe ieysuida concept a as studied widely been has model voter adaptive The eateto hsc,Hmod nvriy etnt.15 Newtonstr. University, Humboldt Physics, of Department .INTRODUCTION I. aclP lme n acWiedermann Marc and Klamser P. Pascal Dtd oebr8 2017) 8, November (Dated: oahnF Donges F. Jonathan ekV Donner V. Reik pno,btti rcs svr nieya oprdto compared their as change unlikely very still nodes. as is can other zealots process form they stubborn this weak work, but fully present a opinion, to the contrast as In in seen studied 6]. be [5, confi- also zealots bounded can of in models extremists opinion might that It dence note [17]. to an state receiving interesting their nodes be pins as which seen signal be input can external the- zealots control of to type analogy (opinion) In this ory, VM. state the net- dynamical of take a evolution their we the of during study, change nodes present never as the zealots that define In work and forma- 16]. viewpoint 15, opinion latter the 4, different [3, of models context exten- tion the been have in one zealotry [4] studied of maintain sively types fully Both or opinion. [3] favor stub- specific either are zealots who Here, agents born system. the to zealots introducing sst isteoealoiin oevr te opinion other Moreover, in opinion. overall 6], the extrem- for bias placements [5, to good in the be ists extremists made to spread found with were been hubs to which models also where have confidence from observations bounded positions Similar good opinion. be were hubs to zealots high-degree found whereby of decreased. 19], placement [16, drastically investigated topological where was was optimal zealots consensus the of reach density Moreover, study to specific latter time a The at the regular, [18]. transition on graphs a zealots random found of number and finite [4] single a a complete with with and lattices regular [3] on zealot studied been already has ntecneto h lsia M h ffc fzealotry of effect the VM, classical the of context the In rkt2,141 tchl,See,EU Sweden, Stockholm, 19 114 2B, triket e erea oprdt riaynodes ordinary to compared as degree sed h ffc fzaos .. oe aiming nodes i.e., zealots, of effect the n pedn nue yzaosdepending zealots by induced spreading dteroiint h hl network. whole the to opinion their ad ain eelta eo h fragmen- the below that reveal lations h rbblt frwrn n/rthe and/or rewiring of probability the nteaatv oe oe sstrongly is model voter adaptive the in 3 41 osa,Gray Uand EU Germany, Potsdam, 14412 03, 2 3 41 osa,Gray Uand EU Germany, Potsdam, 14412 03, 2 uldodnr ieeta equations. differential ordinary oupled ge fnra oe ncomparison in nodes normal of egree ia nig sn enfil ap- mean-field a using findings rical 3 41 osa,Gray EU Germany, Potsdam, 14412 03, 2 iee h oe oe nastatic a on model voter the sidered eo h elt’oiin hr are there opinion, zealots’ the of ce 28 eln emn,EU Germany, Berlin, 12489 , dpientoktplg co- topology network adaptive n a oe o pno formation opinion for model ual del 2 formation models like the majority voter model [15] are dom and possess additional links and therefore an ex- also known to exhibit particularly rich dynamics on in- cess degree compared to ordinary nodes, which is ben- troduction of zealots. eficial for the spreading, as former related studies sug- All aforementioned studies assumed dynamics on a gest [5, 6, 16, 19, 27]. The excess degree is not only fixed network structure. However, the results obtained do motivated by the expected effect on the spreading but not apply to systems with a time-evolving network topol- also from real world examples [28]. For instance, elec- ogy [20]. Adaptive network models like the adaptive tion campaigns aim to reach as many voters (nodes) as voter model (AVM) [21] or other more realistic models possible. Following their mission, campaigners are not of opinion formation [22] generate a time-dependent net- convinced by voters, there is only a unidirectional in- work structure through a feedback mechanism between fluence of campaigners on voters. Additionally, cam- topology and node (agent) states. The AVM extended paigners reach an effectively increased degree in social the classical VM by giving the nodes the possibility to networks due to their professional outreach efforts [28]. break an existing link and reconnect to a like-minded Another example for zealotry are lobbyists intervening node which results in a temporally evolving network. in political processes. The number of zealots and their The standard AVM has been analyzed in terms of mean- excess degree can be interpreted as a measure of the re- field theory which revealed the existence of two absorbing sources that have been invested to pursue the campaign. states, namely the active and frozen (fragmented) state, Another important issue is how easily a system can be where for finite systems, the active state asymptotically controlled and what needs to be changed in order to in- becomes the consensus state (exhibiting a giant compo- crease its resilience against external pressure or corrup- nent with a single opinion only) [23, 24]. Further, it was tion. found that minor changes in the microscopic update rules After the description of the model and methods in can reduce or increase the time necessary to reach a final Sect. II, we focus on the question how the zealot opin- state due to network topology-state feedbacks [25]. ion is spread over clusters of different sizes in Sect. III. In the context of temporal networks (i.e., networks Thereby, we observe the emergence of subgraphs with a with links that are only present intermittently, includ- significantly larger mean degree than that of the whole ing the adaptive networks investigated in the present network. We identify three different parameter regimes, work as a specific class), controllability was investigated, closely related to the phase transition in the AVM, in alongside other studies, by a time-respecting path-based which different effects lead to a significant increase in the method [26] and by an analytical approach combined spreading efficiency of the zealots’ opinion. Our numer- with graphical tools [27]. The latter study revealed a pos- ical microscopic results are further supported by an an- itive relation between the aggregated degree of a node, alytical macroscopic approximation of the AVM includ- the number of interactions during a given time, and the ing zealots (Sect. IV). Ultimately, we explore the conse- size of the subset which is controlled by it. Both stud- quences of two different adaptation rules, in which the ies quantified the controllable subset by the influence of zealots either do not rewire at all or obey heterophilic a single node, assuming linear dynamics and consider- rewiring only (Sect. V). All obtained results are discussed ing networks which are statistically equivalent at dif- and conclusions drawn in Sect. VI. ferent times. However, it has remained unclear so far how the control of the AVM could be best achieved or avoided, because its dynamics are strongly nonlinear and II. METHODS the network is in certain parameter regimes evolving in such a way that it is not statistically equivalent at dif- A. Model description ferent times. Furthermore, previous studies on network controllability commonly addressed only temporal net- We study the AVM in the version originally formulated works without feedback between topology and node state, by Holme and Newman [21], but thoroughly extended by which is a key property of the AVM. Despite the resulting introducing zealots with excess degree. Here, the two differences between the AVM and controllability studies processes governing the opinion dynamics in a network on other temporal networks, we are confident that the with N nodes, also referred to as voters, and G different combination of both aspects is a very promising field of opinions are the change of node opinions and the link research. To our best knowledge, the concept of zealots, rewiring process, the mathematical formulations of which widely studied in the static VM, has not yet been ap- will be presented in the following. As the most crucial plied in the AVM. Since the latter is a relevant concep- parameter of the resulting adaptive network model, the tual model for opinion formation on temporal networks, rewiring probability φ is considered as the fraction of which obeys relatively simple rules allowing for a fair de- cases in which the latter process takes place instead of gree of understanding of the resulting dynamics, we see a the former. Each node i initially possesses an opinion great interest in addressing the issues mentioned above. gi, which is on average the opinion of a number of γ0 = Consequently, this paper addresses the efficiency of N/G nodes. For each node i, we define the set of nodes control by zealot opinion spreading (ZOS) in an exten- Si = Sgi \{{i} ∪ Ni} by excluding from the set Sgi of sion of the AVM. Here, the zealots are chosen at ran- nodes having opinion gi the node i itself and the set of 3 its direct neighbors Ni = {j : Aij =1} (here, Aij denote the average number of nodes initially holding a certain 0 the entries of the network’s adjacency matrix at a given opinion different from gz to γ = γ0(1 − nz). point in time, where we have suppressed the associated time index for brevity). The dynamic update cycle of the model is then B. Perspectives on zealot opinion spreading described as follows: Zealots are nodes with excess degree intending to Step 1: Randomly select a node i. If the degree ki of spread their fixed opinion gz to as many other non-zealot node i is zero do nothing, otherwise randomly select nodes as possible. In this special case, we investigate the a neighbor j ∈ Ni. zealot opinion spreading (ZOS) process characterized by Step 2(a): With probability φ, delete the link to j ∗ the fraction of nodes nz(t) holding the zealots’ opinion, and rewire to a randomly selected node j ∈ Si with ∗ which is a special opinion since it is (unlike the others) j 6= j (rewiring). If |Si| = 0 do nothing. always present due to the zealots that cannot become Step 2(b): With probability (1 − φ), node i imitates convinced by others. The ZOS efficiency is defined as the opinion of node j and, thus, gi → gj (imitation). the fraction of nodes nz(tc) holding opinion gz when the In contrast to the classical formulation of the AVM, final state is reached. in this formulation no multiple links and self-loops are As emphasized above, the considered problem can also be viewed from a controllability perspective. Here, possible due to node rewiring to the set Si (instead of zealots are just normal nodes that are influenced by a rewiring to the set Sgi as in the standard AVM). Note that step 2 is applied regardless of an existing opinion constant control signal b(t)= b which fixes their opinion conflict. The algorithm is iteratively repeated until a to gz. The excess degree kx of zealots can be viewed as a “topological input signal” applied to the network only time tc, where the final state is reached in which only like- once at t = 0. In real world applications, constraints exist minded nodes are connected to each other (Ni ⊆ Sgi ∀i). in terms of resource limitations (campaigners or lobbyists It shall be stressed, that the rules employed in this model variant are based on node selection [21, 24, 29] need to be paid) or ideology (not everyone wants to be a zealot), which motivates us to limit the number of nodes instead of link selection [30–32], because we consider it more realistic for a social network that agents (nodes) which can receive the said control signal b(t). spend on average the same time for communicating with others. The number of nodes N and the total number of III. NUMERICAL RESULTS links M stay constant over time, which implies that the 2M mean degree k0 = N of the network is kept fixed. As already mentioned in the introduction, in this work, The AVM, obeying the discrete-time update rules as we consider zealots as nodes that cannot be convinced. described in Section II A, starts evolving at t = 0 from an For the sake of simplicity, we assume here that all zealots Erd˝os-Rényi random graph, which is known to exhibit a giant component if k0 ≥ 1 [1, 33]. A phase transition oc- carry the same specific opinion gz. The set of zealots Sz is created at the start of each simulation by randomly curs when increasing the rewiring probability φ to a level 0 at which the giant component of the initial network van- declaring a fraction of nodes nz to be zealots. Their key property of having an immutable opinion is ensured by ishes and homogeneous clusters emerge which are formed modifying the last step in the above scheme: by like-minded nodes. This fragmentation transition occurs at a critical rewiring probability φc (for fixed mean Step 2 (b) With probability (1 − φ), if i ∈ Sz do noth- degree k) or at a critical mean degree kc (for fixed φ). ing, otherwise node i adopts the opinion of node j. In what follows, we will mostly follow the strategy of In comparison with all non-zealot nodes, the initial varying φ but keeping the mean degree fixed. mean degree of zealots can be further increased by an excess degree kx to kz(t =0)= k0 +kx, which describes additional links that are randomly connecting each zealot to A. Fragmentation transition non-zealots. This excess degree (which is considered here to be the same for each zealot to simplify the following One way of identifying the fragmentation transition analyses) is motivated by the additional efforts of cam- and the associated parameter value φc is based on crit- paigners or lobbyists to convince as many nodes as pos- ical slowing down [21] indicated here by a maximum of sible. Starting with a specific configuration, the degree the coefficient of variation Vt = σtc /tc of the convergence of individual zealots then changes over time according to time as a characteristic parameter, where tc and σtc de- the considered rules of the update cycle. Let the initial note the empirical mean value and standard deviation of 0 fraction of nodes holding opinion gz be ngz (0) = nz. To the convergence time estimated from a sufficiently large simplify the notation, ngz will be denoted in the follow- ensemble of independent realizations of the considered ing as nz = ngz . Note that introducing the excess degree AVM variant [21]. In the present work, Vt is estimated 4 kx increases the mean degree of the whole network to as a function of φ for N = 800, k0 = 4, γ0 = 10 and 10 0 k = k0 +2nzkx, and the introduction of zealots changes simulation runs (Fig. 1). If no zealots exist (Fig. 1(a)), 4

a) in degree, while for the same φ, some simulations exhibit 0.7  negative degree correlations while others show positive 0.6 ones. The negative degree correlations imply that hubs, 0.5 which are mostly zealots due to their excess degree kx, 0.4 are mainly connected to low degree non-zealot nodes. Be- 0.3 cause randomly chosen nodes mimic the opinion of a randomly chosen neighbor, the opinion of high degree nodes 0.7  b) gets imitated more often, resulting in a faster spread of 0.6 the zealot opinion and a shorter tc. In turn, positive de- 0.5 gree correlations result in a longer t , since in this situa- c 0.4 tion, hubs are frequently connected with other hubs and, 0.3 thus, their opinion gz reaches other hubs (mostly being zealots with the same opinion themselves) with elevated

0.7  c) probability in a short amount of time, but takes much 0.6 longer to aﬀect the rest of the graph. Consequently, the

0.5 secondary peak of Vt indicates the rewiring probability 0.4 that is necessary to compensate for the eﬀect of the intro- 0.3 duced degree heterogeneity. Accordingly, the coexistence

of both types of degree correlations triggers a large σtc 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 and therefore a large V . However, this effect has minor Rewiring probability  t relevance for the ZOS because the secondary maximum of Vt is located below the fragmentation threshold. FIG. 1. (Color online) Coefficient of variation of convergence time Vt, indicating a phase transition at φ = φc where the giant component vanishes. The different parameter settings B. Cluster size distributions 0 0 are (a) no zealots nz = 0, (b) zealot density nz = 0.01 with no 0 excess degree kx = 0, and (c) zealot density nz = 0.01 with an The most straightforward approach to maximize the excess degree of kx = 20. The red lines in (a,b) at φc = 0.45 and in (c) at φc = 0.505 indicate the approximate positions ZOS efficiency nz(tc) is to dominate the largest connected of the phase transition, which is shifted if the excess degree component in the final state. In the following, we will 4 kx is nonzero. Distributions are computed from n = 10 runs refer to connected communities with homogeneous node on graphs with N = 800 nodes and a mean degree (without state (opinion) as clusters. In Fig. 2, the resulting fre- the zealots’ excess degree) of k0 = 4. quency distribution P (s) of cluster sizes s (i.e., the number of nodes in a cluster) is shown for rewiring probabilities below (a,b,c), close to (d,e,f) and above the frag- the phase transition occurs at φc ≈ 0.45, as expected mentation threshold (g,h,i) of the three cases presented from previous studies [21]. If only few zealots are intro- in Fig. 1(a,b,c). 0 duced (nz =0.01) (Fig. 1(b)) the phase transition point The distributions without zealots are characterized by is not altered. However, declaring them as hubs with an a giant component below the critical point φc (Fig. 2(a)), additional excess degree of kx = 20 (Fig. 1(c)) shifts the by a power law behavior close to the transition (Fig. 2(d)) phase transition to φc =0.505. Since the mean degree k and by the absence of a giant component above the frag- has a strong impact on the transition [21], this shift in mentation threshold (Fig. 2(g)), where the cluster size s φc can be explained by its change k0 → k = 4.4. Simu- is distributed around γ, indicated by a vertical line in lations with k0 =4.4 and no zealots (not shown) indeed Fig. 2(g). Since no zealots and therefore no opinion gz perfectly reproduce the observed shift, implying the re- are present the ZOS efficiency is always zero. distribution of the additional links to the whole network Below the fragmentation transition, zealots dominate in t ≤ tc. the giant component in both cases with and without ex- Another phenomenon is observed if the zealots’ excess cess degree (Fig. 2(b,c)) and on average spread their opin- degree is increased, which manifests in a second broad ion to a fraction of nz ≈ 0.98N nodes. Close to φc, the maximum of Vt in Fig. 1(c) at about φ = 0.08. This distribution of clusters having the zealots’ opinion gz (red secondary maximum is shifted towards larger φ if kx is squares) without excess degree (Fig. 2(e)) is similar to the further increased (not shown) and is explained by the total distribution (blue triangles) with the difference that 0 assortativity in degree, the degree correlation of neigh- convinced clusters have a size of at least Nnz, indicated boring nodes, at the final state [34]. The initial state has by a red dashed line. The excess degree (Fig. 2(f)) causes negative assortativity because the hub-zealots are mostly a peak at larger cluster sizes. The size distribution of con- 0 connected to lower degree nodes at t = 0. The secondary vinced clusters starts at a size significantly larger than nz maximum coincides with the peak of the coefficient of and coincides at larger cluster sizes with the total distri- variation of the assortativity (not shown) and low values bution, which shows that the corresponding maximum of its mean. Thus, the network is on average uncorrelated consists solely of convinced clusters. At the fragmen- 5 tation transition the cluster distribution is expected to obey a power law (Fig. 2(d,e)), which is disturbed by the peak if zealots posses an excess degree (Fig. 2(f)). Conse- quently, it can be assumed that the excess degree is split- ting the system into gz-dominated subgraphs which have a larger mean degree than subgraphs that are not influenced by the zealots’ opinion and excess degree. Hence, the gz-dominated subgraphs effectively lie below or close to the local fragmentation transition and, therefore, tend to form larger clusters (Fig. 2(f), red squares). These subgraphs give rise to a five times larger ZOS efficiency as compared to the case of zealots without excess degree (Fig. 2(e)). The gz-dominated subgraphs are also present above the fragmentation threshold, which is indicated by the fact that the previously discussed maximum of the cluster size distribution does not vanish suddenly but is grad- ually shifted towards smaller cluster sizes (Fig. 2(i)). In FIG. 2. (Color online) Frequency distributions of cluster size P (s) at t = tc for all clusters (blue triangles) and for the the case without excess degree (Fig. 2(h)), far above the fraction of convinced clusters (red squares) computed from 4 transition point the size distribution of convinced clusters n = 10 simulation runs on graphs with N = 800 and γ0 = 10. 0 0 exhibits a maximum at nz (red line), while in the case The distributions are computed for no zealots nz = 0 (a,d,g), 0 with excess degree (Fig. 2(i)) this maximum is shifted for zealot density nz = 0.01 and no excess degree kx = 0 0 towards larger cluster sizes, resulting in almost twice as (b,e,h), and for zealot density nz = 0.01 and excess degree high ZOS efficiency as compared to the case with kx = 0. kx = 20 (c,f,i). The rows represent different values of the In summary, the cluster size distributions reveal that rewiring probability φ: φ = 0.04 ≪ φc (a,b,c), φ = φc(k) ZOS close to and above the fragmentation transition is (d,e,f) and φ = 0.96 ≫ φc (g,h,i). Blue solid and red dashed vertical lines indicate γ0 = 10 and initial number of zealots more efficient if hub-zealots are present due to a forma- 0 Nnz tion of convinced subgraphs with a larger local k. Below , respectively. Note the black arrows indicating the giant component. the transition point, the mere existence of zealots is sufficient to reach a maximum ZOS efficiency.

10.25 a) b) c) 0.8 C. Systematic parameter study 9.0 7.75 9.0 0.6 4.4 5.25 0.4 So far, the only parameter settings considered have 2.75 0.2 0 been no zealots and zealots with initial density nz =0.01 0.25 with and without excess degree kx = 20. We now study 10.25 d) e) f) a much wider set of parameters regarding the ZOS at 0.4 7.75 0.3

final state nz(tc) (Fig. 3). The horizontal and vertical 5.25 0.2 axes in Fig. 3 represent the excess degree kx and density  0 2.75 0.1 nz of zealots. The parameters k0 = 4 and N = 800 are 0.25 kept constant. Three different rewiring probabilities are 10 20 30 40 10 20 30 40 10 20 30 40 considered: φ = 0.4 < φc(k0) (Fig. 3(a)), φ = 0.5 > Excess degree φc(k0) (Fig. 3(b)) and φ = 0.75 > φc(k0) (Fig. 3(c)). Note that the fragmentation threshold depends on both kx and φ, and the critical rewiring probability φc(k0) = 0.45 can therefore only be considered as a reference point FIG. 3. (Color online) (a,b,c) Fraction of convinced nodes at the final state nz(tc) and (d,e,f) the former minus nzk(tc) for kx = 0. Thus, for a fixed φ there exists a critical mean 0 in which the kxNnz edges are distributed homogeneously degree kc. If this mean degree is exceeded, the system 0 instead to the Nnz zealots (no zealot-hubs) resulting in is below its fragmentation threshold, and, hence, in a ∆nz(tc)= nz(tc)−nzk(tc), both as a function of excess degree 0 regime where ZOS profits from the presence of a giant kx and zealot density nz. The rewiring probability is increased component. from (a,d) φ = 0.4 via (b,e) φ = 0.5 to (c,f) φ = 0.75. The Below the fragmentation transition (Fig. 3(a)), we find black (b,c,e,f) and white (c,f) lines indicate contours where that an increase in either the excess degree or the num- k = kc and kz = kc, respectively, highlighting the fragmen- ber of zealots strongly increases the ZOS efficiency, which tation transition, with kz being the mean degree of gz nodes. quickly maximizes because a giant component is easily All results have been obtained from 200 simulation runs on dominated. Sufficiently far below the transition point networks with N = 800 nodes and a mean degree without (e.g., for φ = 0.1, not shown), every individual configu- zealots of k0 = 4. 6 ration results in maximum ZOS. does not favor nodes of a specific opinion. In turn, Just above the transition point (Fig. 3(b)), the pres- marked differences emerge if the excess degree is fixed ence of additional zealots increases the ZOS efficiency as to nodes of a single opinion (in our case, the zealots). expected. However, ZOS efficiency quickly saturates to a maximum if together with an elevated excess degree, the critical degree kc = 4.4 is reached where the giant IV. MACROSCOPIC APPROXIMATION component fragments. Finally, far above φc(k0) (Fig. 3(c)), the fragmentation For a model similar to that studied in the present transition is reached at kc = 9 (black contour line) with, work [29], it was recently shown that a mean-field ap- of course, larger amounts of zealots and excess degrees. proximation can be performed to derive analytical results It is remarkable that already below the corresponding by considering solely pairwise interactions and under the transition point, ZOS efficiency increases to multiples of assumptions that the network is large and fully connected 0 nz. This effect can be explained by the convinced clus- at initiation. In the latter work, only two distinct node ters that have larger mean degree than unconvinced ones states (aka opinions) have been present. This assumption (compare Fig. 2(f)). Note that the mean degree of gz is adopted in the following by treating all opinions differ- nodes kz at the final state crosses kc already at a smaller ent from gz as equivalent, i.e., as the opinion of the others zealot density and excess degree than the total mean de- go. Following [29], this simplification reduces the prob- gree (the case kz = kc is highlighted by a white contour lem to three coupled differential equations for the time line in Fig. 3(c)). Those subgraphs with larger mean evolution of three macroscopic properties of the model: degree are a key finding of this study and can be inter- the fraction of nodes nz holding opinion gz and the av- preted as a community in which intense discourses were erage numbers of links per node mzz (moo) among nodes triggered by the excess degrees of the zealots and their holding opinion gz (go): opinion. Note that the critical mean degree kc is estimated by reproducing Fig. 1(a) with a different k0 until the peak of V coincides with the value of φ used in dnz z 0 o t = (1 − φ) noPo − (nz − nz)Pz , (1) Fig. 3(b,c), respectively. dt 0 In order to properly interpret the results discussed dmzz o z nz − nz o = φnzPz + (1 − φ) Po mzo − 2Pz mzz , above, it is important to understand to which extend dt nz the observed emergence of densely connected subgraphs (2) responsible for ZOS efficiency originates from the pres- 0 dmoo z nz − nz o z ence of zeolots with distinct excess degree as opposed to = φnoPo + (1 − φ) Pz mzo − 2Po moo . dt n mere effects of an elevated mean degree k of the whole z (3) network. For this purpose, Fig. 3(d,e,f)) presents the difference ∆nz(tc) = nz(tc) − n (tc) between the ZOS zk Here, mzo = Mzo/N where Mzo is the number of “ac- efficiency nz(tc) of our AVM variant as discussed above tive” links between nodes of distinct opinions gz and go, and the ZOS efficiency n (tc) that would arise if the 0 zk respectively. Note that the excess degree kx enters the additional kxNnz links were distributed homogeneously above equation via the fraction of initially active links among all nodes of the networks instead of assigning them 0 o mzo(t = 0) = n (k0 + kx). P is the probability of a exclusively to the zealots. This difference is always posi- z z gz node to interact with a go node and is given by the tive and largest in the non-fragmented phase close to the heterogeneous mean-field approximation [29] as transition point (Fig. 3(e))). However, at large rewiring probabilities (φ = 0.75, Fig. 3(f)), the parameter range o o kz mzo for which marked differences between both settings are Pz = = . (4) kz 2mzz + mzo found extends over large parts of the parameter subspace o corresponding to the non-fragmented phase. These find- Here, kz is the mean degree of gz nodes and kz is the z ings demonstrate that the emergence of subgraphs with mean number of links from a gz node to go nodes. Po increased mean degree is mainly caused by the presence follows analogously by exchanging the indices. of hub-zealots. Equation (1) implies an increase of nz by go nodes be- We emphasize that qualitatively and quantitatively ing convinced by gz nodes and a decrease through non- similar results are observed (not shown) for situations in zealot gz nodes being convinced by go nodes. The first which the excess degree is distributed among randomly term in Eq. (2) describes an increase of mzz by gz nodes selected nodes of the network (random hubs) or among cutting their link to go nodes. The second term consists randomly chosen non-zealots only (non-zealot hubs). In on the one hand of an increase due to go nodes becom- fact, the results for the random hubs are quantitatively ing convinced, transforming the mean per-node number extremely similar to the case with homogeneously dis- of links between different opinions from mzo/no to mzz, tributed edges (Fig. 3(d,e,f)). Thus, we conclude that it and on the other hand a decrease by non-zealot gz nodes does not matter much if the mean degree is increased ho- changing their opinion and, thus, transforming the mean mogeneously or heterogeneously as long as the increase number of mzz links of one gz node from 2mzz/nz to mzo 7

(analogously for Eq. (3)). Also note that Eqs. (1)–(3) are closed, i.e., mzz + moo + mzo = k/2 and nz + no = 1. The main diﬀerence to the model in [29] is the presence 0.8 a) b) 0.8 of zealots included by reducing the fraction nz in the

0 0.4 0.4 convincing process by nz. For the model considered here, five fixed points can be identified as unstable or outside the regime of inter- 0.0 0.0 0 2 c) d) 6 est (0 ≤ nz ≤ nz ≤ 1, 0 ≤ mzz + moo ≤ k/2). A ODE ODE

=10 =N=400 two-dimensional manifold, which represents the consen- 1 4 =100 =N=800 sus/ﬁnal state, also satisﬁes the stationarity criterion: 2 =800 =N=1600 0 0 ⋆ k mzz = − moo (mzo = 0). (5)

2 1.5 5

Note that this manifold extends over all values of nz and 3 0.5 moo. Its linear stability properties are determined by the e) f) 1 ⋆ eigenvalues of the Jacobian at mzz, which are 0, 0 and 0.0 0.2 0.4 0.6 0.8 0.6 0.7 0.8 0.9 1.0 Rewiring probability  0 f(nz,moo,nz, k, φ)= 0 φnz 2moo(nz − 2nz)(φ − 1)+(nz − 1)nzφ FIG. 4. (Color online) Results of microscopic ensemble simu- + . lations (symbols) and macroscopic approximation (solid lines) 2moo − k 2moonz (6) using Eqs. (1)–(3). Mean values have been estimated from 200 simulation runs. (a) ZOS efficiency nz(tc), (c) density of Since two eigenvalues are zero, the manifold can- links between convinced nodes mzz(tc), (e) density of links not be asymptotically stable but (un-)stable if between unconvinced nodes moo(tc) at convergence time tc 0 with N = 800 and varying γ as indicated in (e). (b,d,f) Same f(nz,moo,n , k, φ) < 0 (> 0). In the following, it is z as for (a,c,e) but with varying N = γ0 and therefore only assumed that the links are, regardless of the initial con- z ⋆ one opinion different from g . The blue dashed vertical line ditions, homogeneously distributed at tc, i.e., moo = marks the approximate phase transition φ˜c (Eq. (9)). For all 0 k(1 − nz)/2, which simplifies the nonzero eigenvalue to runs, nz = 0.01 and kx = 20 have been kept fixed. In (a,c,e), 0 k0 = 4, while in (b,d,f), k0 = 15. Note the different scales ⋆ 0 nz(1 − φ) 2(1 + k)φ along the horizontal axes. f(nz,moo,nz, k, φ)=2 − − . (7) nz k The latter changes its sign at the microscopic model (Fig. 4). Both, the excess degree with k = 20 and the zealot density with n0 = 0.01, ˜ 0 2nz x z φc(nz,nz, k)=1 − 0 . (8) are kept constant across all simulations. A low mean 2(1 + k)nz − knz degree k0 = 4 is considered in Fig. 4(a,c,e). Simula- ˜ tions with this setting show a fragmentation transition Thus, the transition depends on k such that limk→0 φc = ˜ at φc ≈ 0.5, which disagrees with the analytical approx- 0 and limk→∞ φc = 1. The macroscopic approximation ˜ is strictly valid in the case N → ∞; otherwise, finite- imation of Eq. (9), φc ≈ 0.814, marked by a dashed ver- size effects may lead to deviations of the system from the tical line, by a relative error of δφc = 0.38. This large analytically calculated behavior. discrepancy is to be expected since a large mean degree The fragmentation threshold is approximated by is necessary for the mean-field approximation to reach a φ˜ (n = 1), because at the transition weakly connected good agreement with the full numerical results of micro- c z scopic simulations [35]. clusters split up regardless of their size and φ˜c increases In Fig. 4(b,d,f) the mean degree is increased to k0 = with nz. This results in 15. The numerical results indicate φc ≈ 0.85, while the ˜ ˜ 0 1 analytics give φc ≈ 0.94, which reduces the discrepancy φc(1,nz, k)=1 − . (9) 0 to a relative error of δφc = 0.096. The remaining dis- 1+ k(1 − nz/2) agreement can very likely be further reduced by including 0 Note that for nz ≪ 1 the transition can be considered as higher-order terms in the AVM [24], which is a promising 0 independent of nz, which is in perfect agreement with the task for future research. Note that far above and below 0 numerical findings in Fig. 1 with nz = 0.01 ≪ 1, where the transition point, both modeling approaches agree well 0 the transition is not shifted by changing nz (Fig. 1(b)) with each other. 0 but by a change in kx (Fig. 1(c)), since k = k0 +2nzkx. Additionally, in Fig. 4(a,c,e) the validity of the sim- We compare the mean-field model and the microscopic plification of considering all opinions other than gz the model by computing the final state values of nz, mzz and same as go is checked by decreasing γ, which increases moo for Eqs. (1)–(3) by forward integration and by tak- the diversity of opinions. Below but close to the phase ing an average over the outcomes of 200 simulations of transition point, simulations with a higher diversity show 8 a reduced ZOS efficiency due to the higher probability of the latter. Accordingly, the second step in the update small groups to cluster and cut their links to the giant scheme of Sec. II A is then modified as follows: component, cf. the black circles in Fig. 4(a). Above and Step 2(a): With probability φ, if i ∈ Sz and |Si| 6= 0, sufficiently below the fragmentation transition, the ob- delete the link to j and rewire to a randomly selected tained results are robust for different γ. node of the set Sgo \ Ni. If i∈ / Sz and |Si| 6= 0, rewire In Fig. 4(b,d,f) we check for finite size effects by simu- to a node of the set S . Otherwise, do nothing. lating at N = 400, 800and 1, 600 while keeping the diver- i sity of opinions at its minimum (γ = N). The finite size Former studies on AVM variants with heterophily smoothens the transition from the giant component to found that the fragmented phase vanishes [24]. However, the fragmented phase, as can be seen in Fig. 4(b). Here, in contrast to these previous works which introduced het- simulations with the largest number of nodes (green tri- erophily to all nodes, we here only declare zealots as being angles) switch more suddenly to the fragmented phase heterophilic, whereas normal nodes remain homophilic. (compare black circles and green triangles) than such for Thus, the fragmented phase can still be reached if there smaller networks. exist no links to any gz-node, i.e., mzz = mzo = 0. We note that the excess degree is responsible for the In Fig. 5, the coefficient of variation of the conver- observed slow decrease of nz(tc) after the fragmentation 0 gence time, Vt, is analyzed for both types of zealots for transition until it reaches nz at φ = 1 (Fig. 4(a,b)). If the case without excess degree only, since we focus here 0 the kxNnz additional links are instead distributed ho- on the effect of different zealot update schemes rather mogeneously among all nodes, resulting in an elevated than on the effect of hub zealots (cf. Fig. 1(b)). Inter- degree without distinct zealot hubs, the ZOS efficiency estingly, φc ≈ 0.45 remains valid for the case of passive is markedly reduced already at the transition point (not 0 zealots (Fig. 5(a)). Furthermore, the difference ∆nz(tc) shown), which can also be seen in our microscopic simu- between the ZOS efficiency in the case of passive and lations for the ZOS efficiency difference between the two normal zealots, as shown in Fig. 3(d,e,f), does not show cases in Fig. (3(e,f)) close to the transition point. any differences between the two rewiring schemes over a wide parameter range (not shown). A possible explanation for the aforementioned finding V. ALTERNATIVE ZEALOT UPDATE SCHEMES is that, as expected, the unrewired links to go nodes increase the probability to interact with go because of an z increase in mzo, while keeping moo constant increases Po In the previous sections, we have treated zealots as and therefore the first term in Eq. (1), resulting in larger different from normal agents exclusively in terms of their values of n . However, m is reduced in the absence of additional property to be stubborn. Doing so, we have z zz homophilic rewiring of the zealots. Assuming that mzz is been in line with vast parts of the existing literature. decreased by the amount of links which increases mzo re- However, considering that real world zealots can be ex- o sults in an increase of Pz which reduces nz, as can be seen pected to aim to maximize the spread of their respective from the second term of Eq. (1). In fact, P o also increases opinion, it appears unintuitive that they should follow z if mzz stays constant. Thus, these two processes compete almost the same update rules as non-zealots. Therefore, against each other, which could explain the observed ef- in the following we investigate two alternative update fect. However, in this perspective, we have neglected schemes, which should potentially increase the ZOS ef- (besides other aspects) that the considered missing mzz ficiency. For brevity, we leave any in-depth analysis of links connect to zealots. Thus, non-zealot gz nodes are the resulting dynamics along the lines of the previous more likely to be convinced and stay go nodes than in sections as a subject of future research. the setting involving normal zealots. As a first possible scenario, we investigate the case of passive zealots that do not rewire at all and, thus, do In contrast to the aforementioned behavior, the fragmentation transition for heterophilic zealots (Fig. 5(b)) not cut links to go-nodes. This procedure then enhances the probability to convince the latter. In this case, the shifts strongly to φc ≈ 0.77. This shift causes an increase second step in the update cycle as presented in Sec. IIA of ZOS efficiency, since a giant component is present over is modified such that a larger parameter range in φ. This increased stability of the giant component is caused by the mixed behav- Step 2: If i ∈ Sz do nothing, otherwise follow steps ior of the gz nodes: while the non-zealots exhibit ho- 2(a) and (b). mophilic rewiring and therefore stay connected to the Passive zealots can be interpreted as agents who are zealot gz nodes, the zealots themselves reach out for other forced to a trade-off between the large number of con- go nodes, which could already be parts of an isolated com- nections kept to other agents and their ability to rewire ponent consisting only of like-minded nodes. These iso- to new agents. They benefit from their large degree but lated components were not reachable for normal zealots, in turn cannot manage to establish new links. whereas heterophilic zealots can “invade” these groups As a second case, we study heterophilic zealots that of nodes and convert parts of them or even the whole rewire exclusively to go nodes. This restriction is ex- groups. Note that non-zealot gz nodes provide the zealots pected to further enhance the efficiency of convincing with links via their homophilic rewiring. 9

0.7 a) ties and excess degrees, three regimes were identified in  0.6 which different effects allow a maximum or an increased

0.5 ZOS efficiency. Below the fragmentation transition, max- 0.4 imum ZOS is easily achieved by introducing solely a few 0.3 zealots without excess degree due to the emergence of a giant component, which allows for ZOS across the whole 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 network. Shortly above the fragmentation transition, an 1.6 b) increase of zealot density is insufficient to increase ZOS  efficiency. In addition, an excess degree needs to be in- 1.2 troduced and raised to a certain level to allow for the

0.8 system returning to the non-fragmented phase in which ZOS efficiency is quickly maximized. Thereby, either the 0.4 mean degree of a large amount of zealots is slightly in- 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 creased by the excess degree, or a small amount of zealots Rewiring Probability  is declared as hubs. Far above the phase transition, only a large density of hub-zealots is able to push the system to the non-fragmented phase. However, the formation FIG. 5. (Color online) Coefficient of variation of convergence of zealot-dominated subgraphs with increased mean de- time Vt, indicating a phase transition at φ = φc where the gree plays a crucial role and allows for increased ZOS giant component vanishes. Passive (a) and heterophilic (b) zealots are compared. The red line in (a) marks φc = 0.45 efficiency already far below the phase transition point. and in (b) φc = 0.77. Distributions are computed from n = Since these subgraphs emerge in the fragmented phase, 104 runs on graphs with N = 800 nodes, a zealot density of they only include a specific fraction of nodes, determined 0 nz = 0.01 with no excess degree kx = 0 and a mean degree of by the cluster size distribution. k0 = 4. We have macroscopically approximated the model by considering pairwise interactions. An analytical approximation of the phase transition point was found, which In summary, the two alternative zealot rewiring was validated by forward integration. However, the theo- schemes discussed in this section illustrate that there is retically approximated critical point is much larger than a large class of update schemes, which perform equally suggested by the numerical simulation of the microscopic well or even better than the one discussed in the previ- model. This discrepancy was reduced by considering ous sections in terms of ZOS efficiency. Especially the systems with larger mean degree. Previous studies [24] strong effect of heterophilic zealots and how non-zealot suggest the potential for further improvement if the ap- gz nodes provide links to them to enable the zealots to proximation is not only based on pairwise interactions. connect to not like-minded nodes demonstrate that there However, far below and above the phase transition point, can exist interesting feedback effects between two coex- analytics and numerics agree well with each other. isting rewiring schemes, which seems a promising field of Finally, we have studied the effect of two alternative future research. update schemes for the zealots. Here, it was shown that passive zealots, which do not rewire, perform as good as their normal counterparts in terms of ZOS efficiency. VI. CONCLUSIONS In turn, heterophilic zealots, which rewire to not like- minded nodes only, shift the fragmentation transition In this paper, we have introduced zealots with an in- strongly to larger rewiring probabilities and therefore creased mean degree, the excess degree, into the adaptive have a much larger ZOS efficiency. These update schemes voter model (AVM) and investigated how their fixed, uni- provide an interesting starting point for future research, form and new opinion spreads over a social network until which should investigate the effects of a coexistence of a full or fragmented final state is reached. The efficiency different update rules in the same model. of zealot opinion spreading (ZOS) has been quantified by The finding of zealot-dominated subgraphs with a the fraction of nodes holding the zealots’ opinion at the larger mean degree than in the rest of the graph in the asymptotic state. fragmented phase allows drawing the following conclu- After reproducing the results of a previous study [21] sion: Large communities (sub-networks) with agents en- by means of numerical simulation, a detailed comparison gaging in active discourse (larger mean degree than in of the resulting cluster size distributions below, at and far other communities) are likely to be targeted by the in- above the fragmentation transition revealed the existence terests and the resources of an already convinced group. of zealot-dominated subgraphs if the introduced zealots In our model, an active discussion can imply that there is exhibit an excess degree. These subgraphs are charac- an attempt to control the system. In order to avoid being terized by an elevated mean degree and the presence of controlled by an external opinion, each node should keep opinions enforced by the hub-zealots. links to nodes from other subgraphs. By investigating a wide range of initial zealot densi- How to maximize or minimize opinion spreading by 10 zealots is a relevant question in the context of the AVM this case, we have already studied two model variants as well as related models of opinion dynamics. Similar with passive and heterophilic zealots, whereby the lat- studies already focused on the static voter model [16, 19] ter results in a large increase of ZOS efficiency. A for- or on general dynamical systems [11–14, 36] on static mer study investigated how agents maximize their power networks. In the latter research, it was of crucial inter- (represented by a score function increasing with central- est if there exist specific nodes which have topologically ity and decreasing with degree, the diplomats dilemma) favorable or unfavorable positions to spread their opin- by following specific link creation and deletion strategies, ion, or if there exists a minimum set of nodes necessary assuming knowledge up to the second neighborhood [37]. to spread an opinion across the whole network. Along The most common strategy to increase the corresponding these lines, it is of interest to quantify the effect of net- performance was to delete the link to the direct neighbor work topology on ZOS efficiency in the AVM in future with largest centrality and add an edge to the neighbor research. An increase of the rewiring probability changes in the second neighborhood with largest centrality. This the topology faster and increases (below the fragmenta- could be one further update rule, next to many others. In tion transition) the convergence time. Consequently, the this context, an interesting question would be to identify larger the rewiring probability φ, the stronger the ini- specific update schemes which are robust across different tial topology is modified. Hence, topological effects are network types. only expected to play a role at low rewiring probabilities. In summary, we emphasize that the combination of However, as shown by our study for random graphs, in intervening in opinion adoption and rewiring processes, this regime opinion spreading is already maximized by different zealot update schemes, considerations regard- randomly positioning a small number of zealots without ing the role of complex network topology and generaliz- excess degree. For networks with a more complex topol- ing those approaches to more realistic models of social ogy than a random graph, the corresponding effect might network dynamics (e.g., [22, 38]) are promising fields of be quantitatively or even qualitatively different. Above future research. the transition point, especially the topological effects on cluster formation with increased mean degree are of interest. ACKNOWLEDGMENTS More generally, this study has also presented a strong motivation for further investigating the controllability of P.P.K., M.W. and R.V.D. have been financially sup- the AVM as an example for a nonlinear dynamical sys- ported by the German Federal Ministry of Educa- tem on a dynamic network. This is because the ZOS ef- tion and Research (BMBF) via the Young Investigators ficiency, on which we have focused our interest, does not Group CoSy-CC2 (grant no. 01LN1306A). J.F.D. thanks represent exactly the controllable subset of the zealots, the Stordalen Foundation (via the Planetary Bound- but might still be closely related to this concept. If the ary Research Network PB.net) and the Earth League’s ZOS efficiency would represent the subset of nodes con- EarthDoc program for financial support. M.W and trollable by the zealots, we could drive this subset from J.F.D thank the Leibniz Association (project DOMI- any initial state to any desired state within a time tc by NOES) for financial support. The authors gratefully an appropriate input signal. However, the ZOS efficiency acknowledge the European Regional Development Fund only shows that we drove the subset to a specific, not an (ERDF), the German Federal Ministry of Education arbitrary, state. and Research and the Land Brandenburg for support- This study has been based on the approach of intervening this project by providing resources on the high ing in the opinion adoption process in the AVM. Clearly, performance computer system at the Potsdam Insti- a complementary approach would be to interfere with tute for Climate Impact Research. The presented re- the rewiring process, which could imply to declare spe- search was conducted within the scope of the CO- cific links as unbreakable or harder to break. Also, the PAN flagship project on co-evolutionary pathways at creation of links between zealots by rewiring could be ex- the Potsdam Institute for Climate Impact Research cluded, assuming that campaigners or lobbyists have no (http://www.pik-potsdam.de/copan) and the Inter- interest in clustering among themselves. First, it would national Research Training Group IRTG 1740/TRP be easier to identify them as zealots and second, they 2014/50151-0, jointly funded by the German Research would waste their linkage resource needed for influencing Foundation (DFG, Deutsche Forschungsgemeinschaft) other nodes. and the São Paulo Research Foundation (FAPESP, Another way of intervening in the opinion formation Funda¸cão de Amparo à Pesquisa do Estado de São process is to change the update rules for the zealots. In Paulo).

[1] R. Albert and A.-L. Barab´asi, [2] M. E. J. Newman, SIAM Rev. 45, 167 (2003). Rev. Mod. Phys. 74, 47 (2002). [3] M. Mobilia, Phys. Rev. Lett. 91, 028701 (2003). 11

[4] M. Mobilia, A. Petersen, and S. Redner, [22] C.-F. Schleussner, J. F. Donges, D. A. Engemann, and J. Stat. Mech. Theory Exp. 2007, P08029 (2007). A. Levermann, Sci. Rep. 6, 30790 (2016). [5] F. Amblard and G. Deffuant, Physica A 343, 725 (2004). [23] F. Vazquez, V. M. Egu´ıluz, and M. S. Miguel, [6] D. W. Franks, J. Noble, P. Kaufmann, and S. Stagl, Phys. Rev. Lett. 100, 108702 (2008). Adapt. Behav. 16, 264 (2008). [24] D. Kimura and Y. Hayakawa, [7] F. Sélley, A.´ Besenyei, I. Z. Kiss, and P. L. Simon, Phys. Rev. E 78, 016103 (2008). SIAM J. Appl. Dyn. Syst. 14, 168 (2015). [25] C. Nardini, B. Kozma, and A. Barrat, [8] T. Gross, C. J. Dommar D’Lima, and B. Blasius, Phys. Rev. Lett. 100, 158701 (2008). Phys. Rev. Lett. 96, 208701 (2006). [26] M. Pósfai and P. Hövel, New J. Phys. 16, 123055 (2014). [9] R. Pastor-Satorras and A. Vespignani, [27] Y. Pan and X. Li, PLoS One 9, e94998 (2014). Phys. Rev. Lett. 86, 3200 (2001). [28] R. D. Enos and A. Fowler, [10] R. May and A. Lloyd, Phys. Rev. E 64, 066112 (2001). Polit. Sci. Res. Methods , doi:10.1017/psrm.2016.21 (in press). [11] Y.-Y. Liu, J.-J. Slotine, and A.-L. Barabási, [29] M. Wiedermann, J. F. Donges, J. Heitzig, W. Lucht, and Nature 473, 167 (2011). J. Kurths, Phys. Rev. E 91, 052801 (2015). [12] Z. Yuan, C. Zhao, Z. Di, W.-X. Wang, and Y.-C. Lai, [30] H. Silk, G. Demirel, M. Homer, and T. Gross, Nat. Commun. 4, 2447 (2013). New J. Phys. 16, 093051 (2014). [13] N. J. Cowan, E. J. Chastain, D. A. Vilhena, [31] G. A. Böhme and T. Gross, J. S. Freudenberg, and C. T. Bergstrom, Phys. Rev. E 83, 035101 (2011). PLoS One 7, e38398 (2012). [32] G. Demirel, F. Vazquez, G. A. Böhme, and T. Gross, [14] J. C. Nacher and T. Akutsu, Physica D 267, 68 (2014). New J. Phys. 14, 073005 (2012). [33] P. Erdös and A. Rényi, Publ. Math. 6, 290 (1959). [15] S. Galam and F. Jacobs, Physica A 381, 366 (2007). [34] M. E. J. Newman, Phys. Rev. E 67, 026126 (2003), [16] E. Yildiz, D. Acemoglu, A. Ozdaglar, A. Saberi, and 0209450 [cond-mat]. A. Scaglione, ACM Trans. Econ. Comput. 1, 19 (2013). [35] J. P. Gleeson, S. Melnik, J. A. Ward, M. A. Porter, and [17] F. Sorrentino, M. di Bernardo, F. Garofalo, and G. Chen, P. J. Mucha, Phys. Rev. E 85, 026106 (2012). Phys. Rev. E 75, 046103 (2007). [36] J. Sun and A. E. Motter, [18] J. Xie, S. Sreenivasan, G. Korniss, W. Zhang, C. Lim, Phys. Rev. Lett. 110, 208701 (2013). and B. K. Szymanski, Phys. Rev. E 84, 011130 (2011). [37] P. Holme and G. Ghoshal, in Adaptive Networks, edited [19] F. Wu and B. Huberman, (2004), by T. Gross and H. Sayama (Springer, Dordrecht, 2009) arXiv:cond-mat/0407252. pp. 269–288. [20] P. Holme, Eur. Phys. J. B 88, 234 (2015). [38] W. Barfuss, J. F. Donges, M. Wiedermann, and [21] P. Holme and M. E. J. Newman, W. Lucht, Earth Syst. Dyn. 8, 255 (2017). Phys. Rev. E 74, 056108 (2006).