Abstract
Complex networks impact the diffusion of ideas and innovations, the formation of opinions, and the evolution of cooperative behavior. In this context, heterogeneous structures have been shown to generate a coordinationlike dynamics that drives a population towards a monomorphic state. In contrast, homogeneous networks tend to result in a stable coexistence of multiple traits in the population. These conclusions have been reached through the analysis of networks with either very high or very low levels of degree heterogeneity. In this paper, we use methods from Evolutionary Game Theory to explore how different levels of degree heterogeneity impact the fate of cooperation in structured populations whose individuals face the Prisoner’s Dilemma. Our results suggest that in large networks a minimum level of heterogeneity is necessary for a society to become evolutionary viable. Moreover, there is an optimal range of heterogeneity levels that maximize the resilience of the society facing an increasing number of social dilemmas. Finally, as the level of degree heterogeneity increases, the evolutionary dominance of either cooperators or defectors in a society increasingly depends on the initial state of a few influential individuals. Our findings imply that neither very unequal nor very equal societies offer the best evolutionary outcome.
Introduction
It is well known that the structure of complex networks affects the outcomes of dynamical processes^{1,2}, such as epidemic outbreaks^{3,4,5}, diffusion of innovations^{6,7}, opinion formation^{8,9,10} and behavioral evolution^{11,12}. In that context, network heterogeneity is a topological feature of complex networks that is of particular importance^{13,14,15,16,17}. The network heterogeneity refers to the heterogeneity among the number of ties maintained by each individual in a network; it is often measured by the variance of the degree distribution^{18} and called degree heterogeneity. Until now most works have, however, dealt with structures that exhibit opposite levels of degree heterogeneity, such as networks where all individuals have a similar number of network ties or networks with a power law distribution of the individuals’ ties. In this paper, we explore the gap between these two extremes in the context of evolutionary game theory^{19,20} and the problem of the evolution of cooperation^{21,22}. Our goal is to understand if (1) there is a minimum threshold of heterogeneity necessary for the evolution of cooperation, and (2) if is there an optimal level of heterogeneity?
Previous theoretical works have extensively explored the role of the population structure in the evolution of cooperation. Yet, they have mainly dealt with structures that exhibit opposite levels of degree heterogeneity^{23,24,25,26,27,28,29,30,31,32,33,34}. In this regard, regular and highly symmetric structures with low levels of heterogeneity–such as lattices^{11,25,26,27}–have been shown to promote a stable coexistence dynamics between different competing behaviors, such as cooperators and defectors^{35}. This coexistence dynamics stems from the formation of compact clusters that constrain the evolution of cooperation to the boundaries between regions dominated by individuals of a single type. In contrast, structures with high levels of degree heterogeneity (e.g. scalefree networks) prompt a coordinationlike dynamics^{36,37} that result in the population being quickly driven to a monomorphic state which is dominated by one of the existing behaviors (i.e. either only cooperators or defectors). In these heterogeneous networks, the observed coordination dynamics results from the existence of higher degree nodes (i.e. hubs) that act as role models and lead to effective dynamics of competition between hubs. Thus, while both low and high heterogeneous networks promote cooperation, heterogeneous networks are able to promote cooperation in a wider range of social dilemmas.
Few works, though, have explored how a gradual interpolation between extreme levels of topological features in complex networks impact the outcome of dynamical processes. For instance, Santos et al.^{38} explored how the randomization of the connectivity patterns of a network–interpolating from a regular to a randomized homogeneous network–impacts the evolution of cooperative behavior and epidemic spreading. Following this line of research, we explore in this paper how the level of degree heterogeneity impacts the evolution of cooperation. For this purpose we analyze a wide range of networks with different levels of degree heterogeneity to show to which extent degree heterogeneity facilitates the dominance or coexistence of cooperators and/or defectors in a society.
Previous works have linked degree heterogeneity and degree variance also with social diversity, in the sense that it captures diversity among individual’s degrees of popularity and influence^{18}. Social diversity is, however, not limited to popularity and influence, nor is it only able to be modeled through to the interaction structure of the population. Indeed, recent works have explored how diversity among interactions weights^{39}, individuals’ aspirations^{40}, wealth distribution^{41,42,43}, age distribution^{44}, and the type of social dilemmas faced by individuals^{45} play an important role in dictating the fate of cooperative behavior in social systems. Central for the positive impact of diversity in promoting cooperation is the heterogeneous distribution of wealth and influence it generates^{46}. Thus our work also contributes to a better understanding of how social diversity drives the evolution of cooperation in complex networks.
To explore how different levels of degree heterogeneity impact the evolution of cooperation, we make extensive Monte Carlo simulations of the evolutionary dynamics. Moreover, we consider a wide range of populations whose structures smoothly interpolate between low and high levels of heterogeneity (see Model Section). Our findings show that many situations require a minimum threshold level of heterogeneity for cooperation to be sustainable. Moreover, there is an optimal range of heterogeneity that maximizes the population’s resilience across a wide range of social dilemmas. This is the case because an increasingly heterogeneous network structure becomes increasingly correlated with the initial state of a few highly influential individuals (see Discussion section). The findings imply that neither very equal nor extremely unequal societies, in terms of their level of degree heterogeneity, tend to offer the best evolutionary outcome.
Model
Let us consider a population with Z individuals. The structure of interaction between the Z individuals of the population is modeled by means of a complex network, where nodes represent individuals, and the links between the nodes represent social relationships between the individuals. The total number of links of an individual i corresponds to his/her degree z _{ i }. The degree distribution, D(z), describes the fraction of individuals with degree z. To generate networks with different levels of degree heterogeneity D(z), we use an algorithm of network growth and biased preferential attachment by Fortunato et al.^{47}. The algorithm proceeds as follows: it starts from a clique of three connected nodes and then sequentially Z − 3 other nodes are added; each of the newly added nodes gets attached to two preexisting ones sampled proportionally to t ^{α} (α ≥ 0.0) where t corresponds to the age of the nodes. In the limiting case of Z → ∞, the networks exhibit a degree distribution of the form D(z) ≈ z ^{−γ} with γ = (1 + α)/α ^{47,48}. Here α takes the role of the degree heterogeneity parameter exhibiting a onetoone relationship with var(z), so that a low (large) α is associated with a low (large) level of degree variance or heterogeneity.
Let us also consider that individuals behave unconditionally either as Cooperators (C) or Defectors (D). At each moment, there are k Cs and Z − k Ds in the population. The evolutionary process follows the fermi update rule ^{49,50}, thus at each time step a random individual (A) imitates the strategy of a random neighbor, B, with the probability
where β regulates the selection pressure (i.e. the level of randomness of the individuals’ decision to be a cooperator or a defector), while f _{ i } corresponds to the fitness of the individual i. When β ≪ 1, the population evolves under weak selection pressure; when β ≫ 1 the population evolves under strong selection pressure. In the limit of strong selection (large β), the results obtained under the fermi update rule become qualitatively similar to the replicator dynamics^{51}.
Individuals obtain a payoff from each interaction with their neighbors. The value of the payoff depends only on the strategies of both. The four possible outcomes can be summarized in the following payoff matrix:
where R represents the Reward payoff resulting from mutual cooperation; P is the Punishment payoff from mutual defection; S is the Sucker’s payoff obtained by a Cooperator that faced a Defector; and T represents the Temptation payoff that a Defector obtains when facing a Cooperator. Here, we consider a simplified parameter space in which R = 1, P = 0, T = λ and S = 1 − λ with 1 ≤ λ ≤ 2. Hence, λ represents both the temptation for defection and the fear of being cheated, encapsulating the harshness of the dilemma faced by individuals.
The fitness of an individual corresponds to the accumulated payoff over all interactions in which he/she participates and can be computed as
where \({n}_{i}^{C}\) (\({n}_{i}^{D}\)) is the number of Cs (Ds) in the neighborhood of i (\({n}_{i}^{C}+{n}_{i}^{D}={z}_{i}\)) and s _{ i } = 1 (0) if i is a C (D). Moreover, λ defines the strength of the dilemma faced by individuals. When λ > 0, individuals engage in the Prisoner’s Dilemma (PD) and λ captures both the temptation to defect towards a C and the fear of being cheated by a D.
In order to measure the impact of different network structures on the evolution of cooperation, we measure the level of cooperation (η _{ x0}). More details on how the level of cooperation (η _{ x0}) is computed can be found in the Methods section. We estimate η _{ x0} starting from a fraction x _{0} of cooperators and for different levels of selection pressure β, degree heterogeneity α, and harshness of the social dilemma λ. The level of cooperation (η _{ x0}) can thus be understood as the expected number of cooperators in the equilibrium. In order to explore the role that different individuals play in the evolutionary dynamics, we classify the individuals according to their network degree^{18}. Hence let us denote nodes as low degree (LDN) if their degree is bounded by 3z _{ i } ≤ m _{ α }; medium degree (MDN) when m _{ α } < 3z _{ i } ≤ 2m _{ α }; and high degree (HDN) for 3z _{ i } > 2m _{ α }, where m _{ α } corresponds to the observed maximum degree of an individual in a network with a level of degree heterogeneity of α.
Results and Discussion
We start by analyzing how the level of degree heterogeneity (α) and the selection pressure (β) impact the level of cooperation (\({\eta }_{{x}_{0}}\)) under different values of dilemma strength (λ). In the panels A, B and C of Fig. 1, three different regions can be highlighted according to their evolutionary outcome: (i) the domain of parameters strictly dominated by Ds (blue region); (ii) the domain of parameters strictly dominated by Cs (red region) and (iii) the domain of parameters where neither Cs or Ds are strictly dominant. The results in Fig. 1 suggest that cooperation requires a minimum level of degree heterogeneity (α) to be evolutionary viable – i.e. to achieve a state of C dominance with a nonzero likelihood. This result is particularly evident under weaker selection regimes (β < 1). Moreover, there is a region in the α × β domain under which the viability of cooperation is more resilient to variations in the harshness of the social dilemma faced by individuals (λ).
It must be noted that optimal levels of selection pressure have been identified in previous works^{26,27,52}. Here, though, we show that this effect is restricted to a certain interval of network heterogeneity, in our case for α > 0.6. Figure 2A shows how the level of cooperation is affected by degree heterogeneity for different levels of the dilemma strength. While a sharp transition in the levels of cooperation characterizes the evolution under small values of degree heterogeneity, increasing values of degree heterogeneity lead to a slower decay of cooperation with rising levels of dilemma strengths.
Previous works have shown that heterogeneous networks facilitate a coordinationlike dynamics at the populationwide level^{36}. Hence, next we compute the Average Gradient of Selection (Γ_{ g }(k)) in order to inspect of the nature of the populationwide dynamics with different levels of degree heterogeneity. Details on how to compute Γ_{ g }(k) can be found in the methods section. In a nutshell, negative (positive) values of Γ_{ g }(k) indicate that a population with k Cooperators and Z − k Defectors population is more likely to see an decrease (increase) in the number of cooperators. For instance, in a wellmixed population (i.e. a complete graph, where all individuals are connected with each other) and when individuals engage in the Prisoner’s Dilemma, we observe Γ_{ g }(k) to be always negative within the interval 0 < k < Z. Moreover, Γ_{ g }(k) has two trivial solutions, one at k ^{*} = 0 and another at k ^{*} = Z. Hence, we inspect the existence of internal solutions (k ^{*}) that lie inside the domain 0 < k ^{*} < Z.
Figure 2B shows the internal fixed points of the Average Gradient of Selection (Γ_{ g }(k)) at the 150^{th} generation. Gray areas represent regions in which Γ_{ g }(k) < 0 and white areas regions in which Γ_{ g }(k) > 0. Arrows point towards the direction of selection. White circles represent coordination points. These are unstable fixed points, because when the population is close to one of this points it will evolve away from such state. Conversely, black circles represent coexistence points. These are stable fixed points, because the population will converge to such a state. We see that for the interval of heterogeneity under analysis the Γ_{ g }(k) is mostly dominated by a single coordination point (x _{ L }), which implies that the population will be driven to a monomorphic state, depending on the initial fraction of cooperators (x _{0}). Hence, these results reassure us that measuring the level of cooperation (\({\eta }_{{x}_{0}}\)) equates to finding the likelihood with which a population reaches a state of full cooperation.
Figure 2C and D shows the internal fixed points of Γ_{ g }(k) during the first 150 generations when α = 0.35 and α = 1.00, respectively. Although at the beginning the Γ_{ g }(k) is mostly a negative, the evolutionary dynamics will lead to the emergence of an internal coordination point, which shapes the outcome of the population.
Next, we show the level of cooperation (\({\eta }_{{x}_{0}=0.5}\)) as a function of the temptation to defect (λ) and of the level of degree heterogeneity (α) in the case of a weak (β = 0.3) and a strong (β = 3.0) selection pressure (see Fig. 3A and B). The level of heterogeneity largely defines the domain of social dilemmas under which cooperation is evolutionary viable. Although increasing levels of degree heterogeneity diminishes the range of λ that are strictly dominated by Ds, the degree heterogeneity also does not increase the region strictly dominated by Cs. In fact, for low selection pressures there is an optimal level of heterogeneity that maximizes the range of social dilemmas strictly dominated by Cs. For β = 0.3 this is between 0.6 ≤ α ≤ 0.9. This implies that if a population structure evolves to become increasingly heterogeneous, the population will eventually reach a level under which cooperation becomes less viable. Moreover, Fig. 3A and B show that for increasing levels of degree heterogeneity there is an increase in the region in which neither strategy is strictly dominant (white area), a result that is also suggested from the analysis of the α × β domain in Fig. 1A,B and C.
We scrutinize the mechanisms that lead to the prevalence of balanced levels of cooperation (white regions) for such a wide range of parameters. Our hypothesis is that the evolutionary outcome of the population becomes correlated with the initial state of the most wellconnected individuals. With increasing levels of degree heterogeneity, the higher connected individuals (hubs) also receive more and more connections during the network growth. In consequence, the network becomes increasingly more unequal in terms of the distribution of influence and wealth. In other words, the fact that a small number of individuals in the population can participate in the majority of the available social interactions, thus accumulating a proportionally large fitness, results in these individuals dictating the outcome of the population according to their initial state.
To explore this hypothesis, we compute the phi correlation (ϕ _{M}) between the initial strategy of individuals in a degree class (LDN, MDN or HDG) and the final state of the population. This quantity is formally defined as
where n _{M}(X, Y) is the number of times we observe an individual in the degree class M with initial strategy X in a simulation where the population ended in a state dominated by strategy Y, n _{M}(X, •) = n _{M}(X, Y) + n _{M}(X, X) and n _{M}(•, X) = n _{M}(X, X) + n _{M}(Y, X). Hence, when ϕ _{M} = 0, the initial state of the nodes in degree class M are not correlated with the final state of the population, however when ϕ _{M} = 1, the population always ends up in a state dominated by the initial strategy of the nodes of degree class M. Conversely, when ϕ _{M} = −1, the population always ends up in the opposite state of the initial strategy of the nodes of degree class M.
Figure 3C to H show the correlation between the initial state of individuals from different degree classes and the final state of the population. We draw three main conclusions from these results: (i) in the region dominated by a single strategy (blue and red areas in Fig. 4A and B) the final state is not predicted by the elements of any of the degree classes; (ii) Medium degree (MDN) and Low degree (LDN) nodes are weak predictors of the final level of cooperation; and (iii) the most connected, high degree, individuals (HDN) play a fundamental role in the region associated with a balanced levels of cooperation. In this region, we observe a strong correlation between their initial state and the fate of the population, meaning that the evolutionary outcome of the society depends on the probability (set by the initial fraction of cooperators parameter) of the highly connected individuals starting as cooperators or defectors. This confirms the hypothesis that the evolutionary outcome is driven by the initial state of the most influential individuals in the population. Moreover, HDN nodes play a central role in determining the final level of cooperation of the population in domains associated with the phase transition between the two strictly dominance regimes. Indeed, at this point the evolution of cooperation is dictated by the initial state of the most influential and wealthiest individuals.
However, the initial state of the most influential individuals (HDN) does not suffice in explaining what determines the final evolutionary outcome of the society. Without the existence of other Cs in the population, the odds of cooperation being sustainable are very slim. Figure 4 shows how the evolutionary outcome (\({\eta }_{{x}_{0}}\)) depends on the initial abundance of cooperators (x _{0} in the Yaxis) under different heterogeneity levels (α in the Xaxis) and the temptation to defect (λ). Again, it is clear that the evolution of cooperation requires a minimum level of degree heterogeneity and, moreover, a minimum number of cooperators to be evolutionary viable. For increasingly difficult social dilemmas, i.e. larger λ, cooperation viability requires a larger level of degree heterogeneity, hence a more unequal population.
Figure 5 shows the correlation between the initial state of the most influential individuals (HDN) and the final state of the population under similar conditions as considered in Fig. 4. The obtained results suggest, once again, a strong correlation in the domain characterized by a balanced levels of cooperation (nonblue or nonred areas in Fig. 4).
Conclusions
In this work we have explored how different levels of degree heterogeneity impact the evolution of cooperation in the Prisoner’s Dilemma. To that end, we made use of an algorithm of network growth and biased preferential attachment to generate a large set of network structures that interpolate between low and high levels of heterogeneity.
We find evidence of a threshold level of degree heterogeneity above which there is a substantial increase in the range of dilemmas where cooperation is evolutionary dominant. More importantly, we show that increasing levels of degree heterogeneity lead to an evolutionary trade off for populations: although the range of dilemmas in which cooperation can evolve increases, the range of dilemmas under which cooperation is strictly dominant decreases. Indeed, we observe that the evolution of cooperation becomes coupled with the initial state of the most well connected individuals. Hence, the ability of populations to overcome the underlying social dilemmas individuals engage becomes tied to the ability of a few individuals do so. Our findings imply that neither equal nor very unequal societies (in terms of degree heterogeneity) tend to offer the best evolutionary outcome, but instead suggest that there is an optimal level of equality. Although here we have explored static structures, recent works have shown how coevolution of population structure often leads to an increase in degree heterogeneity^{53,54,55,56}. Whenever the degree heterogeneity of a population increases, it also increases the dependence of the population on the decisions of a handful of influential individuals.
Open questions, though, are whether competition between societies can explain the preferential selection of heterogeneous social structures and whether a particular level of heterogeneity is favored under such a competitive environment. Finally, we explored the role of degree heterogeneity in the context of the evolution of cooperation. Our results suggest that there is a need to pay attention to the role that different levels of network heterogeneity play in the context of opinion dynamics and epidemic outbreaks.
Methods
Level of cooperation
(\({\eta }_{{x}_{0}}\)) is estimated by averaging the number of cooperators in the population after a large transient period of up to 10^{5} generations from 10^{4} independent simulations. Each simulation lasts for 5 × 10^{3} generations and starts from a fraction x _{0} of Cs that are randomly distributed across the network. Given the coordination nature of the populationwide dynamics in heterogeneous networks^{36}, \({\eta }_{{x}_{0}}\) provides a good approximation of the likelihood of a population reaching a monomorphic state dominated by Cs. Along the manuscript we compute under \({\eta }_{{x}_{0}}\) different scenarios of degree heterogeneity (α), selection pressure (β) and the temptation to defect (λ).
Average Gradient of Selection
(AGoS) captures the populationwide dynamics on structured populations. This quantity is the numerical counterpart of the drift term in a BirthDeath stochastic processes on finite and wellmixed populations^{57}. We can estimate this quantity conveniently by computing the difference between the probability to increase (\({\xi }_{g}^{+}(k)\)) and the probability to decrease (\({\xi }_{g}^{}(k)\)) the number of Cooperators by one when the population is in a state with k Cooperators and Z − k Defectors. The AGoS is, thus, formally defined as
where \({\xi }_{g}^{\pm }(k)\) is numerically computed according to
where Λ_{ g }(k) accounts for the total number of times the population was observed in a state with k Cooperators at generation g over all Ω simulations and Θ(a,b) is a square discrete function that is equal to 1 if b − 1 ≤ a < b, being 0 otherwise. Finally \({\xi }_{\omega }^{\pm }(k,t)\) is the transition probability at time t of timeseries ω, that is
where ζ _{ j } is the set of closest neighbours of j and k _{ t } is the number of Cs at time t. In this expression, the first sum is calculated over each individual in the population and the second sum is calculated over each individual’s neighbourhood. In order to estimate Γ(k) we let a population evolve for 150 generations. Each simulation of the evolutionary process starts from a possible random state (i.e. starting from a random number of Cooperators k randomly selected from the interval 1 < k < Z). We repeat this for a total of Ω = 2.5 × 10^{7} times. For each iteration we estimate the \({\xi }_{\omega }^{\pm }(k,t)\) as described above.
References
 1.
