Periodic and traveling wave solutions to Volterra-Lotka equations with diffusion

Analytic and numerical solutions to two coupled nonlinear diffusion equations are studied. They are the modified equations of Volterra and Lotka for the spatially stratified predatorprey population model. In a bounded domain with the reflecting boundary, equilibrium, stability, and transition to time-periodic solutions are analyzed. For a wide class of initial states, the solutions to the initial boundary-value problem evolve into their corresponding stable, space-homogeneous, periodic oscillations. In an unbounded domain, a family of traveling wave solutions is found for certain exponential, initial distributions in the limit as the diffusion coefficientv ₁ of the prey tends to zero. In the presence of both diffusions, the results of a numerical simulation to an initial-value problem showed the rapid formation of the Pursuit-Evasion Waves whose speed of propagation and amplitudes increase with the diffusion coefficientv ₁. Presented at the 1974 SIAM Fall Meeting.     

Travelling wave solutions of diffusive Lotka-Volterra equations

We establish the existence of travelling wave solutions for two reaction diffusion systems based on the Lotka-Volterra model for predator and prey interactions. For simplicity, we consider only 1 space dimension. The waves are of transition front type, analogous to the travelling wave solutions discussed by Fisher and Kolmogorov et al. for a scalar reaction diffusion equation. The waves discussed here are not necessarily monotone. For any speed c there is a travelling wave solution of transition front type. For one of the systems discussed here, there is a distinguished speed c^* dividing the waves into two types, waves of speed c < c^* being one type, waves of speed c ? c^* being of the other type. We present numerical evidence that for this system the wave of speed c^* is stable, and that c^* is an asymptotic speed of propagation in some sense. For the other system, waves of all speeds are in some sense stable. The proof of existence uses a shooting argument and a Lyapunov function. We also discuss some possible biological implications of the existence of these waves.     

Minimum domains for spatial patterns in a class of reaction diffusion equations

We study a general class of scalar reaction/interacting population diffusion equations in two space dimensions: convective terms, due to wind, are included. We consider boundary conditions which include a measure of the hostility to the species in the exterior of the domain. The main point of the paper is to obtain estimates for the minimum domain size which can sustain spatially heterogeneous structures and indicate the type of patterns which could appear.     

Traveling wave solutions of a reaction diffusion model for competing pioneer and climax species

Presented is a reaction-diffusion model for the interaction of pioneer and climax species. For certain parameters the system exhibits bistability and traveling wave solutions. Specifically, we show that when the climax species diffuses at a slow rate there are traveling wave solutions which correspond to extinction waves of either the pioneer or climax species. A leading order analysis is used in the one-dimensional spatial case to estimate the wave speed sign that determines which species becomes extinct. Results of these analyses are then compared to numerical simulations of wave front propagation for the model on one and two-dimensional spatial domains. A simple mechanism for harvesting is also introduced.     

Square-root models for the Volterra equations and the explicit solution of these models

Volterra's (1926) equations for competition and predator-prey interactions are modified by introduction of root terms. A critical comparison with the original equations shows that the dynamic properties of the systems remain essentially alike, while the modification allows for explicit solution of the differential equations. Detailed solutions and numerical examples are given.     

Numerical inversion of the Perrin equations for rotational and translational diffusion constants by iterative techniques.

An iterative numerical technique is presented which allows the semiaxes for prolate and oblate ellipsoids to be determined from the Perrin equations for rotational and translational diffusion constants. The use of this inversion technique is illustrated by application to the proteins: lysozyme, bovine serum albumin, human transferrin, and bovine rhodopsin solubilized in digitonin.     

Existence and uniqueness of a sharp travelling wave in degenerate non-linear diffusion Fisher-KPP equations

In this paper we use a dynamical systems approach to prove the existence of a unique critical value c ^* of the speed c for which the degenerate density-dependent diffusion equation u _ct = [D(u)u _x ] _x + g(u) has: 1. no travelling wave solutions for 0 < c < c ^*, 2. a travelling wave solution u(x, t) = (x - c ^* t) of sharp type satisfying (– ) = 1, () = 0 ^*; '(^*–) = – c ^*/D'(0), '(^*+) = 0 and 3. a continuum of travelling wave solutions of monotone decreasing front type for each c > c ^*. These fronts satisfy the boundary conditions (– ) = 1, '(– ) = (+ ) = '(+ ) = 0. We illustrate our analytical results with some numerical solutions.     

On travelling wave solutions in a model for the Belousov-Zhabotinskii reaction.

DNA fragments partially digested by a 3′- or 5′-specific nuclease to produce single chain ends of opposite polarity will form a ring if the ends contain complementary sequences and are allowed to anneal. The frequency of rings can then be used as an assay to determine where and how identical repetitious sequences are arranged in the DNA. Thomas et al. (1973b) showed that all eucaryote chromosomes studied contain similar if not identical repetitious sequences clustered into regions called g-regions. To account for the observed ring frequency under different experimental conditions Thomas, Zimm &; Dancis (1973c) derived equations for two possible models of g-region organization. In the pure tandem model, the repetitious sequences are contiguous and occupy the entire g-region. In the intermittent repetition model, the repetitious sequences are simple copolymers and are irregularly arranged among non-repetitious sequences which are heterogenous in length. In the present paper, the results of Thomas et al. (1973c) are extended to cover the fractional tandem model. In this model, adjacent repetitious sequences are separated by non-repetitious sequences of uniform length. In addition, the equations for both the pure tandem and intermittent repetition models are shown to be special cases of the fractional tandem model but not vice versa.The capabilities and limitations of an analysis of ring formation are demonstrated using data from Drosophila. Although it is not possible to rule out any of the three models, the analysis can limit the ranges of the parameters describing each of the models that are consistent with the data. Previous conclusions that the data could only be explained by a pure tandem model which lacks any intervening unique sequences (Bick, Huang &; Thomas, 1973; Thomas et al., 1973b), are shown to be incorrect, in part because the equations for the fractional tandem model had not then been derived. Thus ring theory equations can be used to show the presence of clusters of similar if not identical sequences from ring-forming experiments, but they may not be able to determine the exact spacing and arrangement of these sequences within the clusters.     

