Three-State Majority-vote Model on Scale-Free Networks and the Unitary Relation for Critical Exponents

We investigate the three-state majority-vote model for opinion dynamics on scale-free and regular networks. In this model, an individual selects an opinion equal to the opinion of the majority of its neighbors with probability 1 − q, and different to it with probability q. The parameter q is called the noise parameter of the model. We build a network of interactions where z neighbors are selected by each added site in the system, a preferential attachment network with degree distribution k−λ, where λ = 3 for a large number of nodes N. In this work, z is called the growth parameter. Using finite-size scaling analysis, we obtain that the critical exponents \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta /\bar{\nu }$$\end{document}β/ν¯ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma /\bar{\nu }$$\end{document}γ/ν¯ associated with the magnetization and the susceptibility, respectively. Using Monte Carlo simulations, we calculate the critical noise parameter qc as a function of z for the scale-free networks and obtain the phase diagram of the model. We find that the critical exponents add up to unity when using a special volumetric scaling, regardless of the dimension of the network of interactions. We verify this result by obtaining the critical noise and the critical exponents for the two and three-state majority-vote model on cubic lattice networks.

temperature T c , and the average number of interacting spins 〈k〉 for the scale-free networks investigated 17 . They also find that the critical temperature of the Ising model is proportional to the logarithm of the finite system size N 18 . This result was also verified for the Potts model with two states on semi-directed Barabási-Albert networks 20 .
The complex network framework allowed physicists to propose models aiming to investigate the critical phenomena in social dynamics [21][22][23][24][25][26][27][28][29][30][31][32][33][34][35] . Although human individuals may be difficult to simulate, the dynamics of people in groups seem to be conceivable to model due to herd effects and other group behaviors. By employing statistical physics techniques, several works confirm that these models present order-disorder transitions and universality among other standard characteristics of condensed matter physical systems. From microscopic interactions among individuals in a social network, several sociophysics models exhibit a robust collective behavior. Concerning opinion dynamics in complex networks, we highlight the majority-vote model with noise 24 . In this model, an individual assumes the same opinion of the majority of its neighbors with probability p and the opposite opinion with probability q = (1 − p). In its two-state version, the opinion of an individual is represented by the spin variable σ, which assumes the value +1 or −1 in a given time. In a regular square lattice of interactions, this model undergoes a nonequilibrium phase transition at a critical probability .  q 0 075 c , and the critical exponents are the same as those of the equilibrium two-dimensional magnetic Ising model.
The three-state majority-vote model with noise is a system of spins, where each one is allowed to be in one of three states, that is σ = 1, 2, or 3 [31][32][33] . As in the two-state model, each spin assumes the state of the majority of its neighboring spins with probability p = (1 − q) and the opposite state with probability q, which is known as the noise parameter of the model. The increase of the parameter q promotes the formation of different opinion configurations in the model, acting as an increase of the social temperature of the system. For a regular square lattice network, the three-state majority-vote model with noise present an order-disorder phase transition at the critical noise value .  q 0 118 c , where the consensus in the system vanishes. In this work, we investigate the influence of a network with preferential attachment, built by using the classical Barabási-Albert model, on the three-state majority-vote model with noise. We perform Monte Carlo simulations to estimate the critical noise parameter as a function of the growth parameter q c (z). Moreover, we use the standard finite-size scaling techniques to obtain the critical exponents for several values of the parameter z for the networks investigated. Based on the results, we propose a unitary relation to examine the criticality of the system. We conjecture that this relation is universal, regardless of the network of interactions. We also perform simulations for the three-state majority-vote model on cubic networks that confirm our results, and we obtain the critical noise for this system and its critical exponents.
The remainder of this paper is divided into four sections. In Section II, we describe the three-state majority-vote model with noise for opinion dynamics, the network construction process, and introduce the relevant quantities used in our simulations. Section III contains our results for complex and regular networks, along with a discussion. In Section IV, we present our conclusions and final remarks.

