Three-component Reaction-diffusion System Research Articles

An innovative Vieta–Fibonacci wavelet collocation method for the numerical solution of three-component Brusselator reaction diffusion system of fractional order

An innovative Vieta–Fibonacci wavelet collocation method for the numerical solution of three-component Brusselator reaction diffusion system of fractional order

Arbitrarily weak head-on collision can induce annihilation: the role of hidden instabilities

AbstractIn this paper, we focus on annihilation dynamics for the head-on collision of traveling patterns. A representative and well-known example of annihilation is the one observed for one-dimensional traveling pulses of the FitzHugh–Nagumo equations. In this paper, we present a new and completely different type of annihilation arising in a class of three-component reaction diffusion system. It is even counterintuitive in the sense that the two traveling spots or pulses come together very slowly but do not merge, keeping some separation, and then they start to repel each other for a certain time. Finally, up and down oscillatory instability emerges and grows enough for patterns to become extinct eventually (see Figs. 1, 2, 3). There is a kind of hidden instability embedded in the traveling patterns, which causes the above annihilation dynamics. The hidden instability here turns out to be a codimension 2 singularity consisting of drift and Hopf (DH) instabilities, and there is a parameter regime emanating from the codimension 2 point in which a new type of annihilation is observed. The above scenario can be proved analytically up to the onset of annihilation by reducing it to a finite-dimensional system. Transition from preservation to annihilation is also discussed in this framework.

Matched asymptotic expansion approach to pulse dynamics for a three-component reaction–diffusion system

We study the existence and stability of standing pulse solutions to a singularly perturbed three-component reaction–diffusion system with one activator and two inhibitors. We apply the MAE (matched asymptotic expansion) method to construct solutions and the SLEP (singular limit eigenvalue problem) method to evaluate their stability. This approach is not just an alternative approach to geometric singular perturbation and the associated Evans function, it also yields two advantages: one is extendability to higher dimensional cases, and the other is that it allows us to obtain more precise information about the behaviors of critical eigenvalues. This implies the existence of codimension-two singularities of drift and Hopf bifurcations for the standing pulse solution. Moreover, it is numerically confirmed that stable standing and traveling breathers emerge around the singularity in a physically acceptable region.

Stability and Hopf Bifurcation of a Reaction–Diffusion System with Weak Allee Effect

A delayed three-component reaction–diffusion system with weak Allee effect and Dirichlet boundary condition is considered. The existence and stability of the positive spatially nonhomogeneous steady-state solution are obtained via the implicit function theorem. Moreover, taking delay as the bifurcation parameter, the Hopf bifurcation near the spatially nonhomogeneous steady-state solution is proved to occur at a critical value. Especially, the direction of Hopf bifurcation is forward and the bifurcated periodic solutions are unstable. Finally, the general results are applied to four types of three-species population models with weak Allee effect in growth.

Traveling waves for a three-component reaction-diffusion model of farmers and hunter-gatherers in the Neolithic transition.

The Neolithic transition began the spread of early agriculture throughout Europe through interactions between farmers and hunter-gatherers about 10,000years ago. Archeological evidences indicate that the expanding velocity of farming into a region occupied by hunter-gatherers is roughly constant all over Europe. In the late twentieth century, from the contribution of the radiocarbon dating, it could be found that there are two types of farmers: one is the original farmer and the other is the converted farmer which is genetically hunter-gatherers but learned agriculture from neighbouring farmers. Then this raises the following questions: Which farming populations play a key role in the expansion of farmer populations in Europe? and what is the fate of hunter-gatherers (e.g., become extinct, or live in lower density, or live in agricultural life-style)? We consider a three-component reaction-diffusion system proposed by Aoki, Shida and Shigesada, which describes the interactions among the original farmers, the converted farmers, and the hunter-gatherers. In order to resolve these two questions, we discuss traveling wave solutions which give the information of the expanding velocity of farmer populations. The main result is that two types of traveling wave solutions exist, depending on the growth rate of the original farmer population and the conversion rate of the hunter-gatherer population to the converted farmer population. The profiles of traveling wave solutions indicate that the expansion of farmer populations is determined by the growth rate of the original farmer and the (maximal) carrying capacity of the converted farmer, and the fate of hunter-gatherers is determined by the growth rate of the hunter-gatherer and the conversion rate of the hunter-gatherer to the converted farmer. Thus, our results provide a partial answer to the above two questions.

Taxis-driven Formation of Singular Hotspots in a May--Nowak Type Model for Virus Infection

A three-component reaction-diffusion system is considered which originates from an extension of the classical May--Nowak model for viral infections to situations in which spatially heterogeneous dy...

Spreading properties of a three-component reaction-diffusion model for the population of farmers and hunter-gatherers

