Startseite On a class of critical elliptic systems in ℝ4
Artikel Open Access

On a class of critical elliptic systems in ℝ4

  • Xin Zhao EMAIL logo und Wenming Zou
Veröffentlicht/Copyright: 13. September 2020

Abstract

In the present paper, we consider the following classes of elliptic systems with Sobolev critical growth:

Δu+λ1u=μ1u3+βuv2+2qpyu2qp1v2inΩ,Δv+λ2v=μ2v3+βu2v+2yu2qpvinΩ,u,v>0inΩ,u,v=0onΩ,

where Ω ⊂ ℝ4 is a smooth bounded domain with smooth boundary ∂Ω; p, q are positive coprime integers with 1 < 2qp < 2; μi > 0 and λi ∈ ℝ are fixed constants, i = 1, 2; β > 0, y > 0 are two parameters. We prove a nonexistence result and the existence of the ground state solution to the above system under proper assumptions on the parameters. It seems that this system has not been explored directly before.

MSC 2010: 35J50; 35J15; 35J60

1 Introduction

In the current paper, we consider the following system:

Δu+λ1u=μ1u3+βuv2+2qpyu2qp1v2inΩ,Δv+λ2v=μ2v3+βu2v+2yu2qpvinΩ,u,v>0inΩ,u,v=0onΩ, (1.1)

where Ω ⊂ ℝ4 is a smooth bounded domain with a smooth boundary ∂Ω; p, q are positive coprime integers with 1 < 2qp < 2, μi > 0 and λi ∈ ℝ are fixed constants, i = 1, 2; β > 0, y > 0 are parameters.

In ℝ4, the Sobolev critical exponent 2* = 4, therefore, the system (1.1) has both critical nonlinearities and subcritical coupling perturbations. The study about the solutions of (1.1) when y = 0 arises in the solitary wave solutions of the following Schrödinger equations, known as the Gross-Pitaevskii equations(cf.[13, 20]):

ιtΨ1=ΔΨ1+μ1|Ψ1|2Ψ1+β|Ψ2|2Ψ1,ιtΨ2=ΔΨ2+μ2|Ψ2|2Ψ2+β|Ψ1|2Ψ2,Ψi=Ψi(t,x)H1(RN;C),i=1,2,N=1,2,3. (1.2)

Here ι is the imaginary unit. This system was developed as a mathematical model for Bose-Einstein condensation in hyperfine spin states. To obtain solitary wave solutions, we set Ψi(t, x) = eιtλiui(x) for i = 1, 2, then ui(x) satisfy the elliptic system:

Δu1+λ1u1=μ1u13+βu22u1inΩ,Δu2+λ2u2=μ2u23+βu12u2inΩ,u1,u2>0inΩ,u1=u2=0onΩ. (1.3)

Since the solitary waves of (1.2) arise in many physical problems, especially in the Hartree-Fock theory and nonlinear optics, the system (1.3) in low dimensions (1 ≤ N ≤ 3) has been well studied. When N ≤ 3, the nonlinearity and the coupling terms in (1.3) are of subcritical growth. For example, we refer to [11, 14, 16, 20, 21, 25] and the references therein. Recently, multiple results of (1.3) in ℝN with N ≥ 4 have been derived (cf. [10, 19]). In these cases, the nonlinearities and coupled terms are Sobolev critical for N = 4 and even super-critical for N ≥ 5. Therefore, the study of the cases in high dimensions is much more complicated due to the lack of compactness. In [10], the authors proved that the 2-component system (1.3) has a positive least-energy solution when N = 4 for negative β, positive small β and positive large β with μ1, μ2 > 0, –α1(Ω) < λ1, λ2 < 0, where α1(Ω) is the first eigenvalue of –Δ with the Dirichlet boundary condition. In [8, 19], the existence of infinitely many sign-changing solutions for the k-component system (1.3) for N ≥ 2, k ≥ 2 was obtained.

From the perspective of mathematics, we also note that the related system as (1.1) was studied by [6] in 1984. Actually, in [6] it was considered the systems of equations

Δui(x)=gi(u(x)),i=1,2,...N0,

where the functions gi are the gradients of some functions. Under some hypotheses on gi, the existence results of the ground state solution were established (see Theorem 2.2 and Theorem 2.3 in [6]).

Furthermore, in [23], the author obtained some results including the existence and nonexistence of the ground states when N = 4 to the following Sobolev critical system, which has both linear perturbations and nonlinear couplings arising from (1.3):

Δu1+λ1u1=μ1u13+βu22u1+yu2inΩ,Δu2+λ2u2=μ2u23+βu12u2+yu1inΩ,u1,u2>0inΩ,u1=u2=0onΩ, (1.4)

here β > 0 (attractive), y > 0. Such kinds of systems arise in nonlinear optics (cf. [1]). There are also some results on the linearly coupled system like (1.4). For example, in [2] the multi-bump solitons were obtained in the case of the dimension N ≤ 3. For more general nonlinear terms and N = 1, we refer to [3] (see also [4] for N ≥ 2).

Our interest in (1.1) mainly comes from the challenge in mathematics since this system has not been explored directly before. An idea of dealing with the lack of compactness contributed by Brezis and Nirenberg in [7], is to establish a threshold for the compactness of the (PS) sequence. Recall that the methods in [9, 10] for the critical systems with either only nonlinear couplings or only linear couplings cannot be applied in the present paper, because of the existence with both critical and subcritical nonlinear couplings in our case (1.1). We will devote to obtain a uniqueness result to the related algebraic equation when β > 0 is properly large, and to estimate the least energy of the related functional by exploring the minimizing problem (3.2) (see Lemma 3.1) below.

Let us present the necessary notations and definitions before stating the main results for (1.1). Set 𝓗 := H01 (Ω) × H01 (Ω). Define the functional I : 𝓗 → ℝ by

I(u,v):=12Ω(|u|2+|v|2+λ1u2+λ2v2)14Ω(μ1u4+μ2v4+2βu2v2)yΩu2qpv2. (1.5)

Then I(u, v) is of class C1 in 𝓗 and the positive solutions of (1.1) correspond to the positive critical points of I(u, v). We define the Nehari manifold

N:={(u,v)H{(0,0)}:F(u,v):=I(u,v)(u,v)=0}, (1.6)

where

F(u,v)=Ω(|u|2+|v|2+λ1u2+λ2v2)Ω(μ1u4+μ2v4+2βu2v2)(2qp+2)yΩu2qpv2. (1.7)

Then any nontrivial positive solution of (1.1) belongs to 𝓝. It is not hard to check that 𝓝 ≠ ∅. Define

B:=inf(u,v)NI(u,v). (1.8)

Since

I|N(u,v)=q2(p+q)Ω(|u|2+|v|2+λ1u2+λ2v2)+pq4(p+q)Ω(μ1u4+μ2v4+2βu2v2)

here μi > 0 and β ≥ 0, then I is bounded from below on 𝓝 with 0 < 2qp < 2.

Definition 1.1

The pair (u, v) is called a nontrivial solution iff u ≠ 0, v ≠ 0; we call (u, v) is positive iff u > 0, v > 0; the nontrivial solution (u0, v0) is called a ground state solution to (1.1) iff I(u0, v0) ≤ I(u, v) for all nontrivial critical points (u, v).

Firstly, we consider a little relaxed system:

Δu+λ1u=μ1u3+βuv2+2qpyu2qp1v2inΩ,Δv+λ2v=μ2v3+βu2v+2yu2qpvinΩ,u=v=0onΩ, (1.9)

where u, v H01 (Ω) and y, λ1, λ2 ∈ ℝ; μ1, μ2 > 0. The first conclusion in this paper is about the nonexistence of solution to the systems (1.1) and (1.9).

Theorem 1.2

Let α1(Ω) > 0 be the first eigenvalue ofΔ in H01 (Ω).

  1. Assume μ1, μ2, β, y > 0. If min(λ1, λ2) ≤ –α1(Ω), then system (1.1) has no solution;

  2. Assume Ω is star-shaped or Ω = ℝ4. If y ≤ 0 and min(λ1, λ2) > 0, then (1.9) has no nontrivial solution.

