# Evolution of Cooperation with Moore Neighborhood and Self-Playing Rule

Authors: Yanbing Yang; Junhu Ruan; Bin Liu; Yi Liu; Yan Shi
Journal: Discrete Dynamics in Nature and Society
 (2015)

---

## Abstract

Evolutionary spatial game is a promising way to unravel the mystery of cooperation, and it has been well recognized that spatial structure could enable cooperation to persist. Schweitzer et al.’s lattice model provides an innovative method to solve the problem. In this paper, we conduct simulations using the same von Neumann neighborhood as in Schweitzer et al.’s study (2002) and observe the effect of initial population and lattice size on the evolution of cooperation. Then, we extend the model with a more complicated Moore neighborhood and self-playing rule for each central player. Simulation results not only provide new evidence for the persistence of cooperation in the evolution with spatial structures, but also investigate critical conditions for the spatial coexistence or the invasion of cooperators and defectors with the more complicated neighborhood.


---

## Body

## 1. Introduction

Cooperation is one crucial driving force for human society to burgeon and flourish. Actually, it is a fundamental factor that impels nature to evolve unceasingly into the amazing world. As a ubiquitous phenomenon, cooperation can be found everywhere, from foraging ant colonies to workers in assembly lines, from evolving chromosomes in evolution of human to communicating nodes in ad hoc networks [1–3]. However, why can the cooperation be sustained and evolved?To research into the evolution of cooperation, evolutionary game is a powerful framework [4, 5]. Since evolutionary stable strategy (ESS) was first proposed by Smith and Price [6], it has drawn a lot of attention from biologists, economists, and game theorists. Its striking advantages originate from the expansion of classic game theory in aspects of both bounded rationality and dynamics. As Axelrod remarked [7], “the evolutionary approach is based on a simple principle: whatever is successful is likely to appear more often in the future”; thus successful strategies spread widely and less successful ones are doomed to be wiped out by virtue of nature selection.Prisoner’s dilemma is a most famous paradigm to dig into the emergence and persistence of cooperation [8]. As a two-person nonzero sum game, it reflects profound conflict between collective rationality and individual rationality. Everyone in this game has two pure strategies: cooperate and defect (resp., denoted by C and D), and they have to make decisions simultaneously. The source generating dilemma situation is that either prisoner can achieve maximum payoff through adopting D. In evolutionary contexts where population is well-mixed, defection is the only ESS. However, the iterated prisoner’s dilemma introduced by Axelrod casts a light on this issue [7]. Axelrod conducted computer tournaments and identified that tit for tat (TFT) strategy, which simply cooperates at the first round and then takes whatever the opponent did in his last round, surprisingly performed best. The force behind this strategy that drives cooperation to evolve is kindness, forgiving, and retaliatory [2]. In fact, TFT is an ESS when the shadow of the future is large enough and can invade the population of defection through forming clusters.In order to take a step towards realistic situations, spatial structure is further introduced into evolutionary prisoner’s dilemma [9]. Axelrod seems to be the first one modeling this [2]. He considered that a territory is subdivided into square cells, with each cell occupied by one player. Everyone only plays games with four nearest neighbors and updates his strategy by imitating the strategy used by the most successful neighbor in the last round. After simulations, he drew a conclusion that strategies in the spatial structure perform at least as good as mean-field structure, and he also revealed the process that one TFT invader spreads to the most population composed by the ALLD (all-defect) strategy.In this paper, we focus on how to demonstrate various simulation results with respect to dynamics of lattice models. Schweitzer et al. [10] put forward an innovative way to comprehend the spatial effect. In Schweitzer et al.’s work [10], the evolution is based on the von Neumann neighborhood where each player has four nearest neighbors. However, different neighborhood structures may lead to different evolutionary patterns of the invasion or the coexistence of cooperator and defectors. Meanwhile, Schweitzer et al. [10] assumed each player could not play with itself. The assumption means each player cannot learn from itself, which is obviously inconsistent with the natural and especially social evolution of cooperation. Motivated by these observations, in this paper we first conduct simulations as Schweitzer et al. did to make further analysis on the evolution of cooperation with the von Neumann neighborhood and more importantly extend the model to a complicated level where the neighborhood turns into the structure of Moore and each cell can play games with itself. Based on our simulations, we analyze the effect of initial population and system size on the evolution of cooperation and investigate critical conditions for the spatial coexistence or the invasion of players with the Moore neighborhood and self-playing rule.In Section2, we give a brief account of current works on evolutionary spatial game. In Section 3, the evolution models with the von Neumann and Moore neighborhoods are, respectively, presented. In Section 4, we conduct simulations with different models and conditions to investigate the corresponding evolution patterns. Conclusions are drawn in Section 5, with recommendations on future studies.


## 2. Literature Review

Spatial structure can enable the persistence of cooperation in the evolutionary prisoner’s dilemma game, which is a very positive and promising result and first derived by Nowak and May [11, 12]. Like Axelrod’s model, players are located in a two-dimensional lattice, and everyone merely communicates with neighbors surrounding them. To observe the spatial effect, they pulled out the repeated effect by assuming that every player only has two choices in one generation: C and D, and they also adopted deterministic update rule: imitate-the-best. By the succinct model, Nowak and May found that C and D persist together indefinitely to reach variety of spatial and temporal dynamics dependent on the value of a polymorphic parameter b which includes stable, quasistable, dynamic fractal, and chaotically varying states. A decade later, this theoretical prediction was confirmed by biological experiments [13]. Huberman and Glance doubted the above model’s usefulness for real-world systems [14] and proposed an asynchronous updating system where players can be updated at any time. Besides the lattice-based models, more and more attention has been paid on network-based evolutionary games. Detailed advances regarding this could refer to recent works by Szabó and Fáth [15], Poncela et al. [16], Perc [17], and Vilone et al. [18].Stochasticity is very important for finite population, so the scores obtained by the cells within the neighborhood do not fully determine which player will take over in the next generation but only specify the probability that a player will take over. Killingback and Doebeli extended their research object by studying what evolutionary dynamics Hawk-Dove game can bring. As a result, it has higher level of cooperation and the spatiotemporal chaos is quite different from prisoner’s dilemma [9, 19]. Ellison concentrated on the coordination games in spatial structure and showed that they quickly converge to the risk-dominant equilibrium [20]. Szabó et al. studied the spatial evolutionary games with myopic players whose payoff interest is tuned from selfishness to other-regarding preference via fraternity [21, 22].Recently, more social behaviors are considered in the evolutionary dilemma. Wang et al. [5] argued that the strategy choice of players should be affected by all their neighbors rather than the most successful neighbor and proposed an adaptive strategy-adoption rule to imitate the real-world social influence. Du and Fu [23] thought noise may be present in individual strategy learning and compared the impact of noise on spatial organization of cooperation. Gruji et al. [24] carried out comparative experiments and found that moody conditional cooperation is more supported than noninnovative game dynamics such as imitate-the-best or pairwise comparison rules. Cong et al. [25] introduced the reputation into the evolution model of the prisoner’s dilemma and analyzed the effects of population density, mobility sensitivity, and reputation memory on the emergence and persistence of cooperation. Li et al. [26] considered the trust between each player and its neighbors to study the coevolution of quantum strategies on scale-free networks.Moreover, a lot of works considering specific evolutionary mechanisms and interaction topologies have been contributed. In some recent reviews, corresponding advances as well as future directions have been well identified. The coevolutionary rules could be regarded as a natural upgrade of evolutionary games, but how kinds of rules influence the evolution of cooperation needs to be investigated. Motivated by this, Perc and Szolnoki [27] presented a detailed review on the effects of coevolutionary rules on dynamical interactions between players, population growth, teaching activity, mobility, and aging of players. Nowak suggested that there are five rules for the evolution of cooperation: kin selection, direct reciprocity, indirect reciprocity, network reciprocity, and group selection [28]. Specifically, Perc et al. [29] reviewed advances regarding the evolutionary dynamics of group interactions on structured populations and made a comparison with the results obtained on well-mixed populations. Interaction topology is another key affecting the evolution of cooperation. As summarized in [27], common interaction topologies include rock-paper-scissors game, the ultimatum game, and public goods game. Szolnoki et al. [30] presented a latest review on cyclic dominance in evolutionary games.Although a lot of simulation models and findings have been achieved, there are few well-accepted theoretical works focusing on how to demonstrate various simulation results with respect to dynamics of lattice models. To the best of our knowledge, Szabó and Fáth [15] introduced the technique of generalized mean-field to analyze lattice system, and another fascinating method is pair approximation which considers the pair correlations between nearest neighbors [31]. Hutson and Vickers [32] used partial differential equations to build reaction-diffusion models and formulate two-dimensional grid dynamics. Schweitzer et al. [10] put forward an innovative way to comprehend the spatial effect. As mentioned in Section 1, in this work we conduct simulations using the same von Neumann neighborhood with [10] to observe the effect of initial population and system size on the evolution of cooperation. Then, we extend the model with a more complicated Moore neighborhood and self-playing rule for each central player, in order to investigate critical conditions for the spatial coexistence or the invasion of players in the extended evolution model.


