Short note on a solution with large amplitude for the limiting system arising from the competition-diffusion system

This paper is concerned with positive solutions of the limiting system arising from the Shigesada–Kawasaki–Teramoto model with large interspecific competition rate. It has previously been suggested that the limiting system has positive solutions with large amplitude for a certain value of parameters. As a first step, the purpose of this paper is to discuss the precise spatial profile of such solutions by employing formal calculations by the singular perturbation technique.


Introduction
In 1979, Shigesada et al. [7] proposed the competition-diffusion system with nonlinear diffusion effect to model the segregation of interacting species.The following is the stationary problem of their proposed system: where the variables w and z mean the population density of two competing species, the parameters ε, d, θ , and γ are positive, the cross-diffusion coefficients α and β are nonnegative, the habitat of two competing species is a bounded domain in R with smooth boundary ∂ , the integer is the space dimension of , the vector ν = ν(x) is the outward unit normal vector on x ∈ ∂ , the notation ∂ ξ f denotes the derivative of the function f with respect to ξ .We remark that the values θ and γ θ correspond to the interspecific competition rate.Since w(x) and z(x) mean the population density, we restrict our discussion to the positive solution of problem (1), where we say that (w, z)(x) is positive if (w, z)(x) is in the positive quadrant for any x in the closure Cl of .Since then, many mathematicians have researched the solution structure of positive solution for problem (1) to understand what kinds of phenomena occur by the cross-diffusion effect (for example, see Jüngel [3] and Yamada [8,9] as surveys of mathematical study related to problem (1)).
Since the change of variables p = (ψ p , ψ q )(w) with w = (w, z), p = (p, q), p = ψ p (w) = (1 + αz)w, q = ψ q (w) = (1 + βw)z is a bijective map from the positive quadrant to itself, it turns out that there exists an inverse map w = ( w , z )(p) of p = (ψ u , ψ v )(w) in the positive quadrant.After simple calculations, we obtain for any p in the positive quadrant, and then we represent problem (1) as which is a stationary problem of reaction-diffusion system with linear diffusion term, where We remark that hold true for any p in the positive quadrant.
In the previous paper [5], it was shown that problem (2) as θ → +∞ has two kinds of limiting systems.One limiting system is represented as by applying the argument in Dancer et al. [2], where = -∪ + ∪ , is the interface which separates two subregions and μ = μ(x) is a unit vector normal to .We note that the solution of limiting system (3) approximates the bounded solution of problem (1) for sufficiently large θ > 0. By employing the change of variables u = γ θp and v = θ q, and taking the limit as θ → +∞, we obtain the other limiting system of problem (2).Since the positive constant solution of problem (1) in the positive quadrant is the corresponding constant solution of problem (4) becomes Hence, it turns out that the solution of limiting system (4) approximates the solution with small amplitude of problem (1) in a neighborhood of w = ŵ for sufficiently large θ > 0. We note that the solution structure of limiting system (3) has been already investigated (for example, see Dancer et al. [2]).However, although limiting system (4) is simple in a class of 2-component reaction-diffusion systems, we have not yet understood enough solution structure of positive solution for limiting system (4).As a first step to approach the solution structure of problem (4), we assume that the habitat is the ball with center origin and radius π , and we restrict our discussion to the radially symmetric positive solution, that is, we study positive solutions of the problem where In the previous paper [5], it was suggested that problem (5) has a solution with very large amplitude for some ε > 0. In this paper, to discuss precise spatial profile of such a solution of problem (5), we construct the solution for some ε > 0.
To do this, we employ formal calculations by the singular perturbation technique because the rigorous proof on the existence theorem of solution requires lengthy arguments (for example, see Ito [4]).

Preliminary
Let (u, v)(r) be an arbitrary bounded positive solution of problem (5) for ε > 0 and d > 0.
When u(r) > 1 (respectively, v(r) > 1) is satisfied for any r ∈ [0, π], we have because of the boundary condition.This contradiction implies that the following holds true.
Lemma 1 Let (u, v)(r) be an arbitrary bounded positive solution of problem (5) for ε > 0 and d > 0. Then there exist r u and r v ∈ [0, π] such that u(r u ) ≤ 1 and v(r v ) ≤ 1 are satisfied. Setting we consider solutions of the problem with We denote by u + (r) (respectively, u -(r)) the solution of the problem 0 = L D(u) + u, r > 0 with u(0) = 1 (respectively, u(0) = -1) and ∂ r u(0) = 0.The phase plane analysis leads to the following.
For each k ∈ N and σ ∈ {-, +}, setting we find out that ūk (r) and ū+ k (r) are solutions of problem ( 6) for ε = ε - k and ε = ε + k , respectively, and that u - k (r) and u + k (r) are quasi-positive solutions of the problem where (0, π) = I -∪ I 0 ∪ I + , I 0 is the interface that separates two subintervals and we say that We remark that problem ( 7) is a linearized one for problem (5) around (u, v) = 0.Moreover, in the previous paper [5], it was suggested that positive solutions with large amplitude for problem (5) appear in a neighborhood of ε = ε - 1 and ε = ε + 1 .