the Model
The Barabási-Albert network. The three-state majority-vote model with noise consists in a set of spin variables {σ i } with i = 1, 2, …, N, where each variable can be assigned to one of the values σ = 1, 2, or 3, representing the opinion for an individual in the community at a given time t. The individuals are distributed in the nodes of a scale-free network of social interactions with N sites. That is, before adding individuals and their opinions in our system, we build a network of interactions. We start with a core of z +1 fully connected nodes, and then we add new nodes -one at a time -with z free links, which will be connected by preferential attachment to the existing nodes of the network. The probability of connecting a new node j to node i, Π k ( ) i , depends on the degree k i of the node i. Thus, for Barabási-Albert networks with linear preferential attachment we write where the summation is equal to the total number of existing links in the scale-free network, and more nodes are added to the network until it reaches a total of N sites. A double connection to the same site is forbidden.
In Fig. 2 we show an illustration of the Barabási-Albert network with N = 100 sites for z = 5 (left). Observe that some nodes have a high number of connections, despite the small average value of links per site in the network. We also present the histogram of the degree of the nodes for networks with N = 20000 nodes and different www.nature.com/scientificreports www.nature.com/scientificreports/ values of the growth parameter z (right), where we obtain the characteristic scale-free degree distribution plot with exponent λ = 3. Note that the decay exponent λ is the same, even for different values of the average degree per node z as expected 2 .
Dynamics and numerical quantities. The dynamics of the system consists of a generalization for three opinions of the two-state majority-vote model [24][25][26][28][29][30][31][32][33][36][37][38][39][40][41][42][43] . For a randomly selected individual σ i we determine the majority opinion of the individuals that are linked to it. With probability 1 − q, the selected agent σ i adopts the dominant opinion of its neighbors (follows the majority), and with probability q, the individual adopts a different opinion (follow a minority). For a tie between the three states, the selected agent σ i changes to any opinion with the same probability equal to 1/3. For the case of a tie between two majorities, σ i assumes one of these tied states with probability (1 − q)/2, and the minority with probability q. Finally, for the case of a single majority, σ i follows any of the two minorities with probability q/2, and the majority with probability 1 − q. That is, if n α is the number of neighbors of the individual σ i in a given state α = 1, 2, 3, then, the probabilities for σ i to assume the opinion α = 1 is: P n n n q P n n n q P n n n q P n n n q n n n These transition rules present the C 3v symmetry for the simultaneous change of all opinions. The probabilities for the other states σ = 2 and 3 are obtained by symmetry operations. The total number of individuals connected to σ i is n = n 1 + n 2 + n 3 . We remark that all probabilities satisfy the relation In our simulations, we choose the consequential update for the opinion dynamics of the agents. To investigate the critical behavior of the three-state majority-vote model, we first calculate the average opinion, defined by where the sum is over all sites in the scale-free network of social interactions and δ(α, σ i ) is the Kronecker delta function. In this way, to investigate the critical behavior of the model, we consider the following numerical quantities: the magnetization M, the magnetic susceptibility χ, and the Binder's fourth-order cumulant U defined by t c t c 2 2 where q is the noise parameter, z is the growth parameter of the network, N is the total number of individuals or agents, 〈…〉 t represents the time average taken in the stationary regime and 〈…〉 c denotes the configurational average. The critical behavior of the model is investigated by performing computer simulations and by using finite-size scaling analysis. The three-state majority-vote dynamics evolves over time according to the probability rules given by Eq. (2). After a transient, the model reaches a steady state that presents complete order, partial order(disorder) or complete disorder. That is, for q = 0 the system exhibits an ordered steady macrostate, characterized by the predominance of individuals with one of the three possible opinions. Assuming for this case that σ = 1 ∀ N, we obtain and m = 1, yielding M = 1 for q = 0. The upper limit for q, that is, the infinite social temperature, is obtained when the probability of agreeing with the majority is equal to the probability of agreeing with any minority, thus − = ⇒ = q q q 1 /2 2/3. In this case, any opinion or individual state σ = 1, 2 or 3 can be found with equal probability in the network of social interactions. Thus, leading to  m m m , , 0 and M(q = 2/3) = 0 in the thermodynamic limit N → ∞.