Nowak, M. A. & May, R. M. Evolutionary games and spatial chaos. Nature 359, 826–829 (1992).
 2.
Barrat, A., Barthelemy, M. & Vespignani, A. Dynamical processes on complex networks (Cambridge university press, 2008).
 3.
PastorSatorras, R. & Vespignani, A. Epidemic spreading in scalefree networks. Physical Review Letters 86, 3200 (2001).
 4.
Moreno, Y., PastorSatorras, R. & Vespignani, A. Epidemic outbreaks in complex heterogeneous networks. The European Physical Journal BCondensed Matter and Complex Systems 26, 521–529 (2002).
 5.
Newman, M. E. Spread of epidemic disease on networks. Physical review E 66, 016128 (2002).
 6.
Kocsis, G. & Kun, F. The effect of network topologies on the spreading of technological developments. Journal of Statistical Mechanics: Theory and Experiment 2008, P10014 (2008).
 7.
Castellano, C., Fortunato, S. & Loreto, V. Statistical physics of social dynamics. Reviews of Modern Physics 81, 591 (2009).
 8.
Fortunato, S. Damage spreading and opinion dynamics on scalefree networks. Physica A: Statistical Mechanics and its Applications 348, 683–690 (2005).
 9.
Acemoglu, D. & Ozdaglar, A. Opinion dynamics and learning in social networks. Dynamic Games and Applications 1, 3–49 (2011).
 10.