Before stating the next main theorem, we recall the Brezis-Nirenberg type critical exponent equation (see [7]):

Δu+λiu=μiu3,u0,uH01(Ω),

it has a positive least energy solution uμiC2(Ω) ∩ C(Ω̄) when –α1(Ω) < λi < 0, N = 4 with energy

Bμi:=12Ω(|uμi|2+λiuμi2)14Ωμiuμi4=14Ω(|uμi|2+λiuμi2)=14Ωμiuμi4>0, (1.10)

and there is no solution when λi ≥ 0 and Ω is a (smooth) starshaped domain.

Remark 1.3

Assume thatα1(Ω) < λi < 0. Let 𝓢 be the sharp constant of D1,2(ℝ4) ↪ L4(ℝ4). Then Bμi mentioned in (1.10) satisfies (See [9]):

1N(α1(Ω)+λiα1(Ω))N2μi2N2SN2Bμi<1Nμi2N2SN2,i=1,2.

Define

Nμi:={uH01(Ω):Ω(|u|2+λiu2μiu4)=0}. (1.11)

Let C(Ω) be the embedding constant of H01 (Ω) to L2(Ω), i.e.,

Ω|f|2(C(Ω))2Ω|f|2. (1.12)

Here, C(Ω)=(1α1(Ω))12.

Theorem 1.4

Let λ1, λ2 > –α1(Ω), μ1, μ2, β, y > 0. Assume that λ1 satisfies one of the following two conditions:

  1. λ1 < 0 and β>(Ω|v0|2+λ22v02)2y(Ωuμ12qpv02)Ωuμ12v02, for some 0 ≠ v0 H01 (Ω);

  2. λ1 > 0 with Ω is a (smooth) starshaped domain,

and that λ2 satisfies one of the following two conditions:

  1. λ2 < 0;

  2. λ2 > 0 with Ω is a (smooth) starshaped domain.

Then the system (1.1) has a ground state solution provided that one of the following cases holds:

  1. min(λ1, λ2) < 0 < max(λ1, λ2) and β > max(μ1, μ2);

  2. min(λ1, λ2) > 0 and β > max(μ1, μ2);

  3. max(λ1, λ2) < 0 and β > μ1μ2.

Remark 1.5

If, for example, (11) and (21) hold, then only the case (c) occurs by default.

The rest of this paper will be denoted to the proofs of these theorems. We denote the norm of Lp(Ω) by up=(Ω|u|pdx)1p, the norm of H01 (Ω) by ∥u∥ = ∥∇u2 and positive constants (possibly different) by C. The paper is organized as follows. In section 2, we study the nonexistence of solutions by establishing a kind of Pohozaev’s type identity. Then in Section 3, we shall show some preliminary results including energy estimation. In section 4, we will use Ekeland variational principle and the results of Section 3 to prove Theorem 1.4.

2 Nonexistence

Lemma 2.1

[Pohozaev identity] Suppose u, v are (smooth) functions satisfying

Δu=g(u,v)inΩ,Δv=f(u,v)inΩ,u,v=0onΩ, (2.1)

where Ω ⊂ ℝn(n ≥ 3); g and f are continuous on2. Then we have

Ω(n(G+F)+(Gvv+Fuu)x)=n2(u2+v2)+12Ω(|u|2+|v|2)(xν), (2.2)

where 2=2nn2,G(u,v)=0ug(t,v)dt,F(u,v)=0vf(u,t)dt and ν⃗ denotes the outward normal to ∂Ω.

Proof

We multiply the two equations of (2.1) by (x⃗ ⋅ ∇u), (x⃗ ⋅ ∇v) respectively and integrate over Ω. Then add them up, we have

Ω((xu)Δu(xv)Δv)=Ω(g(u,v)(xu)+f(u,v)(xv)). (2.3)

By Pohozaev identity (cf. [7]), we have

ΩΔu(xu)=(1n2)Ω|u|212Ω|u|2(xν)ds,ΩΔv(xv)=(1n2)Ω|v|212Ω|v|2(xν)ds (2.4)

and

Ωg(u,v)(xu)=Ωi=1ng(u,v)xiuxi=Ωi=1nGxixii=1nGvvxixi=nΩG(u,v)ΩGv(xv). (2.5)

Similarly, we have

Ωf(u,v)(xv)=nΩF(u,v)ΩFu(xu). (2.6)

Then plug (2.4), (2.5) and (2.6) into (2.3), we obtain (2.2), which completes the proof.□

Proof of Theorem 1.2

  1. Recall system (1.1), we know that β > 0, y > 0, μ1 > 0, μ2 > 0. If λ1 = min(λ1, λ2) ≤ –α1 (Ω) and let φ be the corresponding eigenfunction of α1 (Ω). Suppose that (u, v) is a solution of system (1.1), i.e., u > 0, v > 0. Then we multiply (1.9) with (φ, 0) and integrate by parts, we obtain

    0(λ1+α1(Ω))Ωuφ=Ωuφ+λuφ=Ω(μ1u3+βuv2+2qpyu2qp1v2)φ>0,

    which is impossible. Similarly, if λ2 = min (λ1, λ2) ≤ –α1(Ω) and we multiply (1.1) with (0, φ), we have

    0(λ2+α1(Ω))Ωvφ=Ω(μ2v3+βu2v+2yu2qpv)φ>0,

    which leads to a contradiction. So system (1.1) has no solution.

  2. Let (u, v) be a solution of the system (1.9). We take

    G(u,v)=λ12u2+μ14u4+β2u2v2+yu2qpv2,F(u,v)=λ22v2+μ24v4+β2u2v2+yu2qpv2.

    By the Pohozaev identity in Lemma 2.1, we have

    Ω(|u|2+|v|2+2(λ1u2+λ2v2)(μ1u4+μ2v4+2βu2v2+4yu2qpv2))=12Ω(|u|2+|v|2)(xν)ds. (2.7)

    One the other hand,

    Ω(|u|2+|v|2+λ1u2+λ2v2(μ1u4+μ2v4+2βu2v2)(2qp+2)yu2qpv2)=0. (2.8)

    Therefore,

    Ω(λ1u2+λ2v2(22qp)yu2qpv2)=12Ω(|u|2+|v|2)(xν)ds.

    When one of the following two conditions holds: (a) Ω = ℝ4 or (b) Ω is star-shaped, i.e., (x⃗ν⃗) > 0 a.e. ∂Ω, then we have – 12 ∂Ω(|∇u|2 + |∇v|2)(x⃗ν⃗)ds ≤ 0, it follows that if λ1, λ2 > 0, y ≤ 0, system (1.9) has no nontrivial solution.□

3 Some Preliminaries

Consider the system

Δu=μ1u3+βuv2R4,Δv=μ2v3+βu2vR4,

the corresponding functional is

E(u,v)=12R4(|u|2+|v|2)14R4(μ1u4+μ2v4+2βu2v2) (3.1)

defined in 𝓓 = (D1,2(ℝ4))2, and 𝓓 is a Hilbert space equipped with the inner product

(u1,v1),(u2,v2)R4=R4(u1u2+v1v2).

Define

M={(u,v)D{0}:E(u,v)(u,v)=0}.

We consider the minimizing problem

c:=infME(u,v).

Define

f(t1,t2)=t122+t222μ1t14+μ2t244β2t12t22,P={(t1,t2)R2{(0,0)}:t12+t22(μ1t14+μ2t24)2βt12t22=0}. (3.2)

Repeating the proof of Proposition 3.2 in [23], we have

Lemma 3.1

(see [23]) Let

d:=infPf(t1,t2).

Then c := d𝓢2 is attained by (u, v) if and only if (u, v) = (t1Uε,z, t2Uε,z), where 𝓢 is the best Sobolev embedding constant from D1,2(ℝ4) to L4(ℝ4) and

Uε,z(x)=(8ε2)12ε2+|xz|2

