Kinetic Continuous Opinion Dynamics Model on Two Types of Archimedean Lattices

1. Introduction

In 1986 Galam introduced the use of local majority rules to study voting systems as bottom-up democratic voting in hierarchical structures [1], see also the references [2–4]. Although sociophysics has been rejected by some physicists in the eighties [5], it is has become today an active and promising area of research for interdisciplinary physicists [6, 7].

In this same context, in 1982, de Oliveira [8] proposed a non-equilibrium model called Majority Vote Model (MVM). On two-dimensional lattices it shows critical phenomena with critical exponents ν, β, and γ, as for [8–10] the equilibrium Ising model [11, 12], in agree with a hypothesis of Grinstein et al. [13].

In 2012 Biswas et al. [14] proposed a kinetic model of opinion formation. This model kinetic continuous opinion dynamics (KCOD) presents mutual interactions between the individual i, j that can be both positive and negative. In this model the fraction of negative interactions is represented by a parameter p in order to characterize the different types of distributions for the mutual interactions. The results of the continuous version the KCOD model, obtained through numerical simulations indicate the existence of a universal continuous phase transition at p = p_c with exponents of mean field (νd = 2.00(1), β = 0.50(1), and γ = 1.00(1)).

Similar to this KCOD model is the one of Deffuant et al. [15], where each person i selects another person j to talk to (no lattice) and both move their opinion toward that of the other person. For the Krause-Hegselmann model [16], each person averages over the opinions of the others in the population, again no lattice. While these two models use continuous opinions, those in the Sznajd model [17] usually are Ising-like (+1 or −1) and restricted to a square lattice with four nearest neighbors. Two Sznajd people happening to agree in their opinion convince all their six neighbors of this opinion. In all these models one may start from a random distribution of opinions and then check if the opinions all converge to one consensus, two opposing camps of opinions, or a fragmentation into many opinion clusters [18–20].

Recently, Mukherjee and Chatterjee [21] studied the KCOD model on square and cubic lattices (2D and 3D). Their numerical results indicate that the same critical behavior of the KCOD model and the Ising model in the corresponding dimensions.

In this work, we studied the KCOD on two Archimedean lattices—namely, (4, 6, 12), and (4, 8²)—through extensive Monte Carlo simulations. Pictures of the (4, 6, 12) and (4, 8²) AL are shown in Figure 1. The AL are vertex transitive graphs that can be embedded in a plane such that every face is a regular polygon. Kepler showed that there are exactly 11 such graphs. The AL are labeled according to the sizes of faces incident to a given vertex. The face sizes are sorted, starting from the face for which the list is the smallest in lexicographical order. In this way, the square lattice gets the name (4, 4, 4, 4) (abbreviated to (4⁴)), honeycomb is called (6³), and Kagome is (3, 6, 3, 6). Here, we also compared our results with those of the MVM made on (3, 4, 6, 4) and (3⁴, 6) AL.

FIGURE 1

Figure 1. Pictures of the (4, 6, 12) and (4, 8²) AL.

One of our goals, besides finding the critical exponents of the KCOD model on (4, 6, 12) and (4, 8²) AL, is to verify the Grinstein et al. criterion [13], for non-equilibrium stochastic spin systems with up-down symmetry on (4, 6, 12) and (4, 8²) AL belong to the same universality class as the Ising model on (4⁴) as suggested by Grinstein et al. [13]. Here, we also compared our results with those of the MVM on (4, 6, 12) and (4, 8²) [22].

2. Definition and Simulation

The KCOD model [14] is defined by a set of individuals with continuous opinion variables o_i(t), where the opinion of a person i at time t, takes the values in the interval [−1, +1], is situated on every node of the (4, 6, 12) and (4, 8²) AL with N = 12L² sites for (4, 6, 12) and N = 4L² sites for (4, 8²). In a population of N individuals, opinions change out of pair-wise interactions via mutual influences/couplings μ_ij as:

$\begin{array}{ll}{o}_i(t+1)=o_i(t)+{\mu }_{ij}{o}_j(t)& \left(1\right)\end{array}$

$\begin{array}{ll}{o}_j(t+1)=o_j(t)+{\mu }_{ij}{o}_i(t).& \left(2\right)\end{array}$

