Emergence of complex structures from nonlinear interactions and noise in coevolving networks

We study the joint effect of the non-linearity of interactions and noise on coevolutionary dynamics. We choose the coevolving voter model as a prototype framework for this problem. By numerical simulations and analytical approximations we find three main phases that differ in the absolute magnetization and the size of the largest component: a consensus phase, a coexistence phase, and a dynamical fragmentation phase. More detailed analysis reveals inner differences in these phases, allowing us to divide two of them further. In the consensus phase we can distinguish between a weak or alternating consensus (switching between two opposite consensus states), and a strong consensus, in which the system remains in the same state for the whole realization of the stochastic dynamics. Additionally, weak and strong consensus phases scale differently with the system size. The strong consensus phase exists for superlinear interactions and it is the only consensus phase that survives in the thermodynamic limit. In the coexistence phase we distinguish a fully-mixing phase (both states well mixed in the network) and a structured coexistence phase, where the number of links connecting nodes in different states (active links) drops significantly due to the formation of two homogeneous communities of opposite states connected by a few links. The structured coexistence phase is an example of emergence of community structure from not exclusively topological dynamics, but coevolution. Our numerical observations are supported by an analytical description using a pair approximation approach and an ad-hoc calculation for the transition between the coexistence and dynamical fragmentation phases. Our work shows how simple interaction rules including the joint effect of non-linearity, noise, and coevolution lead to complex structures relevant in the description of social systems.