is the Talanti function satisfyingΔU = U3 in4 and (t1, t2) satisfies

t1=μ1t13+βt1t22,t2=μ2t23+βt12t2,t1,t20,t1+t2>0. (3.3)

Moreover, if β > max{μ1, μ2}, then t1, t2 > 0.

Assume 0 ∈ Ω, ρ > 0 such that B(0, 2ρ) ⊂ Ω. Take a radial symmetric cut-off function Ψ C02 (B(0, 2ρ)), 0 ≤ Ψ(x) ≤ 1, Ψ(x) ≡ 1 in B(0, ρ). Define Vε(x) = Ψ(x) ⋅ Uε,0(x), where

Uε,0(x)=(8ε2)12ε2+|x|2.

Then, it is well know (cf. [10]) that

Vε22=S2+o(ε2),Vε44=S2+o(ε4),Vε22Cε2|lnε|+o(ε2). (3.4)

We have to estimate ℬ, which is given by (1.8). To this purpose, we need to prepare Lemma 3.2, Proposition 3.3 and Lemma 3.4. According to the definition of 𝓝, if ∃ a1, a2 > 0 such that the following equation

((μ1a14+2βa12a22+μ2a24)Vε44)k2+((2qp+2)ya12qpa22Vε2qp+22qp+2)k2qp((a12+a22)Vε22+(λ1a12+λ2a22)Vε22)=0 (3.5)

has a solution k, where k is dependent on ε, a1, a2, then (ka1Vε, ka2Vε) ∈ 𝓝.

Define

Gε(k):=((μ1a14+2βa12a22+μ2a24)Vε44)k2+((2qp+2)ya12qpa22Vε2qp+22qp+2)k2qp((a12+a22)Vε22+(λ1a12+λ2a22)Vε22). (3.6)

Then the equation (3.5) is equivalent to Gε(k) = 0. Moreover, Gε (k) ≥ 0, ∀ k ∈ [0, +∞) as μi, β > 0, which means Gε(k) is monotonically increasing on [0, +∞).

Lemma 3.2

Assume μ1, μ2, β, y > 0.

  1. If min(λ1, λ2) < 0 < max(λ1, λ2), thena1, a2 > 0 such that (3.5) has a unique positive solution k andQ < 0 such that

    λ1a12+λ2a222Vε2qp+22qp+2Vε22ya12qpa22k2qpQ<0, (3.7)

    for ρ > 0 and ε > 0 small enough.

  2. If min(λ1, λ2) > 0, thena1, a2 > 0, (3.5) has a unique positive solution k andQ < 0 such that (3.7) holds forε > 0 small enough.

  3. If max(λ1, λ2) < 0, thena1, a2 > 0, (3.5) has a unique positive solution k andQ < 0 such that (3.7) holds forρ, ε > 0 small enough.

Proof

Since Gε(k) ∈ C1(ℝ) and monotonically increasing on [0, +∞) with

limk+Gε(k)=+.

It is sufficient to verify that the unique solution k > 0 exists as long as

Pε:=Gε(0)=(a12+a22)Vε22+(λ1a12+λ2a22)Vε22>0. (3.8)

  1. Without loss of generality, we assume λ1 < 0 < λ2 in this case. Let

    Ai:=Vε22+λiVε22,i=1,2.

    Then we have A1 < A2 and A2 > 0. If there exists a1, a2 > 0 such that

    (a2a1)2>A1A2, (3.9)

    then Pε=a12(A1+(a2a1)2A2)>0 and (3.8) holds, which means (3.5) has positive solution k. Next we will show that ∀ a > 0, a1 = a, a2 = λ12λ2 a satisfy (3.9). In fact, – A1A2 < 0 for ε > 0 small enough. Recall that λ1 < 0 ≤ λ2, A2 > 0 and

    |Vε2|B(0,2ρ)|14Vε4,

    here |B(0, 2ρ)| denotes the measure of the ball B(0, 2ρ), then |B(0, 2ρ)| = 8π2ρ4 when N = 4, we have

    A1A2=Vε22+λ1Vε22Vε22+λ2Vε22S2+o(ε2)+λ1(|B(0,2ρ)|12Vε42)Vε22+λ2Vε22=S2+22πρ2Sλ1+o(ε2)Vε22+λ2Vε22. (3.10)

    Take ρ > 0 small enough and satisfying

    S+22πρ2λ1>0, (3.11)

    here 𝓢 is the best Sobolev embedding constant from D1,2(ℝN) to L4(ℝN). Thus we have 𝓢2 + 2 2 πρ2 λ1𝓢 + o(ε2) > 12 (𝓢2 + 2 2 πρ2λ1𝓢) > 0 when ε > 0 small enough and

    S2+22πρ2λ1S+o(ε2)Vε22+λ2Vε22<12(S2+22πρ2λ1S)Vε22+λ2Vε2212(S2+22πρ2λ1S)S2+o(ε2)+λ2(|B(0,2ρ)|12Vε42)=12(S2+22πρ2λ1S)S2+22πρ2λ2S+o(ε2)<12(S2+22πρ2λ1S)32(S2+22πρ2λ2S)<0. (3.12)

    Together with (3.10), when ε > 0 and ρ > 0 small enough, we obtain – A1A2 < 0. Then ∀ a > 0, (a1, a2) = (a,λ12λ2a) satisfies (3.9). Therefore (3.5) has a positive solution, and in particular

    λ1a12+λ2a222<0,a1,a2>0.

    So we may choose Q=λ1a12+λ2a222<0 and then (3.7) holds.

  2. For the case λ1, λ2 > 0, then obviously Pε > 0. Therefore, the solution k exists for ∀ a1, a2 > 0. In the following proposition, we will show that if the positive solution k exists, then k has positive lower bound. Note that

    Vε2qp+22qp+2Vε22=B2ρ(0)(Ψ(x)Uε,0(x))2qp+2dxB2ρ(0)(Ψ(x)Uε,0(x))2dx=02ρεr3(1+r2)2qp+2dr+o(1)8qpε2qp02ρεr3(1+r2)2dr+o(1)+,ε0+,

    we can take G=λ1a12+λ2a22ya12qpa22k02qp>0, here k0 will be given by Proposition 3.3. Then we obtain that

    λ1a12+λ2a222Vε2qp+22qp+2Vε22ya12qpa22k2qp<λ1a12+λ2a222Gya1a22k0=λ1a12+λ2a222

    when ε > 0 is small enough, so we may let Q=λ1a12+λ2a222, then Q < 0 and satisfies (3.7).

  3. Recall the definition of Ai in the case (a) and λ1, λ2 < 0, if ρ satisfies

    0<ρ<(S22πmin(λ1,λ2))12,

    then

    Ai=Vε22+λiVε22S2+22πρ2λiS+o(ε2)>12(S2+22πρ2λiS)>0,i=1,2

    for ε > 0 small enough. Similar to the proof of the case (a), we obtain Pε=a12(A1+(a2a1)2A2)>0 for any a1, a2 > 0, so (3.5) has a unique positive solution. Further, if we take Q=λ1a12+λ2a222, then (3.7) holds.□

Proposition 3.3

Assume μ1, μ2, β, y > 0, then for all the three cases of Lemma 3.2, ∃ a1, a2 > 0 such that the positive solution k of (3.5) exists and satisfies:

0<k0<k<1,

for ε > 0 small enough, here k0 is a constant depending on μi, β, y, a1, a2.