The pairs i, j are unrestricted, meaning the original model is defined on a fully-connected graph, giving a mean-field-like limit (infinite range interactions). The pairwise interaction implies no sum over the index j. with real random couplings μ_ij. Agent i updates his/her opinion via Equation (1) by interacting with agent j and is influenced by the coupling μ_ij. The opinions are limited to −1 ≤ o_i(t) ≤ 1. If the opinion value of an individual become higher (lower) than +1 (−1), then it is set equal to +1 (−1). This bound, along with Equation (1), defines the dynamics of the model. Here, μ_ij is a continuous random variable defined in the range [−1, +1], i.e., it takes a random real value in the range [−1, 0] or ([0, 1]) with probability p or (1−p). In other words, the disorder parameter p denotes the fraction of negative pairwise interactions. The average opinion $O=\left|\sum _i{o}_i\right|/N$ measures the ordering in the system. As a function of the fraction p of negative interactions a symmetry breaking transition occurs between an ordered and a disordered phase: below a certain value p_c of the parameter p, the system is ordered (giving a non zero value of the opinion O, defined in the following), while it is disordered above p_c (O = 0).

The critical properties of model we are interested in are the order parameter O, the order parameter fluctuations (susceptibility) OF and the reduced fourth-order cumulant of the O, namely here by O₄, defined as

$\begin{array}{ll}O(p)\equiv \langle O\rangle ,& \left(3\text{a}\right)\end{array}$

$\begin{array}{ll}OF(p)\equiv N(\langle O^2\rangle -{\langle O\rangle }^2),& \left(3\text{b}\right)\end{array}$

$\begin{array}{ll}{O}_4(p)\equiv 1-\frac{\langle O^4\rangle }{3{\langle O^2\rangle }^2},& \left(3\text{c}\right)\end{array}$

where 〈⋯ 〉 stands for time averages, computed at the steady states. N_run independent simulations are averaged over.

The quantities O, OF, and O₄ depend on the disorder parameter p and obey

$\begin{array}{ll}O=L^{-\beta /\nu }{f}_o(x),& \left(4\text{a}\right)\end{array}$

$\begin{array}{ll}OF=L^{\gamma /\nu }{f}_{of}(x),& \left(4\text{b}\right)\end{array}$

$\begin{array}{ll}\frac{d{O}_4}{dp}=L^{1/\nu }{f}_{o4}(x),& \left(4\text{c}\right)\end{array}$

(finite-size scaling) with β, γ, and ν as the usual critical exponents, f_o(x), f_of(x), f_o4(x) as the finite-size scaling functions and

$\begin{array}{ll}x=(p-p_c)L^{1/\nu }& \left(4\text{d}\right)\end{array}$

as the scaling variable. Thus, the size dependence of O and OF gives us the exponents β/ν and γ/ν, respectively. The maximum of susceptibility also scales as L^γ/ν. Moreover, the value of p* for which the susceptibility has a maximum scales with L as

$\begin{array}{ll}{p}^{*}=p_c+b{L}^{-1/\nu }\mathrm{\text{ with }}b\approx 1.& \left(5\right)\end{array}$

Therefore, Equations (4c) and (5) give the exponent 1/ν. The effective dimensionality, D_eff, is given by the hyperscaling hypothesis

$\begin{array}{ll}2\beta /\nu +\gamma /\nu =D_{\mathrm{\text{eff}}}.& \left(6\right)\end{array}$

Monte Carlo simulations were made on (4, 6, 12) and (4, 8²) AL with various systems of size N = 768, 3,072, 12,288, 49,152, and 196,608 for (4, 6, 12) and N = 256, 1,024, 4,096, 16,384, 65,536, and for (4, 8²) AL. We use 2 × 10⁵ Monte Carlo steps (MCS) to make the system reach the steady state, and then the time averages are estimated over the next 3 × 10⁵ MCS. One MCS is accomplished after N attempts to update the opinions of agents i and j, considering the evolution Equations (1) and (2). The results are averaged over N_run (500 ≤ N_run ≤ 2, 000) independent simulation runs for each lattice and for given set of parameters (p, N).

3. Results and Discussion

Figure 2 displays the dependence of the opinion O, OF, and O₄ on the disorder parameter p, obtained from simulations on (4, 6, 12) and (4, 8²) AL with N ranging from N = 256 to 196,608 sites. The shape of O(p), OF, and O₄ curves, for a given value of L, indicate the occurrence of a second-order phase transition in the system. The phase transition occurs at the value of the critical disorder parameter p_c. Such critical disorder parameter p_c is estimated as the point where the curves of the Binder cumulant O₄ for different system sizes N intercept each other [23]. The corresponding value of O₄ is represented by $O_4^{*}$. Then, we obtain p_c = 0.086(3) and $O_4^{*}=0.59(4)$; p_c = 0.109(3) and $O_4^{*}=0.606(5)$ for (4, 6, 12), and (4, 8²) AL, respectively.

FIGURE 2

