Daphnias: from the individual based model to the large population equation

The class of deterministic 'Daphnia' models treated by Diekmann et al. (J Math Biol 61: 277-318, 2010) has a long history going back to Nisbet and Gurney (Theor Pop Biol 23: 114-135, 1983) and Diekmann et al. (Nieuw Archief voor Wiskunde 4: 82-109, 1984). In this note, we formulate the individual based models (IBM) supposedly underlying those deterministic models. The models treat the interaction between a general size-structured consumer population ('Daphnia') and an unstructured resource ('algae'). The discrete, size and age-structured Daphnia population changes through births and deaths of its individuals and throught their aging and growth. The birth and death rates depend on the sizes of the individuals and on the concentration of the algae. The latter is supposed to be a continuous variable with a deterministic dynamics that depends on the Daphnia population. In this model setting we prove that when the Daphnia population is large, the stochastic differential equation describing the IBM can be approximated by the delay equation featured in (Diekmann et al., l.c.)

Evolutionary branching in a stochastic population model with discrete mutational steps

Evolutionary branching is analysed in a stochastic, individual-based population model under mutation and selection. In such models, the common assumption is that individual reproduction and life career are characterised by values of a trait, and also by population sizes, and that mutations lead to small changes in trait value. Then, traditionally, the evolutionary dynamics is studied in the limit of vanishing mutational step sizes. In the present approach, small but non-negligible mutational steps are considered. By means of theoretical analysis in the limit of infinitely large populations, as well as computer simulations, we demonstrate how discrete mutational steps affect the patterns of evolutionary branching. We also argue that the average time to the first branching depends in a sensitive way on both mutational step size and population size.Comment: 12 pages, 8 figures. Revised versio

From supported membranes to tethered vesicles: lipid bilayers destabilisation at the main transition

We report results concerning the destabilisation of supported phospholipid bilayers in a well-defined geometry. When heating up supported phospholipid membranes deposited on highly hydrophilic glass slides from room temperature (i.e. with lipids in the gel phase), unbinding was observed around the main gel to fluid transition temperature of the lipids. It lead to the formation of relatively monodisperse vesicles, of which most remained tethered to the supported bilayer. We interpret these observations in terms of a sharp decrease of the bending rigidity modulus $\kappa$ in the transition region, combined with a weak initial adhesion energy. On the basis of scaling arguments, we show that our experimental findings are consistent with this hypothesis.Comment: 11 pages, 3 figure

Active Membrane Fluctuations Studied by Micropipet Aspiration

We present a detailed analysis of the micropipet experiments recently reported in J-B. Manneville et al., Phys. Rev. Lett. 82, 4356--4359 (1999), including a derivation of the expected behaviour of the membrane tension as a function of the areal strain in the case of an active membrane, i.e., containing a nonequilibrium noise source. We give a general expression, which takes into account the effect of active centers both directly on the membrane, and on the embedding fluid dynamics, keeping track of the coupling between the density of active centers and the membrane curvature. The data of the micropipet experiments are well reproduced by the new expressions. In particular, we show that a natural choice of the parameters quantifying the strength of the active noise explains both the large amplitude of the observed effects and its remarkable insensitivity to the active-center density in the investigated range. [Submitted to Phys Rev E, 22 March 2001]Comment: 14 pages, 5 encapsulated Postscript figure

COMPETITIVE OR WEAK COOPERATIVE STOCHASTIC LOTKA-VOLTERRA SYSTEMS CONDITIONED TO NON-EXTINCTION

International audienceWe are interested in the long time behavior of a two-type density-dependent biological population conditioned to non-extinction, in both cases of competition or weak cooperation between the two species. This population is described by a stochastic Lotka-Volterra system, obtained as limit of renormalized interacting birth and death processes. The weak cooperation assumption allows the system not to blow up. We study the existence and uniqueness of a quasi-stationary distribution, that is convergence to equilibrium conditioned to non extinction. To this aim we generalize in two-dimensions spectral tools developed for one-dimensional generalized Feller diffusion processes. The existence proof of a quasi-stationary distribution is reduced to the one for a $d$-dimensional Kolmogorov diffusion process under a symmetry assumption. The symmetry we need is satisfied under a local balance condition relying the ecological rates. A novelty is the outlined relation between the uniqueness of the quasi-stationary distribution and the ultracontractivity of the killed semi-group. By a comparison between the killing rates for the populations of each type and the one of the global population, we show that the quasi-stationary distribution can be either supported by individuals of one (the strongest one) type or supported by individuals of the two types. We thus highlight two different long time behaviors depending on the parameters of the model: either the model exhibits an intermediary time scale for which only one type (the dominant trait) is surviving, or there is a positive probability to have coexistence of the two species

Polymorphic evolution sequence and evolutionary branching

We are interested in the study of models describing the evolution of a polymorphic population with mutation and selection in the specific scales of the biological framework of adaptive dynamics. The population size is assumed to be large and the mutation rate small. We prove that under a good combination of these two scales, the population process is approximated in the long time scale of mutations by a Markov pure jump process describing the successive trait equilibria of the population. This process, which generalizes the so-called trait substitution sequence, is called polymorphic evolution sequence. Then we introduce a scaling of the size of mutations and we study the polymorphic evolution sequence in the limit of small mutations. From this study in the neighborhood of evolutionary singularities, we obtain a full mathematical justification of a heuristic criterion for the phenomenon of evolutionary branching. To this end we finely analyze the asymptotic behavior of 3-dimensional competitive Lotka-Volterra systems

Fluctuation spectrum of fluid membranes coupled to an elastic meshwork: jump of the effective surface tension at the mesh size

We identify a class of composite membranes: fluid bilayers coupled to an elastic meshwork, that are such that the meshwork's energy is a function $F_\mathrm{el}[A_\xi]$ \textit{not} of the real microscopic membrane area $A$, but of a \textit{smoothed} membrane's area $A_\xi$, which corresponds to the area of the membrane coarse-grained at the mesh size $\xi$. We show that the meshwork modifies the membrane tension $\sigma$ both below and above the scale $\xi$, inducing a tension-jump $\Delta\sigma=dF_\mathrm{el}/dA_\xi$. The predictions of our model account for the fluctuation spectrum of red blood cells membranes coupled to their cytoskeleton. Our results indicate that the cytoskeleton might be under extensional stress, which would provide a means to regulate available membrane area. We also predict an observable tension jump for membranes decorated with polymer "brushes"

A Semi-Lagrangian scheme for a modified version of the Hughes model for pedestrian flow

In this paper we present a Semi-Lagrangian scheme for a regularized version of the Hughes model for pedestrian flow. Hughes originally proposed a coupled nonlinear PDE system describing the evolution of a large pedestrian group trying to exit a domain as fast as possible. The original model corresponds to a system of a conservation law for the pedestrian density and an Eikonal equation to determine the weighted distance to the exit. We consider this model in presence of small diffusion and discuss the numerical analysis of the proposed Semi-Lagrangian scheme. Furthermore we illustrate the effect of small diffusion on the exit time with various numerical experiments

The frequency-dependent Wright-Fisher model: diffusive and non-diffusive approximations

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