Proof

  1. We are going to show that 0 < k < 1.

    1. For the case min(λ1, λ2) < 0 < max (λ1, λ2). According to the proof of (a) in Lemma 3.2, ∀ a > 0, ∃ (a1, a2) = (a,λ12λ2a) such that k exists and λ1a12+λ2a22 < 0. Then

      Gε(1)=((μ1a14+2βa12a22+μ2a24)Vε44)+((2qp+2)ya12qpa22Vε2qp+22qp+2)((a12+a22)Vε22+(λ1a12+λ2a22)Vε22)(μ1a14+2βa12a22+μ2a24)(S2+o(ε4))(a12+a22)(S2+o(ε2))C(λ1a12+λ2a22)(ε2|lnε|+o(ε2))=(μ1a14+2βa12a22+μ2a24(a12+a22))S2+o(ε).

      We can take a > 0 large enough, such that a1 and a2 satisfy

      μ1a14+2βa12a22+μ2a24(a12+a22)>0.

      Therefore, Gε(1) > 0 for ε > 0 small enough, and we observe that ∃ k ∈ (0, 1) such that Gε(k) = 0, here k depends on a1, a2, ε only, therefore it is the unique solution of the algebraic equation (3.5).

    2. For the case λ1, λ2 > 0. It follows from λ1a12+λ2a22 > 0 and

      Vε2|B(0,2ρ)|14Vε4,|B(0,2ρ)|=12π2(2ρ)4=8π2ρ4,

      that

      Gε(1)=((μ1a14+2βa12a22+μ2a24)Vε44)+((2qp+2)ya12qpa22Vε2qp+22qp+2)((a12+a22)Vε22+(λ1a12+λ2a22)Vε22)(μ1a14+2βa12a22+μ2a24)(S2+o(ε2))(a12+a22)(S2+o(ε4))22πρ2(λ1a12+λ2a22)(S+o(ε4))=S[(μ1a14+2βa12a22+μ2a24(a12+a22))S22πρ2(λ1a12+λ2a22)]+o(ε2).

      According to the case (b) of Lemma 3.2, ∀ a1, a2 > 0, the unique solution k of (3.5) exists. So we can take a1, a2 large enough such that

      (μ1a14+2βa12a22+μ2a24(a12+a22))S22πρ2(λ1a12+λ2a22)>0,

      and then Gε(1) > 0 for ε small enough, which means 0 < k < 1 since Gε(0) < 0.

    3. For the case of max(λ1, λ2) < 0, we see that λ1a12+λ2a22 < 0. We can apply the same calculation as that used in i) to show that for a1, a2 large enough and ρ, ε > 0 small enough, there holds 0 < k < 1.

  2. We show that k has a lower bound, i.e., k > k0 > 0.

    1. Assume λ1 < 0 < λ2. As the Step 1, we consider a1, a2 so that λ1a12+λ2a22 < 0. Together with

      Vε2qp+2|B(0,2ρ)|pq4(p+q)Vε4,

      we have

      Gε(k)(μ1a14+2βa12a22+μ2a24)(S2+o(ε4))k2(λ1a12+λ2a22)|B(0,2ρ)|12(S+o(ε4))+(2qp+2)ya12qpa22|B(0,2ρ)|pq2p(Sqp+1+o(ε4))k2qp(a12+a22)(S2+o(ε2))=:G(k)+o(ε),

      here

      G(k)=(μ1a14+2βa12a22+μ2a24)S2k2+(2qp+2)ya12qpa22|B(0,2ρ)|pq2pSqp+1k2qp(a12+a22)S2(λ1a12+λ2a22)|B(0,2ρ)|12S.

      As the proof of Lemma 3.2, we take ρ > 0 such that 𝓢 + 2 2 πρ2λ1 > 0 and (a1, a2) = (a,λ12λ2a) . Then we obtain that – S+22πρ2λ2S+22πρ2λ1 < 0, it follows that (a1a2)2>S+22πρ2λ2S+22πρ2λ1 always holds, which means that

      G(0)=(a12+a22)S2(λ1a12+λ2a22)|B(0,2ρ)|12S<0.

      Together with limk+ G*(k) = +∞, there exists k0 > 0 dependent on μi, β, y, a1, a2 such that G*(k0) < 0. Thus

      Gε(k0)G(k0)+o(ε)<0,

      when ε is small enough and then we obtain k0 < k.

    2. Assume that λi > 0, i = 1, 2. In this case, λ1a12+λ2a22 > 0, it follows that

      Gε(k)(μ1a14+2βa12a22+μ2a24)(S2+o(ε4))k2(λ1a12+λ2a22)(Cε2|lnε|+o(ε2))(a12+a22)(S2+o(ε2))+(2qp+2)ya12qpa22|B(0,2ρ)|pq2p(Sqp+1+o(ε4))k2qp(μ1a14+2βa12a22+μ2a24)S2k2+(8π2ρ4)pq2p(2qp+2)ya12qpa22Sqp+1k2qp(a12+a22)S2+o(ε)=:G(k)+o(ε),

      here

      G(k)=(μ1a14+2βa12a22+μ2a24)S2k2+(8π2ρ4)pq2p(2qp+2)ya12qpa22Sqp+1k2qp(a12+a22)S2.

      Obviously, G*(0) < 0 and limk+ G*(k) = +∞, then ∃ k0 > 0 dependent on μi, β, y, a1, a2 such that G*(k0) < 0. Thus Gε(k0) ≤ G*(k0) + o(ε2) < 0 for ε > 0 small enough. Recall that k is the unique positive solution satisfying Gε(k) = 0, we obtain k0 < k.

    3. For the case max(λ1, λ2) < 0, we have λ1a12+λ2a22 < 0 as that in case 1). According to the proof of Lemma 3.2, we take ρ > 0 satisfying 𝓢 + 2 2 πρ2λi > 0, i = 1, 2. Thus we can take a1, a2 > 0, note that (a1a2)2>S+22πρ2λ2S+22πρ2λ1, then we have

      G(0)=(a12+a22)S2(λ1a12+λ2a22)|B(0,2ρ)|12S<0.

      Therefore, similar to the case 1), we get a constant k0 > 0 independent of ε such that k0 < k.□

Lemma 3.4

Assume μ1, μ2, β, y > 0. If one of the following cases holds:

  1. min(λ1, λ2) < 0 < max(λ1, λ2) and β > max(μ1, μ2);

  2. min(λ1, λ2) > 0 and β > max(μ1, μ2);

  3. max(λ1, λ2) < 0;

then ℬ < c.

Proof

The main point of the proof is constructing a test function (kv1, kv2) ∈ 𝓝.

  1. For the case min(λ1, λ2) < 0 < max(λ1, λ2) and β > max(μ1, μ2), we choose (kv1, kv2) = (ka1 Vε, ka2 Vε). We have to determine the values k, a1, a2. By Lemma 3.2, ∃ a1, a2 > 0 such that for ρ satisfying (3.11) and ε > 0 small enough, there exists a unique k which is only dependent on ε, a1, a2 such that (kv1, kv2) ∈ 𝓝. Moreover, ∃ Q < 0 such that (3.7) holds. Together with Proposition 3.3, k is bounded, which means

    k2(ε)o(ε2)=o(ε2),asε0 (3.13)

    and

    Qk2(ε)<Qk02<0,Q<0. (3.14)

    Then it follows from (3.4) and (3.7) that

    B=infNII(ka1Vε,ka2Vε)=12(a12+a22)k2Vε22+12(λ1a12+λ2a22)k2Vε2214(μ1a14+2βa12a22+μ2a24)k4Vε44ya12qpa22k2qp+2Vε2qp+22qp+212(a12+a22)k2(S2+o(ε2))14(μ1a14+2βa12a22+μ2a24)k4(S2+o(ε4))+Qk2Vε22{(ka1)2+(ka2)2214[μ1(ka1)4+μ2(ka2)4+2β(ka1)2(ka2)2]}S2+Qk2(Cε2|lnε|+o(ε2))+o(ε2). (3.15)

    It is standard (cf.[24] lemma 2.2) to show that f(t1, t2) given by (3.2) has a global maximum point (t1, t2) as (3.3) in (ℝ+)2 when β > max(μ1, μ2). This claims

    {(ka1)2+(ka2)2214[μ1(ka1)4+μ2(ka2)4+2β(ka1)2(ka2)2]}S2dS2=c,

    combining with (3.13)(3.14), it follows that ℬ ≤ c2|ln ε| + o(ε2), which implies that ℬ < c as ε > 0 small enough.

  2. the proof is the same as a) by Lemma 3.2 and Proposition 3.3.

  3. For the case max(λ1, λ2) < 0, we take the test function of ℬ as

    (kv1,kv2)=(kt1Vε,kt2Vε)N,