Introduction
Coevolving or adaptive network models 1 provide a better representation of real-world systems in comparison with static or evolving networks . Most empirical networks display both network dynamics (evolution of the network's topology) as well as dynamics of the state of the nodes 2, 3 . Moreover, a nontrivial feedback loop between these aspects renders a simple sum of effects analyzed separately incomplete. Adaptive mechanisms coupling network and nodes state dynamics give rise to new phenomena absent when coevolution process is not taken into account [4][5][6][7][8][9][10][11][12][13][14] . Coevolution models incorporate microscopic assumptions in better agreement with empirical observations, and they also produce new macroscopic results.
Another essential feature of many real-world systems is the non-linearity associated with non-dyadic interactions. It is often assumed in network models that an interaction occurs pairwise, only between two selected vertices. From a single node point of view it means selecting one of its neighbors at random for the interaction. This leads to a linear relation between the number of neighbors in a given state and the probability of choosing one of them. However, in non-dyadic or group interactions, linearity is lost 15 . In contagion or spreading processes, the difference between these two types of interaction goes under the name of simple vs. complex contagion [16][17][18] A third crucial empirical element in many dynamical processes on networks is noise. This is specially important in social systems where noise is inevitable 19,20 . It can manifest itself on various levels. First, people chose other people to interact with at random. The exact form of this randomness can take different forms, nevertheless the structure's evolution is never hard-coded. But the most fundamental part of randomness lays probably within individual choices. For example, having exactly the same influence on two people's opinions we can not be sure of the outcome. This mechanism is sometimes referred to as non-conformism 21 . It reflects the ability of agents to change state independently of the states of their neighbors. It is often a model parameter that needs to be calibrated to reproduce empirical data 22 .
In this paper we aim at exploring the joint effect of these three important aspects -the coevolution of network structure and  Figure 1. Schematic illustration of update rules in the nonlinear coevolving voter model with noise. Every node is in a state +1 or −1, indicated by orange and blue colors. The active node i is chosen randomly. Then with probability (a i /k i ) q an interaction occurs and one of the active links (the one to the node j) is selected randomly. With probability p the link (i, j) is rewired to the link (i, k), where k is a random node being in the same state as i. With probability 1 − p the focal node i copies the state of the node j. At the end of the time step, regardless of what happened before, the active node draws a random state with probability ε.
node states, the non-linearity of interactions and the noise -on the behavior of the system. As the framework we choose the simple voter model 23,24 . With a binary state it is often used as a model of opinion dynamics, but different forms and extensions of the model have been fruitful in explaining empirical observations in fairly distinct phenomena such as electoral processes 22 , stock market 25 or online communities 26 . The consequences of coevolution [27][28][29] , noise [30][31][32][33] , and non-linearity 15 have been already considered separately in the voter model. The joint effect of these aspects, however, turns out to be more complex than a mere superposition of the results obtained so far. The coevolving voter model (CVM) 27 was among the pioneers in introducing adaptive mechanisms in general. In the standard voter model, node state dynamics follows an imitation rule which is here coupled with link rewiring, introducing the coevolution. This leads to a network fragmentation transition. The effect of noise in the CVM 34 prevents the existence of absorbing configurations so that the different phases of the system are described by dynamically active stationary states. These include a striking new dynamical fragmentation phase. Additionally, a fully-mixing phase is found, as could be expected for large noise levels. Nonlinear interactions have been also considered in the CVM 35,36 . Non-linearity changes the stability of fixed points in the voter model dynamics, leading to a new dynamically trapped coexistence phase. Finally, the joint effects of noise and non-linearity in the voter model have also been considered 37 . It was found that non-linearities transform a finite size transition known from the noisy voter model into a bona fide phase transition that survives in the thermodynamic limit. Here we introduce a CVM in which noise and non-linearities are jointly taken into account.

The model
Our noisy and nonlinear CVM is defined by its time evolution over discrete time steps as follows. First a random graph is generated and every node is assigned a state s i ∈ {−1, +1} at random 1 . In each time step a node i is chosen at random, we call it the active or focal node. Then, with probability (a i /k i ) q ≡ ρ q i an interaction occurs, where k i is the degree of the focal node i, a i is the number of neighbors of the node i being in the opposite state, and q is the non-linearity parameter of the model. If an interaction occurs, one of the a i neighbors in a different state is chosen randomly, call it j. Then, with probability p a link rewiring is performed and with complementary probability 1 − p a state copying. When rewiring, the node i cuts the connection to the node j and creates a new link with a randomly selected node in the same state (if there is no such node, nothing happens). When copying the state, the node i replicates the state of the node j, i.e. s i → s i = s j . At the and of the time step, regardless of what happened before, the active node with probability ε draws a random state. Note that this is equivalent to flipping the current state with probability ε/2. The algorithm of the model is illustrated in Figure 1. Our model has three parameters, namely the noise rate ε, the plasticity p and the non-linearity q. The parameter p is a network plasticity parameter measuring the ratio of time scales of node dynamics and network dynamics. The non-linearity parameter q measures the nonlinear effect of local majorities: q = 1 corresponds to the ordinary voter model with a mechanism of random imitation, while q < 1 indicates a probability of imitation above random imitation and q > 1 a probability below random imitation. The ordinary voter model corresponds to p = ε = 0 and q = 1, while p = ε = 0 corresponds to the nonlinear voter model and p = 0 and q = 1 to the noisy voter model. The CVM is obtained for ε = 0 and q = 1, the noisy CVM is obtained for q = 1 and the nonlinear CVM for ε = 0.
Our simulations are ran from an initial random network with N nodes and M links or average degree µ = ∑ i k i /N = 2M/N, and with random initial conditions for the nodes states until a stationary state or a frozen configuration is reached (the latter being possible only for ε = 0).

Results
We explore the space of possible values of the parameters (p, q, ε) by means of computer simulations and analytical approximations. In order to describe the system, we adapt typical order parameters, such as magnetization m = ∑ i s i /N, size of the largest network component S, and density of active links ρ = ∑ i a i /2M. By an active link we mean a connection between two nodes being in opposite states. All these values are usually normalized to fit the range [0, 1]. Additionally, we define a new indicator, namely community overlap (ov.), in order to be able to distinguish structural changes. The community overlap is a fraction of nodes assigned to the same community by both their state and by the algorithm of community detection in the network 38 . Consequently, if a given node was assigned to the same community in both cases it increases the ov. by 1/N (where the denominator comes from normalization). Table 1. Average values of the order parameters in different phases. The border between phases A and B is given by |m| = 0.5, between phases B1 and B2 is given by ov. = 0.75 which is approximated by ρ = 0.1, the border between phases B and C is given by S = 0.75. The difference between phases A1 and A2 can be observed in the dynamical behavior of the magnetization and on the system size scaling (see Figure 2 and 7).

Phase diagram
We numerically study the p-ε phase diagram for three different values of the q parameter -the sublinear case q = 0.5, the ordinary linear case q = 1, and the superlinear case q = 2. These phase diagrams are shown in Figure 2 for two different network sizes. Obviously, for any finite amount of noise in the system a frozen configuration does not exist, and any phase is described by a characteristic dynamical stationary state. We can distinguish three general phases in the model. Phase A, indicated by the red area in the figure, is a consensus phase. In this range of parameters the system stays in a consensus state for most of the time, i.e. magnetization is close to ±1 and the network is connected having a single large component S = 1 and a small number of active links ρ 0 . If we increase the noise rate ε or the plasticity p sufficiently, we obtain phase B indicated by the white area in Figure 2, and referred to as a fully-mixing phase in previous work for the noisy CVM 34 . In this phase the magnetization drops to zero, m = 0, hence there is no consensus in the system any more. In addition, the network stays connected most of the time, S 1 . We refer to this phase as a coexistence phase. As we will see, phase B is not homogeneous in its whole range of parameters and it can be either fully-mixing (phase B1) or structured (phase B2), what is indicated by different values of the density of active links ρ and community overlap ov. Finally, for values of the rewiring probability above  Table 1).
the critical point p c of the nonlinear CVM 35 , and relatively small noise rates, phase C arises. It is marked by the blue area in the figure. In this region we find dynamical fragmentation -the network consists of two separate components with nodes in each one being in opposite states, so that m 0, S 0.5 and ρ = 0. It is possible, however, that the two network components get connected intermittently in the stationary state due to noise and random rewiring, creating again a single component network with m 0 and a small number of active links ρ 0. Phase C can be described as dynamical switching between these two arrangements. Values of all analyzed quantities for every phase are summarized in Table 1.
For the linear case (q = 1) phases A and C exist only for a finite size of the network, and the size of these phases in the parameter space decreases with growing number of nodes. For the sublinear scenario q < 1, we can see in Figure 2 that the same holds for phase C, while phase A does not exist at all. The only point where the average absolute magnetization slightly raises is at p c and for ε ≈ 0, but its maximal value is only about 0.3. This raise is due to higher fluctuations close to the transition point. On the other hand, phase C prevails for even twice larger noise rate than in the linear scenario.
In the superlinear case (q = 2) phase C is much smaller in parameter space and disappears faster with growing system size. But phase A prevails for much larger noise than in the linear case q = 1. We observe phase A even for ε larger by almost two orders of magnitude. Additionally, the system size scaling is different for the superlinear scenario. Indeed, non-linearity has significant influence on the nature of phase A. For q < 1 it does not exist, for q = 1 it exists only in finite networks, and for q > 1 phase A persists in the thermodynamic limit.
A closer look at the phases In order to better understand the behavior of the system and the differences between phases we analyse horizontal cross-sections of Figure 2 and single-run trajectories which are presented in Figures 3-6 for the linear, sublinear and superlinear cases. In panel (a) of each of them, a phase diagram with respect to rewiring probability p and for a particular value of ε is presented. Values of the noise rate are chosen in such a way that allows to show three phases in one panel. For an horizontal cross-section of the full phase diagram it is difficult to capture all phases. Therefore, areas of the middle phase can be narrow, but still different values of the order parameters can be distinguished. Panels (b)-(d) show typical time traces of the order parameters in different phases.
For the sublinear case (q = 0.5) we can see the differences between phase B1(fully-mixing, Figure 3b) and phase B2 (structured, Figure 3c). Phases B1 and B2 have zero average absolute magnetization |m| ≈ 0, but we can distinguish a region  with high density of active links (B1) and small density of active links (B2). Phase C also has a low density of active links, but the largest network component switches from S = 1 to S = 0.5, giving and average value S ≈ 0.5, whereas phase B2 on the average stays connected.
Results for the linear case (q = 1) are shown in Figure 4. Phase A is characterized by a magnetization which tends to stay at one of the consensus states, but it can switch from −1 to +1, or the other way around, during the time evolution. Therefore, |m| ≈ 1 but m = 0. To distinguish it from the superlinear case where m ≈ ±1 we call this phase A1. We also observe that for q = 1 in phase B2 ( Figure 4c) the network can fragment close to the transition line, however it remains a single component network most of the time. Figure 5 and Figure 6 correspond, respectively, to small and large noise rates in the superlinear case (q = 2). In this scenario the consensus phase A prevails for much larger noise rate. In panels (b) of Figure 5 and Figure 6 we can see how the system behaves in phase A for q = 2. It quickly reaches a consensus state for either m = 1 or m = −1 and remains at this value of magnetization. Therefore, |m| ≈ 1 and in contrast to the linear case in a time average also m ≈ ±1. To account for this difference we call the consensus phase A2 for q > 1 .
The difference between the consensus phase A1 and A2 is also clearly visible from the probability distribution of the magnetization in a given realization of the dynamical process ( Figure 7). In the linear case there is a bimodal distribution for the magnetization with two equal peaks at values +1 and −1 (Figure 7d), while for the superlinear case there is a single peak for a value of the magnetization at either of the boundary values +1 or −1, depending on the run (Figure 7g). For q = 2, once the consensus is reached the system stays there with minor fluctuations (phase A2), while for q = 1 the system goes back and forth between opposite consensus states (phase A1). Furthermore, phase A2 is robust against finite-size fluctuations, while phase A1 disappears in the thermodynamic limit 34 (see Figure 2).
The distribution of the magnetization gives additional insights on the phase diagram: The fact that phase A does not exist in the sublinear case (q < 1), is reflected in a distribution with a single peak at 0 for all values of ε (Figure 7a-c). However, the variance of the distribution takes its maximal value for noise going to zero and p = p c , i.e. close to the transition point between coexistence and fragmentation phase in the nonlinear CVM 35 . A different form of the transition between phases A and B for the linear and superlinear case is also observed. For q = 1 there is a flat distribution at the transition point (Figure 7e) , while a trimodal distribution is found for q = 2 (Figure 7h). A trimodal magnetization distribution was reported before in the noisy voter model on a static network 37 , but only for non-linearity parameter equal 5 or larger. With coevolution, trimodality here is obtained already for q = 2.
Phase C can be defined in terms of the size of the largest component. In phases A and B it is equal to the size of the whole network (S = 1), while the phase C is characterized by a dynamical fragmentation into two components of similar size and opposite state. Due to noise expressed in random changes of nodes states and rewiring the components are constantly being reconnected and disconnected. It can be examined looking at the trajectory or at the probability distribution of the size of the largest component, which is presented in Figure 8.