## 3. The Model

Our model is built on the structure of a two-dimensional lattice. Each cell is occupied by one player, and everyone is only connected with the nearest neighbors. There are two kinds of neighborhood we will investigate. One is von Neumann neighborhood where only four nearest neighbors are considered, as Figure1 shows. The other is Moore neighborhood where eight nearest neighbors are taken into account, as Figure 2 shows.Figure 1
Von Neumann neighborhood where each cell has four nearest neighbors.Figure 2
Moore neighborhood where each cell has eight nearest neighbors.Instead of the mean-field structure in global models, individuals can only make local contacts by playing prisoner’s dilemma with neighbors. Also, there are no iterated games, no complicated algorithms, and no memories. Every individual only plays games with neighbors once in every generation and is restricted to choose from two strategies: cooperate or defect. At the end of each generation, each player sums up its payoff against neighbors as its present fitness and decides what kind of strategy to use in the next generation. Players cannot determine what to do by collecting information to anticipate what others will do in the next generation but simply imitate the player whose strategy gives rise to the highest fitness. If there is a tie, every player always keeps its original strategy because switching strategy requires paying the price.The above updating rule has twofold meanings. From the perspective of biology, strategies with higher fitness can survive and thrive; they will spread gradually to take over the majority of cells, while the ones with lower fitness are doomed to die out. From the perspective of sociology, people tend to adopt strategies proved to be more successful. As we can see, this is a deterministic and synchronous updating rule. Players may be regarded as cellular automata which change their states synchronously at discrete time steps, and the next state of each cell depends on the current states of the neighboring cells according to an update rule. All cells use the same rule, and the rule is applied to all cells at the same time [33]. In the model, the array has periodic boundaries, so the marginal cells on one side of the array are regarded as the neighbors of the cells on the opposite margin of the lattice.In the lattice withn×n cells, we denote a finite set of cells by C={i,i=1,2,…,n×n} and a set of corresponding states by S={sit,i=1,2,…,n×n,t=1,2,…}, where sit means the state of cell i in the generation t. As mentioned above, every cell’s state is only concerned with its neighbors, so, in Moore structure, each cell has a neighborhood with eight cells. We denote neighbors of cell i by Ni={n0i,n1i,…,n8i}, in which n0i represents cell i itself. When cell i begins to play games with his neighborhood, it may receive the payoff Ui=∑j=08‍ui(i,nji), where ui(i,nji) represents the payoff of cell i over its neighbor nji. As a matter of fact, we count the payoff that cell i plays games with itself, which seems to be paradox but is very common in nature. We just need to regard every cell as a group, and then this consideration becomes reasonable. Last but not least, the force at the core in the dynamic system, that is, the updating rule, is F:sit+1=snjit, j=argmax⁡j=0,1,…,8Unji.Specifically speaking, there are only two states that every cell hassi=0, which means the cell defects, or si=1, which means the cell cooperates. The interaction between two cells is playing prisoner’s dilemma game, and either of them does not have any prior knowledge about what the other will do. After all players choose their strategies from C and D, they will receive the payoffs according to the payoff matrix in Table 1. T means temptation payoff for defecting in a cooperative environment, S means sucker’ payoff for cooperating in a defecting environment, R means reward payoff for cooperating in a likewise environment, and P means punishment payoff for defecting in a likewise environment [10].Table 1
Payoff matrix of prisoner’s dilemma.
Playeri
Playerj
C
D
C
(R, R)
(S, T)
D
(T, S)
(P, P)They both get payoffR as reward if they both cooperate, while they receive payoff P as punishment if they both defect. A player using D obtains the payoff T against a C-player, and the counterpart has to suffer the minimum income S as for sucker, which signifies that selfish player can exploit the cooperator. The source leading to the dilemma lies in T>R>P>S. To promote cooperation and maintain this difference, assuming 2R>T+S is also necessary, which implies that the only Pareto efficient result is that they both cooperate.


## 4. The Spatial Dynamics and Analysis