where t1, t2 are constants given by (3.3) and k = k(ε) satisfies

((t12+t22)Vε44)k2+((2qp+2)yt12qpt22Vε2qp+22qp+2)k2qp((t12+t22)Vε22+(λ1t12+λ2t22)Vε22)=0. (3.16)

According to Lemma 3.2, such a k > 0 exists. Now we define

gε(k):=k22Vε22k44Vε44.

A direct calculation shows that

dgε(k)dk=Vε22kVε44k3,

so the function gε(k) is increasing when 0 < k < Vε2Vε42 and decreasing when k > Vε2Vε42 . Recall that

Vε2Vε42=S+o(ε2)S+o(ε4)=1+o(ε2),

we obtain

gε(k)gε(1+o(ε2))=12Vε2214Vε44+o(ε2).

Then

I(kt1Vε,kt2Vε)=k22Ω(t12+t22)|Vε|2k44Ω(μ1t14+μ2t24+2βt12t22)Vε4+k22Ω(λ1t12+λ2t22)Vε2yk2qp+2Ωt12qpt22Vε2qp+2t12+t222k2Vε22t12+t224k4Vε44+k22(λ1t12+λ2t22)Vε22t12+t222(S2+o(ε2))t12+t224(S2+o(ε4))Cε2|lnε|+o(ε2)=14(t12+t22)S2Cε2|lnε|+o(ε2). (3.17)

On the other hand, note that Uε,022=Uε,044=S2, we have

14(t12+t22)S2=12(t12+t22)Uε,02214(μ1t14+μ2t24+2βt12t22)Uε,044=E(t1Uε,0,t2Uε,0)=c. (3.18)

Then we have ℬ ≤ I(kt1Vε, kt2Vε) ≤ c2|ln ε| + o(ε2). Thus we can obtain ℬ < c by taking ε > 0 small enough.□

In the last section, we will use the Nehari manifold 𝓝 given by (1.6) to obtain a (PS)-sequence {(un, vn)} ⊂ 𝓝 at the least energy level ℬ by the Ekeland’s variational principle, here ℬ is given by (1.8). Indeed, we will show that 𝓝 is a C1-manifold, and it is a natural constraint.

Lemma 3.5

The Nehari manifold 𝓝 of I(u, v) is a natural constraint in 𝓗 in the case β ≥ – μ1μ2 , y > 0.

Proof

Assume (u, v) ∈ 𝓝 such that I|N (u, v) = 0. Then there exists a constant k ∈ ℝ which is a Lagrange multiplier, such that

I(u,v)+kF(u,v)=0,

then

I(u,v)(u,v)+kF(u,v)(u,v)=0.

We have I′(u, v)(u, v) = 0, since (u, v) ∈ 𝓝, then

kF(u,v)(u,v)=0.

On the other hand,

F(u,v)(u,v)=2Ω(|u|2+|v|2+λ1u2+λ2v2)4Ω(μ1u4+μ2v4+2βu2v2)(2qp+2)2yΩu2qpv2=2Ω(μ1u4+μ2v4+2βu2v2)2qp(2qp+2)yΩu2qpv22Ω(μ1u2μ2v2)22qp(2qp+2)yΩu2qpv2<0,

since β ≥ – μ1μ2 , y > 0. Thus we have k = 0, which means I′(u, v) = 0.□

The following is a version of the Ekeland’s variational principle (cf. [22, Theorem 8.5]).

Lemma 3.6

Let X be a Banach space and let GC2(X, ℝ) be such that, for every vV := {vX | G(v) = 1}, G′(v) ≠ 0. Let FC1(X, ℝ) be bounded below on V, vV and ε, δ > 0. If

F(v)infVF+ε,

then there exists uV such that

F(u)infVF+2ε,minλRF(u)λG(u)8εδ,uv2δ.

Remark 3.7

(cf. [22, Definition 5.10]) Let X be a Banach space, ψC2(X, ℝ), V := {vX : ψ(v) = 1} and ψ′(v) ≠ 0 for every vV. The norm of the derivative of the restriction of ψ to V at v is defined byφ′(v)∥*. Moreover, (cf. [22, Proposition 5.12]), if φC1(X, ℝ) and uV then

φ(u)=minλRφ(u)λψ(u).

Lemma 3.8

Assume μ1, μ2, β, y > 0, –α1(Ω) < λ1 < 0, –α1(Ω) < λ2, uμ1 is given by (1.10). Ifv0 H01 (Ω) ∖ {0}, such that β, y satisfy

(Ωuμ12v02)β+(2Ωuμ12qpv02)y>(Ω|v0|2+λ2v02), (3.19)

then ℬ < Bμ1.

Proof

∀ (u, v) ∈ 𝓗, assume that ∃ t > 0 such that (tu, tv) ∈ 𝓝, then t satisfies

(Ωμ1u4+μ2v4+2βu2v2)t2+((2qp+2)yΩu2qpv2)t2qp=(Ω|u|2+|v|2+λ1u2+λ2v2). (3.20)

Let (u, v) = (uμ1, s v0), here uμ1 is given by (1.10), 0 ≠ v0 H01 (Ω) and 0 ≠ s ∈ ℝ. By (3.20), there exists t(s) = t(s, uμ1, v0) such that (t(s)uμ1, t(s)sv0) ∈ 𝓝, and t(s) is the positive solution of the following equation of t:

Hs(t):=(Ωμ1uμ14+μ2s4v04+2βuμ12v02s2)t2+((2qp+2)ys2Ωu2qpv02)t2qp(Ω|uμ1|2+s2|v0|2+λ1uμ12+s2λ2v02)=0. (3.21)

In fact, Hs(t) is continuous and monotonically increasing on [0, +∞) with

limk+Hs(k)=+

and Ω|v0|2+λ2v022Bμ2(Ωv04)12>0, when λ2 < 0; Ω|v0|2+λ2v020, when λ2 ≥ 0. So

Hs(0)=Ω(|uμ1|2+λ1uμ12+s2|v0|2+s2λ2v02)=(4Bμ1+s2Ω|v0|2+λ2v02)<0.

Then ∃ unique t = t(s) > 0 such that (3.21) holds. Moreover, t(s) → 1, when s → 0. In order to simplify the calculation process, we define

m1:=Ωμ2v04,m2:=Ω2βuμ12v02,m3:=(2qp+2)yΩuμ12qpv02,m4:=4Bμ1,m5:=Ω|v0|2+λ2v02.

We get t′(s) from (3.21) by implicit function differentiation

t(s)=2m1s3t2+m2st2+m3st2qpm5sm1s4t+m2s2t+qpm3s2t2qp1+m4t,
t(s)s=2m1s2t2+m2t2+m3t2qpm5m1s4t+m2s2t+qpm3s2t2qp1+m4tm2+m3m5m4=(Ω|v0|2+λ2v02)(2qp+2)y(Ωuμ12qpv02)2βΩuμ12v024Bμ1,as s0. (3.22)

Let

k:=(Ω|v0|2+λ2v02)(2qp+2)y(Ωuμ12qpv02)2βΩuμ12v024Bμ1. (3.23)

Then

t(s)=ks+o(s),t(s)=1+k2s2+o(s2),t2(s)=1+ks2+o(s2),t2qp+2(s)=1+(qp+1)ks2+o(s2).

Therefore we have

I(tuμ1,tsv0)=t24Ω(|uμ1|2+λ1uμ12+s2(|v0|2+λ2v02))(pq2p)s2t2qp+2yΩuμ12qpv02=t2Bμ1+14t2s2(Ω|v0|2+λ2v02)(pq2p)s2t2qp+2yΩuμ12qpv02=(1+ks2+o(s2))Bμ1+s24(1+ks2+o(s2))(Ω|v0|2+λ2v02)(pq2p)s2(1+(qp+1)ks2+o(s2))yΩuμ12qpv02=Bμ1+(kBμ1+14Ω(|v0|2+λ2v02)(pq2p)yΩuμ12qpv02)s2+o(s2)<Bμ1, (3.24)

