arxiv: v2 [cs.si] 27 Apr 2017

Transcription

1 Opinion Dynamics in Networks: Conergence, Stability and Lack of Explosion arxi: [cs.si] 27 Apr 2017 Tng Mai Georgia Institte of Technology Vijay V. Vazirani Georgia Institte of Technology Abstract Ioannis Panageas MIT-SUTD Inspired by the work of Kempe et al. [16], we introdce and analyze a model on opinion formation; the pdate rle of or dynamics is a simplified ersion of that of [16]. We assme that the poplation is partitioned into types whose interaction pattern is specified by a graph. Interaction leads to poplation mass moing from types of smaller mass to those of bigger. We show that starting niformly at random oer all poplation ectors on the simplex, or dynamics conerges point-wise with probability one to an independent set. This settles an open problem of [16], as applicable to or dynamics. We beliee that or techniqes can be sed to settle the open problem for the Kempe et al. dynamics as well. Next, we extend the model of Kempe et al. by introdcing the notion of birth and death of types, with the interaction graph eoling appropriately. Birth of types is determined by a Bernolli process and types die when their poplation mass is less than a parameter ɛ. We show that if the births are infreqent, then there are long periods of stability in which there is no poplation mass that moes. Finally we show that een if births are freqent and stability is not attained, the total nmber of types does not explode: it remains logarithmic in 1/ɛ. Tng Mai and Vijay V. Vazirani wold like to acknowledge NSF Grant CCF Ioannis Panageas wold like to acknowledge a MIT-SUTD postdoctoral fellowship. 1

2 1 Introdction The birth, growth and death of political parties, organizations, social commnities and prodct adoption grops (e.g., whether to se Windows, Mac OS or Linx) often follows common patterns, leading to the belief that the dynamics nderlying these processes has mch in common. Understanding this commonality is important for the prposes of predictability and hence has been the sbject of stdy in mathematical social science for many years [4, 7, 8, 14, 27]. In recent years, the growth of social commnities on the Internet, and their increasing economic and social ale, has proided fresh impets to this stdy [1, 2, 5, 18]. In this paper, we contine along these lines by bilding on a natral model proposed by Kempe et al. [16]. Their model consists of an inflence graph G on n ertices (types, parties) into which the entire poplation mass is partitioned. Their main tenet is that indiidals in smaller parties tend to get inflenced by those in bigger parties 1. Indiidals in the two ertices connected by an edge can interact with each other. These interactions reslt in indiidals moing from smaller to bigger in poplation ertices. Kempe et al. characterize stable eqilibria of this dynamics ia the notion of Lyapno stability, and they show that nder any stable eqilibrim, the entire mass lies in an independent set, i.e., the poplation breaks into non-interacting islands. The message of this reslt is clear: a poplation is (Lyapno) stable, in the sense that the system does not change by mch nder small pertrbations, only if people of different opinions do not interact. They also showed conergence to a fixed point, not necessarily an independent set, starting from any initial poplation ector and inflence graph. One of their main open problems was to determine whether starting niformly at random oer all poplation ectors on the nit simplex, their dynamics conerge with probability one to an independent set. We first settle this open problem in the affirmatie for a modification of the dynamics, which howeer is similar as that of Kempe et al. in spirit in that it moes mass from smaller to bigger parties (the dynamics is defined below along with a jstification). We beliee that the ideas behind or analysis can be sed to settle the open problem for the dynamics of Kempe et al. as well, ia a more complicated spectral analysis of the Jacobian of the pdate rle of the dynamics (see Section 3.2). Whereas the model of Kempe et al. captres and stdies the effects of migration of indiidals across types in a ery satisfactory manner, it is qite limited in that it does not inclde the birth and death of types. In this paper, we model birth and death of types. In order to arrie at realistic definitions of these notions, we first condcted case stdies of political parties in seeral contries. We present below a case stdy on Greek politics, bt similar phenomena arise in India, Spain, Italy and Holland (see Wikipedia pages). The Siriza party in Greece proides an excellent example of birth of a party (this information is readily aailable in Wikipedia pages). This party was essentially in a dormant state ntil the first 2012 elections in which it got 16.8% of the ote, mostly taken away from the Pasok party, which dropped from 43.9% to 13.2% in the process (Wikipedia). In the second election in 2012, Siriza increased its ote to 26.9% and Pasok dropped to 12.3%. Finally, in 2015, Siriza increased to 36.3% and Pasok dropped frther to 4.7%. Another party, Potami, was formed in 2015 and got 6.1% of the ote, again mainly from Pasok. Howeer, in a major 2016 poll, it seems to hae collapsed and is 1 Changes in the sizes of political parties and other organizations can occr for a mltitde of possible reasons, sch as changes in economic conditions, immigration flows, wars and terrorism, and drastic changes in technology (sch as the introdction of the Internet, smart phones and social media). Stdying changes de to these mltitde of reasons in a systematic qantitatie manner is nrealistic. For this reason, many athors in compter science and the social sciences hae limited their work to stdying the effects of relatie sizes of the grops, in itself a key factor, e.g. see [16]. Following these works, or paper also takes a similar approach. 2

3 likely be absorbed by other parties. In contrast, the KKE party in Greece, which had almost no interactions with the rest of the parties (and was like a disconnected component), has remained between % of the ote oer the last 26 years. Motiated by these examples, we hae modeled birth and death of types in the following manner. We model poplation as a continm, as is standard in poplation dynamics, and time is discrete. This is the same as arxi Version 1 of [16], which is what we will refer to throghot this paper; the later ersions stdy the continos time analog. The birth of a new type in or model is determined by a Bernolli process, with parameter p. The newly born type absorbs mass from all other types ia a randomized process gien by an arbitrary distribtion with finite spport (see Section 2.2). After birth, the new type is connected to an arbitrary, thogh non-empty, set of other types. Or model has a parameter ɛ, and when the size of a type drops below ɛ, it simply dies, moing its mass eqally among its neighbors. Or rle for migration of mass, which is somewhat different from that of Kempe et al. is motiated by the following considerations. For a type, x will denote the fraction of poplation that is of type. Assme that types and hae an edge, i.e., their poplations interact. If so, we will assme that some indiidals of the smaller type get inflenced by the larger one and moe to the larger one. The qestion is what is a reasonable assmption on the poplation mass that moes. For arriing at the rle proposed in this paper, consider three sitations. If x =.02 and x =.25, i.e., the smaller type is ery small, then clearly not many people will moe. If x =.22 and x =.25, i.e., the types are approximately of the same poplation size, then again we expect not many people to moe. Finally, if x =.15 and x =.25, i.e., both types are reasonably big and their difference is also reasonably big, then we expect seeral people to moe from the smaller to the bigger type. From these considerations, we propose that the amont of poplation mass moing from to, assming x < x, is gien by the rle f (t) = F ( ), where F (x) = F (x) is a fnction that captres the leel of inflence between,. We assme that F : [ 1, 1] [ 1, 1] is continosly differentiable, F (0) = 0 (there is no poplation flow between two neighboring types if they hae the same fraction of poplation), is increasing and finally it is odd, i.e., F ( x) = F (x) (so that f (t) = f ). (t) In this simplified setting we hae made the assmption that the system is closed, i.e., that it does not get inflence from otside factors (e.g., economical crisis, immigrations flows, terrorism etc). 1.1 Or reslts and techniqes We first stdy or migration dynamics withot birth and death and settle the open problem of Kempe et al., as it applies to or dynamics. We show that the dynamics conerges set-wise to a fixed point, i.e., there is a set S containing only fixed points sch that the distance between the trajectory of the dynamics and S goes to zero for all starting poplation ectors. To show this conergence reslt, we se a simple potential fnction of the poplation mass namely, the l 2 2 norm of the poplation ector, and we show that this potential is strictly increasing at each time step (nless the dynamics is at a fixed point). Moreoer, the potential is bonded, hence the reslt follows. Next, we strengthen this reslt by showing point-wise conergence as well. The latter reslt is technically deeper and more difficlt, since it means that eery trajectory conerges to a specific fixed point p. We show point-wise conergence by constrcting a local potential fnction that is decreasing in a small neighborhood of the limit point p. The potential fnction is always non-zero in that small neighborhood and is zero only at p. 3

