Fundamental ingredients for discontinuous phase transitions in the inertial majority vote model.

Discontinuous transitions have received considerable interest due to the uncovering that many phenomena such as catastrophic changes, epidemic outbreaks and synchronization present a behavior signed by abrupt (macroscopic) changes (instead of smooth ones) as a tuning parameter is changed. However, in different cases there are still scarce microscopic models reproducing such above trademarks. With these ideas in mind, we investigate the key ingredients underpinning the discontinuous transition in one of the simplest systems with up-down Z2 symmetry recently ascertained in [Phys. Rev. E 95, 042304 (2017)]. Such system, in the presence of an extra ingredient-the inertia- has its continuous transition being switched to a discontinuous one in complex networks. We scrutinize the role of three central ingredients: inertia, system degree, and the lattice topology. Our analysis has been carried out for regular lattices and random regular networks with different node degrees (interacting neighborhood) through mean-field theory (MFT) treatment and numerical simulations. Our findings reveal that not only the inertia but also the connectivity constitute essential elements for shifting the phase transition. Astoundingly, they also manifest in low-dimensional regular topologies, exposing a scaling behavior entirely different than those from the complex networks case. Therefore, our findings put on firmer bases the essential issues for the manifestation of discontinuous transitions in such relevant class of systems with Z2 symmetry.

Introduction

Spontaneous breaking symmetry manifests in a countless sort of systems besides the classical ferromagnetic-paramagnetic phase transition[1,2]. For example, fishes moving in ordered schools, as a strategy of protecting themselves against predators, can suddenly reverse the direction of their motion due to the emergence of some external factor, such as water turbulence, or opacity[3]. Also, some species of Asian fireflies start (at night) emitting unsynchronized flashes of light but, some time later, the whole swarm is flashing in a coherent way[4]. In social systems as well, order-disorder transitions describe the spontaneous formation of a common language, culture or the emergence of consensus[5].

Systems with Z2 (“up-down”) symmetry constitute ubiquitous models of spontaneous breaking symmetry, and their phase transitions and universality classes have been an active topic of research during the last decades[1,2,6]. Nonetheless, several transitions between the distinct regimes do not follow smooth behaviors[7-9], but instead, they manifest through abrupt shifts. These discontinuous (nonequilibrium) transitions have received much less attention than the critical transitions and a complete understanding of their essential aspects is still lacking. In some system classes, essential mechanisms for their occurrence[10], competition with distinct dynamics[11,12], phenomenological finite-size theory[13] and others[14-17] have been pinpointed.

Heuristically, the occurrence of a continuous transition in systems with Z2 symmetry is described (at a mean field level) by the logistic equation , that exhibit the steady solutions m = 0 and . The first solution is stable for negative values of the tuning parameter a, while the second is stable for positive values of a. For the description of abrupt shifts, on the other hand, one requires the inclusion of an additional term + cm5, where c > 0 ensures finite values of m. In such case, the jump of m yields at , reading . Despite portrayed under the simple above logistic equation, there are scarce (nonequilibrium) microscopic models forecasting discontinuous transitions[18].

Recently, Chen et al.[19] showed that the usual majority vote (MV) model, an emblematic example of nonequilibrium system with Z2 symmetry[20-22], exhibits a discontinuous transition in complex networks, provided relevant strengths of inertia (dependence on the local spin) is incorporated in the dynamics. This results in a stark contrast with the original (non-inertial) MV, whose phase transition is second-order, irrespective the lattice topology and neighborhood. The importance of such results is highlighted by the fact that behavioral inertia is an essential characteristic of human being and animal groups. Therefore, inertia can be a significant ingredient triggering abrupt transitions that arise in social systems[5].

Although inertia plays a central role for changing the nature of the phase transition, their effects allied to other components have not been satisfactorily understood yet[23,24]. More concretely, does the phase transition become discontinuous irrespective of the neighborhood or on the contrary, is it required a minimal neighborhood for (additionally to the inertia) promoting a discontinuous shift? Another important question concerns the topology of the network. Is it a principal ingredient? Do complex and low-dimensional regular structures bring us similar conclusions?