as |s| > 0 small enough, here (3.19) and (3.23) are applied. Then

BI(t(s)uμ1,t(s)sv0)<Bμ1.

This completes the proof.□

Lemma 3.9

Assume μ1, μ2, y, β > 0. Ifα1(Ω) < λ2 < 0, then ℬ < Bμ2.

Proof

According to the proof of Lemma 3.8, we take the test function as (tsu, tuμ2), here 0 < s, t(s) ∈ ℝ, 0 ≠ u H01 (Ω) and uμ2 is given by (1.10). In order to ensure (tsu, tuμ2) ∈ 𝓝, we have

Ls(t):=(Ωμ1s4u4+μ2uμ24+2βu2uμ22s2)t2+((2qp+2)ys2qpΩu2qpuμ22)t2qp(Ω|uμ2|2+s2|u|2+λ2uμ22+s2λ1u2)=0. (3.25)

Then as the same as Lemma 3.8, ∃ unique t = t(s) > 0 such that (3.25) holds and t(s) → 1, when s → 0. We define

w1:=Ωμ1u4,w2:=Ω2βu2uμ22,w3:=(2qp+2)yΩu2qpuμ22,w4:=4Bμ2,w5:=Ω|u|2+λ1u2.

By the implicit function differentiation, we have

t(s)=2w1s3t2+w2st2+qpw3s2qp1t2qpw5sw1s4t+w2s2t+qpw3s2qpt2qp1+w4t,
t(s)s2qp1=2w1s42qpt2+w2s22qpt2+qpw3t2qpw5s22qpw1s4t+w2s2t+qpw3s2qpt2qp1+w4tqw3pw4=qp(2qp+2)yΩu2qpuμ224Bμ2,as s0. (3.26)

Let

d:=qp(2qp+2)yΩu2qpuμ224Bμ2. (3.27)

Then

t(s)=ds2qp1+o(s2qp1),t(s)=1+p2qds2qp+o(s2qp),t2(s)=1+pqds2qp+o(s2qp),t2qp+2(s)=1+(qp+1)pqds2qp+o(s2qp).

Therefore we have

I(tsu,tuμ2)=t24Ω(s2(|u|2+λ1u2)+(|uμ2|2+λ2uμ22))(pq2p)s2qpt2qp+2yΩu2qpuμ22=t2Bμ2+14t2s2(Ω|u|2+λ1u2)(pq2p)s2qpt2qp+2yΩu2qpuμ22=(1+pqds2qp+o(s2qp))Bμ2+s24(1+pqds2qp+o(s2qp))(Ω|u|2+λ1u2)(pq2p)s2qp(1+(qp+1)pdqs2qp+o(s2qp))yΩu2qpuμ22=Bμ2+(pdqBμ2(pq2p)yΩu2qpuμ22)s2qp+o(s2qp)=Bμ2(yΩu2qpuμ22)s2qp+o(s2qp)<Bμ2, (3.28)

as |s| > 0 small enough, here y > 0 and (3.7) are applied. Then

BI(t(s)su,t(s)uμ2)<Bμ2.

This completes the proof.□

Lemma 3.10

Assume β, y, μ1, μ2 > 0. λ1, λ2 ∈ ℝ, λi > –α1(Ω), i = 1, 2. Suppose that λ1 satisfies one of the following two conditions:

  1. λ1 < 0 and β, y satisfy (3.19);

  2. λ1 ≥ 0 with Ω is a (smooth) starshaped domain,

and that λ2 satisfies one of the following two conditions:

  1. λ2 < 0;

  2. λ2 ≥ 0 with Ω is a (smooth) starshaped domain.

If (u1, v1) is the minimizer of I(u, v) on 𝓝, then (u2, v2) = (|u1|, |v1|) is a ground state solution of system (1.1): u2 > 0, v2 > 0.

Proof

Since ∀ u H01 (Ω), |∇|u|| ≤ |∇u| a.e. ∈ ℝ4, we obtain that (u2, v2) ∈ 𝓝 and I(|u1|, |v1|) ≤ I(u1, v1) = minN I. Thus (u2, v2) is also a minimizer of I(u, v) on 𝓝. Since 𝓝 is a natural constraint by Lemma 3.5, (u2, v2) satisfies the system

Δu+λ1u=μ1u3+βuv2+2qpyu2qp1v2inΩ,Δv+λ2v=μ2v3+βu2v+2yu2qpvinΩ,u,v0inΩ,u2,v2=0onΩ.

By the maximum principle and 0 ≠ (u2, v2) ∈ 𝓝, there exists the following three possibilities: u2 = 0, v2 > 0; u2 > 0, v2 > 0 or u2 > 0, v2 = 0. Suppose that (u2, v2) is not a solution of (1.1), then u2 > 0, v2 = 0 or u2 = 0, v2 > 0. First, we suppose u2 > 0, v2 = 0 if λ1 < 0, which leads to a contradiction since B=infNI(u,v)<Bμ1=infNμ1I(u,0) by Lemma 3.8; if λ1 ≥ 0, then u2 satisfies the equation following:

Δu2+λ1u2=μ1u23inΩ,u2>0inΩ,u2=0onΩ,

which leads to a contradiction for that Brezis-Nirenberg equation has no solution when λ1 ≥ 0. Suppose u2 = 0, v2 > 0 and prove as the last assumption by Lemma 3.9, we obtain a contradiction. Thus u2 > 0, v2 > 0 and (u2, v2) is a ground state of system (1.1).□

4 Proof of Theorem 1.4

Define the function

z(x,y):=2qpq(λ1x+λ2y)+μ1x2+μ2y2+2βxy,x,yR. (4.1)

Lemma 4.1

If μ1, μ2 > 0, β > μ1μ2 or β < – μ1μ2 , then z(x, y) ≥ 0.

Proof

Let

zx=2qpqλ1+2μ1x+2βy=0zy=2qpqλ2+2μ2y+2βx=0,

we have critical point (x0, y0) as

(x0,y0)=qpq(λ1μ2βλ2β2μ1μ2,λ2μ1βλ1β2μ1μ2).

Then, 2zx2|(x0,y0)=2μ1,2zy2|(x0,y0)=2μ2,2zxy|(x0,y0)=2β. Since (2β)2 – (2μ1)(2μ2) > 0, z(x, y)min = z(x0, y0) = (qpq)2(λ12μ2+λ22μ12λ1λ2ββ2μ1μ2) ≥ 0. The last inequality holds since the function h(x, y) := μ2x2 + μ1y2 – 2βxy has minimum 0 at (0, 0) by the same calculation as above.□

Lemma 4.2

Assume that β, y > 0 and λ1, λ2 > –α1(Ω), where α1(Ω) > 0 is the first eigenvalue ofΔ in H01 (Ω). Then there exists C > 0 such thatu, v ∈ 𝓝,

u22+v22C.

Proof

By the Hölder and Young inequalities, ∀ ε > 0, ∃ Cε > 0, such that

Ωu4qp(Ωu2)22qp(Ωu4)2qp1ε(Ωu2)+Cε(Ωu4).

For u, v ∈ 𝓝, we have

Ω|u|2+|v|2=Ωμ1u4+μ2v4+2βu2v2+(2qp+2)yu2qpv2λ1u2λ2v2Ω(μ1+β)u4+(μ2+β+(qp+1)y)v4+(qp+1)yu4qpλ1u2λ2v2(μ1+β+Cε(qp+1)y)Ωu4+((qp+1)yελ2)Ωv2+(μ2+β+(qp+1)y)Ωv4λ1Ωu2,

that is,

Ω|u|2+λ1u2+|v|2+(λ2(qp+1)yε)Ωv2CΩu4+v4C(Ω|u|2+|v|2)2. (4.2)

By the Poíncaré inequality

fL2(Ω)1α1(Ω)fL2(Ω),fH01(Ω)

we obtain

Ω|u|2+λ1u2C1Ω|u|2,