Sood, V. & Redner, S. Voter model on heterogeneous graphs. Physical Review Letters 94, 178701 (2005).
 11.
Szabó, G. & Fath, G. Evolutionary games on graphs. Physics Reports 446, 97–216 (2007).
 12.
Rand, D. G., Nowak, M. A., Fowler, J. H. & Christakis, N. A. Static network structure can stabilize human cooperation. Proceedings of the National Academy of Sciences 111, 17093–17098 (2014).
 13.
Liljeros, F., Edling, C. R., Amaral, L. A. N., Stanley, H. E. & Åberg, Y. The web of human sexual contacts. Nature 411, 907–908 (2001).
 14.
Dorogovtsev, S. N. & Mendes, J. F. Evolution of networks: From biological nets to the Internet and WWW (OUP Oxford, 2013).
 15.
Csányi, G. & Szendröi, B. Structure of a large social network. Physical Review E 69, 036131 (2004).
 16.
Mislove, A., Marcon, M., Gummadi, K. P., Druschel, P. & Bhattacharjee, B. Measurement and analysis of online social networks. In Proceedings of the 7th ACM SIGCOMM conference on Internet measurement, 29–42 (ACM, 2007).
 17.
Gjoka, M., Kurant, M., Butts, C. T. & Markopoulou, A. Walking in facebook: A case study of unbiased sampling of osns. In Infocom, 2010 Proceedings IEEE, 1–9 (IEEE, 2010).
 18.