In this paper, we investigate the spreading properties of solutions of farmer and hunter-gatherer model which is a three-component reaction-diffusion system. Ecologically, the model describes the geographical spreading of an initially localized population of farmers into a region occupied by hunter-gatherers. This model was proposed by Aoki, Shida and Shigesada in 1996. By numerical simulations and some formal linearization arguments, they concluded that there are four different types of spreading behaviors depending on the parameter values. Despite such intriguing observations, no mathematically rigorous studies have been made to justify their claims. The main difficulty comes from the fact that the comparison principle does not hold for the entire system. In this paper, we give theoretical justification to all of the four types of spreading behaviors observed by Aoki et al. Furthermore, we show that a logarithmic phase drift of the front position occurs as in the scalar KPP equation. We also investigate the case where the motility of the hunter-gatherers is larger than that of the farmers, which is not discussed in the paper of Aoki et al.

Conditional symmetries and exact solutions of a nonlinear three-component reaction-diffusion model

AbstractQ-conditional (non-classical) symmetries of the known three-component reaction-diffusion (RD) system [K. Aoki et al. Theor. Popul. Biol. 50, 1–17 (1996)] modelling interaction between farmers and hunter-gatherers are constructed for the first time. A wide variety of Q-conditional symmetries are found, and it is shown that these symmetries are not equivalent to the Lie symmetries. Some operators of Q-conditional (non-classical) symmetry are applied for finding exact solutions of the RD system in question. Properties of the exact solutions (in particular, their asymptotic behaviour) are identified and possible biological interpretation is discussed.

Exact Solutions of a Mathematical Model Describing Competition and Co-Existence of Different Language Speakers.

The known three-component reaction–diffusion system modeling competition and co-existence of different language speakers is under study. A modification of this system is proposed, which is examined by the Lie symmetry method; furthermore, exact solutions in the form of traveling fronts are constructed and their properties are identified. Plots of the traveling fronts are presented and the relevant interpretation describing the language shift that has occurred in Ukraine during the Soviet times is suggested.

Cops-on-the-dots: The linear stability of crime hotspots for a 1-D reaction-diffusion model of urban crime

In a singularly perturbed limit, we analyse the existence and linear stability of steady-state hotspot solutions for an extension of the 1-D three-component reaction-diffusion (RD) system formulated and studied numerically in Jones et. al. [Math. Models. Meth. Appl. Sci., 20, Suppl., (2010)], which models urban crime with police intervention. In our extended RD model, the field variables are the attractiveness field for burglary, the criminal population density and the police population density. Our model includes a scalar parameter that determines the strength of the police drift towards maxima of the attractiveness field. For a special choice of this parameter, we recover the ‘cops-on-the-dots’ policing strategy of Jones et. al., where the police mimic the drift of the criminals towards maxima of the attractiveness field. For our extended model, the method of matched asymptotic expansions is used to construct 1-D steady-state hotspot patterns as well as to derive nonlocal eigenvalue problems (NLEPs) that characterise the linear stability of these hotspot steady states to ${\cal O}$(1) timescale instabilities. For a cops-on-the-dots policing strategy, we prove that a multi-hotspot steady state is linearly stable to synchronous perturbations of the hotspot amplitudes. Alternatively, for asynchronous perturbations of the hotspot amplitudes, a hybrid analytical–numerical method is used to construct linear stability phase diagrams in the police vs. criminal diffusivity parameter space. In one particular region of these phase diagrams, the hotspot steady states are shown to be unstable to asynchronous oscillatory instabilities in the hotspot amplitudes that arise from a Hopf bifurcation. Within the context of our model, this provides a parameter range where the effect of a cops-on-the-dots policing strategy is to only displace crime temporally between neighbouring spatial regions. Our hybrid approach to study the NLEPs combines rigorous spectral results with a numerical parameterisation of any Hopf bifurcation threshold. For the cops-on-the-dots policing strategy, our linear stability predictions for steady-state hotspot patterns are confirmed from full numerical PDE simulations of the three-component RD system.

Unfolding Symmetric Bogdanov–Takens Bifurcations for Front Dynamics in a Reaction–Diffusion System

This manuscript extends the analysis of a much studied singularly perturbed three-component reaction-diffusion system for front dynamics in the regime where the essential spectrum is close to the origin. We confirm a conjecture from a preceding paper by proving that the triple multiplicity of the zero eigenvalue gives a Jordan chain of length three. Moreover, we simplify the center manifold reduction and computation of the normal form coefficients by using the Evans function for the eigenvalues. Finally, we prove the unfolding of a Bogdanov-Takens bifurcation with symmetry in the model. This leads to stable periodic front motion, including stable traveling breathers, and these results are illustrated by numerical computations.

Spreading Waves in a Farmers and Hunter-Gatherers Model of the Neolithic Transition in Europe.