and

Ω|v|2+(λ2(qp+1)yε)v2C2Ω|u|2,

for ε small enough, here C1, C2 are positive constants. Together with (4.2), we obtain u22+v22C.

Proof of Theorem 1.4

First, we verify the assumptions of Ekeland’s variational principle, i.e., I|𝓝(u, v) is bounded below on 𝓝.

I|N(u,v)=q2(p+q)Ω(|u|2+|v|2+λ1u2+λ2v2)+pq4(p+q)Ω(μ1u4+μ2v4+2βu2v2)=pq4(p+q)Ω(2qpq(|u|2+|v|2)+z(u2,v2)),

where z(x, y) is defined in (4.1), we obtain that I|𝓝(un, vn) ≥ 0. Then according to the Ekeland’s variational principle, there is a (PS) sequence {(un, vn)} ⊂ 𝓝 of I(u, v) at the least energy level ℬ. Since ℬ = limn0 I|𝓝(un, vn) ≥ 0 and z(u2, v2) ≥ 0, we have

q2(p+q)Ω(|un|2+|vn|2)I|N(un,vn)M.

So {(un, vn)} is bounded in 𝓗, we may assume that (un, vn) ⇀ (u, v) weakly in 𝓗. Next, we will prove that (u, v) ≠ (0, 0). Suppose to the contrary that (u, v) = (0, 0), then up to a subsequence

un0,vn0L2(Ω),un0,vn0L2qp+2(Ω). (4.3)

Therefore, we have

|un|2qp+20,|vn|2qp+20L1(Ω),

as n → ∞ and

Ω|un2qpvn2|(Ω|un|2qp+2)qp+q(Ω|vn|2qp+2)pp+q=o(1), (4.4)
un22+vn22=μ1un44+μ2vn44+2βΩun2vn2+o(1). (4.5)

Take tn as following such that (tnun, tnvn) ∈ 𝓜 and with (4.5), we have

tn=un22+vn22μ1un44+μ2vn44+2βΩun2vn2=un22+vn22(un22+vn22)+o(1)=1+o(1), (4.6)

since un22+vn22C>0 by Lemma 4.2. On the other hand, by (4.3)-(4.4)

I(un,vn)=E(un,vn)+12Ωλ1un2+λ2vn22yun2qpvn2=o(1),

here E(u, v) is given in (3.1). Therefore, together with (4.6) we have

B+o(1)=I(un,vn)=E(tnun,tnvn)+o(1)c+o(1).

Thus we obtain that ℬ ≥ c, which contradicts Lemma 3.4. Thus (u, v) ≠ (0, 0).

Since I|N (un, vn) = o(1) and according to Lemma 3.6 and Remark 3.7, we have

I(un,vn)=o(1). (4.7)

As (un, vn) ⇀ (u, v) weakly in 𝓗, passing to a subsequence, we assume that

unu,vnv, weakly in L4(Ω),unu,vnv,a.e.xΩ,unu,vnv,strongly inL2(Ω).

Thus by (4.7), we have I′(u, v) = 0. Together with (u, v) ≠ (0, 0), we obtain (u, v) ∈ 𝓝. Then by the Fatou Lemma, we have

B=limnI(un,vn)I(u,v)infNI=B.

Thus we obtain that (u, v) is the minimizer of I(u, v) in 𝓝. Then by Lemma 3.10, system (1.1) has a positive ground state solution: u > 0 and v > 0.□

Acknowledgements

The authors wish to thank the referees very much for his/her careful reading the manuscript and valuable comments and suggestions.

This work is supported by NSFC(11771234, 11926323).

References

[1] N. Akhmediev, A. Ankiewicz, Novel soliton states and bifurcation phenomena in nonlinear fiber couplers, Phys. Rev. Lett., 70, 2395-2398, 1993.10.1103/PhysRevLett.70.2395Suche in Google Scholar PubMed

[2] A. Ambrosetti, E. Colorado, D. Ruiz, Multi-bump solitons to linearly coupled systems of nonlinear Schrödinger equations, Calc. Var. Partial Differential Equations, 30, 85-112, 2007.10.1007/s00526-006-0079-0Suche in Google Scholar

[3] A. Ambrosetti, Remarks on some systems of nonlinear Schrödinger equations, Fixed Point Theory Appl., 4, 35-46, 2008.10.1007/s11784-007-0035-4Suche in Google Scholar

[4] A. Ambrosetti, G. Cerami, D. Ruiz, Solitons of linearly coupled systems of semilinear non-autonomous equations on ℕ RN, J. Funct. Anal., 254, 2816-2845, 2008.10.1016/j.jfa.2007.11.013Suche in Google Scholar

[5] T. Bartsch and N. Soave, A natural constraint approach to normalized solutions of nonlinear Schrödinger equations and systems, Arch.Ration.Mech.Anal., 205, no. 2, 515-551, 2012.10.1016/j.jfa.2017.01.025Suche in Google Scholar

[6] H. Brezis, E.H. Lieb, Minimum action solutions of some vector field equations, Comm. Math. Phys., 96, 97-C113, 1984.10.1007/978-3-642-55925-9_45Suche in Google Scholar

[7] H. Brezis and L. Nirenberg, Positive solutions of nonlinear elliptic equations involving critical Sobolev exponents, Commun.Pure Appl.Math., 36, 437-477, 1983.10.1002/cpa.3160360405Suche in Google Scholar

[8] Z. Chen, C. S. Lin, W. Zou, Infinitely many sign-changing and semi-nodal solutions for a nonlinear Schrödinger system. Ann. Sc. Norm. Super. Pisa Cl. Sci., (5) 15(2016), 859-897.10.2422/2036-2145.201401_002Suche in Google Scholar

[9] Z. Chen and W. Zou, Positive least energy solutions and phase separation for coupled Schrd̈inger equations with critical exponent: higher dimensional case, Calc. Var. Partial Differential Equations, 52, 423-467, 2015.10.1007/s00526-014-0717-xSuche in Google Scholar

[10] Z. Chen and W. Zou, Positive least energy solutions and phase separation for coupled Schrödinger equations with critical exponent, Arch. Ration. Mech. Anal., 205, 515-551, 2012.10.1007/s00205-012-0513-8Suche in Google Scholar

[11] E. Dancer and J. Wei, Spike solutions in coupled nonlinear Schrödinger equations with attractive interaction, Trans. Ammer. Math. Soc., 361, 1189-1208, 2009.10.1090/S0002-9947-08-04735-1Suche in Google Scholar

[12] D. Gilbarg and N. Trudinger, Elliptic Partial Differential Equations of Second Order, 2nd edn. Springer, Berlin, 1998.Suche in Google Scholar

[13] D. Hall, M. Matthews, J. Ensher, C. Wieman and E. Cornell, Dynamics of component separation in a binary mixture of Bose-Einstein condensates, Phys. Rev. Lett., 81, 1539-1542, 1998.10.1142/9789812813787_0071Suche in Google Scholar

[14] T. C. Lin and J. Wei, Ground state of N coupled nonlinear Schrödinger equations in ℝn, n ≤ 3, Commun. PDEs, 18, 1043-1054, 2005.Suche in Google Scholar

[15] S. Peng, S. Wei and Q. Wang, Multiple positive solutions for linearly coupled nonlinear elliptic systems with critical exponent, J. Differential Equations., 263, no.1, 709-731, 2017.10.1016/j.jde.2017.02.053Suche in Google Scholar

[16] B. Sirakov, Least energy solitary waves for a system of nonlinear Schrödinger equations in ℝN, Commun.Math.Phys., 271, 199-221, 2009.10.1007/s00220-006-0179-xSuche in Google Scholar

[17] M. Struwe, Variational Methods.Applications to Nonlinear Partial Differential Equations and Hamiltonian Systems, 2nd ed., Ergeb. Math. Grenzgeb.(3)34, Springer, Berlin, 1996.Suche in Google Scholar

[18] A. Szulkin and T. Weth, Ground state solutions for some indefinite variational problems, J. Funct. Anal., 257(12), 3802-3822, 2009.10.1016/j.jfa.2009.09.013Suche in Google Scholar

