arxiv: v2 [physics.soc-ph] 25 Apr 2011

Transcription

1 Social consensus through the influence of committed minorities J. Xie, 1 S. Sreenivasan, 1, G. Korniss, W. Zhang, C. Lim, and B. K. Szymanski 1 1 Dept. of Computer Science, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY 1180 Dept. of Physics, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY 1180 Dept. of Mathematics, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY 1180 arxiv:110.91v [physics.soc-ph] 5 Apr 011 We show how the prevailing majority opinion in a population can be rapidly reversed by a small fraction p of randomly distributed committed agents who consistently proselytize the opposing opinion and are immune to influence. Specifically, we show that when the committed fraction grows beyond a critical value p c 10%, there is a dramatic decrease in the time, T c, taken for the entire population to adopt the committed opinion. In particular, for complete graphs we show that when p < p c, T c exp(α(p)n), while for p > p c, T c lnn. We conclude with simulation results for Erdős-Rényi random graphs and scale-free networks which show qualitatively similar behavior. PACS numbers: 87..Ge, Fb I. INTRODUCTION Human behavior is profoundly affected by the influenceability of individuals and the social networks that link them together. Well before the proliferation of online social networking, offline or interpersonal social networks have been acknowledged as a major factor in determining how societies move towards consensus in the adoption of ideologies, traditions and attitudes [1, ]. As a result, the dynamics of social influence has been heavily studied in sociological, physics and computer science literature [ 7]. In the sociological context, work on diffusion of innovations has emphasized how individuals adopt new states in behavior, opinion or consumption through the influence of their neighbors. Commonly used models for this process include the threshold model [8] and the Bass model[9]. Akeyfeatureinboththesemodelsisthatonce an individual adopts the new state, his state remains unchanged at all subsequent times. While appropriate for modeling the diffusion of innovation where investment in a new idea comes at a cost, these models are less suited to studying the dynamics of competing opinions where switching one s state has little overhead. Here we address the latter case. From among the vast repertoire of models in statistical physics and mathematicalsociology,wefocus ononewhich isa-opinionvariant [10] of the Naming Game (NG) [11 15] and that we refer to as the binary agreement model. The evolution of the system in this model takes place through the usual NG dynamics, wherein at each simulation time step, a randomly chosen speaker voices a random opinion from his list to a randomly chosen neighbor, designated the listener. If the listener has the spoken opinion in his list, both speaker and listener retain only that opinion, else the listener adds the spoken opinion to his list (see Table 1). The order of selecting speakers and listeners is known to influence the dynamics, and we stick to choos- Corresponding author: sreens@rpi.edu ing the speaker first, followed by the listener. It serves to point out that an important difference between the binary agreement model and the predominantly used opinion dynamics models [4, 6, 16 18] is that an agent is allowed to possess both opinions simultaneously in the former, and this significantly alters the time required to attain consensus starting from uniform initial conditions. Numerical studies in [10] have shown that for the binary agreement model on a complete graph, starting from an initial condition where each agent randomly adopts one of the two opinions with equal probability, the system reaches consensus in time T c lnn (in contrast, for example, to T c N for the voter model). Here N is the number of nodes in the network, and unit time consists of N speaker-listener interactions. The binary agreement model is well suited to understanding how opinions, perceptions or behaviors of individuals are altered through social interactions specifically in situations where the cost associated with changing one s opinion is low [19], or where changes in state are not deliberate or calculated, but unconscious [0]. Furthermore, by its very definition, the binary agreement model may be applicable to situations where agents while trying to influence others, simultaneously also have a desire to reach global consensus [1]. Another merit of the binary agreement dynamics in modeling social opinion change seems worth mentioning. Two state epidemic-like models of social contagion (examples in []) suffer from the drawback that the rules governing the conversion of a node from a given state to the other are not symmetric for the two states. In contrast, in the binary agreement model, both singular opinion states are treated symmetrically in their susceptibility to change. Here we study the evolution of opinions in the binary agreement model starting from an initial state where all agents adopt a given opinion B, except for a finite fraction p of the total number of agents who are committed agents and have state A. Committed agents, introduced previously in [], are defined as nodes that can influence other nodes to alter their state through the usual prescribed rules, but which themselves are immune to