Discussion and Results
Monte carlo simulations. We perform numerical simulations on scale-free networks with sizes ranging from N = 1000 to 20000 using the Monte Carlo method. For each pair of values for the parameters q and z, we set a fraction of sites f 0 = 0.8 to point to one opinion, i.e., σ i = 1, and the remaining fraction 1 − f 0 to point to the other two options, equally distributed. We next select a randomly chosen individual and update its opinion with the probabilities given by Eq. (2). This process is repeated N times to allow that the opinion of each agent is updated once (on average) for each Monte Carlo Step (MCS). We skip 10 5 MCS in the simulation to overcome the initial transients and allow the system to reach a steady-state, characterized by the set of parameters q and z. Next, we perform the time averages in the following 2 × 10 5 MCS. For each realization, we generate at least 100 independent random samples to obtain the configurational averages. Different values for f 0 yield the same final steady-state for the system, which becomes ordered or disordered, depending on the values for q, z, and N. In the ordered phase, the dominant opinion is found in one of the possible states 1, 2, or 3. In the disordered phase, the three opinions are equally distributed in the network of social interactions.
In Fig. 3, we illustrate the effect of the network of interactions with preferential attachment in the consensus (order) of the system. We show the plot for the magnetization M(q, z, N), for the susceptibility χ(q, z, N) and for the Binder's fourth-order cumulant U(q, z, N) as a function of the noise q for N = 20000 and z = 2, 3, 4, 5, …, 10. We find that for small values of the noise parameter q, the system presents an ordered macrostate, or phase, where M(q, z, N) = O (1). In this phase, the society becomes ordered, with a preference for a dominant opinion. By increasing the social temperature q, the magnetization will continuously decrease to zero near a critical value q c , denoting the second-order phase transition of the system. In this region (q > q c ), there is no prevailing opinion in the society, and every state can be found with the same probability. In Fig. 3(b) the magnetic susceptibility χ(q, z, N) exhibits a maximum near some critical value q c , where the transition order-disorder occurs. This behavior is also denoted by the rapid decrease of the Binder's fourth-order cumulant U(q, z, N), showed in Fig. 3(c). From www.nature.com/scientificreports www.nature.com/scientificreports/ the results, we obtain that the critical noise is an increasing function of the growth parameter of the network z, indicating that the consensus is stronger if there are more connections to an individual in the society.
Next, we consider the finite-size effects on our measured quantities. In Fig. 4 we show the (a) magnetization  M(q, z, N), (b) the susceptibility χ(q, z, N) and (c) the Binder's fourth-order cumulant U(q, z, N) versus the noise parameter q for z = 10, and several system sizes N. We remark that M(q, z, N) ≠ 0 for high values of the parameter q due to finite-size effects. The susceptibility χ(q, z, N) exhibits a sharper peak as we increase the system size, and the position of its maximum in the horizontal axis depends on N. Thus, we write the pseudocritical noise as q c (z, N), and the values for q c (z, N) are our first estimative for the critical noises of the system. To obtain the critical noise for each z in the thermodynamic limit, q c (z), we calculate the Binder cumulant of the system U(q, z, N). In Fig. 4(c) we show the Binder's fourth-order cumulant as a function of the social temperature q for different system sizes. The critical noise parameter q c (z) can be estimated for the value of q where the curves of U(q, z, N) for different system sizes intercept each other. In this figure, we estimate q c = 0.513 (1). Figure 5(a) illustrates the dependence of the Binder's fourth-order cumulant on the noise q for z = 5 and different system sizes. The curves for different values of N intercept in the region 0.43 < q < 0.44. Figure 5(b) shows the details of the Binder cumulant data, along with a polynomial fit for the lines near the interception. At this point, the critical noise does not depend on the system size, and for this case, shown in figure q c = 0.4326 (4). By calculating the cumulant U(q, z, N) for other values of z, we obtain the phase diagram for the three-state majority-vote  www.nature.com/scientificreports www.nature.com/scientificreports/ model on Barabási-Albert scale-free networks shown in Fig. 5(c). The orange region denotes the ordered phase of the system, where one of the three opinions is the majority state of the system. In this result, the error bars are smaller than the thickness of the line. Note that the critical noise increases with the growth parameter z, which controls the average number of connections 〈k〉 of a given individual in the social system. Here, we obtain q c = 2/3 as z → ∞, where q c = 2/3 is the infinite temperature equivalent for the social system of the three-state majority-vote model. This result agrees with the critical temperature obtained in other studies for ferromagnetic spin systems in scale-free networks 17,18 . the unitary relation and scaling results. To obtain the critical exponents in complex networks, we propose that near the critical noise q c the correlation length ξ scales with the actual volume of the system 44,45 as ξ ∼ .
N (9) Thus, the pseudocritical noise q c (N), the magnetization M(q, z, N), the susceptibility χ (q, z, N), and the Binder cumulant U(q, z, N) satisfy the finite-size scaling relations where ε = q − q c is the distance to the critical noise, b is a constant, and M , χ , and ∼ U are scaling functions that only depend on the scaled variable ε = ν x N 1/ . For regular networks, we recall that N = L d , where d is the effective dimension of the network, and L is an effective linear size of the system. In this case, we obtain for the magnetization and the magnetic susceptibility˜ε d / 1 /  We use the notation ν instead of v since we changed the correlation length scaling relation from the usual linear scaling ξ ~ L to a "volumetric scaling" ξ ~ L d . In this case, the hyperscaling relation now reads β ν γ ν regardless of the effective dimension d of the network. This result allows us to remark that the hyperscaling relation cannot be used to estimate the dimension of these networks when using the volumetric scaling ξ ~ L d , in contrast to the results of previous studies 25,28,33,38,39 . Nevertheless, the unitary relation (Eq. (16)) was verified in these works for random graphs and scale-free networks. In this context, we rewrite the unitary relation by denoting a new exponent upsilon υ, defined as where we conjecture that the exponent υ = 1 for any network under the condition of the volumetric scaling of Eq. (9). In this work, we denote the equation υ = 1 as the unitary relation for critical exponents. We validate the consistency of this result according to the comparison with the numerical findings for the critical exponents β ν / and γ ν / for regular and complex networks. By calculating the logarithm of Eqs. (10), (11) and (12) at the critical point q c , we obtain an explicit relation involving the critical exponents, the measured quantities and the system volume N ln , (20) and we use the Eqs. (18), (19) and (20) to obtain the critical exponents of the system. (2020) 10:8255 | https://doi.org/10.1038/s41598-020-63929-1 www.nature.com/scientificreports www.nature.com/scientificreports/  Table 1 provides critical noise, critical exponents, and the unitary relation values for each growth parameter investigated in the model. Note that the critical exponent for the magnetization (susceptibility) is a decreasing (an increasing) function of the growth parameter z. The critical noise q c increases with z, while the critical exponent ν 1/ decrease with z. For all values of the critical exponents, we obtain υ ∼ 1 as expected.
By relating the critical exponents, we obtain the characteristic unitary line of the model showed in Fig. 7. Here, we plot the values of the critical exponents γ ν / versus β ν 2 / . The linear fit of the data yield = − . + . y x 0 96(1) 1 02(1), and an averaged unitary exponent υ = .
1 02(1) for the three-state majority-vote model on Barabási-Albert networks. We conjecture that this line is universal, regardless of the geometric structure of the network used in the model. Thus, one can use the volumetric scaling ξ ~ N with the unitary relation Eq. (17), and the Eqs. (18), (19) and (20) to obtain the critical exponents and the unitary line for any spin model under consideration, with or without a system size length clearly defined. Figure 8 shows the plot of the rescaled (a) magnetization / 0300, γ ν = . / 044 and ν = .
1/ 045 with q c = 0.5282. Other values for the growth parameter z exhibit the same features and the same qualitative results for the data collapse of the magnetization, susceptibility, and Binder cumulant. From our simulation results and analysis, we conclude that the three-state majority-vote  susceptibility ln χ(q c , z, N), and (c) ln [q c (N) − q c ] versus ln N using z = 2, 5, 14, 20 and 50, with g = 1, 2, 3, 4 and 5. Here, g is an integer used to better display the data. In all plots the error bars are smaller than the symbol size. www.nature.com/scientificreports www.nature.com/scientificreports/ model defined on Barabási-Albert networks and on the Erdös-Rényi random graphs belong to different universality classes when the volumetric scaling ξ ~ N is used 33 .