In Schweitzer et al.’s work [10], they consider a cell and its nearest neighbors as a group, and what happens among them is regarded as an n-person game even though the actual behavior among each other in this group is independent. The group game can be decomposed as a combination of two-person games happening n-1 times, and the payoff of the central cell in each generation is the sum of its payoffs by playing with every neighbor. From this perspective, the fitness of each cell is impacted intensely by the number of cooperators in the neighborhood. For instance, cell i has k neighbors who are cooperators under the structure of Moore neighborhood. Then it can get the payoff:(1)Ui=∑j=08‍uii,nji=k+1R+8-kS=R-Sk+R+8S.If its own statesi=1, then it can get the payoff:
(2)Ui=∑j=08‍uii,nji=kT+9-kP=T-Pk+9P.If its own statesi=0, the game of cell i ends in the generation.Due toT>R>P>S, payoff Ui is an increasing function of the number of cooperators. Thus, no matter what state a cell is in, the more cooperators there are in its neighborhood, the higher fitness it will receive. Since the updating rule is determined by the cell with the highest score in the group, it is essential to take into account the whole configuration. Let the configuration Qij indicate that the state of the central cell in its neighborhood is i and the number of cooperators in this group is j rather than standalone cell configuration. Apparently, the value of i can be 1 or 0, and 0≤j≤8 in the structure of Moore neighborhood. Each Qij corresponds to uij that indicates the payoff of cell i against the jth cooperator in its group. We can get some simple summaries from the definition:(1)
Ifj>k, then u0j>u0k and u1j>u1k (as we mentioned above);(2)
u0j+1>u1j, j=0,1,2,…,8, because
(3)u0j+1-u1j=jT+9-jP-jR-9-jS=jT-R+9-jP-S>0.Based on this theoretical foundation, Schweitzer et al. [10] did some simulations in the model of von Neumann structure. We extend the model to the Moore neighborhood and add the self-playing for every cell. The aim of this paper is to verify the accuracy of this theory and dig deeper into the application environment of the theory.The flash point that makes the theory so special is that it transforms invasion from a macroscopic view into a microstructural one. There is a concept of invasion for every group when the value ofi in Qij changes. If i turns into 1 from 0, it can be regarded that configuration Q0j invades configuration Q1j. In our model, there are 18 different kinds of configurations, that is, {Q00,Q01,Q02,Q03,Q04,Q05,Q06,Q07,Q08} and {Q10,Q11,Q12,Q13,Q14,Q15,Q16,Q17,Q18}. Since there is no room for invasion in the configuration Q00, we do not consider it. Configuration Q10 is extremely unstable which rapidly turns into Q00 after one generation, so we also rule it out.Through the above discussion, we have(4)u01<u02<u03<u04<u05<u06<u07<u08,u11<u12<u13<u14<u15<u16<u17<u18.So if the following condition is met,u18<u01, the highest score a cooperator can obtain is still lower than the lowest score a defector can receive. Then, cooperators are doomed to be invaded, and defectors are able to invade cooperators completely no matter what initial distribution they have. However, there is no way by which cooperators completely invade defectors because u11<u17<u08. Meanwhile, u1j<u0j+1, j=1,2,…,7 implies that structure Q1j must be invaded by structure Q0j+1. Since there is no limitation on u1j+1 and u0j+1, if u1j<u0j+1<u1j+1 is satisfied, the power of aggression of Q0j+1 will be crippled, and if u0j+1<u1j+1 is not satisfied, u0j+1<u1j+2 can serve as the second strict restrictive condition, and so forth. We may get several different patterns in lattice system with the relaxation of limitation on Q0j’s invasion.The payoff in classic prisoner’s dilemma isT=5, R=3, P=1, S=0. To investigate varying degrees of limitation on the invasion of Q0j, we set the related parameters as follows: R=3, S=0, u0j being a binary variable of T and P, and u1j as a constant. As a result, we could partition the P,T-plane according to the different rank ordering of R, S, T, and P.


### 4.1. The Spatial Dynamics of von Neumann Neighborhood without Self-Playing