4 Using the latter reslt and one of the most important theorems in dynamical systems, the Center Stable Manifold Theorem, we proe that with probability one, nder an initial poplation ector picked niformly at random from the nit simplex, or dynamics conerges point-wise to a fixed point p, where the actie types w in p, i.e., w V (G) so that p w > 0, form an independent set of G. This inoles characterization of the linearly stable (see Section 2.3 for definition) fixed points and proing that the pdate rle of the dynamics is a local diffeomorphism 2. This settles the open problem of Kempe et al., mentioned in the Introdction, for or dynamics. This reslt is important becase it allows s to perform a long-term aerage case analysis of the behaior of or dynamics and make predictions. Next, we introdce birth and death in or model. Clearly there will be no conergence in this case since new parties are created all the time. Instead we define and stdy a notion of stability which is different from the classical notions that appear in dynamical systems (see Section 2.3 for the definition of the classical notion and Definition 4.2 for or notion). A dynamics is (T, d)-stable if and only if t : T t T + d, no poplation mass moes at step t. We show that despite birth and death, there are arbitrarily long periods of stability with high probability, for a sfficiently small p. Finally, we show that in the long rn, with high probability, for a sfficiently large p, the nmber of types in the poplation will be O(log(1/ɛ)). This may seem conter-intitie, since with a large p new types will be created often; howeer, since new types absorb mass from old types, the old types die freqently. In contrast, in the short term (from the definition of ɛ) we can hae p to Θ(1/ɛ) types. Let s gie an interpretation of the reslts of the preios paragraph in terms of political parties of certain contries (information obtained from Wikipedia). Contries do hae periods of political stability, e.g., dring , , no new major (with more than 1% of the ote) parties were formed in Greece, moreoer there was no sbstantial change in the percentage of otes won by parties in sccessie elections. The parameter ɛ can be interpreted as the fraction of people that can form a party that participates in elections. The minimm size of a party arises for organizational and legal reasons, and is Θ(1/Q), where Q is the poplation of the contry and therefore ɛ is inersely related to poplation. The message of the latter theorem is that the nmber of political parties grows at most as the logarithm of the poplation of the contry, i.e. O(log Q). The following data spports this fact. The poplation of Greece, Spain and India in 2015 was 1.1e7, 4.6e7 and 1.2e9, respectiely, and the nmber of parties that participated in the general elections was 20, 32 and 50, respectiely. 1.2 Related work As stated aboe, we bild on the work of Kempe et al. [16]. They model their dynamics in a similar way, i.e., there is a flow of poplation for eery interacting pair of types,. The flow goes from smaller to bigger types; in or case the mass is jst the poplation of a type. One ery interesting common trait between the two dynamics is that the fixed points hae similar description: all types with positie mass belonging to the same connected component C hae the same mass. Stable fixed points also hae the same properties in both dynamics, namely they are independent sets. The pdate rles of the two dynamics are somewhat different; or simpler dynamics helps s in proing stronger reslts. One of the most stdied models is the following: there is a graph G in which each ertex denotes an indiidal haing two possible opinions. At each time step, an indiidal is chosen at random who next chooses his opinion according to the majority (best response) opinion among his neighbors. 2 Continosly differentiable, the inerse exists and is also continosly differentiable. 4

5 This has been introdced by Galam [10] and appeared in [9, 23], where they address the qestion: in which classes of graphs do indiidals reach consenss. The same dynamics, bt with each agent choosing his opinion according to noisy best response (the dynamics is a Marko chain) has been stdied in [17, 22] and many other papers referenced therein. They gie bonds for the hitting time and expected time of the consenss state (risk dominant) respectiely. Another well-known model for the dynamics of opinion formation in mltiagent systems is Hegselmann-Krase [11]. Indiidals are discrete entities and are modeled as points in some opinion space (e.g., real line). At eery time step, each indiidal moes to the mass center of all the indiidals within nit distance. Typical qestions are related to the rate of conergence (see [6] and references therein). Finally, another classic model is the oter model, where there is a fixed graph G among the indiidals, and at eery time step, a random indiidal selects a random neighbor and adopts his opinion [12]. For more information on opinion formation dynamics of an indiidal sing information learned from his neighbors, see [13] for a srey. Other works, inclding dynamical systems that show conergence to fixed points, are [19 21, 24, 25, 28]. [28] focses on qadratic dynamics and they show conergence in the limit. On the other hand [3] shows that sampling from the distribtion this dynamics indces at a gien time step is PSPACE-complete. In [20, 24], it is shown that replicator dynamics in linear congestion and 2-player coordination games conerges to pre Nash eqilibria, and in [19, 25] it is shown that gradient descent conerges to local minima, aoiding saddle points een in the case where the fixed points are ncontably many. Organization: In Section 2 we describe or dynamics formally and gie the necessary definitions abot dynamical systems. In Section 3 we show that or dynamics withot births/deaths conerges with probability one to fixed points p so that the set of types with positie poplation, i.e., actie types, form an independent set of G. Finally, in Section 4 we first show that there is no explosion in the nmber of types (i.e., the order neer becomes Θ(1/ɛ)) and also we perform stability analysis sing or notion. 2 Preliminaries Notation: We denote the probability simplex on a set of size n as n. Vectors in R n are denoted in boldface and x j denotes the jth coordinate of a gien ector x. Time indices are denoted by sperscripts. Ths, a time indexed ector x at time t is denoted as. We se the letters J, J to denote the Jacobian of a fnction and finally we se f t to denote the composition of f by itself t times. 2.1 Migration dynamics Let G = (V, E) be an ndirected graph on n ertices (which we also call types), and let N denote the set of neighbors of in G. Dring the whole dynamical process, each ertex has a non-negatie poplation mass representing the fraction of the poplation of type. We consider a discrete-time process and let denote the mass of time step t. It follows that the condition V (G) = 1, mst be maintained for all t, i.e., n 3 for all t N. 3 Recall that n denotes the simplex of size V (G) = n. 5

