Research Article | Open Access

Bartłomiej Nowak, Katarzyna Sznajd-Weron, "Homogeneous Symmetrical Threshold Model with Nonconformity: Independence versus Anticonformity", *Complexity*, vol. 2019, Article ID 5150825, 14 pages, 2019. https://doi.org/10.1155/2019/5150825

# Homogeneous Symmetrical Threshold Model with Nonconformity: Independence versus Anticonformity

**Academic Editor:**Yong Xu

#### Abstract

We study two variants of the modified Watts threshold model with a noise (with nonconformity, in the terminology of social psychology) on a complete graph. Within the first version, a noise is introduced via so-called independence, whereas in the second version anticonformity plays the role of a noise, which destroys the order. The modified Watts threshold model, studied here, is homogeneous and possesses an up-down symmetry, which makes it similar to other binary opinion models with a single-flip dynamics, such as the majority-vote and the -voter models. Because within the majority-vote model with independence only continuous phase transitions are observed, whereas within the -voter model with independence also discontinuous phase transitions are possible, we ask the question about the factor, which could be responsible for discontinuity of the order parameter. We investigate the model via the mean-field approach, which gives the exact result in the case of a complete graph, as well as via Monte Carlo simulations. Additionally, we provide a heuristic reasoning, which explains observed phenomena. We show that indeed if the threshold , which corresponds to the majority-vote model, an order-disorder transition is continuous. Moreover, results obtained for both versions of the model (one with independence and the second one with anticonformity) give the same results, only rescaled by the factor of 2. However, for the jump of the order parameter and the hysteresis is observed for the model with independence, and both versions of the model give qualitatively different results.

#### 1. Introduction

Models of opinion dynamics are among the most studied models of complex systems [1–4]. This is not surprising, because they can be treated as a zero-level approach to various more complex social processes, including polarization of opinion [5–7], diffusion of innovation [8–10], or political voting [11–13]. In most of these models, public opinion is formed as an outcome from individual opinions of mutually interacting agents. Particularly interesting is a subset of the binary opinion models, including the voter model [14], the majority-vote model [14–16], the Galam model [17], the Sznajd model [18], the Watts threshold model [19], the -voter model [20], or the threshold -voter model [21, 22]. All these models belong to the broader class of binary-state dynamics, likewise the kinetic Ising models [23]. The binary decision/opinion framework is not only attractive from physicist’s point of view but also natural in the social sciences [19, 24].

In this paper, we will focus on a particular subclass of the binary opinion models with a single-flip dynamics, which means that one agent at most can change her/his state in a single update [25]. Such an updating scheme is used within the voter model, the majority-vote model, the Watts threshold model, and the -voter model. There are several common features possessed by all these models, at least in their original formulations. Within these models, we have the following:(1)We consider individuals that are tied to the nodes of some graph. Each node of a graph is occupied by exactly one agent.(2)Each individual is described by the dynamical binary variable , that represents an opinion on a given subject (yes/no, agree/disagree, etc.) at given time . Such a variable reminds an Ising spin and therefore wording, “individual”, “agent”, “voter” and “spin”, is used interchangeably; is often represented by , whereas is by .(3)Interactions between agents are local; i.e., they take place only if two agents are directly linked.(4)At each elementary update a single spin is randomly chosen and it can flip to the opposite direction with a probability that depends on the model’s details.(5)Conformity, i.e., an act of matching opinions, attitudes, beliefs, and/or behaviors to the certain group of influence, is the main type (often the only type) of the social response.(6)Agents are memoryless, which means that opinion of a given agent at time depends only on her/his own opinion at the previous time step and opinions of her/his neighbors also at the previous time step .

The differences between models depend mainly on the condition under which the conformity takes place. In some models all neighbors of a given target agent influence her/him (this applies to the majority-vote model [14, 16] or the Watts threshold model [19]), whereas in others only a certain group of influence is chosen from the neighborhood (e.g., in the -voter model [20, 21, 26]). In the -voter model unanimity of opinion is needed to influence a voter [20], whereas in the majority-vote model absolute majority is sufficient [16]. The last and particularly important difference consists of the presence or the absence of the up-down symmetry. The linear voter, the -voter model, and majority-vote are symmetrical; i.e., they are invariant to the swap of state labels, as the Ising model without the external field. However, the Watts threshold model is not symmetrical in its original formulation, because it was introduced as a model of innovation diffusion and flips from to were forbidden [19, 27]. Moreover, the Watts threshold model is heterogeneous, even on homogeneous graphs; i.e., each agent is characterized by its individual threshold needed for conformity. However, modification of the model to make it symmetrical and homogeneous (in a sense that all agents are characterized by the same threshold) is straightforward and will be investigated in this work.

One may ask what is the motivation to modify the Watts threshold model into the symmetric case? First of all, opinion dynamics concerns not only asymmetrical problems, as the diffusion of innovations, but many symmetric or almost symmetric issues, such as voting to one of two political parties and choosing one of two products on the duopoly market. Although the Watts threshold model was originally introduced to model the diffusion of innovation, the main idea of the threshold is also very natural in the broader context. It has been shown in many social experiments that simple majority () may not be sufficient to convince people; for short review see [28]. The second reason for the modification of the Watts threshold model comes from the basic research. Most of the binary opinion models within sociophysics have the yes-no (up-down) symmetry, so to make a comparison with them it is necessary to deal with the symmetric version of the Watts threshold model. In fact, the second reason was our main motivation. We wanted to understand the nature of the phase transitions observed within models of binary opinions with a single-flip dynamics and up-down symmetry.

