Receptor-based models with hysteresis for pattern formation in hydra

In this paper, we propose a new receptor-based model for pattern formation and regulation in a fresh-water polyp, namely hydra. The model is defined in the form of a system of reaction-diffusion equations with zero-flux boundary conditions coupled with a system of ordinary differential equations. The production of diffusible biochemical molecules has a hysteretic dependence on the density of these molecules and is modeled by additional ordinary differential equations. We study the hysteresis-driven mechanism of pattern formation and we demonstrate the advantages and constraints of its ability to explain different aspects of pattern formation and regulation in hydra. The properties of the model demonstrate a range of stationary and oscillatory spatially heterogeneous patterns, arising from multiple spatially homogeneous steady states and switches in the production rates.     

Stochastic differential equation models for spatially distributed neurons and propagation of chaos for interacting systems.

Distribution or nuclear space-valued stochastic differential equations (SDEs) (diffusions as well as discontinuous equations) are discussed as stochastic models for the behavior of voltage potentials of spatially distributed neurons. A propagation of chaos result is obtained for an interacting system of Hilbert space-valued SDEs.     

Multiscale Modelling of Fluid and Drug Transport in Vascular Tumours

A model for fluid and drug transport through the leaky neovasculature and porous interstitium of a solid tumour is developed. The transport problems are posed on a micro-scale characterized by the inter-capillary distance, and the method of multiple scales is used to derive the continuum equations describing fluid and drug transport on the length scale of the tumour (under the assumption of a spatially periodic microstructure). The fluid equations comprise a double porous medium, with coupled Darcy flow through the interstitium and vasculature, whereas the drug equations comprise advection–reaction equations; in each case the dependence of the transport coefficients on the vascular geometry is determined by solving micro-scale cell problems.     

On the effects of spatial heterogeneity on the persistence of interacting species

 The dynamics of two interacting theoretical populations inhabiting a heterogeneous environment are modelled by a system of two weakly coupled reaction–diffusion equations having spatially dependent reaction terms. Longterm persistence of both populations is guaranteed by an invasibility condition, which is itself expressed via the signs of certain eigenvalues of related linear elliptic operators with spatially dependent lowest order coefficients. The effects of change in these coefficients upon the eigenvalues are here exploited to study the effects of spatial heterogeneity on the persistence of interacting species through two particular ecological topics of interest. The first concerns when the location of favorable hunting grounds within the overall environment does or does not affect the success of a predator in predator–prey models, while the second concerns cases of competition models in which the outcome of competition in a spatially varying environment differs from that which would be expected in a spatially homogeneous environment. Received: 9 June 1997     

A mathematical model for weed dispersal and control

Mathematical models for weed dispersal and control are developed, analyzed and numerically simulated. A model incorporating periodic control, e.g. herbicide application, is derived for a plant population in a spatially homogeneous setting. The model is extended to a spatially heterogeneous population where plant dispersal is incorporated. The dispersal and control model involves integrodifference equations, discrete in time and continuous in space. The models are analyzed to determine values of the control parameter that prevent weed spread. The effects of the control on travelling wave solutions are investigated numerically.     

Spatially Extended Relativistic Particles Out of Traveling Front Solutions of Sine-Gordon Equation in (1+2) Dimensions

Slower-than-light multi-front solutions of the Sine-Gordon in (1+2) dimensions, constructed through the Hirota algorithm, are mapped onto spatially localized structures, which emulate free, spatially extended, massive relativistic particles. A localized structure is an image of the junctions at which the fronts intersect. It propagates together with the multi-front solution at the velocity of the latter. The profile of the localized structure obeys the linear wave equation in (1+2) dimensions, to which a term that represents interaction with a slower-than-light, Sine-Gordon-multi-front solution has been added. This result can be also formulated in terms of a (1+2)-dimensional Lagrangian system, in which the Sine-Gordon and wave equations are coupled. Expanding the Euler-Lagrange equations in powers of the coupling constant, the zero-order part of the solution reproduces the (1+2)-dimensional Sine-Gordon fronts. The first-order part is the spatially localized structure. PACS: 02.30.Ik, 03.65.Pm, 05.45.Yv, 02.30.Ik.     

