Cooperation has been vital to the evolution of all living things, including single-celled organisms (Velicer, 2005, 2003; Velicer & Stredwick, 2002; Crespi, 2001; Velicer et al., 2000; Boorman & Levitt, 1980), fish (Brosnan et al., 2003; Dugatkin, 1991, 1992, 1997; Milinski, 1987), birds (Brown & Brown, 1996; Faaborg et al., 1995), canines (Creel & Creel, 2002; Courchamp & Macdonald, 2001; Fentress & Ryon, 1986), felines (Caro, 1994; Packer & Pusey, 1982), non-human primates (Brosnan & de Waal, 2003; de Waal, 1996, 1982; Harcourt & de Waal, 1992; Chapais, 1992), and humans (Ostrom et al., 1999; Fehr & Fischbacher, 2003; Johnson et al., 2003).

Even so, the evolution of cooperative, pro-social behavior under circumstances in which individual interests are at odds with common interests, (circumstances characterized as social dilemmas (Gotts et al., 2003; Dawes & Messick, 2000)), remains a largely unsolved, multidisciplinary puzzle (Hammerstein, 2003). Approaches to these types of problems have, for the most part, been applications of evolutionary game theory (Gintis, 2000; Hofbauer & Sigmund, 1998; Maynard-Smith, 1982; Maynard-Smith & Price, 1973; Trivers, 1971; Hamilton, 1967; von Neumann & Morgenstern, 1944) and due to their importance as generalized models of many important socio-economic situations (Tomassini, 2006), iconic games such as the Prisoner’s Dilemma have been widely employed as metaphors (Doebeli & Hauert, 2005; Axelrod & Hamilton, 1981; Axelrod, 1984; Nowak & Sigmund, 1992, 2004; Nowak & May, 1992; Maynard-Smith, 1982; Sugden, 1986).

At the same time, the study of networks, complex systems, and nonlinear dynamics has pervaded all of science (Strogatz, 2001)1. Indeed, E.O. Wilson, who once characterized the evolution of cooperation as one of the greatest challenges for modern biology (Wilson, 2000), more recently made a more emphatic appeal for research on complex systems. “The greatest challenge today, not just in cell biology and ecology, but in all of science, is the accurate and complete description of complex systems. Scientists have broken down many kinds of systems. They think they know most of the elements and forces.

Fig. 1: Figure 1

A complex network in which each agent has N — 1 possible connections, a full range of degree and heterogeneity of degree.

The next task is to reassemble them, at least in mathematical models that capture the key properties of the entire ensembles.” (Wilson, 1998: 85). The application of complex systems tools and network analysis methodologies to the study of social dilemmas represents a very new, but extremely promising means of shedding light on the quandary of cooperation (Hanaki et al., 2007; Fu et al., 2007; Ohtsuki et al., 2006; Santos et al., 2006a, 2006b; Szabo & Fath, 2006; Tomassini et al., 2006; Santos & Pacheco, 2005; Santos et al., 2005; Lieberman et al., 2005; Abramson & Kuperman, 2001; Skyrms & Pemantle, 2000).

In this work, we situate players engaging in social dilemma games on complex social networks, thereby avoiding artificial limitations placed on the number of partners each player may have (a measure known as “degree”) by the more structured, lattice-like architectures common to the earlier literature on such games. These more generalized network architectures offer a full range of degree (from 0 to N — 1), which allows us to more fully investigate the impact of average degree on the evolution of pro-social behavior. In addition, large numbers of simulations subjected to statistical analysis help to address the primary concerns of this research: the effect of degree variance, what we call heterogeneity of degree, on cooperative behavior and the extent to which this effect can be separated from other explanatory factors. We begin with detailed descriptions of the models under study. Next, we describe results of simulations conducted in this framework, systematically varying certain aspects of network architecture and measuring the effect on the evolution of cooperative behavior. Finally, we discuss the implications of these findings for conflict management and briefly describe two examples where they may find practical application.

Model Descriptions

We examine a population of N players each engaging in a repeated prisoner’s dilemma game with a neighborhood of other players defined by particular network architectures. See Figure Two. The set of players with whom player i interacts in period t is denoted by Ωi,t. In each generation, which is comprised of g games, each player accumulates an adaptive score based upon a standard payoff matrix described in more detail below. At the end of each generation, each player observes the payoffs and strategies of each neighbor and stochastically updates2 their strategy with probability ρ ∈ [0,1] by imitating

Fig. 2: Figure 2