Aimed at addressing questions mentioned above, here we examine separately, the role of three key ingredients: inertia, system degree, and the lattice topology. For instance, we consider regular lattice and random regular (RR) networks for different system degrees through mean-field theory (MFT) treatment and numerical simulations. Our findings point out that a minimal neighborhood is also an essential element for promoting an abrupt transition. Astonishing, a discontinuous transition is also observed in low-dimensional regular structures, whose scaling behavior is wholly different from that presented in complex networks[13]. Therefore, our upshots put on firmer bases the minimum and essential issues for the manifestation of “up-down” discontinuous transitions.

Model and Results

In the original MV, with probability 1 − f each node i tends to align itself with its local neighborhood majority and, with complementary probability f, the majority rule is not followed. By increasing the misalignment parameter f, a continuous order-disorder phase transition takes place, irrespective the lattice topology[20-22]. Chen et al.[19] included in the original model a term proportional to the local spin σ, with strength θ, given bywhere sign(X) = ±1 and 0, according to X > 0, <0 and X = 0, respectively. Note that one recovers the original rules as θ = 0.

MFT results

In several cases, a mean field treatment affords a good description of the model properties. By following the main steps from refs[19,21,24,25], we derive relations for evaluating the order parameter m for fixed f, θ and k [see Methods, Eqs (3–8)]. Figure 1 shows the main results for k = 4, 8 and 12. Note that MFT predicts a continuous phase transition for k = 4 irrespective the value of θ [see panels (a) and (b)], in which m is a decreasing monotonic function of the misalignment parameter f. An opposite scenario is drawn for k = 8 and 12, where phase coexistence stems as θ increases [see panels (c–f)]. They are signed by the presence of a spinodal curve, emerging at f [see, e.g., panel (c) and (e)] and meeting the monotonic decreasing branch at f. For k = 8, the coexistence line arises only when θ > 1/3 and is very tiny (f − f is about 2.10−4), but they are more pronounced for θ > 3/7. Analogous phase coexistence hallmarks also appear for k = 12 (panel (e)) and k = 20 (Fig. 6 and ref.[19]). Thus, MFT insights us that large θ and k (k ≥ 6) are core ingredients for the appearance of a discontinuous phase transition. A remarkable feature concerning the phase diagrams is the existence of plateaus, in which the transition points present identical values within a range of inertia values. As it will be explained further, that is a consequence of the regular topology. Also, the number of plateaus increase by raising k.

Figure 1

From top to bottom, mean-field results for regular networks for k = 4, k = 8 and k = 12, respectively. The left panels show the behavior of m versus f for distinct θ’s, whereas the right ones show the respective phase diagrams. ORD and DIS correspond to the ordered and disordered phases, respectively. Location of forward and backward transitions are exemplified by arrows in panel (c).

Figure 6

The phase diagram θ versus f for the MV with k = 20 in a bidimensional lattice. Continuous and dashed lines correspond to critical and discontinuous phase transitions, respectively. Inset: The same, but for the RR topology. Circles (×) correspond to the increase (decrease) of f starting from an ordered (disordered) phase.

Numerical results

Numerical simulations furnish more realistic outcomes than the MFT ones, since the dynamic fluctuations are taken into account. The actual simulational protocol is described in [Methods]. Starting with the random topology, Fig. 2 shows the phase diagrams for k = 4, 8, and k = 12, respectively.

Figure 2

RR Networks: From the top to bottom, numerical results for k = 4, k = 8 and k = 12, respectively. The left panels exemplify the behavior of 〈m〉 versus f for θ = 0.33 (k = 4) and 0.35 (k = 8 and 12), whereas right ones show the phase diagrams. Inset: Reduced cumulant U4 vs. f for θ = 0.05 and k = 8. Circles (times) correspond to the increase (decrease) of f starting from an ordered (disordered) phase.