6 Additionally, we consider a dynamical migration rle where the poplation can moe along edges of G at each step. The moement at step t is determined by. Specifically, for E(G), the amont of mass moing from to at step t is gien by f (t) = F ( ). For all E(G) we assme that F : [ 1, 1] [ 1, 1] is a continosly differentiable fnction sch that: 1. F (0) = 0 (there is no poplation flow between two neighboring types if they hae the same fraction of poplation), 2. F is increasing (the larger x x, the more poplation moing from to ), 3. F is odd i.e., F ( x) = F (x) (so that f (t) = f (t) ). It can be easily deried from the assmptions that F (x) F (0) = 0 for x 0 and F ( x) = F (x) for x [ 1, 1], where F denotes the deriatie of F. Note that f (t) > 0 implies that poplation is moing from to, and f (t) < 0 implies that poplation is moing in the other direction. The pdate rle for the poplation of type can be written as x (t+1) = = + f (t) (1) N + N F ( We denote the pdate rle of the dynamics as g : n n, i.e., we hae that x (t+1) = g( ). ). (2) Therefore it holds that = g t (x (0) ), where g t denotes the composition of g by itself t times. It is easy to see g is well-defined for sp x [ 1,1] F (x) 1 for all E(G), in the sense that if n then x (t+1) n. This is tre since for all we get (sing indction, i.e., n ) x (t+1) = x (t) + F ( ) N N (1 ) 0, moreoer it holds x (t+1) = + N + N + (1 ) + 1 = 1, F ( ) and also x (t+1) = = 1 (the other terms cancel ot). 6

7 2.2 Birth and death of types Political parties or social commnities don t tend to srie once their size becomes small and hence there is a need to incorporate death of parties in or model. We will define a global parameter ɛ in or model. When the poplation mass of a type becomes smaller than some fixed ale ɛ, we consider it to be dead and moe its mass arbitrarily to existing types. Formally, if ɛ then 0 and + / N for all N. Also, ertex is remoed and edges are added arbitrarily on its neighbors to ensre connectiity of the reslting graph. Remark 2.1. It is not hard to see that the maximm nmber of types is 1/ɛ (by definition). We say that we hae explosion in the nmber of types if they are of Θ(1/ɛ). In Theorem 4.6 we show that in the long rn, the nmber of types is mch smaller, it is O(log(1/ɛ))) with high probability. Eery so often, new political opinions emerge and like-minded people moe from the existing parties to create a new party, which then follows the normal dynamics to either srie or die ot. To model birth of new types, at each time step, with probability p, we create a new type sch that takes a portion of mass from each existing type independently. The amont of mass going to from each follows an arbitrary distribtion in the range [β min, β max ]. Specifically, let Z D where D is a distribtion with spport [β min, β max ], the amont of mass going from to is Z x. We connect to the existing graph arbitrarily sch that it remains connected. Additionally, we make a small change to the migration dynamics defined in Section 2.1 to make it more realistic. Or tenet is that poplation mass migrates from smaller to bigger types becase of inflence. Howeer, if the two types are of approximately the same size, the difference is size is not discernible and hence migration shold not happen. To incorporate this, we introdce a new parameter δ > 0 and if x x δ, we assme that no poplation moes from to. Finally, each step of the dynamics consists of there phases in the following order: 1. Migration: the dynamics follows the pdate rle from Section Birth: with probability p, a new type is created and takes mass from the existing types. 3. Death: a type with mass smaller than ɛ dies ot and moe its mass to the existing types. Remark 2.2. For any different order of phases, all proofs in the paper still go throgh with minimal changes. 2.3 Definitions and basics A recrrence relation of the form x (t+1) = f( ) is a discrete time dynamical system, with pdate rle f : S S (for or prposes, the set S is n ). The point z is called a fixed point or eqilibrim of f if f(z ) = z. A fixed point z is called Lyapno stable (or jst stable) if for eery ε > 0, there exists a ζ = ζ(ε) > 0 sch that for all x with x z < ζ we hae that f k (x ) z < ε for eery k 0. We call a fixed point z linearly stable if, for the Jacobian J(z ) of f, it holds that its spectral radis is at most one. It is tre that if a fixed point z is stable then it is linearly stable bt the conerse does not hold in general [26]. A seqence (f t (x (0) )) t N is called a trajectory of the dynamics with x (0) as starting point. A common techniqe to show that a dynamical system conerges to a fixed point is to constrct a fnction P : m R sch that P (f(x )) > P (x ) nless x is a fixed point. We call P a potential or Lyapno fnction. 7

8 3 Conergence to independent sets almost srely In this section we proe that the deterministic dynamics (assming no death/birth of types, namely the graph G remains fixed) conerges point-wise to fixed points p where { : p > 0} (set of actie types) is an independent set of the graph G, with probability one assming that the starting point x (0) follows an atomless distribtion with spport in n. To do that, we show that for all starting points x (0), the dynamics conerges point-wise to fixed points. Moreoer we proe that the pdate rle of the dynamics is a diffeomorphism and that the linearly stable fixed points p of the dynamics satisfy the fact that the set of actie types in p is an independent set of G. Finally, or main claim of the section follows by sing a well-known theorem in dynamical systems, called Center-Stable Manifold theorem. Strctre of fixed points. The fixed points of the dynamics (1) are ectors p sch that for each E(G), at least one of the following conditions mst hold: 1. p = p, 2. p = 0, 3. p = 0. Therefore, for each fixed point p, the set of actie types (types with non-zero poplation mass) with respect to p mst form a set of connected components sch that all types in each component hae the same poplation mass. We first proe that the dynamics conerges point-wise to fixed points. 3.1 Point-wise conergence First we consider the following fnction Φ(x ) = x 2 and state the following lemma on Φ. Lemma 3.1 (Lyapno (potential) fnction). Let x be a ector with x > x. Let y be another ector sch that y = x d, y = x + d for some 0 < d x and y z = x z for all z,. Then Φ(x ) < Φ(y). Proof. By the definition of Φ, Φ(y) = y 2 + y 2 + z, y 2 z = (x d) 2 + (x + d) 2 + z, y 2 z = x 2 2dx + d 2 + x 2 2dx + d 2 + = Φ(x ) + 2d(x x ) + 2d 2 > Φ(x ). The ineqality follows becase d > 0 and x > x. If we think of x as a poplation ector, Lemma 3.1 implies that Φ(x ) increases if poplation is moing from a smaller type to a bigger type. z, x 2 z 8

