On some theory of monostable and bistable pure birth-jump integro-differential equations

We study an integral-differential equation that models a pure birth-jump process, where birth and dispersal cannot be decoupled. A case has been made that these processes are more suitable for phenomena such as plant dynamics, fire propagation, and cancer cell dynamics. We contrast the dynamics of this equation with those of the classical reaction-diffusion equation, where the reaction term models either logistic growth or a strong Allee effect. Recent evidence of an Allee effect has been found in plant dynamics during the germination process (due to seed predation) but not in the generation of seeds. This motivates where the Allee effect is included in our model. We prove the global existence and uniqueness of solutions with bounded initial data and analyze some properties of the solutions. Additionally, we prove results related to the persistence or extinction of a species, which are analogous to those of the classical reaction-diffusion equation. A key finding is that in some cases a population which is initially below the Allee threshold in some area, even if small, will actually survive. This is in contrast to solutions of the classical reaction-diffusion with the same initial data. Another difference of note is the lack of regularization and an infinite number of discontinuous equilibrium solutions to the birth-jump model.     

Density-dependent selection migration model with non-monotone fitness functions

This paper studies the classical single locus, diallelic selection model with diffusion for a continuously reproducing population. The phase variables are population density and allele frequency (or allele density). The genotype fitness depend only on population density but include one-hump functions of the density variable. With mild assumptions on genotype fitnesses, we study the geometry of the nullclines and the asymptotic behavior of solutions of the selection model without diffusion. For the diffusion model with zero Neumann boundary conditions, we use this geometric information to show that if the initial data satisfy certain conditions then the corresponding solution to the reaction-diffusion equation converges to the spatially constant stable equilibrium which is closest to the initial data.Research partially supported by NSF grant DMS-8920597Research supported by funds provided by the USDA-Forest Service, Southeastern Forest Experiment Station, Pioneering (Population Genetics of Forest Trees) Research Unit, Raleigh, North Carolina     

一类具时滞的禽流感模型

针对具扩散和时滞的SI-SIR传染病模型,用特征分析和Lyapunov泛函方法研究了相应的具齐次Neumann边界条件反应扩散方程组解的渐近性质.最后给出数值模拟来说明如果染病鸟类的接触率和染病人类的接触率小,那么全系统的无病平衡点是全局渐近稳定的;但当染病鸟类的接触率大或者和染病人类的接触率大时,变异的禽流感将在人类中扩散.     

Computational methods for diffusion-influenced biochemical reactions

MOTIVATION: We compare stochastic computational methods accounting for space and discrete nature of reactants in biochemical systems. Implementations based on Brownian dynamics (BD) and the reaction-diffusion master equation are applied to a simplified gene expression model and to a signal transduction pathway in Escherichia coli. RESULTS: In the regime where the number of molecules is small and reactions are diffusion-limited predicted fluctuations in the product number vary between the methods, while the average is the same. Computational approaches at the level of the reaction-diffusion master equation compute the same fluctuations as the reference result obtained from the particle-based method if the size of the sub-volumes is comparable to the diameter of reactants. Using numerical simulations of reversible binding of a pair of molecules we argue that the disagreement in predicted fluctuations is due to different modeling of inter-arrival times between reaction events. Simulations for a more complex biological study show that the different approaches lead to different results due to modeling issues. Finally, we present the physical assumptions behind the mesoscopic models for the reaction-diffusion systems. AVAILABILITY: Input files for the simulations and the source code of GMP can be found under the following address: http://www.cwi.nl/projects/sic/bioinformatics2007/     

A new 3D mass diffusion�Creaction model in the neuromuscular junction

A three-dimensional model of the reaction-diffusion processes of a neurotransmitter and its ligand receptor in a disk shaped volume is proposed which represents the transmission process of acetylcholine in the synaptic cleft in the neuromuscular junction. The behavior of the reaction-diffusion system is described by a three-dimensional diffusion equation with nonlinear reaction terms due to the rate processes of acetylcholine with the receptor. A new stable and accurate numerical method is used to solve the equations with Neumann boundaries in cylindrical coordinates. The simulation analysis agrees with experimental measurements of end-plate current, and agrees well with the results of the conformational state of the acetylcholine receptor as a function of time and acetylcholine concentration of earlier investigations with a smaller error compared to experiments. Asymmetric emission of acetylcholine in the synaptic cleft and the subsequent effects on open receptor population is simulated. Sensitivity of the open receptor dynamics to the changes in the diffusion parameters and neuromuscular junction volume is investigated. The effects of anisotropic diffusion and non-symmetric emission of transmitter at the presynaptic membrane is simulated.     