Unitary relation on regular networks. To confirm the validity of our statements, we performed Monte
Carlo simulations for the majority-vote model with two and three states on regular square lattices and cubic networks. From our simulations of the three-state majority vote model on cubic networks, we obtain the Fig. 9 that shows the (a) magnetization M(q, L), the (b) susceptibility χ(q, L) and the (c) Binder cumulant U(q, L) versus the noise parameter q, where N = L 3 . We observe some familiar results such as M(q, L) → 0 for q > q c , with L → ∞, and χ(q, L) that also exhibits a sharper peak as we increase L. From the Binder parameter, we find the critical noise of the model q c = 0.25230 (2).
By performing numerical simulations for the majority-vote model with two and three states on regular square lattices and on cubic networks, we build the Table 2 with the critical exponents for the magnetization and susceptibility. We also calculate the unitary relation υ and effective dimension d obtained for each model with the volumetric (ξ ~ N) and linear (ξ ~ L) scalings, respectively. We conclude that the critical exponents with the linear and volumetric scalings relate by β ν β ν = d / / and γ ν γ ν = d / / , as expected by recalling that ξ ~ L d . In Fig. 10, we plot the logarithm of the magnetization, and magnetic susceptibility used to obtain the critical exponents for the majority-vote model with two and three states on a square lattice and cubic networks, where we used the volumetric scaling. Our results confirm that the unitary relation holds for this model on these networks, and it points that the effective dimension obtained by previous works with the majority-vote model on random graphs and on Barabási-Albert networks might be not equal to unity 25   / 03, γ ν = . / 044 and ν = .
1/ 045. The error bars are smaller than the symbol size.