First of all, we observe that the positions of plateaus are identical than those predicted from the MFT. Also, the phase transition is continuous for k = 4, irrespective the inertia value. For θ > 1/3, no phase transition is displayed and the system is constrained into the disordered phase. In all cases (see, e.g., Fig. 2(a) for k = 8 and θ = 0.05), the critical transitions are absent of hysteresis and U4 curves cross at with .

Opposite to the low k, discontinuous transitions are manifested for k = 8 and 12 in the regime of pronounced θ. More specifically, the crossovers take place at θ = 1/3 and θ = 1/4 for the former and latter k, respectively. Notwithstanding, there are some differences between approaches. As expected, MFT predicts overestimated transition points than numerical simulations. Although MFT predicts a continuous phase transition in the interval (k = 12), numerical simulations suggest that it is actually discontinuous.

In Fig. 3, the previous analysis is extended for regular (bidimensional) versions. In order to mimic the increase of connectivity, the cases k = 4, 8 and 12 cases are undertaken by restricting the interaction between the first, first and second, first to third next neighbors, as exemplified in panels (a–c) in Fig. 4, respectively.

Figure 3

Bidimensional regular lattices for distinct system sizes N = L × L: Left panels show the reduced cumulant U4 vs f for the nearest neighbor (a) second-neighbor (c) and third-neighbor (e) versions, respectively. Inset: the same but for the variance χ. Right panels show their correspondent phase diagrams. In all cases, continuous lines correspond to critical phase transitions.

Figure 4

Local configuration for the (bidimensional) versions with interactions between the first (a) first and second (b) first to third (c) and first to fourth (d) next neighbors.

The positions of plateaus are identical than both previous cases, but with lower f’s. This is roughly understood by recalling that homogeneous complex networks exhibit a mean-field structure, whose correspondent transition points are thus larger than those from regular lattices. Similarly, all critical points are obtained from the crossing among U4 curves, but the value is different from the RR case, following to the Ising universality class value [1,6,20]. Thereby, there is an important difference between random and regular structures: The phase transitions are continuous irrespective the inertia value for k from k = 4 to k = 12.

An absolutely different scenario is unveiled by extending interactions range up to the fourth next neighbors spins (mimicking the case k = 20) [see, e.g., Fig. 4(d)] and large inertia values, in which the phase transition becomes discontinuous (see, e.g., Fig. 5 for θ = 0.35). Contrary to the random complex case, hysteresis is absent [panel (a)], on the other hand the order-parameter distribution exhibits a bimodal shape [panel (b) shows the value of f in which the peaks describing the ordered and disordered phases present equal areas]. Complementary, U4 behaves quite differently from continuous transitions and presents a minimum whose value decreases with N [panel (c)]. Also, the maximum of χ increases with N (inset). These features, analogous than discontinuous absorbing phase transitions[13], also discloses the emergence of discontinuous transitions in regular lattices as k is increased. The f’s (estimated from the equal area position, maximum of χ and minimum of U4) scales with N−1 [panel (d)], from which one obtains the estimates f0 = 0.0687(1) (equal area and maximum of χ) and f0 = 0.0689(1) (minimum of U4) [see Methods for obtaining the finite-size scaling relation].

Figure 5

Results for k = 20 and θ = 0.35: Panel (a) compares the order parameter 〈m〉 versus f for the RR network (circles and stars) and regular lattice (symbol×) for a system with N = 104 sites. Regular lattice case: Panels (b) and (c) show the non normalized equal area probability distribution and the U4 vs. f for distinct L’s (N = L × L), respectively. Inset: The variance χ versus f. In (d), the positions of maxima of χ, minima of U4 and equal area versus 1/N.

In Fig. 6, the phase diagram is presented. As in previous cases, the positions of the plateaus are identical to the RR for k = 20 (see inset and ref.[19]). The phase coexistence occurs for θ > 1/3, larger than θ > 3/13 (RR structure). For θ < 1/3, the phase transition is continuous, although U4 presents a value different from in the interval 2/7 < θ < 1/3.

Origin of Plateaus