2 influence. In the presence of committed agents adopting state A, the only absorbing fixed point of the system is the consensus state where all influenceable nodes adopt opinion A - the opinion of the committed agents. The question that we specifically ask is: how does the consensus time vary with the size of the committed fraction? More generally, our work addresses the conditions under which an inflexible set of minority opinion holders can win over the rest of the population. The effect of having un-influencable agents has been considered to some extent in prior studies. Biswas et al. [4] considered for two state opinion dynamics models in one dimension, the case where some individuals are rigid in both segments of the population, and studied the time evolution of the magnetization and the fraction of domain walls in the system. Mobilia et. al. [5] considered the case of the voter model with some fraction of spins representing zealots who never change their state, and studied the magnetization distribution of the system on the complete graph, and in one and two dimensions. Our study differs from these not only due to the particular model of opinion dynamics, but also in its explicit consideration of different network topologies and of finite size networks, specifically deriving how consensus times scale with network size for the particular case of the complete graph. Furthermore, the above mentioned studies do not explicitly consider the initial state that we care about - one where the entire minority set is un-influencable. A notable exception to the latter is the study by Galam and Jacobs [6] in which the authors consideredthecaseof inflexibles inatwostatemodelof opinion dynamics with opinion updates obeying a majority rule. While this study provides several useful insights and is certainly the seminal quantitive attempt at understanding the effect of committed minorities, its analysis is restricted to the mean-field case, and has no explicit consideration of consensus times for finite systems. II. COMPLETE GRAPHS A. Infinite network size limit We start along similar lines as [6] by considering the case where the social network connecting agents is a complete graph with the size of the network N. We designate the densities of uncommitted nodes in states A,B as n A,n B respectively. Consequently, the density of nodes in the mixed state AB is n AB = 1 p n A n B, where p is the fraction of the total number of nodes, that are committed. Neglecting correlations between nodes, and fluctuations, one can write the following rate equations for the evolution of densities: dn A dt dn B dt = n A n B +n AB +n ABn A + pn AB = n A n B +n AB +n ABn B pn B (1) Before interaction After interaction A A A A A B A A AB B B A B B B B B AB AB A A AB A B AB A AB AB B A AB B B AB B AB A - A A - AB A - A B - AB B - B B - B A - A AB - AB A - A AB - AB B - B B - B TABLE I: Shown here are the possible interactions in the binary agreement model. Nodes can possess opinion A, B or AB, and opinion updates occur through repeated selection of speaker-listener pairs. Shown in the left column are the opinions of the speaker (first)and listener (second) before the interaction, and the opinion voiced by the speaker during the interaction is shown above the arrow. The column on right shows the states of the speaker-listener pair after the interaction. The terms in these equations are obtained by considering all interactions which increase (decrease) the density of agents in a particular state and computing the probability of that interaction occurring. Table 1 lists all possible interactions. As an example, the probability of the interaction listed in row eight is equal to the probability that a node in state AB is chosen as speaker and a node in state B is chosen as listener (n AB n A ) times the probability that the speaker voices opinion A ( 1 ). The fixed-point and stability analyses (see Appendix) ofthesemean-fieldequationsshowthatforanyvalueofp, theconsensusstateinthecommittedopinion(n A = 1 p, n B = 0) is a stable fixed point of the mean-field dynamics. However, below p = p c = 5 ) ( ( ), two additional fixed points appear: one of these is an unstable fixed point (saddle point), whereas the second is stable and represents an active steady state where n A,n B and n AB are all non-zero (except in the trivial case where p = 0). Fig. 1(a) shows (asterisks) the steady state density of nodes in state B obtained by numerically integrating the mean-field equations at different values of the committed fraction p and with initial condition n A = 0,n B = 1 p. As p is increased, the stable density ofb nodes n B abruptly jumps from to zero at the critical committed fraction p c. A similar abrupt jump also occurs for the stable density of A nodes from a value very close to zero below p c, to a value of 1, indicating consensus in the A state (not shown). In the study of phase transitions, an order parameter is a suitable quantity changing(either contin-

3 uously or discontinuously) from zero to a non zero-value at the critical point. Following this convention, we use n B - the density of uncommitted nodes in state B - as the order parameter appropriate for our case, characterizing the transition from an active steady state to the absorbing consensus state. In practice, for a complete graph of any finite size, consensus is always reached. However, we can still probe how the system evolves, conditioned on the system not having reached consensus. Fig. 1(a) shows the results of simulating the binary agreement model on a complete graph for different system sizes (solid lines). For p < p c, in each realization of agreement dynamics, neglecting the initial transient, the density ofnodes in state B, n B, fluctuates around a non-zero steady state value, until a large fluctuation causes the system to escape from this active steady state to the consensus state. Fig. 1(a) shows these steadystatevaluesofn B conditionedonsurvival, forseveral values of p. As expected, agreement of simulation results with the mean-field curve improves with increasing system size, since Eq. (1) represents the true evolution of the system in the asymptotic network size limit. Accordingly the critical value of the committed fraction obtained from the mean-field equations is designated as p c ( ), although, for brevity, we refer to it simply as p c throughout this paper. The existence of the transition as p is varied and when the initial condition for densities is (n A = 0,n B = 1 p) can be further understood by observing the motion of the fixed points in phase space. Fig. 1(b) shows how the stablefixedpointandthe unstablefixedpointmoveinphase space as p is varied from 0 to p c. The active steady state moves downward and right while the saddle point moves upwards and left. At the critical value p c the two meet and the only remaining stable fixed point is the consensus fixed point. A similar observation was made in the model studied in [6]. The fact that the value of n B converges to 0.65 and does not smoothly approach zero as the stable fixed point and the saddle point approach each other, explains the origin of the first-order nature of the phase transition. Fig shows the representative trajectories obtained by integrating the mean-field equations for the cases where p = 0.05 (< p c ) and p = 0.1 (> p c ). B. Finite network size: Scaling results for consensus times Even though consensus is always reached for finite N, limits on computation time prohibit the investigation of the consensus time, T c, for values of p below or very close top c. Wethereforeadoptasemi-analyticalapproachprescribed in [7] that allows us to estimate the consensus times for different N for an appreciablerange of p including values below p c. We start with the master equation which describes the evolution of the probability that the network of size N has n (m) uncommited nodes in state A (B). We denote by c, the number of committed nodes, and by l(= N n m c), the number of uncommitted nodes in state AB. dp nm dt 1 N = 1 ( N p nm (ln+ lc+nm+l(l 1) (l+1)(n 1+c) + lm+mc)+p n 1,m (n+1)(m+l 1) + p n+1,m (l+)(l+1) + p n,m (l+1)(m 1) + p n,m 1 (m+1)(n+c+l 1) + p n,m+1 ) (l+)(l+1) + p n,m () The factor of 1/N in the LHS comes from the fact that a transition between states takes place in an interval of time1/n. Thetransitionratesineachtermaretheproduct of two densities, which is responsible for the overall factor of 1/N in the RHS. The probabilities are defined over all allowed states of the system ( i.e. 0 n N c, and 0 m N c n for given n ) and the allowed transitions from any point {n, m} in the interior of this state-space are {n,m} {n,m ± 1},{n,m} {n±1,m},{n,m} {n,m+},{n,m} {n+,m}. We know from the mean field equations that in the asymptotic limit, and below a critical fraction of committed agents, there exists a stable fixed point. For finite stochastic systems, escape from this fixed point is always possible, and therefore it is termed metastable. For a finite system, the probability of having escaped to the metastable fixed point as a function of time is P e (t) = 1 P s (t) where P s (t) is the survival probability. The surviving fraction is constrained to be in the allowed region of n, m quadrant excluding the true fixed point {N c,0}. If the number of committed agents is far lower than p c N we expect that this surviving fraction will occupy configurations around the metastable fixed point, and the occupation probabilities p n,m will be peaked around the metastable fixed point. In systems which exhibit such long lived metastable states in addition to an absorbing fixed point, applying a quasistationary (QS) approximation has been found to be useful in computing quantities of interest [7 9]. This approximation assumes that after a short transient, the occupation probability, conditioned on survival, of allowed states excluding the consensus state, is stationary. Following this approximation, the distribution of occupation probabilities conditioned on survival can be written as, p nm = p nm (t)/p s (t) [7] and using this form in the mas-

4 4 FIG. 1: (a) The steady state density n B of nodes in state B as a function of committed fraction p for complete graphs of different sizes, conditioned on survival of the system. Simulation results are from 100 realizations of the binary agreement dynamics. (b) Movement of the stable fixed point and the saddle point in phase space as a function of committed fraction p (see text). The location of these points in phase space was obtained through a stability analysis (not shown) of the mean field equations Eq. (1). The inset shows the density of nodes in state B at the stable (red) and unstable(blue) fixed points as p is varied; the expressions from which these numerical values are obtained are given in the Appendix. ter equation (Eq. ()), we get: dp s (t) p nm = P ( s(t) p nm (ln+ dt N lc+nm + l(l 1)+lm+mc) (l+1)(n 1+c) p n 1,m (n+1)(m+l 1) p n+1,m (l+)(l+1) p n,m (l+1)(m 1) p n,m 1 (m+1)(n+c+l 1) p n,m+1 ) (l+)(l+1) p n,m () Considering transitions from states {N c 1,0} and {N c,0}to the absorbingstate {N c,0}, we obtain the decay rate of the survival probability dp s (t)/dt: dp s (t) dt [ ( ) (N 1) = P s (t) p N c 1,0 + p N c,0 N N (4) Substituting Eq. (4) into Eq. (), we finally obtain a condition that the occupation probabilities conditioned on survival must satisfy [0] : p nm = Q nm W nm Q 0 (5) where, Q nm = Q nm (t)/p s (t) is obtained through explicit ] consideration of the terms in the master equation: (l+1)(n 1+c) Q nm = p n 1,m (n+1)(m+l 1) (l+) + p n+1,m + p n,m (l+1)(m 1) (l +) + p n,m 1 + p n,m (m+1)(n+c+l 1) + p n,m+1 (6) [ )+ p N c,0 ] and Q ( 0 = p (N 1) N c 1,0 is the term arising from the decay of the survival probability (Eq. (4)). W nm is the coefficient of p nm ( p nm ) within the brackets on the right hand side of Eq. () (Eq. ()) and is equal to the transition rate out of state {n,m} times N. Eq. (4) indicates that the survival probability decays exponentiallywith arateλ = Q 0 /N. Sincethe meanlifetime of an exponentially decaying process is the inverse ofthe decayrate, it follows that the mean consensus time (neglecting the short transient before the QS distribution is attained) is T c 1 [ ( ) ] (N 1) λ = 1/ p N c 1,0 + p N c,0 N N (7) Thus, knowledge of p nm s (in particular, p N c 1,0 and p N c,0 ) wouldallow us to calculate T c througheq. (7). In order to obtain p nm (for all 0 n,m), we adopt the iterative procedure proposed in [0]. Following this procedure, we start with an arbitrary initial distribution p 0 nm, and obtain a new distribution using: pi+1 nm = α p i nm +(1 α) Qi nm, where 0 α 1 is an arbitrary Wnm i Q i 0

5 5 FIG. : Trajectories (obtained from integration of the mean-field equations Eq. (1)) in the phase-plane show the nature of flows from different regions of the phase-plane into existing fixed points for a) p = 0.05 (< p c) and b) p = 0.1 (> p c) parameter, and Q i nm,wi nm and Q i 0 are all obtained using the probability distribution at the current iteration, p i nm. With sufficient number of iterations, this procedure is expected to converge to a distribution that satisfies Eq. (5) and which is thus, the desired QS distribution. In our case, we obtained acceptable convergence with a choice of α = 0.5 and 0000 iterations. Following the above method, we obtain the QS distribution, and consequently the mean consensus times T c for different values of committed fraction p and system size N. Fig. (a) shows how the consensus time grows as p is decreased beyond the asymptotic critical point p c for finite N. For p < p c, the growth of T c is exponential in N (Fig. (b), consistent with what is known regarding escape times from metastable states. For p > p c, the QS approximation does not reliably provide information on mean consensus times, since consensus times themselves are small and comparable to transient times required to establish a QS state. However, simulation results show that above p c the scaling of mean consensus time with N is logarithmic (Fig. (c)). A snapshot of the QS distribution (Fig. 4) near p c (p = 0.09) for a system of size N = 100 shows clearly the bimodal nature of the distribution, with the two modes centered around the stable fixed point, and the consensus fixed point. The precise dependence of consensus times on p can also be obtained for p < p c by considering the rate of exponential growthof T c with N. In other words, assuming T c exp(α(p)n), we can obtain α(p) as a function of p. Figure (d) shows that α(p) p p c ν where ν Thus, below p c, we have: T c (p < p c ) exp((p c p) ν N) (8) This exponential growth is presumably modulated by factors of logn which become dominant only when p = p c. Above p c, the dependance of T c on p as seen from simulations is negligible (not shown). III. SPARSE NETWORKS Next, we present simulation results for the case when the underlying network topology is chosen from an ensemble of Erdős-Rényi (ER) random graphs with given size N and given average degree k. The qualitative features of the evolution of the system in this case are the same as that of the complete graph, although the critical fraction p c displays some dependence on k. For small k and fixed N, the drop in consensus times occursslightlyearlierinpforergraphsthanforacomplete graph of the same size, as shown in Fig. 5(a). However for p > p c, a complete graph has shorter consensus times (on average) than an ER graph of the same size. Above p c, the difference between consensus times for a graph with average degree k and the complete graph, T c, decays approximately as power law with increasing k (Fig. 5(b)). The deviation from a perfect power law is likely due to other weaker k dependent terms, presumably logarithmic in k. We also performed simulations of the binary agreement model on Barabasi-Albert (BA) networks (Fig. c)), and found similar qualitative behavior as observed for ER networks including the difference from mean-field behavior. We leave a detailed analysis of the dependence of the critical fraction p c and the consensustimes T c onthe averagedegree k ofsparsenetworks, for future work. IV. SUMMARY In closing, we have demonstrated here the existence of a tipping point at which the initial majority opinion of a network switches quickly to that of a consistent and inflexible minority. There are several historical precedents for such events, for example, the suffragette movement in the early 0th century and the rise of the American civil-rights movement that started shortly after the size of