conclusion and final Remarks
We have investigated the dynamics of the three-state majority-vote model for opinion dynamics on Barabási-Albert networks. We obtained the phase diagram and the critical exponents of the model, and we verified that the second-order phase transition occurs for networks with growth parameter z > 1. The critical exponents γ ν / and ν 1/ decrease with z, while β ν / increases. We also find that the critical noise q c is an increasing function of the growth parameter z, which converges to the infinite social temperature value 2/3 as z → ∞. In    for Q = 2. For a cubic network (b), we obtain β ν = .
Scientific RepoRtS | (2020) 10:8255 | https://doi.org/10.1038/s41598-020-63929-1 www.nature.com/scientificreports www.nature.com/scientificreports/ other words, if z is big enough, the three-state opinion society remains ordered, with a consensus, even for a high value of the social disorder (or preference for dissensus q). From Fig. 5 we estimate that this occurs for z ≥ 20, where q c ~ 2/3, indicating that it is harder to destroy consensus in a more connected society.
By assuming that near criticality, the correlation length ξ scales with the actual volume of the system ξ ~ N, we found that the new hyperscaling relation is equal to 1, regardless of the effective dimension of the network of interactions. Nevertheless, by performing Monte Carlo simulations we verified that the unitary relation υ β ν γ ν ≡ + = 2 / / 1 for all values of the growth parameter z investigated. The unitary relation was also verified for the majority-vote model with two and three states on regular square lattices and cubic networks, with well defined effective dimensions. Our analysis sheds light on a curious result of several works, where authors find that the effective dimension of different complex networks is equal to one 25,28,33,38,39 . We show that this result is a consequence of the scaling used and that this is not the real dimension of those networks, excluding eventual coincidences. Furthermore, the unitary relation here defined for a sociophysics model holds for other spin-like models in condensed matter systems, regarding the geometric dimensions in which these systems are embedded.
We remark that obtaining the effective dimension for some complex networks -such as the random graphs and the Barabási-Albert networks -using the finite-size scaling analysis near the critical points of such spin systems, and the relations between its critical exponents, remains a task to perform. Although other methods and scalings have been successfully developed 46,47 . Our results may suggest future research finding and using other finite-size scaling relations, such as some power-law with a characteristic length instead. We also recommend the use of the unitary relation (Eq. (17)) and the unitary line (Fig. 7), defined and discussed in this work, as a form investigation of the criticality of systems with complex interactions of unknown effective dimension, or with a linear size not precisely defined.