The Neolithic transition began the spread of early agriculture throughout Europe through interactions between farmers and hunter-gatherers about 10,000 years ago. Archeological evidence produced by radiocarbon dating indicates that the expanding velocity of farming is roughly constant all over Europe. Theoretical understanding of such evidence has been performed from mathematical modeling viewpoint. However, the expanding velocity determined by existing modeling approaches is faster than the observed velocity. For understanding this difference, we propose a three-component reaction-diffusion system which consists of two different types of farmers (sedentary and migratory) and hunter-gatherers from the viewpoint of the influence of farming technology. Our purpose is to study the relation between the expanding velocity of farmers and the farming technology parameter (say, [Formula: see text]). In this paper, we mainly focus on the one-dimensional traveling wave solution with minimal velocity and show that the minimal velocity decreases, as [Formula: see text] increases. This can be compatible with the observed velocity when farming technology is developed. Our results suggest that the reason for the slowdown of the Neolithic transition might be related to the increase in the development of farming technology.

Traveling wave solutions for a bacteria system with density-suppressed motility

<p style='text-indent:20px;'>In 2011, Liu et. al. proposed a three-component reaction-diffusion system to model the spread of bacteria and its signaling molecules (AHL) in an expanding cell population. At high AHL levels the bacteria are immotile, but diffuse with a positive diffusion constant at low distributions of AHL. In 2012, Fu et. al. studied a reduced system without considering nutrition and made heuristic arguments about the existence of traveling wave solutions. In this paper we provide rigorous proofs of the existence of traveling wave solutions for the reduced system under some simple conditions of the model parameters.

On a Variational Problem Arising from the Three-component FitzHugh-Nagumo Type Reaction-Diffusion Systems

We study a variational problem arising from the three-component Fitzhugh-Nagumo type reaction diffusion systems and its shadow systems. In [15], Oshita studied the two-component systems. He revealed that a minimizer of energy corresponding to the problem oscillates under an appropriate condition and also studied its stability. Moreover, he mentioned its energy estimate without a proof. We investigate the behavior of a minimizer corresponding to the three-component problem, its stability and its energy estimate and extend some results of Oshita to the three-component systems and its shadow systems. In particular, we give a necessary and sufficient condition that the minimizer highly oscillates as $ \epsilon \to 0 $. Also, we establish a precise order of the energy estimate of the minimizer as $ \epsilon \to 0 $. In the proof of the energy estimate, we propose a new interpolation inequality.

A hunter-gatherer–farmer population model: Lie symmetries, exact solutions and their interpretation

The Lie symmetry classification of the known three-component reaction–diffusion system modelling the spread of an initially localized population of farmers into a region occupied by hunter-gatherers is derived. The Lie symmetries obtained for reducing the system in question to systems of ordinary differential equations (ODEs) and constructing exact solutions are applied. Several exact solutions of travelling front type are also found, their properties are identified and biological interpretation is discussed.

Spiral wave chimeras in locally coupled oscillator systems.

The recently discovered chimera state involves the coexistence of synchronized and desynchronized states for a group of identical oscillators. In this work, we show the existence of (inwardly) rotating spiral wave chimeras in the three-component reaction-diffusion systems where each element is locally coupled by diffusion. A transition from spiral waves with the smooth core to spiral wave chimeras is found as we change the local dynamics of the system or as we gradually increase the diffusion coefficient of the activator. Our findings on the spiral wave chimera in the reaction-diffusion systems suggest that spiral chimera states may be found in chemical and biological systems that can be modeled by a large population of oscillators indirectly coupled via a diffusive environment.

Butterfly Catastrophe for Fronts in a Three-Component Reaction–Diffusion System

We study the dynamics of front solutions in a three-component reaction–diffusion system via a combination of geometric singular perturbation theory, Evans function analysis, and center manifold reduction. The reduced system exhibits a surprisingly complicated bifurcation structure including a butterfly catastrophe. Our results shed light on numerically observed accelerations and oscillations and pave the way for the analysis of front interactions in a parameter regime where the essential spectrum of a single front approaches the imaginary axis asymptotically.

Time-delayed feedback control of breathing localized structures in a three-component reaction-diffusion system.

The dynamics of a single breathing localized structure in a three-component reaction-diffusion system subjected to time-delayed feedback is investigated. It is shown that variation of the delay time and the feedback strength can lead either to stabilization of the breathing or to delay-induced periodic or quasi-periodic oscillations of the localized structure. A bifurcation analysis of the system in question is provided and an order parameter equation is derived that describes the dynamics of the localized structure in the vicinity of the Andronov-Hopf bifurcation. With the aid of this equation, the boundaries of the stabilization domains as well as the dependence of the oscillation radius on delay parameters can be explicitly derived, providing a robust mechanism to control the behaviour of the breathing localized structure in a straightforward manner.

Singular limit analysis of a reaction-diffusion system with precipitation and dissolution in a porous medium

In this paper we consider a three-component reaction-diffusion system with a fast precipitation and dissolution reaction term. We investigate its singular limit as the reaction rate tends to infinity. The limit problem is described by a combination of a Stefan problem and a linear heat equation. The rate of convergence with respect to the reaction rate is established in a specific case.

A Hopf bifurcation in a three-component reaction-diffusion system with a chemorepulsion

A Hopf bifurcation in a three-component reaction-diffusion system with a chemorepulsion