Santos, F. C., Pinheiro, F. L., Lenaerts, T. & Pacheco, J. M. The role of diversity in the evolution of cooperation. Journal of Theoretical Biology 299, 88–96 (2012).
 19.
Perc, M. et al. Statistical physics of human cooperation. Physics Reports (2017).
 20.
Nowak, M. A. Evolutionary Dynamics (Harvard University Press, 2006).
 21.
Hofbauer, J. & Sigmund, K. Evolutionary games and population dynamics (Cambridge university press, 1998).
 22.
Nowak, M. A. Evolving cooperation. Journal of Theoretical Biology 299, 1–8 (2012).
 23.
Nakamaru, M., Matsuda, H. & Iwasa, Y. The evolution of cooperation in a latticestructured population. Journal of Theoretical Biology 184, 65–81 (1997).
 24.
Szolnoki, A. & Szabó, G. Cooperation enhanced by inhomogeneous activity of teaching for evolutionary prisoner’s dilemma games. Europhysics Letters 77, 30004 (2007).
 25.
Szabó, G. & Töke, C. Evolutionary prisoner’s dilemma game on a square lattice. Physical Review E 58, 69 (1998).
 26.
Szabó, G., Vukov, J. & Szolnoki, A. Phase diagrams for an evolutionary prisoner’s dilemma game on twodimensional lattices. Physical Review E 72, 047107 (2005).
 27.