Because all binary models, mentioned above, have been extensively investigated for years, many modifications and extensions of their original formulations have been proposed; a short review on modifications of the majority-vote model can be found in [29], on the Watts threshold model in [30], on the -voter model in [31], and on the Galam model comprehensive review written by the author of the model [32, 33]. Among many extensions, going into different directions, the introduction of an additional type of the social response was particularly interesting from the point of view of social/psychological sciences, as well as the theory of nonequilibrium phase transitions. This new type of social response is called nonconformity and can take one of two possible forms: independence (resisting influence) or anticonformity (rebelling against influence) [21, 34, 35]. In the first case, the situation is evaluated independently of the group norm, which means that a state of a given spin is not affected by its neighborhood. In the second case, a voter is influenced by the others but takes the position that is opposite to the group of influence. Therefore it is said that anticonformity and conformity are opposites at the operational level but at the same time similar at the conceptual level, because both indicate behavior that has been influenced by the source [35].

It is clear that conformity increases agreement (ferromagnetic order) in the system, whereas both types of nonconformity act against consensus. As a result of this competition, an order-disorder phase transition emerges. Interestingly, the type of the phase transition (continuous or discontinuous) may depend on the type of nonconformity. For example, it has been shown that within the -voter model with anticonformity only continuous phase transitions are possible, whereas for the -voter model with independence tricriticality (a switch between continuous and discontinuous phase transition) appears [21, 26, 36–38]. The majority-vote model contains conformity and anticonformity in its original formulation, but recently an additional noise, in the form of independence, has been introduced [29, 39]. It has been shown that the presence of an additional noise does not affect the type of the phase transition, which remains continuous independently of the network structure.

The question that naturally arises here is “which factor is responsible for the discontinuous phase transition within the -voter model, since we do not observe the analogous phenomenon within the majority-vote model?” We believe that investigating the modified version of the threshold model on the complete graph could help to understand this phenomenon. Due to our knowledge, the role of a noise has not been explored yet within the Watts threshold model [19, 30]. Therefore in this paper we will introduce two versions of the model, analogously as it was done for the -voter model, one with independence and the second one with anticonformity.

The paper is organized as follows. In the next section, we describe original versions of binary opinion models with single-flip dynamics. Then we present model’s extensions, which consist of introducing the noise into the models and we describe briefly the results that show how this noise impacts phase transitions. In the following subsection, we modify the original Watts threshold model to make it symmetrical and homogeneous, which makes it comparable to other binary opinion models with single-flip dynamics. Subsequently, we propose two versions of the symmetrical, homogeneous Watts threshold model, one with independence and the second one with anticonformity. Then, we analyze the model on the complete graph, which corresponds to the mean-field approach. We compare results obtained within Monte Carlo simulations results with those obtained via analytical treatment. Moreover, we provide a heuristic explanation of the obtained results. Finally, we discuss results in the context of other binary opinion models with single-flip dynamics.

#### 2. Methods

##### 2.1. Binary Opinion Models with a Single-Flip Dynamics

The general framework of all binary opinion models with a single-flip dynamics has been described above so we will not repeat it here. Instead, we present updating rules that define the dynamics of models within this class. The most extensively studied among all is the linear voter model [14]. On the other hand, the linear voter model is a special case of the more general -voter model [1] and thus we will not discuss it separately. The dynamics of the original -voter model is the following:(1)At a given time , choose one voter at random, located at site .(2)Choose randomly neighbors of site from its neighbors, where is degree of a node . In the original formulation and in some later versions repetitions were allowed to make the model universal (arbitrary value of on the arbitrary graph is possible) [20, 40–42]. However, in many other papers repetitions were forbidden [21, 26, 28, 36].(3)If all neighbors have the same opinion, the spin at site takes the same state as neighbors.(4)Otherwise, i.e., in lack of unanimity, spin at site can flip to the opposite direction with probability . In most of later modifications [21, 26, 28, 36, 40–42] and here we also refer to this case.(5)Time is updated .

In [43] two types of noise (interpreted as nonconformity) have been introduced to the model, but not simultaneously. Initially two versions of the model have been introduced: one with independence and the second one with anticonformity. In each of these models nonconformity (independence or anticonformity) takes place with probability , whereas with complementary probability agent conforms. Summarizing, the algorithm of a single step is as follows:(1) At a given time , choose one voter at random, located at site .(2) Update the opinion :(i)Model* I* (with Independence)(a)With probability , an agent changes opinion independently, i.e., she/he changes opinion to the opposite one with probability .(b)With probability , an agent conforms, i.e., if all agents, randomly chosen from all neighbors of site , are in the same state, then the voter at site takes the same position as those agents.(ii)Model* A* (with Anticonformity)(a)With probability , an agent anticonforms, i.e., acts against a group of influence, i.e., if all agents, randomly chosen from all neighbors of site , are in the same state, then the voter at site takes the opposite position to those agents.(b)With probability , an agent changes opinion as in Model I.(3) Time is updated .

The generalized versions of the model that consists of a threshold [21, 22] and two types of nonconformity simultaneously [28] have been also introduced. In the case of the threshold -voter model, only among neighbors have to share the same opinion in order to influence a voter.

The -voter model with independence was studied on the complete graph [21, 28, 43], as well as on various complex networks [36, 37], whereas the -voter model with anticonformity only on the complete graph. It has been shown that within Model A only continuous phase transitions are possible for all , whereas within Model I both types of phase transitions appear; for there is continuous phase transition and for transition is discontinuous. It has occurred that the tricritical point , even if can take noninteger values [38].