6 6 FIG. : (a) Mean consensus time T c for p < p c obtained by using the QS approximation. (b) Exponential scaling of T c with N, for p < p c; mean consensus times (circles, squares) are obtained using the QS approximation. The lines are guides to the eye. (c) Logarithmic scaling of T c with N for p = 0. > p c; mean consensus times are obtained from simulations. The line shows the best linear fit to the data. (d) The rate, α(p), of exponential growth of the consensus time with N as a function of p p c (see text). Circles show the values of α(p) obtained for p = 0.,0.,0.4,0.5,0.6,0.7,0.8 by considering the scaling of T c with N, for these values of p. The straight line shows a linear fit to the data plotted on a log-log scale. the African-American population crossed the 10% mark. Such processes have received some attention in sociological literature under the term minority influence [6, 1]. Our motivation here has been to study this process in more detail through semi-analytical methods and simulations for finite-sized and sparse networks, within the realm of a particular social influence model - the binary agreement model. There are several open questions and extensions of this work that are worth studying, in our opinion: for example, given a network with non-trivial community structure, what is the optimal scheme for selecting committed agents (for a given committed fraction) that would minimize consensus times, and reduce p c? Secondly, extensions of the model to include utilitydriven opinion switching by agents may be useful in designing optimal incentive schemes for opinion spreading. V. ACKNOWLEDGMENTS This work was supported in part by the Army Research Laboratory under Cooperative Agreement Number W911NF , by the Army Research Office grant W911NF and by the Office of Naval Research Grant No. N The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the official policies either expressed or implied of the Army Research Laboratory or the U.S. Government. S.S. thanks R. Dickman and M. M. de Oliveira for useful discussions.