Figure 2. (Color on-line) O, OF, and O₄ vs. p, for sizes L = 2³, 2⁴, 2⁵, 2⁶ and 2⁷ and N = 12L² (for 4, 6, 12) AL (first line) and N = 4L² sites for (4, 8²) AL (second line).

The results obtained from Figures 3–6 and used in Figure 7 are summarized in Table 2.

FIGURE 3

Figure 3. Log-log plot of the dependence of the opinion O* = O(p_c) on the linear system size L. Fitting data, we obtained the estimation for the critical ratio β/ν.

FIGURE 4

Figure 4. Log-log plot of the OF at p_c vs. L for (4, 6, 12), and (4, 8²) AL. Fitting data, we obtained the estimation for the critical ratio γ/ν.

FIGURE 5

Figure 5. OF at p_χmax(N) vs. L for (4, 6, 12) and (8²), AL. Fitting data, we obtained another estimation for the critical ratio γ/ν.

FIGURE 6

Figure 6. Plot of ln |p_c(L) − p_c| vs. the linear system size L for (4, 6, 12) and (4, 8²) AL. Fitting data, we obtained the estimation for the critical exponent 1/ν.

FIGURE 7

Figure 7. (Color on-line) Data collapse of the opinion O, OF, and O₄ shown in Figure 2 for L = 32, 64, and 128 (4, 6, 12) and (4, 8²) AL. The exponent ratios used here were β/ν = 0.23(7), γ/ν = 1.43(5), and 1/ν = 0.60(3) for (4, 6, 12), and β/ν = 0.149(4), γ/ν = 1.56(4), and 1/ν = 0.94(4) for (4, 8²) AL.

The excellent curve collapses Figure 7 for distinct system sizes corroborates our estimated values for p_c and exponent ratios β/ν, γ/ν and 1/ν.

The resulting critical exponents and disorder parameters are collected in Table 2. One can also see that the exponent ratios β/ν, γ/ν, 1/ν are similar to MVM, Table 1, and are different from the Ising model contrary to the Grinstein's hypothesis [13]. They are different from β/ν = 0.125 and γ/ν = 1.75 obtained for a d = 2 lattices, but obey hyperscaling relation (within the error bars). Equation (6) yields effective dimensionality of systems D_eff = 1.89(6) for (4, 6, 12) and D_eff = 1.76(7) for (4, 8²). The effective dimensionalities of KCOD on our two AL are close to those for MVM on (4, 6, 12) AL (D_eff = 1.78(7)) and on (4, 8²) AL (D_eff = 1.83(6)). The results of simulations are collected in Table 2.

TABLE 1

Table 1. Critical temperatures, exponents and effective dimensionalities for MVM on (4, 6, 12) and (4, 8²) AL [22].

TABLE 2

Table 2. Critical parameter, exponents and effective dimension for KCOD model on (4, 6, 12) and (4, 8²).

4. Conclusion

We studied a non-equilibrium KCOD model through extensive Monte Carlo simulations on (4, 6, 12) and (4, 8²) AL. On these lattices, the KCOD model shows a second-order phase transition. Our Monte Carlo simulations suggest that the effective dimensionality D_eff is close to two, i.e., that hyperscaling relation 2β/ν + γ/ν ≈ 2 may be valid.

Finally, we remark that the critical exponents γ/ν, β/ν, and 1/ν for KCOD on (4, 6, 12) and (4, 8²) AL are similar to the MVM model on (4, 6, 12) and (4, 8²) AL [22], see Tables 1, 2. Therefore, this model not belong to the Ising universality class [24] and the hypothesis of Grinstein has been disproved.

Author Contributions

The author confirms being the sole contributor of this work and approved it for publication.

Conflict of Interest Statement

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

FL is grateful to Dietrich Stauffer for stimulating discussions and for critical reading of the manuscript. He also acknowledges the support the system SGI Altix 1350 the computational park CENAPAD, UNICAMP-USP, SP-BRASIL.

References

1. Galam S. Majority rule, hierarchical structures, and democratic totalitarianism: a statistical approach. J Math Psychol. (1986) 30:426–34; earlier Galam S, Gefen Y, Shapir Y. Sociophysics: a new approach of sociological collective behaviour. I. mean-behaviour description of a strike. J Math Sociol. (1982) 9:1–13. doi: 10.1080/0022250X.1982.9989929

CrossRef Full Text | Google Scholar

2. Galam S. Social paradoxes of majority rule voting and renormalization group. J Stat Phys. (1990) 61: 943–51. doi: 10.1007/BF01027314

CrossRef Full Text | Google Scholar

3. Fortunato S, Macy M, Redner S. Editorial. J Stat Phys. (2013) 151:1–8. doi: 10.1007/s10955-013-0703-2

CrossRef Full Text | Google Scholar