It has been also shown that in the case of the threshold -voter model there is a critical threshold that decreases with , above which discontinuous phase transitions are possible [21, 28]. In the most general case, when independence and anticonformity are introduced simultaneously, it occurs that is monotonically decreasing function of and is always greater than , even for very large [28]. Although the analytical form of has not been found, it was predicted that for the critical value . It means that the absolute majority is not sufficient for the discontinuous phase transition, even if for very large .

Another model with a single-flip dynamics, which has been analyzed in the presence of a noise, is the majority-vote model [14]. The dynamics of the original majority-vote model is as follows:(1)At a given time , choose one voter at random, located at site .(2)With probability a spin at site adopts the minority sign of neighboring spins.(3)With complementary probability a spin at site adopts the majority sign of neighboring spins.(4)Time is updated .

As seen from the above description, within the original majority-vote model anticonformity takes place with probability and conformity with , similarly as within the -voter model with anticonformity. Recently an additional noise has been introduced into the majority-vote model: with probability a voter acts independently, i.e., analogously like within the -voter model with independence, and flips randomly to the opposite direction with probability 1/2. With complementary probability the original rule is applied [29, 39]. The model has been investigated on the square lattice [29], as well as on several graphs, including homogeneous and heterogeneous structures, through the mean-field calculations and the Monte Carlo simulations [39]. It has been shown that only continuous phase transitions are possible within this model, similarly as for the model without additional noise. Discontinuous phase transitions do not appear even for highly connected networks. This result is consistent with the result obtained for the generalized threshold -voter model [21, 28], as described above.

##### 2.2. Symmetrical, Homogeneous Watts Threshold Model

In computational sociology, a particularly popular class of models describing the spread of innovation/idea/behavior are threshold models, based on the idea introduced by Granovetter [24]. Probably the simplest among them is the Watts threshold model [19]. The updating rule within its original formulation is as follows:(1)At a given time , choose one voter at random, located at site .(2)An agent at site is influenced by its neighbors. If at least a threshold fraction of its neighbors are in state , then an agent adopts this state; otherwise, nothing happens.(3)Time is updated .

There are two characteristic features that make the model different from other models described in the previous subsection. The first visible difference is heterogeneity. In other models it is introduced only by the heterogeneity of a graph; here, agents possess individual thresholds. Originally, each agent is assigned a threshold that is drawn at random from a given probability distribution function (PDF). Of course, as a special case, we can choose a one-point PDF, which takes the value equal to one at , and zero otherwise. Within such a formulation is an external (control) parameter of the model.

Another difference between Watts threshold model and other models, presented above, is the lack of the up-down symmetry. An agent who ones adopted cannot go back to an unadopted state. However, we can easily modify the model to make it symmetrical, in the following way:(1) At a given time , choose one voter at random, located at site .(2) An agent at site is influenced by all neighbors:(a)if at least a threshold fraction of its neighbors are in state , then an agent takes state 1, else(b)if at least a threshold fraction of its neighbors are in state , then an agent takes state -1, else(c)an agent remains in its old state.(3) Time is updated .

It should be noticed that within the above definition, the rules are defined unambiguously only for . To clarify this, let us give here an example. Imagine that and at a given time step . It means that the ratio of the positive opinions is equal to and the ratio of the negative opinions is equal to . Because we have both conditions: (a) a threshold fraction of neighbors are in state and (b) a threshold fraction of neighbors in the state are fulfilled, which means that there is no unambiguous choice. Of course one could think about another model in which this ambiguity could be solved by introducing the probabilistic rule. However, in such a case interpretation of what is conformity and what is anticonformity would be far less clear. Therefore here we consider only . This can be interpreted also as a supermajority (or a qualified majority), whereas corresponds to the simple majority rule.

Now we are ready to introduce two versions of the model with nonconformity, one with Independence (Model* I*) and the second one with Anticonformity (Model* A*), analogously as it was done for -voter model. Within Model A an agent conforms or anticonforms if the fraction of neighboring spins having the same state is larger than a fixed threshold . In the conformity case (which takes place with probability ), an agent follows the opinion of the group of influence, whereas in the case of anticonformity (which takes place with probability ) she/he takes the opposite opinion to the group, as in the majority-vote or in the -voter model with anticonformity. Within Model I, instead of anticonformity, independence takes place with probability : an agent flips to the opposite state with probability .

The algorithm of a single update is as follows:

At a given time , choose one voter at random, located at site .

Update the opinion :(i)Model I (with Independence)(a)With probability , an agent acts independently; i.e., she/he changes an opinion to the opposite one with probability .(b)With probability , an agent conforms to its neighbors:(i)if at least a threshold fraction of its neighbors are in state , then an agent takes state , else(ii)if at least a threshold fraction of its neighbors are in state , then an agent takes state , else(iii)an agent remains in its old state.(ii)Model A (with Anticonformity)(a)With probability , an agent anticonforms to its neighbors:(i)if at least a threshold fraction of its neighbors are in state , then an agent takes state , else(ii)if at least a threshold fraction of its neighbors are in state , then an agent takes state , else(iii)an agent remains in its old state.(b)With probability , an agent conforms to its neighbors, analogously as in Model I.

Time is updated .

##### 2.3. The Mean-Field Approach

