Modeling and Simulation of Synchronous Threshold in Vent Collective Behavior

With the strengthening of the social contradiction, the outbreak of vent collective behavior tends to be frequent. The essence of vent collective behavior is emergence of synchronization. In order to explore the threshold of consensus synchronization in vent collective behavior, a mathematic model and a corresponding simulation model based on multi-agent are proposed. The results of analysis by mean field theory and simulation experiments show the following. (1) There is a threshold K c for consensus synchronization in global-coupling and homogeneous group, and when the system parameter K is greater than K c , consensus synchronization emerge. Otherwise the system cannot achieve synchronization. The conclusion is verified by further study of multiagent simulation. (2) Comparedwith the global-coupling situation, the process of synchronization is delayed in local-coupling and homogeneous group. (3) For local-coupling and heterogeneous group, consensus dissemination can achieve synchronization only when the effects of the parameters meet the threshold requirements of consensus synchronization.


Introduction
Synchronous behaviors formed by self-organization are ubiquitous in nature [1,2].A large number of physical and biological systems can be thought of as a synchronous system.The most striking visual example is given by the synchronous flashing of fireflies observed in summer season in some regions of South Asia [3], and more examples include chorusing crickets [4] and coupled laser arrays [5].One interesting common phenomenon in all these systems is their ability to synchronize: a large number of the individuals lock onto a common frequency.
In fact, the synchronous behaviors are widespread not only in the biological world, but also in society, such as the synchronization applause in music hall [6].Sometimes, some burst factors cause people to vent their emotions.Under some invisible "cohesion, " they participate in a temporary group by self-organization, and, then, synchronous behavior may be formed.We call the synchronous behavior as vent collective behavior [7].As a social phenomenon, how does the infection of emotions and consensus work in the group composed of independent individuals?And, through individuals' consensus transmission and behavior imitation, how does the synchronous behavior emerge?Is there a synchronous threshold in vent collective behavior?So far, there has been very little work on the above questions.This is due to the complexity of vent collective behavior.It is difficult to establish a model corresponding with the actual situation by quantitative analysis.In recent years, with the in-depth study of social complex systems, related topics begin to arouse numerous scholars' interests in different fields, such as sociology, psychology, computer science, and physics, but what they care about is different.For sociologists [8] and psychologists [9], they qualitatively analyze the problem from the political angle, computer scientists [10] are concerned about the simulation of collective escape behavior under the emergency, and physicists [11] are more willing to study the synchronization in complex network.
With the development of complex system science, simulation based on multiagents is more and more popular as a new research method [12][13][14].Some scholars begin to use it to study the consensus transmission in collective behaviors.The early simulation researches made individuals in collective behavior homogeneous [15].After that, researches on the (a)  0 = 0.3;  2 = 0.17 (b)  0 = 0.7;  2 = 0.17 (c)  0 = 0.5;  2 = 0.17 effects of various characteristics on collective behavior were developed [16].Recently, some scholars used continuous view dynamic model [17] to study the formation dynamics of consensus in collective behavior.They took the real number from 0 to 1 to represent the individual's opinions on social events.Through the simulation experiments, they observed rich evolution dynamic phenomenon.In fact, in the collective behavior, individuals learn from each other and then make the strategy choice in the next time step.In this process, each individual just can be seen as an agent who has ability to learn from others and "think" of its future strategy.And then, the group can be regarded as a complex system with multiagents.
Motivated by the above work [6,[15][16][17] and the discussions earlier, in this paper, we are concerned with the problem of synchronous threshold in vent collective behavior.In Section 2, we set up a model of consensus dissemination of vent collective behavior.The model is analyzed by mean field method [18] in Section 3, and the common conclusions of the synchronous threshold in vent collective behavior are given.For more complex situation of local-coupling and homogeneous group, multiagent simulation model based on sociology and psychology is proposed in Section 4. We will put all individuals in a two-dimensional grid space and define that each individual has heterogeneous conformity feature, attitude, and influence ability.The attitude of the individual to the social events is influenced by the whole group consensus or near neighbors' opinion.Through a series of simulation experiments in Section 5, we will prove that there is a threshold which induces the synchronous consensus in collective behavior, and the synchronous consensus is delayed in the local-coupling environment than global-coupling with the same parameters.For local-coupling and heterogeneous group, consensus dissemination can achieve synchronization only when the effects of the parameters meet the threshold requirements of consensus synchronization.

Consensus Dissemination Model
In vent collective behavior, most participants are not direct gainers who participate in collective behavior not due to benefits but just venting their emotions.Different from direct gainers who have direct interest in the event, they have no direct relationship with event, and they also do not know each other, so they are relatively independent.Nevertheless, each individual has been infected by group consensus to adjust its attitude of whether to participate in the behavior.Because vent collective behavior bursts out generally and has a short duration, we assume that the number of participants is , and each individual  has relatively stable Ω  neighbors for the great crowd density.Normally, when a person is faced with an event, he has the attitude which will decide his behavior strategies.So we assume that individual i's attitude is   () at time .Clearly the evolution of individual's attitude is influenced by its attitude in the previous stage and coupling effect between individuals, where the coupling effect between individuals contains two aspects: external impact on the individual and individual's acceptance of external impact.Firstly, we consider the external impact on the individual.Because each individual has a different number of neighbors in vent collective behavior, we set that the individual  has   neighbors.Since each individual's ability of influencing others to change attitudes is different, we set that   represents individual j's influence on others, which is different caused by the difference of abilities, such as expressive ability.Take that  obey the normal distribution function ℎ().Considering that the opinion of individuals is affected not only by their neighbors but also by the opinion of themselves, we express the individual j's influence on neighbor i as where  is the attitudes influence function.[  () −   ()] reflects the differences between individuals' opinion and their neighbors' opinion; the greater the difference the greater the impact on the individual and the converse is also true.When   >   , individual  is more inclined to participate in vent collective behavior than individual , and the converse is also true.Individual j's influence on individual  will be greater when   is greater.We use the mean value of all neighbors' influence to represent the combined influence on individual , namely, Secondly, we consider individual's acceptance for external impact.Although the individual will be influenced by group consensus to adjust its own attitude, the change degree of each individual is different.To characterize this feature, we set that   shows individual's ability of accepting influence.  is a random value of a normal distribution function ().
In addition, individual's conformity, namely, frequency of changing its attitude, is different for the external impact.We set that  denotes individual's conformity coefficient.The  obeys the normal distribution function ().Therefore, individual's attitude in vent collective behavior can update as the following formula: where   ,   , and   are all real numbers.Formula (3) shows that individual's attitude is influenced by its attitude in the previous stage, its conformity coefficient, and the average influence of its neighbors.Individual's attitude update rule considers not only individual's receptivity, but also individual's infective ability, so it can represent the coupling influence among individuals in vent collective behavior.

Synchronous Threshold in Vent Collective Behavior
In fact, because there are heterogeneous neighbors (localcoupling) and heterogeneous parameters in model (3), we cannot get the analytic solution of this complex model; to find the threshold of synchronization in vent collective behavior, let us consider one special case firstly.
In global-coupling group (  = ), we set  =     , where  is a constant.Consider that, in vent collective behavior, individuals maybe change their attitudes frequently.For convenience, we assume that  is a periodic function, and [  () −   ()] = sin[  () −   ()].So formula (3) is essentially the Kuramoto model which shows that there is a coupling strength threshold which can lead to synchronization of a set of oscillators with fixed amplitude mutually coupled by a 2 periodic interaction [19].In recent years, various extensions and generalizations models based on Kuramoto model have been analyzed more deeply [20].However, in the original form of this model, the individuals are global-coupling and parameter  is homogeneous.
By constructing the order parameter as [21] then formula (3) can be rewritten as where  is the average value of all individuals' attitudes.
Supposing the density function of  is (, , ) in a limited population, the density function (, , ) accords with the conditions [22] Formula (5) shows that when individuals' attitudes incline to synchronization,  is close to 1, and when individuals' attitudes incline to dissimilarity,  trails off.
Inspired by the above analysis, some people consider whether there is a threshold   which will make attitudes achieve synchronization.According to the mean field theory [18], they find that when  → 0 + , the threshold   = 2/(0) will make individuals' attitudes achieve synchronization [20,22].This is to say, in global-coupling and homogeneous parameters case, when  >   , consensus of vent collective behavior will achieve partial or complete synchronization.When  ≤   , consensus cannot achieve synchronization.Then, can the same conclusion be got for more complex situation just as local-coupling case?And when we consider the heterogeneous conformity and influence ability of individuals, what laws will formula (3) show?
Because there is no analytic solution with local-coupling and heterogeneity parameters case, in the rest of the sections, we will set up multiagent simulation model in Section 4 and verify the conclusions of Section 3 in Section 5.1; then, we will study the more complex case of the model by multiagent simulation experiments in Sections 5.2 and 5.3.

Multiagent Simulation Model
We use agent to represent an independent person with the ability to make decisions and make two-dimensional grid the agent's environment which generally take  line  column, such as 20 × 30 or 50 × 50.Every agent occupies one node.The individual's attitude is related to both its conformity coefficient and consensus of vent collective behavior.influence of neighbors.Generally,   is a constant.In the initial state,   randomly comes forth according to normal distribution function () and does not alter with time; (4) Φ = { 1 ,  2 , . . .,   } is influence ability feature space.It is used to reflect the individual's ability of influencing others.Generally,   is a constant.In the initial state,   randomly comes forth according to normal distribution function ℎ() and does not alter with time; (5)  = { 1 (),  2 (), . . .,   ()} is individual's attitude feature space, and () denotes individual's attitude for whether to participate in vent collective behavior, and it alters with time; (6)  = {1, 2, 3 . ..} is system clock, which is the basis of the evolution of simulation system.
In the process of multiagent modeling, how to set up the effective individual interaction rule is a key problem.In social physics, what behavior strategy individual takes in vent collective behavior is a kind of group imitation action, and it is also a kind of gradual learning process.At present, there are a lot of evolutionary algorithms, such as the local evolutionary stable strategy and replicated dynamic [23] and comparison study by using the social comparison theory [24].Actually, in vent collective behavior, how individual's attitude evolves is restricted by its conformity coefficient and neighbors' attitude.This restriction includes not only interacting imitation or learning but also other nonlinear couplings, such as emotion abreaction.This paper designs a kind of blended learning rules in terms of formula (3).In addition, to estimate whether the system achieves synchronization, the criterion is given as where Ẋ() is the average of ẋ  ().Whether  is close or equal to 1 explains whether the system achieves partial or complete synchronization.
Another notice is that () is a normal distribution function, namely,  ∼ ( 0 ,  2 ).When  becomes bigger, the conformity coefficient will be bigger, and individual is more inclined to change its attitude by influence of consensus.When  2 becomes bigger, individual's attitude goes to heterogeneity more clearly (Figure 1).In this paper, ℎ() and () have the same law with ().

Simulation and Results
In this paper, we use the NetLogo software to make a twodimensional lattice whose size is 15 × 15.Let   ∼ ( 0 ,  2 ); by taking  0 = 0.5 and  = 0.17, we guarantee that the most values of are between 0 and 1.Furthermore, we only take   ∈ [0, 1] in this paper, so  is almost always a positive value.When  = 0, the individual is very stubborn; when  = 1, the individual will change his opinion according to others completely.Assume that the initial attitude   (0) obeys uniform distribution   (0)∼Uniform (0, 1).The coupling effects between individuals' attitude evolve in accordance with (3).
In this paper, we will study three different cases which are global-coupling and homogeneous group, local-coupling and homogeneous group, and local-coupling and heterogeneous group.For global-coupling group, every individual is linked with all of the other individuals.However, the individual is only linked with its neighbors (e.g., nearest neighbors) in local-coupling group.We will simplify the parameters of   and   as constant in homogeneous group but set them to obey uniform distribution in heterogeneous group.
To verify the effect of established model and the result of Section 3, we will reveal the process of consensus synchronization in global-coupling and homogeneous group of vent collective behavior firstly.And, then, we change the condition of global-coupling to local coupling to find the laws of synchronous consensus.At last, we will find the deep laws for more complex situation of heterogeneous parameters.

Threshold of Consensus Synchronization in Global-
Coupling and Homogeneous Group.Figures 2(a) and 2(b) are the subfigures which show that synchronization criterion r and collective attitudes () change with time, respectively; Figures 2(c) and 2(d) are the histograms of the distribution of attitude () and the distribution of attitude changing ratio ẋ  (), respectively.From the four subfigures of Figure 2, we can see when K = 0.5, neither the attitude nor the attitude changing ratio of groups has not reached the synchronization status.
As can be seen from Figure 3, when K = 1, the attitude changing ratio of vent collective behavior has almost reached the synchronization status; then  ≈ 1.The attitude changing ratio is stabilized at about 0.5 to 0.6.But from Figure 3(c) we can see that there are still differences in group opinion; that is, group attitude is not consistent.
Setting  = 2, we find that not only the attitude but also the attitude changing ratio of groups is close to the synchronization status (Figure 4).With increasing the value of , we find that as  increases, the attitude and the attitude changing ratio tend to be more synchronous (Figure 5).
Figures 2-5 show that, with the increasing of , the system's capability of synchronization also increases.Figure 6 gives the results of  changing with the , where the simulation time is 1500, and every dot in the figure is the average value of 5 experiments.Corresponding with the conclusion in Section 3, there is a threshold of consensus synchronization in vent collective behavior   ≈ 0.6.When  >   , the system reaches synchronization.

Threshold of Consensus Synchronization in Local-Coupling and Homogeneous
Group.In reality, the individual is often not influenced by the whole group to change its attitude but is influenced only by their adjacent small groups.To reflect this fact, we further study the threshold of consensus synchronization in local-coupling group by simulation experiments.
We can see from Figure 7 that when K = 1, neither the attitude nor the attitude changing ratio of groups has not reached the synchronization status. Figure 8 shows the results of simulation system with  = 1.2.Although the attitude changing ratio of vent collective behavior has reached the synchronization status, the group attitude is not consistent.Compared with Section 5.1, we find that, with the same parameters, the process of consensus synchronization is delayed to the global-coupling situation.
Figure 9 further reveals the above law.In global-coupling case, the system can basically be synchronized when K > 0.7, but, in the case of local-coupling, only when K > 1.2, the system can achieve synchronization.Obviously, this result is not corresponding with [25][26][27].This is because there are many factors that may influence the synchronization of consensus, such as the coupling strength K and the number of neighbors.Though Strogatz and Mirollo proved that synchronization was not possible in the thermodynamic limit [26] and similar conclusions were reached by Daido [27], they also point out that the system will be an attainable synchronization for a finite  provided  is sufficiently large.Different from the previous studies [25][26][27], in this paper, the system is two-dimensional lattice but not one-dimensional chain, the probability distribution of natural frequencies  is also somehow different, and the local-coupling is not the nearest neighbors but is more widely ranging.We think all of these differences may result in the different results.

Threshold of Consensus Synchronization in Local-Coupling and Heterogeneous
Group.Finally, we study the threshold of consensus synchronization in the conditions of local coupling, and the individual's influence ability  and acceptance influence ability  are heterogeneous.
Figure 10 shows the results of heterogeneous parameters for individuals.Set that  and  satisfy the normal distribution whose mathematical expectation is 0.5 and variance is 0.17.Obviously, as can be seen from Figure 10, the system cannot achieve synchronization.The reason is that  and  are too small to affect the system to reach the threshold of consensus synchronization.
Adjusting the mathematical expectation to 1.5 and maintaining variance to be unchanged, we can get the results shown as in Figure 11.The system achieves synchronization now.The reason for this result is that, with the increasing of the mean value of  and , their effects are also increased.When the value of  and  reaches the threshold which can achieve synchronization, the consensus synchronization emerges.

Conclusions
When the social contradiction and public animosity reach some degree, an emergency usually causes the outbreak of vent collective behavior.The essence of vent collective behavior is emergence of synchronization.Usually, synchronous behavior will occur with a certain condition.So, what is the synchronous condition of vent collective behavior?To answer this question, we propose a mathematic model corresponding to the reality.Using the mean field theory and multiagent simulation method, we study three situations which are global-coupling and homogeneous group, local-coupling and heterogeneous group, and local-coupling and homogeneous group.Along with the increasing of complexity of model, we carry out series simulation experiments.The results show that there is a threshold   for consensus synchronization in global-coupling and homogeneous group and the threshold   is related to the original distribution of conformity coefficient.When K > K c , consensus of vent collective behavior will achieve partial or complete synchronization.
When  ≤   , consensus cannot achieve synchronization.For local-coupling and homogeneous group, the process of consensus synchronization is delayed compared with the global-coupling situation for the same parameters.For localcoupling and heterogeneous group, with the increasing of the mean value of  and , their effects are also increased; when the value of  and  reaches the threshold which can achieve synchronization, the consensus synchronization emerges.

Figure 1 :Figure 2 :Figure 3 :
Figure 1: Influence of parameters in normal distribution on individual's attitude; 0 ≤  ≤ 1.Light color expresses individual is a standpatter (more inclined to keep its attitude).Dark color expresses individual is capricious (more inclined to change its attitude).The influence of mathematical expectation on individual's initial conformity is shown in (a) and (b).The influence of variance on individual's conformity is shown in (c) and (d).

Figure 4 :
Figure 4: Results of simulation system with  = 2. (a) is the subfigure which shows synchronization criterion  changes with time; (b) is the subfigure which shows group attitudes () versus time; (c) is the histogram of the final distribution of agents versus group attitudes (); (d) is the histogram of the final distribution of agents versus attitude changing ratio ẋ  ().

Figure 5 :
Figure 5: Results of simulation system with  = 5.(a) is the subfigure which shows synchronization criterion  changes with time; (b) is the subfigure which shows group attitudes () versus time; (c) is the histogram of the final distribution of agents versus group attitudes (); (d) is the histogram of the final distribution of agents versus attitude changing ratio ẋ  ().

Figure 6 :
Figure 6: When group is global-coupling, the threshold   is about 0.6 for consensus synchronization in vent collective behavior.

Definition 1 .Figure 7 :
Figure 7: Results of simulation system with  = 1.(a) is the subfigure which shows synchronization criterion  changes with time; (b) is the subfigure which shows group attitudes () versus time; (c) is the histogram of the final distribution of agents versus group attitudes (); (d) is the histogram of the final distribution of agents versus attitude changing ratio ẋ  ().

Figure 8 :
Figure 8: Results of simulation system with  = 1.2.(a) is the subfigure which shows synchronization criterion  changes with time; (b) is the subfigure which shows group attitudes () versus time; (c) is the histogram of the final distribution of agents versus group attitudes (); (d) is the histogram of the final distribution of agents versus attitude changing ratio ẋ  ().

Figure 9 :
Figure 9: When group is local-coupling, the threshold   is about 1.2 for consensus synchronization in vent collective behavior.

Figure 10 :
Figure 10: The simulation results with heterogeneous parameters, where  and  satisfy the normal distribution whose mathematical expectation is 0.5 and variance is 0.17.(a) is the subfigure which shows synchronization criterion  changes with time; (b) is the subfigure which shows group attitudes () versus time; (c) is the histogram of the final distribution of agents versus group attitudes (); (d) is the histogram of the final distribution of agents versus attitude changing ratio ẋ  ().

Figure 11 :
Figure 11: The simulation results with heterogeneous parameters, where  and  satisfy the normal distribution whose mathematical expectation is 1.5 and variance is 0.17.(a) is the subfigure which shows synchronization criterion  changes with time; (b) is the subfigure which shows group attitudes () versus time; (c) is the histogram of the final distribution of agents versus group attitudes (); (d) is the histogram of the final distribution of agents versus attitude changing ratio ẋ  ().