Isolated periodic solutions of the Hodgkin-Huxley equations

Periodic solutions of the current clamped Hodgkin-Huxley equations (Hodgkin & Huxley, 1952 J. Physiol. 117, 500) that arise by degenerate Hopf bifurcation were studied recently by Labouriau (1985 SIAM J. Math. Anal. 16, 1121, 1987 Degenerate Hopf Bifurcation and Nerve Impulse (Part II), in press). Two parameters, temperature T and sodium conductance gNa were varied from the original values obtained by Hodgkin & Huxley. Labouriau's work proved the existence of small amplitude periodic solution branches that do not connect locally to the stationary solution branch, and had not been previously computed. In this paper we compute these solution branches globally. We find families of isolas of periodic solutions (i.e. branches not connected to the stationary branch). For values of gNa in the range measured by Hodgkin & Huxley, and for physically reasonable temperatures, there are isolas containing orbitally asymptotically stable solutions. The presence of isolas of periodic solutions suggests that in certain current space clamped membrane experiments, action potentials could be observed even though the stationary state is stable for all current stimuli. Once produced, such action potentials will disappear suddenly if the current stimulus is either increased or decreased past certain values. Under some conditions, "jumping" between action potentials of different amplitudes might be observed.     

A relaxation wave solution of the FitzHugh-Nagumo equations

For the system of FitzHugh-Nagumo equations and certain values of parameters, there exists a relaxation wave solution: a moving spatially periodic structure for which one of the dependent variables has intervals of slow and fast changes within its period. An asymptotic expansion for the solution is constructed (as a power series in the small parameter ) consisting of regular terms and transition layers. Along with the asymptotic expansion, the velocity of the structure, its period, and the distance between jumps of one of the functions are determined.Visiting from the Moscow State University, Russia.     

Generalization of the Spiegler-Kedem-Katchalsky frictional model equations of the transmembrane transport for multicomponent non-electrolyte solutions.

The Spiegler-Kedem-Katchalsky frictional model equations of the transmembrane transport for systems containing n-component, non-ionic solutions is presented. The frictional interpretation of the phenomenological coefficients of membrane and the expressions connecting the practical coefficients (Lp, sigma i, omega ij) with frictional coefficients (fij) are presented.     

The uniqueness of wave solutions for the deterministic non-reducible n-type epidemic

In a recent paper, [8], we investigated the existence of wave solutions for a model of the deterministic non-reducible n-type epidemic. In this paper we first prove two properties left as an open question in that paper. The uniqueness of the wave solutions at all speeds for which a wave solution exists is then established. Only an exceptional case is not covered.     

Polarized fluorescence photobleaching recovery for measuring rotational diffusion in solutions and membranes.

A variation of fluorescence photobleaching recovery (FPR) suitable for measuring the rate of rotational molecular diffusion in solution and cell membranes is presented in theory and experimental practice for epi-illumination microscopy. In this technique, a brief flash of polarized laser light creates an anisotropic distribution of unbleached fluorophores which relaxes by rotational diffusion, leading to a time-dependent postbleach fluorescence. Polarized FPR (PFPR) is applicable to any time scales from seconds to microseconds. However, at fast (microsecond) time scales, a partial recovery independent of molecular orientation tends to obscure rotational effects. The theory here presents a method for overcoming this reversible photobleaching, and includes explicit results for practical geometries, fast wobble of fluorophores, and arbitrary bleaching depth. This variation of a polarized luminescence "pump-and-probe" technique is compared with prior ones and with "pump-only" time-resolved luminescence anisotropy decay methods. The technique is experimentally verified on small latex beads with a variety of diameters, common fluorophore labels, and solvent viscosities. Preliminary measurements on a protein (acetylcholine receptor) in the membrane of nondeoxygenated cells in live culture (rat myotubes) show a difference in rotational diffusion between clustered and nonclustered receptors. In most experiments, signal averaging, high laser power, and automated sample translation must be employed to achieve adequate statistical accuracy.     

Reaction rates of oxygen with hemoglobin measured by non-equilibrium facilitated oxygen diffusion through hemoglobin solutions

The purpose of this study was to verify the concept of non-equilibrium facilitated oxygen diffusion. This work succeeds our previous study, where facilitated oxygen diffusion by hemoglobin was measured at conditions of chemical equilibrium, and which yielded diffusion coefficients of hemoglobin and of oxygen. In the present work chemical non-equilibrium was induced using very thin diffusion layers. As a result, facilitation was decreased as predicted by theory. Thus, this work presents the first experimental demonstration of non-equilibrium facilitated oxygen diffusion. In addition, association and dissociation rate parameters of the reaction between oxygen and bovine and human hemoglobin were calculated and the effect of the homotropic and heterotropic interactions on each rate parameter was demonstrated. The results indicate that the homotropic interaction--which leads to increasing oxygen affinity with increasing oxygenation--is predominantly due to an increase in the association rate. The heterotropic interaction--which leads to decreasing oxygen affinity by anionic ligands--appears to be effected in two ways. Cl- increases the dissociation rate. In contrast, 2,3-diphosphoglycerate decreases the association rate.