In this work, analogously as in [26], we analyze the model on a complete graph, which means that, for each agent, all other agents in the system are neighbors. On the one hand, one can argue that in the case of the symmetric homogeneous threshold model such a structure will give trivial, easily predictable results. On the other hand, only for the complete graph the mean-field approach is exact. Moreover, most of the results for -voter model with nonconformity were obtained on the complete graph, so this structure is adequate for comparison between models. Finally, we hope that such an approach will be helpful in understanding the nature of the phase transitions within the -voter model and the majority-vote model, what will be discussed in the last section of this paper.

As an aggregated quantity, which fully describes the system in case of the complete graph, we choose an average concentration of agents with positive opinions:where denotes the number of agents in the state at time . Alternatively, we could choose an average opinion (magnetization), which is a natural order parameter [44]:However, for the analytical treatment is more convenient, because within the mean-field approach it gives the probability that a randomly chosen agent, at any site, is positive.

As for all other binary opinion models with a single-flip dynamics, in an elementary time step the number of up spins can increase by one, decrease by one, or remain the same. It means that the concentration of up spins can change only by . We follow the notation from [21, 26]: where , as usual. Of course with the complementary probability the state of the system will not change.

Using above probabilities, we obtain a recursive formula for the concentration of up spins:which for gives the rate equation [45, 46]:

Probabilities depend of course on model’s details and can be derived looking at all possible changes that may occur in a single update. Within the homogeneous symmetrical threshold model with independence (Model I), the following changes are possible: where and denote states of a target agent, and is the floor function of , which follows from the restriction that the number of agents can take only integer value, whereas is a real number.

Within the homogeneous symmetrical threshold model with anticonformity (Model A):In all other situations, the state of the system will not change. Therefore, where are probabilities related to conformity. They are the same for both versions of the model:Within conformity, an agent accepts the opinion of her/his neighbors; i.e., it takes the same state as the majority above the threshold . Thus the probability that an agent changes her/his opinion from to is equal to the probability that a randomly selected agent is at state multiplied by the probability that at least of her/his neighbors are in the state . Analogously, the probability that an agent changes her/his opinion from to is equal to the probability that a randomly selected agent is at state multiplied by the probability that at least of her/his neighbors are in the state .

On the other hand, , are related to nonconformity; i.e., they depend on the model’s version. In the model with independence an agent changes her/his state to the opposite one with probability . Thus is equal to the probability that a randomly selected agent is at state multiplied by . Analogously, is equal to the probability that a randomly selected agent is at state multiplied by and thus Whereas for anticonformity, Within anticonformity, an agent takes the opposite state to her/his neighbors, which are in the majority above the threshold . Thus the probability that an agent changes her/his opinion from to is equal to the probability that a randomly selected agent is at state multiplied by the probability that at least of her/his neighbors are in the state . Analogously, the probability that an agent changes her/his opinion from to is equal to the probability that a randomly selected agent is at state multiplied by the probability that at least of her/his neighbors are in the state .

For the large systems we use the following approximation, analogously as in [21]: where are positive, finite constants. In such a case equations for and take much simpler forms. For Model I,For Model A,Summations in above formulas can be calculated using the cumulative distribution function of the binomial distribution and thus we obtain for Model Iand for Model AWe use the following notation:where is a binomially distributed random variable with number of trials and success probability in each trial equal to , and is a binomially distributed random variable with number of trials and success probability in each trial : For the system size large enough, we can approximate binomial distribution by the normal one, which should simplify calculations. Unfortunately, even within such an approximation, we are not able to derive an analytical formula for , described by the rate equation (5). However, we can solve the equation numerically or, in a case of the finite system, calculate by iterating (4).

However, usually, we are more interested in the stationary state than in the time evolution. Especially, the aim of this work is to understand the nature of phase transitions induced by the noise and thus we are interested in the dependence between the stationary value of the concentration of up spins and the level of noise (probability of nonconformity) .

From (4) we see that in the stationary state, the probability and should be equal. Thus to calculate stationary values of concentration we should simply solve the equationSolving analytically (19), i.e., finding, as a function of is impossible, but we can easily derive the opposite relations satisfying (19), analogously as it was done for the -voter model [26]. If we use formulas (15) and (16), then we obtain for Model Iwhereas for Model AWe have used above formulas to plot the dependency between and probability for several values of the threshold on Figures 1 and 2. Although the relation is unknown, only the relation was calculated; we can plot diagrams by rotating the figure, analogously as in [26].

#### 3. Results

We have investigated the model within the mean-field approach, described in the previous section, as well as via the Monte Carlo simulations. We have conducted simulations for several system sizes varying from to and for we have obtained satisfying agreement with formulas obtained within MFA for the large system, what can be seen in Figures 1 and 2. Results were averaged only over samples, but this was enough for this size of the system to get the good statistics—error bars in Figures 1 and 2 are of the same size or even smaller than the symbols. Solid lines indicate stable (attracting) steady values of concentration, whereas dashed lines indicate unstable (repelling) steady states.

It is seen that, generally, the dependence between and is quite trivial (linear) for both variants of the model. For , which corresponds to the simple majority, results for Model I and Model A are almost identical, only rescaled. For both modelswhere for Model I, whereas for model A. It means that for there are two stable solutions: for any initial value the system eventually reaches , whereas for any initial value the system eventually reaches . For there is only one stable solution with up-down symmetry, i.e., . This behavior could be identified as the continuous phase transition, although the dependence is trivial.