Modeling of Multifactor Taxis in a Predator–Prey System

Biophysics - A model of predator–prey dynamics in a spatially heterogeneous range is considered using a system of two nonlinear equations of the diffusion–advection reaction. The...     

Mechanisms for stabilisation and destabilisation of systems of reaction-diffusion equations

Potential mechanisms for stabilising and destabilising the spatially uniform steady states of systems of reaction-diffusion equations are examined. In the first instance the effect of introducing small periodic perturbations of the diffusion coefficients in a general system of reaction-diffusion equations is studied. Analytical results are proved for the case where the uniform steady state is marginally stable and demonstrate that the effect on the original system of such perturbations is one of stabilisation. Numerical simulations carried out on an ecological model of Levin and Segel (1976) confirm the analysis as well as extending it to the case where the perturbations are no longer small. Spatio-temporal delay is then introduced into the model. Analytical and numerical results are presented which show that the effect of the delay is to destabilise the original system.     

From infanticide to parental care: why spatial structure can help adults be good parents

We investigate the evolution of parental care and cannibalism in a spatially structured population where adults can either help or kill juveniles in their neighborhood. We show that spatial structure can reverse the selective pressures on adult behavior, leading to the evolution of parental care, whereas the nonspatial model predicts that cannibalism is the sole evolutionary outcome. Our analysis emphasizes that evolution of such spatially structured populations is best understood at the level of the cluster of invading mutants, and we define invasion fitness as the growth rate of that cluster. We derive an analytical expression for the selective pressures on the trait and show that relatedness and Hamilton's rule are recovered as emergent properties of the spatial ecological dynamics. When adults can also help other adults, the benefits to each class of recipients are weighted by the class reproductive value, a result consistent with that of other models of kin selection. Finally, we advocate a different approach to moment equations and argue that even though the development of moment closure approximations is a necessary line of research, much-needed ecological and evolutionary insight can be gained by studying the unclosed moment equations.     

Integro-differential equations and the stability of neural networks with dendritic structure

 We analyse the effects of dendritic structure on the stability of a recurrent neural network in terms of a set of coupled, non-linear Volterra integro-differential equations. These, which describe the dynamics of the somatic membrane potentials, are obtained by eliminating the dendritic potentials from the underlying compartmental model or cable equations. We then derive conditions for Turing-like instability as a precursor for pattern formation in a spatially organized network. These conditions depend on the spatial distribution of axo-dendritic connections across the network. Received: 2 November 1994/Accepted in revised form: 7 March 1995     

A singular perturbation approach to diffusion reaction equations containing a point source,with application to the hemolytic plaque assay

Summary Many cells secrete factors which diffuse and bind to receptors on neighboring cells. These processes can be described by a nonlinear diffusion equation with a point source and a spatially distributed binding reaction. We show via perturbation analysis how approximate solutions can be obtained for such equations when the binding reaction is fast compared to diffusive transport. We base our analysis on an example which is of great practical importance in immunology, the hemolytic plaque technique.     

Nonuniformity of the chemical composition of a capillary discharge plasma

A steady-state distribution of the concentration of two ion species in a capillary discharge plasma is studied using MHD equations for a plasma with a spatially nonuniform, time-dependent chemical composition. In our case, the set of equations is significantly simplified because of the steady-state character and symmetry of the problem. Even with such simplification, however, some results could be obtained only by numerical integration. The factors affecting the distribution of heavy ions are studied. It is shown that the distribution of the heavy impurity over the discharge cross section can be much more nonuniform than the distribution of the main component (hydrogen). A simple criterion for such a nonuniformity is obtained.     

Solitary wave and spatially locked solitary pattern in a chemical reaction system