Community structure
Phase B is generally defined by zero average magnetization, also zero average absolute magnetization, and by the existence of one large network component. Nonetheless, this description leaves room for different possible configurations. Analysis of the trajectories showed that the density of active links can vary within phase B, but the question is weather this is a sign of a topological change. In the linear case (q = 1) only the fully-mixing phase was reported 34 , with nodes of states +1 and −1 well mixed inside a random graph. We refer to this configuration as phase B1. On the other hand, we can satisfy conditions for the phase B having two evident communities, highly connected internally and of opposite states, with only a few links bridging them. There is still zero magnetization and one large network component in such configuration characterized by a small number of active links. We call this phase B2. The difference between phases B1 and B2 is clearly seen in Figure 9.
Although the difference between phases B1 and B2 can be seen in the density of active links, a closer look at Figure 9 suggests that phase B2 has well defined topological communities. Therefore, we propose an alternative quantitative measure for the difference between phases B1 and B2, the overlap between state communities -defined by the state of the nodes -and structural communities found by a community detection algorithm. We use a classical algorithm from 38 , but the result does not differ much when using other algorithms. Each node is assigned to the state community by its state and to a structural community by the algorithm's result. The relative overlap between these two communities is a new quantitative indicator of the phase of the system 2 . For a random assignment or no community structure the overlap will be close to 0.5. This is the situation in phase B1. For phase B2 the overlap will be close to 1. This means that our dynamical coevolving model generates clear topological communities emerging from local interactions involving only state of nodes. This result may potentially explain process of formation of communities in social networks, where such structures are especially common 39 .