Vukov, J., Szabó, G. & Szolnoki, A. Cooperation in the noisy case: Prisoner’s dilemma game on two types of regular random graphs. Physical Review E 73, 067103 (2006).
 28.
Santos, F. C. & Pacheco, J. M. Scalefree networks provide a unifying framework for the emergence of cooperation. Physical Review Letters 95, 098104 (2005).
 29.
Santos, F. C., Santos, M. D. & Pacheco, J. M. Social diversity promotes the emergence of cooperation in public goods games. Nature 454, 213–216 (2008).
 30.
Assenza, S., GómezGardeñes, J. & Latora, V. Enhancement of cooperation in highly clustered scalefree networks. Physical Review E 78, 017101 (2008).
 31.
Poncela, J., GómezGardenes, J., Flora, L. & Moreno, Y. Robustness of cooperation in the evolutionary prisoner’s dilemma on complex networks. New Journal of Physics 9, 184 (2007).
 32.
Poncela, J., GómezGardenes, J. & Moreno, Y. Cooperation in scalefree networks with limited associative capacities. Physical Review E 83, 057101 (2011).
 33.
Xia, C.Y., Meng, X.K. & Wang, Z. Heterogeneous coupling between interdependent lattices promotes the cooperation in the prisoner’s dilemma game. PLoS ONE 10, e0129542 (2015).
 34.