7 7 FIG. 4: The quasistationary distribution p nm for p = 0.09 and N = 100. Appendix: Fixed points of the mean-field equations Here, we analyze the mean-field equations for the existence of fixed points. To simplify notation we use the notation x = n A and y = n B. Thus for a fixed point of the evolution given by Eq. (1): xy + (1 x y p) +x(1 x y p) + 1.5p(1 x y p) = 0 xy + (1 x y p) +y(1 x y p) yp = 0 (9) which can be reduced further to: x = [(1 y p/4) 9p /16]/(p/+1) y = (1 x p) (10) Substituting the expression of x into the expression for y, and denoting z = y we get z(z ( p/)z +p/+1) = 0 (11) For any value of p, z = z 0 = 0 is a solution to the above equation. In other words, for any value of p, the meanfield equations admit a stable fixed point, n A = x 1 = 1 p, n B = y 1 = 0 which represents the network having reached a consensus state where all nodes have adopted the opinion of the committed agents. The remaining fixed points are roots of f(z) = z ( p/)z +p/+1 = 0. (1) In order to find the criterion which has to be satisfied for valid roots (i.e. 0 z ((1 p))) of the above equation to exist, we analyze the extrema of the function f(z) which are given by: FIG. 5: (a) The steady state density n B of nodes in state B as a function of committed fraction p for Erdős-Rényi graphs of different sizes with k = 10, conditioned on survival of the system. Lines are mean values of n B obtained from simulations; asterisks show mean-field consensus times obtained by integrating Eq. 1 (b) Scaling of T c( k ) (defined in text) with k ; the line shows the best linear fit to the data. (c) The steady state density n B of nodes in state B as a function of committed fraction p for Barabasi-Albert (BA) networks of different sizes with k = 10, conditioned on survival of the system f (z) = z +p/ = 0 (1)

