1. Introduction
In recent decades, statistical physics has expanded its scope to venture into the field of sociology, giving rise to a discipline called sociophysics [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10]. A commonly studied phenomenon is the dynamics of opinion formation, by means of simple mathematical models. In these models, individuals are called agents, and each of them is characterized by the value of a variable that represents its opinion on a particular topic—such as the intention to vote for a candidate on a ballot—which, for simplicity, can take one of two possible values (+1 or −1). The opinion of each agent can change after interacting with other agents following simple rules. One of the most implemented interaction rules is that introduced in a model by Galam [
11] and extensively studied later on [
6,
12,
13,
14], to which we refer as the Galam Majority Model (GMM), in which all agents in a group chosen at random adopt the opinion of the majority in that group. This local dynamic drives a steady increase in the initial global majority opinion (provided the system’s symmetry is not broken by ties for even-size groups), which eventually ends with a consensus, i.e., an absorbing state where all agents share the same opinion. Multiple extensions of the GMM have been studied in the literature, including the possibility of contrarian behavior, that is, all members of a chosen group taking the minority opinion [
7]. This work studied the effects of introducing a fixed fraction
a of contrarian agents to the original GMM, where it was found that, instead of a frozen consensus, as in the model with no contrarians, the system reaches an ordered stationary state for
and a disordered stationary state for
. The transition value
separates an ordered phase, where a large majority of agents hold the same opinion, from a disordered phase, in which both opinions are equally represented in the population.
Many other opinion formation models with contrarians were also studied in [
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34]. In particular, the effects of contrarian behavior were also investigated in the voter model (VM) for opinion formation [
24], where agents interact in pairs and one adopts the opinion of the other, with a probability of
(imitation), or the opposite opinion, with a probability of
(contrarian). It was shown that the model displays a transition from order to disorder when the probability of having a contrarian behavior overcomes the threshold
in a system of
N agents. The contrarian voter model [
24] was recently studied under the presence of a mass media propaganda that influences agents’ decisions [
34]. The propaganda was implemented in the form of an external oscillating field that tends to align agents’ opinions in the direction of the field. A stochastic resonance (SR) phenomena within an oscillatory regime was found, that is, there is an optimal level of noise for which the population effectively responds to the modulation induced by the external field [
35,
36].
In order to expand our knowledge on the combined effects of contrarians and mass media propaganda on opinion models, we study in this article the GMM with contrarian behavior under the presence of an external oscillating field. Each agent in the population can either follow a majority rule that increases similarity with its neighbors or behave as a contrarian by adopting the opposite opinion, with respective probabilities
and
. The majority probability
varies periodically in time according to an external field or propaganda, based on a mathematical form introduced in [
21,
22] for the Sznajd model and implemented in [
34,
37] for the VM, so that agents tend to follow the majority when it is aligned with the field. Although we understand that this model is relatively far from being adequate to describe real social systems where mass media does not necessarily change sign periodically, it provides analytical insight on features that could be actually relevant in the real word. We also notice that, unlike models where a predetermined group of agents always behaves as a contrarian (adopt the minority opinion) and the rest of agents always behave as a follower (adopt the majority opinion), like in [
7], in our model each agent can behave as both, as a contrarian or as a follower at each iteration step of the dynamics, with a frequency of choosing the contrarian behavior given by
. However, we note that these two models lead to identical mean-field (MF) equations for the system’s evolution in an all-to-all interaction setup. By exploring the dynamics of the GMM under the influence of an oscillating external field and the presence of contrarians, we aim to gain deeper insights into the manifestation of the SR phenomenon in opinion dynamics models. We show that this model exhibits unimodal and bimodal oscillatory regimes, as well as an SR that is observed close to the transition between the two regimes.
It is worth mentioning that the GMM studied in this article belongs to a class of “non-linear” models, while the VM with contrarians described above and studied in [
34] belongs to a completely different class. That is, in the VM an agent switches opinion (
,
) with a probability equal to the fraction of neighbors that hold the opposite opinion
(linear in
), while in the GMM the switch happens with a probability that is proportional to a combination of
and
, i.e., non-linear in
. Then, the equation that governs the time evolution of the fraction of agents in a given opinion can be written as a Ginzburg–Landau equation with an associated potential that has two symmetric wells in the GMM, as in the Ising model for ferromagnetism, while in the VM the potential is zero. A main consequence of this difference in the potential is in the type of domain coarsening observed in spatial systems: while in the GMM the coarsening is driven by curvature, in the VM it is without surface tension [
38] (driven by noise). Another consequence is that, in the version of these models with contrarians, the order–disorder transition in the thermodynamic limit (
) takes place at a finite fraction of contrarian agents
in the GMM, while in the VM the transition happens at a vanishing contrarian probability (
).
We also need to mention that the SR effect has also been observed in other opinion models. For instance, in [
21,
22] the authors found SR in a variation of the Sznajd model with stochastic driving and a periodic signal. The work in [
14] analyzed a majority rule dynamics under the action of noise and an external modulation, and found an SR that depends on the randomness of the small-world network. There are also other works [
39,
40,
41,
42,
43] that explored the combined effects of stochastic driving and an external signal on a majority rule dynamics. However, none of these works have incorporated a contrarian behavior in the dynamics.
The rest of this article is organized as follows. We introduce the model in
Section 2. In
Section 3, we present numerical simulation results for the evolution of the system and the behavior of different magnitudes that characterize the SR phenomena. In
Section 4, we develop an MF approach that gives an insight into the system’s evolution and the relation between the SR and the transition between different regimes. Finally, in
Section 5 we summarize our findings and discuss the results.
2. The Model
We consider a population of
N interacting agents where a given agent
i (
) can hold one of two possible opinion states,
. We denote by
and
the fraction of nodes with respective states
and
at time
t, such that
for all
. In a time step
of the dynamics, we follow the basic GMM using groups of size three to update individual opinions. However, here, for our purpose of investigating the effects of propaganda at the level of a single individual, we implement the rule in a different setting, where only one individual can update its opinion state in an iteration step. Instead of selecting three agents randomly to update all of them at once, we pick one agent
i with state
and a group of three other different agents
(
), all randomly chosen. In the
limit, their respective states are
with probability
. A majority of + choices is thus obtained for the configurations
,
,
and
, yielding an overall probability
Similarly, a majority of − occurs for
,
,
and
, with the overall probability
Then, agent
i updates its state in two steps. (i) First, the update follows the basic GMM, where agent
i simply adopts the majority state of the group of the three agents
and
l. We thus have
with probability
, or
with probability
. (ii) Second, agent
i can either preserve this majority state (
) with probability
, or change to the opposite (minority) state (
) with the complementary probability
, where
is defined below. The implication of this second step is that each agent can behave as a “contrarian” by adopting the state opposed to the majority (minority state) with probability
, or as a “majority follower” with probability
. Thus, there is no fixed fraction of contrarian agents, as in [
7].
At this point, we introduce the effect of an external field
H on agent
i in state
within a Boltzmann scheme by assuming that the probability
to preserve the majority state is larger when
is aligned with
H [i.e.,
],
where
is a parameter that plays the role of a
social temperature analogous to the contrarian feature of the GMM. The related probability to oppose the field is
. We assume that
H is an oscillating periodic field
with amplitude
(
), frequency
and period
, which represents an external propaganda. Thus, according to Equation (
3), agents are more likely to keep the opinion that is aligned with the propaganda. In addition to the external field, we introduce an individual “inertia” parameter
I, which provides an agent with a weight to preserve its current state against a field favoring the opposite state. It is a self-interaction
which modifies Equation (
3) as
which can be rewritten as
where
and
T have been rescaled as
and
using
.
At this stage, we combine the GMM with the inertia and field effects by taking
for the probability of agent
i to keep the majority state
, and
for the probability to adopt the opposite (minority) state
, which can be interpreted as a noise. Finally, combining Equations (
1), (
2) and (
6), the probability
for a randomly selected agent
i to adopt the state + in a single time step
is given by
where the first term comes from following a local majority + among the three selected agents, which happens with probability
, while the second term corresponds to opposing the state − in case of a majority of − among the three selected agents, which happens with probability
. Analogously, the state − is selected with probability
.
As noted above, in a single iteration only one agent (the “focal agent”
i) updates its state, unlike in the original Galam’s model where all three agents in the chosen group update their states at once [
11]. However, as these two settings use the same majority rule and they differ only in the number of agents updated in an iteration, they turn out to be equivalent at the MF level, i.e., in an all-to-all interaction setup. That is, they both have the same rate equations for the evolution of macroscopic quantities (see
Section 4). Therefore, the model introduced in this section is the same as the original Galam’s model, with the novelty of the addition of an external oscillating field that is coupled to agents’ opinions, which extends the parameter space of Galam’s model and leads to new phenomena, as we shall see.
Equation (
6) shows that individuals are more prone to adopt the opinion of the majority when it is aligned with the propaganda. In addition,
and
approach the value one as
, which makes this case equivalent to the original GMM, with neither contrarians nor an external field. In the opposite limit
,
and
approach the value
, which corresponds to the pure noise case where agents take one of the two opinions at random, independent of the field.
In the Glossary, we display a table with a list of the symbols of variables and their names that we use throughout the text.
4. Mean-Field Approach
In this section, we analyze the behavior of the model within an MF approach by deriving a rate equation for the evolution of
m that corresponds to the dynamics introduced in
Section 2. Let us write the fractions of + and − agents in terms of the magnetization
m,
and
. As we described in
Section 2, in a time step
, a random agent
i with state
is chosen with probability
, and then adopts the state + (
flip) with probability
, which corresponds to adopt either the majority state + of a selected + majority, or the minority state + of a selected − majority, where
and
are given by Equations (
1) and (
2), respectively. This flip
leads to an overall change
in
m. Conversely, with probability
, the chosen agent
i has state
, and flips to
(
flip) with probability
, leading to a change
. Assembling these factors, the mean change of
m in a time step can be written as
which becomes, in the
limit, the rate equation
after replacing the expressions for
and
and doing some algebra. Here,
are the probabilities of adopting the state
and
of a majority, respectively, as defined in Equation (
6).
We note that if we consider the original Galam’s model where, in an iteration step, a group of three agents are chosen at random and all adopt the majority state, we can derive a rate equation for
m that has the same form as Equation (
8), with a prefactor of three multiplying the right-hand-side of the equation. Therefore, both models turn out to be equivalent in MF; they have the same stationary states and the same behavior. They only differ by a factor of three in the time scale associated with the relaxation to the steady state.
For the zero field case, (
) is
, and, thus, Equation (
8) is reduced to the simple equation
Equation (
10) has three fixed points corresponding to the possible stationary states of the agent-based model. The fixed point
is stable for
and corresponds to a disordered active state with equal fractions of + and − agents (
), whereas the two fixed points
are stable for
, and they represent asymmetric active states of coexistence of + and − agents, with stationary fractions
and
. The stable fixed points are plotted with a solid line in
Figure 1b, where we observe a good agreement with MC simulation results (solid circles). Equation (
11) shows the existence of a transition from order to disorder as
T overcomes the value
(
), as we already mentioned in
Section 3.1. Notice that the probability of behaving as a contrarian
is identical to the critical proportion of contrarians
obtained in the GMM for groups of size three [
7]. Given that Equation (
10) can be rewritten as a Ginzburg–Landau equation with an associated double-well potential with two minima at
, we expect a bistable regime for
, where, in a single realization,
m jumps between
and
.
For a field that is constant in time (
), the fixed points of Equation (
8) are given by the roots of a cubic polynomial, and
is no longer a root. Only one root is real, and corresponds to the stationary state of the agent-based model. As the analytical expression for the real root is large and not very useful, we integrated Equation (
8) numerically to find the stationary value
, which we plot with a dashed line in
Figure 1b for
and
. We observe a good agreement with MC simulations (symbols). A positive field
breaks the symmetry in favor of the + state, given that
, leading to a positive stationary value
that increases monotonically with
H.
For an oscillating field
, we have that
and
oscillate in time according to
, which in turn leads to oscillations in
. In order to explore, within the MF approach, the behavior of
m in the different regimes described in
Section 3.1, we plot in
Figure 5a the evolution of
m obtained from the numerical integration of Equation (
8) for
,
and various temperatures. For low temperatures, we see that
m oscillates around a positive value (it could also be a negative value for other initial conditions), but when the temperature is increased beyond a threshold value, oscillations turn to be around
. At first sight, this transition that happens in the oscillatory regime of
m, already reported in
Section 3.1 from MC simulations, appears to be quite sharp, where the center of oscillations seems to jump from a large value to zero after a small increment of
T. To better characterize the transition, we plot in
Figure 5b the temporal average of
m from
to
, called
, as a function of
T and for several periods
. The value of
can be seen as an order parameter, which takes a positive or negative value in the bimodal regime and a value close to zero in the unimodal regime. We can see that
decreases continuously with
T for low
(see curve for
), and that the transition becomes more abrupt as
increases (see curves for
). The inset shows a more detailed view of the transition in the value of
.
In
Figure 2a, we compare the evolution of
m obtained from the MF approach (dashed lines) with that from MC simulations, for
,
and various temperatures. We observe a good agreement with single realizations of the dynamics, except for the temperature
, which is close to the transition value
, estimated from
Figure 5b as the point where
becomes zero. This discrepancy is due to the fact that the MF approach cannot reproduce the random jumps of
from the value
in the bimodal regime to
in the unimodal regime. These jumps are induced by finite-size fluctuations, and are more frequent when the control parameter
T is close to the transition point
.
An insight into the behavior of the resonance temperature
with the period
can be obtained from the MF approach, assuming that the response
reaches a maximum value at a temperature similar to the transition point
, that is, we expect,
. This is because, in the bimodal regime
, the magnetization
m oscillates around a positive or a negative value and eventually crosses
around times
,
, etc., via finite-size fluctuations, leading to multiple peaks in the residence time distribution. Then, at
, oscillations start to be centered at
, and, thus, we expect that the
shows a single peak at
. By increasing
T beyond
, we expect that the height of the peak for
is reduced by the presence of a higher noise that induces another maximum of the
at
, as explained in
Section 3.2, leading to a smaller
. Therefore, we expect that
is at its maximum at
.
Figure 4b shows in diamonds the value of
obtained from
Figure 5b for various periods
. We see that
decreases with
, as it happens with
(circles), although discrepancies between
and
increase as
decreases.
5. Summary and Discussion
In this article, we studied the dynamics of the binary-state majority rule model introduced by Galam for opinion formation under the presence of an external propaganda and contrarian behavior. When an agent has to update its opinion, it can either follow the majority opinion among three random neighbors, similarly to the original GMM, or adopt the opposite (contrary) opinion, i.e., the minority opinion. The probability to adopt the majority opinion is coupled to an external field that oscillates periodically in time (propaganda), in a way that agents are more likely to adopt the majority opinion when it is align with the field. This rule tries to reproduce a reinforcing mechanism by which individuals have a tendency to follow the majority opinion when it is in line with mass media propaganda. Additionally, the majority probability depends on a parameter T (temperature), which acts as an external source of noise, in such a way that by increasing T from zero the system goes from following the majority opinion only ( for ) to adopting a random opinion for large temperatures ( for ).
We explored the model in a complete graph (all-to-all interactions) and found different phenomena associated with different regimes as T is varied. For T below a threshold value , the system is in a bimodal regime, where the mean opinion m oscillates in time around a positive or negative value, , and performs jumps between and due to finite-size fluctuations, similarly to what happens in a bistable system. As the temperature is increased beyond , there is a transition to an unimodal regime in which m oscillates around zero, where the amplitude of oscillations decreases with T and eventually vanishes in the limit that corresponds to pure noise. The transition at becomes more abrupt as the period of the field increases. We also studied the response of the system to the external field by means of the distribution of residence times, i.e., the time interval between two consecutive changes of the sign of m. We found that there is an optimal temperature for which the response is maximum, that is, a stochastic resonance phenomenon induced by the external noise controlled by T. Furthermore, we developed an MF approach that lead to a non-linear rate equation for the time evolution of m in the thermodynamic limit, whose numerical solution agrees very well with MC simulations of the model. We used this equation to give a numerical estimate of , and found that the behavior of with the period is qualitatively similar to that of . Although the transition temperature is similar to the resonance temperature only for large , this analysis shows that they are related.
A possible interpretation of these results in a social context is the following. Reacting with a contrarian attitude occasionally (small T/low noise) on a given issue, that is, adopting an opposite position to that of the majority of our acquaintances, leads to a state of collective agreement in a population, which can be reversed completely after some time by means of a collective decision, independently of the external propaganda. This alternating behavior between opposite opinions might be seen as more “socially healthy” than a frozen full consensus in one of the two alternatives, which happens in populations with a total absence of contrarian attitudes (). However, having a contrarian behavior more often could induce a collective state where the mean opinion oscillates in time following the external propaganda, which can be interpreted as a society whose opinions are manipulated optimally by the mass media, in opposition to collective freedom. Finally, in the extreme case of having a very frequent contrarian attitude (), the population falls into a state of opinion bipolarization, where there are two groups of similar size with opposite opinions.
Given that the opinion formation model introduced and studied in this article is simple enough, it allowed for a mathematical treatment that led to numerical and analytical results. However, due to its simplicity, the model is not very adequate to describe the complexity observed in real social systems, and, thus, some extensions are necessary to cope with real life scenarios. For instance, one can possibly collect real data to infer how mass media propaganda evolves over time. We could expect that real propaganda oscillates with a frequency that is not constant like in the model, but varies in time, and is centered at an arbitrary positive or negative number in the interval . Nevertheless, we have seen that the model provides insights on features that could be actually relevant in the real word, like the conditions for which the opinions of a society are optimally manipulated by mass media, as described in the paragraph above, related to the stochastic resonance effect.
The results presented in this article correspond to a fully connected network. Although we expect that the conclusions remain valid qualitatively for other interaction topologies, it might be worthwhile to study the model in complex networks like scale-free or Erdös Renyi networks, which better-represent social interactions. This version of the model would be more realistic, not only for the interaction topology, but also because the size of the majority group would be different for each agent or node in the network (its degree or number of neighbors), allowing for the case of ties if the group has an even number of agents. It might also be interesting to explore how the stochastic resonance effect depends on the topology of the network. Finally, the model can be easily adapted to include some of the possible real time features of mass media propaganda described above. We leave all these ideas for future work.