The existence and stability of spike equilibria in the one-dimensional Gray–Scott model on a finite domain

Results for the existence, stability, and pulse-splitting behavior of spike patterns in the one-dimensional Gray–Scott model on a finite domain in the semi-strong spike interaction regime are summarized. Conditions on the parameters for the existence of competition instabilities, synchronous oscillatory instabilities, or pulse-splitting behavior of spike patterns are given.     

The Existence and Stability of Spike Equilibria in the One-Dimensional Gray–Scott Model: The Low Feed-Rate Regime

In a singularly perturbed limit of small diffusivity ɛ of one of the two chemical species, equilibrium spike solutions to the Gray–Scott (GS) model on a bounded one-dimensional domain are constructed asymptotically using the method of matched asymptotic expansions. The equilibria that are constructed are symmetric k -spike patterns where the spikes have equal heights. Two distinguished limits in terms of a dimensionless parameter in the reaction-diffusion system are considered: the low feed-rate regime and the intermediate regime. In the low feed-rate regime, the solution branches of k -spike equilibria are found to have a saddle-node bifurcation structure. The stability properties of these branches of solutions are analyzed with respect to the large eigenvalues λ in the spectrum of the linearization. These eigenvalues, which have the property that  λ= O (1)  as  ɛ→ 0  , govern the stability of the solution on an O (1) time scale. Precise conditions, in terms of the nondimensional parameters, for the stability of symmetric k -spike equilibrium solutions with respect to this class of eigenvalues are obtained. In the low feed-rate regime, it is shown that a large eigenvalue instability leads either to a competition instability, whereby certain spikes in a sequence are annihilated, or to an oscillatory instability (typically synchronous) of the spike amplitudes as a result of a Hopf bifurcation. In the intermediate regime, it is shown that only oscillatory instabilities are possible, and a scaling-law determining the onset of such instabilities is derived. Detailed numerical simulations are performed to confirm the results of the stability theory. It is also shown that there is an equivalence principle between spectral properties of the GS model in the low feed-rate regime and the Gierer–Meinhardt model of morphogenesis. Finally, our results are compared with previous analytical work on the GS model.     

Zigzag and Breakup Instabilities of Stripes and Rings in the Two-Dimensional Gray–Scott Model

Two different types of instabilities of equilibrium stripe and ring solutions are studied for the singularly perturbed two‐component Gray–Scott (GS) model in a two‐dimensional domain. The analysis is performed in the semi‐strong interaction limit where the ratio O(?^?2) of the two diffusion coefficients is asymptotically large. For ?→ 0 , an equilibrium stripe solution is constructed where the singularly perturbed component concentrates along the mid‐line of a rectangular domain. An equilibrium ring solution occurs when this component concentrates on some circle that lies concentrically within a circular cylindrical domain. For both the stripe and the ring, the spectrum of the linearized problem is studied with respect to transverse (zigzag) and varicose (breakup) instabilities. Zigzag instabilities are associated with eigenvalues that are asymptotically small as ?→ 0 . Breakup instabilities, associated with eigenvalues that are O(1) as ?→ 0 , are shown to lead to the disintegration of a stripe or a ring into spots. For both the stripe and the ring, a combination of asymptotic and numerical methods are used to determine precise instability bands of wavenumbers for both types of instabilities. The instability bands depend on the relative magnitude, with respect to ?, of a nondimensional feed‐rate parameter A of the GS model. Both the high feed‐rate regime A=O(1) , where self‐replication phenomena occurs, and the intermediate regime O(?^1/2) ?A?O(1) are studied. In both regimes, it is shown that the instability bands for zigzag and breakup instabilities overlap, but that a zigzag instability is always accompanied by a breakup instability. The stability results are confirmed by full numerical simulations. Finally, in the weak interaction regime, where both components of the GS model are singularly perturbed, it is shown from a numerical computation of an eigenvalue problem that there is a parameter set where a zigzag instability can occur with no breakup instability. From full‐scale numerical computations of the GS, it is shown that this instability leads to a large‐scale labyrinthine pattern.     