[19] S. Tavares and S. Terracini, Sign-changing solutions of competition-diffusion elliptic systems and optimal partition problem, Ann.Inst.H.Poincaré Anal.Non Linéaire, 92, 279-300, 2012.10.1016/j.anihpc.2011.10.006Suche in Google Scholar

[20] S. Terracini and G. Verzini, Multipulse phases in k-mixtures of Bose-Einstein condensates, Arch.Ration.Mech.Anal., 194, 717-741, 2009.10.1007/s00205-008-0172-ySuche in Google Scholar

[21] J. Wei and T. Weth, Radial solutions and phase separation in a system of two coupled Schrödinger equations, Arch.Ration.Mech.Anal., 190(3), 83-106, 2008.10.1007/s00205-008-0121-9Suche in Google Scholar

[22] M. Willem, Minimax Theorems. Birkhäuser, Boston, 1996.10.1007/978-1-4612-4146-1Suche in Google Scholar

[23] Y. Wu, Ground states of a K-component critical system with linear and nonlinear couplings: the attractive case, Adv. Nonlinear Stud, 19(3), 595-623, 2019.10.1515/ans-2019-2049Suche in Google Scholar

[24] Yu. Wu, On a K-component elliptic system with the Sobolev critical exponent in high dimensions: the repulsive case, Calc. Var. Partial Differential Equations, 56, Art. 151, 51 pp., 2017.10.1007/s00526-017-1252-3Suche in Google Scholar

[25] Y. Wu, T. F. Wu and W. Zou, On a two-component Bose-Einstein condensate with steep potential wells, ANN.Mat., 196, 1695-1737, 2017.10.1007/s10231-017-0635-6Suche in Google Scholar

[26] Y. Wu and W. Zou, Spikes of the two-component elliptic system in ℝ4 with the critical Sobolev exponent. Calc. Var. Partial Differential Equations, 58, Art.24, 44 pp., 2019.10.1007/s00526-018-1479-7Suche in Google Scholar

Received: 2020-05-17
Accepted: 2020-07-30
Published Online: 2020-09-13

© 2021 Xin Zhao and Wenming Zou, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Artikel in diesem Heft

  1. Editorial
  2. Editorial to Volume 10 of ANA
  3. Regular Articles
  4. Convergence Results for Elliptic Variational-Hemivariational Inequalities
  5. Weak and stationary solutions to a Cahn–Hilliard–Brinkman model with singular potentials and source terms
  6. Single peaked traveling wave solutions to a generalized μ-Novikov Equation
  7. Constant sign and nodal solutions for superlinear (p, q)–equations with indefinite potential and a concave boundary term
  8. On isolated singularities of Kirchhoff equations
  9. On the existence of periodic oscillations for pendulum-type equations
  10. Multiplicity of concentrating solutions for a class of magnetic Schrödinger-Poisson type equation
  11. Nehari-type ground state solutions for a Choquard equation with doubly critical exponents
  12. Gradient estimate of a variable power for nonlinear elliptic equations with Orlicz growth
  13. The structure of 𝓐-free measures revisited
  14. Solvability of an infinite system of integral equations on the real half-axis
  15. Positive Solutions for Resonant (p, q)-equations with convection
  16. Concentration behavior of semiclassical solutions for Hamiltonian elliptic system
  17. Global existence and finite time blowup for a nonlocal semilinear pseudo-parabolic equation
  18. On variational nonlinear equations with monotone operators
  19. Existence results for nonlinear degenerate elliptic equations with lower order terms
  20. Blow-up criteria and instability of normalized standing waves for the fractional Schrödinger-Choquard equation
  21. Ground states and multiple solutions for Hamiltonian elliptic system with gradient term
  22. Positivity of solutions to the Cauchy problem for linear and semilinear biharmonic heat equations
  23. Convex solutions of Monge-Ampère equations and systems: Existence, uniqueness and asymptotic behavior
  24. Multiple solutions for critical Choquard-Kirchhoff type equations
  25. Regularity for sub-elliptic systems with VMO-coefficients in the Heisenberg group: the sub-quadratic structure case
  26. Inducing strong convergence of trajectories in dynamical systems associated to monotone inclusions with composite structure
  27. A posteriori analysis of the spectral element discretization of a non linear heat equation
  28. Liouville property of fractional Lane-Emden equation in general unbounded domain
  29. Large data existence theory for three-dimensional unsteady flows of rate-type viscoelastic fluids with stress diffusion
  30. On some classes of generalized Schrödinger equations
  31. Variational formulations of steady rotational equatorial waves
  32. On a class of critical elliptic systems in ℝ4
  33. Exponential stability of the nonlinear Schrödinger equation with locally distributed damping on compact Riemannian manifold
  34. On a degenerate hyperbolic problem for the 3-D steady full Euler equations with axial-symmetry
  35. Existence, multiplicity and nonexistence results for Kirchhoff type equations
  36. Combined effects of Choquard and singular nonlinearities in fractional Kirchhoff problems
  37. Convergence analysis for double phase obstacle problems with multivalued convection term
  38. Multiple solutions for weighted Kirchhoff equations involving critical Hardy-Sobolev exponent
  39. Boundary value problems associated with singular strongly nonlinear equations with functional terms
  40. Global solvability in a three-dimensional Keller-Segel-Stokes system involving arbitrary superlinear logistic degradation
  41. Multiplicity and concentration behaviour of solutions for a fractional Choquard equation with critical growth
  42. Concentration results for a magnetic Schrödinger-Poisson system with critical growth
  43. Periodic solutions for a differential inclusion problem involving the p(t)-Laplacian
  44. The concentration-compactness principles for Ws,p(·,·)(ℝN) and application
  45. Regularity for commutators of the local multilinear fractional maximal operators
  46. An improvement to the John-Nirenberg inequality for functions in critical Sobolev spaces
  47. Local versus nonlocal elliptic equations: short-long range field interactions
  48. Analysis of a diffusive host-pathogen model with standard incidence and distinct dispersal rates
  49. Blowing-up solutions of the time-fractional dispersive equations
  50. Fixed point of some Markov operator of Frobenius-Perron type generated by a random family of point-transformations in ℝd
  51. Non-stationary Navier–Stokes equations in 2D power cusp domain
  52. Non-stationary Navier–Stokes equations in 2D power cusp domain
  53. Nontrivial solutions to the p-harmonic equation with nonlinearity asymptotic to |t|p–2t at infinity
  54. Iterative methods for monotone nonexpansive mappings in uniformly convex spaces
  55. Optimality of Serrin type extension criteria to the Navier-Stokes equations
  56. Fractional Hardy-Sobolev equations with nonhomogeneous terms
  57. New class of sixth-order nonhomogeneous p(x)-Kirchhoff problems with sign-changing weight functions
  58. On the set of positive solutions for resonant Robin (p, q)-equations
  59. Solving Composite Fixed Point Problems with Block Updates
  60. Lions-type theorem of the p-Laplacian and applications
  61. Half-space Gaussian symmetrization: applications to semilinear elliptic problems
  62. Positive radial symmetric solutions for a class of elliptic problems with critical exponent and -1 growth
  63. Global well-posedness of the full compressible Hall-MHD equations
  64. Σ-Shaped Bifurcation Curves
  65. On the critical behavior for inhomogeneous wave inequalities with Hardy potential in an exterior domain
  66. On singular quasilinear elliptic equations with data measures
  67. On the sub–diffusion fractional initial value problem with time variable order
  68. Partial regularity of stable solutions to the fractional Geľfand-Liouville equation
  69. Ground state solutions for a class of fractional Schrodinger-Poisson system with critical growth and vanishing potentials
  70. Initial boundary value problems for the three-dimensional compressible elastic Navier-Stokes-Poisson equations
Heruntergeladen am 8.9.2025 von https://www.degruyterbrill.com/document/doi/10.1515/anona-2020-0136/html?lang=de
Button zum nach oben scrollen