Evolution of a semilinear parabolic system for migration and selection without dominance

The semilinear parabolic system that describes the evolution of the gene frequencies in the diffusion approximation for migration and selection at a multiallelic locus without dominance is investigated. The population occupies a finite habitat of arbitrary dimensionality and shape (i.e., a bounded, open domain in Rd). The selection coefficients depend on position; the drift and diffusion coefficients may depend on position. The primary focus of this paper is the dependence of the evolution of the gene frequencies on λ, the strength of selection relative to that of migration. It is proved that if migration is sufficiently strong (i.e., λ is sufficiently small) and the migration operator is in divergence form, then the allele with the greatest spatially averaged selection coefficient is ultimately fixed. The stability of each vertex (i.e., an equilibrium with exactly one allele present) is completely specified. The stability of each edge equilibrium (i.e., one with exactly two alleles present) is fully described when either (i) migration is sufficiently weak (i.e., λ is sufficiently large) or (ii) the equilibrium has just appeared as λ increases. The existence of unexpected, complex phenomena is established: even if there are only three alleles and migration is homogeneous and isotropic (corresponding to the Laplacian), (i) as λ increases, arbitrarily many changes of stability of the edge equilibria and corresponding appearance of an internal equilibrium can occur and (ii) the conditions for protection or loss of an allele can both depend nonmonotonically on λ. Neither of these phenomena can occur in the diallelic case.     

Homoclinic bifurcations at the onset of pulse self-replication

We establish a series of properties of symmetric, N-pulse, homoclinic solutions of the reduced Gray-Scott system: u^″=uv², v^″=v−uv², which play a pivotal role in questions concerning the existence and self-replication of pulse solutions of the full Gray-Scott model. Specifically, we establish the existence, and study properties, of solution branches in the (α,β)-plane that represent multi-pulse homoclinic orbits, where α and β are the central values of u(x) and v(x), respectively. We prove bounds for these solution branches, study their behavior as α→∞, and establish a series of geometric properties of these branches which are valid throughout the (α,β)-plane. We also establish qualitative properties of multi-pulse solutions and study how they bifurcate, i.e., how they change along the solution branches.     

Numerical detection of instability regions for delay models with delay-dependent parameters

In this paper we are interested in gaining local stability insights about the interior equilibria of delay models arising in biomathematics. The models share the property that the corresponding characteristic equations involve delay-dependent coefficients. The presence of such dependence requires the use of suitable criteria which usually makes the analytical work harder so that numerical techniques must be used. Most existing methods for studying stability switching of equilibria fail when applied to such a class of delay models. To this aim, an efficient criterion for stability switches was recently introduced in [E. Beretta, Y. Kuang, Geometric stability switch criteria in delay differential systems with delay dependent parameters, SIAM J. Math. Anal. 33 (2002) 1144–1165] and extended [E. Beretta, Y. Tang, Extension of a geometric stability switch criterion, Funkcial Ekvac 46(3) (2003) 337–361]. We describe how to numerically detect the instability regions of positive equilibria by using such a criterion, considering both discrete and distributed delay models.     

Dynamic behaviors of a delay differential equation model of plankton allelopathy

In this paper, we consider a modified delay differential equation model of the growth of n-species of plankton having competitive and allelopathic effects on each other. We first obtain the sufficient conditions which guarantee the permanence of the system. As a corollary, for periodic case, we obtain a set of delay-dependent condition which ensures the existence of at least one positive periodic solution of the system. After that, by means of a suitable Lyapunov functional, sufficient conditions are derived for the global attractivity of the system. For the two-dimensional case, under some suitable assumptions, we prove that one of the components will be driven to extinction while the other will stabilize at a certain solution of a logistic equation. Examples show the feasibility of the main results.     

A comparison between two distinct continuous models in projective cluster theory: The median and the tight-span construction

It is possible to consider two variants of cluster theory: Inaffine cluster theory, one considers collections ofsubsets of a given setX of objects or states, whereas inprojective cluster theory, one considers collections ofsplits (orbipartitions) of that set. In both contexts, it can be desirable to produce acontinuous model, that is, a spaceT encompassing the given setX which represents in a well-specified and more or less parsimonious way all possibleintermediate objects ortransition states compatible with certain restrictions derived from the given collection of subsets or splits. We investigate an interesting and intriguing relationship between two such constructions that appear in the context of projective cluster theory: TheBuneman construction and thetight-span (or justT)construction.     

Boundedness vs. blow-up in a chemotaxis system

We determine the critical blow-up exponent for a Keller-Segel-type chemotaxis model, where the chemotactic sensitivity equals some nonlinear function of the particle density. Assuming some growth conditions for the chemotactic sensitivity function we establish an a priori estimate for the solution of the problem considered and conclude the global existence and boundedness of the solution. Furthermore, we prove the existence of solutions that become unbounded in finite or infinite time in that situation where this a priori estimate fails.     

Proton–neutron structure of the N=52 nucleus Zr

Following the successful identification of mixed-symmetric one- and two-phonon states in the N=52 nuclei ⁹⁴Mo and ⁹⁶Ru, we have performed a photon scattering experiment on the N=52 isotone ⁹²Zr. Experimental data and shell model calculations show that both, single particle and collective degrees of freedom are present in the low-lying levels of ⁹²Zr. The second excited quadrupole state shows the signatures of the one-phonon mixed-symmetric 2⁺ state, while calculations and data indicate an almost pure neutron configuration for the 2⁺₁ state, in contradiction with the F-spin symmetric limit. Furthermore, two strong dipole excitations, which are candidates for the two-phonon quadrupole–octupole coupled E1 excitation and for the mixed-symmetric 1⁺ two-phonon state, were observed.     

On the geometry of wave solutions of a delayed reaction-diffusion equation

The aim of this paper is to study the existence and the geometry of positive bounded wave solutions to a non-local delayed reaction-diffusion equation of the monostable type.     

Stationary probability distributions for multiresponse linear learning models

A necessary and sufficient condition is given for the existence of stationary probability distributions of a non-Markovian model with linear transition rule. Similar to the Markovian case, stationary probability distributions are characterized as eigenvectors of nonnegative matrices. The model studied includes as special cases the Markovian model as well as the linear learning model and has applications in psychological and biological research, in control theory, and in adaption theory.