Identifying phase boundaries
So far, we gave a description of different phases with different qualitative behavior. Transition lines between these different types of behavior are not clearly or unambiguously defined because every analyzed phase indicator changes value significantly across the phase diagram. We do not focus on properties of phase transitions in this work, but rather on the properties of different emerging structures. Therefore we follow here a simple and pragmatic way to identify boundaries between different phases, employing previous approaches 34 . These boundaries have to be understood as an arbitrary way to deal with crossover system behavior.
Each of the order parameters or phase indicators analyzed has a continuous range of values (for a large N) and different phases are described by the extreme values of these quantities, allowed in this range. Therefore, a straightforward way of dividing the phase diagram into separate phases is to use the middle value, i.e. the value in the middle between the maximum and minimum of a given range. For example, the absolute magnetization is defined in the range [0, 1], taking a value close to 1 in phase A and a value close to zero in phase B. Hence, we identify the border between the two phases with |m| = 0.5.
Likewise, the size of the normalized largest network component takes values in the range [0.5, 1]. In phases A and B there is a single component network (S = 1), while phase C is characterized by the dynamical fragmentation into two components of similar size and opposite state, so that S = 0.5. In this phase, due to noise expressed in random changes of nodes states and rewiring, the components are constantly being reconnected and disconnected. We then identify 34 the border between phase B and phase C by the middle value S = 0.75. This is a line at which the network is half of the time fragmented and half of the time contains only one big component. This phase boundary can also be obtained from the probability distribution of the size of the largest component (Figure 8). At the transition line (panel b) two peaks have the same area.
Finally the boundary between phases B1 and B2 can be identified in terms of the the overlap (ov.), which takes values in the range [0.5, 1]. There is a transition from values around 0.5 in B1, up to 1 in B2, and so we define the transition line at ov. = 0.75. However, this parameter is computationally very demanding. Alternatively, we can approximate the identification of the transition line by a small value of the density of active links which we arbitrarily fix at ρ = 0.1