Since the transition rate depends only on the signal of resulting argument in Eq. (1), the phase diagrams will present plateaus provided the number of neighbors is held fixed. Generically, let us take a lattice of degree k with the central site σ0 with and nearest neighbors with spins +1 and −1, respectively (obviously ). Taking for instance σ0 = −1 (similar conclusions are earned for σ0 = 1). In such case, the argument of sign(X) reads , implying that for all the transition rate −1 → +1 will be performed with the same rate 1 − f and thus the transition points are equal. Only for θ > θ the transition −1 → +1 is performed with probability f. Table 1 lists the plateau points θ for k = 8 and distinct ′s. For example, for and , all transition rates are equal, implying the same f for such above set of inertia. For the second local configuration becomes different and thereby f is different from the value for θ < θ. Keeping so on with other values of , the next plateau positions are located. It is worth mentioning that leads to negative θ’s, that not have been examined here.

Table 1

For the central site σ0 = −1 and connectivity k = 8, the signal function for distinct local configurations.

σ 0	nk+	nk−	X > 0	θp
−1	4+	4−	−θ > 0	0
−1	5+	3−	1−5θ4>0	15
−1	6+	2−	1−3θ2>0	13
−1	7+	1−	3−7θ4>0	37
−1	8+	0−	1 − 2θ > 0	12

X is the value of resulting expression (see main Text) and θ denotes the plateaus positions.

Methods

We consider a class of systems in which each site i can take only two values ±1, according to its “local spin” (opinion) σ, is “up” or “down”, respectively. The time evolution of the probability P(σ) associated to a local configuration σ ≡ (σ1, .., σ, σ) is ruled by the master equationwhere the sum runs over the N sites of the system and σ ≡ (σ1, .., −σ, σ) differs from σ by the local spin of the i–th site. From the above, the time evolution of the magnetization of a local site, defined by m = 〈σ〉, is given bywhere w−1→1 and w1→−1 denote the transition rates to states with opposite spin. In the steady state, one has that

By following the formalism from refs[19,24,25], the transition rates w−1→1 and w1→−1 in Eq. (4) are decomposed asandwhere denote the probabilities that the node i of degree k, with spin σ = −1 (σ = 1) changes its state according to the majority (minority) rules, respectively. Such probabilities can be written according towith p±1 being the probability that a nearest neighbor is ±1 and and corresponding to the lower limit of the ceiling function, reading and .

Since we are dealing with uncorrelated structures with the same degree k, p± is simply (1 ± m)/2, from which Eq. (4) readswith being evaluated from Eq. (7). Thus, the solution(s) of Eq. (8) grant the steady values of m.

An alternative way of deriving the MFT expressions consists in writing down the transition rates as the sum of products of the local spins ,where σ is the product of spins belonging to the cluster of k sites, and c is a real coefficient. For example, for k = 8 and θ = 0, we have that , yielding the critical point , in full equivalency with f obtained from Eq. (8).

The numerical simulations will be grouped into two parts: Random regular (RR) network and bidimensional lattices.

The former is generated through a configuration model scheme[26] described as follows: For a system with N nodes and connectivity k, we first start with a set of Nk points, distributed in N groups, in which each one contains exactly k points. Next, one chooses a random pairing of the points between groups and then creates a network linking the nodes i and j if there is a pair containing points in the i-th and j-th sets until Nk/2 pairs (links) are obtained. In the case the resulting network configuration presents a loop or duplicate links, the above process is restarted. The bidimensional topologies also present connectivity k, but it forms a regular arrangement. Note that both structures are quenched, i.e., they do not change during the simulation of the model.

For a given network topology, and with N, f, and θ held fixed, a site i is randomly chosen, and its spin value σ is updated (σ → −σ) according to Eq. (1). With complementary probability, the local spin remains unchanged. A Monte Carlo (MC) step corresponds to N updating spin trials. After repeating the above dynamics a sufficient number of MC steps (in order of 106 MC steps), the system attains a nonequilibrium steady state. Then, appropriate quantities, including the mean magnetization , its variance χ = N[〈m2〉 − 〈m〉2] and the fourth-order reduced cumulant are evaluated, in order to locate the transition point and to classify the phase transition. We have evaluated above averages for a total of 5.106 MC steps.

