Social Conflicts Studied by Statistical Physics Approach and Monte Carlo Simulations

Statistical physics models of social systems with a large number of members, each interacting with a subset of others, have been used in very diverse domains such as culture dynamics, crowd behavior, information dissemination and social conflicts. We observe that such models rely on the fact that large societal groups display surprising regularities despite individual agency. Unlike physics phenomena that obey Newton’s third law, in the world of humans the magnitudes of action and reaction are not necessarily equal. The effect of the actions of group n on group m can differ from the effect of group m on group n. We thus use the spin language to describe humans with this observation in mind. Note that particular individual behaviors do not survive in statistical averages. Only common characteristics remain. We have studied two-group conflicts as well as three-group conflicts. We have used time-dependent Mean-Field Theory and Monte Carlo simulations. Each group is defined by two parameters which express the intra-group strength of interaction among members and its attitude toward negotiations. The interaction with the other group is parameterized by a constant which expresses an attraction or a repulsion to other group average attitude. The model includes a social temperature T which acts on each group and quantifies the social noise. One of the most striking features is the periodic oscillation of the attitudes toward negotiation or conflict for certain ranges of parameter values. Other striking results include chaotic behavior, namely intractable, unpredictable conflict outcomes.


Introduction
Statistical physics studies large systems of interacting particles using their microscopic properties [1]. Statistical physics models are suitable for investigating social systems which are composed of a large number of individuals who interact with each other under a social ambiance which can be associated with the temperature in physics. These models have been used to study complex social systems such as culture dynamics, crowd behavior, information dissemination and social conflict. This is possible [2,3] because the average over large societal groups washes away individual particularities and retains only common characteristics at the end. We see surprising behavior regularities in opinion surveys despite individual uniqueness. Many discussions have been found in the literature on the critique of oversimplification of social dynamics of sociophysics [4]. 2 of 9 In particular, social conflicts have been subject of investigations over the years [5][6][7][8][9]. Many applications to concrete situations such as the Brexit referendum [10,11], the US election in 2016 [10][11][12], the Serbia-Herzegovina-Croatia election in 2018 [13] have been discussed. This paper is organized as follows. Section 2 is devoted to the two-group dynamics with mean-field analysis and Monte Carlo simulations. Section 3 shows the case of three-group dynamics. Conclusion is given in the last section.

Two-Group Dynamics
We use a spin S i to represent the opinion of an individual i in a group. As said in the Introduction, this spin can be discrete or continuous, depending on the problem and the model we wish to make for that problem.

Model: Time-Dependent Mean-Field Theory
On the Renyi-Erdos (equivalent neighbor) network, the mean of preferences s of each group is proportional to the exponential of the intensity of interactions (negative energy): We introduce a lag time as we assume the preference s at time t + 1 interacts with the averages S A and S B at an earlier time t. The time is measured in units of the delay time. The sums on the right hand sites of Equations (1) and (2) involve the Brillouin function [1]: Equations (2)-(9) can be written as: An analysis of the linearized Equations (4) and (5), valid for small S A and S B , gives the region of the parameter space where an ordered phase (|S A | > 0, |S B | > 0) can exist beside the disordered phase (S A = S B = 0): On the manifold of Equation (6) both eigenvalues are equal to 1. The eigenvalues can become complex provided K 12 * K 21 < 0. In such cases the average S A and S B exhibit oscillations as function of time. The period of such oscillations is given by where λ R and λ I are the real and imaginary parts of the eigenvalue. When the eigenvalue absolute value is larger than unity the average S evolves in time to a non-zero value while when the eigenvalue absolute value is less than unity the average S evolves in time to zero. We show some results in the following with a few typical sets of parameters ( Figure 1).

Monte Carlo Simulation
We use the Metropolis algorithm [14] to simulate the model presented in the previous subsection: We calculate the energy of an individual at the time t using the state of the other individuals of his/her group and the average field of the other group at time t. We update his/her state according to the Metropolis algorithm, before taking another individual for updating until all individuals of two groups are updated. This achieves on Monte Carlo step. Unlike problems of statistical physics at equilibrium, here we would like to follow the individuals with time evolution. We need just to equilibrate the groups at T separately (this take times), and turn on the interaction between them. We follow how each group evolves with t.
We change the notations here: In the mean-field theory above, the "social temperature" has been included in the definitions of J 1 , J 2 , K 12 and K 21 . Here we explicitly express it by J 1 = J A /T, J 2 = J B /T, The size of N = 1600 individuals is used for each group with periodic boundary conditions so that individuals being at the borders of the group have the same number of neighbors. This order of size is what is used in opinion surveys so that a precision of a few percents since the uncertainty is proportional to 1/ √ N near the change of phase.
As an example we take the case with an antisymmetric −K AB = K BA . The two groups in the absence of interaction and with interaction at long-time average are shown in Figure 2.
The time evolution of the two groups' stances is shown in Figure 3 at three social temperatures: At low T two groups are stable with opposite stances, at intermediate T their stances oscillate (periodic change of stance), and at high T they oscillate chaotically. These results have been discussed in relation with the Brexit referendum [10,11].