Analytical predictions
The magnetization m and the density of active links ρ obey Equations 3 describing the dynamics of the system, as derived in the Methods section. Several fixed points (m * , ρ * ) of these ordinary differential equations can be found depending on parameter values. However, not all of them are stable, therefore not all of them are observed in numerical simulations. To analyze the stability of these fixed points, we consider flow diagrams of the dynamics in Figure 10. Note, that from panel (a) to (b) only the value of q changes, emphasizing the difference between sublinear and superlinear cases, while from panel (b) to (c) only the value of ε changes, emphasizing the noise effect in the superlinear case. A change in the non-linearity parameter q can reverse the stability of fixed points when going over the boundary value of 1. A change in the noise rate ε can additionally shift the position of fixed points allowing for fixed points different than m = −1, 0, 1. Since the analytical description is derived in the thermodynamic limit, we don't observe stable fixed points at non-zero magnetization for q ≤ 1. This finding is consistent with the scaling behavior of numerical results indicating existence only of phase B in the limit of a large number of nodes N. For the superlinear case, although the fixed points are placed at the same values of the magnetization m = −1, 0, 1 (for p < p c , ε < ε c ), their stability is inverted -now only the solutions of |m| = 1 are stable, corresponding to phase A2. This is clearly seen in the analytical prediction for the phase diagram in Figure 11. These results, obtained in the thermodynamic limit, indicate that phase A2 should be observed for any N when q = 2, which is in agreement with our numerical results in Figure 2. Separate mean-field prediction of the disappearance of phase A1 in the thermodynamic limit was given by Diakonova et. al. 34 for the linear CVM with noise. Additionally, phase C is not obtained in the thermodynamic limit. The non-linear noisy voter model in a fixed network (p = 0) has been thoroughly analytically studied 37 , showing that q = 1 is a bordering value between a unimodal and bimodal distribution of the magnetization m. In other words, it is a transition line between existence and nonexistence of phase A. The agreement of our results with previous studies can be seen when analyzing the extreme value of p = 0 in the phase diagrams (Figure 11a and b). It is also separately presented in the Figure 11c. The transition to phase A, characterized by nonzero absolute magnetization |m|, exists for finite values of ε only when q > 1. In the Methods section we derive a formula for the critical value of the noise rate ε c (p = 0) at which the system looses the consensus state, that is, the transition form phase A to phase B (Equation 7). This result is in agreement with the numerical solution of Equations 3 presented in Figure 11c, giving a critical value of the density of active links ρ c = 1 3 ≈ 0.33 and ε c = 2 11 ≈ 0.18 for the parameters values used in the figure. The analytical solution from Equation 7 also predicts disappearance of the transition at q = 1. In the reference 37 a similar prediction of the critical noise rate was given for a complete graph: ε c (p = 0) = 2 −q (q − 1) which would give for the case of Figure 11c ε c = 0.25. Therefore, a complete graph gives only a first approximation to the value found here.
There has been no previous attempt of an analytical approximation describing the transition from phase B to the dynamical fragmentation phase C already found in the noisy linear CVM. We propose here a simple description of this phenomena. Having two separate, but internally homogeneous network components (clusters) with nodes in opposite states, the only way of connecting them is by a random change of node's state and link rewiring to the second component. The probability of the first event is independent of q and is simply given by ε/2. When the node changing state is selected as the active node the probability of an interaction is ρ q i . Since ρ i ∈ [0, 1], for smaller q the probability of an interaction is higher, except for boundary cases with ρ i = 0, 1. To reconnect the two clusters, rewiring must occur, but this happens always with probability p, despite the value of q. Therefore, for a single node in a state opposite to the whole cluster, the probability of connecting to the other cluster is constant (since ρ i = 1). However, once the two clusters are connected, the higher probability of an interaction for lower q means a higher probability of rewiring causing fragmentation again. Consequently, we expect phase C to persist for larger noise when q is smaller. More detailed description of this process is given in the Methods section, where we derive the following formula for the transition line: Based on this approximation we predict phase C to fade with growing non-linearity parameter q or with growing system size N, as shown in Figure 12. Both predictions are consistent with our numerical results.