*For any threshold *, the situation is slightly more complicated. We can still identify a value below which decays or increases monotonically with , depending on the initial condition . However, this time it is not reached from arbitrary value of . In order to reach one of the ordered state, denoted by the solid lines in Figures 1 and 2, the initial concentration of up spins or alternatively initial concentration of down spins has to be larger than a threshold ; i.e., or . Having this in mind, for we can actually rewrite (22):and again for Model I, whereas for model A for arbitrary value of . The threshold can be easily derived from the condition:We see that for , we obtain , as expected, and (23) becomes almost identical as (22).

However, there is one crucial difference between the case and . For , independently on the version of the model, disordered steady state is unstable for and stable for . For each version of the model behaves differently. To better illustrate differences between Model I and Model A, we present trajectories for on Figure 3, as well as the flow diagrams on Figure 4 for a fixed value of .

**(a)**

**(b)**

**(a)**

**(b)**

*Let us first discuss results for Model I,* shown in Figure 1 and on the left panels of Figures 3 and 4. In this case, is stable for any value of and for two additional symmetrical unstable steady states appear: and . This means that from the initial state the system will be attracted to the disordered state for any value of . In result there is jump of size at , where from (24)Moreover, one could also identify hysteresis, because for the stationary concentration of up spins depends on the initial state . So, looking naively at Figure 1 one could argue that for the model with independence for any value of there is a discontinuous phase transition: the jump of the order parameter increases, whereas hysteresis decreases with growing . Of course, again this result is trivial after a short moment of reflection, but we will discuss it later.

*For the model A*, as shown in Figure 2 and on the right panels of Figures 3 and 4, there is no jump and even for the system remains ordered; i.e., if we start from the initial condition or if we start from the initial condition . For system does not evolve and .

Probably, some of the readers wonder why we use such a formal approach, which eventually led to (20)-(21), as well as the time-consuming Monte Carlo simulations, if nearly all results could be easily deduced directly from the model’s assumptions, without any calculations. Indeed, in the case of this model it would be possible and will be discussed in the next subsection. However, usually the standard approach that we have presented above is the only possible way to obtain results. That was the case for -voter model with independence for which the Monte Carlo simulations, as well as analytical results have been obtained on the complete graph [26] and on various random graphs [36]. Still, heuristically it was not explained why discontinuous phase transition appears in such a model.

Before we proceed to the heuristic explanation of the observed phenomena,* let us summarize the most important results.* Within the model with independence: a discontinuous phase transition is observed, from the completely ordered initial state the jump of the order takes place at and the jump size is equal , and consequently for the transition is continuous. Within the model with anticonformity: for the stationary concentration of the positive opinion decays with and for it has a constant value or depending on the initial value of and consequently for a continuous phase transition is observed at .

##### 3.1. Heuristic Explanation of the Obtained Results

We start with the* heuristic explanation for independence*. If the fraction of up spins is smaller than threshold , i.e., at any time and simultaneously the fraction of down spins is smaller than threshold , i.e., , then, due to the model’s assumptions, conformity cannot take place and changes are caused exclusively by independence. Therefore, for there are only random flips in both directions and as a result we should obtain stationary value , independently on . Indeed we obtain such a result from simulations, as well as from analytical reasoning. For only up spins can cause conformity and therefore with probability spins flip randomly and with probability they flip up. This means that for the system will always reach ordered state with all spins up and thus . On the other hand, for only random flips will occur, which leads to . Expecting the linear dependence between and we can immediately draw the line, which crosses these two points and ; i.e., . Analogous reasoning can be provided for . In such a case the line should cross points and , which gives . Alternative way to deduce the above relation is to realize that in general a fraction of spins flips with probability competing with the complete order , if only and thusWe obtain the value of from the condition: because, as written above, for the system will reach disorder, i.e., . It should become clearer looking at Figure 4. Exactly the same results were obtained from the Monte Carlo simulations, as well as the mean-field approach. However, this simple explanation helps to realize that the jump of at is caused just because conformity cannot take place if and thus in this range the system is attracted to the disordered state with up-down symmetry.

*For the model with anticonformity* heuristic explanation is equally simple. If the fraction of up spins is smaller than threshold , i.e., and simultaneously the fraction of down spins is smaller than threshold , i.e., , then there is not enough social pressure for conformity, as well as for anticonformity. Therefore, for there are no changes in the system and . For only up spins drive changes in the system: with probability a spin flips from to or with probability a spin flips from to , and thuswhere can be obtained analogously like for independence: We can perform analogous reasoning for ; the only difference is that then only down spins drive changes in the system, so again the results will be symmetrical with respect to the line .

#### 4. Discussion

The modification of the Watts threshold model that we have proposed here may be treated as a destruction of the model from the social point of view. However, our aim was not to propose a model describing properly some social phenomena but to understand the nature of the phase transitions observed within models of binary opinions with a single-flip dynamics and up-down symmetry.

It should be noticed that the Model A, proposed here, can be treated as the generalization of the original majority-vote model, which corresponds to . Within such a model only continuous phase transitions are observed, even in the presence of an additional noise [29, 39]. In [29] the model was studied via the mean-field approximation on the square lattice and it was suggested that maybe for larger number of neighbors discontinuous phase transitions would be observed. However, the same model was recently examined on various graphs and it was shown that the presence of independence does not change the type of the phase transition, which remains continuous even for highly connected networks [39]. Indeed we show here that for transition is always continuous even in the limiting (the most connected) case of a complete graph. However, we have shown that discontinuous phase transition may appear within the model with independence for . This result agrees with those obtained for the threshold -voter model, for which also discontinuous phase transition is possible only if [21].