4. Galam S. The drastic outcomes from voting alliances in three-party democratic voting (1990 –> 2013). J Stat Phys. (2013) 151:46–68. doi: 10.1007/s10955-012-0641-4

CrossRef Full Text | Google Scholar

5. Galam S. Sociophysics: a personal testimony. Phys A (2004) 336:49–55. doi: 10.1016/j.physa.2004.01.009

CrossRef Full Text | Google Scholar

6. Stauffer D. A biased review of sociophysics. J Stat Phys. (2013) 151:9–20. doi: 10.1007/s10955-012-0604-9

CrossRef Full Text | Google Scholar

7. Galam S. Sociophysics: A Physicist's Modeling of Psycho-political Phenomena. New York, NY: Springer (2012). doi: 10.1007/978-1-4614-2032-3

CrossRef Full Text | Google Scholar

8. de Oliveira MJ. Isotropic majority-vote model on a square lattice. J Stat Phys. (1992) 66:273–81. doi: 10.1007/BF01060069

CrossRef Full Text | Google Scholar

9. Santos MA, Teixeira S. Anisotropic voter model. J Stat Phys. (1995) 78:963–70. doi: 10.1007/BF02183696

CrossRef Full Text | Google Scholar

10. Crochik L, Tomé T. Entropy production in the majority vote model. Phys Rev E (2005) 72:057103–7. doi: 10.1103/PhysRevE.72.057103

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Hasenbusch M. Monte Carlo studies of the three-dimensional Ising model in equilibrium. Int J Mod Phys C (2001) 12:911–1009. doi: 10.1142/S0129183101002383

CrossRef Full Text | Google Scholar

12. Binney JJ, Dowrick NJ, Fisher AJ, Newman MEJ. A Theory of Critical Phenomena. An Introduction to the Renormalization Group. Oxford: Clarendon Press (1992).

Google Scholar

13. Grinstein G, Jayaprakash C, He Y. Statistical mechanics of probabilistic cellular automata. Phys Rev Lett. (1985) 55:2527–30. doi: 10.1103/PhysRevLett.55.2527

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Biswas S, Chatterjee A, Sen P. Disorder induced phase transition in kinetic models of opinion dynamics. Phys A (2012) 391:3257–65. doi: 10.1016/j.physa.2012.01.046

CrossRef Full Text | Google Scholar

15. Deffuant G, Amblard F, Weisbuch G, Faure T. How can extremism prevail? A study based on the relative agreement interaction model. J Artif Soc Soc Simul. (2002) 5:1–26.

Google Scholar

16. Hegselmann R, Krause U. Opinion dynamics and bounded confidence models, analysis and simulation. J Artif Soc Soc Simul. (2002) 5:1–32.

Google Scholar

17. Sznajd-Weron K, Sznajd J. Opinion evolution in closed community. Int J Mod Phys C (2000) 11:1157–65. doi: 10.1142/S0129183100000936

CrossRef Full Text | Google Scholar

18. Costa Filho RN, Almeida MP, Andrade JS Jr, Moreira JE. Scaling behavior in a proportional voting process. Phys Rev E (1999) 60:1067–8. doi: 10.1103/PhysRevE.60.1067

CrossRef Full Text | Google Scholar

19. Araripe LE, Costa Filho RN, Herrmann HJ, Andrade JS Jr. Pluraty voting: the statistical laws of democracy in Brazil. Int J Mod Phys C (2006) 17:1809–13. doi: 10.1142/S0129183106010200

CrossRef Full Text | Google Scholar

20. Araripe LE, Costa Filho RN. Role of parties in the vote distribution of proportional elections. Phys A (2009) 388:4167–70. doi: 10.1016/j.physa.2009.06.023

CrossRef Full Text | Google Scholar

21. Mukherjee S, Chatterjee A. Disorder induced phase transition in an opinion dynamics model: results in 2 and 3 dimensions. Phys Rev E (2016) 94:062317–22. doi: 10.1103/PhysRevE.94.062317

CrossRef Full Text | Google Scholar

22. Welington F, Lima S. Nonequilibrium model on Archimedean lattices. Cent Eur J Phys. (2014) 12:185–91. doi: 10.2478/s11534-014-0435-1

CrossRef Full Text | Google Scholar

23. Binder K, Heermann DW. Monte Carlo Simulation in Statistical Physics. Berlin; Heidelberg: Springer Verlag (1988). doi: 10.1007/978-3-662-08854-8

CrossRef Full Text | Google Scholar

24. Temperley HNV. Two-dimensional ising models, in: Domb C, Green MS, editors. Phase Transitions and Critical Phenomena. Vol. 1 London: Academic Press (1972).