Discussion
In this paper we analyzed the nonlinear coevolving voter model with noise. Depending on the values of the three main parameters -the rewiring probability p, the noise intensity ε, and the non-linearity parameter q -we observed three distinct phases: a consensus phase A, a coexistence phase B, and a dynamical fragmentation phase C. We observed, however, significant internal differences within phases A and B. The first one can be further divided into phase A1 and A2. Phase A1, for q = 1, is a consensus phase with absolute magnetization equal 1 on average, but real magnetization switching between −1 and +1 states, giving rise to a bimodal magnetization distribution within one realization of the stochastic dynamics. In phase A2, observed for q > 1, there is a stable consensus, i.e. global magnetization states −1 and +1 are stable. Consequently, during one realization of the stochastic dynamics the system remains in a given consensus state, producing a unimodal magnetization distribution with a peak at the maximal (minimal) magnetization. Additionally, phases A1 and A2 have different system size scaling -phase A1 disappears in the thermodynamic limit, while phase A2 is stable against finite size fluctuations. Finally, phase A does not exist for q < 1. Phase B can be similarly divided into phases B1 and B2. Phase B1 is a fully-mixing phase with random network structure and random states of the nodes, giving zero magnetization. But for larger p and low noise intensity we observe phase B2, which has the same vanishing average magnetization, but a different network structure. In the structured phase B2 one can easily distinguish two communities of opposite states connected by just a few inter-community links. This structural difference is confirmed by community detection algorithms.
Phase C is associated with a dynamical fragmentation of the network -two components in opposite magnetization states are being constantly connected and disconnected. We derive an analytical description of this behavior and an approximated value for the transition line between phase B2 and phase C.
The only phases surviving in the thermodynamic limit are phases A2, B1 and B2. The transition line form phase A2 to phase B is largely independent of finite size fluctuations. We have also presented an analytical pair approximation able to describe these findings and the main features of phases A and B in the thermodynamic limit.
Our work fills a gap in the studies of the CVM. It provides a binding bridge between studies on the CVM with noise 34 and studies on the nonlinear CVM 35 . It reduces to the nonlinear noisy voter model 37 and the ordinary CVM 27 for a proper configuration of parameters values. We obtain full consistency with those limiting cases and we explore new parameter domains. Our work brings the analysis of the voter model to a greater complexity by taking into account the joint effect of noise, coevolution and non-linearity which turns out not to be a mere superposition of them. It may provide a tool for the evaluation of the relevance of different mechanisms in the description of opinion dynamics, but can be also a reference point in the study of coevolving network models. We also show how nonlinear vs linear interactions can change the stability of a consensus state in the network and how topological communities can arise from non-topological interactions. These results are of relevant value in the description of social networks.

