Skip to main content

Existence, uniqueness and stability of positive solutions to a general sublinear elliptic systems

Abstract

In this paper, we make use of a new stability result and bifurcation theory to study the existence and uniqueness of positive solutions to semilinear elliptic systems with some general sublinear conditions. Moreover, we obtain the precise global bifurcation diagrams of the system in a single monotone solution curve.

MSC:35J55, 35B32.

1 Introduction

We consider positive solutions of a semilinear elliptic system with n equations (n2)

{ Δ u 1 + λ f 1 ( u 1 , u 2 , , u n ) = 0 , x Ω , Δ u 2 + λ f 2 ( u 1 , u 2 , , u n ) = 0 , x Ω , Δ u n + λ f n ( u 1 , u 2 , , u n ) = 0 , x Ω , u 1 ( x ) = = u n ( x ) = 0 , x Ω ,
(1.1)

where λ>0 is a positive parameter, Ω is a bounded smooth domain in R N , and f i (i=1,2,,n) are smooth real-valued functions defined on R + n ={( u 1 , u 2 ,, u n ) R n : u i >0,1in} satisfying the following.

Cooperativeness Define the Jacobian of the vector field ( f 1 ,, f n ) as

H( u 1 ,, u n )= ( f 1 u 1 f 1 u n f n u 1 f n u n ) = ( f 11 f 1 n f n 1 f n n ) .
(1.2)

Then f i / u j 0 (ij) for ( u 1 ,, u n ) R + n .

The purpose of this paper is to study the existence, uniqueness and stability of positive solutions of such cooperative system (1.1) under certain conditions of f i .

The existence, uniqueness and stability of positive solutions to sublinear semilinear elliptic systems with two equations have been recently studied in [14]. The sublinear condition plays an important role. In this paper, we continue the effort in [3] to prove the stability of a positive solution to (1.1) under some reasonable sublinear conditions, and the stability implies the uniqueness of the positive solution. We also prove corresponding existence results using bifurcation theory and the continuation method. This is motivated by the existence study of exact multiplicity (and uniqueness) of positive solutions to the scalar semilinear elliptic equation:

{ Δ u + λ f ( u ) = 0 , x Ω , u ( x ) = 0 , x Ω ,
(1.3)

starting from Korman et al. [5, 6]. In Ouyang and Shi [7, 8], their result classified the different global exact multiplicity of (1.3) for more general nonlinearity f. There are more results on the existence and uniqueness of solution to the semilinear cyclic elliptic system

{ Δ u + λ f ( v ) = 0 , x Ω , Δ v + λ g ( u ) = 0 , x Ω , u ( x ) = v ( x ) = 0 , x Ω .
(1.4)

The notation of sublinearity and superlinearity of the nonlinear vector field or the ones in higher dimension was considered in Sirakov [9]. Our definition of sublinear nonlinearity is quite different, and ours is similar to the one in Ouyang and Shi [8] for the scalar case. Dalmasso [10] obtained the existence and uniqueness result for a more special sublinear system, and it was extended by Shi and Shivaji [4]. The uniqueness of a positive solution for large λ was proved in Hai [11, 12], Hai and Shivaji [13]. If Ω is a finite ball or the whole space, then the positive solutions of systems are radically symmetric and decreasing in radical direction by the result of Troy [14]. Hence the system can be converted into a system of ODEs. Several authors have taken that approach for the existence of the solutions, see Serrin and Zou [15, 16], and much success has been achieved for Lane-Emden systems. Using the scaling invariant in the Lane-Emden system, the uniqueness of the radial positive solution was shown in Dalmasso [10], Korman and Shi [17]. Cui et al. [18, 19] considered cyclic systems with three equations, and the uniqueness and existence of positive solutions were obtained. For the Lane-Emden systems with n equations, Maniwa [20] obtained the uniqueness and existence of positive solutions to systems under the sublinear conditions.

We organize the rest of this paper in the following way. In Section 2, we recall the maximum principle and prove the main stability result. In Section 3, we use the stability result and bifurcation theory to prove the existence and uniqueness of a positive solution. We also obtain the precise global bifurcation diagrams of the system (the bifurcation diagram is a single monotone solution curve in all cases) and give some examples. In Section 4, we consider the similar question for merely Hölder continuous nonlinearities, and we use monotone methods for existence. We use W 2 , p (Ω) and W loc 2 , p (Ω) for the standard Sobolev space, C( Ω ¯ ) for the space of continuous functions defined on Ω ¯ , and C 0 ( Ω ¯ )={uC( Ω ¯ ):u(x)=0,xΩ}. We use N(L) and R(L) to denote the null space and the range space of a linear operator L.

2 Stability and linearized equations

In this section, we study the stability result about a positive solution. Let U=( u 1 ,, u n ) be a solution of (1.1). We shall denote the partial derivative of f i (U) with respect to u j by f i j (U) or f i j . The stability of U is determined by the eigenvalue equation

{ Δ ξ 1 + λ f 11 ξ 1 + λ f 12 ξ 2 + + λ f 1 n ξ n = μ ξ 1 , x Ω , Δ ξ 2 + λ f 21 ξ 1 + λ f 22 ξ 2 + + λ f 2 n ξ n = μ ξ 2 , x Ω , Δ ξ n + λ f n 1 ξ 1 + λ f n 2 ξ 2 + + λ f n n ξ n = μ ξ n , x Ω , ξ 1 ( x ) = ξ 2 ( x ) = = ξ n ( x ) = 0 , x Ω ,
(2.1)

which can be written as

Lu=Hu+μu,
(2.2)

where

u= ( ξ 1 ξ n ) ,Lu= ( Δ ξ 1 Δ ξ n ) ,andH= ( f 11 f 1 n f n 1 f n n ) .
(2.3)

Definition 2.1 [16]

An n×n matrix A is reducible if for some permutation matrix Q,

QA Q T = ( B 0 C D ) ,

where B and C are square matrices and Q T is the transpose of Q. Otherwise, A is irreducible.

Throughout this paper, H is assumed to be irreducible, since if not the case, the linearized system (2.1) can be reduced to two subsystems with one being not coupled with the other. If we assume that H is cooperative and irreducible, then the maximum principles hold for such systems. Before stating our results, we recall some known results as required.

Lemma 2.2 Let X= [ W loc 2 , p ( Ω ) C 0 ( Ω ¯ ) ] n , and let Y= [ L p ( Ω ) ] n for p>n. Suppose that L, H are given as in (2.3), and H=( f i j ) is irreducible and satisfies f i j 0 (ij) for ( u 1 ,, u n ) R + n . Then we have the following.

  1. 1.

    μ 1 =inf{μspt(LH)} is a real eigenvalue of LH, where spt(LH) is the spectrum of LH.

  2. 2.

    For μ= μ 1 , there exists a unique (up a constant multiple) eigenfunction u 1 X, and u 1 >0 in Ω.

  3. 3.

    For μ< μ 1 , the equation Lu=Hu+μu+f is uniquely solvable for any fY, and u>0 as long as f0.

  4. 4.

    (Maximum principle) For μ μ 1 , suppose that u [ W loc 2 , p ( Ω ) C ( Ω ¯ ) ] n satisfies LuHu+μu in Ω, u0 on Ω, then u0 in Ω.

  5. 5.

    If there exists u [ W loc 2 , p ( Ω ) C ( Ω ¯ ) ] n satisfying LuHu in Ω, u0 on Ω, and either u0 on Ω or LuHu in Ω, then μ 1 >0.

For the result and proofs, see Sweers [21], Proposition 3.1 and Theorem 1.1. Moreover, from a standard compactness argument, there are countably many eigenvalues { μ i } of LH, and | μ i μ 1 | as i. We notice that μ i is not necessarily real-valued. We call a solution U stable if μ 1 >0; and otherwise, it is unstable ( μ 1 0).

For our purpose, in this section, we also need to consider the adjoint operator of LH. Let the transpose matrix of the Jacobian be

H T = ( f 11 f n 1 f 1 n f n n ) .
(2.4)

Then, evidently, the results in Lemma 2.2 also hold for the eigenvalue problem

L u = H T u +μ u ,
(2.5)

which is

{ Δ ξ 1 + λ f 11 ξ 1 + λ f 21 ξ 2 + + λ f n 1 ξ n = μ ξ 1 , x Ω , Δ ξ 2 + λ f 12 ξ 1 + λ f 22 ξ 2 + + λ f n 2 ξ n = μ ξ 2 , x Ω , Δ ξ n + λ f 1 n ξ 1 + λ f 2 n ξ 2 + + λ f n n ξ n = μ ξ n , x Ω , ξ 1 ( x ) = ξ 2 ( x ) = = ξ n ( x ) = 0 , x Ω ,
(2.6)

where u = ( ξ 1 , ξ 2 , , ξ n ) T . It is easy to verify that L H T is the adjoint operator of LH, while both are considered as operators defined on subspaces of [ L 2 ( Ω ) ] n . By using the well-known functional analytic techniques (see [21, 22]), one can show the following.

Lemma 2.3 Let X, Y, L and H be the same as in Lemma  2.2. Then the principal eigenvalue μ 1 of LH is also a real eigenvalue of L H T , μ 1 =inf{μspt(L H T )}, and for μ= μ 1 , there exists a unique eigenfunction u 1 [ W loc 2 , p ( Ω ) C 0 ( Ω ¯ ) ] n of L H T (up a constant multiple), and u 1 >0 in Ω.

Cui et al. [3] obtained the stability result of a positive solution for the system with two equations. We give the following stability result about a positive solution of (1.1).

Theorem 2.4 Suppose that U=( u 1 ,, u n ) is a positive solution of (1.1), f i is cooperative and H=( f i j ) is irreducible, then U is stable if f i satisfies one of the following conditions: for any U R + n ,

( H 1 )

f i (U)> j = 1 n f i j (U) u j ,i=1,,n;or
(2.7)

( H 2 )

f i (U)> j = 1 n f j i (U) u j ,i=1,,n.
(2.8)

Proof Let U=( u 1 ,, u n ) be a positive solution of (1.1), and let ( μ 1 , u ) be the corresponding principal eigen-pair of (2.6) such that ξ i >0 in Ω for any 1in.

We assume that f i satisfies ( H 1 ). Multiplying the system (1.1) by u , the system (2.6) by U, integrating over Ω and subtracting, we obtain

μ 1 Ω i = 1 n u i ξ i d x = λ Ω i = 1 n f i ξ i d x λ Ω i = 1 n ξ i j = 1 n f i j u j d x = λ Ω i = 1 n ξ i ( f i j = 1 n f i j u j ) d x .
(2.9)

Hence μ 1 >0 if ( H 1 ) is satisfied.

Similar to the proof above, let ( μ 1 ,u) be the corresponding principal eigen-pair of (2.1) such that ξ i >0 in Ω for any i=1,,n. We assume that f i satisfies ( H 2 ). Multiplying the system (1.1) by u, the system (2.1) by U, integrating over Ω and subtracting, we can get

μ 1 Ω i = 1 n u i ξ i d x = λ Ω i = 1 n f i ξ i d x λ Ω i = 1 n ξ i j = 1 n f j i u j d x = λ Ω i = 1 n ξ i ( f i j = 1 n f j i u j ) d x .
(2.10)

Hence μ 1 >0 if ( H 2 ) is satisfied. □

On the other hand, the same proof also implies the following instability result under the opposite condition of ( H 1 ) and ( H 2 ).

Theorem 2.5 Suppose that U=( u 1 ,, u n ) is a positive solution of (1.1), f i is cooperative and H T =( f j i ) is irreducible, then U is unstable if f i satisfies one of the following conditions: for any U R + n ,

( H 1 )

f i (U)< j = 1 n f i j (U) u j ,i=1,,n;or
(2.11)

( H 2 )

f i (U)< j = 1 n f j i (U) u j ,i=1,,n.
(2.12)

Remark 2.6

  1. 1.

    In [8], for positive solutions of the scalar equation

    Δu+λh(u)=0,xΩ,u(x)=0,xΩ,

    the function h(u) is called a sublinear function if h(u)>u h (u), and it is superlinear if h(u)<u h (u). It was proved in Proposition 3.14 of [8] that a positive solution u is stable if h is sublinear, and u is unstable if h is superlinear. Now our conclusions, Theorems 2.4 and 2.5, are generalizations of the corresponding results. The condition ( H i ) for i=1,2 is the generalization of sublinearity (or superlinearity) to n-variable vector fields.

  2. 2.

    The conditions ( H 1 ) and ( H 2 ) can be written in a vector form F(U)>H(U)U and F(U)> H T (U)U, where F(U)= ( f 1 ( U ) , , f n ( U ) ) T , U= ( u 1 , , u n ) T , and H is the original Jacobian matrix of the vector field ( f 1 (U),, f n (U)), and H T is the transpose matrix of the Jacobian H. The condition ( H 1 ) is clearly more natural as the conditions for f i are separate. Hence the sublinearity can be defined for a single n-variable function. The condition ( H 2 ) is defined for the whole vector field ( f 1 (U),, f n (U)).

  3. 3.

    If a solution U is stable, then it is necessarily a non-degenerate solution. That is, any eigenvalue μ i of (2.1) has a positive real part. But when a solution is proved to be unstable, it can be a degenerate one with zero or pure imaginary eigenvalues.

3 Existence and uniqueness

In this section, we consider the uniqueness and existence of positive solutions for the following problem:

{ Δ u 1 + λ ( g 11 ( u 1 ) + g 12 ( u 2 ) + + g 1 n ( u n ) ) = 0 , x Ω , Δ u 2 + λ ( g 21 ( u 1 ) + g 22 ( u 2 ) + + g 2 n ( u n ) ) = 0 , x Ω , Δ u n + λ ( g n 1 ( u 1 ) + g n 2 ( u 2 ) + + g n n ( u n ) ) = 0 , x Ω , u 1 ( x ) = u 2 ( x ) = = u n ( x ) = 0 , x Ω .
(3.1)

Suppose that each g i j (x) (i,j=1,2,,n) is a smooth real function defined on R + satisfying

( A 1 ) g i j (0)0;

( A 2 ) g i j (x)0, ( g i j ( x ) / x ) 0, for all x0.

The Perron-Frobenius theorem plays a critical role in our main result.

Lemma 3.1 (Perron-Frobenius theorem: strong form [[23], Theorem 5.3.1])

Let n×n matrix A be a nonnegative irreducible matrix. Then ρ(A) is a simple eigenvalue of A, associated to a positive eigenvector, where ρ(A) denotes the spectral radius of A. Moreover, ρ(A)>0.

Here let ( λ 1 , φ 1 ) be the principal eigen-pair of

Δφ=λφ,xΩ,φ(x)=0,xΩ,
(3.2)

such that φ 1 (x)>0 in Ω and φ 1 =1. We have the following existence and uniqueness result for this sublinear problem.

Theorem 3.2 Assume that each of g i j (x) satisfies ( A 1 ), ( A 2 ) and

( A 3 ) lim x g i j ( x ) x =0.

  1. 1.

    If at least one of g i j (0) (i=1,2,,n) is positive and matrix G=( g i j (0)) is irreducible, then (3.1) has a unique positive solution U(λ)=( u 1 (λ), u 2 (λ),, u n (λ)) for all λ>0;

  2. 2.

    If g i j (0)=0, g i j (0)0 for each i,j=1,2,,n and matrix G =( g i j (0)) is irreducible, then for some λ >0, (3.1) has no positive solution when λ λ , and (3.1) has a unique positive solution U(λ)=( u 1 (λ), u 2 (λ),, u n (λ)) for λ> λ .

Moreover, {(λ, u 1 (λ), u 2 (λ),, u n (λ)):λ> λ } (in the first case, we assume λ =0) is a smooth curve so that u i (λ) is strictly increasing in λ, and u i (λ)0 as λ λ + .

Proof Our proof follows that of Theorem 6.1 in [4]. Firstly, we extend g i j to be defined on R and they are continuously differentiable on R. From ( A 2 ), g i j (x)> g i j (x)x implies that j = 1 n g i j ( u j )> j = 1 n g i j ( u j ) u j , so f i (U)= f i ( u 1 ,, u n )= j = 1 n g i j ( u j ) satisfies ( H 1 ). Hence from Theorem 2.4, any positive solution of (3.1) is stable.

Let us define

F(λ,U)= ( Δ u 1 + λ [ g 11 ( u 1 ) + g 12 ( u 2 ) + + g 1 n ( u n ) ] Δ u 2 + λ [ g 21 ( u 1 ) + g 22 ( u 2 ) + + g 2 n ( u n ) ] Δ u n + λ [ g n 1 ( u 1 ) + g n 2 ( u 2 ) + + g n n ( u n ) ] ) ,
(3.3)

where λR and u 1 , u 2 ,, u n C 0 2 , α ( Ω ¯ ). Here g i j are at least C 1 , then F:R×XY is continuously differentiable, where X= [ C 0 2 , α ( Ω ¯ ) ] n and Y= [ C α ( Ω ¯ ) ] n . For weak solutions U=( u 1 ,, u n ), one can also consider X= [ W 2 , p ( Ω ) W 0 1 , p ( Ω ) ] n and Y= [ L p ( Ω ) ] n where p>1 is properly chosen.

It is easy to see that (λ,U)=(0,0,,0) is a solution of (3.1). We apply the implicit function theorem at (λ,U)=(0,0,,0). The Fréchet derivative of F is given by

F U (λ,U) ( ϕ 1 ϕ 2 ϕ n ) = ( Δ ϕ 1 + λ [ g 11 ( u 1 ) ϕ 1 + g 12 ( u 2 ) ϕ 2 + + g 1 n ( u n ) ϕ n ] Δ ϕ 2 + λ [ g 21 ( u 1 ) ϕ 1 + g 22 ( u 2 ) ϕ 2 + + g 2 n ( u n ) ϕ n ] Δ ϕ n + λ [ g n 1 ( u 1 ) ϕ 1 + g n 2 ( u 2 ) ϕ 2 + + g n n ( u n ) ϕ n ] ) .
(3.4)

Then F U (0,0,,0) ( ϕ 1 , , ϕ n ) T = ( Δ ϕ 1 , , Δ ϕ n ) T is an isomorphism from X to Y, and the implicit function theorem implies that F(λ,U)=0 has a unique solution (λ,U(λ)) for λ(0,δ) for some small δ>0, and ( u 1 (0),, u n (0)) is the unique solution of

{ Δ ϕ 1 + g 11 ( 0 ) + g 12 ( 0 ) + + g 1 n ( 0 ) = 0 , x Ω , Δ ϕ 2 + g 21 ( 0 ) + g 22 ( 0 ) + + g 2 n ( 0 ) = 0 , x Ω , Δ ϕ n + g n 1 ( 0 ) + g n 2 ( 0 ) + + g n n ( 0 ) = 0 , x Ω , ϕ 1 ( x ) = = ϕ n ( x ) = 0 , x Ω .
(3.5)

Therefore,

{ u 1 ( 0 ) = ( g 11 ( 0 ) + g 12 ( 0 ) + + g 1 n ( 0 ) ) e u 2 ( 0 ) = ( g 21 ( 0 ) + g 22 ( 0 ) + + g 2 n ( 0 ) ) e u n ( 0 ) = ( g n 1 ( 0 ) + g n 2 ( 0 ) + + g n n ( 0 ) ) e ,

where e is the unique positive solution of

Δe+1=0,xΩ,e(x)=0,xΩ.
(3.6)

If j = 1 n g i j (0)>0 for any i=1,2,,n, then U(λ)=( u 1 (λ),, u n (λ)) is positive for λ(0,δ). If there exists i such that g i 1 (0)+ g i 2 (0)++ g i n (0)>0, then u i (λ)>0 for λ(0,δ). For ki, Δ u k (λ)=λ[ g k 1 (λ)++ g k n (λ)] and g i j is positive, hence u k (λ)>0 as well.

Next we assume that g i j (0)=0 and g i j (0)>0 (i,j=1,,n). The linearized operator is

F U ( λ , 0 , , 0 ) ( Φ 1 Φ 2 Φ n ) = ( Δ Φ 1 + λ [ g 11 ( 0 ) Φ 1 + g 12 ( 0 ) Φ 2 + + g 1 n ( 0 ) Φ n ] Δ Φ 2 + λ [ g 21 ( 0 ) Φ 1 + g 22 ( 0 ) Φ 2 + + g 2 n ( 0 ) Φ n ] Δ Φ n + λ [ g n 1 ( 0 ) Φ 1 + g n 2 ( 0 ) Φ 2 + + g n n ( 0 ) Φ n ] ) = ( Δ Φ 1 Δ Φ 2 Δ Φ n ) + λ ( g 11 ( 0 ) g 12 ( 0 ) g 1 n ( 0 ) g 21 ( 0 ) g 22 ( 0 ) g 2 n ( 0 ) g n 1 ( 0 ) g n 2 ( 0 ) g n n ( 0 ) ) ( Φ 1 Φ 2 Φ n ) = ( Δ Φ 1 Δ Φ 2 Δ Φ n ) + λ J ( Φ 1 Φ 2 Φ n ) ,
(3.7)

where

J= ( g 11 ( 0 ) g 1 n ( 0 ) g n 1 ( 0 ) g n n ( 0 ) ) .

Since g i j (0)>0, all entries of matrix J are positive. Therefore, by using Lemma 3.1, there exist a positive principal eigenvalue χ J and the corresponding eigenvector ( k 1 , k 2 , , k n ) T of J for some k i >0 such that ( k 1 φ 1 , k 2 φ 1 , , k n φ 1 ) T is a positive eigenvector of F ( u 1 , , u n ) ( λ ,0,,0), where λ = λ 1 / χ J . Similarly, J T has the same principal eigenvalue χ J , and the corresponding eigenvector ( k 1 φ 1 , k 2 φ 1 , , k n φ 1 ) T , where k i (i=1,,n) is a positive constant.

Hence when λ= λ = λ 1 / χ J , F U (λ,0,,0) is not invertible and λ= λ is a potential bifurcation point. More precisely, the null space N( F U ( λ ,0,,0))=span{( k 1 φ 1 ,, k n φ n )} is one-dimensional.

Suppose that ( h 1 , , h n ) T R( F U ( λ ,0,,0)), then there exist ( ψ 1 ,, ψ n )X such that

F U ( λ , 0 , , 0 ) ( ψ 1 ψ n ) = ( Δ ψ 1 Δ ψ n ) + λ ( g 11 ( 0 ) g 1 n ( 0 ) g n 1 ( 0 ) g n n ( 0 ) ) ( ψ 1 ψ n ) = ( h 1 h n ) .
(3.8)

Consider the adjoint eigenvalue equation

( Δ ψ 1 Δ ψ n ) + λ J T ( ψ 1 ψ n ) = ( Δ ψ 1 Δ ψ n ) + λ ( g 11 ( 0 ) g n 1 ( 0 ) g 1 n ( 0 ) g n n ( 0 ) ) ( ψ 1 ψ n ) =0,
(3.9)

where ( ψ 1 , , ψ n ) T = ( k 1 φ 1 , , k n φ n ) T .

Multiplying the system (3.8) by ( ψ 1 ,, ψ n ), the system (3.9) by ( ψ 1 ,, ψ n ), integrating over Ω and subtracting, we get

Ω ( h 1 ψ 1 + + h n ψ n ) dx= Ω ( k 1 h 1 φ 1 + + k n h n φ 1 ) dx=0.
(3.10)

Hence ( h 1 , , h n ) T R( F U ( λ ,0,,0)) if and only if (3.9) holds, which implies that the codimension of R( F U ( λ ,0,,0)) is one.

Next we verify that F λ U ( λ ,0,,0) [ k 1 φ 1 , , k n φ 1 ] T R( F U ( λ ,0,,0)). Indeed,

F λ U ( λ ,0,,0) ( k 1 φ 1 k n φ 1 ) =J ( k 1 k n ) φ 1 = χ J ( k 1 k n ) .
(3.11)

But this contradicts with

0= χ J Ω ( k 1 k 1 + + k n k n ) φ 1 2 dx>0.
(3.12)

Hence

F λ U ( λ ,0,,0) [ k 1 φ 1 , , k n φ 1 ] T R ( F U ( λ , 0 , , 0 ) ) .

By using a bifurcation from a simple eigenvalue theorem of Crandall-Rabinowitz [24], we conclude that ( λ ,0,,0) is a bifurcation point. The nontrivial solutions of F(λ, u 1 ,, u n )=(0,,0) near ( λ ,0,,0) are in the form of {(λ(s), u 1 (s),, u n (s))} for s(δ,δ), where u i (s)= k i s φ 1 +o(s) (i=1,,n). From the stability of positive solutions, each positive solution is non-degenerate.

Next we claim that (1.1) has no positive solution when λ λ . Let ( u 1 ,, u n ) is a positive solution of (1.1) and ( k 1 φ 1 ,, k n φ 1 ) satisfy

( Δ k 1 φ 1 Δ k n φ 1 ) + λ ( g 11 ( 0 ) g n 1 ( 0 ) g 1 n ( 0 ) g n n ( 0 ) ) ( k 1 φ 1 k n φ 1 ) =0.
(3.13)

Multiplying the system (3.1) by ( k 1 φ 1 , k 2 φ 1 ,, k n φ n ), the system (3.13) by ( u 1 ,, u n ), integrating over Ω and subtracting, by ( A 2 ) and j = 1 n g i j (0)=0, we obtain

λ Ω i = 1 n [ j = 1 n k i g i j ( 0 ) u j φ 1 ] d x = λ Ω j = 1 n [ i = 1 n k i g i j ( u j ) φ 1 ] d x < λ Ω j = 1 n [ i = 1 n k i g i j ( 0 ) u j φ 1 ] d x .
(3.14)

Hence (3.1) has no positive solution when λ λ . And the solution {(λ(s), u 1 (s),, u n (s))} can also be parameterized as (λ, u 1 (λ),, u n (λ)) for λ( λ , λ +δ). With the implicit function theorem, we can extend this curve to the largest λ .

Let Γ={(λ, u 1 (λ),, u n (λ)): λ <λ< λ }. We show that ( u 1 (λ),, u n (λ)) is strictly increasing in λ for λ( λ , λ ). In fact, ( u 1 (λ)/λ,, u n (λ)/λ) satisfies the equation

F U ( λ , U ) ( u 1 ( λ ) λ u n ( λ ) λ ) = ( Δ u 1 ( λ ) λ Δ u n ( λ ) λ ) + λ ( g 11 ( u 1 ( λ ) ) g 1 n ( u n ( λ ) ) g n 1 ( u 1 ( λ ) ) g n n ( u n ( λ ) ) ) ( u 1 ( λ ) λ u n ( λ ) λ ) = ( g 11 ( u 1 ( λ ) ) + + g 1 n ( u n ( λ ) ) g n 1 ( u 1 ( λ ) ) + + g n n ( u n ( λ ) ) ) .
(3.15)

Hence ( u 1 (λ)/λ,, u n (λ)/λ) from the maximum principle (Lemma 2.2 part 3) and the fact that μ 1 (( u 1 (λ),, u n (λ)))>0 from the stability of positive solutions. We claim that λ =. Suppose not, then λ < and ( u 1 ( λ ) , , u n ( λ ) ) X <. Then one can show that the curve Γ can be extended to λ= λ from some standard elliptic estimates, then from the implicit function theorem, Γ can be extended beyond λ= λ , which is a contradiction; if λ < and ( u 1 ( λ ) , , u n ( λ ) ) X =, a contradiction can be derived with the solution curve which cannot blow up at finite λ (see similar arguments for the scalar equation in [25]). Hence we must have λ =.

Finally, we claim the uniqueness. If there is another positive solution for some λ> λ , then the arguments above show that this solution also belongs to a solution curve defined for λ( λ ,), and the solutions on the curve are increasing in λ, but the nonexistence of positive solutions for λ< λ and the local bifurcation at λ= λ exclude the possibility of another solution curve. Hence the positive solution is unique for all λ> λ . This completes the proof. □

Example 3.3

We consider the following cyclic system:

{ Δ u 1 + λ f 1 ( u 2 ) = 0 , x Ω , Δ u 2 + λ f 2 ( u 3 ) = 0 , x Ω , Δ u n + λ f n ( u 1 ) = 0 , x Ω , u 1 ( x ) = = u n ( x ) = 0 , x Ω ,
(3.16)

where f i are smooth real functions defined on R + satisfying ( A 1 ), ( A 2 ) and ( A 3 ). Then Theorem 3.2 implies the existence and uniqueness of a positive solution of problem (3.16), and the bifurcation diagram is a single monotone solution curve. The n-dimensional cyclic positone and semipositone system was considered in [26] and [27]. They got the existence and multiplicity of positive solutions result for some combined sublinear condition by the method of sub-super solutions.

Example 3.4

Consider the general Lane-Emden system

{ Δ u 1 + λ b 1 ( a 11 + u 1 ) p 11 ( a 1 n + u n ) p 1 n = 0 , x Ω , Δ u i + λ b i ( a i 1 + u 1 ) p i 1 ( a i n + u n ) p i n = 0 , x Ω , Δ u i + λ b n ( a n 1 + u 1 ) p n 1 ( a n n + u n ) p n n = 0 , x Ω , u 1 ( x ) = = u n ( x ) = 0 , x Ω ,
(3.17)

where λ>0 is a positive parameter, b i , a i j (1i,jn) are positive constants, p i j (1i,jn) are nonnegative constants satisfying i = 1 n p i j <1 and Ω R N denotes a bounded domain of class C 2 , α , α(0,1). For the case N=2, (3.17) has been studied by many authors. Especially, Dalmasso [10] proved the uniqueness and existence of positive solutions to (3.17) for the case N=2, p 11 = p 22 =0, 0< p 12 p 21 <1.

In this section, we show the uniqueness and existence of positive solutions for (3.17) by using super-subsolution methods and the stability of positive solutions by using Theorem 2.4.

Let e(x) be the unique positive solution of

Δe+1=0,xΩ,e(x)=0,xΩ.

We construct a super-solution ( u ¯ 1 ,, u ¯ n )=M(e(x),,e(x)). There exists a suitable positive constant M such that

Δ ( M e ) + λ b i ( a i 1 + M e ) p i 1 ( a i n + M e ) p i n M + λ b i ( a ¯ i + M e ) p i 1 + + p i n 0 ,

where a ¯ i =max{ a i 1 ,, a i n }. We construct a sub-solution in the form of ( u ̲ 1 ,, u ̲ n )=(ε φ 1 ,,ε φ 1 ), where ε will be specified later. Recall that φ 1 is the positive principal eigenfunction with φ 1 =1. Now, for the i th equation of (3.17), we have

Δ ( ε φ 1 ) + λ b i ( a i 1 + ε φ 1 ) p i 1 ( a i n + ε φ 1 ) p i n ε λ 1 φ 1 + λ b i ( a ̲ i + ε φ 1 ) p i 1 + + p i n ε λ 1 φ 1 + λ b i ( ε φ 1 ) p i 1 + + p i n ,

where a ̲ i =min{ a i 1 ,, a i n }. Hence ( u ̲ 1 ,, u ̲ n ) is a subsolution of (3.17), if we choose ε smaller, and it satisfies ε< ( λ b i / λ 1 ) 1 1 ( p i 1 + + p i n ) φ 1 1 . Hence if we choose ε smaller so that ( u ̲ 1 ,, u ̲ n )( u ¯ 1 ,, u ¯ n ), then ( u ̲ , v ̲ ) is a subsolution of (4.1). Therefore, there exists a positive solution ( u 1 ,, u n ) of (3.17) between the supersolution and the subsolution.

Next, we show the solution is stable. Letting f i (U)= b i j = 1 n ( a i j + u j ) p i j , by simple calculation, we get

f i j (U)= b i p i j l = 1 , l j n ( a i l + u l ) p i l ( a i j + u i j ) p i j 1

and

f i ( U ) j = 1 n f i j ( U ) u j = b i j = 1 n ( a i j + u j ) p i j j = 1 n b i p i j l = 1 , l j n ( a i l + u l ) p i l ( a i j + u j ) p i j 1 = b i l = 1 n ( a i l + u l ) p i l [ 1 j = 1 n p i j a i j + u j u j ] > b i l = 1 n ( a i l + u l ) p i l [ 1 j = 1 n p i j ] > 0 .

So, f i satisfies ( H 1 ), hence from Theorem 2.4, any positive solution of (3.17) is stable.

4 Application: Hölder continuous case

In this section, we consider that f i C α ( R ¯ + n ) C 1 ( R + n ) for α(0,1).

Example 4.1

Consider

{ Δ u + λ ( a u u 2 + u v ) = 0 , x Ω , Δ v + λ ( b v v 2 + u v ) = 0 , x Ω , u ( x ) = v ( x ) = 0 , x Ω ,
(4.1)

where λ>0 is a positive parameter, a,b0.

We use the monotone method to prove the existence of a solution. For the supersolution, we choose ( u ¯ , v ¯ )=(a+b+1,a+b+1). Then

Δ u ¯ + λ ( a u ¯ u ¯ 2 + u ¯ v ¯ ) λ [ a ( a + b + 1 ) ( a + b + 1 ) 2 + ( a + b + 1 ) ] λ ( a + b + 1 ) [ ( a + 1 ) ( a + b + 1 ) ] 0 .

Similarly, we have 0<v<(a+b+1). Then it is clear that ( u ¯ , v ¯ ) is a supersolution of (4.1). We construct a sub-solution in the form of ( u ̲ , v ̲ )=(ε φ 1 ,ε φ 1 ), where ε will be specified later. Recall that φ 1 is the positive principal eigenfunction with φ 1 =1. Now, for the equation of u or the equation of v, we have

Δ ( ε φ 1 ) + λ ( a ε φ 1 ( ε φ 1 ) 2 + ε φ 1 ) = ε λ 1 φ 1 + λ ( a ε φ 1 ( ε φ 1 ) 2 + ε φ 1 ) = ε φ 1 [ λ ( a + 1 ) λ 1 λ ε φ 1 ] .

Hence if we choose ε smaller and it satisfies ε<min{ λ ( a + 1 ) λ 1 λ φ 1 , λ ( b + 1 ) λ 1 λ φ 1 }, so that ( u ̲ , v ̲ )( u ¯ , v ¯ ), then ( u ̲ , v ̲ ) is a subsolution of (4.1). Therefore, there exists a positive solution (u,v) of (4.1) between the supersolution and the subsolution when λ>max{ λ 1 /(a+1), λ 1 /(b+1)}. We remark that the stability defined in Section 2 can still be established for a nonlinear function f i to become ∞ near Ω by using Remark 3.1 in [28]. Thus the positive solution (u,v) of (4.1) is stable.

Remark 4.2 Since (4.1) is a cooperative model from ecology with logistic growth rate and sublinear interaction term, we can get the stable result. When the interaction terms are uv (Lotka-Volterra type) and u v as proposed here, and they do not satisfy the conditions of Theorem 3.2, thus, the solution may not be unique.

References

  1. Chen ZY, Chern JL, Shi J, Tang YL: On the uniqueness and structure of solutions to a coupled elliptic system. J. Differ. Equ. 2010, 249: 3419-3442.

    Article  MathSciNet  Google Scholar 

  2. Chern JL, Tang YL, Lin CS, Shi J: Existence, uniqueness and stability of positive solutions to sublinear elliptic systems. Proc. R. Soc. Edinb. A 2011, 141: 45-64.

    Article  MathSciNet  Google Scholar 

  3. Cui, RH, Li, P, Shi, J, Wang, YW: Existence, uniqueness and stability of positive solutions for a class of semilinear elliptic systems. Topol. Methods Nonlinear Anal. (to appear)

  4. Shi, J, Shivaji, R: Exact multiplicity of positive solutions to cooperative elliptic systems. Preprint (2012)

    Google Scholar 

  5. Korman P, Li Y, Ouyang T: Exact multiplicity results for boundary value problems with nonlinearities generalising cubic. Proc. R. Soc. Edinb. A 1996, 126: 599-616.

    Article  MathSciNet  Google Scholar 

  6. Korman P, Li Y, Ouyang T: An exact multiplicity result for a class of semilinear equations. Commun. Partial Differ. Equ. 1997, 22: 661-684.

    Article  MathSciNet  Google Scholar 

  7. Ouyang T, Shi J: Exact multiplicity of positive solutions for a class of semilinear problem. J. Differ. Equ. 1998, 146: 121-156.

    Article  MathSciNet  Google Scholar 

  8. Ouyang T, Shi J: Exact multiplicity of positive solutions for a class of semilinear problem, II. J. Differ. Equ. 1999, 158: 94-151.

    Article  MathSciNet  Google Scholar 

  9. Sirakov B: Notions of sublinearity and superlinearity for nonvariational elliptic systems. Discrete Contin. Dyn. Syst. 2005, 13: 163-174.

    Article  MathSciNet  Google Scholar 

  10. Dalmasso R: Existence and uniqueness of positive solutions of semilinear elliptic systems. Nonlinear Anal. 2000, 39: 559-568.

    Article  MathSciNet  Google Scholar 

  11. Hai DD: Existence and uniqueness of solutions for quasilinear elliptic systems. Proc. Am. Math. Soc. 2005, 133: 223-228.

    Article  MathSciNet  Google Scholar 

  12. Hai DD: Uniqueness of positive solutions for semilinear elliptic systems. J. Math. Anal. Appl. 2006, 313: 761-767.

    Article  MathSciNet  Google Scholar 

  13. Hai DD, Shivaji R: An existence result on positive solutions for a class of semilinear elliptic systems. Proc. R. Soc. Edinb. A 2004, 134: 137-141.

    Article  MathSciNet  Google Scholar 

  14. Troy WC: Symmetry properties in systems of semilinear elliptic equations. J. Differ. Equ. 1981, 42: 400-413.

    Article  MathSciNet  Google Scholar 

  15. Serrin J, Zou H: Existence of positive solutions of the Lane-Emden system. Atti Semin. Mat. Fis. Univ. Modena 1998, 46: 369-380. suppl.

    MathSciNet  Google Scholar 

  16. Serrin J, Zou H: Existence of positive entire solutions of elliptic Hamiltonian systems. Commun. Partial Differ. Equ. 1998, 23: 577-599.

    MathSciNet  Google Scholar 

  17. Korman P, Shi J: On Lane-Emden type systems. Discrete Contin. Dyn. Syst. 2005, suppl.: 510-517.

    MathSciNet  Google Scholar 

  18. Cui RH, Wang YW, Shi J: Uniqueness of the positive solution for a class of semilinear elliptic systems. Nonlinear Anal. 2007, 67: 1710-1714.

    Article  MathSciNet  Google Scholar 

  19. Cui RH, Shi J, Wang YW: Existence and uniqueness of positive solutions for a class of semilinear elliptic systems. Acta Math. Sin. Engl. Ser. 2011, 27: 1079-1090.

    Article  MathSciNet  Google Scholar 

  20. Masaaki M: Uniqueness and existence of positive solutions for some semilinear elliptic systems. Nonlinear Anal. 2004, 59: 993-999.

    Article  MathSciNet  Google Scholar 

  21. Sweers G:Strong positivity in C( Ω ¯ ) for elliptic systems. Math. Z. 1992, 209: 251-271.

    Article  MathSciNet  Google Scholar 

  22. Kato T Classics in Mathematics. In Perturbation Theory for Linear Operators. Springer, Berlin; 1995. Reprint of the 1980 edition.

    Google Scholar 

  23. Serre D Graduate Texts in Mathematics 216. In Matrices: Theory and Applications. Springer, Berlin; 2002.

    Google Scholar 

  24. Crandall MG, Rabinowitz PH: Bifurcation from simple eigenvalues. J. Funct. Anal. 1971, 8: 321-340.

    Article  MathSciNet  Google Scholar 

  25. Shi J: Blow up points of solution curves for a semilinear problem. Topol. Methods Nonlinear Anal. 2000, 15: 251-266.

    MathSciNet  Google Scholar 

  26. Ali J, Brown K, Shivaji R:Positive solutions for n×n elliptic systems with combined nonlinear effects. Differ. Integral Equ. 2011, 24(3-4):307-324.

    MathSciNet  Google Scholar 

  27. Lee EK, Shivaji R, Ye J:Classes of infinite semipositone n×n systems. Differ. Integral Equ. 2011, 24(3-4):361-370.

    MathSciNet  Google Scholar 

  28. Ambrosetti A, Brézis H, Cerami G: Combined effects of concave and convex nonlinearities in some elliptic problems. J. Funct. Anal. 1994, 122: 519-543.

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

Partially supported by the National Natural Science Foundation of China (No. 11071051, 11271100 and 11101110), the Aerospace Supported Fund, China, under Contract (Grant 2011-HT-HGD-06), Science Research Foundation of the Education Department of Heilongjiang Province (Grant No. 12521153), Science Foundation of Heilongjiang Province (Grant A201009) and Harbin Normal University advanced research Foundation (11xyg-02).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Renhao Cui.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

BW and RC carried out the proof of the main part of this article, RC corrected the manuscript and participated in its design and coordination. All authors have read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wu, B., Cui, R. Existence, uniqueness and stability of positive solutions to a general sublinear elliptic systems. Bound Value Probl 2013, 74 (2013). https://doi.org/10.1186/1687-2770-2013-74

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-74

Keywords