9 Theorem 3.2 (Set-wise conergence). Φ( ) is strictly increasing along eery nontriial trajectory, i.e., Φ(x t+1 ) = Φ(g( )) Φ( ) with eqality only when is a fixed point. As a corollary, the dynamics conerges to fixed points (set-wise conergence). Proof. First we proe that the dynamical process conerges to a set of fixed points by showing that Φ( ) is strictly increasing as t grows. The idea is breaking a migration step from to x (t+1) into mltiple steps sch that each small step only inoles migration between two types. Moreoer, in each small step, poplation is moing from a smaller type to a bigger type. Lemma 3.1 garantees that Φ is strictly increasing in eery small step, and ths strictly increasing in the combined step from to x (t+1). Let D be the directed graph representing the migration moement from to x (t+1). Formally, for each edge E(G) we direct in both directions and let f = max(0, f ). (t) Define the following process on D: 1. If there exists a directed path z of length 3 in D sch that f and f z are both positie, we make the following modification to the flow in D. Let = min(f, f z ), and f f f z f z f z f z Keep repeating the preios step ntil there is no path of length 3 carrying positie flow. The aboe process mst terminate since the fnction f ( ) 2 strictly increases in each modification and is bonded aboe. At the end of the process, there is no path of length 3 carrying positie flow. In other words, each type can not hae both flows coming into it and flows coming ot of it. Frthermore, it is easy to see that the net flow V (G) f at each type is presered and f > 0 only if >. We can break a migration step from to x (t+1) into mltiple migrations sch that each small migration corresponds to a flow f > 0 at the end of the process. It follows that in each small migration, poplation is moing from one smaller type to one bigger type. To finish the proof we proceed in a standard manner as follows: Let Ω n be the set of limit points of a trajectory (orbit). Φ( ) is increasing with respect to time t by aboe and so, becase Φ is bonded on n, Φ( ) conerges as t to Φ = sp t {Φ( )}. By continity of Φ we get that Φ(p) = lim t Φ( ) = Φ for all p Ω, namely Φ is constant on Ω. Also p (t) = lim k x (t k+t) as k for some seqence of times {t i }, with p (0) = lim k x (t k) and p (0) Ω. Therefore p (t) lies in Ω for all t N, i.e., Ω is inariant. Ths, since p (0) Ω and the orbit p (t) lies in Ω, we get that Φ(p (t) ) = Φ on the orbit. Bt Φ is strictly increasing except on fixed points and so Ω consists entirely of fixed points. Theorem 3.3 (Point-wise conergence). The dynamics conerges point-wise to fixed points. Proof. We first constrct a local potential fnction Ψ sch that Ψ(x, p) is strictly decreasing in some small neighborhood of a limit point (fixed point) p. Formally we initially proe the following: 9

10 Claim 3.4 (Local Lyapno fnction). Let p be a limit point (which will be a fixed point of the dynamics by Theorem 3.2) of trajectory y (t), and Ψ(x, p) be the following fnction Ψ(x, p) = :p >0 (p x ). There exists a small ε > 0 sch that if y (t) p 1 ε then Ψ(y (t+1), p) Ψ(y (t), p). Proof of Claim 3.4. We know that the set of actie types in p mst form a set of connected components sch that all types in each component hae the same poplation mass with respect to p. Let C be one sch component and let δ(c) = { : C, C, E(G)}. We choose ε to be so small so that y (t) > y (t) for all C and δ(c), becase y (t) < ε since p = 0 (ths is arbitrarily close to zero) for δ(c). Therefore, the net flow into C mst be non-negatie, ths Ψ(g(y (t) ), p) Ψ(y (t), p). Additionally, we hae that y (t) p for all C. Sppose otherwise, then there exists a type w def = w(t) C with y (t) w > p w, which is a contradiction. To see why, consider w(t ) to be arg max z C y (t ) z, then it shold hold that y (t ) w(t ) p w(t ) + s for some constant s > 0 independent of t and t t, becase y (t ) w(t ) is increasing and therefore p cannot be a limit point of the trajectory y (t). Hence, Ψ(y (t), p) mst be non-negatie and only eqal to zero when y (t) = p (i). To finish the proof of the theorem, if p is a limit point of y (t), there exists an increasing seqence of times {t i }, with t n and y (tn) p. We consider ε sch that the set S = {x : x p and Ψ(x, p) ε } is inside B = x p 1 < ε where ε is from Claim 3.4 abot the local potential. Since y (tn) p, consider a time t N where y (t N ) is inside S. From Claim 3.4 we get that if y (t) B then Ψ(y (t+1), p) Ψ(y (t), p) (and also y (t) S), ths Ψ(y (t) ) Ψ(y (t N ), p) ε and y (t) p for all t t N (namely the orbit remains in S; we se Claim 3.4 indctiely). Therefore Ψ(, p) is decreasing in S and since Ψ(y (tn), p) Ψ(p, p) = 0, it follows that Ψ(y (t), p) 0 as t. Hence y (t) p as t sing (i). 3.2 Diffeomorphism and stability analysis ia Jacobian In this section we compte the Jacobian J of g and then perform spectral analysis on J. The Jacobian of g is the following: g = J, = 1 + x N x [ F (x x ) + x F (x x ) ], g x = J, = x [ F (x x ) x F (x x ) ] if E(G) else 0. Lemma 3.5 (Local Diffeomorphism). The Jacobian is inertible on the sbspace x = 1, for sp x [ 1,1] F (x) < 1 2 for each E(G). Moreoer, g is a local diffeomorphism in a neighborhood of n. Proof. First we hae that N x F (x x ) + N x x F (x x ) N x 1 and hence J, > 0 for all and x n. The first ineqality comes from the fact that F (x) 0 10

11 (F is increasing) and F (x) > 1/2 > 1. Additionally, we get that J, J, = 1 + = N x > 1 N,J,>0 N,J,>0 N x 0, [ F (x x ) + x F (x x ) ] J,>0 x [ F (x x ) + x F (x x ) ] x F (x x ) J + J,<0 where we sed the fact that 1 2 > F (x) > 1 2, F ( x) = F (x) and that N x 1. Therefore we conclde that J is diagonally dominant. Finally, assme that J is not inertible, then there exists a nonzero ector y so that J y = 0 (ii). We consider the index type with the maximm absolte ale in y, say w. Hence we hae y w y for all V (G). Finally, sing (ii) we hae that J w,w y w = w J,wy, ths J w,w w J w, y y w w J w, (first ineqality is triangle ineqality and second comes from assmption on w). We reach a contradiction becase we showed before that J is diagonally dominant. Therefore, J(x ) is inertible for all x n. Moreoer, from inerse fnction theorem we get the claim that the pdate rle of the dynamics is a local diffeomorphism in a neighborhood of n. Lemma 3.6 (Linearly stable fixed point independent set). Let p be a fixed point so that there exists a connected component C, C > 1 with p > 0 same poplation mass for all C. Then the Jacobian at p has an eigenale with absolte ale greater than one. Proof. The Jacobian at p has eqations: 1. Assme p = 0 then J, = 1 + N p F ( p ) = 1 J, = 0 for all. N p F (p ), J 2. Assme p > 0 then J, = 1 + F (0)p 2. N,p >0 For all V (G) so that p = 0, it follows that J, is nonzero only when = (diagonal entry) and ths the corresponding eigenale of J is 1 N p F (p ) 1 with left eigenector (0,..., 0, }{{} 1, 0,..., 0) (it is clear that p F (p ) 0 since F (x) F (0) = 0 for x 0). Hence -th the characteristic polynomial of J at p is eqal to λ 1 p F (p ) det(λi J), :p =0 N 11