Schematic illustration of various network architectures, all with 25 nodes, roughly in ascending order of heterogeneity. (a) Fully connected network. (b) Ring lattice with all nodes connected to it neighbors out to some range k (here k = 3). (c) Small world network starting with ring lattice and adding shortcut links between random pairs of nodes. (d) Random network constructed with connection probability, p = .15. (e) Scale-free network constructed by attaching nodes at random to previously existing nodes, where the probability of attachment is proportional to the degree of the target node, i.e., “the rich get richer.”

the strategy of the neighbor with the highest adaptive score (including themselves). Ties in high scores are broken at random.

Strategic Dynamics

For each period t, players choose to either cooperate C or defect D with each of its neighbors Ωi,t and the strategic decision for each neighbor is independent of the decisions with regard to other neighbors, that is, a player can choose to cooperate with some neighbors and defect with others (c.f. Hanaki, 2007). Each neighbor j ∈ Ωi,t faces a symmetrical decision giving rise to the to a standard payoff matrix

i, j C D
C R, R S, T
D T, S P,P
Where π(si ,sj) is the payoff for player i choosing strategy si when neighbor j chooses strategy sj ,

π(C,C) = R, π(C,D) = S, π(D,C) = T

and π(D,D) = P.

In keeping with the standard structure of a social dilemma, T > R > P > S, which makes defection a dominant strategy, that is, defection results in a higher payoff as compared to cooperation regardless of what strategy the opponent neighbor chooses, and 2R > (T + S), which insures that mutual cooperation is preferred over all other strategy sets in the sense that it produces maximum aggregate outcomes. The unique equilibrium for the game, mutual defection, thus leads to a Pareto-suboptimal solution.

For each generation, each player accumulates an adaptive score for g games for all neighbors. Following the logic that the maintenance of networks with more neighbors would involve more cost than networks with fewer neighbors, we reduce adaptive scores by θ(k), the total cost of interaction with a network of k neighbors. Thus, the net payoff for each player i accumulated in a time period t is:

where θ(k) is an increasing function of k with the specific form θ(k) = ckα , where α≥1 and 0 ≤ cP (Hanaki, 2007).

Imitation Dynamic

After each generation, each player examines the accumulated adaptive scores of each of its neighbors, Ωi,t, and its own accumulated adaptive score, and either adopts by imitation the strategy of the most successful neighbor, or keeps it own strategy if it has been most successful, to be employed in the next generation, formally

Fig. 3: Table 1

Model Variable Specification

If more than one player in the neighborhood shares the highest accumulated adaptive score, ties are broken at random.

Network measures

For each run of the simulation, which is comprised of a large number of generations sufficient to arrive at equilibrium in the strategy population, a number of variables are recorded: population, average degree, heterogeneity of degree, network architecture, and cooperation. Network architecture is recorded as lattice, small world, random, or scale-free (fully connected is a special case of lattice). Cooperation is measured as a ratio of player decisions to cooperate to the total number of cooperation/defection decisions. See Table 1.

Simulation Results

We ran the simulation as described above 1,000 times creating stochastic networks by drawing network architecture uniformly from lattice, small world, random, or scale-free; drawing population uniformly from a range of 10 to 100; and drawing average degree uniformly from a range of 2 to 10. Heterogeneity of degree ranged from 0 to 4 as a function largely of network architecture. Each run was for 1,000 generations with the cooperation ratio measured in the last 100.

First, in Model 1, we regressed cooperation on population, average degree, heterogeneity of degree and three indicator variables representing four network architectures, lattice (as the base case), small world, random and scale free. We included the indicator variables to capture any variation resulting from network architectural differences not captured by the other independent variables.

We hypothesized that population size would have a positive effect on cooperation, and that both average degree and heterogeneity of degree would have a negative effect. While Model 1 bore out the first two hypotheses (population coefficient = .245, p <= .000 and average degree coefficient = -8.170, p <= .000), heterogeneity of degree showed a significant positive effect (heterogeneity of degree coefficient = 31.847, p <= .000). However, collinearity diagnostics indicated that the two most heterogeneous network architectures were highly collinear with the heterogeneity

Fig. 4: Model 1

Cooperation ( “Coop”) regressed on Population ( “Pop”), Average Degree ( “Degree”), Heterogeneity of Degree ( “Hetero”) and three indicator variables representing four network architectures, Lattice (base case), Small World ( “SWDum”), Random ( “RNDum”), and Scale Free ( “SFDum”).

Fig. 5: Model 2

Cooperation ( “Coop”) regressed on Population ( “Pop”), Average Degree ( “Degree”), and Heterogeneity of Degree ( “Hetero”)

of degree variable (variance proportions on dimension 7: random network = .91, scale-free network = .94, and heterogeneity of degree = .91). This collinearity made model coefficients uninterpretable.3 Further, the small-world network indicator variable failed to achieve statistical significance (p = .660). These results offered confidence that the indicator variables were not adding significant additional explanatory power.

