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

On the wellposedness of some McKean models with moderated or singular diffusion coefficient

We investigate the well-posedness problem related to two models of nonlinear McKean Stochastic Differential Equations with some local interaction in the diffusion term. First, we revisit the case of the McKean-Vlasov dynamics with moderate interaction, previously studied by Meleard and Jourdain in [16], under slightly weaker assumptions, by showing the existence and uniqueness of a weak solution using a Sobolev regularity framework instead of a Holder one. Second, we study the construction of a Lagrangian Stochastic model endowed with a conditional McKean diffusion term in the velocity dynamics and a nondegenerate diffusion term in the position dynamics

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