Numerical Challenges for Resolving Spike Dynamics for Two One-Dimensional Reaction-Diffusion Systems

Asymptotic and numerical methods are used to highlight different types of dynamical behaviors that occur for the motion of a localized spike-type solution to the singularly perturbed Gierer–Meinhardt and Schnakenberg reaction-diffusion models in a one-dimensional spatial domain. Depending on the parameter range in these models, there can either be a slow evolution of a spike toward the midpoint of the domain, a sudden oscillatory instability triggered by a Hopf bifurcation leading to an intricate temporal oscillation in the height of the spike, or a pulse-splitting instability leading to the creation of new spikes in the domain. Criteria for the onset of these oscillatory and pulse-splitting instabilities are obtained through asymptotic and numerical techniques. A moving-mesh numerical method is introduced to compute these different behaviors numerically, and results are compared with corresponding results computed using a method of lines based software package.     

Weak solutions and supersolutions in $ L^1 $ for reaction-diffusion systems

We prove here that limits of nonnegative solutions to reaction-diffusion systems whose nonlinearities are bounded in always converge to supersolutions of the system. The motivation comes from the question of global existence in time of solutions for the wide class of systems preserving positivity and for which the total mass of the solution is uniformly bounded. We prove that, for a large subclass of these systems, weak solutions exist globally. RID="h1" ID="h1"A Philippe, mon ma?tre et ami     

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.     

基于C-STGARCH模型的大庆原油现货价格波动研究

引用Dueker等(2011)提出的同期门槛平滑转换广义自回归条件异方差(C-STGARCH)模型对我国大庆原油现货价格的波动状态进行了实证分析,以求对大庆原油现货价格的波动有一个新的、更深刻的量度.研究显示:第一,大庆原油现货价格的波动是不稳定的,并且存在显著的非对称和非线性现象;第二,CSTGARCH模型能很好地刻画大庆原油现货价格波动的这些现象,并且发现油价的波动以3.738%为门槛点存在高波动区和低波动区两种状态,低波动区的波动持续性比高波动区强,然而,对平滑转换持续性的影响方面,高波动区要略大于低波动区.     

Asymptotic behavior in chemical reaction-diffusion systems with boundary equilibria

We consider the asymptotic behavior for large time of solutions to reaction-diffusion systems modeling reversible chemical reactions. We focus on the case where multiple equilibria exist. In this case, due to the existence of so-called "boundary equilibria", the analysis of the asymptotic behavior is not obvious. The solution is understood in a weak sense as a limit of adequate approximate solutions. We prove that this solution converges in L^1 toward an equilibrium as time goes to infinity and that the convergence is exponential if the limit is strictly positive.     

The effect of strongly anisotropic turbulent mixing on critical behavior: Renormalization group analysis of two nonstandard systems

We consider the effect of strongly anisotropic turbulent mixing on the critical behavior of two systems: a φ ³ critical dynamics model describing universal properties of metastable states in the vicinity of a firstorder phase transition and a reaction-diffusion system near the point of a second-order transition between fluctuation and absorption states (a simple epidemic process or the Gribov process). In both cases, we demonstrate the existence of a new strongly nonequilibrium, anisotropic scaling regime (universality class) for which both the mixing and the nonlinearity in the order parameter are relevant. We evaluate the corresponding critical dimensions in the one-loop approximation of the renormalization group.     

Second Order Semiclassics with Self-Generated Magnetic Fields