12 where J corresponds to J at p by deleting rows corresponding to types sch that p = 0 and I the identity matrix. So it sffices to proe that J has an eigenale with absolte ale greater than one. Assme J has size l l, in other words the nmber of actie types is l in p. Eery type that has no actie neighbors satisfies J, = 1 and eery type that has at least one actie neighbor satisfies J, = 1 + N,p >0 F (0)p 2 > 1. Therefore trace(j) > l by assmption on p. Hence the sm of the eigenales of J is greater than l, ths there exists an eigenale with absolte size greater than Center-stable manifold and aerage case analysis In this section we proe or first main reslt, Corollary 3.11, which is a conseqence of the following theorem: Theorem 3.7. Assme that max x [ 1,1] F (x) < 1/2 for all E(G). The set of points x n sch that dynamics 1 starting at x conerges to a fixed point p whose actie types do not form an independent set of G has measre zero. To proe Theorem 3.7, we are going to se argably one of the most important theorems in dynamical systems, called Center Stable Manifold Theorem: Theorem 3.8 (Center-stable Manifold Theorem [29]). Let p be a fixed point for the C r local diffeomorphism f : U R m where U R m is an open neighborhood of p in R m and r 1. Let E s E c E be the inariant splitting of R m into generalized eigenspaces of the Jacobian of g, J(p) corresponding to eigenales of absolte ale less than one, eqal to one, and greater than one. To the J(p) inariant sbspace E s E c there is an associated local f inariant C r embedded disc Wloc sc tangent to the linear sbspace at p and a ball B arond p sch that: f(w sc loc ) B W sc loc. If f m (x ) B for all m 0, then x W sc loc. (3) Since an n-dimensional simplex n in R n has dimension n 1, we need to take a projection of the domain space ( x = 1) and accordingly redefine the map g. Let x be a point mass in n. Let be a fixed type and define h : R n R n 1 so that we exclde the ariable x from x, i.e., h(x) = x. We sbstitte the ariable x with 1 x and let g be the reslting pdate rle of the dynamics g (x ) = g(x ). The following lemma gies a relation between the eigenales of the Jacobians of fnctions g and g. Lemma 3.9. Let J, J be the Jacobian of g, g respectiely. Let λ be an eigenale of J so that λ does not correspond to left eigenector (1,..., 1) (with eigenale 1). Then J has also λ as an eigenale. Proof. By chain rle, the eqations of J are as follows: g = J x,w(x ) = J,w (x, 1 x ) J, (x, 1 x ). w Assme λ is associated with left eigenector r def = (r 1,..., r n 1, r n ) (we label the types with nmbers 1,..., n with taking index n). We claim that λ is an eigenale of J with right eigenector 12

13 r def = (r 1 r n,..., r n 1 r n ). First it is easy to see that n 1 n 1 J ji(r j r n ) = (J ji J jn )(r j r n ) j=1 j=1 n 1 = j=1 n 1 J ji r j j=1 n 1 J jn r j j=1 n 1 J ji r n + j=1 J jn r n. Since r is a left eigenector, we get that n j=1 J jir j = λr i for all i [n] and also it holds that n j=1 J ji = 1 for all i [n]. Therefore n 1 J ji(r j r n ) = (λr i J ni r n ) (λr n J nn r n ) (1 J ni )r n + (1 J nn )r n j=1 = λ(r i r n ), namely r J = λr and the lemma follows. Before we proceed with the proof of Theorem 3.7, we state the following which is a corollary of Lemmas 3.5, 3.6 and 3.9 and also ses classic properties for determinants of matrices. Corollary Let p be a fixed so that the actie types are not an independent set in G, then J at h(p) has an eigenale with absolte ale greater than one. Additionally, The Jacobian J of g is inertible in h( n ) and as a reslt g is a local diffeomorphism in a neighborhood of h( n ). Proof. If p is a fixed point where the actie types are not an independent set in G, then by Lemma 3.6 we get that J at p has an eigenale with absolte ale greater than one, hence sing Lemma 3.9 it follows that J at h(p) has an eigenale with absolte ale greater than one. Let B be the reslting matrix if we add all the first n 1 rows to the n-th row and then sbtract the n-th colmn from all other colmns in matrix J. It is clear that det(b) = det(j) 0 (determinant not zero since J is inertible from Lemma 3.5). Additionally, the last row of B is all 0 s and B nn = 1, so det(b) = det(b ) where B is the reslting matrix if we delete from B last row,colmn. Bt B = J, hence 0 det(j) = det(b) = det(j ) and ths J is inertible, therefore g is a local diffeomorphism in a neighborhood of h( n ) (by Inerse fnction theorem). Proof of Theorem 3.7. Let p be a fixed point of fnction g(x ) so that the set of actie types is not an independent set. We consider the projected fixed point p def = h(p) of fnction g. Then p is a linearly nstable fixed point. Let B p be the (open) ball (in the set R n 1 ) that is deried from center-stable manifold theorem. We consider the nion of these balls A = p B p. De to Lindelőf s lemma A.1 stated in the appendix, we can find a contable sbcoer for A, i.e., there exist fixed points p 1, p 2,... sch that A = m=1 B p. Starting from a point x def = h( m n ), there mst exist a t 0 and m so that g t (x) B p m for all t t 0 becase of Theorem 3.3, i.e., becase the dynamics (1) conerges point-wise. From center-stable manifold theorem we get that g t (p ) Wloc sc (p m) where we sed the fact that g ( ) (the poplation ector is always in simplex, see Section 2.1), namely the trajectory remains in for all times 4. 4 W sc loc(p m ) denotes the center stable manifold of fixed point p m 13

14 By setting D 1 (p m) = g 1 (W sc loc (p m) ) and D i+1 (p m) = g 1 (D i (p m) ) we get that x D t (p m) for all t t 0. Hence the set of initial points in so that dynamics conerges to a fixed point p so that the set of actie types is not an independent set of G is a sbset of P = m=1 t=0 D t (p m)). (4) Since p m is linearly nstable fixed point (for all m), the Jacobian J has an eigenale greater than 1, and therefore the dimension of Wloc sc (p m) is at most n 2. Ths, the set Wloc sc (p m) has Lebesge measre zero in R n 1. Finally since g is a local diffeomorphism, g 1 is locally Lipschitz (see [26] p.71). g 1 preseres the nll-sets (by Lemma A.2 that appears in the appendix) and hence (by indction) D i (p m) has measre zero for all i. Thereby we get that P is a contable nion of measre zero sets, i.e., is measre zero as well. Corollary 3.11 (Conergence to independent sets). Sppose that max x [ 1,1] F (x) < 1/2 for all E(G). If the initial mass ector x (0) n is chosen from an atomless distribtion, then the dynamics conerges point-wise with probability 1 to a point p so that the actie types form an independent set in G. Proof. The proof comes from Theorem 3.3 and Theorem 3.7. From Theorem 3.3 we hae that lim t exists for all x (0) n and from Theorem 3.7 we get the probability that the dynamics conerges to fixed points where the actie types are not an independent set is zero. Hence the dynamics conerges to fixed points where the actie types is an independent set, with probability one. Corollary 3.11 is illstrated in Figre 1 for the case of a 3-path and a triangle. As shown in the figre, if the initial condition is chosen niformly at random from a point in the simplex, the dynamics conerges to an independent set with probability one. 4 Stability and bond on the nmber of types In this section we consider dynamical systems with migration, death and birth and proe two probabilistic statements on stability and nmber of types. Recall that in this settings, the dynamics at each step contains three phases in order: migration birth death. The following direct application of Chernoff s bond is sed intensiely to attain probabilistic garantees. Lemma 4.1. Consider a period of t steps. 1. There are at least tp/2 births with probability at least 1 e tp/8. 2. There are at most 3tp/2 births with probability at least 1 e tp/6. Proof. Let X i = 1 if there is a birth at step i, and X i = 0 otherwise. The nmber of births in k step is k X = It follows that E[X] = tp. From Chernoff s bond, i=1 X i Pr (X tp/2) e tp/8. In other words, Pr (X tp/2) 1 e tp/8. 14

