### Msg 0.3 error 011105 - thanks for

## Eco-evolutionary dynamics of multigames with mutations

^{ 1 }Sayantan Nag Chowdhury, Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Writing &#x; original draft,

^{ 2 }Prakash Chandra Mali, Validation, Writing &#x; review &#x; editing,

^{ 1 }Matja Perc, Funding acquisition, Validation, Visualization, Writing &#x; review &#x; editing,

^{ 3 ,}

^{ 4 ,}

^{ 5 ,}

^{ 6 }and Dibakar Ghosh, Project administration, Supervision, Validation, Visualization, Writing &#x; review &#x; editing

^{}

^{ 2 ,}

^{*}

### Sourav Roy

^{1} Department of Mathematics, Jadavpur University, Kolkata, West Bengal, India

Find articles by Sourav Roy

### Sayantan Nag Chowdhury

^{2} Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata, West Bengal, India

Find articles by Sayantan Nag Chowdhury

### Prakash Chandra Mali

^{1} Department of Mathematics, Jadavpur University, Kolkata, West Bengal, India

Find articles by Prakash Chandra Mali

### Matja Perc

^{3} Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia

^{4} Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan

^{5} Alma Mater Europaea, Maribor, Slovenia

^{6} Complexity Science Hub Vienna, Vienna, Austria

Find articles by Matja Perc

### Dibakar Ghosh

^{2} Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata, West Bengal, India

Find articles by Dibakar Ghosh

The Anh Han, Editor

Author informationArticle notesCopyright and License informationDisclaimer

^{1} Department of Mathematics, Jadavpur University, Kolkata, West Bengal, India

^{2} Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata, West Bengal, India

^{3} Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia

^{4} Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan

^{5} Alma Mater Europaea, Maribor, Slovenia

^{6} Complexity Science Hub Vienna, Vienna, Austria

Teesside University, UNITED KINGDOM

^{}Corresponding author.

**Competing Interests: **We have declared that no competing interests exist.

* E-mail: [email protected]

Received Jun 20; Accepted Jul

Copyright © Roy et al

This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

### Abstract

Most environments favor defection over cooperation due to natural selection. Nonetheless, the emergence of cooperation is omnipresent in many biological, social, and economic systems, quite contrary to the well-celebrated Darwinian theory of evolution. Much research has been devoted to better understanding how and why cooperation persists among self-interested individuals despite their competition for limited resources. Here we go beyond a single social dilemma since individuals usually encounter various social challenges. In particular, we propose and study a mathematical model incorporating both the prisoner&#x;s dilemma and the snowdrift game. We further extend this model by considering ecological signatures like mutation and selfless one-sided contribution of altruist free space. The nonlinear evolutionary dynamics that results from these upgrades offer a broader range of equilibrium outcomes, and it also often favors cooperation over defection. With the help of analytical and numerical calculations, our theoretical model sheds light on the mechanisms that maintain biodiversity, and it helps to explain the evolution of social order in human societies.

### 1 Introduction