In this section, like Schweitzer et al.’s work [10], each cell only communicates with its four nearest neighbors in von Neumann neighborhood and does not play game with itself. Here, u0j(j=1,2,3,4) provides four kinds of lines with different gradients, u1j(j=2,3,4) stands for three different intercepts, and each gradient is matched with several intercepts according to the rule u0j>u1j. In the P,T-plane, there are six lines intersecting with each other and cutting apart the plane with eleven areas. The lines are, respectively, represented by
(5)u01=u14⟹T+3P=12;u01=u13⟹T+3P=9;u01=u12⟹T+3P=6;u02=u14⟹T+P=6;u02=u13⟹T+P=4.5;u03=u14⟹3T+P=12.The inside plane formed by two crossover points is regarded as a classification basis for six lines. Obviously, lines intersecting at the same point can be classified into the same category. Therefore we have three sets:E1={l12,l23,l34}, E2={l13,l24}, and E3={l14}, in which lij represents straight line u0i=u1j. Each set also corresponds to a level of degree e1 that indicates how intense it influences various configurations, and e1 is the most intense class of all and so forth. Then we discuss distinct regions in accordance with simulation of spatial prisoner’s dilemma.Simulation I. The most stringent limitation for configuration Q0j which makes defectors too weak to recklessly invade cooperators is u0i<u1i+1, i=1,2,3, namely,
(6)u01<u12,u02<u13,u03<u14.From another point of view, this group of conditions enables cooperators with various sorts of configuration to invade the defectors of configuration in the process of evolution. Cooperation occupies the majority of domain in this situation as we expect, stabilizing above the average percentage of 78%. Defectors are clustered in linked curves, carving up the sea of cooperation as massive regions and sometimes accompanied by periodic oscillating cross, as presented in Figure3.Figure 3
Simulation results on a50×50 square lattice with von Neumann neighborhood and random initial distribution of players. This pattern is performed by the parameters region: u01<u12, u02<u13, u03<u14.In Figures3, 4, and 5 red cells represent players who defect in this generation and in the previous generation, blue is cooperation following cooperation, green is cooperation following defection, and yellow is defection following cooperation.Simulation results on a50×50 square lattice with von Neumann neighborhood and random initial distribution of players. The pattern in (a) is performed by the parameters region: u01>u12, u02<u13, u03<u14, and the pattern in (b) is by the region: u01<u12, u02<u13, u03>u14.
(a)
(b)Simulation results on a50×50 square lattice with von Neumann neighborhood and random initial distribution of players. The pattern in (a) is performed by the parameters region u01>u12, u02>u13, u03<u14, and the pattern in (b) is by the region u01<u12, u02>u13, u03>u14.
(a)
(b)Simulation II. We loosen the restriction for configuration Q0j to change the sign of inequality of one condition, which produces two different feasible regions:
(7)u01>u12,u02<u13,u03<u14,(8)u01<u12,u02<u13,u03>u14.In the situation formed by (7), the number of defectors increases sharply and captures the most territory under both circumstances. Cooperation is found to exist in the form of big clusters, making up about 16% of players, and it is hard to find any periodically oscillating pattern, as Figure 4(a) shows.The third situation is presented in Figure4(b) where a peculiar phenomenon reveals itself: there is no terminal point for the evolution, and the pattern is changing all the time with high level of oscillation. This phenomenon can be referred to as spatial chaos: at first, defection dominates most parts, but cooperators cluster to compete with rivals; every time any clusters of cooperation are close enough to each other, the defectors between them will take advantage of them to grow; as the number of their alliances reaches to a point, all of them are facing an unfavorable payoff resulting in that cooperation has the better hand again; repetition of devouring of both sides brings about the phenomenon we have seen. The defectors outnumber cooperators in Figure 4(b) in the ratio of three to two.Simulation III. We make further efforts to relax the restriction of defectors. The sign of inequality in two conditions has been changed, resulting in two different feasible regions:
(9)u01>u12,u02>u13,u03<u14,(10)u01<u12,u02>u13,u03>u14.Figure5(a) shows the final state under conditions (9). Cooperation is suffering to survive as small clusters of crosses reaching the average percentage of almost 7.5%, and the pattern is perfectly stable without oscillation. When we examine the matter from different angles, group of conditions (9) is apparently further relaxation for conditions (7), and the clusters have been transformed from big to small and from many to rare, as Figures 4(a) and 5(a) show.Group of conditions (10), which is the relaxation of the third, yields a sort of special pattern for von Neumann neighborhood. Cooperators grow as crosses, following the trace of the four directions, and there are always some separately oscillating crosses, as Figure 5(b) shows. In a nutshell, the ratio ofT toP determines the stability level of the ultimate state. The larger the ratio is, the more oscillating the crosses will appear.Simulation IV. Removing all the limitations for defectors in the first class, we have
(11)u01>u12,u02>u13,u03>u14.The simulated pattern shows the domination of defection and the average percentage of cooperation does not exceed 0.6%. There is absolutely no cooperator under the assumption of initial random distribution especially whenu02>u14, u01<u14.Discussions. From the above simulation results, we have observed the phenomena of steady or oscillating coexistence between cooperators and defectors. These results verify that cooperators and defectors could coexist in the evolution if the local spatial structure is considered. Meanwhile, different combinations of conditions on players have resulted in different evolutions. What the results indicate is somewhat similar with the underlying evolution of real-world nature and society. The neighborhood of each cell is like the relationship of each living being, and the parameters region is like the role of local environment setting.We already know that cooperators are inevitable to die out whileu01>u14 no matter what the initial distribution is. As a matter of fact, we have observed huge differences of results from those in Schweitzer et al.’s work [10], especially in the regions {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. They showed that these regions stand for the case of unstable coexistence, which means the dynamics will always lead to an invasion of the defectors except for some special case. However, our simulations show that there are always stable clusters in the regions and sometimes even very big clusters.After making a detailed comparison with [10], we find that the main differences of our simulations from Schweitzer et al.’s are the initial population of players and the lattice size. In [10], the initial distribution was randomly generated such that the number of defectors was equal to that of cooperators, and the lattice size was 40×40; In this work, we randomly assign each player on a 50×50 lattice as 0 (to defect) or 1 (to cooperate), so the number of initial defectors is basically equivalent to that of cooperators. Thus, the inconsistent results from [10] may be owing to the different initial structures of players and the lattice size.In order to verify the impact of the randomness of initial structures on evolutionary coexistence, we conducted additional simulations, respectively, with regions{u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. Like [10], we use the global frequency of cooperating players to show the evolution of these simulations, as Figure 6 shows.Simulation results on a50×50 square lattice with von Neumann neighborhood and different random initial distributions of players. The x-axis represents the number of generation, and y-axis represents the frequency of cooperators in the evolution. The results in (a) and (b) are, respectively, performed by the parameters region {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. f(0) means the frequency of cooperators in the random initial distribution.
(a)
(b)From the results in Figure6, we could find (i) for the two regions {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}, all the conducted evolutions although with different random initial distributions finally came into the stable coexistence of cooperators and defectors; (ii) the frequency of cooperators in the initial distribution has direct impact on the ratio of players in the final coexistence.In order to investigate the effect of population size on the evolution, we extend the simulations to square lattices in different sizes. As Figure7 shows, we obtain the final frequencies of cooperators with lattice sizes ranging from 10×10 to 80×80. Due to space limitations, we just give the results performed by the parameters region {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. From the results in Figure 7, we could have the following observations:(1)
For the parameters region{u01>u12,u02<u13,u03<u14}, the final states on all the lattices evolve into steady coexistence of cooperators and defectors, with the slightly fluctuated frequencies of cooperators in the final patterns.(2)
For the parameters region{u01>u12,u02>u13,u03<u14}, the pattern turns into the complete invasion by defectors when the square lattice is 10×10, but all other patterns evolve into steady coexistence of cooperators and defectors.Simulation results on square lattices in different sizes with von Neumann neighborhood andf(0)=0.50. The x-axis represents the lattice size, and y-axis represents the frequency of cooperators in the final pattern. The results in (a) are performed by the parameters region {u01>u12,u02<u13,u03<u14}, and the results in (b) are by the region {u01>u12,u02>u13,u03<u14}.
(a)
(b)To make more comprehensive investigation of the validity of our results, in the following subsection we extend the model to a more complicated level: the neighborhood turns into the formation of Moore and each cell can play games with itself. The self-playing rule is inspired by the inheritance of species in nature, which would have the side effect to improve the performance of cooperators.


### 4.2. The Spatial Dynamics of Moore Neighborhood with Self-Playing

As mentioned above, with Moore neighborhood, there are eight sorts of configurationQ0j which represent eight sorts of gradients of straight lines and eight sorts of configuration Q1j which imply eight sorts of intersects. Synthesizing various conditions presents thirty-six lines. According to whether they are intersecting at the same crossover point, these lines can be categorized into eight subsets:
(12)E1=l11,l22,l33,l44,l55,l66,l77,l88,E2=l12,l23,l34,l45,l56,l67,l78,E3=l13,l24,l35,l46,l57,l68,E4=l14,l25,l36,l47,l58,E5=l15,l26,l37,l48,E6=l16,l27,l38,E7=l17,l28,E8=l18.E
1 embodies the relationship between two entities with homogeneous level (the same amount of cooperators). Once there is a limitation up to someone in this category, it will have the greatest strength to affect outcomes. Following the case of von Neumann, we create a degree ek to measure the influence that Ek could exert on any configuration. We can easily deduce ek≥ej, k≥j. In contrast with absence of self-playing, the significant parts lie in {u1i<u0i+1,i=1,2,…,7} and the fact that Q08 is no longer the invincible structure. Q18 is reasonable to be superior to Q08 here.Since overwhelming lines need to be analyzed individually, in Figure8 we partition the P,T-plane according to the different patterns of spatial dynamics.Figure 8
Partitioning of theP,T-plane according to the different patterns of spatial dynamics.Simulation V. Cooperation holds the trump card in region A formed by u0i<u1i+2(i=1,2,…,6) with the average percentage above 90%, as shown in Figure 9. Defectors survive as single dot or oscillating cluster with the size of 3×3, and sometimes there are a few inconsistent short lines.Figure 9
Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern is performed by the parameters region u0i<u1i+2(i=1,2,…,6).Simulation VI. While we make modifications to one of the inequalities u0i<u1i+2(i=1,2,…,6), only two groups of inequalities turn out to be feasible:
(13)u01>u13,u0i<u1i+2i=2,3,…,6,(14)u06>u18,u0i<u1i+2i=1,2,…,5.Under above two circumstances, defectors join up as integrated lines to partition the whole territory of cooperation into pieces, which manifests the enhancement of their strength of incursion. However, great distinctions still remain in the two regions. When conditions (13) are met, only steady clusters of defectors emerge, leaving oscillating clusters no place to live. Nevertheless, conditions (14) incubate various states with high level of oscillation, and every line or cluster formed by defectors is surrounded by oscillating cells with two periods. Figures 10(a) and 10(b) show the corresponding patterns.Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern in (a) is performed by the parameters region u01>u13, u0i<u1i+2(i=2,3,…,6), and the pattern in (b) is by the region u06>u18, u0i<u1i+2(i=1,2,…,5).
(a)
(b)Simulation VII. To make detailed exploration on how the dynamics changes gradually along with the relaxation of restriction of defection from u0i<u1i+2(i=1,2,…,6), we modify two inequalities to achieve two feasible regions:
(15)u01>u13,u02>u14,u0i<u1i+2i=3,4,…,6,(16)u05>u17,u06>u18,u0i<u1i+2i=1,2,…,4.These two regions are the expansion from conditions (13) and (14), respectively, and the dynamics are also in the wake of invasion degree of defection becoming worse in terms of cooperation, whose patterns are shown in Figures 11(a) and 11(b). It is reflected by the fact that there are more stable lines of defectors. Conditions (13) combined with conditions (15) form region B in Figure 8, and conditions (14) and conditions (16) form the region C.Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern in (a) is performed by the parameters region u01>u13, u02>u14, u0i<u1i+2(i=3,4,…,6), and the pattern in (b) is by the region u05>u17, u06>u18, u0i<u1i+2(i=1,2,…,4).
(a)
(b)Simulation VIII. In the light of the law of development, three conditions need to be changed. Region u0i>u1i+2i=4,5,6, u0i<u1i+2(i=1,2,3) assumes hybrid dynamics of spatial chaos and majority of cooperation. u05<u18 (region D) holds the predominance of cooperation: clusters as the shape of tentacles are formed inside the picture and hardly touch the edges, as Figure 12(a) shows.Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern in (a) is performed by the parameters region u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05<u18, and the pattern in (b) is by the region u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05>u18.
(a)
(b)When it comes tou05>u18 (region E), there comes the chaos, as presented in Figure 12(b), and we notice that the chaos pattern resulting from Moore neighborhood is quite different from that from von Neumann neighborhood. The former starts with a local fascinating pattern continuing migrating and slowly diffuses from small plots of cooperation in the territory of defectors until the whole is in the chaotic state that is never able to be restored again. Meanwhile, every plot of cooperation is engulfed repeatedly by defectors from inside. While the latter is relatively stable, as soon as cooperators get a firm foothold, they only stretch and shrink on the boundary and the percentage of cooperation stabilizes rapidly.Simulation IX. The conditions are as follows:
(17)u0i>u1i+2i=1,2,3,u0i<u1i+2i=4,5,6,u0i>u1i+2i=1,2,…,4,u0i<u1i+2i=5,6,u0i>u1i+2i=1,2,…,5,u06<u18.These conditions (region F) are the watershed between the domination of cooperation and defection, forming the patterns similar with that in Figure13. The force of cooperation has been deteriorated by opponents, reaching the state where cooperators survive in the shape of small rectangle clusters. This reflects another distinction between two structures of neighborhood. A crisscross pattern formed by cooperators is another choice for von Neumann neighborhood to cluster.Figure 13
Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern is performed by the parameters region u0i>u1i+2(i=1,2,3), u0i<u1i+2(i=4,5,6); u0i>u1i+2(i=1,2,…,4), u0i<u1i+2(i=5,6); u0i>u1i+2(i=1,2,…,5), u06<u18.Other Simulations. The conditions are as follows:
(18)u01<u13,u02<u14,u0i>u1i+2i=3,4,…,6,u01<u13,u0i>u1i+2i=2,3,…,6.Defectors totally control the whole situation and fill up all the cells except whenu05<u18. u05>u18 sometimes causes very few rectangle clusters of cooperators but always brings about the complete domination of defectors if the initial distribution is random. For some special initial distributions such that there is only one defector in the central cell, defection is not able to invade them all.u
0
i
>
u
1
i
+
2
(
i
=
1,2
,
…
,
6
) releases the force of incursion thoroughly on the basis of level of e3. This region is very much like the last one, which has few clusters when extra condition u05<u18 is satisfied. Together with inequations (18), the region G is regarded as almost complete defected area.Onceu18<u01 (region H) is met, defection will completely dominate cooperation.Discussions. Every group of straight lines in Ek divides the picture with similar presentation and the number of lines in Ek also decreases gradually. The area corresponding to E1 is
(19)P1=x,y∣y<278-18xory<6-8x
which can be divided into three pieces:
(20)P11=x,y∣y<278-18x,y<6-8x,P12=x,y∣y<278-18x,y>6-8x,P13=x,y∣y>278-18x,y<6-8x.The level of cooperation inP11 certainly is highest among all the three. P12 turns out to be more oscillating and the clusters formed in this region are generally in the shape of curve, and P13 intends to be more stable without oscillation and the clusters are found in the shape of lump. The patterns forming the similar clusters like P13 are also different, and the lines in E1 separate the level of cooperation in P13 into secondary classes changing along with the direction of radiation. The same law can be applied to other regions Pk(k=2,3). The extent of emerging of cooperation is diminishing from Pk to Pk+1(k=1,2,…,7).The above discussion is subject to the dynamics on the basis ofe3, rather than the case of von Neumann which we start with e1. The reason is that the dynamics of Moore structure are far more complicated than that of von Neumann structure, and the influence of e3 sufficiently brings the majority of cooperation to the system. Since each group of lines in Ek corresponds to a set of constraint conditions, e1 obviously has the strictest constraint on configuration Q0j which is {u0i<u1i,i=1,2,…,6}. Just as we expected, defectors suffer as independent dots and the number is rare. Every time we convert an inequality to unleash a level of constraint for Q0j, e2 would put on some opposite strength to prevent Q0j from inflating excessively. These forces of expansion and restriction result in various amazing patterns.Our modification is made from both sides ofEk(k=1,2,3), so we come to the conclusion that u01<u11+(k-1) has less power to undermine the defectors than u06<u16+(k-1). The same rule works for {u01<u11+(k-1),u02<u12+(k-1)}, {u06<u16+(k-1),u05<u15+(k-1)}, and so forth. Simulation results show that the regions dominated by cooperators in Moore neighborhood are larger than those in von Neumann, and the clusters of defectors in the former structure are in diversified shapes which are much more than those in the later. All the phenomena demonstrate the fact that the more neighbors surround the central cell, the more cooperators would be produced.With regard to the state of occupation of cooperation, the procedure of evolution is relatively slow and regular: at the beginning defectors rule the most part and cooperators form the clusters of local and scattered dots; then cooperators develop gradually. The cooperated state with oscillation is generated by major and small clusters; the less cooperated state with stable defectors is usually formed by minor and big clusters. The spatial chaos reveals the range that promotes the situation in both neighborhood structures is generally the same. When the value ofT/P is over a threshold, the chaos appears; after it reaches some limit, the chaos will vanish again.In order to test the above results against the population size, we select six typical parameters regions to conduct simulations with Moore neighborhood and self-playing on square lattices in different sizes. As Figure14 shows, we obtain the final frequencies of cooperators with lattice sizes ranging from 10×10 to 80×80. Note that the two curves in the oscillating patterns like Figure 14(a) and chaotic patterns like Figure 14(e), respectively, represent the maximal and minimal frequencies of cooperators. From the results in Figure 14, we can find the following.(1)
The steady patterns performed by parameters regionsu01>u13, u0i<u1i+2(i=2,3,…,6) and u0i>u1i+2(i=1,2,3), u0i<u1i+2(i=4,5,6) almost keep the coexistence of cooperators and defectors on square lattices in different sizes with Moore neighborhood and self-playing, as Figures 14(b) and 14(f) show. However, the evolution with Moore neighborhood may also turn into the complete invasion of defectors if the population size is not big enough, as Figure 14(f) shows. Together with the extinction phenomenon in Figure 7(b), we find that it is easy to produce complete invasion of defectors in the evolution with loose constraints on defectors and small population size.(2)
The population size has the impact on the evolution of oscillating patterns such as Figures14(a), 14(c), and 14(d), but the impact is not big enough to turn the oscillating patterns into complete invasion of cooperators or defectors. The similar impact is observed for the chaotic patterns, as Figure 14(e) shows.(3)
Comparing the results with those on the50×50 lattice with Moore neighborhood and self-playing, we basically test the effectiveness of the results against population size. The results on lattices in different sizes provide further evidence on that the Moore neighborhood and self-playing rule could promote the cooperation in evolutionary games.Simulation results on square lattices in different sizes with Moore neighborhood andf(0)=0.50. The x-axis represents the lattice size, and y-axis represents the frequency of cooperators in the final pattern. The results in (a)–(f) are, respectively, performed by the parameters region u0i<u1i+2(i=1,2,…,6); u01>u13, u0i<u1i+2(i=2,3,…,6); u06>u18, u0i<u1i+2(i=1,2,…,5); u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05<u18; u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05>u18, and u0i>u1i+2(i=1,2,3), u0i<u1i+2(i=4,5,6).
(a)
(b)
(c)
(d)
(e)
(f)


## 4.1. The Spatial Dynamics of von Neumann Neighborhood without Self-Playing

In this section, like Schweitzer et al.’s work [10], each cell only communicates with its four nearest neighbors in von Neumann neighborhood and does not play game with itself. Here, u0j(j=1,2,3,4) provides four kinds of lines with different gradients, u1j(j=2,3,4) stands for three different intercepts, and each gradient is matched with several intercepts according to the rule u0j>u1j. In the P,T-plane, there are six lines intersecting with each other and cutting apart the plane with eleven areas. The lines are, respectively, represented by
(5)u01=u14⟹T+3P=12;u01=u13⟹T+3P=9;u01=u12⟹T+3P=6;u02=u14⟹T+P=6;u02=u13⟹T+P=4.5;u03=u14⟹3T+P=12.The inside plane formed by two crossover points is regarded as a classification basis for six lines. Obviously, lines intersecting at the same point can be classified into the same category. Therefore we have three sets:E1={l12,l23,l34}, E2={l13,l24}, and E3={l14}, in which lij represents straight line u0i=u1j. Each set also corresponds to a level of degree e1 that indicates how intense it influences various configurations, and e1 is the most intense class of all and so forth. Then we discuss distinct regions in accordance with simulation of spatial prisoner’s dilemma.Simulation I. The most stringent limitation for configuration Q0j which makes defectors too weak to recklessly invade cooperators is u0i<u1i+1, i=1,2,3, namely,
(6)u01<u12,u02<u13,u03<u14.From another point of view, this group of conditions enables cooperators with various sorts of configuration to invade the defectors of configuration in the process of evolution. Cooperation occupies the majority of domain in this situation as we expect, stabilizing above the average percentage of 78%. Defectors are clustered in linked curves, carving up the sea of cooperation as massive regions and sometimes accompanied by periodic oscillating cross, as presented in Figure3.Figure 3
Simulation results on a50×50 square lattice with von Neumann neighborhood and random initial distribution of players. This pattern is performed by the parameters region: u01<u12, u02<u13, u03<u14.In Figures3, 4, and 5 red cells represent players who defect in this generation and in the previous generation, blue is cooperation following cooperation, green is cooperation following defection, and yellow is defection following cooperation.Simulation results on a50×50 square lattice with von Neumann neighborhood and random initial distribution of players. The pattern in (a) is performed by the parameters region: u01>u12, u02<u13, u03<u14, and the pattern in (b) is by the region: u01<u12, u02<u13, u03>u14.
(a)
(b)Simulation results on a50×50 square lattice with von Neumann neighborhood and random initial distribution of players. The pattern in (a) is performed by the parameters region u01>u12, u02>u13, u03<u14, and the pattern in (b) is by the region u01<u12, u02>u13, u03>u14.
(a)
(b)Simulation II. We loosen the restriction for configuration Q0j to change the sign of inequality of one condition, which produces two different feasible regions:
(7)u01>u12,u02<u13,u03<u14,(8)u01<u12,u02<u13,u03>u14.In the situation formed by (7), the number of defectors increases sharply and captures the most territory under both circumstances. Cooperation is found to exist in the form of big clusters, making up about 16% of players, and it is hard to find any periodically oscillating pattern, as Figure 4(a) shows.The third situation is presented in Figure4(b) where a peculiar phenomenon reveals itself: there is no terminal point for the evolution, and the pattern is changing all the time with high level of oscillation. This phenomenon can be referred to as spatial chaos: at first, defection dominates most parts, but cooperators cluster to compete with rivals; every time any clusters of cooperation are close enough to each other, the defectors between them will take advantage of them to grow; as the number of their alliances reaches to a point, all of them are facing an unfavorable payoff resulting in that cooperation has the better hand again; repetition of devouring of both sides brings about the phenomenon we have seen. The defectors outnumber cooperators in Figure 4(b) in the ratio of three to two.Simulation III. We make further efforts to relax the restriction of defectors. The sign of inequality in two conditions has been changed, resulting in two different feasible regions:
(9)u01>u12,u02>u13,u03<u14,(10)u01<u12,u02>u13,u03>u14.Figure5(a) shows the final state under conditions (9). Cooperation is suffering to survive as small clusters of crosses reaching the average percentage of almost 7.5%, and the pattern is perfectly stable without oscillation. When we examine the matter from different angles, group of conditions (9) is apparently further relaxation for conditions (7), and the clusters have been transformed from big to small and from many to rare, as Figures 4(a) and 5(a) show.Group of conditions (10), which is the relaxation of the third, yields a sort of special pattern for von Neumann neighborhood. Cooperators grow as crosses, following the trace of the four directions, and there are always some separately oscillating crosses, as Figure 5(b) shows. In a nutshell, the ratio ofT toP determines the stability level of the ultimate state. The larger the ratio is, the more oscillating the crosses will appear.Simulation IV. Removing all the limitations for defectors in the first class, we have
(11)u01>u12,u02>u13,u03>u14.The simulated pattern shows the domination of defection and the average percentage of cooperation does not exceed 0.6%. There is absolutely no cooperator under the assumption of initial random distribution especially whenu02>u14, u01<u14.Discussions. From the above simulation results, we have observed the phenomena of steady or oscillating coexistence between cooperators and defectors. These results verify that cooperators and defectors could coexist in the evolution if the local spatial structure is considered. Meanwhile, different combinations of conditions on players have resulted in different evolutions. What the results indicate is somewhat similar with the underlying evolution of real-world nature and society. The neighborhood of each cell is like the relationship of each living being, and the parameters region is like the role of local environment setting.We already know that cooperators are inevitable to die out whileu01>u14 no matter what the initial distribution is. As a matter of fact, we have observed huge differences of results from those in Schweitzer et al.’s work [10], especially in the regions {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. They showed that these regions stand for the case of unstable coexistence, which means the dynamics will always lead to an invasion of the defectors except for some special case. However, our simulations show that there are always stable clusters in the regions and sometimes even very big clusters.After making a detailed comparison with [10], we find that the main differences of our simulations from Schweitzer et al.’s are the initial population of players and the lattice size. In [10], the initial distribution was randomly generated such that the number of defectors was equal to that of cooperators, and the lattice size was 40×40; In this work, we randomly assign each player on a 50×50 lattice as 0 (to defect) or 1 (to cooperate), so the number of initial defectors is basically equivalent to that of cooperators. Thus, the inconsistent results from [10] may be owing to the different initial structures of players and the lattice size.In order to verify the impact of the randomness of initial structures on evolutionary coexistence, we conducted additional simulations, respectively, with regions{u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. Like [10], we use the global frequency of cooperating players to show the evolution of these simulations, as Figure 6 shows.Simulation results on a50×50 square lattice with von Neumann neighborhood and different random initial distributions of players. The x-axis represents the number of generation, and y-axis represents the frequency of cooperators in the evolution. The results in (a) and (b) are, respectively, performed by the parameters region {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. f(0) means the frequency of cooperators in the random initial distribution.
(a)
(b)From the results in Figure6, we could find (i) for the two regions {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}, all the conducted evolutions although with different random initial distributions finally came into the stable coexistence of cooperators and defectors; (ii) the frequency of cooperators in the initial distribution has direct impact on the ratio of players in the final coexistence.In order to investigate the effect of population size on the evolution, we extend the simulations to square lattices in different sizes. As Figure7 shows, we obtain the final frequencies of cooperators with lattice sizes ranging from 10×10 to 80×80. Due to space limitations, we just give the results performed by the parameters region {u01>u12,u02<u13,u03<u14} and {u01>u12,u02>u13,u03<u14}. From the results in Figure 7, we could have the following observations:(1)
For the parameters region{u01>u12,u02<u13,u03<u14}, the final states on all the lattices evolve into steady coexistence of cooperators and defectors, with the slightly fluctuated frequencies of cooperators in the final patterns.(2)
For the parameters region{u01>u12,u02>u13,u03<u14}, the pattern turns into the complete invasion by defectors when the square lattice is 10×10, but all other patterns evolve into steady coexistence of cooperators and defectors.Simulation results on square lattices in different sizes with von Neumann neighborhood andf(0)=0.50. The x-axis represents the lattice size, and y-axis represents the frequency of cooperators in the final pattern. The results in (a) are performed by the parameters region {u01>u12,u02<u13,u03<u14}, and the results in (b) are by the region {u01>u12,u02>u13,u03<u14}.
(a)
(b)To make more comprehensive investigation of the validity of our results, in the following subsection we extend the model to a more complicated level: the neighborhood turns into the formation of Moore and each cell can play games with itself. The self-playing rule is inspired by the inheritance of species in nature, which would have the side effect to improve the performance of cooperators.


## 4.2. The Spatial Dynamics of Moore Neighborhood with Self-Playing

As mentioned above, with Moore neighborhood, there are eight sorts of configurationQ0j which represent eight sorts of gradients of straight lines and eight sorts of configuration Q1j which imply eight sorts of intersects. Synthesizing various conditions presents thirty-six lines. According to whether they are intersecting at the same crossover point, these lines can be categorized into eight subsets:
(12)E1=l11,l22,l33,l44,l55,l66,l77,l88,E2=l12,l23,l34,l45,l56,l67,l78,E3=l13,l24,l35,l46,l57,l68,E4=l14,l25,l36,l47,l58,E5=l15,l26,l37,l48,E6=l16,l27,l38,E7=l17,l28,E8=l18.E
1 embodies the relationship between two entities with homogeneous level (the same amount of cooperators). Once there is a limitation up to someone in this category, it will have the greatest strength to affect outcomes. Following the case of von Neumann, we create a degree ek to measure the influence that Ek could exert on any configuration. We can easily deduce ek≥ej, k≥j. In contrast with absence of self-playing, the significant parts lie in {u1i<u0i+1,i=1,2,…,7} and the fact that Q08 is no longer the invincible structure. Q18 is reasonable to be superior to Q08 here.Since overwhelming lines need to be analyzed individually, in Figure8 we partition the P,T-plane according to the different patterns of spatial dynamics.Figure 8
Partitioning of theP,T-plane according to the different patterns of spatial dynamics.Simulation V. Cooperation holds the trump card in region A formed by u0i<u1i+2(i=1,2,…,6) with the average percentage above 90%, as shown in Figure 9. Defectors survive as single dot or oscillating cluster with the size of 3×3, and sometimes there are a few inconsistent short lines.Figure 9
Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern is performed by the parameters region u0i<u1i+2(i=1,2,…,6).Simulation VI. While we make modifications to one of the inequalities u0i<u1i+2(i=1,2,…,6), only two groups of inequalities turn out to be feasible:
(13)u01>u13,u0i<u1i+2i=2,3,…,6,(14)u06>u18,u0i<u1i+2i=1,2,…,5.Under above two circumstances, defectors join up as integrated lines to partition the whole territory of cooperation into pieces, which manifests the enhancement of their strength of incursion. However, great distinctions still remain in the two regions. When conditions (13) are met, only steady clusters of defectors emerge, leaving oscillating clusters no place to live. Nevertheless, conditions (14) incubate various states with high level of oscillation, and every line or cluster formed by defectors is surrounded by oscillating cells with two periods. Figures 10(a) and 10(b) show the corresponding patterns.Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern in (a) is performed by the parameters region u01>u13, u0i<u1i+2(i=2,3,…,6), and the pattern in (b) is by the region u06>u18, u0i<u1i+2(i=1,2,…,5).
(a)
(b)Simulation VII. To make detailed exploration on how the dynamics changes gradually along with the relaxation of restriction of defection from u0i<u1i+2(i=1,2,…,6), we modify two inequalities to achieve two feasible regions:
(15)u01>u13,u02>u14,u0i<u1i+2i=3,4,…,6,(16)u05>u17,u06>u18,u0i<u1i+2i=1,2,…,4.These two regions are the expansion from conditions (13) and (14), respectively, and the dynamics are also in the wake of invasion degree of defection becoming worse in terms of cooperation, whose patterns are shown in Figures 11(a) and 11(b). It is reflected by the fact that there are more stable lines of defectors. Conditions (13) combined with conditions (15) form region B in Figure 8, and conditions (14) and conditions (16) form the region C.Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern in (a) is performed by the parameters region u01>u13, u02>u14, u0i<u1i+2(i=3,4,…,6), and the pattern in (b) is by the region u05>u17, u06>u18, u0i<u1i+2(i=1,2,…,4).
(a)
(b)Simulation VIII. In the light of the law of development, three conditions need to be changed. Region u0i>u1i+2i=4,5,6, u0i<u1i+2(i=1,2,3) assumes hybrid dynamics of spatial chaos and majority of cooperation. u05<u18 (region D) holds the predominance of cooperation: clusters as the shape of tentacles are formed inside the picture and hardly touch the edges, as Figure 12(a) shows.Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern in (a) is performed by the parameters region u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05<u18, and the pattern in (b) is by the region u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05>u18.
(a)
(b)When it comes tou05>u18 (region E), there comes the chaos, as presented in Figure 12(b), and we notice that the chaos pattern resulting from Moore neighborhood is quite different from that from von Neumann neighborhood. The former starts with a local fascinating pattern continuing migrating and slowly diffuses from small plots of cooperation in the territory of defectors until the whole is in the chaotic state that is never able to be restored again. Meanwhile, every plot of cooperation is engulfed repeatedly by defectors from inside. While the latter is relatively stable, as soon as cooperators get a firm foothold, they only stretch and shrink on the boundary and the percentage of cooperation stabilizes rapidly.Simulation IX. The conditions are as follows:
(17)u0i>u1i+2i=1,2,3,u0i<u1i+2i=4,5,6,u0i>u1i+2i=1,2,…,4,u0i<u1i+2i=5,6,u0i>u1i+2i=1,2,…,5,u06<u18.These conditions (region F) are the watershed between the domination of cooperation and defection, forming the patterns similar with that in Figure13. The force of cooperation has been deteriorated by opponents, reaching the state where cooperators survive in the shape of small rectangle clusters. This reflects another distinction between two structures of neighborhood. A crisscross pattern formed by cooperators is another choice for von Neumann neighborhood to cluster.Figure 13
Simulation results on a50×50 square lattice with Moore neighborhood and self-playing. The initial distribution of cooperators and defectors is randomly generated. The pattern is performed by the parameters region u0i>u1i+2(i=1,2,3), u0i<u1i+2(i=4,5,6); u0i>u1i+2(i=1,2,…,4), u0i<u1i+2(i=5,6); u0i>u1i+2(i=1,2,…,5), u06<u18.Other Simulations. The conditions are as follows:
(18)u01<u13,u02<u14,u0i>u1i+2i=3,4,…,6,u01<u13,u0i>u1i+2i=2,3,…,6.Defectors totally control the whole situation and fill up all the cells except whenu05<u18. u05>u18 sometimes causes very few rectangle clusters of cooperators but always brings about the complete domination of defectors if the initial distribution is random. For some special initial distributions such that there is only one defector in the central cell, defection is not able to invade them all.u
0
i
>
u
1
i
+
2
(
i
=
1,2
,
…
,
6
) releases the force of incursion thoroughly on the basis of level of e3. This region is very much like the last one, which has few clusters when extra condition u05<u18 is satisfied. Together with inequations (18), the region G is regarded as almost complete defected area.Onceu18<u01 (region H) is met, defection will completely dominate cooperation.Discussions. Every group of straight lines in Ek divides the picture with similar presentation and the number of lines in Ek also decreases gradually. The area corresponding to E1 is
(19)P1=x,y∣y<278-18xory<6-8x
which can be divided into three pieces:
(20)P11=x,y∣y<278-18x,y<6-8x,P12=x,y∣y<278-18x,y>6-8x,P13=x,y∣y>278-18x,y<6-8x.The level of cooperation inP11 certainly is highest among all the three. P12 turns out to be more oscillating and the clusters formed in this region are generally in the shape of curve, and P13 intends to be more stable without oscillation and the clusters are found in the shape of lump. The patterns forming the similar clusters like P13 are also different, and the lines in E1 separate the level of cooperation in P13 into secondary classes changing along with the direction of radiation. The same law can be applied to other regions Pk(k=2,3). The extent of emerging of cooperation is diminishing from Pk to Pk+1(k=1,2,…,7).The above discussion is subject to the dynamics on the basis ofe3, rather than the case of von Neumann which we start with e1. The reason is that the dynamics of Moore structure are far more complicated than that of von Neumann structure, and the influence of e3 sufficiently brings the majority of cooperation to the system. Since each group of lines in Ek corresponds to a set of constraint conditions, e1 obviously has the strictest constraint on configuration Q0j which is {u0i<u1i,i=1,2,…,6}. Just as we expected, defectors suffer as independent dots and the number is rare. Every time we convert an inequality to unleash a level of constraint for Q0j, e2 would put on some opposite strength to prevent Q0j from inflating excessively. These forces of expansion and restriction result in various amazing patterns.Our modification is made from both sides ofEk(k=1,2,3), so we come to the conclusion that u01<u11+(k-1) has less power to undermine the defectors than u06<u16+(k-1). The same rule works for {u01<u11+(k-1),u02<u12+(k-1)}, {u06<u16+(k-1),u05<u15+(k-1)}, and so forth. Simulation results show that the regions dominated by cooperators in Moore neighborhood are larger than those in von Neumann, and the clusters of defectors in the former structure are in diversified shapes which are much more than those in the later. All the phenomena demonstrate the fact that the more neighbors surround the central cell, the more cooperators would be produced.With regard to the state of occupation of cooperation, the procedure of evolution is relatively slow and regular: at the beginning defectors rule the most part and cooperators form the clusters of local and scattered dots; then cooperators develop gradually. The cooperated state with oscillation is generated by major and small clusters; the less cooperated state with stable defectors is usually formed by minor and big clusters. The spatial chaos reveals the range that promotes the situation in both neighborhood structures is generally the same. When the value ofT/P is over a threshold, the chaos appears; after it reaches some limit, the chaos will vanish again.In order to test the above results against the population size, we select six typical parameters regions to conduct simulations with Moore neighborhood and self-playing on square lattices in different sizes. As Figure14 shows, we obtain the final frequencies of cooperators with lattice sizes ranging from 10×10 to 80×80. Note that the two curves in the oscillating patterns like Figure 14(a) and chaotic patterns like Figure 14(e), respectively, represent the maximal and minimal frequencies of cooperators. From the results in Figure 14, we can find the following.(1)
The steady patterns performed by parameters regionsu01>u13, u0i<u1i+2(i=2,3,…,6) and u0i>u1i+2(i=1,2,3), u0i<u1i+2(i=4,5,6) almost keep the coexistence of cooperators and defectors on square lattices in different sizes with Moore neighborhood and self-playing, as Figures 14(b) and 14(f) show. However, the evolution with Moore neighborhood may also turn into the complete invasion of defectors if the population size is not big enough, as Figure 14(f) shows. Together with the extinction phenomenon in Figure 7(b), we find that it is easy to produce complete invasion of defectors in the evolution with loose constraints on defectors and small population size.(2)
The population size has the impact on the evolution of oscillating patterns such as Figures14(a), 14(c), and 14(d), but the impact is not big enough to turn the oscillating patterns into complete invasion of cooperators or defectors. The similar impact is observed for the chaotic patterns, as Figure 14(e) shows.(3)
Comparing the results with those on the50×50 lattice with Moore neighborhood and self-playing, we basically test the effectiveness of the results against population size. The results on lattices in different sizes provide further evidence on that the Moore neighborhood and self-playing rule could promote the cooperation in evolutionary games.Simulation results on square lattices in different sizes with Moore neighborhood andf(0)=0.50. The x-axis represents the lattice size, and y-axis represents the frequency of cooperators in the final pattern. The results in (a)–(f) are, respectively, performed by the parameters region u0i<u1i+2(i=1,2,…,6); u01>u13, u0i<u1i+2(i=2,3,…,6); u06>u18, u0i<u1i+2(i=1,2,…,5); u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05<u18; u0i>u1i+2(i=4,5,6), u0i<u1i+2(i=1,2,3), u05>u18, and u0i>u1i+2(i=1,2,3), u0i<u1i+2(i=4,5,6).
(a)
(b)
(c)
(d)
(e)
(f)


## 5. Conclusions

Evolutionary game is a powerful framework to study the evolution of cooperation, and different spatial structures may bring about different patterns of the invasion or the coexistence of cooperator and defectors. In this work, we conducted lattice-based simulations with von Neumann neighborhood as Schweitzer et al. did and found some inconsistent results. In order to find the cause, the impact of initial population distribution and the lattice size on the evolution has been investigated. What’s more, we have extended the model to a complicated level where the neighborhood turns into the structure of Moore and each cell can play games with itself.Simulation results provided new evidence for the persistence of cooperation in the evolution with spatial structures and obtain some observations: (i) both the randomness of initial population distribution and the size of lattice have direct impact on the ratio of players in the final patterns; (ii) it is easy to produce complete invasion of defectors in the evolution with loose constraints on defectors and small population size; and (iii) comparing with the model with von Neumann neighborhood, our extended model with Moore neighborhood and self-playing rule could promote the cooperation in evolutionary games.There are still further works which we will perform. Considerations will be given to conduct simulations in heterogeneous structures where the initial distribution of players is not random. Specific real-world evolution situations need be considered in order to observe more practical rules for promoting the cooperation behavior.



---

*Source: 989024-2015-01-28.xml*