Speeds of invasion in a model with strong or weak Allee effects

We study an invasion model based on a reaction-diffusion equation with an Allee effect. We use a special, piecewise-linear, population growth rate. This function allows us to obtain traveling wave solutions and to compute wave speeds for a full range of Allee effects, including weak Allee effects. Some investigators claim that linearization fails to give the correct speed of invasion if there is an Allee effect. We show that the minimum speed for a sufficiently weak Allee may, in fact, be the same as that derived by means of linearization.     

Simple and complex spatiotemporal structures in a glycolytic allosteric enzyme model

Pattern formation in glycolysis is studied with a classical reaction-diffusion allosteric enzyme model. It is found that, similar to recent experimental reports in the yeast extracts, a small magnitude local perturbation can induce transient target waves in a two dimensional oscillatory medium. An above threshold stimulation generates target waves which eventually evolve into spatiotemporal chaos upon collisions with the boundary or other wave activities. Detailed simulation studies show that the studied simple glycolytic reaction-diffusion model can support three types of spatiotemporal behaviors which are independent of the boundary conditions: (1) a spatially uniform stable steady state, (2) periodic global oscillations and (3) spatiotemporal chaos.     

Modelling the impact of an invasive insect via reaction-diffusion

An exotic, specialist seed chalcid, Megastigmus schimitscheki, has been introduced along with its cedar host seeds from Turkey to southeastern France during the early 1990s. It is now expanding in plantations of Atlas Cedar (Cedrus atlantica). We propose a model to predict the expansion and impact of this insect. This model couples a time-discrete equation for the ovo-larval stage with a two-dimensional reaction-diffusion equation for the adult stage, through a formula linking the solution of the reaction-diffusion equation to a seed attack rate. Two main diffusion operators, of Fokker-Planck and Fickian types, are tested. We show that taking account of the dependence of the insect mobility with respect to spatial heterogeneity, and choosing the appropriate diffusion operator, are critical factors for obtaining good predictions.     

Persistence in reaction diffusion models with weak allee effect

We study the positive steady state distributions and dynamical behavior of reaction-diffusion equation with weak Allee effect type growth, in which the growth rate per capita is not monotonic as in logistic type, and the habitat is assumed to be a heterogeneous bounded region. The existence of multiple steady states is shown, and the global bifurcation diagrams are obtained. Results are applied to a reaction-diffusion model with type II functional response, and also a model with density-dependent diffusion of animal aggregation. J. S. is partially supported by United States NSF grants DMS-0314736 and EF-0436318, College of William and Mary summer grants, and a grant from Science Council of Heilongjiang Province, China.     

Spontaneous separation of bi-stable biochemical systems into spatial domains of opposite phases

Bi-stable chemical systems are the basic building blocks for intracellular memory and cell fate decision circuits. These circuits are built from molecules, which are present at low copy numbers and are slowly diffusing in complex intracellular geometries. The stochastic reaction-diffusion kinetics of a double-negative feedback system and a MAPK phosphorylation-dephosphorylation system is analysed with Monte-Carlo simulations of the reaction-diffusion master equation. The results show the geometry of intracellular reaction compartments to be important both for the duration and the locality of biochemical memory. Rules for when the systems lose global hysteresis by spontaneous separation into spatial domains in opposite phases are formulated in terms of geometrical constraints, diffusion rates and attractor escape times. The analysis is facilitated by a new efficient algorithm for exact sampling of the Markov process corresponding to the reaction-diffusion master equation.     

Non-existence of wave solutions for the class of reaction--diffusion equations given by the Volterra interacting-population equations with diffusion.

The system of equations is reduced to a single nonlinear parabolic equation on which a maximum principle can be used. It is then shown that the effect of uniform diffusion on the Volterra equations for any even number of interacting populations which have non-zero equilibrium values, is to damp out all spatial variations. The inclusion of population saturation terms is shown to enhance the damping process, as would be expected. The main consequence of the results is that such reaction-diffusion equations (given in section 5) cannot have physically realistic wave-like solutions, that is stable solutions, with non-negative values of the concentrations, which evolve from a time dependent solution.     

Drosophila segmentation: supercomputer simulation of prepattern hierarchy

