The KCOD Model on ( 3 , 4 , 6 , 4 ) and ( 3 4 , 6 ) Archimedean Lattices

Through Monte Carlo simulations, we studied the critical properties of kinetic models of continuous opinion dynamics on (3, 4, 6, 4) and (34, 6) Archimedean lattices. We obtain pc and the critical exponents’ ratio from extensive Monte Carlo studies and finite size scaling. The calculated values of the critical points and Binder cumulant are pc = 0.085(6) and O∗ 4 = 0.605(9); and pc = 0.146(5) and O∗ 4 = 0.606(3) for (3, 4, 6, 4) and (3 4, 6) lattices, respectively, while the exponent ratios β/ν, γ/ν and 1/ν are, respectively: 0.126(1), 1.50(7), and 0.90(5) for (3, 4, 6, 4); and 0.125(3), 1.54(6), and 0.99(3) for (34, 6) lattices. Our new results agree with majority-vote model on previously studied regular lattices and disagree with the Ising model on square-lattice.


Introduction
The study of the behavior of individuals in a society by physicists is known as sociophysics, having as the main contributor in this new research area Serge Galam who introduced the use of local majority rules to study voting systems as bottom-up democratic voting in hierarchical structures [1][2][3][4].Although sociophysics was rejected by some physicists in the eighties [5], it has today become an active field of research among physicists all over the world [3,6,7].
In this same context and based on the criterion of Grinstein et al. [8] (where a nonequilibrium model presenting up-down symmetry in two-state dynamic systems implies the same critical behavior (same universality class) as the equilibrium Ising model), Oliveira [9] proposed a nonequilibrium version of Ising model called majority vote model (MVM).On two-dimensional regular lattices, this presents a second-order phase transition with critical exponents β, γ, ν, as for [10,11] the equilibrium Ising model [12,13].
A multiagent model for opinion formation in society by modifying kinetic exchange dynamics studied in the context of income, money, or wealth distributions in a society where a spontaneous symmetry-breaking transition to polarized opinion states starting from nonpolarized opinion states was proposed by M. Lallouache et al. [15].
A model of continuous opinion dynamics (KCOD) was proposed by Biswas et al. [16] in 2012.In the KCOD model, the mutual interactions can be both positive and negative and a single parameter p denoting the fraction of negative interactions was considered in order to characterize the different types of distributions for the mutual interactions.Numerical simulations of the continuous version of this model 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)) (see also [17]).(3,4,6,4) and (3 4 , 6)  [14].For completeness, we cite data for Ising model on (4 4 ) as well [18].

MVM
(3, The KCOD model on square and cubic lattices (2D and 3D) was studied by Mukherjee and Chatterjee [19].Their numerical results indicate that the critical behavior of the KCOD model is the same as that of the Ising model in the corresponding dimensions.
Recently, C. Anteneodo and N. Crokidakis [20] studied a model of like KCOD model in the presence of a social temperature.The critical behavior of this model showed three different kinds of collective states (symmetric, asymmetric, and neutral) and nonequilibrium transitions between them (see also [21,22]).
In this work, we studied the KCOD on two Archimedean lattices-namely, (3,4,6,4) and (3 4 , 6)-through extensive Monte Carlo simulations.The topologies of (3,4,6,4), and (3 4 , 6) AL are presented 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 eleven 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 4 )), honeycomb is called (6 3 ), 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 4 , 6) AL.

Model and Simulations
The KCOD [16] model is defined as follows: A set of agents (individuals) with continuous opinion variables o i (t) is situated on every node of the (3,4,6,4) where the i, j interactions are pair-wise interactions between nearest neighbors, which implies no sum over the index j, and µ ij = µ ji are real random variables.In the above dynamics (Equation (1)), an agent i updates his opinion by interacting with agent j and is influenced by the mutual influence term µ ij .Here, j is selected randomly from one of the nearest neighbors.Unlike other models (such as Ising model and MVM) that present up-down symmetry [23], in the KCOD model the opinions are bounded (i.e., −1 ≤ o i (t) ≤ 1).If the opinion value of an agent becomes higher (lower) than +1 (−1), then it is made equal to +1 (−1) to preserve this bound.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].The ordering in the system is measured by the quantity O = | ∑ i o i |/N , the average opinion.
Changing the fraction p of negative interactions, one can observe a symmetry breaking transition between an ordered and a disordered phase below a particular value p c of the parameter p, the system orders (giving a non-zero, finite value of the order parameter O (opinion), defined in the following), while a disordered phase exists above p c (O = 0).
To study the critical behavior of the model, we are interested in the average opinion O, order parameter fluctuations OF, and the reduced fourth-order cumulant of the O (herein named as O 4 ), defined as where • • • stands for time averages, computed at the steady states.The results are averaged over the N run independent simulations.The above-mentioned quantities are functions of the disorder parameter p, and obey the finite-size scaling relations where ν, β, and γ are the usual critical exponents, f o , f o f , f o 4 (x) are the finite-size scaling functions with being the scaling variable.Therefore, from the size dependence of O and OF, we obtained the exponent ratios β/ν (O) and γ/ν (OF).The maximum value of susceptibility also scales as L γ/ν .Moreover, the value of p * for which OF has a maximum is expected to scale with the system size as Therefore, the relations (3c) and ( 4) may be used to get the exponent 1/ν.We also evaluate the effective dimensionality, D eff , from the hyperscaling hypothesis 2β/ν + γ/ν = D eff . ( Monte Carlo simulations were performed on (3,4,6,4) and (3 4 , 6) AL with various systems of size N = 384, 1536, 6144, 24, 576, and 98, 304 for (3,4,6,4) and (3 4 , 6) AL.It takes 2 × 10 5 Monte Carlo steps (MCS) to let the system reach the steady state, and then the time averages are calculated over the next 3 × 10 5 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 (1000 ≤ N run ≤ 2000) independent simulation runs for each lattice and for given set of parameters (p, N).

Conclusions
We studied a nonequilibrium KCOD model through extensive Monte Carlo simulations on (3,4,6,4) and (3 4 , 6) AL.On these lattices, the KCOD 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.
and (3 4 , 6) AL with N = 6L 2 sites.The opinion of an individual i at time t takes the values in the range [−1, +1], in a system of N agents.Here, the opinions change out of pair-wise interactions via mutual influences/couplings µ ij as:

Figure 2
displays the dependence of the opinion O, OF, and O 4 on the disorder parameter p, obtained from simulations on(3,4,6,4) and(3 4 , 6) AL with L ranging from L = 8 to L = 128.The shape of O(p), OF, and O 4 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 .This critical disorder parameter p c is estimated as the point where the curves of the Binder cumulant O 4 for different system sizes N intercept each other [24].The corresponding value of O 4 is represented by O * 4 .Then, we obtained p c = 0.085(6) and O * 4 = 0.605(9); p c = 0.146(5) and O * 4 = 0.606(4) for (3, 4, 6, 4), and (3 4 , 6) AL, respectively.

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

Table 1 .
Critical parameter (p c ), exponents, and effective dimension for majority vote model (MVM) on