Jin, Q., Wang, L., Xia, C.Y. & Wang, Z. Spontaneous symmetry breaking in interdependent networked game. Scientific reports 4 (2014).
 35.
Hauert, C. & Doebeli, M. Spatial structure often inhibits the evolution of cooperation in the snowdrift game. Nature 428, 643–646 (2004).
 36.
Pinheiro, F. L., Pacheco, J. M. & Santos, F. C. From local to global dilemmas in social networks. PLoS ONE 7, e32114 (2012).
 37.
Skyrms, B. The stag hunt and the evolution of social structure (Cambridge University Press, 2004).
 38.
Santos, F. C., Rodrigues, J. & Pacheco, J. M. Epidemic spreading and cooperation dynamics on homogeneous smallworld networks. Physical Review E 72, 056128 (2005).
 39.
Ma, Z.Q. et al. Heterogeneous link weight promotes the cooperation in spatial prisoner’s dilemma. International Journal of Modern physics C 22, 1257–1268 (2011).
 40.
Perc, M. & Wang, Z. Heterogeneous aspirations promote cooperation in the prisoner’s dilemma game. PLoS ONE 5, e15117 (2010).
 41.
Wang, J., Fu, F. & Wang, L. Effects of heterogeneous wealth distribution on public cooperation with collective risk. Physical Review E 82, 016102 (2010).
 42.