Subsequently, in Model 2, we removed the indicator variables, and regressed cooperation on population, average degree, and heterogeneity of degree. In this more parsimonious model, population had a significant positive effect (population coefficient = .428, p <= .000), average degree had a significant negative effect (average degree coefficient = -4.78, p <= .000), and heterogeneity of degree had a significant negative effect (heterogeneity of degree coefficient = -4.208, p <= .001). Additionally, collinearity diagnostics showed that each of the independent variables was loading highly on its own dimension (population variance proportion on dimension 4 = .75, average degree variance proportion on dimension 3 = .86, heterogeneity of degree variance proportion on dimension 2 = .93). Based on this evidence, we concluded that heterogeneity of degree has a significant negative effect on the evolution of cooperation and that this effect is independent of the negative effect of average degree.


Generalized models of social dilemmas have been employed to increase understanding related to many important socio-economic circumstances (Tomassini, 2006), providing insights into the importance of social architecture that may even dominate individual preference (Schelling, 1978) and the ability of local interaction or “correlated association” to drive cooperative outcomes contrary to those predicted by Nash equilibria (Skyrms, 1996, 2003). These insights have been shown to have a high level of universality, with significance for numerous human social contexts as well as many non-human social species.

Our work places these social dilemma dynamics on complex social networks and confirms the importance of local interaction. As average degree increases and the relative importance of correlated association is decreased, our models demonstrate that pro-social behavior is less likely to evolve. Our models also offer a new insight. As the variance associated with degree within the network, or the heterogeneity of degree increases, the likelihood of pro-social strategy evolution is again reduced, and this negative effect is shown to be independent of the negative effect of average degree.

In Bowling Alone, Robert Putnam worries that the decline of social capital that he sees in the declining memberships in civic organizations may undermine the civil engagement that according to him is necessary for a strong democracy (Putnam, 2001). The results of our study suggest that the problem may be more nuanced. It may not be, in fact, the mere magnitude of social connections, but the nature of these connections that should concern us most.4 Merely promoting the development of dense social networks may lead us down a path to social decline. More important may be the design of institutions that promote homogeneity in social connectedness — increasing homogeneity of degree produces network effects that make the emergence of pro-social behaviors more likely thereby increasing overall social welfare. As such, homogeneity of degree is properly thought of as a collective good.

This prescription offers those involved in conflict management a deceptively simple means of promoting the evolution of cooperation within any given group — by increasing homogeneity in social connectedness. Astute readers may point out, with appropriate concern, that such homogeneity can be obtained by suppressing highly connected individuals as well as taking steps to increase connectedness among relative social isolates. To this objection we would respond by suggesting that cooperation is certainly not the only goal for would-be social engineers and it is, in our view, not the most important goal. Freedom of association may trump cooperation in these considerations and counsel against attempting to restrict the density of any individual’s social capital. On the other hand, helping to integrate social isolates into beneficial social networks would, if we are right, benefit us all. The European Union has made more progress along these lines. Where as in the United States, we largely treat poverty as a welfare problem and provide financial resources, the EU views poverty as a lack of access to social capital and offers targeted transportation, Internet access, and job training and education — all directed at social inclusion.

One other example, that offers a living laboratory of sorts where social architecture could be quite easily manipulated, are the plethora of on-line social network communities that have sprung up on the Internet in recent years. In these virtual communities, membership categories and differential incentive structures could allow operators to experimentally influence the structure of networks and sub-networks and facilitate the capture of metrics for cooperation. Such experiments would not only help to refine interventions that conflict managers may employ in the “real” world, but would increasingly find direct application as the “real” world and the virtual world of social interaction become less and less distinguishable.


  1. For example, these ideas have recently been employed to study the architecture of the Internet (Faloutsos et al., 1999), the topology of food webs (Williams & Martinez, 2000), and the metabolic network of the bacterium Escherichia coli (Jeong et al., 2000).

  2. Note that where ρ<1, updates are asynchronous (Huberman & Glance, 1993).

  3. This is not to say that Model 1 as a whole is unreliable. Highly correlated independent variables can bring about dramatic, unpredictable changes in coefficient estimates, sometimes even changing the direction, or sign. Additionally, in this particular case, the collinearity of two of the indicator variables with the heterogeneity variable and the lack of significance related to the third indicator variable suggests that these indicator variables do not explain enough additional variance to warrant their inclusion.

  4. Putnam’s work also suggests that local cohesiveness, or “clumpiness” may have a determinative effect—this is a network measure not included here, but planned for future studies.