We consider the semiclassical asymptotics of the sum of negative eigenvalues of the three-dimensional Pauli operator with an external potential and a self-generated magnetic field B. We also add the field energy bòB²{\beta \int B^{2}} and we minimize over all magnetic fields. The parameter β effectively determines the strength of the field. We consider the weak field regime with β h ² ≥ const > 0, where h is the semiclassical parameter. For smooth potentials we prove that the semiclassical asymptotics of the total energy is given by the non-magnetic Weyl term to leading order with an error bound that is smaller by a factor h^1+e{h^{1+\varepsilon}} , i.e. the subleading term vanishes. However for potentials with a Coulomb singularity, the subleading term does not vanish due to the non-semiclassical effect of the singularity. Combined with a multiscale technique, this refined estimate is used in the companion paper (Erdős et al. in Scott correction for large molecules with a self-generated magnetic field, Preprint, 2011) to prove the second order Scott correction to the ground state energy of large atoms and molecules.     

Hopf bifurcations in a reaction-diffusion population model with delay effect

A reaction-diffusion population model with a general time-delayed growth rate per capita is considered. The growth rate per capita can be logistic or weak Allee effect type. From a careful analysis of the characteristic equation, the stability of the positive steady state solution and the existence of forward Hopf bifurcation from the positive steady state solution are obtained via the implicit function theorem, where the time delay is used as the bifurcation parameter. The general results are applied to a “food-limited” population model with diffusion and delay effects as well as a weak Allee effect population model.     

Saddle-point behavior for monotone semiflows and reaction-diffusion models

The saddle-point behavior is established for monotone semiflows with weak bistability structure and then these results are applied to three reaction-diffusion systems modeling man-environment-man epidemics, single-loop positive feedback and two species competition, respectively.     

Numerical simulation of reaction-diffusion systems by modified cubic B-spline differential quadrature method

In this paper, we have applied modified cubic B-spline based differential quadrature method to get numerical solutions of one dimensional reaction-diffusion systems such as linear reaction-diffusion system, Brusselator system, Isothermal system and Gray-Scott system. The models represented by these systems have important applications in different areas of science and engineering. The most striking and interesting part of the work is the solution patterns obtained for Gray Scott model, reminiscent of which are often seen in nature. We have used cubic B-spline functions for space discretization to get a system of ordinary differential equations. This system of ODE’s is solved by highly stable SSP-RK43 method to get solution at the knots. The computed results are very accurate and shown to be better than those available in the literature. Method is easy and simple to apply and gives solutions with less computational efforts.     

Asymmetric Spotty Patterns for the Gray–Scott Model in R2

In this paper, we rigorously prove the existence and stability of asymmetric spotty patterns for the Gray–Scott model in a bounded two-dimensional domain. We show that given any two positive integers   k ₁, k ₂  , there are asymmetric solutions with   k ₁  large spots (type A) and   k ₂  small spots (type B). We also give conditions for their location and calculate their heights. Most of these asymmetric solutions are shown to be unstable. However, in a narrow range of parameters, asymmetric solutions may be stable.     

Explicit constructions and properties of generalized shift-invariant systems in $L^{2}(\mathbb {R})$

Generalized shift-invariant (GSI) systems, originally introduced by Hernández et al. and Ron and Shen, provide a common frame work for analysis of Gabor systems, wavelet systems, wave packet systems, and other types of structured function systems. In this paper we analyze three important aspects of such systems. First, in contrast to the known cases of Gabor frames and wavelet frames, we show that for a GSI system forming a frame, the Calderón sum is not necessarily bounded by the lower frame bound. We identify a technical condition implying that the Calderón sum is bounded by the lower frame bound and show that under a weak assumption the condition is equivalent with the local integrability condition introduced by Hernández et al. Second, we provide explicit and general constructions of frames and dual pairs of frames having the GSI-structure. In particular, the setup applies to wave packet systems and in contrast to the constructions in the literature, these constructions are not based on characteristic functions in the Fourier domain. Third, our results provide insight into the local integrability condition (LIC).     

Smoothing and Rothe's method for Stefan problems in enthalpy form