15 Applying Chernoff s bond again, we hae Pr (X 3tp/2) e tp/6, and Pr (X 3tp/2) 1 e tp/ Stability We define the notion of stability and gie a stability reslt for a dynamical system inoling migration, death and birth. For the rest of the paper, we denote by α min = min E(G),x [ 1,1] F (x) and α max = max E(G),x [ 1,1] F (x). It can be seen easily that for each E(G), α min x F (x) F (0) α max x. Definition 4.2 ((T, d)-stable dynamics). A dynamics is (T, d)-stable if and only if T t T +d, no poplation mass moes in the migration phase at step t. We state the following two lemmas whose proofs come from the definition of Φ. Lemma 4.3. If the dynamics is not (t, 0)-stable, the migration phase at time t increases Φ by at least 2α min ɛδ 3. Proof. By Theorem 3.2, we know that Φ is strictly increasing in each migration when the dynamics is not at a fixed point. Moreoer, it is easy to see that the increase is minimized when there is only one edge carrying poplation flow. From the proof of Lemma 3.1, the additie increase in Φ is at least 2d x x, where d is the amont of mass moing along. Withot loss of generality, we may assme that x > x. Since x ɛ and x x δ, d = x x F (x x ) α min x x (x x ) α min ɛδ 2. It follows that 2d x x 2α min ɛδ 3. Lemma 4.4. Each birth can decrease Φ by at most 2β max. Proof. Recall that Φ(x ) = x 2. The potential after a new type is created is x 2 (1 Z ) 2 + Z 2 x 2. 15

16 Therefore, the net decrease in potential is Φ = x 2 x 2 (1 Z ) 2 Z 2 x 2 = x 2 x 2 (1 2Z + Z 2 ) Z 2 x 2 = 2 2 Z x 2 Z 2 x 2 Z x 2 2β max x 2 Z 2 x 2 2β max ( x ) 2 = 2β max. With the two aboe lemmas, we can gie a theorem on the stability of the dynamics. At a high leel, it says that if the probability of a new type emerging is small enogh, then as time goes on, there mst be a long period period sch that there is no migration. ( Theorem 4.5 ( Stable for long enogh). Let p < min ɛδ 3 α min ( With probability at least 1 e tp/6, the dynamics is T, 1 3p 3β max, 2 3 ) ) -stable for some T t. and t > 1 ɛδ 3 α min 3pβ max. Proof. Consider a period of t steps. By Lemma 4.1, there are at most 3tp/2 births in the period with probability at least 1 e tp/6. Note that p < 2/3 garantees that 3tp/2 < t. In the migration phases of the period, Φ can either increase if there is a migration or remain nchanged otherwise. Assme that Φ increases in more than t/2 migration phases. By Lemma 4.3, the amont of increase in potential de to migrations is at least t 2 2α minɛδ 3 = tα min ɛδ 3. Since there are at most 3tp/2 births, Lemma 4.4 garantees that the amont of decrease in potential de to births is at most 3tp 2 2β max = 3tpβ max. Therefore, the net increase of Φ is least tα min ɛδ 3 3tpβ max = t(α min ɛδ 3 3pβ max ). 1 Since t > ɛδ 3 α min 3pβ max, the net increase in Φ is greater than 1, which is a contradiction. It follows that Φ cannot increase in more than t/2 migration phases, and mst remain nchanged in at least t/2 migration phases. Note that if the dynamics is (t, 0)-stable for some t, it will be (t, d)-stable ntil the next birth at t = t + d + 1. Since there are at most 3tp/2 births, there mst be no migration in a period of t/2 3tp/2 = 1 3p consectie steps. 16

17 4.2 Bond on the nmber of types In this section we inestigate a behaior of the dynamics following a long period of time. Specifically, we show that after a large nmber of steps, the nmber of types can not be too high. Or goal is to proe the following theorem: Theorem 4.6 (Lack of explosion). Let α max p/512 and t (16/p) log 2 (1/ɛ). The dynamics at step t has at most 72 log(1/ɛ) types with probability at least 1 3ɛ. First we gie the following lemma, which says that if the nmber of types is large enogh, then after a fixed period of time, it will decrease by a factor of roghly 2. Lemma 4.7. Let α max p/512 and k be the nmber of types at step t 0. If k 48 log(1/ɛ), with probability at least 1 2ɛ 2, the nmber of types at step t 0 +(16/p) log(1/ɛ) is at most k/2+24 log(1/ɛ). Proof. Withot loss of generality, we may assme t 0 = 0. Or goal is to show that with high probability, at least half of the original types (types at step 0) will die from step 0 to step (16/p) log(1/ɛ). From Lemma 4.1, with probability at least 1 2ɛ 2, there are at least 8 log(1/ɛ) and at most 24 log(1/ɛ) births in that period. Therefore, we may assme that the nmber of births in the period is between 8 log(1/ɛ) and 24 log(1/ɛ). Since the nmber of new types created is at most 24 log(1/ɛ), if at least half of the original types die ot from step 0 to step (16/p) log(1/ɛ), the lemma will immediately follow. Assme that the nmber of original types dying in the period is less than k/2 for the sake of contradiction. It follows that the nmber of remaining original types at the end of the period is at least k/2. Since 24 log(1/ɛ) k/2, the total nmber of types dying ot in the period is at most k/ log(1/ɛ) k. Therefore, the total mass that can be added to the remaining types from the dead types is at most kɛ. It follows that on aerage, a remaining original type receie at most 2ɛ from the dead types. Marko s ineqality says that at most half of the remaining original types can receie more than 4ɛ from the dead types. Therefore, at most k/4 original types receie more than 4ɛ from the dead types. Since the aerage mass of an original type is 1/k, Marko s ineqality again garantees that the nmber of original types haing mass greater than 4/k is at most k/4. Combining the aboe two reasons, we can conclde that at least k/2 original types hae initial mass less than 4/k and receie less than 4ɛ from the dead types in the whole period. We will proe that those types can not remain at the end of the period. The idea is to bond the total increase in the masses of those types, and arge that after at least 8 log(1/ɛ) births, their masses will all be less than ɛ. Consider an original type sch that x (0) 4/k and receies less than 4ɛ from the dead types from step 0 to step (16/p) log(1/ɛ). Recall that the change in mass of at step t is = N α max N ( N α max ( F ) 3 ( (k + 24 log(1/ɛ)) α max ( < 2kα max ( ) 3. ) ) ) 3 17