Solutions with large amplitude
Setting we represent problem (5) as follows: We remark that for each k ∈ N, u - k (r) and u + k (r) are quasi-positive solutions of problem ( 8) We only discuss the construction of positive solution for problem (8) in a neighborhood of ε = ε + 1 by formal calculations employing the singular perturbation technique, because we can similarly discuss that in a neighborhood of ε = ε - 1 .We set Since our purpose of this paper is to seek classical solutions of problem (8), we may impose on the condition at r = τ .In an analogous way to the argument in Ito [4], we introduce δ as an auxiliary parameter and assume the following estimate: as δ → 0, where each function is bounded and continuous.Setting for each k ∈ N and γ 0 = v + 0,1 , we have as δ → 0 for any compact subset of R. The definition of u 0 (r) and v 0 (r) implies that hold true.Substituting the above estimate (10) into problem (8), we have as δ → 0 for any compact subset of R ± , where R -= (-∞, 0) and R + = (0, +∞).Moreover, we see from condition (9) that are satisfied.
Employing v 1 (±∞) = 0, we calculate the asymptotic expansion of φ -(ξ ) and φ + (ξ ) as |ξ | → +∞, and then we have where Bi(t) is the Airy function of the second kind, and C - ± and C + ± are suitable nonzero constants.Employing the asymptotic expansion of Ai(t) and Bi(t) as t → +∞, we have as t → +∞, and then we obtain Let f (ξ ) be an arbitrary bounded function from R to itself.After simple calculations, we see that the bounded solution W (ξ ) of the problem for any ξ ∈ R, where We denote by K -1 f (ξ ) the solution W (ξ ) given in equation (15).Since holds true because of φ -(ξ ) = φ + (-ξ ) for any ξ ∈ R, we obtain W (-ξ ) = W (ξ ) (respectively, W (-ξ ) = -W (ξ ))for any ξ ∈ R if f (ξ ) is an even (respectively, odd) function in ξ ∈ R Let z(ξ ) be a solution of problem (13) with z(±∞) = 0 for (ε 1 , q 2 , q2 ) ∈ R 3 .Since A(ξ ) is an even function in ξ ∈ R, it follows that Since F -(ξ ) and F + (ξ ) are bounded in ξ ∈ R, we can represent the solution z(ξ ) of problem (13) as where By integration by parts, we have which implies that q2 and Ẑ+ (0) are determined by the above equations.Hence, we see from q 2 = Ẑ+ (0)ε 1 that the positive solution of problem (12) is determined for any fixed ε 1 ∈ R. Because rigorous proofs are not given in the above arguments, the summary of this paper is written down as a conjecture.

Conclusion and future works
Since we can study the rough solution structure and the fine solution structure of problem (1) through limiting systems (3) and ( 4), respectively, we might guess that these structures lead to the global solution structure for problem (1).In this paper, we formally constructed a family of solutions for problem (8) arising from limiting system (4).To do this, we employed the singular perturbation technique such as in Ito [4].
Recently in [6] it was shown that 0 < ε - 1 < ε + 1 holds true for any d > 1.Furthermore, employing the argument as in Dancer et al. [2], it follows that limiting system (3) has no nontrivial solutions for any ε > ε + 1 and has quasi-positive solutions for any 0 < ε < ε + 1 .Although we have not yet determined the value of ε 1 in estimate (10), it is possible that we could determine its value by studying the higher order asymptotic expansion of the solution as δ 0. If the validity of ε(δ) > ε + 1 can be proved, then we obtain a positive solution of problem (1), which does not appear in limiting system (3), by employing limiting system (4).
Since problem (1) has nonlinear diffusion terms, only the linearized stability is not sufficient for the stability analysis of solution for the system of evolution equations whose stationary problem is represented as problem (1), so that the stability analysis is still open also for the solution corresponding to the positive solution of limiting systems (3) and (4).Since in this paper we have obtained the asymptotic expansion of positive solutions with large amplitude for limiting system (4), it may be possible to establish the stability of the solution by using the expansion of the solution constructed in this paper.
Since the existence of the solutions constructed in this paper seems to be proved in an analogous way to the argument in Ito [4], one of our future works is that we establish a rigorous proof.