8 8 Thus, the extrema occur at : z 1, = ± / p/6 (14) It can be seen from Eq. (1) that f(z) is increasing, decreasing and increasing again in the intervals (,z 1 ),(z 1,z ),(z,+ ) respectively. Consequently, f(z) achieves a maximum at 1 < z 1 = / p/6 < 0 and a minimum at 0 < z = / p/6 < 1. Furthermore, since f( ) = p/ < 0 and f( 1) = > 0, one root of f(z) = 0 occurs in the interval < z < 1. Sincef(z)ispositiveatz 1, decreasingfromz 1 toz where it achieves a minimum, and increasing thereafter, it follows that a necessary and sufficient condition for more roots of f(z) = 0 to exist, is that f(z ) be less than zero: f(z ) = z ( p/)z +p/+1 < 0 Denoting z = q and p = 4 6q (from Eq. (14)) yields thefollowinginequalityforq asaconditionformoreroots of f(z) = 0 to exist: f(q) = q +1.5q 1.5 > 0 (15) (Note that z is itself a function of p). Analyzing the derivative of f(q) enables us to glean that the inequality Eq. (15) is satisfied for q > q 0 where q 0 is the solution of the cubic equation f(q) = 0 and is given by: [ q 0 = ] 4 1 / Thus, the original fixed point equation Eq. (9) has at least one valid root besides z = 0, so long as p is less than or equal to: p c = 5 ( 5+ ) 4 1 ( 5 ) 4 1 (16) which using standard computer algebra software is evaluated to be p c = Using, z = y = q 0 and Eq. (10), we obtain the state of the system at p c to be {n A,n B } = {0.0957,0.6504}. It also follows from the expression for f(z), that f(0) > 0 and therefore if f(z ) is negative, Eq. (1) has two rootson the positive line when p < p c. Thus there are two fixed points of Eq. (9) when p < p c. The exact expressions for these fixed points (that can also be obtained numerically), obscure their dependence on p. We therefore adopt an approximation which exhibits a much clearer dependence of the fixed point values on p, and numerically yield values close to those obtained from the exact expressions. Substituting z = t p in Eq. (1) reduces it to: 1+p/ t t+r = 0 (17) where r = ( p/). Clearly, r is a monotonically 1 increasing function of p, and therefore r < 1+p c/ = ( p c/) for 0 p < p c, our range of interest. Function g(t) = t t is monotonically decreasing for t < 1 and g( 1) = 0, while g( ) < 1 < r. Hence, there is a real root t 1 (, 1) to Eq. (17) which is a monotonically decreasing function of r, but which clearly does not yield a valid fixed point. This root can be expressed as t 1 (r) = +ǫ(r), where ǫ(r) is monotonically decreasing from less than to 0 over the range of our interest for r. Substituting this expression back into Eq. 17 and neglecting powers of ǫ(r) higher than unity, we get an approximation of ǫ in terms of r, and consequently an approximation for t 1 : t 1 (r) 16 9 r (18) with relative error of less than 0.01%. Now, we can factorize the LHS of Eq. (17) and write it as (t +bt+c)(t t 1 ). Equating this factorized expression with t t+r, gives us b and c in terms of r. Thus, two more roots of Eq. (17) are obtained in terms of r by solving the quadratic equation t +bt+c = 0 which yields: t, = r 6 ± r 9 r 1 (19) Finally, we can obtain the values of z associated with the above roots, and therefore the values of x and y written in terms of these roots are derived as: y, = t 4 p, x, = (4 4y, p) 9p 8p+16 (0) The stability of these fixed points can be checked via linear stability analysis. Following the standard procedure, the stability matrix is given by: [ 1 p S = +y + p +x +p 1 where(x,y )isthefixedpointunderconsideration. The eigenvalues of the stability matrix are given by: λ = 1 4 ( 4 p±(17p + 64(x 1)(y 1) ] + 16p(x +4y 5)) 1 ) (1) From the expression for the eigenvalues we numerically determine that the real part of both the eigenvalues is negative for (x,y ) over the range 0 p < p c indicating that (x,y ) that is a stable fixed point. This is however, not the case for (x,y ), making it unstable. Similarly, the consensus fixed point (x 1,y 1 ) is found to be stable for 0 p 1. Finally, we note that as p 0, the stable fixed point converges to n A = 0,n B = 1, while the unstable fixed point converges to n A = n B 0.8.

