In the animal world, collective action to shelter, protect and nourish requires the cooperation of group members. Among humans, many situations require the cooperation of more than two individuals simultaneously. Most of the relevant literature has focused on an extreme case, the N-person Prisoner's Dilemma. Here we introduce a model in which a threshold less than the total group is required to produce benefits, with increasing participation leading to increasing productivity. This model constitutes a generalization of the two-person stag hunt game to an N-person game. Both finite and infinite population models are studied. In infinite populations this leads to a rich dynamics that admits multiple equilibria. Scenarios of defector dominance, pure coordination or coexistence may arise simultaneously. On the other hand, whenever one takes into account that populations are finite and when their size is of the same order of magnitude as the group size, the evolutionary dynamics is profoundly affected: it may ultimately invert the direction of natural selection, compared with the infinite population limit.

In the animal world, performing a given task which is beneficial to an entire group requires the cooperation of several individuals of that group who often share the workload required to perform the task. The mathematical framework to study the dynamics of collective action is game theory. Here we study the evolutionary dynamics of cooperators and defectors in a population in which groups of individuals engage in N-person, non-excludable public goods games. We explore an N-person generalization of the well-known two-person snowdrift game. We discuss both the case of infinite and finite populations, taking explicitly into consideration the possible existence of a threshold above which collective action is materialized. Whereas in infinite populations, an N-person snowdrift game (NSG) leads to a stable coexistence between cooperators and defectors, the introduction of a threshold leads to the appearance of a new interior fixed point associated with a coordination threshold. The fingerprints of the stable and unstable interior fixed points still affect the evolutionary dynamics in finite populations, despite evolution leading the population inexorably to a monomorphic end-state. However, when the group size and population size become comparable, we find that spite sets in, rendering cooperation unfeasible.

No abstract

We study the large population limit of the Moran process, assuming weak-selection, and for different scalings. Depending on the particular choice of scalings, we obtain a continuous model that may highlight the genetic-drift (neutral evolution) or natural selection; for one precise scaling, both effects are present. For the scalings that take the genetic-drift into account, the continuous model is given by a singular diffusion equation, together with two conservation laws that are already present at the discrete level. For scalings that take into account only natural selection, we obtain a hyperbolic singular equation that embeds the Replicator Dynamics and satisfies only one conservation law. The derivation is made in two steps: a formal one, where the candidate limit model is obtained, and a rigorous one, where convergence of the probability density is proved. Additional results on the fixation probabilities are also presented.

ABSTRACT. We study a class of processes that are akin to the Wright-Fisher model, with transition probabilities weighted in terms of the frequency-dependent fitness of the population types. By considering an approximate weak formulation of the discrete problem, we are able to derive a corresponding continuous weak formulation for the probability density. Therefore, we obtain a family of partial differential equations (PDE) for the evolution of the probability density, and which will be an approximation of the discrete process in the joint large population, small time-steps and weak selection limit. If the fitness functions are sufficiently regular, we can recast the weak formulation in a more standard formulation, without any boundary conditions, but supplemented by a number of conservation laws. The equations in this family can be purely diffusive, purely hyperbolic or of convection-diffusion type, with frequency dependent convection. The particular outcome will depend on the assumed scalings. The diffusive equations are of the degenerate type; using a duality approach, we also obtain a frequency dependent version of the Kimura equation without any further assumptions. We also show that the convective approximation is related to the replicator dynamics and provide some estimate of how accurate is the convective approximation, with respect to the convective-diffusion approximation. In particular, we show that the mode, but not the expected value, of the probability distribution is modelled by the replicator dynamics. Some numerical simulations that illustrate the results are also presented. Wright-Fisher process and diffusion approximations and continuous limits and replicator equation 92D15 and 92D25 and 35K57 and 35K67 and 35L65

In vector-borne diseases, the human hosts' epidemiology often acts on a much slower time scales than the one of the mosquitos transmitting as a vector from human to human, due to their vastly different life cycles. We investigate in how far the fast time scale of the mosquito epidemiology can be slaved by the slower human epidemiology, so that for the understanding of human disease data mainly the dynamics of the human time scale is essential and only slightly perturbed by the mosquito dynamics.

We study the evolution of the probability density of an asexual, one locus population under natural selection and random evolution. This evolution is governed by a Fokker-Planck equation with degenerate coefficients on the boundaries, supplemented by a pair of conservation laws. It is readily shown that no classical or standard weak solution definition yields solvability of the problem. We provide an appropriate definition of weak solution for the problem, for which we show existence and uniqueness. The solution displays a very distinctive structure and, for large time, we show convergence to a unique stationary solution that turns out to be a singular measure supported at the endpoints. An exponential rate of convergence to this steady state is also proved.2000 Mathematics Subject Classification. Primary 95D15; Secondary 35K65. Key words and phrases. gene fixation, evolutionary dynamics, degenerate parabolic equations, boundarycoupled weak solutions.The authors want to thank Peter Markowich for helpful comments, and an anonymous referee for pointing out a mistake in an earlier version of Proposition 4 and for general comments that improved the presentation.

We present a derivation of the classical SIR model through a mean-field approximation from a discrete version of SIR. We then obtain a hyperbolic forward Kolmogorov equation, and show that its projected characteristics recover the standard SIR model. Moreover, we show that the long time limit of the evolution will be a Dirac measure. The exact position will depend on the well-know $R_0$ parameter, and it will be supported on the corresponding stable SIR equilibrium