However, results obtained here help not only to understand the difference between the -voter model and the majority-vote model, but additionally they help to understand what the difference between anticonformity and independence is. For both variants of the model proposed here, conformity cannot take place for and thus ordering cannot occur. However, there is big difference between Model I and Model A, which cannot be visible for : independence can take place for , whereas anticonformity cannot. In result for Model I there is an attracting point for any value of , whereas freezing of the system, i.e., within Model A. This causes jump between partially ordered state to the disordered state within Model I whereas freezing (and in consequence a lack of jump) within Model A.

Threshold plays one more role, namely, introduces a kind inertia on the microscopic level. For a spin always takes a position of majority, independently of its state, whereas for it may happen that there is no majority above the threshold and in such a case neighborhood does not influence a target spin. This point of view is particularly interesting if we recall recent results obtained in [47, 48]. It has been shown that inertia introduced on the microscopic level into the majority-vote model, which causes that spin-flip probability of a given spin depends not only on the states of its neighbors, but also on its own state, can change the type of the phase transition from continuous to discontinuous. Interestingly, in [48] the dependence between the order parameter, i.e., average magnetization , and the control parameter, i.e., the probability of anticonformity (denoted by ), was also nearly linear for . Moreover was stable for , similarly as in our model.

Another possibility of obtaining a discontinuous phase transition has been recently suggested within the generalized threshold -voter model, in which two types of nonconformity were introduced simultaneously [28]. It has been shown that discontinuous phase transition could be obtained even without independence if only the threshold for anticonformity would be smaller than for conformity. In the face of the reasoning carried out here, this result also becomes obvious, because the “forbidden” range for conformity is larger than for anticonformity and the jump will increase with the difference between these two thresholds.

We are aware of the fact that many people may wonder why to care about the type of the phase transition. Is there any reason, other than academic, to distinguish between continuous and discontinuous phase transitions? In the face of the social observations, and more recently also laboratory experiments, it seems that discontinuous phase transitions are particularly important, mainly because of the notion of the social hysteresis and the critical mass [49–53]. Both phenomena are strictly related to discontinuous phase transitions. There is no hysteresis within continuous phase transitions, which means that the state of the system is fully determined by the external conditions and does not depend on the history of the system. On the other hand, hysteresis, which is observed within discontinuous phase transitions, means that under the same external conditions the system can be in different states depending on its previous states (history). The social hysteresis was observed in animal [49, 51, 52] as well as in human societies [50, 54, 55]. The second phenomenon, the so-called critical mass, which was recently observed experimentally in social convention [53], is also strongly related to the discontinuous phase transitions. Within the continuous phase transitions, there is no phase coexistence and the transition appears strictly at a given critical point due to the fluctuations and the infinite range of the correlations between them. Because of these correlations, the transition takes place immediately and simultaneously in the entire system. On the other hand, a seed (a “critical mass”) initiating the transition is needed to change the phase under discontinuous phase transition.

We realize that more interesting results could be obtained for different homogeneous and heterogeneous networks. Moreover, both noises could be introduced simultaneously as in [28, 29, 39]. However, the aim of this work was different—we wanted to understand the reason for the discontinuous phase transition in the simplest possible settings. Still, we believe that further studies of the model proposed here would be desirable task for the future, since it is a simple generalization of the majority-vote rule.

#### Data Availability

No data were used to support this study (no empirical data were used; only analytical calculations and Monte Carlo simulations were conducted).

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work was supported by funds from the National Science Center (NCN, Poland) through grant no. 2016/21/B/HS6/01256.

#### References