Mean-Field Results
where j n = J n /T and k n,m = K n,m /T for n,m = 1, 2, 3. We use units such that k B = 1. Note that in Equations (9)-(11), s at time t + 1 interacts with the averages s 1 , s 2 and s 3 evaluated at an earlier time t. Here time is measured in units of the delay time. The sums on the right hand sides of Equations (9)-(11) give This model depends on 3 values M, 3 values of the intra-network couplings J, and 6 values of the inter-network couplings K. We consider next only M ≡ M 1 = M 2 = M 3 = 3 as in the above two-group case.
At lower temperatures T = 1.0 the oscillations are sustained and the three groups' attitudes are synchronized exhibiting the same period, (see Figure 4). The trajectory in the (s 1 , s 2 , s 3 ) space evolves in the long run to a closed-loop attractor shown in Figure 4. Lowering the temperatures increases the period of oscillations, and the closed-loop attractor begins to fragment (Figures 6 and 7).  At lower temperatures T = 1.0 the oscillations are sustained and the three groups' attitudes are synchronized exhibiting the same period, (see Figure 5). The trajectory in the (s 1 , s 2 , s 3 ) space evolves in the long run to a closed-loop attractor shown in Figure 5. Lowering the temperatures increases the period of oscillations, and the closed-loop attractor begins to fragment (Figures 6 and 7). Lowering the temperatures increases the period of oscillations, and the closed-loop attractor begins to fragment (Figures 6 and 7). We argue that the higher temperature attractor represents the essence of acute intractability: there is no single point at which the conflict settles, but rather a never-ending (non-sequential) cycling occurs among possible outcomes on the attractor.
At lower temperatures the attractor fragments into a discreet number of fixed points. The system still cycles among them but the discrete configuration corresponds to a lower degree of intractability than the continuous attractor ( Figure 4).

Monte Carlo Simulation
Monte Carlo simulations have been carried out as described in Section 2 but with three groups. An example is shown in Figure 8. For each group, there is a social temperature T C beyond which the stance of a group is lost, namely s = 0. We see that T is not the same for three groups: the higher the We argue that the higher temperature attractor represents the essence of acute intractability: There is no single point at which the conflict settles, but rather a never-ending (non-sequential) cycling occurs among possible outcomes on the attractor.
At lower temperatures the attractor fragments into a discreet number of fixed points. The system still cycles among them but the discrete configuration corresponds to a lower degree of intractability than the continuous attractor ( Figure 5).
At very low temperatures the attractor collapses into a single point that corresponds to the ordered phase of the static model. Because of the symmetry S → −S there is another static solution where all the attitude values are replaced by their negatives.

Monte Carlo Simulation
Monte Carlo simulations have been carried out as described in Section 2 but with three groups. An example is shown in Figure 8. For each group, there is a social temperature T C beyond which the stance of a group is lost, namely s = 0. We see that T C is not the same for three groups: The higher the intra-group interaction J, the higher T C threshold. At very low temperatures the attractor collapses into a single point that corresponds to the ordered phase of the static model. Because of the symmetry S → −S there is another static solution where all the attitude values are replaced by their negatives.

Monte Carlo Simulation
Monte Carlo simulations have been carried out as described in Section 2 but with three groups. An example is shown in Figure 8. For each group, there is a social temperature T C beyond which the stance of a group is lost, namely s = 0. We see that T C is not the same for three groups: the higher the intra-group interaction J, the higher T C threshold. As in the mean-field calculation above, an individual in a given group interacts at the time t + 1 with the average of the action field created by the other groups at the earlier time t. The only difference from the mean-field calculation is the short-range interaction considered in the Monte Carlo simulation. We will see that the results differ in some important aspects. As in the mean-field calculation above, an individual in a given group interacts at the time t + 1 with the average of the action field created by the other groups at the earlier time t. The only difference from the mean-field calculation is the short-range interaction considered in the Monte Carlo simulation. We will see that the results differ in some important aspects.
Once the equilibrium is reached for each group, we turn on the interactions between groups at time t.
An example at low T is shown in Figure 8 where the inter-group interactions may or not destroy the order of a group. We have chosen the interaction strengths and signs in the example below to illustrate the case of Serbia-Herzegovina-Croatia 2018 election [13].
At higher T, the order of each group is weakened. The inter-group interactions cause the groups' stances to oscillate widely without periodicity as also seen in the long-range mean-field results above. We observe that at times the stronger group 2 dominates the other two. This pattern reflects the level of intractability of the three-group conflict simulated here, consistent with the longer-term mean-field results. While the conflict is intractable at all the temperatures of Figure 9, at the lower temperature (corresponding to a stable context) the groups are 'stuck' in predictable ways (see Figure 9a); as the context gets heated, the three-group system cycles unpredictably through various stages (Figure 9b,c).
These results have been discussed in the case of Serbia-Herzegovina-Croatia election in 2018. The parameters used in Figure 9 correspond to the situation of this three-group competition [13].
(corresponding to a stable context) the groups are 'stuck' in predictable ways (see Figure 9a); as the context gets heated, the three-group system cycles unpredictably through various stages (Figure 9b and c).
These results have been discussed in the case of Serbia-Herzegovina-Croatia election in 2018. The parameters used in Figure 9 correspond to the situation of this three-group competition [13].

Conclusions
We have shown in this paper that statistical physics is a suitable tool to study social conflicts. In particular, we have outlined the Mean-Field Theory and the Monte Carlo simulation for the two-and three-group conflicts. We have shown the importance of the inter-group interaction which is at the origin of the fluctuation of the group stances if the social temperature is in a region of social instability. Intractability of conflicts is also seen with an appropriate choice of those parameters.