Xu, B. & Lan, Y. The distribution of wealth and the effect of extortion in structured populations. Chaos, Solitons & Fractals 87, 276–280 (2016).
 43.
Hartmann, D. Economic complexity and human development: how economic diversification and social networks affect human agency and welfare, vol. 110 (Routledge, 2014).
 44.
Wang, Z., Zhu, X. & Arenzon, J. J. Cooperation and age structure in spatial games. Physical Review E 85, 011149 (2012).
 45.
Amaral, M. A., Wardil, L., Perc, M. & da Silva, J. K. Evolutionary mixed games in structured populations: Cooperation and the benefits of heterogeneity. Physical Review E 93, 042304 (2016).
 46.
Perc, M. & Szolnoki, A. Social diversity and promotion of cooperation in the spatial prisoner’s dilemma game. Physical Review E 77, 011904 (2008).
 47.
Fortunato, S., Flammini, A. & Menczer, F. Scalefree network growth by ranking. Physical Review Letters 96, 218701 (2006).
 48.
Goh, K.I., Kahng, B. & Kim, D. Universal behavior of load distribution in scalefree networks. Physical Review Letters 87, 278701 (2001).
 49.
Traulsen, A., Nowak, M. A. & Pacheco, J. M. Stochastic dynamics of invasion and fixation. Physical Review E 74, 011909 (2006).
 50.