We investigate the behavior of a one-dimensional two component dynamical system. The dynamical equations are obtained by extracting an essence out of equations which describe the behavior of a biochemical reaction catalyzed by an allosteric protein. The obtained dynamical equations are similar to van der Pol equations. The dynamical equations are solved numerically. In the continuous system, a solitary wave is found to occur in certain ranges of the parameter space. The condition of occurrence of the solitary wave is investigated. The solitary wave can be induced by various initial perturbations, including rectangular ones with space-wise length longer than a certain critical value. The property of the solitary wave is similar to that of the impulses in nervous systems. In the discrete system, a spatially locked solitary pattern is found to occur in certain ranges of the parameter space.     

The evolution of juvenile-adult interactions in populations structured in age and space

We study the evolution of a spatially structured population with two age classes using spatial moment equations. In the model, adults can either help juveniles by increasing their survival, or adopt a cannibalistic behaviour and consume juveniles. While cannibalism is the sole evolutionary outcome when the population is well-mixed, both cannibalism and parental care can be evolutionarily stable if the population is viscous. Our analysis allows us to make two main technical points. First, we present a method to define invasion fitness in class-structured viscous populations, which allows us to apply adaptive dynamics methodology. Second, we show that ordinary pair approximation introduces an important quantitative bias in the evolutionary model, even on random networks. We propose a correction to the ordinary pair approximation that yields quantitative accuracy, and discuss how the bias associated with this approach is precisely what allows us to identify subtle aspects associated with the evolutionary dynamics of spatially structured populations.     

Evolution of complex dynamics in spatially structured populations

Dynamics of populations depend on demographic parameters which may change during evolution. In simple ecological models given by one-dimensional difference equations, the evolution of demographic parameters generally leads to equilibrium population dynamics. Here we show that this is not true in spatially structured ecological models. Using a multi-patch metapopulation model, we study the evolutionary dynamics of phenotypes that differ both in their response to local crowding, i.e. in their competitive behaviour within a habitat, and in their rate of dispersal between habitats. Our simulation results show that evolution can favour phenotypes that have the intrinsic potential for very complex dynamics provided that the environment is spatially structured and temporally variable. These phenotypes owe their evolutionary persistence to their large dispersal rates. They typically coexist with phenotypes that have low dispersal rates and that exhibit equilibrium dynamics when alone. This coexistence is brought about through the phenomenon of evolutionary branching, during which an initially uniform population splits into the two phenotypic classes.     

Spatial patterns in population conflicts

We consider a dynamical model for evolutionary games, and enquire how the introduction of diffusion may lead to the formation of stationary spatially inhomogeneous solutions, that is patterns. For the model equations being used it is already known that if there is an evolutionarily stable strategy (ESS), then it is stable. Equilibrium solutions which are not ESS's and which are stable with respect to spatially constant perturbations may be unstable for certain choices of the dispersal rates. We prove by a bifurcation technique that under appropriate conditions the instability leads to patterns. Computations using a curve-following technique show that the bifurcations exhibit a rich structure with loops joined by symmetry-breaking branches.     

Genetic Drift in a Cline

A model is developed of genetic drift in a cline maintained by spatially varying natural selection and local dispersal of individuals. The model is analyzed by an approximation scheme which is valid for weak selection and small migration rates. The results, which are based on numerical iterations of the approximate equations, are that the cline is less steep than predicted on the basis of the deterministic theory but that for weak selection the correlation between random fluctuations in neighboring colonies is approximately the same as in models of migration and drift in the absence of selection.     

Instability of non-constant equilibrium solutions of a system of competition-diffusion equations

The system of interaction-diffusion equations describing competition between two species is investigated. By using a version of the Perron-Frobenius theorem of positive matrices generalized to function spaces, it is proved that any non-constant equilibrium solution of the system is unstable both under Neumann boundary conditions (for the rectangular parallelepiped domain) and under periodic conditions. It is conjectured that this result extends to convex domains, and that the simple interaction-diffusion model cannot explain spatially segregated distributions of two competing species in such domains.