- C. Castellano, S. Fortunato, and V. Loreto, “Statistical physics of social dynamics,”
*Reviews of Modern Physics*, vol. 81, no. 2, pp. 591–646, 2009. View at: Publisher Site | Google Scholar - S. Galam,
*Sociophysics: A Physicist’s Modeling of Psycho-Political Phenomena*, Springer-Verlag, New York, NY, USA, 2012. View at: Publisher Site | MathSciNet - P. Sen and B. K. Chakrabarti,
*Sociophysics: An Introduction*, Oxford University Press, 2013. - A. Srbu, V. Loreto, V. D. P. Servedio, and F. Tria, “Opinion dynamics: models, extensions and external effects,” in
*Participatory Sensing, Opinions and Collective Awareness*, pp. 363–401, Springer, Cham, Switzerland, 2017. View at: Google Scholar - G. Deffuant, D. Neau, F. Amblard, and G. Weisbuch, “Mixing beliefs among interacting agents,”
*Advances in Complex Systems (ACS)*, vol. 3, pp. 87–98, 2000. View at: Publisher Site | Google Scholar - R. Hegselmann and U. Krause, “Opinion dynamics and bounded confidence: models, analysis and simulation,”
*Journal of Artificial Societies and Social Simulation*, vol. 5, no. 3, 2002. View at: Google Scholar - T. Krueger, J. Szwabiński, and T. Weron, “Conformity, anticonformity and polarization of opinions: insights from a mathematical model of opinion dynamics,”
*Entropy*, vol. 19, no. 7, p. 371, 2017. View at: Publisher Site | Google Scholar - K. Sznajd-Weron, J. Szwabiński, R. Weron, and T. Weron, “Rewiring the network. What helps an innovation to diffuse?”
*Journal of Statistical Mechanics: Theory and Experiment*, vol. 2014, no. 3, p. P03007, 2014. View at: Publisher Site | Google Scholar - K. Byrka, A. Jȩdrzejewski, K. Sznajd-Weron, and R. Weron, “Difficulty is critical: The importance of social factors in modeling diffusion of green products and practices,”
*Renewable & Sustainable Energy Reviews*, vol. 62, pp. 723–735, 2016. View at: Publisher Site | Google Scholar - T. Weron, A. Kowalska-Pyzalska, and R. Weron, “The role of educational trainings in the diffusion of smart metering platforms: An agent-based modeling approach,”
*Physica A: Statistical Mechanics and its Applications*, vol. 505, pp. 591–600, 2018. View at: Publisher Site | Google Scholar - S. Galam, “From 2000 Bush–Gore to 2006 Italian elections: voting at fifty-fifty and the contrarian effect,”
*Quality & Quantity*, vol. 41, no. 4, pp. 579–589, 2007. View at: Publisher Site | Google Scholar - S. Galam, “The Trump phenomenon: an explanation from sociophysics,”
*International Journal of Modern Physics B*, vol. 31, no. 10, 17 pages, 2017. View at: Publisher Site | Google Scholar | MathSciNet - S. Galam, “Unavowed abstention can overturn poll predictions,”
*Frontiers in Physics*, vol. 6, 24 pages, 2018, https://www.frontiersin.org/article/10.3389/fphy.2018.00024. View at: Publisher Site | Google Scholar - T. M. Liggett,
*Interacting Particle Systems*, Springer, 1985. - T. Tome, M. J. Oliveira, and M. A. Santos, “Non-equilibrium Ising model with competing Glauber dynamics,”
*Journal of Physics A: Mathematical and General*, vol. 24, no. 15, pp. 3677–3686, 1991. View at: Publisher Site | Google Scholar - M. J. de Oliveira, “Isotropic majority-vote model on a square lattice,”
*Journal of Statistical Physics*, vol. 66, no. 1-2, pp. 273–281, 1992. View at: Publisher Site | Google Scholar - S. Galam, “Social paradoxes of majority rule voting and renormalization group,”
*Journal of Statistical Physics*, vol. 61, no. 3-4, pp. 943–951, 1990. View at: Publisher Site | Google Scholar - K. Sznajd-Weron and J. Sznajd, “Opinion evolution in closed community,”
*International Journal of Modern Physics C*, vol. 11, no. 6, pp. 1157–1165, 2000. View at: Publisher Site | Google Scholar - D. J. Watts, “A simple model of global cascades on random networks,”
*Proceedings of the National Acadamy of Sciences of the United States of America*, vol. 99, no. 9, pp. 5766–5771, 2002. View at: Publisher Site | Google Scholar | MathSciNet - C. Castellano, M. A. Muñoz, and R. Pastor-Satorras, “Pastor-Satorras, Nonlinear q-voter model,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 80, no. 4, 2009. View at: Publisher Site | Google Scholar - P. Nyczka and K. Sznajd-Weron, “Anticonformity or independence?—insights from statistical physics,”
*Journal of Statistical Physics*, vol. 151, no. 1-2, pp. 174–202, 2013. View at: Publisher Site | Google Scholar | MathSciNet - A. R. Vieira and C. Anteneodo, “Threshold q-voter model,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 97, no. 5, 2018. View at: Publisher Site | Google Scholar - J. P. Gleeson, “Binary-state dynamics on complex networks: pair approximation and beyond,”
*Physical Review X*, vol. 3, no. 2, 2013. View at: Publisher Site | Google Scholar - M. Granovetter, “Threshold models of collective behavior,”
*American Journal of Sociology*, vol. 83, no. 6, pp. 1420–1443, 1978. View at: Publisher Site | Google Scholar - A. Jędrzejewski and K. Sznajd-Weron, “Impact of memory on opinion dynamics,”
*Physica A: Statistical Mechanics and its Applications*, vol. 505, pp. 306–315, 2018. View at: Publisher Site | Google Scholar - P. Nyczka, K. Sznajd-Weron, and J. Cisło, “Phase transitions in the q-voter model with two types of stochastic driving,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 86, no. 1, 2012. View at: Publisher Site | Google Scholar - D. J. Watts and P. S. Dodds, “Inf luentials, networks, and public opinion formation,”
*Journal of Consumer Research*, vol. 34, no. 4, pp. 441–458, 2007. View at: Publisher Site | Google Scholar - P. Nyczka, K. Byrka, P. R. Nail, K. Sznajd-Weron, and N. Fytas, “Conformity in numbers—Does criticality in social responses exist?”
*PLoS ONE*, vol. 13, no. 12, p. e0209620, 2018. View at: Publisher Site | Google Scholar - A. R. Vieira and N. Crokidakis, “Phase transitions in the majority-vote model with two types of noises,”
*Physica A: Statistical Mechanics and its Applications*, vol. 450, pp. 30–36, 2016. View at: Publisher Site | Google Scholar - S. Oh and M. A. Porter, “Complex contagions with timers,”
*Chaos: An Interdisciplinary Journal of Nonlinear Science*, vol. 28, no. 3, p. 033101, 2018. View at: Publisher Site | Google Scholar - A. Jedrzejewski and K. Sznajd-Weron, “Statistical physics of opinion formation: is it a spoof?” 2019, https://arxiv.org/abs/1903.04786. View at: Google Scholar
- S. Galam, “Sociophysics: a review of galam models,”
*International Journal of Modern Physics C*, vol. 19, no. 3, pp. 409–440, 2008. View at: Publisher Site | Google Scholar - S. Galam, “The drastic outcomes from voting alliances in three-party democratic voting (1990 → 2013),”
*Journal of Statistical Physics*, vol. 151, no. 1-2, pp. 46–68, 2013. View at: Publisher Site | Google Scholar - P. R. Nail, G. MacDonald, and D. A. Levy, “Proposal of a four-dimensional model of social response.,”
*Psychological Bulletin*, vol. 126, no. 3, pp. 454–470, 2000. View at: Publisher Site | Google Scholar - P. R. Nail, S. I. Di Domenico, and G. MacDonald, “Proposal of a double diamond model of social response.,”
*Review of General Psychology*, vol. 17, no. 1, pp. 1–19, 2013. View at: Publisher Site | Google Scholar - A. Jędrzejewski, “Pair approximation for the q-voter model with independence on complex networks,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 95, no. 1, 2017. View at: Publisher Site | Google Scholar - A. Chmiel and K. Sznajd-Weron, “Phase transitions in the q-voter model with noise on a duplex clique,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 92, no. 5, 2015. View at: Publisher Site | Google Scholar - A. F. Peralta, A. Carro, M. San Miguel, and R. Toral, “Analytical and numerical study of the non-linear noisy voter model on complex networks,”
*Chaos: An Interdisciplinary Journal of Nonlinear Science*, vol. 28, no. 7, p. 075516, 2018. View at: Publisher Site | Google Scholar - J. Encinas, H. Chen, M. M. de Oliveira, and C. E. Fiore, “Majority vote model with ancillary noise in complex networks,”
*Physica A: Statistical Mechanics and its Applications*, vol. 516, pp. 563–570, 2019. View at: Publisher Site | Google Scholar - M. Mobilia, “Nonlinear q-voter model with inflexible zealots,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 92, no. 1, 2015. View at: Publisher Site | Google Scholar - A. Mellor, M. Mobilia, and R. K. Zia, “Characterization of the nonequilibrium steady state of a heterogeneous nonlinear q-voter model with zealotry,”
*EPL (Europhysics Letters)*, vol. 113, no. 4, p. 48001, 2016. View at: Publisher Site | Google Scholar - A. Mellor, M. Mobilia, and R. K. Zia, “Heterogeneous out-of-equilibrium nonlinear q-voter model with zealotry,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 95, no. 1, 2017. View at: Publisher Site | Google Scholar - P. Nyczka, J. Cisło, and K. Sznajd-Weron, “Opinion dynamics as a movement in a bistable potential,”
*Physica A: Statistical Mechanics and its Applications*, vol. 391, no. 1-2, pp. 317–327, 2012. View at: Publisher Site | Google Scholar - L. D. Landau, “On the theory of phase transitions,”
*Journal of Experimental and Theoretical Physics*, vol. 7, pp. 19–32, 1937. View at: Google Scholar - H. Hinrichsen, “Non-equilibrium critical phenomena and phase transitions into absorbing states,”
*Advances in Physics*, vol. 49, no. 7, pp. 815–958, 2000. View at: Publisher Site | Google Scholar - P. L. Krapivsky, S. Redner, and E. Ben-Naim,
*A Kinetic View of Statistical Physics*, Cambridge University Press, 2010. - H. Chen, C. Shen, H. Zhang, G. Li, Z. Hou, and J. Kurths, “First-order phase transition in a majority-vote model with inertia,”
*Physical Review E: Statistical, Nonlinear, and Soft Matter Physics*, vol. 95, no. 4, 2017. View at: Publisher Site | Google Scholar - J. M. Encinas, P. E. Harunari, M. M. de Oliveira, and C. E. Fiore, “Fundamental ingredients for discontinuous phase transitions in the inertial majority vote model,”
*Scientific Reports*, vol. 8, no. 1, 2018. View at: Publisher Site | Google Scholar - M. Beekman, D. J. Sumpter, and F. L. Ratnieks, “Phase transition between disordered and ordered foraging in Pharaoh's ants,”
*Proceedings of the National Acadamy of Sciences of the United States of America*, vol. 98, no. 17, pp. 9703–9706, 2001. View at: Publisher Site | Google Scholar - M. Scheffer, F. Westley, and W. Brock, “Slow Response of Societies to New Problems: Causes and Costs,”
*Ecosystems*, vol. 6, no. 5, pp. 493–502, 2003. View at: Publisher Site | Google Scholar - J. N. Pruitt, A. Berdahl, C. Riehl et al., “Social tipping points in animal societies,”
*Proceedings of the Royal Society B Biological Science*, vol. 285, no. 1887, p. 20181282, 2018. View at: Publisher Site | Google Scholar - G. N. Doering, I. Scharf, H. V. Moeller, and J. N. Pruitt, “Social tipping points in animal societies in response to heat stress,”
*Nature Ecology & Evolution*, vol. 2, no. 8, pp. 1298–1305, 2018. View at: Publisher Site | Google Scholar - D. Centola, J. Becker, D. Brackbill, and A. Baronchelli, “Experimental evidence for tipping points in social convention,”
*Science*, vol. 360, no. 6393, pp. 1116–1119, 2018. View at: Publisher Site | Google Scholar - J. Elster, “A note on hysteresis in the social sciences,”
*Synthese*, vol. 33, no. 1, pp. 371–391, 1976. View at: Publisher Site | Google Scholar - A. Clark, “Unemployment as a social norm: psychological evidence from panel data,”
*Journal of Labor Economics*, vol. 21, no. 2, pp. 323–351, 2003. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2019 Bartłomiej Nowak and Katarzyna Sznajd-Weron. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.