Traulsen, A., Pacheco, J. M. & Nowak, M. A. Pairwise comparison and selection temperature in evolutionary game dynamics. Journal of theoretical biology 246, 522–529 (2007).
 51.
Xia, C.Y., Meloni, S., Perc, M. & Moreno, Y. Dynamic instability of cooperation due to diverse activity patterns in evolutionary social dilemmas. EPL (Europhysics letters) 109, 58002 (2015).
 52.
Pinheiro, F. L., Santos, F. C. & Pacheco, J. M. How selection pressure changes the nature of social dilemmas in structured populations. New Journal of Physics 14, 073035 (2012).
 53.
Santos, F. C., Pacheco, J. M. & Lenaerts, T. Cooperation prevails when individuals adjust their social ties. PLoS Comput Biol 2, e140 (2006).
 54.
Pinheiro, F. L., Santos, F. C. & Pacheco, J. M. Linking individual and collective behavior in adaptive social networks. Physical Review Letters 116, 128702 (2016).
 55.
Poncela, J., GómezGardenes, J., Flora, L. M., Sánchez, A. & Moreno, Y. Complex cooperative networks from evolutionary preferential attachment. PLoS ONE 3, e2449 (2008).
 56.
Poncela, J., GómezGardeñes, J., Traulsen, A. & Moreno, Y. Evolutionary game dynamics in a growing structured population. New Journal of Physics 11, 083031 (2009).
 57.
Traulsen, A., Claussen, J. C. & Hauert, C. Coevolutionary dynamics in large, but finite populations. Physical Review E 74, 011901 (2006).
Acknowledgements
This work was funded by the Cooperative Agreement between the Masdar Institute of Science and Technology and the MIT Media Lab Consortia. The authors would like to thank Pedro Neves, Fernando P. Santos, Vítor V. Vasconcelos, Cristian CandiaCastroVallejos, Francisco C. Santos and Jorge M. Pacheco for helpful discussions and insights and the two anonymous referees for the helpful feedback.
Author information
Affiliations
Contributions
F.L.P. and D.H. conceived the experiment(s), conducted the experiment(s) and analyzed the results.
Corresponding authors
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Pinheiro, F.L., Hartmann, D. Intermediate Levels of Network Heterogeneity Provide the Best Evolutionary Outcomes. Sci Rep 7, 15242 (2017). https://doi.org/10.1038/s41598017155557
Received:
Accepted:
Published:
Further reading

Pancancer mapping of differential proteinprotein interactions
Scientific Reports (2020)

The network effects of NGOs on social capital and innovation among smallholder farmers: a case study in Peru
The Annals of Regional Science (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.