18 It follows that x (t+1) + 2kα max ( Ignoring the effect of births from step 0 to step t < (16/p) log(1/ɛ), we claim that 8 ( ) 8 3 k + 16α maxkt. k We will proe the aboe claim by indction. We may assme that receies all the mass from dead types at the beginning, since this assmption can only increase. We hae x (0) 4 k + 4ɛ 8 k, where the last ineqality follows since k 1/ɛ. Hence, the base case is satisfied. Sppose that the claim is tre for t, then we hae ( ) 3 x (t+1) + 2kα max 8 ( 8 k + 16α maxkt k 8 ( 8 k + 16α maxkt k = 8 k + 16α maxk(t + 1) The last ineqality follows becase and ths, ) 3 + 2kα max ( 8 ) 3 + 2kα max ( 2 8 k ( ) 8 3. k ) 3. k + 16α maxkt k 48 log(1/ɛ) 48 log(1/ɛ)(512α max /p) > 1024α max t ) 3 8 k = k 8 k 2 > 1024α maxt 8 k 2 = 16α maxtk ( ) 8 3. k ( ) ) k Now the mass of decreases by at least a mltiplicatie factor of 1 β min at each birth. We may assme that the decreases on x happen after the increases since this assmption can only increase the bond on x. We hae ( 8 k + 16α maxkt ( ) ) 8 3 (1 β min ) B, k where B is the nmber of births in the period of t steps. Setting t = (16/p) log(1/ɛ) and B = 8 log(1/ɛ) gies ( 8 Therefore, k + 16α maxk(16/p) log(1/ɛ) ( 8 k log(1/ɛ) 1 ) k 2 ɛ = 16ɛ k. < ɛ as desired. 18 ( ) ) 8 3 (1 β min ) 8 log(1/ɛ) k

19 Now we can proe Theorem 4.6. Proof of Theorem 4.6. We only consider the last (16/p) log 2 (1/ɛ) steps and assme that t = (16/p) log 2 (1/ɛ). We call a period of (16/p) log(1/ɛ) steps a decreasing period if it satisfies the condition in Lemma 4.7, i.e., if the nmber of types k at the beginning of the period is at least 48 log(1/ɛ), and the nmber of types at the end of the period is at most k/ log(1/ɛ). Constrct a set P of periods of length (16/p) log(1/ɛ) as follows: 1. Start with t = Repeat the following step ntil t = t: (a) If t + (16/p) log(1/ɛ) t and the nmber of types at t is at least 48 log(1/ɛ), let i be the period from t to t + (16/p) log(1/ɛ), and add i to P. Update t t + (16/p) log(1/ɛ). (b) Else pdate t t + 1. Assme that all periods in P are decreasing periods. By Lemma 4.7, the probability of sch an otcome occrring is at least (1 2ɛ 2 ) log(1/ɛ) 1 2ɛ 2 log(1/ɛ) 1 2ɛ. First, we proe that the nmber of types mst become smaller than 48 log(1/ɛ) at some step of the dynamics. Assme that the nmber of types is at least 48 log(1/ɛ) throghot the dynamics. Since the dynamics has (16/p) log 2 (1/ɛ) steps and each period has (16/p) log(1/ɛ) steps, there are log(1/ɛ) periods in P. Let T (n) be the nmber of types after the n-th period is added to P. Since all periods in P are decreasing periods, we hae the recrrence T (n) T (n 1)/ log(1/ɛ). Note that T (0) 1/ɛ in the base case. The aboe recrrence gies T (log(1/ɛ)) < 48 log(1/ɛ). In other words, the nmber of types becomes smaller than 48 log(1/ɛ) at the end of the dynamics, which is a contradiction. Now, we know that the nmber of types becomes smaller than 48 log(1/ɛ) at some step. If it remains smaller than 48 log(1/ɛ), the theorem holds triially. Therefore, we may assme that it reaches 48 log(1/ɛ) at some later step. We consider two cases: 1. There are at least (16/p) log(1/ɛ) sbseqent steps. According to or assmption, that period of (16/p) log(1/ɛ) steps mst be a decreasing period. Therefore, the nmber of types is at most 48 log(1/ɛ) after the period. 2. There are less than (16/p) log(1/ɛ) sbseqent steps after the nmber of types reaches 48 log(1/ɛ). By Lemma 4.1, the probability that in those remaining steps, there are at most 24 log(1/ɛ) births is at least 16 log(1/ɛ) 1 e 6 1 ɛ. By nion bond, the probability that there are at most 72 log(1/ɛ) at the end of the dynamics is at least 1 (2ɛ + ɛ) = 1 3ɛ. Theorems 4.5 and 4.6 are smmarized in Figre 2. 19

20 5 Conclsion In this paper, we introdce and analyze a model on opinion formation. In the first part, the dynamics is deterministic and we don t hae either deaths or births of types. We show that the dynamics in this case conerges point-wise to fixed points p, where the set of actie types in p forms an independent set of G. After introdcing births and deaths of types, we show that with high probability in the long rn we reach a state in which there is no moement of poplation mass for a long period of time (aka stable ). We also show that the nmber of types is logarithmic in 1/ɛ, where ɛ is the size of a type at which it dies. A host of noel qestions arise from this model and there is mch scope for ftre work: Rate of conergence (withot births and deaths): How fast does or migration dynamics conerge point-wise to fixed points p for different choices of fnctions F? How does the strctre of G inflence the time needed for conergence? Assming F (z) = a z (linear fnctions), do the ales of α s affect the speed of conergence? Aerage case analysis: Theorem 3.7 gies qalitatie information for the behaior of the dynamics assming no births and deaths of types. Howeer, it is not clear which independent sets are more likely to occr if we start at random in the simplex. Assming F (z) = a z (linear fnctions) or F (z) = a z 3 etc (cbic), do the ales of α s affect the likelihood of the linearly stable fixed points 5? Understanding the behaior of the dynamics: From Figre 2 we see that when p lies in the regime (Θ(ɛ), Θ(1)), we don t nderstand the behaior of the system, e.g., we don t know if we hae explosion in the nmber of types (i.e., haing Θ(1/ɛ) types) in the long rn. Moreoer, we don t know if the system reaches stability (in or notion). Relaxing the notion of stability: If we relax the notion of stability so that γ-fraction of the poplation is allowed to moe (1 γ does not moe for some 0 < γ < 1), can we gie better garantees in Theorem 4.5? References [1] A. Anagnostopolos, R. Kmar, and M. Mahdian. Inflence and correlation in social networks. 14th ACM SIGKDD International Conference on Knowledge Discoery and Data Mining, pages 7 15, [2] S. Aral, L. Mchnik, and A. Sndararajan. Distingishing inflence-based contagion from homophily-drien diffsion in dynamic networks. Proceedings of the National Academy of Sciences (PNAS), [3] S. Arora, Y. Rabani, and U. Vazirani. Simlating qadratic dynamical systems is pspacecomplete (preliminary ersion). Proceedings of the Twenty-sixth Annal ACM Symposim on Theory of Compting (STOC), pages , [4] R. Axelrod. The dissemination of cltre. Jornal of Conflict Resoltion, pages , [5] E. Bakshy, I. Rosenn, C. A. Marlow, and L. A. Adamic. The role of social networks in information diffsion. 21st International World Wide Web Conference, pages , this likelihood of a fixed point is called region of attraction. 20