The classical two-phase Stefan problem as well as its weak variational formulation model the connection between the different phases of the considered material by interface conditions at the occurring free boundary or by a jump of the enthalpy. One way to treat the corresponding discontinuous variational problems consists in its embedding into a family of continuous ones and applying some standard techniques to the chosen approximation problems. The aim of the present paper is to analyze a semi-discretization via Rothe's method and its convergence behavior in dependence of the smoothing parameter. While in Grossmann et al. (Optimization, in preparation) the treatment of the Stefan problem is based on the given variable, i.e. the temperature, here first a transformation via the smoothed enthalpy is applied. Numerical experiments indicate a higher stability of the discretization by Rothe's method. In addition, to avoid inner iterations a frozen coefficient approach as common in literature is used.     

Stability of planar shock fronts for multidimensional systems of relaxation equations

We investigate stability of multidimensional planar shock profiles of a general hyperbolic relaxation system whose equilibrium model is a system, under the necessary assumption of spectral stability and a standard set of structural conditions that are known to hold for many physical systems. Our main result, generalizing the work of Kwon and Zumbrun in the scalar relaxation case, is to establish the bounds on the Green?s function for the linearized equation and obtain nonlinear L² asymptotic behavior/sharp decay rate of perturbed weak shock profiles. To establish Green?s function bounds, we use the semigroup approach in the low-frequency regime, and use the energy method for the high-frequency bounds, separately. For the system equilibrium case, the analysis of the linearized equation is complicated due to glancing phenomena. We treat this difficulty similarly as in the inviscid and viscous systems, under the constant multiplicity condition.     

Global existence and regularity of solutions to a system of nonlinear Maxwell equations

We consider the model that has been suggested by Greenberg et al. (Physica D 134 (1999) 362-383) for the ferroelectric behavior of materials. In this model, the usual (linear) Maxwell's equations are supplemented with a constitutive relation in which the electric displacement equals a constant times the electric field plus an internal polarization variable which evolves according to an internal set of nonlinear Maxwell's equations. For such model we provide rigorous proofs of global existence, uniqueness, and regularity of solutions. We also provide some preliminary results on the long-time behavior of solutions. The main difficulties in this study are due to the loss of compactness in the system of Maxwell's equations. These results generalize those of Greenberg et al., where only solutions with TM (transverse magnetic) symmetry were considered.     

Discrete infinite-dimensional type-K monotone dynamical systems and time-periodic reaction-diffusion systems

The asymptotic behavior of discrete type-K monotone dynamical systems and reaction-diffusion equations is investigated. The studying content includes the index theory for fixed points, permanence, global stability, convergence everywhere and coexistence. It is shown that the system has a globally asymptotically stable fixed point if every fixed point is locally asymptotically stable with respect to the face it belongs to and at this point the principal eigenvalue of the diagonal partial derivative about any component not belonging to the face is not one. A nice result presented is the sufficient and necessary conditions for the system to have a globally asymptotically stable positive fixed point. It can be used to establish the sufficient conditions for the system to persist uniformly and the convergent result for all orbits. Applications are made to time-periodic Lotka-Volterra systems with diffusion, and sufficient conditions for such systems to have a unique positive periodic solution attracting all positive initial value functions are given. For more general time-periodic type-K monotone reaction-diffusion systems with spatial homogeneity, a simple condition is given to guarantee the convergence of all positive solutions.     

A remark on triviality for the two-dimensional stochastic nonlinear wave equation

We consider the two-dimensional stochastic damped nonlinear wave equation (SdNLW) with the cubic nonlinearity, forced by a space-time white noise. In particular, we investigate the limiting behavior of solutions to SdNLW with regularized noises and establish triviality results in the spirit of the work by Hairer et al. (2012). More precisely, without renormalization of the nonlinearity, we establish the following two limiting behaviors; (i) in the strong noise regime, we show that solutions to SdNLW with regularized noises tend to 0 as the regularization is removed and (ii) in the weak noise regime, we show that solutions to SdNLW with regularized noises converge to a solution to a deterministic damped nonlinear wave equation with an additional mass term.