Numerous interdisciplinary researchers have long sought a way to understand how the tremendous biodiversity among species persists in nature despite the significant differences between them in terms of competitive capability. Darwinian evolution [1] always challenges the emergence of spontaneous social cooperation, as cooperators act entail an inherent cost for displaying altruism [2]. Self-interested defectors generally exploit these cooperators and hinder the maintenance of cooperation. This trivially raises the question of why cooperators act selflessly if only the fittest succeeds and how numerous forms of cooperative behavior sustain in nature ranging from a low-level microbial biological system to a high-level complex social system. We resort to the evolutionary game dynamical interaction given by a 2 × 2 pay-off matrix to solve this riddle. These classical games are capable of analyzing the economic and strategic decisions of rational individuals [3]. An unprecedented spectrum of researchers have already focused overwhelmingly on the evolutionary game theory [4&#x;6] to seek out the mechanisms that sustain and promote cooperation [7&#x;14].

We construct a mathematical model on evolutionary multigames [15&#x;27] by adopting two simple two-player games, viz. the prisoner&#x;s dilemma [28] (PD) and snowdrift game [29, 30] (SD), as metaphors for cooperation between unrelated rational individuals. Our choice of combining two distinct games is motivated by the fact that people may face different social complexities rather than a single one. People take care of problems based on their own perceptions. This real-life behavioral heterogeneity among individuals inspires us to examine the effect of multigames for understanding social diversity. Our theoretical approach is another effective way to explore the impact of the interaction between different cultural backgrounds on large-scale social ordering [31].

While defection is the ultimate profitable strategy in the classical PD game, irrespective of the co-player&#x;s decision, a player&#x;s best strategy in the SD game depends on the co-player&#x;s decision. Thus, we select two different games with two contrasting outcomes. The stable coexistence of both cooperators and defectors is the expected consequence in the SD game, resulting in the persistence of cooperative behavior. In comparison, cheaters are always encouraged to exploit the cooperative individuals in the set-up of the PD game. In our constructed model, unrelated individuals have their probabilistic selection to choose which game they want to play. We also consider the mutation [32&#x;42] as an evolutionary mechanism that enables people to switch their respective strategies. We assume that it is hard to behave rationally under every circumstance. Hence, in our mathematical model, we incorporate a simplistic assumption that cooperators&#x; and defectors&#x; subpopulations can interchange their strategies with a fixed bidirectional mutation rate. To the best of our knowledge, the majority of the studies on multigames have been examined in the absence of mutation.

Apart from that, to reveal the interplay between ecological and evolutionary dynamics, we inspect the influence of free space in our mathematical model. In the context of complex systems, free space [43&#x;49] proves to be a promising factor for offering significant consequences on diverse collective dynamics. We assume that each individual&#x;s birth rate depends not only on the respective average payoffs but are also proportional to the available free space. This free space will give them reproductive opportunities. The population of unrelated players gains a payoff due to the interaction of the multigames. We treat these payoffs as reproductive successes. Individuals with a higher payoff leave more offspring if sufficient empty spaces are available and are able to outcompete less successful ones. Free space will also contribute to the individual&#x;s average payoffs. We consider free space as an ecological variable donating selflessly to all players&#x; fitness. However, free space never anticipates any benefit from others. This generous nature is surprising in a society of self-interested individuals. Still, free space provides reproductive opportunities to each player and loses their own identity to improve the fitness of other individuals. We consider such a selfless act of free space because helping others is a common practice among human beings and other animals. This tendency of unselfish concern for other people is a unique recipe for promoting cooperation and favors the survival of various species. The inclusion of such ecological dynamics into evolutionary multigames gives rise to a fascinating way to reveal the influence of eco-evolutionary dynamics [50&#x;60] on decision-making and the evolution of cooperation.

The concurrence of ecological changes and rational players&#x; evolution, depicted in the present article, prevails through the same time scale, rather than the typical assumption that ecological processes are much faster than evolutionary processes [61]. Our theoretical model uncovers how the combination of ecological dynamics and game dynamics is beneficial for maintaining cooperative behaviors under the influence of mutation and leads to the stable coexistence of interacting competitors. Earlier, Nag Chowdhury et al. [8, 16] examined the influence of eco-evolutionary dynamics on cooperative behavior&#x;s emergence and evolution under a cooperation-supporting mechanism, viz. punishment. Nevertheless, our model exhibits the coexistence of cooperators and defectors in the absence of any such supporting mechanism. Moreover, individuals&#x; birth rates in those Refs. [8, 16] depend only on their respective average payoffs. Apart from the addition of free space&#x;s altruistic behavior in the payoff matrix, we furthermore include the importance of available free space on the birth rates by assuming their birth rates to be proportional to the available space. The rest of the article is structured as follows: In Sec. (2), we discuss how we construct our mathematical model in detail. The following section (3) deals with the existence, uniqueness, and boundedness of the model, along with extensive numerical simulations and discussions. The section (4) presents the summary of our findings. Following that, we provide a brief discussion in the last section (5) on the challenges in this work&#x;s context that need to be addressed and are worth studying in the future.

### 2 The model

To start with, we consider a simplistic assumption that each individual has two distinct choices, viz. (i) cooperation (**C**) and (ii) defection (**D**). Even they can play any of the two possible games (a) PD game and (b) SD game. They can adopt the PD game with probability *p*, and alternatively, they interact with other individuals by playing the SD game with the complimentary probability (1 &#x; *p*). Both of these two-person games can be given by the following two payoff matrices *A* and *B* respectively, where

in which the entries portray the payoff accumulated by the players in the left.

Here, *R*_{PD} and *R*_{SD} contemplate the reward for mutual cooperation among two players in the respective PD and SD games. Similarly, both unrelated players receive the punishment *P*_{PD} and *P*_{SD} for mutual defection in the games PD and SD, respectively. An exploited cooperator gains the sucker&#x;s payoff *S*_{PD} and *S*_{SD}, respectively in the PD and SD games when confronted by a defector. The mixed choice yields the defector temptations *T*_{PD} and *T*_{SD} to exploit a cooperator in the PD and SD games, respectively. The payoff ranking of these four-game parameters determines the two-person games. The conventional relative ordering for the PD game is *T*_{PD} *R*_{PD} *P*_{PD} *S*_{PD} [8, 16, 31] and 2*R*_{PD} *S*_{PD} + *T*_{PD} [62]. Without loss of any generality, we choose *T*_{PD} = *²* 1, *R*_{PD} = 1, *P*_{PD} = *·* &#x; (0, 1), and *S*_{PD} = 0. Similarly, we choose *T*_{SD} = *²* 1, *R*_{SD} = 1, *S*_{SD} = 0, and *P*_{SD} = &#x;*·* &#x; (&#x;1, 0), maintaining the relative ordering *T*_{SD} *R*_{SD} *S*_{SD} *P*_{SD} for the standard SD game [8, 16, 31]. Thus, the payoff matrices *A* and *B* become

Note that, in both of these games, mutual cooperation leads to the payoff *R*_{PD} and *R*_{SD}, which is relatively higher than *P*_{PD} and *P*_{SD}, which one defector receives when playing with a defector. Thus, cooperation always promises higher income than defection if both the rational players choose the same strategy. The difference between these two games&#x; relative ordering leads to a contrasting scenario. In the SD game, the interaction between the cooperator and defector always promises a better income in terms of payoff than the interaction between two defectors. A reverse reflection is observed in the case of the PD game thanks to the choice of such relative ranking of game parameters in the PD game. The interaction between two defectors in the PD game allows them to earn more than a cooperator encountering a defector. The switching between *P* and *S* in the relative ordering of both games thus produces noticeable unexpected consequences on the evolution of cooperation.

Since a player can decide which game they want to play, thus the final payoff matrix for the multigames looks like

Here, *p* &#x; [0, 1] is the probability of playing the PD game. Moreover, we consider free space **F** as an ecological variable, which contributes altruistically by helping others. Nevertheless, free space does not get any benefits by giving them reproductive opportunities. We incorporate this charitable role of free space by extending the 2 × 2 payoff matrix *E* to the 3 × 3 payoff matrix *G* as follows

The matrix *G* clearly reveals that the free space never earns any payoff for their selfless charitable act; however, it contributes a positive payoff *Ã*_{1} and *Ã*_{2} to the cooperators and the defectors, respectively. Since most of the game parameters (not all) lies within the closed interval [0, 1], we assume, for the sake of feasible comparison, *Ã*_{1} and *Ã*_{2} both lies within the interval [0, 1]. When *Ã*_{1} and *Ã*_{2} are equal to zero, free space will not contribute anything to anyone. However, whenever *Ã*_{1} and *Ã*_{2} attain positive values, individuals gain an additional payoff from free space.

Inspired by the Malthusianism, we consider the following set of differential equations governing the changes in frequencies of cooperators and defectors as a function of time *t*

(1)

where

Here, *x* and *y* are the normalized densities of cooperators and defectors, respectively. Let *z* be the available free space. Thus, we have

Relation (2) assures that by studying *x* and *y* alone, one can easily capture the dynamics of the two strategies. We assume the birth rates of each individual depends crucially on the available free space as well as on their respective average fitness. Thus, we consider

(3)

where *f*_{C} and *f*_{D} are average fitness of the cooperators and the defectors, respectively. The average payoff of cooperators and defectors can be determined using the payoff matrix G of the multigames, and the relation (2) as follows,

(4)

Note that the average fitness of free space is

This is expected as free space does not gain anything for its benevolent nature. For simplicity, we further assume that all individuals die at a uniform and constant mortality rate *¾* &#x; (0, 1]. Hence using the relations (3) and (4), our constructed model (1) transforms into

(6)

Now, we introduce a constant probability *¼* as a rate with which each individual mutates from one strategy to the others, in a continuous manner,

(7)

Relation (7) reflects that the mutation probability *¼* from the cooperators to the defectors is identical to the mutation rate from the defectors to the cooperators. So using the system (6), the well-mixed population under the influence of bidirectional mutation gives rise to the differential equations

(8)

The system (8) contains seven different parameters. We summarize the necessary information about these parameters in (Table 1).

### Table 1

**Parameters with their physical significance and domains: The first column represents the symbols used in the present manuscript to describe several parameters.**

At the same time, the second column depicts what these parameters exactly mean. The set of possible input values from which the parameters can assume their values is given in the third column of the table. For further details, please see the main text.

Parameters | Physical Significance | Domain |
---|---|---|

¾ | Death Rate | (0, 1] |

² | The gain of a defector while interacting with a cooperator | (1, 2) |

· | The payoff for mutual defection | (0, 1) |

p | Probability of playing the PD game | [0, 1] |

¼ | Mutation probability | [0, 1] |

Ã_{1} | The altruistic incentive of free space towards the cooperators | [0, 1] |

Ã_{2} | The altruistic incentive of free space towards the defectors | [0, 1] |

Open in a separate window

### 3 Results and discussions

### Existence, uniqueness and positive invariance

Before investigating the model (8) using numerical simulations, we first prove the positive invariance of the proposed system (8). Clearly, the functions on the right-hand side of the differential Eq (8) are continuously differentiable and at the same time, locally Lipschitz in the first quadrant of . This ensures the existence and uniqueness of solutions for the model (8) with suitable initial values. Note that the initial conditions (*x*_{0}, *y*_{0}) must lie within the domain [0, 1] × [0, 1] maintaining the inequality 0 &#x; *x*_{0} + *y*_{0} &#x; 1, as they represent the frequencies of the two strategies. To determine the positivity of the proposed model (8), we write the system as follows,

(9)

where and are two integrable functions in the Riemannian sense. Solving Eq (9), we get

(10)

where *c*_{1} and *c*_{2} are the integrating constants depending on the initial densities *x*_{0} and *y*_{0}. This proves that both *x* and *y* are non-negative. Now, to calculate the upper bound of *x* + *y*, we proceed as follows. The dynamical Eq (8) yield

(11)

Since as per our previous analysis, *x* + *y* &#x; 0 and the mortality rate *¾* is always positive, we have *¾*(*x* + *y*)&#x;0. Thus, (11) reduces to

(12)

Integrating both sides, we get

(13)

where *c*_{3} is the initial density dependent constant. Hence, we find that

(14)

i.e., the overall species density *x* + *y* eventually remains bounded within the region [0, 1]. This boundedness within the closed interval [0, 1] allows us to relate the possible emerging dynamics of the system (8) to physically implementable scenarios with biological relevance. When the sum of the population density (*x* + *y*) is precisely one, the available free space is zero, as *z* = 1 &#x; (*x* + *y*). i.e., there is no reproductive opportunity accessible to any individual in that situation with *z* = 0. When (*x* + *y*) = 0, then individually *x* = 0 and *y* = 0. Hence, all individuals die, and *z* is equal to one. Thus despite the presence of ample free space, all the individuals are extinct under that circumstances.

### Coexistence of different stationary points depending on the initial conditions

Next, we point out the multistable dynamics of our model (8), resulting in the system&#x;s vulnerability to small perturbations. Initially, we set the parameter values at *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = , *Ã*_{2} = and *¼* = 0 in Fig 1(a). We vary the initial conditions (*x*_{0}, *y*_{0}) within the interval [0, 1] × [0, 1] maintaining the inequality (*x*_{0} + *y*_{0})&#x;[0, 1]. We find that the dynamics switching between two stationary points, viz. *E*_{0} = (0, 0) and *E*_{1} = (*x**, 0). We analytically calculate the stationary points of the system (8) in the absence of the mutation, i.e., with *¼* = 0. We trace out four different stationary points,

Open in a separate window

Fig 1

**Alternations between multiple co-existing steady states depending on initial conditions.**

The coexistence of various stable states is portrayed here by varying the initial conditions maintaining 0 &#x; (*x*_{0} + *y*_{0})&#x;1. Red points signify the extinction equilibrium *E*_{0}. Violet points represent the defector-free steady state *E*_{1}. The cooperator-free stationary points *E*_{2} are shown by sea blue points, and the coexistence equilibrium *E*_{3} is plotted using yellow points. The mutation-free model in subfigures (a-c) allows four stable steady states to coexist. However, subfigure (d) supports only bistability for the chosen parameters&#x; values. The parameter values for each of these subfigures are (a) *¾* = , *²* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. (b) *¾* = , *²* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. (c) *¾* = , *²* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. (d) *¾* = , *p* = , *²* = , *¼* = , *Ã*_{1} = , and *Ã*_{2} = Other parameter is *·* =

*E*_{0}= (0, 0) reveals all individuals die. This point is locally stable when*E*_{1}= (*x**, 0) exhibits a society free from any defectors. Here, and*x** &#x; (0, 1]. The stability criteria is given by*E*_{2}= (0,*y**) represents that we have only left with defectors. Here, and*y** &#x; (0, 1]. This stationary state is stable if*E*_{3}= (*x***,*y***) indicates the coexistence equilibrium offering the survival of cooperation and defection simultaneously. Here, , and*x*** satisfies the equation, , where, .*x*** and*y*** both should lie within the interval (0, 1). The local stability yields the conditions for the stability of (*x**,*y**) are

Clearly, the chosen parameters satisfy the local stability criteria of both stationary points *E*_{0} and *E*_{1}. Note that *E*_{1} leads to two different values for the selected values of the parameters, out of which (, 0) always remains locally stable, and (, 0) is unstable. Thus, we see the appearance of two different stationary points *E*_{0} (shown by red points) and *E*_{1} (displayed by violet points) in the basin of attraction portrayed in Fig 1(a). Such toggling between alternate stable stationary points is one of the generic features in some biological systems involving the fundamental processes of life [63&#x;65] and in a few nonlinear dynamical systems [66&#x;69]. The initial condition (*x*_{0}, *y*_{0}) = (0, 0) always helps the system (8) to converge to the stationary point *E*_{0}. The system is never able to give rise to the survival of any cooperators and defectors without the presence of any individual at the beginning. Since the chosen initial point (0, 0) itself is the stationary state, the system always stabilizes in the (0, 0) stationary point irrespective of the parameter values. Fig 1(b) is drawn with *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. Similarly, we find the system (8) in the absence of mutation (i.e., *¼* = 0) converges to *E*_{0} for the single initial condition (*x*_{0}, *y*_{0}) = (0, 0). It is anticipated that the choice of *x*_{0} = 0 always leads to the cooperator-free steady state. The initial absence of cooperators in the mutation-free model will not entertain any exposure for the cooperators in the long run. The line of initial conditions *x*_{0} = 0 and *y*_{0} &#x; 0 produces the stationary state *E*_{2} (sea blue points in Fig 1(b)). Apart from these two steady states, the mutation-free system also switches between *E*_{1} (violet points) and *E*_{3} (yellow points) depending on the suitable choice of initial conditions (See Fig 1(b)). Thus, for the same choices of parameters&#x; values, the system flips between four alternate steady states depending on the initial densities of cooperators and defectors. Similarly, we find the system (8) with the parameters&#x; values *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = , *Ã*_{2} = and *¼* = 0 converges to all these four stationary points. All population extincts for the initial conditions ranging from (0, 0) to (, 0) (See red points in Fig 1(c)). The system (8) can be solved analytically with *Ã*_{1} = *Ã*_{2} = *¼* = 0 in the absence of defectors (i.e., *y* = 0) as follows,

where *c*_{4} is the integrating constant.

Similarly, the system (8) can be solved analytically with *Ã*_{1} = *Ã*_{2} = *¼* = 0 in the absence of cooperators (i.e., *x* = 0) as follows

where *c*_{5} is the integrating constant.

Fig 1 is plotted by solving the differential Eq (8) by varying the initial conditions within [0, 1] × [0, 1] with fixed step-length *´x*_{0} = *´y*_{0} = and maintaining 0 &#x; *x*_{0} + *y*_{0} &#x; 1. To solve our proposed system (8) numerically, we use the 4th order Runge-Kutta (RK4) method with 20 × 10^{5} iterations with fixed integration step length *´t* = The final point (*x*, *y*) is stored to decide the asymptotic dynamics of the system. Other initial conditions with *y*_{0} = 0 stabilize the dynamics in the defector-free state *E*_{1} (violet points). We also track a fair portion of the basin in Fig 1(c), where the system (8) with *¼* = 0 converges to either *E*_{2} (sea blue points) or *E*_{3} (yellow points) depending on the choice of initial conditions.

Although we choose the mutation-free model for the Fig 1(a)&#x;1(c), we consider the contribution of *¼* in Fig 1(d). We generate Fig 1(d) with *¾* = , *p* = , *²* = , *·* = , *¼* = , *Ã*_{1} = , and *Ã*_{2} = This non-zero *¼* leads to the disappearance of two stationary points, (i) the defector-free steady state *E*_{1}, and (ii) the cooperator-free steady state *E*_{2}. The symmetric mutation from one species to another species always gives them two feasible opportunities. Either both strategies survive or all the individuals perish for *¼* &#x; 0. We also mathematically derive two possible stationary points as follows,

- The extinction equilibrium
*E*_{0}= (0, 0), which is stable under the conditions, - The interior equilibrium
*E*_{3}= (*x**,*y**) becomes stable, if

The compelling evidence of bistability under the same choice of parameters&#x; values is recognized in Fig 1(d).

### Emergent dynamics in absence of mutation

We examine the impact of the parameters *p*, *²*, *·*, and *¾* on the system (8) with *¼* = 0 in Fig 2. For the comparability, we choose *Ã*_{1} = *Ã*_{2} = 0 in subfigures (a), (c) and (e). Also, we set *Ã*_{1} &#x; 0 and *Ã*_{2} &#x; 0 for the subfigures (b), (d) and (f). In the absence of mutation (*¼* = 0) and free space induced benefits (*Ã*_{1} = *Ã*_{2} = 0), the defector-free is unable to stabilize. Although *E*_{1} exists for , but two eigenvalues of the Jacobian of the linearized system are

Open in a separate window

Fig 2

**The interplay of different parameters in the mutation-free model (8).**

*p* is the probability of playing the PD game. We investigate the role of *p* in our proposed model (8) in the absence of the mutation. The figures (a), (c), and (e) of the left-handed column are drawn with *Ã*_{1} = *Ã*_{2} = 0. The remaining figures in the right-hand column are plotted with non-zero values of *Ã*_{1} and *Ã*_{2}. The advantages provided by the free space help to promote cooperation, as portrayed through the broad area of yellow and violet regions in figures (b), (d), and (f). Red, sea blue, violet, and yellow depict the stable extinct state, cooperator-free state, defector-free state, and co-existence state, respectively. We also plot our calculated curves from the local stability analysis of the stationary points. The white line represents that of *E*_{0}, whereas the blue line portrays that of the coexistence equilibrium *E*_{3}. The brown dotted line is the stability curve for the cooperator-free stationary state *E*_{2}, while the blue dashed line represents that of the defector-free steady state *E*_{1}. The mismatches with the stability curve in a few places are due to the multistable behavior of our constructed model (8). Other parameter values are kept fixed at (a-d) *¾* = , (c-f) *²* = , (a,b,e,f) *·* = , (b,d,f) *Ã*_{1} = , (b,d,f) *Ã*_{2} = and *¼* = 0.

Since *²* 1 and *¾* 0, »_{2} is always positive. Consequently, *E*_{1} is always unstable. This result is physically meaningful, as both the chosen games (PD and SD) do not encourage a defection-free society without any supporting mechanism. However, we find three different steady states in Fig 2(a), 2(c) and 2(e). These stationary points and their corresponding stability analysis are given below,

The extinction equilibrium

*E*_{0}= (0, 0) (red points in subfigures (a), (c) and (e) of Fig 2) always exists and is always locally stable.*E*_{2}= (0,*y***) is the cooperator-free stationary point (sea blue points in subfigures (a), (c) and (e) of Fig 2), where . This state exists under the condition 0 &#x; 4*¾·*(2*p*&#x; 1)&#x;*·*^{2}(2*p*&#x; 1)^{2}, and becomes stable if*¾**·*(2*p*&#x; 1)(2*y*** &#x; 3*y***^{2}).*E*_{3}= (*x**,*y**) allows the coexistence of both cooperators and defectors (yellow points in subfigures (a), (c) and (e) of Fig 2), where , and , with . The stationary state is stable under the following conditions,

With increasing *p*, players are prone to play the PD game. Hence, we can not anticipate the survival of cooperators for large *p*. Thus, for large *p*, either the extinction equilibrium *E*_{0} or the cooperator-free stationary point *E*_{2} stabilizes in Fig 2(a), 2(c) and 2(e). However, the free-space induced benefits in Fig 2(b), 2(d) and 2(f) facilitate the emergence of cooperation and stabilize the interior equilibrium *E*_{3}. Even, Fig 2(f) provides a range in the *p* &#x; *¾* parameter space, where the defector-free stationary point *E*_{1} stabilizes. The comparative study between the left and right column of Fig 2 ensures the encouraging role of free space in the promotion of cooperation under suitable circumstances.

With increasing *²*, the defectors are getting extra aid, and thus, we get the stabilization of the extinction equilibrium *E*_{0} (red zone) in Fig 2(a). However, the presence of free space assistance converts that area into a coexistence zone (yellow zone) in Fig 2(b). The increment of *p* in both subfigures (a-b) permits only the existence of defectors, as people are playing mostly the PD game in such circumstances. Defectors are getting a favorable atmosphere in the PD game for our chosen PD game parameters&#x; values. A similar sort of stabilization of cooperator-free steady state (sea blue area) is observed in subfigures (c-d) for larger values of *p*. However, lower values of *p* provide the opportunity for playing the SD game. Hence, coexistence equilibrium (yellow region) stabilizes in Fig 2(c) and 2(d) for smaller *p*. Nevertheless, the inclusion of free space induced benefits helps to broaden the region of coexistence in the *p* &#x; *·* parameter plane, as shown in Fig 2(d). The comparison between the Fig 2(e) and 2(f) suggests that appropriate non-zero values of *Ã*_{1} and *Ã*_{2} entertain the stabilization of a defector-free steady state (violet region). Hence, we can trace a fair portion of violet points in the *p* &#x; *¾* parameter space of the subfigure (f). However, too large a mortality rate reduces the opportunity of survivability of any individual, resulting in the stabilization of the extinction equilibrium *E*_{0}. We focus on the effect of mortality rate *¾* more elaborately in Fig 3.

Open in a separate window

Fig 3

**Influence of death rate**

*¾*in the mutation-free model (8).The influence of game parameters *²* and *·* and the mortality rate *¾* on the emergent dynamics of the model (8) in the absence of mutation is examined here. Subfigures (a) and (c) are drawn with *Ã*_{1} = *Ã*_{2} = 0. While subfigures (b) and (d) are generated with *Ã*_{1} = and *Ã*_{2} = Since, *Ã*_{1} *Ã*_{2}, we notice a fair portion of defector-free region (violet points) in subfigures (b) and (d). In the absence of free space induced benefits, we are unable to trace such defector-free regions in subfigures (a) and (c). The red region reflects the disappearance of all individuals for relatively high values of death rate *¾*. Note that we keep the value of *p* fixed at Thus, the system gets more opportunities to play the SD game, which generally encourages stabilizing the coexistence equilibrium. Therefore, we notice lower values of *¾* will lead to the convergence towards the coexistence equilibrium *E*_{3} (yellow points). *²* is kept fixed at for the subfigures (a-b) and *·* is set at for subfigures (c-d). We iterate the system (8) with *¼* = 0 for 20 × 10^{5} iterations with fixed integration step size *´t* = and fixed initial condition (*x*_{0}, *y*_{0}) = (, ).

As expected, higher values of *¾* constantly enlarge the chances of extinction. Thus, we observe a fair portion of the red region in Fig 3. Nevertheless, the amount of this red area is considerably lesser in the right column of Fig 3 compared to the left column. We introduce the non-zero values of *Ã*_{1} and *Ã*_{2} in the right column of this figure. These free space induced benefits encourage maintaining a defector-free society, as shown in Fig 3(b) and 3(d). We choose *Ã*_{1} = = *Ã*_{2} for Fig 3(b) and 3(d). i.e., the free space will provide an additional advantage for the cooperators, and thus depending on the other parameters, the defectors are vanished in the long run, as shown in Fig 3(b) and 3(d). In all of these subfigures of Fig 3, we trace a portion of yellow points depicting the survival of both cooperators as well as defectors simultaneously. We choose *p* = in Fig 3. Thus, people get more chances to play the SD game, which facilitates the concurrence of both cooperation and defection. Thus, smaller values of *¾* provide an opportunity to coexist for all strategies, which is observed in Fig 3 with initial condition (*x*_{0}, *y*_{0}) = (, ). We plot the stability curve of all stationary points in Figs 2&#x;5 as follows,

Open in a separate window

Fig 5

**Investigating the impact of**

*Ã*_{2}on the evolution of cooperation.For the smaller values of *Ã*_{2}, the mutation-free system gives rise to the defector-free (violet region) society. While larger values of *Ã*_{2} encourage the defectors and stabilize the cooperator-free stationary point *E*_{2} (sea blue region). Suitable choices of parameters may favor the co-existence of cooperators as well as defectors. Even for a higher mortality rate *¾*&#x;s values, all individuals die in the red portion of the parameter space *Ã*_{2} &#x; *¾* in subfigure (d). All the parameters are set at *¾* = , *·* = , *²* = , *p* = , *Ã*_{1} = , and *¼* = 0, unless it is varied. The simulations are computed for 20 × 10^{5} iterations with fixed initial condition (*x*_{0}, *y*_{0}) = (, ). The color codes are as follows: (i) red represents the extinct state *E*_{0}, (ii) yellow represents the co-existence state *E*_{3}, (iii) violet represents the defector-free state *E*_{1}, and (iv) sea blue represents cooperator free state *E*_{2}, respectively. We plot the stability curves corresponding to each stationary point in each subfigure.

(i) The brown dotted line for the cooperator-free steady state

*E*_{2},(ii) the blue dashed line for the defector-free steady state

*E*_{1},(iii) the solid blue line for the interior equilibrium

*E*_{3}, and(iv) the solid white line for the extinction equilibrium

*E*_{0}.

In Fig 4, we fix the value of *Ã*_{2} at and examine the role of *Ã*_{1}. Since *Ã*_{1} represents the free space induced benefits towards the cooperators, thus enhancement of its (*Ã*_{1}) value will help in promoting cooperation in the society. Fig 4 reflects the same scenario. A defector-free society (violet region) is noticed in all these subfigures. The increment of temptation parameter *²*&#x;s value challenges the prevalence of cooperation. Thus, for a small *Ã*_{1}, we find a defector dominated society (the sea blue region in Fig 4(a)). *Ã*_{1} proves to be a cooperator facilitating parameter as we detect a wide range of yellow regions in Fig 4(a), where both cooperators and defectors can coexist. Similarly, irrespective of the choice of *·* in Fig 4(b), higher values of *Ã*_{1} provide all cooperators an extra benefit for survival, and thus, the stationary point *E*_{1} (violet points) stabilizes. For the intermediate choice of *Ã*_{1}, the stationary point *E*_{3} (yellow) yields the stable coexistence of both strategies. However, cooperators strive to keep in existence for the lower values of *Ã*_{1}, and we find the sea blue region of cooperator-free steady state. The parameter *p* indicates the probability of playing the PD game. Thus, *p* &#x; 1&#x; always gives defectors a more favorable environment to survive. That&#x;s why we spot a sea blue portion in Fig 4(c) for higher values of *p* and *Ã*_{1}. When *p* is small, people are prone to play the SD game; and thus, we get the coexistence of both strategies (the yellow region in Fig 4(c)) for smaller *p* and *Ã*_{1}. Nevertheless, a larger value of *Ã*_{1} with a moderate value of *p* always provides a reasonable scope for the cooperators to survive, and we discover a healthy portion of the stationary point *E*_{1} (the violet region) in Fig 4(c). A higher mortality rate *¾* always results in the extinction of both cooperators and defectors, and we locate a huge red region in the *Ã*_{1} &#x; *¾* parameter plane of Fig 4(d). We find a very tiny sea blue region in Fig 4(d), where defectors are only able to survive. But, as expected, a higher value of *Ã*_{1} always promotes the cooperation strategy, and we obtain a violet zone of defector-free stationary point and a yellow region of interior equilibrium *E*_{3} in Fig 4(d).

Open in a separate window

Fig 4

**Importance of**

*Ã*_{1}in the enhancement of cooperation.Four different parameter spaces (a) *Ã*_{1} &#x; *²*, (b) *Ã*_{1} &#x; *·*, (c) *Ã*_{1} &#x; *p*, and (d) *Ã*_{1} &#x; *¾* are contemplated here with fixed initial condition *x*_{0} = and *y*_{0} = *Ã*_{1} is varied within [0, 1], and the other parameters&#x; values are for the subfigures (a) *¾* = , *·* = , *p* = , *Ã*_{2} = , and *¼* = 0. *²* is varied within the open interval (1, 2), (b) *¾* = , *²* = , *p* = , *Ã*_{2} = and *¼* = 0. *·* is varied within (0, 1), (c) *¾* = , *²* = , *·* = , *Ã*_{2} = , and *¼* = 0. *p* is varied within the closed interval [0, 1], and (d) *²* = , *·* = , *p* = , *Ã*_{2} = and *¼* = 0. *¾* is varied within the interval (0,1]. The color code represents the following: (i) red represents the extinct state, (ii) yellow portrays the co-existence state, (iii) violet displays the defector-free state, and (iv) sea blue depicts cooperator-free state, respectively. We have run the numerical simulations for 20 × 10^{5} iterations for each point and store the final value for determining the final asymptotic state. Increment of *Ã*_{1} contributes more to the cooperators&#x; payoff, and hence, we observe a defector-free region for higher values of *Ã*_{1} depending on the other parameters. We draw the stability curves for *E*_{0} (the solid white line), *E*_{1} (the blue dashed line), *E*_{2} (the brown dotted line), and *E*_{3} (the solid blue line).

Fig 5(a) shows that except for a smaller portion of the defector-free region (violet zone), the whole *Ã*_{2} &#x; *²* parameter space produces the coexistence of both cooperators and defectors depending on other parameters&#x; values. Despite the increment of *Ã*_{2} and *²*, the cooperators are able to survive along with the defectors due to our choice of other parameters&#x; values. The *Ã*_{2} &#x; *·* parameter space portrays that smaller values of *Ã*_{2} can not provide any benefit to the defectors, and stabilize the defector-free stationary point *E*_{1} (violet region) irrespective of *·*&#x;s value. However, if *Ã*_{2} increases, it will yield a window of opportunity for the defectors to thrive. We observe a wide coexistence region (yellow region) and a small area of cooperator-free stationary point (sea blue region)in Fig 5(b). *Ã*_{2} indicates the free-space induced benefits towards the defectors. So, it is expected that larger values of *Ã*_{2} always enhance the chances of defectors&#x; survivability. Thus, we notice a cooperator-free region (sea blue zone) in both Fig 5(c) and 5(d). Nevertheless, larger values of *p* enhance the chances of playing the PD game, where defectors get a favorable environment to survive. Thus, the cooperator-free sea blue region is found in Fig 5(c). We are also able to detect a small violet region of a defector-free environment in the *Ã*_{2} &#x; *p* parameter space. However, we trace a healthy portion of coexistence (yellow) too in Fig 5(c) due to our chosen parameters&#x; values.

### The influence of bidirectional mutation

Now we investigate the influence of bidirectional mutation on the long-term behavior of the nonlinear differential Eq (8). Since every species can mutate into the other at a specific uniform rate *¼* &#x; (0, 1], thus cooperators and defectors can not remain alive alone. Either all populations will die; otherwise, the dynamics will lead to the coexistence of all two species. Fig 4 already portrays the influence of free space-induced benefits on the cooperators in the absence of mutation. We scrutinize the impact of *Ã*_{1} under the influence of 50% mutation (i.e., *¼* = ) in Fig 6. As *Ã*_{1} increases, the cooperators are getting a better environment for survival. We find a portion of stable coexistence equilibrium in each two-dimensional parameter space for large *Ã*_{1} in Fig 6. Thanks to the mutation, the cooperators can not live alone. The white line in Fig 6 is the stability curve corresponding to the extinction equilibrium. The local stability analysis fits almost exactly with the numerical simulations in Fig 6 with fixed initial condition (*x*_{0}, *y*_{0}) = (, ). There are a few places where the stability analysis fails to predict the stabilization of the extinction stationary point (0, 0). This is mainly due to the multistable behavior of the proposed model 8.

Open in a separate window

Fig 6

**The effect of altruist free space on the nonlinear dynamics of multigames with mutation.**

The influence of free space-induced benefits by varying *Ã*_{1} within the closed interval [0, 1] is established. The whole population goes extinct in the red region, and the yellow area reflects the system&#x;s stable interior point, corresponding to the coexistence of all two strategies. All the parameters are kept fixed at *¾* = , *·* = , *p* = , *Ã*_{2} = , *²* = , *¼* = , unless it is varied. The initial condition is kept fixed at (, ). A notable difference is observed from (Fig 4), which is drawn in the absence of mutation. Larger values of *Ã*_{1} always facilitate the maintenance of cooperation, and the bidirectional mutation reinforces the system&#x;s inherent tendency to flow from cooperators to defectors. The solid white line in each figure is the analytically computed stability curve corresponding to the extinction equilibrium.

As discussed, the increasing values of the parameters *²*, *·*, and *p* always provide the defectors a favorable environment to dominate the cooperators. However, beyond a critical value of *Ã*_{1}, both the cooperators and defectors can coexist, as shown in Fig 6(a)&#x;6(c). The larger values of *¾* always hinder the evolution of cooperators as well as of defectors. However, an intermediate choice of *Ã*_{1} &#x; *¾* favors the successful evolution of both strategies, as portrayed through Fig 6(d). The same feature is also noticeable in Fig 7. The complex evolutionary dynamics switch between two stationary points depending on the choices of parameters&#x; values in Fig 7. For smaller values of *Ã*_{2}, the defectors are not getting enough advantages to survive, and thus the extinction equilibrium (0, 0) stabilizes in the red region of all subfigures of Fig 7. The choice of other parameters&#x; values is also crucial for obtaining these stationary states. The parameter *Ã*_{2} benefits the defectors; hence the defectors can survive beyond a certain threshold of *Ã*_{2}. The employed 50% mutation rate helps to flow a certain fraction of defectors into cooperators, and we have a moderate portion of coexistence state (yellow region) in Fig 7. The observed results may vary for different choices of initial conditions, as the system is multistable. We plot the boundary separating solid white lines in all subfigures by analyzing the local stability analysis of the extinction equilibrium. Clearly, this stability curve agrees well with our numerical simulations, and the places, where they don&#x;t agree with the numerical simulations, is due to the multistable behavior of our proposed model 8. All the simulations are done by iterating for 20 × 10^{5} times with fixed integrating step length *´t* = The last point is gathered to finalize the asymptotic state. All the codes to generate these figures are publicly available at Ref. [70].

Open in a separate window

Fig 7

**The impact of free space-induced benefits on the defectors in the presence of the mutation.**

Either all populations die (red region), or the symmetric mutation preserves the coexistence state (yellow area) in all these subfigures. Beyond a critical value of *Ã*_{2}, both strategies mutate one into another, providing the opportunities for the coexistence of both species. The solid white line is the analytically derived stability curve to stabilize the extinction equilibrium. All the parameters are kept fixed at *¾* = , *·* = , *p* = , *Ã*_{1} = , *¼* = , *²* = , unless they are varied. The initial condition is kept fixed at (, ). The slight mismatch of the stability of the extinction equilibrium (0, 0) in subfigures (a) and (c) is due to the multistable dynamics of our proposed model (8).

Fig 8 demonstrates the importance of the parameters *Ã*_{1}, *Ã*_{2} and *¾* for the enhancement of cooperation under the presence of mutation. The larger values of *Ã*_{1} and *Ã*_{2} facilitate the evolution of at least one species, and the positive mutation rate *¼* &#x; (0, 1] assures that species should mutate into the other. This mechanism will lead to the species&#x; coexistence in a major portion (yellow region) of the two-dimensional parameter spaces represented in Fig 8(a) and 8(b). The results are further validated by plotting the stability curve (white solid lines) below which the stationary point (0, 0) is locally stable. A higher mortality rate never entertains the evolution of both strategies; thus, we obtain the red region in Fig 8(c). This red region indicates the extinction equilibrium (0, 0). Once again, we plot the stability curve (solid white line) in Fig 8(c), above which both the species should be extinct as per our local stability analysis. All the subfigures are drawn with fixed initial condition (, ).

Open in a separate window

Fig 8

**Effect of bidirectional mutation on the emergent dynamics.**

The positive mutation rate *¼* &#x; (0, 1] allows the system (8) to settle into two stationary points depending on the choice of other parameters&#x; values. When the free space-induced benefits are small, all species are extinct, as depicted through subfigures (a-b). However, all strategies can coexist for suitable choices of other parameters. The red region reflects the extinction equilibrium, and the yellow region indicates the stable coexistence of cooperators and defectors. The solid white line represents the analytically derived stability curve below which the extinction equilibrium (0, 0) is locally stable in subfigures (a-b). The initial condition is chosen here as (*x*_{0}, *y*_{0}) = (, ). Other parameters&#x; values are kept fixed at *¾* = , *²* = , *·* = , *p* = , *Ã*_{2} = , *Ã*_{1} = , unless they are varied. The moderate choice of mortality rate *¾* &#x; (0, 1] allows the species&#x; coexistence in subfigure (c). Beyond a critical value of *¾*, both the cooperators as well as the defectors die, as reflected through the red region of subfigure (c).

We inspect the influence of different parameters on the constructed model (8) in Fig 9. We keep fixed all parameters&#x; values at *¾* = , *²* = , *¼* = , *p* = , *·* = , *Ã*_{1} = and *Ã*_{2} = , unless they are varied. Since free space provides additional advantage to the cooperators compared to the defectors as we choose *Ã*_{1} = *Ã*_{2} = , we have a defector-free society at least initially with *¼* = 0 in Fig 9(a). However, whenever the mutation rate *¼* becomes positive, each strategy can mutate into other. In this way, the density of the cooperators (red line) decreases, and the defectors&#x; density (blue line) increases. Eventually, both densities almost become identical for *¼* &#x; 1 &#x;. In of Fig 9(b), the vital role of *p* is investigated in the one-dimensional bifurcation diagram. As *p* &#x; 1 &#x;, the defectors are getting the upper hand over the cooperators as *p* indicates the probability of playing the PD game. PD game always provides additional assistance to the defectors. All these results obtained in Fig 9 are consistent with the social dilemmas considered for constructing the model (8). Fig 9(b) points out that lower values of *p* are always better for the maintenance of cooperation as small values of *p* indicate more rational people are playing the SD game and the SD game always favors the coexistence of both strategies. Increasing the temptation parameter *²* always uplifts the defectors&#x; fraction *y* (blue line). Thus, one needs to choose the value of *²* wisely so that we obtain a moderate range of *²* in Fig 9(c) allowing the coexistence of both strategies. When *²* is small, we found the cooperators&#x; density *x* (red line) dominates the defectors&#x; fraction *y* (blue line). However for *²* , *y* is larger than *x*. Note that the results may alter for a different choice of initial condition as the system 8 is multistable. We plot Fig 9 with fixed initial condition (, ) and the code to generate this figure is freely available at [70]. Increasing the parameter *·* &#x; (0, 1) can provide extra benefits to the defectors. Thus, the rate of increment of *y* is slightly better than that of *x* in subfigure (d) of Fig 9. However, other parameters&#x; values are also crucial for stabilizing the competitive communities. Hence, *y* (blue line) remains always lower than *x* (red line) in Fig 9(d) for our chosen parameters&#x; values.

Open in a separate window

Fig 9

**Effects of various parameters in spreading of cooperative behavior.**

Subfigure (a) reveals two contrasting scenarios, where the cooperators decrease with the increment of the mutation rate *¼*. However, the defectors&#x; fraction gradually increases with increasing *¼*. In the absence of mutation (i.e., for *¼* = 0), all defectors are extinct as the parameters values chosen to produce this subfigure is *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = and *Ã*_{2} = Thus, the free space induced benefit towards cooperator *Ã*_{1} is greater than that of defector *Ã*_{2}. Subfigure (b) shows that both species&#x; densities gradually increase with the probability *p* of playing the PD game. As the probability of playing the PD game increases, the increment of the fraction of defectors (blue line) gets over in margin with respect to the fraction of the cooperators (red line) after *p* = depending on other parameters&#x; values. This fact is obvious from the view that a higher possibility of playing the PD game gets a more beneficial ambiance for defectors. Subfigure (c) indicates the same phenomena as the value of the temptation payoff parameter *²* increases, both the fractions escalate. Whereas at *²* = , both these fractions get the same value. After that value of *²*, the value of the fraction of defector (blue line) becomes more than that of the cooperators (red line). Subfigure (d) also highlights the increment of both the fractions as the punishment payoff parameter *·* increases. But up to the highest received payoff value of *·*, the value of the fraction of defectors can never exceed the cooperators&#x; fraction because the free space provides more benefit to the cooperators than the defectors. Subfigure (e) demonstrates that when *Ã*_{1} &#x; [0, ], both these fractions are extinct, as the mortality rate *¾* has more value than these free space-induced benefits. But as the value of *Ã*_{1} increases, the cooperator gets more advantages than defectors. Subfigure (f) shows that both the fractions increase with the enhancement of *Ã*_{2}. Up to *Ã*_{2} = , the cooperators&#x; density acquires more value than *y*, as the value of *Ã*_{1} = is taken sufficiently high. After *Ã*_{2} = , *y* becomes higher than *x*. All the subfigures are drawn by iterating the system (8) with 20 × 10^{5} times with fixed integrating step length *´t* = and fixed initial condition (*x*_{0}, *y*_{0}) = (, ). The last point at 20 × 10^{5} iteration is collected to identify the asymptotic state of the proposed model. Other parameters are kept fixed at *¾* = , *²* = , *¼* = , *p* = , *·* = , *Ã*_{1} = and *Ã*_{2} = , unless they are varied.

Since the mutation rate *¼* is 50%, we can obtain only two stationary points of the model (8). Thus, initially, for a smaller choice of *Ã*_{1}, extinction of both species prevails in Fig 9(e). Nevertheless, with an increment of *Ã*_{1}, cooperators are getting further assistance from free space, and hence *x* becomes positive beyond *Ã*_{1} = Since *¼* = , 50% of these cooperators mutate bidirectionally into the defectors, and thus, we have positive *y* too in that range of *Ã*_{1}. Since *Ã*_{2} = *¾* = , thus we get the stabilization of (0, 0) within the interval *Ã*_{1} &#x; [0, ]. In fact, the cooperators&#x; density *x* is slightly better than that of *y* for the higher values of *Ã*_{1} in Fig 9(e). The choice of the initial condition and other parameters&#x; values are vital in obtaining all these results. Similar observation can be found in Fig 9(f), where we examine the role of *Ã*_{2} in our proposed model (8). With increasing *Ã*_{2}, the defectors will be benefited from free space, and hence, *y* (blue line) dominates *x* (red line) for larger values of *Ã*_{2}. But since *Ã*_{1} is taken sufficiently large in this Fig 9(f), we initially have a small portion for smaller values of *Ã*_{2} where *x* *y*. Fig 9 discloses suitable choices of all parameters&#x; values not only entertain the coexistence of both strategies but also may promote the evolution of cooperation.

### 4 Concluding remarks

Our mutation-induced model provides widespread coexistence of both strategies under suitable choices of parameters&#x; values. Such stable persistence of ecological communities strengthens the theory of concurrency. Our thorough analysis with several numerical simulations enhances our understanding of the mechanisms that drive the survival of cooperative behavior in the multigames consisting of both the Prisoner&#x;s Dilemma and the Snowdrift games. The inclusion of altruistic behavior of free space along with the mutation in our proposed evolutionary model seems to be a natural course of action as observed in many realistic settings. Our findings clearly demonstrate that there exists an optimal probability of playing each game so that people are more likely to cooperate in such circumstances. Throughout the study, we have pointed out the positive impact of diverse factors (parameters) on significant improvement of cooperation. We have shown that the selfless contribution of free space promotes the coexistence of all strategies efficiently, even in the absence of mutation.

In summary, our research indicates the proposed model (8) may possess four different stationary points in the absence of mutation. The exciting feature of this model is that the system never allows settling into a cooperator-free state in the lack of free space-induced benefits and mutation rate. This precise result is consistent with the chosen games as both the PD and SD games never encourage a cooperator-free society in the usual scenario. The merging of two games with different outcomes provides a more realistic representation of the concept of opinion formation. However, the numerical results presented here are highly sensitive to the variation of initial conditions. The positive invariance and boundedness of the model are analyzed too in this article. We have shown the viable choice of bidirectional mutation allows the system to switch between only two stationary states. Either all people will die, or both the strategies coexist in the eco-evolutionary model with mutation. This is an interesting angle of our research as our model brings forth stable biodiversity in the form of a heterogeneous population (mixed cooperator-defector state). Note that we have only focused on the equilibria of the corresponding dynamical system (8) throughout the article so that we can relate those stationary states from the game-theoretical point of view. Our simple model with mutation sheds light on how cooperation emerges in a complex society. The presented insightful results attest that altruistic behavior and mutation are advantageous for the spontaneous maintenance of biodiversity. We believe the presented results may help us understand the mechanism behind the coexistence of competing species through the co-evolution of both strategies.

### 5 Discussion: Limitation &#x; future perspectives

Lastly, we discuss challenges related to our work that merits further investigation. Even though it is almost impossible to incorporate all complex, realistic relationships among social creatures using minimal modeling, it is nevertheless absolutely essential to investigate such models using elementary mathematical principles. In fact, the literature already provides a number of excellent models that offer a stimulating starting point for exploring many practical situations. Chen et al. [71] proposed an elegant model aimed at understanding the co-evolutionary outcome of the strategies of the traffic management department, drivers, and pedestrians. Reference [72] proposes a hybrid machine learning model to predict traffic accidents. A novel autonomous system without any stationary points may offer hidden attractors, as revealed in Ref. [73]. It is almost impossible to incorporate even a partial list of relevant references to emphasize the importance of minimal modeling here. Motivated by all these works, we present a new model that can offer several valuable insights into the evolution of cooperation and altruism. However, our model is far from perfect. Instead of considering two-person interactions, the introduction of group interactions may yield a deeper understanding of decision-making. Interdependent networks may be an excellent choice to study the evolution of cooperation, as reflected through Refs. [20, 31, 62, 74]. This remains a promising future research generalization, including interdependent networks which may yield several complex dynamical behaviors other than stable equilibrium states.

One should observe that our approach works absolutely fine with other social dilemmas. For any two two-person games of the form , the system constructed with the help of our policy always remains bounded within the closed interval [0, 1]. However, the stationary points alter due to the change in the mathematical model. We have verified this boundedness by considering the snowdrift and the harmony game [75, 76] (results are not shown here). In fact, on the generality of our results, we can comment that a similar impact of mutation and free space can also be expected for other two-person games. It will be interesting to examine how the results may vary with other games like public goods games, rock paper scissor, etc. Besides introducing other realistic scenarios like apology and forgiveness [77], intention recognition [78], delay [79] etc., one may advance our understanding of the roots of cooperation in social and biological systems. Moreover, the impact of higher-order interactions [80] among the agents on the coevolution [81] of cooperation and synchronization [82&#x;84] in a coupled network remains largely unexplored. The interdisciplinary researchers of complex systems can pay attention to this exciting topic for future research.

### Acknowledgments

The authors gratefully acknowledge the anonymous referees for their careful reading, insightful suggestions, and valuable remarks that helped considerably improve and readability of the manuscript. S.R. wants to convey his sincere gratitude to Gourab Kumar Sar, Md Sayeed Anwar, Srilena Kundu, and Subrata Ghosh of the Indian Statistical Institute, Kolkata, for their helpful discussion and suggestions.

### Funding Statement

M.P. was supported by the Javna Agencija za Raziskovalno Dejavnost RS (Grant Nos. P and J).

### Data Availability

All relevant data are within the paper.

### References

1. Darwin C. On the origin of species: A facsimile of the first edition. Harvard University Press; [Google Scholar]

2. West SA, Griffin AS, Gardner A. Social semantics: altruism, cooperation, mutualism, strong reciprocity and group selection. Journal of Evolutionary Biology. ;20(2)&#x; doi: /jx [PubMed] [CrossRef] [Google Scholar]

3. Weibull JW. Evolutionary game theory. MIT press; [Google Scholar]

4. Nowak MA. Evolutionary dynamics: exploring the equations of life. Harvard university press; [Google Scholar]

5. Sigmund K. The calculus of selfishness. Princeton University Press; [Google Scholar]

6. Hofbauer J, Sigmund K, et al. Evolutionary games and population dynamics. Cambridge university press; [Google Scholar]

7. Allen B, Lippner G, Chen YT, Fotouhi B, Momeni N, Yau ST, et al. Evolutionary dynamics on any population structure. Nature. ;()&#x; doi: /nature [PubMed] [CrossRef] [Google Scholar]

8. Nag Chowdhury S, Kundu S, Banerjee J, Perc M, Ghosh D. Eco-evolutionary dynamics of cooperation in the presence of policing. Journal of Theoretical Biology. ; doi: /sprers.eu [PubMed] [CrossRef] [Google Scholar]

9. Pennisi E. How did cooperative behavior evolve?Science. ;()&#x; doi: /science [PubMed] [CrossRef] [Google Scholar]

Fotouhi B, Momeni N, Allen B, Nowak MA. Evolution of cooperation on large networks with community structure. Journal of the Royal Society Interface. ;16() doi: /rsif

initMPS with licenseOption 0

mental ray: got 8 satellite CPUs.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

MSG error cannot start thread Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

## Importing Objects - mental ray error

Thanks Joojaa,

Not sure how to reset global render settings..

The architectural file is only 7mb, and renders fine in mental ray, although after importing a car that is only mb, I get these error messages and a partial render in mental ray, sayin 'out of memory'. This pc has 4gb of RAM. and the car has no attached map files.

Somethins kooky..

mental ray for Maya

mental ray: version , 22 Nov , revision

mental ray: mental ray for Maya - startup done

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

API warn type conflict: "input" is type vector, "maya_ramp" returns type color

API warn type conflict: "specularColor" is type color, "maya_multdiv" returns type vector

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

API warn type conflict: "input" is type vector, "maya_ramp" returns type color

API warn type conflict: "specularColor" is type color, "maya_multdiv" returns type vector

API warn type conflict: "specularColor" is type color, "maya_multdiv" returns type vector

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

MSG error cannot start thread 6: Not enough storage is available to process this command.

MSG error cannot start thread 5: Not enough storage is available to process this command.

MSG error cannot start thread 7: Not enough storage is available to process this command.

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

MEM fatal can't allocate bytes.

MEM fatal can't allocate bytes.

If i try the maya renderer, it locks up and then gives this below

Rendering using 2 threads.

Automatic near/far clipping values: ,

Starting Rendering C:/Users/David/Documents/maya/projects/default/images/tmp/Genesis iff.

Constructing shading groups.

Rendering current frame.

Failed to allocate space for image.

====================================

Cause of memory exception

====================================

Mb Free Memory

Mb Free Swap

Mb Size of alloc

Mb Low Memory Threshold

====================================

Memory use when exception was thrown

====================================

====================================

Page faults

Mb Max resident size

Mb Peak total size(Estimated)

Mb Peak arena size

====================================

Mb Current

Mb POLY_DRAW_CACHE_DATA

Mb Pixel Map

Mb Render Geometry Arena

Mb arguments

Mb AG

Mb MEL

Mb File Texture Mipmaps

Mb Render Cache

Mb Arrays

Mb Transforms

Mb Data Blocks

Mb Ray Tracer

Mb NURBS Surface Shapes

Mb NURBS AG

Mb Nurbs Tesselator Vertex

Mb POLY_DRAW_CACHE_STATIC_DATA

Mb NURBS Geometry Cache

====================================

I went back to an earlier iteration without the animation and tried merging in my character again, but Max crashed instantly, and was unable to autosave the file, siting a 'Memory Error'.

I did a quick search and the only other posting with similar a question is in the Maya forum.

Is the problem simply that my computer is not powerful enough?

Edit:

When I open the Material Editor now the following messages:

SCEN warn switching ambient occlusion off, since trace is off

MSG error cannot start thread 4: Not enough storage is available to process this command.

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

Also, I've just noticed that after adding the animation this file has gone from 9, KB to 14, KB

Max , Operating System: Windows 7 bit

Video Card: 1 Gigabyte ATI Radeon HD

CPU: Intel Core 2 Duo E GHZ, Ram: 6GB

Hard Drive: 1 Terrabyte Western Digital 32mb cache

### Msg 0.3 error 011105 - think

initMPS with licenseOption 0

mental ray: got 8 satellite CPUs.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

MSG error cannot start thread Für diesen Befehl ist nicht genügend Speicher verfügbar.

DISP error failed to memory map virtual frame buffer: Für diesen Befehl ist nicht genügend Speicher verfügbar.

## Eco-evolutionary dynamics of multigames with mutations

^{ 1 }Sayantan Nag Chowdhury, Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Writing &#x; original draft,

^{ 2 }Prakash Chandra Mali, Validation, Writing &#x; review &#x; editing,

^{ 1 }Matja Perc, Funding acquisition, Validation, Visualization, Writing &#x; review &#x; editing,

^{ 3 ,}

^{ 4 ,}

^{ 5 ,}

^{ 6 }and Dibakar Ghosh, Project administration, Supervision, Validation, Visualization, Writing &#x; review &#x; editing

^{}

^{ 2 ,}

^{*}

### Sourav Roy

^{1} Department of Mathematics, Jadavpur University, Kolkata, West Bengal, India

Find articles by Sourav Roy

### Sayantan Nag Chowdhury

^{2} Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata, West Bengal, India

Find articles by Sayantan Nag Chowdhury

### Prakash Chandra Mali

^{1} Department of Mathematics, Jadavpur University, Kolkata, West Bengal, India

Find articles by Prakash Chandra Mali

### Matja Perc

^{3} Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia

^{4} Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan

^{5} Alma Mater Europaea, Maribor, Slovenia

^{6} Complexity Science Hub Vienna, Vienna, Austria

Find articles by Matja Perc

### Dibakar Ghosh

^{2} Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata, West Bengal, India

Find articles by Dibakar Ghosh

The Anh Han, Editor

Author informationArticle notesCopyright and License informationDisclaimer

^{1} Department of Mathematics, Jadavpur University, Kolkata, West Bengal, India

^{2} Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata, West Bengal, India

^{3} Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia

^{4} Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan

^{5} Alma Mater Europaea, Maribor, Slovenia

^{6} Complexity Science Hub Vienna, Vienna, Austria

Teesside University, UNITED KINGDOM

^{}Corresponding author.

**Competing Interests: **We have declared that no competing interests exist.

* E-mail: [email protected]

Received Jun 20; Accepted Jul

Copyright © Roy et al

This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

### Abstract

Most environments favor defection over cooperation due to natural selection. Nonetheless, the emergence of cooperation is omnipresent in many biological, social, and economic systems, quite contrary to the well-celebrated Darwinian theory of evolution. Much research has been devoted to better understanding how and why cooperation persists among self-interested individuals despite their competition for limited resources. Here we go beyond a single social dilemma since individuals usually encounter various social challenges. In particular, we propose and study a mathematical model incorporating both the prisoner&#x;s dilemma and the snowdrift game. We further extend this model by considering ecological signatures like mutation and selfless one-sided contribution of altruist free space. The nonlinear evolutionary dynamics that results from these upgrades offer a broader range of equilibrium outcomes, and it also often favors cooperation over defection. With the help of analytical and numerical calculations, our theoretical model sheds light on the mechanisms that maintain biodiversity, and it helps to explain the evolution of social order in human societies.

### 1 Introduction

Numerous interdisciplinary researchers have long sought a way to understand how the tremendous biodiversity among species persists in nature despite the significant differences between them in terms of competitive capability. Darwinian evolution [1] always challenges the emergence of spontaneous social cooperation, as cooperators act entail an inherent cost for displaying altruism [2]. Self-interested defectors generally exploit these cooperators and hinder the maintenance of cooperation. This trivially raises the question of why cooperators act selflessly if only the fittest succeeds and how numerous forms of cooperative behavior sustain in nature ranging from a low-level microbial biological system to a high-level complex social system. We resort to the evolutionary game dynamical interaction given by a 2 × 2 pay-off matrix to solve this riddle. These classical games are capable of analyzing the economic and strategic decisions of rational individuals [3]. An unprecedented spectrum of researchers have already focused overwhelmingly on the evolutionary game theory [4&#x;6] to seek out the mechanisms that sustain and promote cooperation [7&#x;14].

We construct a mathematical model on evolutionary multigames [15&#x;27] by adopting two simple two-player games, viz. the prisoner&#x;s dilemma [28] (PD) and snowdrift game [29, 30] (SD), as metaphors for cooperation between unrelated rational individuals. Our choice of combining two distinct games is motivated by the fact that people may face different social complexities rather than a single one. People take care of problems based on their own perceptions. This real-life behavioral heterogeneity among individuals inspires us to examine the effect of multigames for understanding social diversity. Our theoretical approach is another effective way to explore the impact of the interaction between different cultural backgrounds on large-scale social ordering [31].

While defection is the ultimate profitable strategy in the classical PD game, irrespective of the co-player&#x;s decision, a player&#x;s best strategy in the SD game depends on the co-player&#x;s decision. Thus, we select two different games with two contrasting outcomes. The stable coexistence of both cooperators and defectors is the expected consequence in the SD game, resulting in the persistence of cooperative behavior. In comparison, cheaters are always encouraged to exploit the cooperative individuals in the set-up of the PD game. In our constructed model, unrelated individuals have their probabilistic selection to choose which game they want to play. We also consider the mutation [32&#x;42] as an evolutionary mechanism that enables people to switch their respective strategies. We assume that it is hard to behave rationally under every circumstance. Hence, in our mathematical model, we incorporate a simplistic assumption that cooperators&#x; and defectors&#x; subpopulations can interchange their strategies with a fixed bidirectional mutation rate. To the best of our knowledge, the majority of the studies on multigames have been examined in the absence of mutation.

Apart from that, to reveal the interplay between ecological and evolutionary dynamics, we inspect the influence of free space in our mathematical model. In the context of complex systems, free space [43&#x;49] proves to be a promising factor for offering significant consequences on diverse collective dynamics. We assume that each individual&#x;s birth rate depends not only on the respective average payoffs but are also proportional to the available free space. This free space will give them reproductive opportunities. The population of unrelated players gains a payoff due to the interaction of the multigames. We treat these payoffs as reproductive successes. Individuals with a higher payoff leave more offspring if sufficient empty spaces are available and are able to outcompete less successful ones. Free space will also contribute to the individual&#x;s average payoffs. We consider free space as an ecological variable donating selflessly to all players&#x; fitness. However, free space never anticipates any benefit from others. This generous nature is surprising in a society of self-interested individuals. Still, free space provides reproductive opportunities to each player and loses their own identity to improve the fitness of other individuals. We consider such a selfless act of free space because helping others is a common practice among human beings and other animals. This tendency of unselfish concern for other people is a unique recipe for promoting cooperation and favors the survival of various species. The inclusion of such ecological dynamics into evolutionary multigames gives rise to a fascinating way to reveal the influence of eco-evolutionary dynamics [50&#x;60] on decision-making and the evolution of cooperation.

The concurrence of ecological changes and rational players&#x; evolution, depicted in the present article, prevails through the same time scale, rather than the typical assumption that ecological processes are much faster than evolutionary processes [61]. Our theoretical model uncovers how the combination of ecological dynamics and game dynamics is beneficial for maintaining cooperative behaviors under the influence of mutation and leads to the stable coexistence of interacting competitors. Earlier, Nag Chowdhury et al. [8, 16] examined the influence of eco-evolutionary dynamics on cooperative behavior&#x;s emergence and evolution under a cooperation-supporting mechanism, viz. punishment. Nevertheless, our model exhibits the coexistence of cooperators and defectors in the absence of any such supporting mechanism. Moreover, individuals&#x; birth rates in those Refs. [8, 16] depend only on their respective average payoffs. Apart from the addition of free space&#x;s altruistic behavior in the payoff matrix, we furthermore include the importance of available free space on the birth rates by assuming their birth rates to be proportional to the available space. The rest of the article is structured as follows: In Sec. (2), we discuss how we construct our mathematical model in detail. The following section (3) deals with the existence, uniqueness, and boundedness of the model, along with extensive numerical simulations and discussions. The section (4) presents the summary of our findings. Following that, we provide a brief discussion in the last section (5) on the challenges in this work&#x;s context that need to be addressed and are worth studying in the future.

### 2 The model

To start with, we consider a simplistic assumption that each individual has two distinct choices, viz. (i) cooperation (**C**) and (ii) defection (**D**). Even they can play any of the two possible games (a) PD game and (b) SD game. They can adopt the PD game with probability *p*, and alternatively, they interact with other individuals by playing the SD game with the complimentary probability (1 &#x; *p*). Both of these two-person games can be given by the following two payoff matrices *A* and *B* respectively, where

in which the entries portray the payoff accumulated by the players in the left.

Here, *R*_{PD} and *R*_{SD} contemplate the reward for mutual cooperation among two players in the respective PD and SD games. Similarly, both unrelated players receive the punishment *P*_{PD} and *P*_{SD} for mutual defection in the games PD and SD, respectively. An exploited cooperator gains the sucker&#x;s payoff *S*_{PD} and *S*_{SD}, respectively in the PD and SD games when confronted by a defector. The mixed choice yields the defector temptations *T*_{PD} and *T*_{SD} to exploit a cooperator in the PD and SD games, respectively. The payoff ranking of these four-game parameters determines the two-person games. The conventional relative ordering for the PD game is *T*_{PD} *R*_{PD} *P*_{PD} *S*_{PD} [8, 16, 31] and 2*R*_{PD} *S*_{PD} + *T*_{PD} [62]. Without loss of any generality, we choose *T*_{PD} = *²* 1, *R*_{PD} = 1, *P*_{PD} = *·* &#x; (0, 1), and *S*_{PD} = 0. Similarly, we choose *T*_{SD} = *²* 1, *R*_{SD} = 1, *S*_{SD} = 0, and *P*_{SD} = &#x;*·* &#x; (&#x;1, 0), maintaining the relative ordering *T*_{SD} *R*_{SD} *S*_{SD} *P*_{SD} for the standard SD game [8, 16, 31]. Thus, the payoff matrices *A* and *B* become

Note that, in both of these games, mutual cooperation leads to the payoff *R*_{PD} and *R*_{SD}, which is relatively higher than *P*_{PD} and *P*_{SD}, which one defector receives when playing with a defector. Thus, cooperation always promises higher income than defection if both the rational players choose the same strategy. The difference between these two games&#x; relative ordering leads to a contrasting scenario. In the SD game, the interaction between the cooperator and defector always promises a better income in terms of payoff than the interaction between two defectors. A reverse reflection is observed in the case of the PD game thanks to the choice of such relative ranking of game parameters in the PD game. The interaction between two defectors in the PD game allows them to earn more than a cooperator encountering a defector. The switching between *P* and *S* in the relative ordering of both games thus produces noticeable unexpected consequences on the evolution of cooperation.

Since a player can decide which game they want to play, thus the final payoff matrix for the multigames looks like

Here, *p* &#x; [0, 1] is the probability of playing the PD game. Moreover, we consider free space **F** as an ecological variable, which contributes altruistically by helping others. Nevertheless, free space does not get any benefits by giving them reproductive opportunities. We incorporate this charitable role of free space by extending the 2 × 2 payoff matrix *E* to the 3 × 3 payoff matrix *G* as follows

The matrix *G* clearly reveals that the free space never earns any payoff for their selfless charitable act; however, it contributes a positive payoff *Ã*_{1} and *Ã*_{2} to the cooperators and the defectors, respectively. Since most of the game parameters (not all) lies within the closed interval [0, 1], we assume, for the sake of feasible comparison, *Ã*_{1} and *Ã*_{2} both lies within the interval [0, 1]. When *Ã*_{1} and *Ã*_{2} are equal to zero, free space will not contribute anything to anyone. However, whenever *Ã*_{1} and *Ã*_{2} attain positive values, individuals gain an additional payoff from free space.

Inspired by the Malthusianism, we consider the following set of differential equations governing the changes in frequencies of cooperators and defectors as a function of time *t*

(1)

where

Here, *x* and *y* are the normalized densities of cooperators and defectors, respectively. Let *z* be the available free space. Thus, we have

Relation (2) assures that by studying *x* and *y* alone, one can easily capture the dynamics of the two strategies. We assume the birth rates of each individual depends crucially on the available free space as well as on their respective average fitness. Thus, we consider

(3)

where *f*_{C} and *f*_{D} are average fitness of the cooperators and the defectors, respectively. The average payoff of cooperators and defectors can be determined using the payoff matrix G of the multigames, and the relation (2) as follows,

(4)

Note that the average fitness of free space is

This is expected as free space does not gain anything for its benevolent nature. For simplicity, we further assume that all individuals die at a uniform and constant mortality rate *¾* &#x; (0, 1]. Hence using the relations (3) and (4), our constructed model (1) transforms into

(6)

Now, we introduce a constant probability *¼* as a rate with which each individual mutates from one strategy to the others, in a continuous manner,

(7)

Relation (7) reflects that the mutation probability *¼* from the cooperators to the defectors is identical to the mutation rate from the defectors to the cooperators. So using the system (6), the well-mixed population under the influence of bidirectional mutation gives rise to the differential equations

(8)

The system (8) contains seven different parameters. We summarize the necessary information about these parameters in (Table 1).

### Table 1

**Parameters with their physical significance and domains: The first column represents the symbols used in the present manuscript to describe several parameters.**

At the same time, the second column depicts what these parameters exactly mean. The set of possible input values from which the parameters can assume their values is given in the third column of the table. For further details, please see the main text.

Parameters | Physical Significance | Domain |
---|---|---|

¾ | Death Rate | (0, 1] |

² | The gain of a defector while interacting with a cooperator | (1, 2) |

· | The payoff for mutual defection | (0, 1) |

p | Probability of playing the PD game | [0, 1] |

¼ | Mutation probability | [0, 1] |

Ã_{1} | The altruistic incentive of free space towards the cooperators | [0, 1] |

Ã_{2} | The altruistic incentive of free space towards the defectors | [0, 1] |

Open in a separate window

### 3 Results and discussions

### Existence, uniqueness and positive invariance

Before investigating the model (8) using numerical simulations, we first prove the positive invariance of the proposed system (8). Clearly, the functions on the right-hand side of the differential Eq (8) are continuously differentiable and at the same time, locally Lipschitz in the first quadrant of . This ensures the existence and uniqueness of solutions for the model (8) with suitable initial values. Note that the initial conditions (*x*_{0}, *y*_{0}) must lie within the domain [0, 1] × [0, 1] maintaining the inequality 0 &#x; *x*_{0} + *y*_{0} &#x; 1, as they represent the frequencies of the two strategies. To determine the positivity of the proposed model (8), we write the system as follows,

(9)

where and are two integrable functions in the Riemannian sense. Solving Eq (9), we get

(10)

where *c*_{1} and *c*_{2} are the integrating constants depending on the initial densities *x*_{0} and *y*_{0}. This proves that both *x* and *y* are non-negative. Now, to calculate the upper bound of *x* + *y*, we proceed as follows. The dynamical Eq (8) yield

(11)

Since as per our previous analysis, *x* + *y* &#x; 0 and the mortality rate *¾* is always positive, we have *¾*(*x* + *y*)&#x;0. Thus, (11) reduces to

(12)

Integrating both sides, we get

(13)

where *c*_{3} is the initial density dependent constant. Hence, we find that

(14)

i.e., the overall species density *x* + *y* eventually remains bounded within the region [0, 1]. This boundedness within the closed interval [0, 1] allows us to relate the possible emerging dynamics of the system (8) to physically implementable scenarios with biological relevance. When the sum of the population density (*x* + *y*) is precisely one, the available free space is zero, as *z* = 1 &#x; (*x* + *y*). i.e., there is no reproductive opportunity accessible to any individual in that situation with *z* = 0. When (*x* + *y*) = 0, then individually *x* = 0 and *y* = 0. Hence, all individuals die, and *z* is equal to one. Thus despite the presence of ample free space, all the individuals are extinct under that circumstances.

### Coexistence of different stationary points depending on the initial conditions

Next, we point out the multistable dynamics of our model (8), resulting in the system&#x;s vulnerability to small perturbations. Initially, we set the parameter values at *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = , *Ã*_{2} = and *¼* = 0 in Fig 1(a). We vary the initial conditions (*x*_{0}, *y*_{0}) within the interval [0, 1] × [0, 1] maintaining the inequality (*x*_{0} + *y*_{0})&#x;[0, 1]. We find that the dynamics switching between two stationary points, viz. *E*_{0} = (0, 0) and *E*_{1} = (*x**, 0). We analytically calculate the stationary points of the system (8) in the absence of the mutation, i.e., with *¼* = 0. We trace out four different stationary points,

Open in a separate window

Fig 1

**Alternations between multiple co-existing steady states depending on initial conditions.**

The coexistence of various stable states is portrayed here by varying the initial conditions maintaining 0 &#x; (*x*_{0} + *y*_{0})&#x;1. Red points signify the extinction equilibrium *E*_{0}. Violet points represent the defector-free steady state *E*_{1}. The cooperator-free stationary points *E*_{2} are shown by sea blue points, and the coexistence equilibrium *E*_{3} is plotted using yellow points. The mutation-free model in subfigures (a-c) allows four stable steady states to coexist. However, subfigure (d) supports only bistability for the chosen parameters&#x; values. The parameter values for each of these subfigures are (a) *¾* = , *²* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. (b) *¾* = , *²* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. (c) *¾* = , *²* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. (d) *¾* = , *p* = , *²* = , *¼* = , *Ã*_{1} = , and *Ã*_{2} = Other parameter is *·* =

*E*_{0}= (0, 0) reveals all individuals die. This point is locally stable when*E*_{1}= (*x**, 0) exhibits a society free from any defectors. Here, and*x** &#x; (0, 1]. The stability criteria is given by*E*_{2}= (0,*y**) represents that we have only left with defectors. Here, and*y** &#x; (0, 1]. This stationary state is stable if*E*_{3}= (*x***,*y***) indicates the coexistence equilibrium offering the survival of cooperation and defection simultaneously. Here, , and*x*** satisfies the equation, , where, .*x*** and*y*** both should lie within the interval (0, 1). The local stability yields the conditions for the stability of (*x**,*y**) are

Clearly, the chosen parameters satisfy the local stability criteria of both stationary points *E*_{0} and *E*_{1}. Note that *E*_{1} leads to two different values for the selected values of the parameters, out of which (, 0) always remains locally stable, and (, 0) is unstable. Thus, we see the appearance of two different stationary points *E*_{0} (shown by red points) and *E*_{1} (displayed by violet points) in the basin of attraction portrayed in Fig 1(a). Such toggling between alternate stable stationary points is one of the generic features in some biological systems involving the fundamental processes of life [63&#x;65] and in a few nonlinear dynamical systems [66&#x;69]. The initial condition (*x*_{0}, *y*_{0}) = (0, 0) always helps the system (8) to converge to the stationary point *E*_{0}. The system is never able to give rise to the survival of any cooperators and defectors without the presence of any individual at the beginning. Since the chosen initial point (0, 0) itself is the stationary state, the system always stabilizes in the (0, 0) stationary point irrespective of the parameter values. Fig 1(b) is drawn with *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = , *Ã*_{2} = , and *¼* = 0. Similarly, we find the system (8) in the absence of mutation (i.e., *¼* = 0) converges to *E*_{0} for the single initial condition (*x*_{0}, *y*_{0}) = (0, 0). It is anticipated that the choice of *x*_{0} = 0 always leads to the cooperator-free steady state. The initial absence of cooperators in the mutation-free model will not entertain any exposure for the cooperators in the long run. The line of initial conditions *x*_{0} = 0 and *y*_{0} &#x; 0 produces the stationary state *E*_{2} (sea blue points in Fig 1(b)). Apart from these two steady states, the mutation-free system also switches between *E*_{1} (violet points) and *E*_{3} (yellow points) depending on the suitable choice of initial conditions (See Fig 1(b)). Thus, for the same choices of parameters&#x; values, the system flips between four alternate steady states depending on the initial densities of cooperators and defectors. Similarly, we find the system (8) with the parameters&#x; values *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = , *Ã*_{2} = and *¼* = 0 converges to all these four stationary points. All population extincts for the initial conditions ranging from (0, 0) to (, 0) (See red points in Fig 1(c)). The system (8) can be solved analytically with *Ã*_{1} = *Ã*_{2} = *¼* = 0 in the absence of defectors (i.e., *y* = 0) as follows,

where *c*_{4} is the integrating constant.

Similarly, the system (8) can be solved analytically with *Ã*_{1} = *Ã*_{2} = *¼* = 0 in the absence of cooperators (i.e., *x* = 0) as follows

where *c*_{5} is the integrating constant.

Fig 1 is plotted by solving the differential Eq (8) by varying the initial conditions within [0, 1] × [0, 1] with fixed step-length *´x*_{0} = *´y*_{0} = and maintaining 0 &#x; *x*_{0} + *y*_{0} &#x; 1. To solve our proposed system (8) numerically, we use the 4th order Runge-Kutta (RK4) method with 20 × 10^{5} iterations with fixed integration step length *´t* = The final point (*x*, *y*) is stored to decide the asymptotic dynamics of the system. Other initial conditions with *y*_{0} = 0 stabilize the dynamics in the defector-free state *E*_{1} (violet points). We also track a fair portion of the basin in Fig 1(c), where the system (8) with *¼* = 0 converges to either *E*_{2} (sea blue points) or *E*_{3} (yellow points) depending on the choice of initial conditions.

Although we choose the mutation-free model for the Fig 1(a)&#x;1(c), we consider the contribution of *¼* in Fig 1(d). We generate Fig 1(d) with *¾* = , *p* = , *²* = , *·* = , *¼* = , *Ã*_{1} = , and *Ã*_{2} = This non-zero *¼* leads to the disappearance of two stationary points, (i) the defector-free steady state *E*_{1}, and (ii) the cooperator-free steady state *E*_{2}. The symmetric mutation from one species to another species always gives them two feasible opportunities. Either both strategies survive or all the individuals perish for *¼* &#x; 0. We also mathematically derive two possible stationary points as follows,

- The extinction equilibrium
*E*_{0}= (0, 0), which is stable under the conditions, - The interior equilibrium
*E*_{3}= (*x**,*y**) becomes stable, if

The compelling evidence of bistability under the same choice of parameters&#x; values is recognized in Fig 1(d).

### Emergent dynamics in absence of mutation

We examine the impact of the parameters *p*, *²*, *·*, and *¾* on the system (8) with *¼* = 0 in Fig 2. For the comparability, we choose *Ã*_{1} = *Ã*_{2} = 0 in subfigures (a), (c) and (e). Also, we set *Ã*_{1} &#x; 0 and *Ã*_{2} &#x; 0 for the subfigures (b), (d) and (f). In the absence of mutation (*¼* = 0) and free space induced benefits (*Ã*_{1} = *Ã*_{2} = 0), the defector-free is unable to stabilize. Although *E*_{1} exists for , but two eigenvalues of the Jacobian of the linearized system are

Open in a separate window

Fig 2

**The interplay of different parameters in the mutation-free model (8).**

*p* is the probability of playing the PD game. We investigate the role of *p* in our proposed model (8) in the absence of the mutation. The figures (a), (c), and (e) of the left-handed column are drawn with *Ã*_{1} = *Ã*_{2} = 0. The remaining figures in the right-hand column are plotted with non-zero values of *Ã*_{1} and *Ã*_{2}. The advantages provided by the free space help to promote cooperation, as portrayed through the broad area of yellow and violet regions in figures (b), (d), and (f). Red, sea blue, violet, and yellow depict the stable extinct state, cooperator-free state, defector-free state, and co-existence state, respectively. We also plot our calculated curves from the local stability analysis of the stationary points. The white line represents that of *E*_{0}, whereas the blue line portrays that of the coexistence equilibrium *E*_{3}. The brown dotted line is the stability curve for the cooperator-free stationary state *E*_{2}, while the blue dashed line represents that of the defector-free steady state *E*_{1}. The mismatches with the stability curve in a few places are due to the multistable behavior of our constructed model (8). Other parameter values are kept fixed at (a-d) *¾* = , (c-f) *²* = , (a,b,e,f) *·* = , (b,d,f) *Ã*_{1} = , (b,d,f) *Ã*_{2} = and *¼* = 0.

Since *²* 1 and *¾* 0, »_{2} is always positive. Consequently, *E*_{1} is always unstable. This result is physically meaningful, as both the chosen games (PD and SD) do not encourage a defection-free society without any supporting mechanism. However, we find three different steady states in Fig 2(a), 2(c) and 2(e). These stationary points and their corresponding stability analysis are given below,

The extinction equilibrium

*E*_{0}= (0, 0) (red points in subfigures (a), (c) and (e) of Fig 2) always exists and is always locally stable.*E*_{2}= (0,*y***) is the cooperator-free stationary point (sea blue points in subfigures (a), (c) and (e) of Fig 2), where . This state exists under the condition 0 &#x; 4*¾·*(2*p*&#x; 1)&#x;*·*^{2}(2*p*&#x; 1)^{2}, and becomes stable if*¾**·*(2*p*&#x; 1)(2*y*** &#x; 3*y***^{2}).*E*_{3}= (*x**,*y**) allows the coexistence of both cooperators and defectors (yellow points in subfigures (a), (c) and (e) of Fig 2), where , and , with . The stationary state is stable under the following conditions,

With increasing *p*, players are prone to play the PD game. Hence, we can not anticipate the survival of cooperators for large *p*. Thus, for large *p*, either the extinction equilibrium *E*_{0} or the cooperator-free stationary point *E*_{2} stabilizes in Fig 2(a), 2(c) and 2(e). However, the free-space induced benefits in Fig 2(b), 2(d) and 2(f) facilitate the emergence of cooperation and stabilize the interior equilibrium *E*_{3}. Even, Fig 2(f) provides a range in the *p* &#x; *¾* parameter space, where the defector-free stationary point *E*_{1} stabilizes. The comparative study between the left and right column of Fig 2 ensures the encouraging role of free space in the promotion of cooperation under suitable circumstances.

With increasing *²*, the defectors are getting extra aid, and thus, we get the stabilization of the extinction equilibrium *E*_{0} (red zone) in Fig 2(a). However, the presence of free space assistance converts that area into a coexistence zone (yellow zone) in Fig 2(b). The increment of *p* in both subfigures (a-b) permits only the existence of defectors, as people are playing mostly the PD game in such circumstances. Defectors are getting a favorable atmosphere in the PD game for our chosen PD game parameters&#x; values. A similar sort of stabilization of cooperator-free steady state (sea blue area) is observed in subfigures (c-d) for larger values of *p*. However, lower values of *p* provide the opportunity for playing the SD game. Hence, coexistence equilibrium (yellow region) stabilizes in Fig 2(c) and 2(d) for smaller *p*. Nevertheless, the inclusion of free space induced benefits helps to broaden the region of coexistence in the *p* &#x; *·* parameter plane, as shown in Fig 2(d). The comparison between the Fig 2(e) and 2(f) suggests that appropriate non-zero values of *Ã*_{1} and *Ã*_{2} entertain the stabilization of a defector-free steady state (violet region). Hence, we can trace a fair portion of violet points in the *p* &#x; *¾* parameter space of the subfigure (f). However, too large a mortality rate reduces the opportunity of survivability of any individual, resulting in the stabilization of the extinction equilibrium *E*_{0}. We focus on the effect of mortality rate *¾* more elaborately in Fig 3.

Open in a separate window

Fig 3

**Influence of death rate**

*¾*in the mutation-free model (8).The influence of game parameters *²* and *·* and the mortality rate *¾* on the emergent dynamics of the model (8) in the absence of mutation is examined here. Subfigures (a) and (c) are drawn with *Ã*_{1} = *Ã*_{2} = 0. While subfigures (b) and (d) are generated with *Ã*_{1} = and *Ã*_{2} = Since, *Ã*_{1} *Ã*_{2}, we notice a fair portion of defector-free region (violet points) in subfigures (b) and (d). In the absence of free space induced benefits, we are unable to trace such defector-free regions in subfigures (a) and (c). The red region reflects the disappearance of all individuals for relatively high values of death rate *¾*. Note that we keep the value of *p* fixed at Thus, the system gets more opportunities to play the SD game, which generally encourages stabilizing the coexistence equilibrium. Therefore, we notice lower values of *¾* will lead to the convergence towards the coexistence equilibrium *E*_{3} (yellow points). *²* is kept fixed at for the subfigures (a-b) and *·* is set at for subfigures (c-d). We iterate the system (8) with *¼* = 0 for 20 × 10^{5} iterations with fixed integration step size *´t* = and fixed initial condition (*x*_{0}, *y*_{0}) = (, ).

As expected, higher values of *¾* constantly enlarge the chances of extinction. Thus, we observe a fair portion of the red region in Fig 3. Nevertheless, the amount of this red area is considerably lesser in the right column of Fig 3 compared to the left column. We introduce the non-zero values of *Ã*_{1} and *Ã*_{2} in the right column of this figure. These free space induced benefits encourage maintaining a defector-free society, as shown in Fig 3(b) and 3(d). We choose *Ã*_{1} = = *Ã*_{2} for Fig 3(b) and 3(d). i.e., the free space will provide an additional advantage for the cooperators, and thus depending on the other parameters, the defectors are vanished in the long run, as shown in Fig 3(b) and 3(d). In all of these subfigures of Fig 3, we trace a portion of yellow points depicting the survival of both cooperators as well as defectors simultaneously. We choose *p* = in Fig 3. Thus, people get more chances to play the SD game, which facilitates the concurrence of both cooperation and defection. Thus, smaller values of *¾* provide an opportunity to coexist for all strategies, which is observed in Fig 3 with initial condition (*x*_{0}, *y*_{0}) = (, ). We plot the stability curve of all stationary points in Figs 2&#x;5 as follows,

Open in a separate window

Fig 5

**Investigating the impact of**

*Ã*_{2}on the evolution of cooperation.For the smaller values of *Ã*_{2}, the mutation-free system gives rise to the defector-free (violet region) society. While larger values of *Ã*_{2} encourage the defectors and stabilize the cooperator-free stationary point *E*_{2} (sea blue region). Suitable choices of parameters may favor the co-existence of cooperators as well as defectors. Even for a higher mortality rate *¾*&#x;s values, all individuals die in the red portion of the parameter space *Ã*_{2} &#x; *¾* in subfigure (d). All the parameters are set at *¾* = , *·* = , *²* = , *p* = , *Ã*_{1} = , and *¼* = 0, unless it is varied. The simulations are computed for 20 × 10^{5} iterations with fixed initial condition (*x*_{0}, *y*_{0}) = (, ). The color codes are as follows: (i) red represents the extinct state *E*_{0}, (ii) yellow represents the co-existence state *E*_{3}, (iii) violet represents the defector-free state *E*_{1}, and (iv) sea blue represents cooperator free state *E*_{2}, respectively. We plot the stability curves corresponding to each stationary point in each subfigure.

(i) The brown dotted line for the cooperator-free steady state

*E*_{2},(ii) the blue dashed line for the defector-free steady state

*E*_{1},(iii) the solid blue line for the interior equilibrium

*E*_{3}, and(iv) the solid white line for the extinction equilibrium

*E*_{0}.

In Fig 4, we fix the value of *Ã*_{2} at and examine the role of *Ã*_{1}. Since *Ã*_{1} represents the free space induced benefits towards the cooperators, thus enhancement of its (*Ã*_{1}) value will help in promoting cooperation in the society. Fig 4 reflects the same scenario. A defector-free society (violet region) is noticed in all these subfigures. The increment of temptation parameter *²*&#x;s value challenges the prevalence of cooperation. Thus, for a small *Ã*_{1}, we find a defector dominated society (the sea blue region in Fig 4(a)). *Ã*_{1} proves to be a cooperator facilitating parameter as we detect a wide range of yellow regions in Fig 4(a), where both cooperators and defectors can coexist. Similarly, irrespective of the choice of *·* in Fig 4(b), higher values of *Ã*_{1} provide all cooperators an extra benefit for survival, and thus, the stationary point *E*_{1} (violet points) stabilizes. For the intermediate choice of *Ã*_{1}, the stationary point *E*_{3} (yellow) yields the stable coexistence of both strategies. However, cooperators strive to keep in existence for the lower values of *Ã*_{1}, and we find the sea blue region of cooperator-free steady state. The parameter *p* indicates the probability of playing the PD game. Thus, *p* &#x; 1&#x; always gives defectors a more favorable environment to survive. That&#x;s why we spot a sea blue portion in Fig 4(c) for higher values of *p* and *Ã*_{1}. When *p* is small, people are prone to play the SD game; and thus, we get the coexistence of both strategies (the yellow region in Fig 4(c)) for smaller *p* and *Ã*_{1}. Nevertheless, a larger value of *Ã*_{1} with a moderate value of *p* always provides a reasonable scope for the cooperators to survive, and we discover a healthy portion of the stationary point *E*_{1} (the violet region) in Fig 4(c). A higher mortality rate *¾* always results in the extinction of both cooperators and defectors, and we locate a huge red region in the *Ã*_{1} &#x; *¾* parameter plane of Fig 4(d). We find a very tiny sea blue region in Fig 4(d), where defectors are only able to survive. But, as expected, a higher value of *Ã*_{1} always promotes the cooperation strategy, and we obtain a violet zone of defector-free stationary point and a yellow region of interior equilibrium *E*_{3} in Fig 4(d).

Open in a separate window

Fig 4

**Importance of**

*Ã*_{1}in the enhancement of cooperation.Four different parameter spaces (a) *Ã*_{1} &#x; *²*, (b) *Ã*_{1} &#x; *·*, (c) *Ã*_{1} &#x; *p*, and (d) *Ã*_{1} &#x; *¾* are contemplated here with fixed initial condition *x*_{0} = and *y*_{0} = *Ã*_{1} is varied within [0, 1], and the other parameters&#x; values are for the subfigures (a) *¾* = , *·* = , *p* = , *Ã*_{2} = , and *¼* = 0. *²* is varied within the open interval (1, 2), (b) *¾* = , *²* = , *p* = , *Ã*_{2} = and *¼* = 0. *·* is varied within (0, 1), (c) *¾* = , *²* = , *·* = , *Ã*_{2} = , and *¼* = 0. *p* is varied within the closed interval [0, 1], and (d) *²* = , *·* = , *p* = , *Ã*_{2} = and *¼* = 0. *¾* is varied within the interval (0,1]. The color code represents the following: (i) red represents the extinct state, (ii) yellow portrays the co-existence state, (iii) violet displays the defector-free state, and (iv) sea blue depicts cooperator-free state, respectively. We have run the numerical simulations for 20 × 10^{5} iterations for each point and store the final value for determining the final asymptotic state. Increment of *Ã*_{1} contributes more to the cooperators&#x; payoff, and hence, we observe a defector-free region for higher values of *Ã*_{1} depending on the other parameters. We draw the stability curves for *E*_{0} (the solid white line), *E*_{1} (the blue dashed line), *E*_{2} (the brown dotted line), and *E*_{3} (the solid blue line).

Fig 5(a) shows that except for a smaller portion of the defector-free region (violet zone), the whole *Ã*_{2} &#x; *²* parameter space produces the coexistence of both cooperators and defectors depending on other parameters&#x; values. Despite the increment of *Ã*_{2} and *²*, the cooperators are able to survive along with the defectors due to our choice of other parameters&#x; values. The *Ã*_{2} &#x; *·* parameter space portrays that smaller values of *Ã*_{2} can not provide any benefit to the defectors, and stabilize the defector-free stationary point *E*_{1} (violet region) irrespective of *·*&#x;s value. However, if *Ã*_{2} increases, it will yield a window of opportunity for the defectors to thrive. We observe a wide coexistence region (yellow region) and a small area of cooperator-free stationary point (sea blue region)in Fig 5(b). *Ã*_{2} indicates the free-space induced benefits towards the defectors. So, it is expected that larger values of *Ã*_{2} always enhance the chances of defectors&#x; survivability. Thus, we notice a cooperator-free region (sea blue zone) in both Fig 5(c) and 5(d). Nevertheless, larger values of *p* enhance the chances of playing the PD game, where defectors get a favorable environment to survive. Thus, the cooperator-free sea blue region is found in Fig 5(c). We are also able to detect a small violet region of a defector-free environment in the *Ã*_{2} &#x; *p* parameter space. However, we trace a healthy portion of coexistence (yellow) too in Fig 5(c) due to our chosen parameters&#x; values.

### The influence of bidirectional mutation

Now we investigate the influence of bidirectional mutation on the long-term behavior of the nonlinear differential Eq (8). Since every species can mutate into the other at a specific uniform rate *¼* &#x; (0, 1], thus cooperators and defectors can not remain alive alone. Either all populations will die; otherwise, the dynamics will lead to the coexistence of all two species. Fig 4 already portrays the influence of free space-induced benefits on the cooperators in the absence of mutation. We scrutinize the impact of *Ã*_{1} under the influence of 50% mutation (i.e., *¼* = ) in Fig 6. As *Ã*_{1} increases, the cooperators are getting a better environment for survival. We find a portion of stable coexistence equilibrium in each two-dimensional parameter space for large *Ã*_{1} in Fig 6. Thanks to the mutation, the cooperators can not live alone. The white line in Fig 6 is the stability curve corresponding to the extinction equilibrium. The local stability analysis fits almost exactly with the numerical simulations in Fig 6 with fixed initial condition (*x*_{0}, *y*_{0}) = (, ). There are a few places where the stability analysis fails to predict the stabilization of the extinction stationary point (0, 0). This is mainly due to the multistable behavior of the proposed model 8.

Open in a separate window

Fig 6

**The effect of altruist free space on the nonlinear dynamics of multigames with mutation.**

The influence of free space-induced benefits by varying *Ã*_{1} within the closed interval [0, 1] is established. The whole population goes extinct in the red region, and the yellow area reflects the system&#x;s stable interior point, corresponding to the coexistence of all two strategies. All the parameters are kept fixed at *¾* = , *·* = , *p* = , *Ã*_{2} = , *²* = , *¼* = , unless it is varied. The initial condition is kept fixed at (, ). A notable difference is observed from (Fig 4), which is drawn in the absence of mutation. Larger values of *Ã*_{1} always facilitate the maintenance of cooperation, and the bidirectional mutation reinforces the system&#x;s inherent tendency to flow from cooperators to defectors. The solid white line in each figure is the analytically computed stability curve corresponding to the extinction equilibrium.

As discussed, the increasing values of the parameters *²*, *·*, and *p* always provide the defectors a favorable environment to dominate the cooperators. However, beyond a critical value of *Ã*_{1}, both the cooperators and defectors can coexist, as shown in Fig 6(a)&#x;6(c). The larger values of *¾* always hinder the evolution of cooperators as well as of defectors. However, an intermediate choice of *Ã*_{1} &#x; *¾* favors the successful evolution of both strategies, as portrayed through Fig 6(d). The same feature is also noticeable in Fig 7. The complex evolutionary dynamics switch between two stationary points depending on the choices of parameters&#x; values in Fig 7. For smaller values of *Ã*_{2}, the defectors are not getting enough advantages to survive, and thus the extinction equilibrium (0, 0) stabilizes in the red region of all subfigures of Fig 7. The choice of other parameters&#x; values is also crucial for obtaining these stationary states. The parameter *Ã*_{2} benefits the defectors; hence the defectors can survive beyond a certain threshold of *Ã*_{2}. The employed 50% mutation rate helps to flow a certain fraction of defectors into cooperators, and we have a moderate portion of coexistence state (yellow region) in Fig 7. The observed results may vary for different choices of initial conditions, as the system is multistable. We plot the boundary separating solid white lines in all subfigures by analyzing the local stability analysis of the extinction equilibrium. Clearly, this stability curve agrees well with our numerical simulations, and the places, where they don&#x;t agree with the numerical simulations, is due to the multistable behavior of our proposed model 8. All the simulations are done by iterating for 20 × 10^{5} times with fixed integrating step length *´t* = The last point is gathered to finalize the asymptotic state. All the codes to generate these figures are publicly available at Ref. [70].

Open in a separate window

Fig 7

**The impact of free space-induced benefits on the defectors in the presence of the mutation.**

Either all populations die (red region), or the symmetric mutation preserves the coexistence state (yellow area) in all these subfigures. Beyond a critical value of *Ã*_{2}, both strategies mutate one into another, providing the opportunities for the coexistence of both species. The solid white line is the analytically derived stability curve to stabilize the extinction equilibrium. All the parameters are kept fixed at *¾* = , *·* = , *p* = , *Ã*_{1} = , *¼* = , *²* = , unless they are varied. The initial condition is kept fixed at (, ). The slight mismatch of the stability of the extinction equilibrium (0, 0) in subfigures (a) and (c) is due to the multistable dynamics of our proposed model (8).

Fig 8 demonstrates the importance of the parameters *Ã*_{1}, *Ã*_{2} and *¾* for the enhancement of cooperation under the presence of mutation. The larger values of *Ã*_{1} and *Ã*_{2} facilitate the evolution of at least one species, and the positive mutation rate *¼* &#x; (0, 1] assures that species should mutate into the other. This mechanism will lead to the species&#x; coexistence in a major portion (yellow region) of the two-dimensional parameter spaces represented in Fig 8(a) and 8(b). The results are further validated by plotting the stability curve (white solid lines) below which the stationary point (0, 0) is locally stable. A higher mortality rate never entertains the evolution of both strategies; thus, we obtain the red region in Fig 8(c). This red region indicates the extinction equilibrium (0, 0). Once again, we plot the stability curve (solid white line) in Fig 8(c), above which both the species should be extinct as per our local stability analysis. All the subfigures are drawn with fixed initial condition (, ).

Open in a separate window

Fig 8

**Effect of bidirectional mutation on the emergent dynamics.**

The positive mutation rate *¼* &#x; (0, 1] allows the system (8) to settle into two stationary points depending on the choice of other parameters&#x; values. When the free space-induced benefits are small, all species are extinct, as depicted through subfigures (a-b). However, all strategies can coexist for suitable choices of other parameters. The red region reflects the extinction equilibrium, and the yellow region indicates the stable coexistence of cooperators and defectors. The solid white line represents the analytically derived stability curve below which the extinction equilibrium (0, 0) is locally stable in subfigures (a-b). The initial condition is chosen here as (*x*_{0}, *y*_{0}) = (, ). Other parameters&#x; values are kept fixed at *¾* = , *²* = , *·* = , *p* = , *Ã*_{2} = , *Ã*_{1} = , unless they are varied. The moderate choice of mortality rate *¾* &#x; (0, 1] allows the species&#x; coexistence in subfigure (c). Beyond a critical value of *¾*, both the cooperators as well as the defectors die, as reflected through the red region of subfigure (c).

We inspect the influence of different parameters on the constructed model (8) in Fig 9. We keep fixed all parameters&#x; values at *¾* = , *²* = , *¼* = , *p* = , *·* = , *Ã*_{1} = and *Ã*_{2} = , unless they are varied. Since free space provides additional advantage to the cooperators compared to the defectors as we choose *Ã*_{1} = *Ã*_{2} = , we have a defector-free society at least initially with *¼* = 0 in Fig 9(a). However, whenever the mutation rate *¼* becomes positive, each strategy can mutate into other. In this way, the density of the cooperators (red line) decreases, and the defectors&#x; density (blue line) increases. Eventually, both densities almost become identical for *¼* &#x; 1 &#x;. In of Fig 9(b), the vital role of *p* is investigated in the one-dimensional bifurcation diagram. As *p* &#x; 1 &#x;, the defectors are getting the upper hand over the cooperators as *p* indicates the probability of playing the PD game. PD game always provides additional assistance to the defectors. All these results obtained in Fig 9 are consistent with the social dilemmas considered for constructing the model (8). Fig 9(b) points out that lower values of *p* are always better for the maintenance of cooperation as small values of *p* indicate more rational people are playing the SD game and the SD game always favors the coexistence of both strategies. Increasing the temptation parameter *²* always uplifts the defectors&#x; fraction *y* (blue line). Thus, one needs to choose the value of *²* wisely so that we obtain a moderate range of *²* in Fig 9(c) allowing the coexistence of both strategies. When *²* is small, we found the cooperators&#x; density *x* (red line) dominates the defectors&#x; fraction *y* (blue line). However for *²* , *y* is larger than *x*. Note that the results may alter for a different choice of initial condition as the system 8 is multistable. We plot Fig 9 with fixed initial condition (, ) and the code to generate this figure is freely available at [70]. Increasing the parameter *·* &#x; (0, 1) can provide extra benefits to the defectors. Thus, the rate of increment of *y* is slightly better than that of *x* in subfigure (d) of Fig 9. However, other parameters&#x; values are also crucial for stabilizing the competitive communities. Hence, *y* (blue line) remains always lower than *x* (red line) in Fig 9(d) for our chosen parameters&#x; values.

Open in a separate window

Fig 9

**Effects of various parameters in spreading of cooperative behavior.**

Subfigure (a) reveals two contrasting scenarios, where the cooperators decrease with the increment of the mutation rate *¼*. However, the defectors&#x; fraction gradually increases with increasing *¼*. In the absence of mutation (i.e., for *¼* = 0), all defectors are extinct as the parameters values chosen to produce this subfigure is *¾* = , *²* = , *·* = , *p* = , *Ã*_{1} = and *Ã*_{2} = Thus, the free space induced benefit towards cooperator *Ã*_{1} is greater than that of defector *Ã*_{2}. Subfigure (b) shows that both species&#x; densities gradually increase with the probability *p* of playing the PD game. As the probability of playing the PD game increases, the increment of the fraction of defectors (blue line) gets over in margin with respect to the fraction of the cooperators (red line) after *p* = depending on other parameters&#x; values. This fact is obvious from the view that a higher possibility of playing the PD game gets a more beneficial ambiance for defectors. Subfigure (c) indicates the same phenomena as the value of the temptation payoff parameter *²* increases, both the fractions escalate. Whereas at *²* = , both these fractions get the same value. After that value of *²*, the value of the fraction of defector (blue line) becomes more than that of the cooperators (red line). Subfigure (d) also highlights the increment of both the fractions as the punishment payoff parameter *·* increases. But up to the highest received payoff value of *·*, the value of the fraction of defectors can never exceed the cooperators&#x; fraction because the free space provides more benefit to the cooperators than the defectors. Subfigure (e) demonstrates that when *Ã*_{1} &#x; [0, ], both these fractions are extinct, as the mortality rate *¾* has more value than these free space-induced benefits. But as the value of *Ã*_{1} increases, the cooperator gets more advantages than defectors. Subfigure (f) shows that both the fractions increase with the enhancement of *Ã*_{2}. Up to *Ã*_{2} = , the cooperators&#x; density acquires more value than *y*, as the value of *Ã*_{1} = is taken sufficiently high. After *Ã*_{2} = , *y* becomes higher than *x*. All the subfigures are drawn by iterating the system (8) with 20 × 10^{5} times with fixed integrating step length *´t* = and fixed initial condition (*x*_{0}, *y*_{0}) = (, ). The last point at 20 × 10^{5} iteration is collected to identify the asymptotic state of the proposed model. Other parameters are kept fixed at *¾* = , *²* = , *¼* = , *p* = , *·* = , *Ã*_{1} = and *Ã*_{2} = , unless they are varied.

Since the mutation rate *¼* is 50%, we can obtain only two stationary points of the model (8). Thus, initially, for a smaller choice of *Ã*_{1}, extinction of both species prevails in Fig 9(e). Nevertheless, with an increment of *Ã*_{1}, cooperators are getting further assistance from free space, and hence *x* becomes positive beyond *Ã*_{1} = Since *¼* = , 50% of these cooperators mutate bidirectionally into the defectors, and thus, we have positive *y* too in that range of *Ã*_{1}. Since *Ã*_{2} = *¾* = , thus we get the stabilization of (0, 0) within the interval *Ã*_{1} &#x; [0, ]. In fact, the cooperators&#x; density *x* is slightly better than that of *y* for the higher values of *Ã*_{1} in Fig 9(e). The choice of the initial condition and other parameters&#x; values are vital in obtaining all these results. Similar observation can be found in Fig 9(f), where we examine the role of *Ã*_{2} in our proposed model (8). With increasing *Ã*_{2}, the defectors will be benefited from free space, and hence, *y* (blue line) dominates *x* (red line) for larger values of *Ã*_{2}. But since *Ã*_{1} is taken sufficiently large in this Fig 9(f), we initially have a small portion for smaller values of *Ã*_{2} where *x* *y*. Fig 9 discloses suitable choices of all parameters&#x; values not only entertain the coexistence of both strategies but also may promote the evolution of cooperation.

### 4 Concluding remarks

Our mutation-induced model provides widespread coexistence of both strategies under suitable choices of parameters&#x; values. Such stable persistence of ecological communities strengthens the theory of concurrency. Our thorough analysis with several numerical simulations enhances our understanding of the mechanisms that drive the survival of cooperative behavior in the multigames consisting of both the Prisoner&#x;s Dilemma and the Snowdrift games. The inclusion of altruistic behavior of free space along with the mutation in our proposed evolutionary model seems to be a natural course of action as observed in many realistic settings. Our findings clearly demonstrate that there exists an optimal probability of playing each game so that people are more likely to cooperate in such circumstances. Throughout the study, we have pointed out the positive impact of diverse factors (parameters) on significant improvement of cooperation. We have shown that the selfless contribution of free space promotes the coexistence of all strategies efficiently, even in the absence of mutation.

In summary, our research indicates the proposed model (8) may possess four different stationary points in the absence of mutation. The exciting feature of this model is that the system never allows settling into a cooperator-free state in the lack of free space-induced benefits and mutation rate. This precise result is consistent with the chosen games as both the PD and SD games never encourage a cooperator-free society in the usual scenario. The merging of two games with different outcomes provides a more realistic representation of the concept of opinion formation. However, the numerical results presented here are highly sensitive to the variation of initial conditions. The positive invariance and boundedness of the model are analyzed too in this article. We have shown the viable choice of bidirectional mutation allows the system to switch between only two stationary states. Either all people will die, or both the strategies coexist in the eco-evolutionary model with mutation. This is an interesting angle of our research as our model brings forth stable biodiversity in the form of a heterogeneous population (mixed cooperator-defector state). Note that we have only focused on the equilibria of the corresponding dynamical system (8) throughout the article so that we can relate those stationary states from the game-theoretical point of view. Our simple model with mutation sheds light on how cooperation emerges in a complex society. The presented insightful results attest that altruistic behavior and mutation are advantageous for the spontaneous maintenance of biodiversity. We believe the presented results may help us understand the mechanism behind the coexistence of competing species through the co-evolution of both strategies.

### 5 Discussion: Limitation &#x; future perspectives

Lastly, we discuss challenges related to our work that merits further investigation. Even though it is almost impossible to incorporate all complex, realistic relationships among social creatures using minimal modeling, it is nevertheless absolutely essential to investigate such models using elementary mathematical principles. In fact, the literature already provides a number of excellent models that offer a stimulating starting point for exploring many practical situations. Chen et al. [71] proposed an elegant model aimed at understanding the co-evolutionary outcome of the strategies of the traffic management department, drivers, and pedestrians. Reference [72] proposes a hybrid machine learning model to predict traffic accidents. A novel autonomous system without any stationary points may offer hidden attractors, as revealed in Ref. [73]. It is almost impossible to incorporate even a partial list of relevant references to emphasize the importance of minimal modeling here. Motivated by all these works, we present a new model that can offer several valuable insights into the evolution of cooperation and altruism. However, our model is far from perfect. Instead of considering two-person interactions, the introduction of group interactions may yield a deeper understanding of decision-making. Interdependent networks may be an excellent choice to study the evolution of cooperation, as reflected through Refs. [20, 31, 62, 74]. This remains a promising future research generalization, including interdependent networks which may yield several complex dynamical behaviors other than stable equilibrium states.

One should observe that our approach works absolutely fine with other social dilemmas. For any two two-person games of the form , the system constructed with the help of our policy always remains bounded within the closed interval [0, 1]. However, the stationary points alter due to the change in the mathematical model. We have verified this boundedness by considering the snowdrift and the harmony game [75, 76] (results are not shown here). In fact, on the generality of our results, we can comment that a similar impact of mutation and free space can also be expected for other two-person games. It will be interesting to examine how the results may vary with other games like public goods games, rock paper scissor, etc. Besides introducing other realistic scenarios like apology and forgiveness [77], intention recognition [78], delay [79] etc., one may advance our understanding of the roots of cooperation in social and biological systems. Moreover, the impact of higher-order interactions [80] among the agents on the coevolution [81] of cooperation and synchronization [82&#x;84] in a coupled network remains largely unexplored. The interdisciplinary researchers of complex systems can pay attention to this exciting topic for future research.

### Acknowledgments

The authors gratefully acknowledge the anonymous referees for their careful reading, insightful suggestions, and valuable remarks that helped considerably improve and readability of the manuscript. S.R. wants to convey his sincere gratitude to Gourab Kumar Sar, Md Sayeed Anwar, Srilena Kundu, and Subrata Ghosh of the Indian Statistical Institute, Kolkata, for their helpful discussion and suggestions.

### Funding Statement

M.P. was supported by the Javna Agencija za Raziskovalno Dejavnost RS (Grant Nos. P and J).

### Data Availability

All relevant data are within the paper.

### References

1. Darwin C. On the origin of species: A facsimile of the first edition. Harvard University Press; [Google Scholar]

2. West SA, Griffin AS, Gardner A. Social semantics: altruism, cooperation, mutualism, strong reciprocity and group selection. Journal of Evolutionary Biology. ;20(2)&#x; doi: /jx [PubMed] [CrossRef] [Google Scholar]

3. Weibull JW. Evolutionary game theory. MIT press; [Google Scholar]

4. Nowak MA. Evolutionary dynamics: exploring the equations of life. Harvard university press; [Google Scholar]

5. Sigmund K. The calculus of selfishness. Princeton University Press; [Google Scholar]

6. Hofbauer J, Sigmund K, et al. Evolutionary games and population dynamics. Cambridge university press; [Google Scholar]

7. Allen B, Lippner G, Chen YT, Fotouhi B, Momeni N, Yau ST, et al. Evolutionary dynamics on any population structure. Nature. ;()&#x; doi: /nature [PubMed] [CrossRef] [Google Scholar]

8. Nag Chowdhury S, Kundu S, Banerjee J, Perc M, Ghosh D. Eco-evolutionary dynamics of cooperation in the presence of policing. Journal of Theoretical Biology. ; doi: /sprers.eu [PubMed] [CrossRef] [Google Scholar]

9. Pennisi E. How did cooperative behavior evolve?Science. ;()&#x; doi: /science [PubMed] [CrossRef] [Google Scholar]

Fotouhi B, Momeni N, Allen B, Nowak MA. Evolution of cooperation on large networks with community structure. Journal of the Royal Society Interface. ;16() doi: /rsif

I went back to an earlier iteration without the animation and tried merging in my character again, but Max crashed instantly, and was unable to autosave the file, siting a 'Memory Error'.

I did a quick search and the only other posting with similar a question is in the Maya forum.

Is the problem simply that my computer is not powerful enough?

Edit:

When I open the Material Editor now the following messages:

SCEN warn switching ambient occlusion off, since trace is off

MSG error cannot start thread 4: Not enough storage is available to process this command.

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

SCEN warn switching ambient occlusion off, since trace is off

RC warn switching ambient occlusion off, since trace is off

Also, I've just noticed that after adding the animation this file has gone from 9, KB to 14, KB

Max , Operating System: Windows 7 bit

Video Card: 1 Gigabyte ATI Radeon HD

CPU: Intel Core 2 Duo E GHZ, Ram: 6GB

Hard Drive: 1 Terrabyte Western Digital 32mb cache

## Importing Objects - mental ray error

Thanks Joojaa,

Not sure how to reset global render settings..

The architectural file is only 7mb, and renders fine in mental ray, although after importing a car that is only mb, I get these error messages and a partial render in mental ray, sayin 'out of memory'. This pc has 4gb of RAM. and the car has no attached map files.

Somethins kooky..

mental ray for Maya

mental ray: version , 22 Nov , revision

mental ray: mental ray for Maya - startup done

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

API warn type conflict: "input" is type vector, "maya_ramp" returns type color

API warn type conflict: "specularColor" is type color, "maya_multdiv" returns type vector

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

API warn type conflict: "input" is type vector, "maya_ramp" returns type color

API warn type conflict: "specularColor" is type color, "maya_multdiv" returns type vector

API warn type conflict: "specularColor" is type color, "maya_multdiv" returns type vector

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

mental ray: got 8 satellite CPUs.

API warn while defining light "sunDirectionShape": light has no light shader

MSG error cannot start thread 6: Not enough storage is available to process this command.

MSG error cannot start thread 5: Not enough storage is available to process this command.

MSG error cannot start thread 7: Not enough storage is available to process this command.

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

mental ray: out of memory

MEM fatal can't allocate bytes.

MEM fatal can't allocate bytes.

If i try the maya renderer, it locks up and then gives this below

Rendering using 2 threads.

Automatic near/far clipping values: ,

Starting Rendering C:/Users/David/Documents/maya/projects/default/images/tmp/Genesis iff.

Constructing shading groups.

Rendering current frame.

Failed to allocate space for image.

====================================

Cause of memory exception

====================================

Mb Free Memory

Mb Free Swap

Mb Size of alloc

Mb Low Memory Threshold

====================================

Memory use when exception was thrown

====================================

====================================

Page faults

Mb Max resident size

Mb Peak total size(Estimated)

Mb Peak arena size

====================================

Mb Current

Mb POLY_DRAW_CACHE_DATA

Mb Pixel Map

Mb Render Geometry Arena

Mb arguments

Mb AG

Mb MEL

Mb File Texture Mipmaps

Mb Render Cache

Mb Arrays

Mb Transforms

Mb Data Blocks

Mb Ray Tracer

Mb NURBS Surface Shapes

Mb NURBS AG

Mb Nurbs Tesselator Vertex

Mb POLY_DRAW_CACHE_STATIC_DATA

Mb NURBS Geometry Cache

====================================

## 0 Comments