21 [6] A. Bhattacharyya and K. Shiragr. How friends and non-determinism affect opinion dynamics. In 54th IEEE Conference on Decision and Control (CDC), pages , [7] J. M. Cohen. Sorces of peer grop homogeneity. Sociology in Edcation, pages , [8] G. Deffant, D. Nea, F. Amblard, and G. Weisbch. Mixing beliefs among interacting agents. Jornal of Conflict Resoltion, pages 87 98, [9] M. Feldman, N. Immorlica, B. Lcier, and S. M. Weinberg. Reaching consenss ia non-bayesian asynchronos learning in social networks. In APPROX/RANDOM, pages , [10] S. Galam. Sociophysics: a reiew of Galam models. International Jornal of Modern Physics C, [11] R. Hegselmann and U. Krase. Opinion dynamics and bonded confidence: models, analysis and simlation. Jornal of Artificial Societies and Social Simlation, [12] R. A. Holley and T. M. Liggett. Ergodic theorems for weakly interacting infinite systems and the oter model. Annals of Probability, [13] M. O. Jackson. Social and economic networks. Princeton Uniersity Press, [14] D. B. Kandel. Homophily, selection, and socialization in adolescent friendships. American Jornal of Sociology, pages , [15] J. L. Kelley. General Topology. Springer, [16] D. Kempe, J. M. Kleinberg, S. Oren, and A. Slikins. Selection and inflence in cltral dynamics. In ACM Conference on Electronic Commerce (EC), pages , [17] G. E. Kreindlera and H. P. Yong. The spread of innoations in social networks. Proceedings of the National Academy of Sciences (PNAS), [18] T. LaFond and J. Neille. Randomization tests for distingishing social inflence and homophily effects. 19th International World Wide Web Conference, pages , [19] J. D. Lee, M. Simchowitz, M. I. Jordan, and B. Recht. Gradient descent only conerges to minimizers. Conference on Learning Theory (COLT), [20] R. Mehta, I. Panageas, and G. Pilioras. Natral selection as an inhibitor of genetic diersity: Mltiplicatie weights pdates algorithm and a conjectre of haploid genetics. Innoations in Theoretical Compter Science (ITCS), [21] R. Mehta, I. Panageas, G. Pilioras, P. Tetali, and V. V. Vazirani. Mtation, sexal reprodction and srial in dynamic enironments. Innoations in Theoretical Compter Science (ITCS), [22] A. Montanari and A. Saberi. The spread of innoations in social networks. Proceedings of the National Academy of Sciences (PNAS), [23] E. Mossel, J. Neeman, and O. Tamz. Majority dynamics and aggregation of information in social networks. In Atonomos Agents and Mlti-Agent Systems (AAMAS),

22 [24] I. Panageas and G. Pilioras. Aerage case performance of replicator dynamics in potential games ia compting regions of attraction. ACM Conference on Economics and Comptation (EC), [25] I. Panageas and G. Pilioras. Gradient descent only conerges to minimizers: Non-isolated critical points and inariant regions. Innoations in Theoretical Compter Science (ITCS), [26] L. Perko. Differential Eqations and Dynamical Systems. Springer, 3nd. edition, [27] K. Poole and H. Rosenthall. Patterns of congressional oting. American Jornal of Political Science, pages , [28] Y. Rabinoich, A. Sinclair, and A. Widgerson. Qadratic dynamical systems. Proc 23rd IEEE Symp Fondations of Compter Science, pages , [29] M. Shb. Global Stability of Dynamical Systems. Springer-Verlag,

23 A Appendix The following theorem holds for eery separable metric space, i.e., eery metric space that contains a contable, dense sbset. In particlar, we se this theorem for R n 1 where n is the nmber of types in the proof of Theorem 3.7. Theorem A.1 (Lindelőf s lemma [15]). For eery open coer there is a contable sbcoer. The following lemma is sed in Theorem 3.7 when we arge that the set of initial poplation ectors so that the dynamics conerges to fixed points with an nstable direction, has measre zero. It roghly states that if a fnction h is locally Lipschitz, then it preseres the measre zero sets (measre zero sets are mapped to measre zero sets). Lemma A.2 (Nll-set presering, Appendix of [20]). Let h : S R m be a locally Lipschitz fnction with S R m, then h is nll-set presering, i.e., for E S if E has measre zero then h(e) has also measre zero. Proof. Let B γ be an open ball sch that h( y) h( x) K γ y x for all x, y B γ. We consider the nion γ B γ which coer R m by the assmption that h is locally Lipschitz. By Lindelőf s lemma we hae a contable sbcoer, i.e., i=1 B i. Let E i = E B i. We will proe that h(e i ) has measre zero. Fix an ɛ > 0. Since E i E, we hae that E i has measre zero, hence we can find a contable coer of open balls C 1, C 2,... for E i, namely E i j=1 C j so that C j B i for all j and also j=1 µ(c j) < ɛ K m i. Since E i j=1 C j we get that h(e i ) j=1 h(c j), namely h(c 1 ), h(c 2 ),... coer h(e i ) and also h(c j ) h(b i ) for all j. Assming that ball C j B( x, r) (center x and radis r) then it is clear that h(c j ) B(h( x), K i r) (h maps the center x to h( x) and the radis r to K i r becase of Lipschitz assmption). Bt µ(b(h( x), K i r)) = K m i µ(b( x, r)) = K m i µ(c j ), therefore µ(h(c j )) K m i µ(c j ) and so we conclde that µ(h(e i )) µ(h(c j )) Ki m j=1 µ(c j ) < ɛ j=1 Since ɛ was arbitrary, it follows that µ(h(e i )) = 0. To finish the proof, obsere that h(e) = i=1 h(e i) therefore µ(h(e)) i=1 µ(h(e i)) = 0. 23

24 B Figres (a) The region with C corresponds to the initial poplation ectors so that the dynamics conerges to the fixed point where all the poplation is of type C. The region A+B corresponds to the initial poplation masses so that the dynamics conerges to a fixed point where part of the poplation is of type A and the rest of type B. (b) Each region A, B, C corresponds to the initial poplation ectors so that the dynamics conerges to all the poplation being of type A, B, C respectiely. It is easy to see that an initial ector (x A, x B, x C) conerges to the fixed point where all poplation is of type arg max i {A,B,C} x i. In case of ties, the limit poplation is split eqally among the tied types (symmetry). Figre 1: Migration dynamics phase portrait for path and triangle of 3 types A, B, C respectiely and for F (x) = 0.5x for all E(G). The black points and the line correspond to the fixed points. x A, x B correspond to the freqencies of people that are of type A, B. We omit x C since x C = 1 x A x B. Figre 2: In the red interal where p is O(ɛ) we hae that in the long rn the system reaches a state where there is no migration of poplation mass for a long period of time. In the green interal where p is Θ(1) we hae that there is no explosion in the nmber of types, namely the nmber is at most O(log(1/ɛ)) in the long term. In the ble interal, we don t hae a qalitatie or qantitatie characterization of the system. 24