For both topologies, continuous phase transitions are trademarked by the algebraic behaviors of and , where β/ν and γ/ν are their associated critical exponents. Another feature of continuous transitions is that U4, evaluated for distinct N's, intersect at . Notwithstanding, the set of critical exponents as well as the crossing value depend on the lattice topology[20,22]. Off the critical point, U4 reads U4 → 2/3 and 0 for the ordered and disordered phases, respectively, when N → ∞.

In similarity with the MFT, numerical analysis of discontinuous transitions in complex networks is commonly identified through the presence of order parameter hysteresis. For the “forward transition”, simulations are started at f = 0 from a full ordered phase (|m| = 1), in which the tuning parameter f is increased by an amount Δf and the system (final) configuration at f used as the initial condition at f + Δf. Then, the system will jump to a disordered phase (|m| = 0) at a threshold value f. Conversely, the system evolution now starts in the full disordered phase and one decreases f (also by the increment Δf), so that the ordered phase (|m| ≠ 0) will be reached at f. Both “forward” and “backward” curves are not expected to coincide themselves at the phase coexistence. The value of Δf to be considered will depend on θ, but in all cases it is constrained between 0.005 and 0.02.

In contrast to complex structures, the behavior of discontinuous transitions is less understood in regular lattices. Recently, a phenomenological finite-size theory for discontinuous absorbing phase transitions was proposed[13], in which no hysteretic nature is conferred, but instead one observes a scaling with the inverse of the system size N−1. Here, we extend it for Z2 up-down phase transitions. Such relation can be understood by assuming that close to the coexistence point, the order-parameter distribution is (nearly) composed of a sum of two independent Gaussians, with each phase σ [σ = o (ordered) and d (disordered)] described by its order parameter value m in such a way thatwhere each term readswhere χ is the variance of the σ–gaussian distribution, Δf = f − f0 denotes the “distance” to the coexistence point f0 and each normalization factor F reads

Note that (10) leads to the probability distribution being a sum of two Dirac delta functions centered at m = m and m = m at f = f0 for N → ∞. For f − f0 → 0+(−), one has a single Dirac delta peak at m = m (m ≠ 0). The pseudo-transition points can be estimated under different ways, such as the value of f in which both phases present equal weight (areas). In such case, from Eq. (10) it follows that for

Since N is supposed to be large, the right side of Eq. (12) becomes small and thus (f − f0) is also small. By neglecting terms of superior order (f − f0)2, we have thatimplying that the difference f − f0 scales with the inverse of the system size N. Evaluation of the position of peak of variance χ provides the same dependence on N−1, whose slope is the same that Eq. (13) [see, e.g., panel (d) in Fig. 5].

In this work, we have classified the (discontinuous) phase transition in regular lattices based on the behavior of the three aforementioned quantities.

Conclusions

A discontinuous phase transition in the standard majority vote model has been recently discovered in the presence of an extra ingredient: the inertia. Results for distinct network topologies revealed the robustness of such phase coexistence trademarked by hysteresis, bimodal probability distribution and others features[19]. Here, we advanced by tackling the essential ingredients for its occurrence. A central conclusion has been ascertained: discontinuous transitions in the MV also manifest in low dimensional regular topologies. Also, its finite size behavior (entirely different from the network cases), is identical to that exhibited by discontinuous phase transitions into absorbing states[13]. This suggests the existence of a common and general behavior for first-order transitions in regular structures. In addition, low connectivity leads to the suppression of the phase coexistence, insighting us that not only the inertia is a central ingredient, but also the connectivity. Numerical simulations reveal that for random regular networks, the minimum neighborhood is k = 7, whereas about k = 20 is required for changing the order of the transition in bidimensional lattices. Summing up, the present contribution aimed not only stemming the key ingredients for the emergence of discontinuous transitions in an arbitrary structure, but also put on firmer basis their scaling behavior in regular topologies.