Spontaneous prepattern formation in a two level hierarchy of reaction-diffusion systems is simulated in three space co-ordinates and time, mimicking gap gene and primary pair-rule gene expression. The model rests on the idea of Turing systems of the second kind, in which one prepattern generates position dependent rate constants for a subsequent reaction-diffusion system. Maternal genes are assumed responsible for setting up gradients from the anterior and posterior ends, one of which is needed to stabilize a double period prepattern suggested to underly the read out of the gap genes. The resulting double period pattern in turn stabilizes the next prepattern in the hierarchy, which has a short wavelength with many characteristics of the stripes seen in actual primary pair-rule gene expression. Without such hierarchical stabilization, reaction-diffusion mechanisms yield highly patchy short wave length patterns, and thus unreliable stripes. The model yields seven stable stripes located in the middle of the embryo, with the potential for additional expression near the poles, as observed experimentally. The model does not rely on specific chemical reaction kinetics, rather the effect is general to many such kinetic schemes. This makes it robust to parameter changes, and it has good potential for adapting to size and shape changes as well. The study thus suggests that the crucial organizing principle in early Drosophila embryogenesis is based on global field mechanisms, not on particular local interactions.     

Travelling front solutions of a nonlocal Fisher equation

We consider a scalar reaction-diffusion equation containing a nonlocal term (an integral convolution in space) of which Fisher‘s equation is a particular case. We consider travelling wavefront solutions connecting the two uniform states of the equation. We show that if the nonlocality is sufficiently weak in a certain sense then such travelling fronts exist. We also construct expressions for the front and its evolution from initial data, showing that the main difference between our front and that of Fisher‘s equation is that for sufficiently strong nonlocality our front is non-monotone and has a very prominent hump. Received: 8 August 1999 / Revised: 3 March 2000 / Published online: 14 September 2000     

Robust exponential stabilization of a class of delayed neural networks with reaction-diffusion terms

In this paper, the problem of global robust exponential stabilization for a class of neural networks with reaction-diffusion terms and time-varying delays which covers the Hopfield neural networks and cellular neural networks is investigated. A feedback control gain matrix is derived to achieve the global robust exponential stabilization of the neural networks by using the Lyapunov stability theory, and the stabilization condition can be verified if a certain Hamiltonian matrix with no eigenvalues on the imaginary axis. This condition can avoid solving an algebraic Riccati equation. Finally, a numerical simulation illustrates the effectiveness of the results.     

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.     

Transition from polar to duplicate patterns

The existence of symmetric nonuniform solutions in nonlinear reaction-diffusion systems is examined. In the first part of the paper, we establish systematically the bifurcation diagram of small amplitude solutions in the vicinity of the two first bifurcation points. It is shown that:

1. The system can adopt a stable symmetric solution (basic wave number 2) if the value of the bifurcation parameter is changed or if the initial polar structure (basic wave number 1) is sufficiently perturbed.
2. This behavior is independent of the particular reaction-diffusion model proposed and of the number of intermediate components (?2) involved.

In the second part of the paper, analogies are established between the possibilities offered by the bifurcation diagrams, involving only the two first primary branches, and the observation that in the early development of different organisms, appropriate experimental manipulations may switch the normal (polar) developmental pattern to a duplicate structure.     

Dynamics induced by delay in a nutrient–phytoplankton model with diffusion

In this paper, a nutrient–phytoplankton model described by a couple of reaction-diffusion equations with delay is studied analytically and numerically. The aim of this research is to provide an understanding of the impact of delay on the nutrient–phytoplankton dynamics. Significantly, the delay can not only induce instability of a positive equilibrium, but also promote the formation of patchiness (an irregular pattern) via Hopf bifurcation. However, if the delay does not exist, the positive equilibrium is always globally asymptotically stable when it exists. In addition, the numerical analysis indicates that the input rate and the loss rate of nutrient also play an important role in the growth of phytoplankton, which supports that eutrophic conditions may be a significant reason inducing phytoplankton blooms. Numerical results are consistent with the analytical results.     

一类自催化反应扩散模型正解的唯一性与稳定性

主要考虑了一类三分子自催化反应扩散系统．在齐次Dirichlet和Robin边界条件下,当反应率c适当小,系统没有共存态;当c适当大,系统至少有一个共存态;当c充分大,系统有唯一渐近稳定的共存态．特别地,在一维空间上共存态是唯一的．在齐次Neumann边界条件下系统是一个简单系统．     

Dissipative structures associated to certain reaction-diffusion equations in mathematical ecology

Sufficient conditions are determined for the existence of stable spatially heterogeneous solutions of the Lotka-Volterra many species equations extended to include spatial diffusion. Some properties of such solutions are obtained. A general definition is proposed for the concept of a dissipative structure associated with a given reaction-diffusion equation. Finally, an approximate solution is presented for two interacting species in one spatial dimension, although the question of stability for this example is left open. Supported in part by NSERC A-7667 to P. L. Antonelli and by a University of Alberta President's Fund grant to J. R. Royce.