9 9 [1] F. Harary, in Studies in Social Power (University of Michigan Press, Ann Arbor, 1959). [] N. E. Friedkin and E. C. Johnson, Journal of Mathematical Sociology 15, 19 (1990). [] T. C. Schelling, Micromotives and Macrobehavior (W. W. Norton, 1978). [4] C. Castellano, S. Fortunato, and V. Loreto, Rev. Mod. Phys. 81, 591 (009). [5] D. Kempe, J. Kleinberg, and E. Tardos, in Proc. 9th ACM SIGKDD international conference on Knowledge discovery and data mining (ACM, New York, NY, USA, 00), KDD 0, p. 17. [6] S. Galam, Physica A 74, 1 (1999). [7] S. Galam, Internat. J. Modern Phys. C 19, 409 (008). [8] M. Granovetter, Am. J. Sociol. 8, 140 (1978). [9] F. M. Bass, Management Science 15, 15 (1969). [10] X. Castelló, A. Baronchelli, and V. Loreto, Eur. Phys. J. B 71, 557 (009). [11] L. Steels, Artificial Life, 19 (1995). [1] A. Baronchelli, M. Felici, E. Caglioti, V. Loreto, and L. Steels, J. Stat. Mech.:Theory Exp. P06014, (006). [1] L. Dall Asta, A. Baronchelli, A. Barrat, and V. Loreto, Phys. Rev. E 74, (006). [14] Q. Lu, G. Korniss, and B. K. Szymanski, Phys. Rev. E 77, (008). [15] A. Baronchelli, Phys. Rev. E 8, (011). [16] D. Stauffer, Comput. Phys. Comm. 146, 9 (00). [17] P. L. Krapivsky and S. Redner, Phys. Rev. Lett. 90, 8701 (00). [18] V. Sood and S. Redner, Phys. Rev. Lett. 94, (005). [19] B. Uzzi, S. Soderstrom, and D. Diermeier (011), submitted to Am. J. Sociol. [0] N. A. Christakis and J. H. Fowler, New Engl. J. Med. 57, 70 (007). [1] M. Kearns, S. Judd, J. Tan, and J. Wortman, Proc. Natl. Acad. Sci. U. S. A 106, 147 (009). [] D. J. Watts and P. S. Dodds, J. Cons. Res. 4, 441 (007). [] Q. Lu, G. Korniss, and B. Szymanski, J. Econ. Interact. Coord. 4, 1 (009). [4] S. Biswas and P. Sen, Phys. Rev. E 80, (009). [5] M. Mobilia, A. Petersen, and S. Redner, J. Stat. Mech.:Theory Exp. 007, P0809 (007). [6] S. Galam and F. Jacobs, Physica A 81, 66 (007). [7] R. Dickman and R. Vidigal, J. Phys. A 5, 1147 (00). [8] R. Dickman, Phys. Rev. E 65, (00). [9] L. Dall Asta and A. Baronchelli, J. Phys. A: Math. Gen 9, (006). [0] M. M. de Oliveira and R. Dickman, Physica A 4, 55 (004). [1] S. Moscovici, E. Lage, and M. Naffrechoux, Sociometry, 65 (1969).