Pair approximation
We use the same approach as used for the nonlinear coevolving voter model 35 to describe the dynamics of magnetization m and the density of active links ρ. Given the network homogeneity due to the random rewiring, we assume each node to have the same average degree µ = 2M/N. Let us denote by n + = (1 + m)/2 and n − = (1 − m)/2 the fraction of nodes in the state +1 and −1 respectively. Then, when we pick a node in the state ±1 as the active node, the probability of choosing a neighbor in the opposite state is given by ρ/2n ± . In other words, ρ 2n ± gives the density of active links ρ i for a node i being in the state s i = ±1. Therefore, the probability of an interaction is given by ρ q i = (ρ/2n ± ) q ≡ n ± q . When q takes integer values this can be also interpreted as the probability of choosing a neighbor in the opposite state q times. Hence, when an interaction occurs with this probability, we can make the approximation that there is at least q neighbors in the opposite state and for the rest of them the probability of being in a different state than the focal node is ρ/2n ± . All together this implies that 3 a i ≈ q + (µ − q) ρ 2n ± . To approximate the evolution of the density of active links ρ we must estimate the contributions of different events that can result in a change of ρ. These events are: (i) rewiring, followed by a change of state through noise, (ii) rewiring, without a change of the node's state due to noise, (iii) changing the state of the node through state copying with no further change due to noise, and (iv) changing the state of the node only as a result of noise, with no previous state copying or rewiring. Let δ ± be the change in the total number of active links given that a node i such that s i = ±1 changed state. The total change in the number of active links in the four possible events above is: (i) 1 + δ ± , (ii) −1, and for events (iii) and (iv) just δ ± . Magnetization changes only when the state of the node is changed via copying or as a result of noise in three possible scenarios: state copying with no noise effect, link rewiring followed by noise effect, and no interaction -neither state copying nor link rewiring -but noise acting alone. When the focal node having a state s i = ±1 changes its state, the total change in the magnetization is ∓2. Hence, in the thermodynamic limit we have: With a few simple algebraic transformations these equations can be rewritten as: When node i changes state, all its a i active links become inactive and all other µ − a i inactive links become active, therefore the total change in the number of active links is δ ± = µ − 2a i . Using the previous approximation for a i we can write δ ± = µ − 2q − 2(µ − q) ρ 2n ± . The simplest stationary solution (fixed point) of Equations 3 is given by taking the magnetization m = 0, which leads to an equation for the stationary value of ρ: A fixed point solution with m = ±1 does not exist for any finite noise rate ε, whilst for ε = 0 a stationary solution is ρ = 0. Setting the noise rate to zero together with m = 0 we obtain the stationary solution of the nonlinear CVM 35 : For ε = 0 and q = 1 we recover the solution of the standard CVM 27 : To compare our results with the nonlinear noisy voter model on static networks we analyze our approximation for the particular case p = 0. Putting ρ = ρ c (1 − m 2 ) in the first of Equations 3 and performing a stability analysis of the fixed point solution m = 0 we can find a critical noise value which depends on the critical value of the density of active links ρ c . The latter can be obtained from Equation 4 as ρ c = µ−2q 2µ−2q . This formula is shown to be in full agreement with numerical solutions of Equations 3 shown in Figure 11c.

Phase C: Finite size scaling
In order to describe the behavior of the dynamical fragmentation phase we first look for an approximation for the probabilities of reconnecting two separate clusters and of disconnecting two clusters sharing at most two links. In this approach we omit events of probability proportional to (1/N) 3 and to ε 2 , or of higher order.
Imagine two separate and internally homogeneous components of opposite states, as it happens in phase C. The simplest way of connecting them under the rules of the nonlinear noisy CVM involves two steps. First, one of the nodes, call it i, must change it's state. This is possible only due to noise and occurs with probability ε/2. Second, node i that has changed its state must rewire one of its links to a node in the opposite cluster, having the same state. This can happen with probability pρ q i /N, because we need to select this particular node as the active node (1/N), an interaction has to occur (ρ q i ), and rewiring must be performed (p). Note, that since node i is the only node in a different state than its cluster ρ q i = 1. Finally, it gives the probability of reconnecting two components equal: Approaching the transition line between phases B and C now from phase B, so that a fragmentation event occurs, we consider one single component network disconnecting into two equal clusters. As done before, imagine a situation two time steps before a possible fragmentation -network has two internally homogeneous components in opposite. One of the nodes i is part of a bridge, i.e. it is connected to two nodes in the opposite cluster. Now, for the fragmentation to occur we need both of the links between the components to be rewired. The probability to rewire the first one is 1 N (2/µ) q p(1 − ε 2 ) + 2 N (1/µ) q p. We have to select the node i (1/N) or one of its neighbors (2/N). An interaction must occur, what happens with probability (a j /µ) q , where the number of active links is 2 for node i and 1 for each of its neighbors in the opposite cluster. Finally, a rewiring must occur with probability p. Additionally if node i was selected, it can not change its state due to noise (1 − ε 2 ), otherwise fragmentation could not be achieved in two steps. The transition occurs, however, for very small values of noise and therefore we can approximate 1 − ε 2 ≈ 1. To rewire the second link we have to select one of the two nodes (2/N) connecting the link, an interaction must occur (1/µ q ), which must be a rewiring event (p). Therefore, the probability of losing the last link between the two clusters is 2 N (1/µ) q p. Finally, we obtain the probability of disconnecting two clusters sharing only two links: Between phases B and C a continuous fragmenting and reconnecting of the network is observed. We define the transition between the two phases when connection and fragmentation happens at such a rate that half of the time the system consists of two separate components and half of the time the network is connected. Therefore, at the transition line we expect P r = P d , which leads to the equation for the critical density of noise given in the main text (Equation 1).