Skip to main content

Multiple solutions to nonlinear Schrödinger equations with critical growth

Abstract

In 2000, Cingolani and Lazzo (J. Differ. Equ. 160:118-138, 2000) studied nonlinear Schrödinger equations with competing potential functions and considered only the subcritical growth. They related the number of solutions with the topology of the global minima set of a suitable ground energy function. In the present paper, we establish these results in the critical case. In particular, we remove the condition c 0 < c , which is a key condition in their paper. In the proofs we apply variational methods and Ljusternik-Schnirelmann theory.

MSC:35J60, 35Q55.

1 Introduction and main result

We investigate the following nonlinear Schrödinger equation:

iħ ψ t = ħ 2 2 m Δψ+W(x)ψg(x,ψ),
(1.1)

which arises in quantum mechanics and provides a description of the dynamics of the particle in a non-relativistic setting. ħ is the Planck’s constant, m>0 denotes the mass of the particle, W: R N R is the electric potential, g is the nonlinear coupling, and ψ is the wave function representing the state of the particle. A standing wave solution of equation (1.1) is a solution of the form ψ(x,t)=u(x) e i E t h . It is clear that ψ(x,t) solves (1.1) if and only if u(x) solves the following stationary equation:

ħ 2 2 m Δu+ ( W ( x ) E ) u=g(x,u).
(1.2)

For simplicity and without loss of generality, we set ε=ħ, V(x)=2m(W(x)E) and g ˜ =2mg, then equation (1.2) is equivalent to

ε 2 Δu+V(x)u= g ˜ (x,u).
(1.3)

A considerable amount of work has been devoted to investigating solutions of (1.3). The existence, multiplicity and qualitative property of such solutions have been extensively studied. For single interior spikes solutions in the whole space R N , please see [19]etc. For multiple interior spikes, please see [10, 11]etc. For single boundary spike solutions with Neumann boundary condition, please see [6, 1215]etc. For multiple boundary spikes, please see [1618]etc. In particular, Wang and Zeng [9] studied the existence and concentration behavior of solutions for NLS with competing potential functions. Cingolani and Lazzo in [19] obtained the multiple solutions for the similar equation. In those papers only the subcritical growth was considered. In the present paper, we complete these studies by considering a class of nonlinearities with the critical growth. In particular, we remove the condition c 0 < c , which is a key condition in [19].

In the sequel, we restrict ourselves to the critical case in which g ˜ (x,u)= | u | 2 2 u+f(x,u). More specifically, we study the following problem:

{ ε 2 Δ u + V ( x ) u = | u | 2 2 u + f ( x , u ) in  R N , u > 0 in  R N , lim | x | u ( x ) = 0 ,
(1.4)

where 2 :=2N/(N2) if N3, and 2 := if N=1,2. fC( R N ×R,R) satisfies

(f1) f(x,t)=0 for each t0;

(f2) lim t 0 + f ( x , t ) t =0;

(f3) there exists q(2, 2 ) such that

lim sup t f ( x , t ) t q 1 <;

(f4) there exists 2<θ< 2 such that

0<θF(x,t)f(x,t)t,t>0,

where F(t):= 0 t f(τ)dτ;

(f5) the function f ( x , t ) t is strictly increasing in t0 for any xΩ.

Our main results are the following theorem.

Theorem 1.1 Let N4. Suppose that f satisfies (f1)-(f5), V is a continuous function in R N and satisfies inf x R N V(x)>0. Then when ε is sufficiently small, the problem (1.4) has at least cat(Σ, Σ δ ) distinct nontrivial solutions.

Here cat(Σ, Σ δ ) denotes the Ljusternik-Schnirelmann category of Σ in Σ δ . By definition (e.g., [20]), the category of A with respect to M, denoted by cat(A,M), is the least integer k such that A A 1 A k , with A i (i=1,,k) closed and contractible in M. We set cat(,M)=0 and cat(A,M)=+ if there are no integers with the above property. We will use the notation cat(M) for cat(M,M).

To prove Theorem 1.1, we mainly use the idea of [15, 19, 21]. More precisely, we can show that the ( PS ) c -condition holds in the subset N ˜ ε (see (4.6)). Hence the standard Ljusternik-Schnirelmann category theory can be applied in N ˜ ε to yield the existence of at least cat( N ˜ ε ) critical points of I ε . And then we construct two continuous mappings

ϕ ε :Σ N ˜ ε
(1.5)

and

β: N ˜ ε Σ δ ,
(1.6)

where

Σ δ = { x R N : dist ( x , Σ ) δ } ,δ>0.
(1.7)

Then a topological argument asserts that

cat( N ˜ ε )2cat(Σ, Σ δ ).

We will also prove that if u is a critical point of I ε satisfying I ε (u) ε N ( c 0 +h(ε)), then u cannot change sign. Hence we obtain at least cat(Σ, Σ δ ) nontrivial critical points of I ε .

The paper is organized as follows. In Section 2, we collect some notations and preliminaries. A compactness result is given in Section 3, which is a key step in our proof. Finally, in Section 4, we prove Theorem 1.1.

2 Notations and preliminaries

H 1 ( R N ) is the usual Sobolev space of real-valued functions defined by

H 1 ( R N ) := { u : u L 2 ( R N )  and  u L 2 ( R N ) }

with the normal

u 2 := R N ( | u | 2 + V ( a ) u 2 ) dx.

Let H ε be the subspace of a Hilbert space H 1 ( R N ) with respect to the norm

u ε 2 := R N ( ε 2 | u | 2 + V ( x ) u 2 ) dx<.

We denote by S the Sobolev constant for the embedding D 1 , 2 ( R N ) L 2 ( R N ), namely

S= inf 0 u D 1 , 2 R N | u | 2 d x ( R N | u | 2 d x ) 2 / 2 ,
(2.1)

where D 1 , 2 ( R N ) is the usual Sobolev space of real-valued functions defined by

D 1 , 2 ( R N ) := { u : u L 2 ( R N )  and  u L 2 ( R N ) } .

We say that a function u H ε is a weak solution of the problem (1.4) if

R N ( ε 2 u v + V ( x ) u v | u | 2 2 u v f ( x , u ) v ) dx=0,v H ε .

In view of (f2) and (f3), we have that the associated functional I ε : H ε R given by

I ε (u)= 1 2 R N ( ε 2 | u | 2 + V ( x ) | u | 2 ) dx 1 2 R N | u | 2 dx R N F(x,u)dx

is well defined. Moreover, I ε C 1 ( H ε ) with the following derivative:

I ε ( u ) , v = R N ( ε 2 u v + V ( x ) u v | u | 2 2 u v f ( x , u ) v ) dx.

Hence, the weak solutions of (1.4) are exactly the critical points of I ε .

Let us recall some known facts about the limiting problem, namely the problem

Δu+V(a)u= | u | 2 2 u+f(a,u)in  R N ,
(2.2)

here a R N acts as a parameter instead of an independent variable. Solutions of (2.2) will be sought in the Sobolev space H 1 ( R N ) as critical points of the functional

J a (u)= 1 2 R N ( | u | 2 + V ( a ) | u | 2 ) dx 1 2 R N | u | 2 dx R N F(a,u)dx.

The least positive critical value G(a) can be characterized as

G(a):= inf u M a J a (u),

where

M a := { v H 1 ( R N ) { 0 } : J a ( u ) , u = 0 } .
(2.3)

An associated critical point w actually solves equation (2.2) and is called a ground state solution or the least energy solution, i.e., w satisfies

J a (w)= inf u M a J a (u).

Moreover, there exist C>0 and δ>0 such that

w(x) e δ | x | ,w(x) e δ | x | for all x R N .
(2.4)

For more details, please see [22, 23].

Set

c 0 := inf a R N G(a),Σ:= { a R N : G ( a ) = c 0 } .

For any δ>0, we denote Σ δ ={x R N :dist(x,Σ)δ}. We need to estimate the super bound of c 0 . In order to do this, we estimate G(a). We shall use a family of radial function defined by

U ε (x)= ( N ( N 2 ) ) N 2 4 ( ε ε 2 + | x | 2 ) N 2 2 .

It is known [20] that

Δ U ε = U ε ( N + 2 ) / ( N 2 ) in  R N .

Moreover, we have

R N | U ε | 2 dx= R N | U ε | 2 dx= S N / 2 .

Set u ε (x)=ϕ(x) U ε (x), where ϕ C 1 is a cut-off function satisfying ϕ(x)1 if |x|δ/2, ϕ(x)0 if |x|δ and 0ϕ(x)1. After a detailed calculation, we have the following estimates:

R N | ( ϕ U ε ) | 2 dx= S N / 2 +O ( ε N 2 ) ,
(2.5)
R N | ϕ U ε | 2 dx= S N / 2 +O ( ε N ) ,
(2.6)
R N | ϕ U ε | 2 dx=α(ε):={ C ε 2 | log ε | + O ( ε 2 ) if  N = 4 , C ε 2 + O ( ε N 2 ) if  N 5 .
(2.7)

Since F(t)0, from (2.5)-(2.7), we conclude

G ( a ) max t > 0 J a ( t u ε ) max t > 0 ( t 2 2 R N ( | u ε | 2 + V ( a ) | u ε | 2 ) d x t 2 2 R N | u ε | 2 d x ) ,
(2.8)

then the maximum value of the right-hand side is achieved at

τ= ( R N ( | u ε | 2 + V ( a ) | u ε | 2 ) d x R N | u ε | 2 d x ) N 2 4
(2.9)

and

max t > 0 J a ( t u ε ) 1 N ( R N ( | u ε | 2 + V ( a ) | u ε | 2 ) d x ) N / 2 ( R N | u ε | 2 d x ) ( N 2 ) / 2 = 1 N ( R N ( | u ε | 2 + V ( a ) | u ε | 2 ) d x ( R N | u ε | 2 d x ) 2 / 2 ) N / 2 1 N ( S N / 2 + V ( a ) α ( ε ) + O ( ε N 2 ) S ( N 2 ) / 2 + O ( ε N 2 ) ) N / 2 < 1 N S N / 2 .
(2.10)

Hence we have

c 0 < 1 N S N / 2 .
(2.11)

We denote the Nehari manifold of I ε by

N ε = { u H ε , A { 0 } : I ε ( u ) , u = 0 } .

3 Compactness result

Proposition 3.1 Let ε n 0 as n. Assume that ( u n ) N ε n satisfies ε n N I ε n ( u n ) c 0 as n. Then uniformly in aΣ, there exist a subsequence of v n (y):= u n ( ε n y+a) (still denoted by v n ), and t n >0 such that w n := t n v n M a . Furthermore, w n converges strongly in H 1 ( R N ) to w, the positive ground state solution of equation (2.2).

Proof Let ( u n ) N ε n be such that ε n N I ε n ( u n ) c 0 . Then, by a change of variable x= ε n y+a, we have

c 0 + 1 ε n N ( I ε n ( u n ) 1 θ I ε n ( u n ) , u n ) ε n N ( 1 2 1 θ ) u n ε n 2 = ( 1 2 1 θ ) R N ( | u n ( ε n y + a ) | 2 + V ( ε n y + a ) | u n ( ε n y + a ) | 2 ) d y .
(3.1)

This implies that ( u n ) is bounded in H 1 ( R N ). Noting that

c 0 + o ( 1 ) = 1 2 R N ( | u n + V ( ε n y + a ) | u n | 2 ) d y 1 2 R N | u n | 2 d y R N F ( x , u n ) d y ,
(3.2)

hence

m:= lim n R N ( | u n + V ( ε n y + a ) | u n | 2 ) dy>0,
(3.3)

since c 0 >0. Now we prove that there exists a sequence ( z n ) R N and constants R,γ>0 such that

lim inf n B R ( z n ) | u n | 2 dyγ>0.
(3.4)

Indeed, if this is not true, then the boundedness of ( u n ) in H 1 ( R N ) and a lemma due to Lions [[24], Lemma I.1] imply that u n 0 in L s ( R N ) for all 2<s< 2 . Given δ>0, we can use (f2), (f3) and u n N ε n to get

R N f(x, u n ) u n dyδ R N | u n | 2 dy+ C δ R N | u n | q dy.

Moreover,

R N ( | u n + V ( ε n y + a ) | u n | 2 ) dy= R N | u n | 2 dy

as n. Therefore

m= lim n R N | u n | 2 dy>0,
(3.5)

and consequently (3.2) yields

c 0 = m 2 m 2 = m N ,

i.e.,

m=N c 0 < S N / 2 (see (2.11)).
(3.6)

However, recall the definition of S in (2.1),

m= lim n R N ( | u n + V ( ε n y + a ) | u n | 2 ) dyS lim n ( R N | u n | 2 d y ) 2 / 2 =S m 2 / 2 ,

equivalent to m S N / 2 , contradicting (3.6). Thus, (3.4) holds. Using the idea of [21, 25], along a subsequence as n, we may assume that

v n (y):= u n ( ε n y+a)v0weakly in  H 1 ( R N ) .

We now consider t n >0 such that w n := t n v n M a (see (2.3)). By a change of variable x= ε n y+a, it follows that

c 0 J a ( w n ) = J a ( t n v n ) sup t > 0 J a ( t v n ) sup t > 0 ε n N I ε n ( t u n ) + o ( 1 ) = ε n N I ε n ( u n ) + o ( 1 ) = c 0 + o ( 1 ) .
(3.7)

Hence J a ( w n ) c 0 , from which it follows that w n 0 in H 1 ( R N ).

Since ( v n ) and ( w n ) are bounded in H 1 ( R N ) and v n 0 in H 1 ( R N ), the sequence ( t n ) is bounded. Thus, up to a subsequence, t n t 0 0. If t 0 =0, then w n H 1 ( R N ) 0, which does not occur. Hence t 0 >0, and therefore the sequence ( w n ) satisfies

J a ( w n ) c 0 , w n w:= t 0 v0weakly in  H 1 ( R N ) .
(3.8)

For fixed v H 1 ( R N ), define

b(w)= R N wvdyfor all w H 1 ( R N ) .

By the Hölder inequality,

| b ( w ) | | w | 2 | v | 2 w H 1 ( R N ) v H 1 ( R N ) .

Hence b H 1 , the dual space of H 1 ( R N ). Consequently, as n, w n w in H 1 ( R N ) implies b( w n )b(w), i.e.,

R N w n vdx= R N wvdx+o(1).
(3.9)

Since w n converges weakly to w in H 1 ( R N ), w n is bounded in L 2 ( R N ). Thus | w n | 2 2 w n is bounded in L ( 2 ) ( R N ). It then follows that there is a subsequence of ( w n ), still denoted by ( w n ), such that | w n | 2 2 w n converges weakly to some w ˜ in L ( 2 ) ( R N ). Next we will show w ˜ = | w | 2 2 w. Choose a sequence ( K m ) m 1 of open relatively compact subsets, with regular boundaries, of R N covering R N , i.e., R N = m 1 K m . It is easy to see that, by compact embedding, w n w in L q ( K m ) for any q< 2 . Hence w n w a.e. on K m . Hence | w n | 2 2 w n | w | 2 2 w a.e. on K m . By the Brezis and Lieb lemma [26], we conclude that | w n | 2 2 w n | w | 2 2 w strongly in L ( 2 ) ( K m ). Thus w ˜ = | w | 2 2 w a.e. on each K m , and then the diagonal rule implies a.e. on R N . Hence

R N | w n | 2 2 w n vdx= R N | w | 2 2 wvdx+o(1).
(3.10)

Similarly, we have

R N V(a) w n vdx= R N V(a)wvdx+o(1).
(3.11)

By (f2) and (f3),

R N |f(x, w n )v|dxδ R N | w n ||v|dx+ C δ R N | w n | q 1 |v|dxδ | w n | 2 | v | 2 + C δ | w n | q q 1 | v | q .

Hence when R is large enough, we get

B R c ( z n ) R N |f(a, w n )v|dx=o(1).

Noting that w n w in L q ( B R ( z n )), 2q< 2 . Therefore we have

B R ( z n ) R N f(a, w n )vdx= B R ( z n ) R N f(a,w)vdx+o(1).

Hence

R N f(a, w n )vdx= R N f(a,w)vdx+o(1).
(3.12)

By (3.9)-(3.12), we derive that

Δw+V(a)w= | w | 2 2 w+f(a,w)in  R N ,
(3.13)

i.e., J a (w)=0.

For any nN let us consider the measure sequence μ n defined by

R N μ n (dy)= R N ( | w n | 2 + | w n | 2 + | w n | q ) dy.

We assume

R N μ n (dy)l.

By the concentration-compactness lemma [24], there exists a subsequence of ( μ n ) (denoted in the same way) satisfying one of the three following possibilities.

Compactness: There exists a sequence z n R N such that for any δ>0 there is a radius R>0 with the property that

lim n B R ( z n ) μ n (dy)lδ.

Vanishing: For all R>0,

lim n ( sup z R N B R ( z ) μ n ( d y ) ) =0.

Dichotomy: There exists a number a ˜ , 0< a ˜ <l, such that for any δ>0 there is a number R>0 and a sequence ( z n ) with the following property: Given R >R there are non-negative measures μ n 1 , μ n 2 such that

  1. (i)

    0 μ n 1 + μ n 2 μ n ,

  2. (ii)

    supp( μ n 1 ) B R ( z n ), supp( μ n 2 ) R N B R ( z n ),

  3. (iii)

    lim sup n (| a ˜ R N μ n 1 (dy)|+|(l a ˜ ) R N μ n 2 (dy)|)δ.

We are going to rule out the last two possibilities so that compactness holds. Our first goal is to show that vanishing cannot occur. Otherwise,

w n H 1 ( R N ) = R N | w n | 2 d y + R N f ( x , w n ) w n d y R N | w n | 2 d y + δ R N | w n | 2 d y + C δ R N | w n | q d y 0 .

Hence J( w n )0, contradicting c 0 >0.

Now for the harder part. Let η be a smooth nonincreasing cut-off function, defined in [0,), such that η=1 if 0t1; η=0 if t2; 0η1 and | η (t)|2. Also, let η r ()=η( r ). We define

ξ(t)=1η(t),

a nondecreasing function on [0,). Denote by ξ r ()=ξ( r ). We show now that dichotomy does not occur. Otherwise there exists a ˜ (0,l) such that for some R >R and z n R N the function μ n splits into μ n 1 and μ n 2 with the following properties:

R N μ n 1 (dy) a ˜ , R N μ n 2 (dy)l a ˜ .
(3.14)

If we denote

w n 1 = η R (x z n ) w n (x), w n 2 = ξ R (x z n ) w n (x),

(3.14) becomes

R N ( | w n 1 | 2 + | w n 1 | 2 + | w n 1 | q ) d y a ˜ , R N ( | w n 2 | 2 + | w n 2 | 2 + | w n 2 | q ) d y l a ˜ .
(3.15)

Denote by Ω := B R ( z n ) B R ( z n ), then

0= J a ( w n ) , χ Ω w n = Ω ( | w n | 2 + V ( a ) | w n | 2 ) dy Ω | w n | 2 dy Ω f(a, w n ) w n dy.

Using Dichotomy (iii), we get

Ω μ n ( d y ) = R N μ n ( d y ) B R ( z n ) μ n ( d y ) B R c ( z n ) μ n ( d y ) R N μ n ( d y ) B R ( z n ) μ n 1 ( d y ) B R c ( z n ) μ n 2 ( d y ) 0 ,

which implies

Ω | w n | 2 dy0, Ω | w n | 2 dy0, Ω | w n | q dy0.

Hence

Ω ( | w n | 2 + V ( a ) | w n | 2 ) d y = Ω | w n | 2 d y + Ω f ( a , w n ) w n d y Ω | w n | 2 d y + δ Ω | w n | 2 d y + C δ Ω | w n | q d y 0 .

Now we observe that supp  w n 1 supp  w n 2 =, therefore

R N ( | w n | 2 + V ( a ) | w n | 2 ) d y = R N ( | w n 1 | 2 + V ( a ) | w n 1 | 2 ) d y + R N ( | w n 2 | 2 + V ( a ) | w n 2 | 2 ) d y + o ( 1 ) ,
(3.16)
R N | w n | 2 dy= R N | w n 1 | 2 dy+ R N | w n 2 | 2 dy+o(1),
(3.17)
R N F(a, w n )dx= R N F ( a , w n 1 ) dy+ R N F ( a , w n 2 ) dy+o(1),
(3.18)

and

R N f(a, w n ) w n dy= R N f ( a , w n 1 ) w n 1 dy+ R N f ( a , w n 2 ) w n 2 dy+o(1),
(3.19)

where o(1)0 as n.

Recall that ( w n ) M a (see (2.3)), which implies

J a ( w n ) , w n = R N ( | w n | 2 + V ( a ) | w n | 2 ) d y R N | w n | 2 d y R N f ( a , w n ) w n d y = o ( 1 ) .
(3.20)

Then using w n 1 and w n 2 in place of w n , respectively, we get

J a ( w n ) , w n 1 = R N ( | w n 1 | 2 + V ( a ) | w n 1 | 2 ) d y R N | w n 1 | 2 d y R N f ( a , w n 1 ) w n 1 d y J ( w n ) , w n 1 = o ( 1 ) , J a ( w n ) , w n 2 = R N ( | w n 2 | 2 + V ( a ) | w n 2 | 2 ) d y R N | w n 2 | 2 d y R N f ( a , w n 2 ) w n 2 d y J ( w n ) , w n 2 = o ( 1 ) .
(3.21)

There exists t n 1 >0 such that t n 1 w n 1 M 0 , i.e.,

( t n 1 ) 2 R N ( | w n 1 | 2 + V ( a ) | w n 1 | 2 ) d y R N | t n 1 w n 1 | 2 d y R N f ( x , t n 1 w n 1 ) t n 1 w n 1 d y = 0 .
(3.22)

By (f2) and (f3), f(a, w n ) w n δ | w n | 2 + C δ | w n | q , we see t n 1 cannot go zero, that is, t n 1 t 0 1 >0. If t n 1 , by (3.21), (3.22) and (f4), we get

R N | w n 1 | 2 d y + R N f ( x , w n 1 ) w n 1 d y = R N ( | w n 1 | 2 + V ( a ) | w n 1 | 2 ) d y = R N | t n 1 w n 1 | 2 + f ( x , t n 1 w n 1 ) t n 1 w n 1 ( t n 1 ) 2 d y ,
(3.23)

since (f5). By (3.15),

R N | w n 1 | 2 d y + R N f ( x , w n 1 ) w n 1 d y R N | w n 1 | 2 d y + δ R N | w n 1 | 2 d y + C δ R N | w n 1 | q d y C R N ( | w n 1 | 2 + | w n 1 | 2 + | w n 1 | q ) d y C ( a ˜ + o ( 1 ) ) ,
(3.24)

a contradiction. Thus 0< t 0 1 t n 1 C. Assume that t n 1 t 1 , we will show t 1 =1. By (3.21) and (3.22), we have

o ( 1 ) = ( t n 1 ) 2 2 R N | w n 1 | 2 d y + R N f ( x , t n 1 w n 1 ) w n 1 t n 1 d y R N | w n 1 | 2 d y R N f ( x , w n 1 ) w n 1 d y = ( ( t n 1 ) 2 2 1 ) R N | w n 1 | 2 d y + R N ( f ( x , t n 1 w n 1 ) t n 1 f ( x , w n 1 ) ) w n 1 d y .

Hence by the Lebesgue dominated convergence theorem, we get

( ( t 1 ) 2 2 1 ) R N | w 1 | 2 dy+ R N ( f ( x , t 1 w 1 ) t 1 f ( x , w 1 ) ) w 1 dy=0.

By (f5), we have t 1 =1. Similarly, t n 2 t 2 =1. Using this together with (3.16), (3.17), (3.18) and (3.19), we obtain

c 0 +o(1)= J a ( w n )= J a ( w n 1 ) + J a ( w n 2 ) +o(1)= J a ( t n 1 w n 1 ) + J a ( t n 2 w n 2 ) +o(1)2 c 0 +o(1).

Contradiction! Thus dichotomy does not occur.

With vanishing and dichotomy ruled out, we obtain the compactness of a sequence μ n , i.e., there exist z n R N and for each δ>0, there exists R>0 such that

B R c ( z n ) ( | w n | 2 + | w n | 2 + | w n | q ) dyδ.
(3.25)

Then ( z n ) must be bounded, for otherwise (3.25) would imply, in the limit n,

R N ( | w | 2 + | w | 2 + | w | q ) dy C 1 δ
(3.26)

for some positive constants C 1 , independent of δ, which implies w0, contrary to (3.8).

From the foregoing, it follows that there exist bounded nonnegative measures μ ˜ , ν ˜ on R N such that | w n | 2 μ ˜ weakly and | w n | 2 ν ˜ tightly as n. Lemma I.1 in [27] declares that there exist sequences ( x j ) R + N , ( μ ˜ j ),( ν ˜ j )(0,) such that

( 1 ) μ ˜ | w | 2 + j J ˜ μ ˜ j δ x j , ( 2 ) ν ˜ = | w | 2 + j J ˜ ν ˜ j δ x j , ( 3 ) μ ˜ j S ν ˜ j 2 / 2 ,
(3.27)

where δ x j denotes a Dirac measure, j J ˜ . Take x j R + N in the support of the singular part of μ ˜ , ν ˜ . We consider ϕ C c ( R N ) such that

ϕ=1on B( x j ,ε),ϕ=0on B ( x j , 2 ε ) c ,|ϕ|2/ε.
(3.28)

Choosing the test function ϕ w n , from I ε ( w n ),ϕ w n 0, we have

R N ϕ μ ˜ ( d x ) + R N ϕ V ( a ) | w | 2 d x R N ϕ ν ˜ ( d x ) R N ϕ f ( x , w ) w d x = lim n R N w n w n ϕ d x C 1 lim n ( B ( x j , 2 ε ) | w n | 2 d x ) 1 / 2 ε 0 0 .
(3.29)

This reduces to

μ ˜ j = ν ˜ j ,

hence (3.27)(3) states

S ν ˜ j 2 / 2 ν ˜ j ,

i.e.,

ν ˜ j =0or ν ˜ j S N / 2 .

Consequently,

ν ˜ | w | 2 + S N / 2 j J δ x j ,

and hence

R N | w n | 2 dx R N ν ˜ (dx) R N | w | 2 dx+ S N / 2 CardJ,
(3.30)

which implies that the set J is at most finite. Here CardJ is the cardinal numbers of set J. Hence

J a ( w n ) 1 2 J a ( w n ) , w n = ( 1 2 1 2 ) R N | w n | 2 d y + R N ( f ( a , w n ) w n F ( a , w n ) ) d y 1 N R N | w n | 2 d y ,
(3.31)

since

f(x,t)tθF(t)F(t)>0.

When n is large enough, recall c 0 < 1 N S N / 2 (see (2.11)), together with (3.30) and (3.31), we obtain

1 N S N / 2 > c 0 1 N Ω | w | 2 dx+ 1 N j J ν j 1 N S N / 2 ,

a contradiction. Therefore J is empty, that is, | w n | 2 2 | w | 2 2 as n. By the Brezis and Lieb lemma [26] again, we get

| w n w | 2 2 0.
(3.32)

Equation (3.25) and compact embedding theorem imply

R N f(a, w n ) w n dx= R N f(a,w)wdx+o(1).
(3.33)

This together with (3.13), (3.20) and (3.32) allows us to deduce easily

w n H 1 ( R N ) = w H 1 ( R N ) +o(1).

Since H 1 ( R + N ) is a uniformly convex Banach space, hence

w n w H 1 ( R N ) 0.
(3.34)

From (3.32), (3.33) and (3.34), we can obtain

c 0 = lim n J a ( w n )= J a (w),
(3.35)

i.e., w is the ground state solution of (2.2) in view of (3.13). The proof of Proposition 3.1 is complete. □

4 Proof of Theorem 1.1

Proposition 4.1 Suppose f satisfies (f2)-(f4). Then I ε satisfies the ( PS ) c -condition for all c< ε N S N / 2 /N, that is, every sequence ( u n ) in H ε such that I ε ( u n )c, I ε ( u n )0, as n, possesses a convergent subsequence.

Proof Suppose that ( u n ) is a sequence in H ε such that I ε ( u n )c< ε N S N / 2 /N, I ε ( u n )0, as n. Using (f4), by a change of variable x=εy, we obtain that

c + o ( 1 ) u n ε I ε ( u n ) 1 θ I ε ( u n ) , u n = ( 1 2 1 θ ) R N ( ε 2 | u n | 2 + V ( x ) | u n | 2 ) d x + ( 1 θ 1 2 ) R N | u n | 2 d x + 1 θ R N ( f ( x , u n ) u n θ F ( x , u n ) ) d x ( 1 2 1 θ ) u n ε 2 = ( 1 2 1 θ ) u n 2 .
(4.1)

This implies that ( u n ) is bounded in H 1 . Therefore we may assume u n u in H 1 and u n u a.e. Let u n = v n +u. Then

R N ( | u n | 2 + V ( ε y ) | u n | 2 ) d y = R N ( | v n | 2 + V ( ε y ) | v n | 2 ) d y + R N ( | u | 2 + V ( ε y ) | u | 2 ) d y + o ( 1 ) , R N f ( x , u n ) u n d y = R N f ( x , v n ) v n d y + R N f ( x , u ) u d y + o ( 1 ) , R N F ( x , u n ) d y = R N F ( x , v n ) d y + R N F ( x , u ) d y + o ( 1 ) ,

and by the Brezis-Lieb lemma [26],

R N | u n | 2 dy= R N | v n | 2 dy+ R N | u | 2 dy+o(1).

For convenience, we denote by

I(u)= 1 2 R N ( | u | 2 + V ( ε y ) | u | 2 ) dy 1 2 R N | u | 2 dy R N F(x,u)dy
(4.2)

and

I ( u ) , v = R N ( u v + V ( ε y ) u v | u | 2 2 u v f ( x , u ) v ) dy.
(4.3)

It is clear that

I ε (u)= ε N I(u), I ε ( u ) , v = ε N I ( u ) , v .
(4.4)

It is easy to verify that I (u),u=0. Hence we have

o(1)= I ( u n ) , u n = I ( v n ) , v n + I ( u ) , u +o(1)= I ( v n ) , v n +o(1),

and thus

lim n R N ( | v n | 2 + V ( x ) | v n | 2 ) dy= lim n R N | v n | 2 dy=:,

since R N f( | v n | 2 ) | v n | 2 dy0 by (f2) and (f3). If =0, then v n 2 0, hence v n ε 2 = ε N v n 2 0 as n, and we can obtain the desired conclusion. Hence it remains to show that =0. By a change of variable, from

I( v n )=I( v n ) 1 2 I ( v n ) , v n 1 N R N | v n | 2 dy0

and

c=I( u n )+o(1)=I( v n )+I(u)+o(1)I( v n )+o(1),

we get

ε N N c< ε N S N / 2 N ,

i.e.,

< S N / 2 .
(4.5)

By the Sobolev inequalities,

( R N | v n | 2 d x ) 2 / 2 S 1 R N ( | v n | 2 + V ( y ) | v n | 2 ) dy.

Letting n, we get 2 / 2 S 1 , so either S N / 2 which contradicts (4.5) or =0. □

Let ι>0 be fixed. Let η be a smooth nonincreasing cut-off function, defined in [0,), such that η=1 if 0tι; η=0 if t2ι; 0η1 and | η (t)|C for some C>0. For any aΣ, let

ψ ε (a)(x)=η ( | x a | ) w ( x a ε ) ,

where w is the positive ground state of (2.2). We may assume that t ε >0 is the unique positive number such that

max t 0 I ε ( t ψ ε ( ξ ) ( x ) ) = I ε ( t ε ψ ε ( ξ ) ( x ) ) .

Let h(ε) be any positive function tending to 0 as ε0, we define the sublevel

N ˜ ε := { u N ε : I ε ( u ) ε N ( c 0 + h ( ε ) ) } .
(4.6)

By Lemma 4.2 below, N ˜ ε is not empty for ε sufficiently small. By noticing that t ε ψ ε (ξ) N ε , we can define ϕ ε :Ω N ˜ ε as

ϕ ε := t ε ψ ε (ξ).

Lemma 4.2 Uniformly in aΣ, we have

lim ε 0 ε N I ε ( ϕ ε ( a ) ) = c 0 .
(4.7)

Proof Let aΣ. Computing directly, we have

ψ ε ( a ) ( x ) ε 2 = R N ( ε 2 | ψ ε ( a ) ( x ) | 2 + V ( x ) | ψ ε ( a ) ( x ) | 2 ) d x = R N ( ε 2 | w ( x a ε ) η ( | x a | ) + 1 ε η ( | x a | ) w ( x a ε ) | 2 + V ( x ) | η ( | x a | ) w ( x a ε ) | 2 ) d x = R N ( η 2 ( | x a | ) ( | w ( x a ε ) | 2 + V ( x ) | w ( x a ε ) | 2 ) ) d x + R N ( ε 2 | η ( | x a | ) | 2 | w ( x a ε ) | 2 + 2 ε η ( | x a | ) w ( x a ε ) η ( | x a | ) w ( x a ε ) ) d x = : I 1 + I 2 .
(4.8)

By a change of variable y= x a ε , we obtain

I 1 = ε N R N η 2 ( ε y ρ ) ( | w ( y ) | 2 + V ( ε y + a ) | w ( y ) | 2 ) d y = ε N ( R N ( | w | 2 + V ( a ) w 2 ) d y + o ( 1 ) ) as  ε 0 ,
(4.9)

uniformly for aΣ.

I 2 = ε N R N ( ε 2 | η ( | ε y | ρ ) | 2 ω 2 ( y ) + 2 ε η ( | ε y | ρ ) ω ( y ) η ( | ε y | ρ ) ω ( y ) ) d y ε N + 2 { y : ρ / ε | y | 2 ρ / ε } C | ω ( y ) | 2 d y + 2 ε N + 1 { y : ρ / ε | y | 2 ρ / ε } C 1 | ω ( y ) | | ω ( y ) | d y .
(4.10)

By the exponential decay of ω, we get

I 2 = ε N ( o ( 1 ) ) as ε0
(4.11)

uniformly for aΣ. Therefore, in the limit that ε is very small, thanks to (4.8) (4.9) and (4.11), we find

ψ ε ( a ) ( x ) ε 2 = ε N ( R N ( | w | 2 + V ( a ) w 2 ) d y + o ( 1 ) ) .
(4.12)

On the other hand, following the idea of [21, 25], from I ε ( t ε ψ ε (ξ)(x)), t ε ψ ε (ξ)(x)=0, by the change of variables y:=(xξ)/ε, we get

t ε ψ ε ( ξ ) ( x ) ε 2 = R N | t ε ψ ε ( ξ ) ( x ) | 2 d x + R N f ( x , t ε ψ ε ( ξ ) ( x ) ) t ε ψ ε ( ξ ) ( x ) d x = ε N ( R N | t ε η ρ ( ε y ) w ( y ) | 2 d y + R N f ( x , t ε η ρ ( ε y ) w ( y ) ) t ε η ρ ( ε y ) w ( y ) d y ) ,
(4.13)
ε N ψ ε ( ξ ) ( x ) ε 2 = | t ε | 2 2 R N | η ρ ( ε y ) w ( y ) | 2 d y + R N f ( x , t ε η ρ ( ε y ) w ( y ) ) t ε η ρ ( ε y ) w ( y ) d y as  t ε ,
(4.14)

which contradicts (4.12). Thus, up to a subsequence, t ε t 0 0.

Since f has subcritical growth and t ε ψ ε (ξ) N ε , it follows that t 0 >0. Thus, we can take the limit in (4.13) to obtain

R N | ( t 0 w ) | 2 + | t 0 w | 2 dy= R N | t 0 w | 2 dy+ R N f(a, t 0 w) t 0 wdy,
(4.15)

from which it follows that t 0 w M a . Since w also belongs to M a , we conclude that t 0 =1. This and Lebesgue’s theorem imply that

R N | t ε ψ ε ( ξ ) | 2 dx= ε N ( R N | w | 2 d y + o ( 1 ) )
(4.16)

and

Ω F ( x , t ε ψ ε ( a ) ) dx= ε N ( R N F ( a , w ) d y + o ( 1 ) )
(4.17)

uniformly for aΣ. Noting t ε 1, from (4.12), (4.16) and (4.17), we have

ε N I ε ( ϕ ε ( ξ ) ) = 1 2 R N ( | w | 2 + V ( a ) w 2 ) d y + 1 2 R N | w | 2 d y + R N F ( a , w ) d y + o ( 1 ) = c 0 + o ( 1 ) .

Thus (4.7) is proved. □

Let β(u) be the center of mass of u N ε in terms of the L 2 norm:

β(u):= R N x | u | 2 d x R N | u | 2 d x ,u N ε .

Lemma 4.3 Let ε n 0 as n. Then for u n N ˜ ε n ,

β( u n )a,as n,

uniformly for aΣ.

Proof By change of variable x= ε n y+a, we have

β ( u n ) = R N x | u n | 2 d x R N | u n | 2 d x = R N ( ε n y + a ) | v n | 2 d y R N | v n | 2 d y = a + ε n R N y | v n | 2 d y R N | v n | 2 d y = a + ε n R N y | t n v n | 2 d y R N | t n v n | 2 d y .

By Proposition 3.1, w n = t n v n converges strongly in H 1 ( R N ) to w, which is a positive ground state solution of equation (2.2). Thanks to the exponential decay of w (see (2.4)), we obtain β( u n )aΣ as n. This completes the proof of Lemma 4.3. □

Proof of Theorem 1.1 By Proposition 4.1, I ε satisfies the ( PS ) c -condition for all c< ε N S N / 2 /N. Now let us choose a function h(ε)>0 such that h(ε)0 as ε0 and such that ε N ( c 0 +h(ε))< ε N S N / 2 /N is not a critical level for I ε . For such h(ε), let us introduce the set N ˜ ε N ε as in (4.6). Then the standard Ljusternik-Schnirelmann theory implies that I ε has at least cat( N ˜ ε ) critical points on N ˜ ε (also see [21]).

By Lemma 4.3, we can assume that for any δ>0, there exists ε δ >0 such that β( N ˜ ε ) Σ δ for any ε< ε δ . For such ε, by Lemma 4.2, we have I ε ( ϕ ε (a)) ε N ( c 0 +h(ε)) uniformly for aΣ, thus ϕ ε (Σ) N ˜ ε . Recall t ε 1, calculating directly, we get

β ( ϕ ε ( a ) ) = R N x | t ε ψ ε ( a ) ( x ) | 2 d x R N | t ε ψ ε ( a ) ( x ) | 2 d x = R N x | η ( | x a | ) w ( x a ε ) | 2 d x R N | η ( | x a | ) w ( x a ε ) | 2 d x + o ( 1 ) = R N ( ε y + a ) | η ( | ε y | ) w ( y ) | 2 d x R N | η ( | ε y | ) w ( y ) | 2 d x + o ( 1 ) = a + R N ε y | η ( | ε y | ) w ( y ) | 2 d x R N | η ( | ε y | ) w ( y ) | 2 d x = a + o ( 1 ) ,

as ε0 uniformly for aΣ. Hence the map β ϕ ε is homotopical equivalence to the inclusion i:Σ Σ δ for ε small enough. We denote N ˜ ε + = N ˜ ε {u N ε :u0 in  R N }. It is easy to verify that cat( N ˜ ε + , N ˜ ε )cat(Σ, Σ δ ) and cat( N ˜ ε , N ˜ ε )cat(Σ, Σ δ ) (cf. [[19], Lemma 2.2]). Hence we have

cat( N ˜ ε )=cat ( N ˜ ε + , N ˜ ε ) +cat ( N ˜ ε , N ˜ ε ) 2cat(Σ, Σ δ ).

Next we show that if u is a critical point of I ε satisfying I ε (u) ε N ( c 0 +h(ε)), then u cannot change sign. Indeed, if u= u + + u with u + 0 and u 0, then from I ε (u),u=0, we have

I ε ( u + ) , u + =0, I ε ( u ) , u =0.

By change of variable x= ε n y+a, we get

J a ( u + ) , u + =0, J a ( u ) , u =0,as  ε n 0,

i.e.,

u + , u M a .

Also, noting

I ε (u)= ε N J a (u),as  ε n 0.

Hence

c 0 +h(ε) J a (u)= J a ( u + ) + J a ( u ) 2G(a)2 c 0 ,

which is a contradiction. Therefore there exist at least cat(Σ, Σ δ ) nonzero critical points of I ε and thus cat(Σ, Σ δ ) solutions of equation (1.4). □

References

  1. Floer A, Weinstein A: Nonspreading wave packets for the cubic Schrödinger equation with a bounded potential. J. Funct. Anal. 1986, 69: 397-408. 10.1016/0022-1236(86)90096-0

    Article  MathSciNet  MATH  Google Scholar 

  2. Oh Y-G:Existence of semi-classical bound states of nonlinear Schrödinger equations with potentials of the class ( V ) a . Commun. Partial Differ. Equ. 1988, 13: 1499-1519. 10.1080/03605308808820585

    Article  MATH  Google Scholar 

  3. Oh Y-G:Corrections to ‘Existence of semi-classical bound state of nonlinear Schrödinger equations with potentials of the class ( V ) a ’. Commun. Partial Differ. Equ. 1989, 14: 833-834.

    MATH  Google Scholar 

  4. Ambrosetti A, Badiale M, Cingolani S: Semiclassical states of nonlinear Schrödinger equations. Arch. Ration. Mech. Anal. 1997, 140: 285-300. 10.1007/s002050050067

    Article  MathSciNet  MATH  Google Scholar 

  5. Del Pino M, Felmer P: Local mountain passes for semilinear elliptic problems in unbounded domains. Calc. Var. 1996, 4: 121-137. 10.1007/BF01189950

    Article  MathSciNet  MATH  Google Scholar 

  6. Del Pino M, Felmer P: Spike-layered solutions of singularly perturbed elliptic problems in a degenerate setting. Indiana Univ. Math. J. 1999, 48(3):883-898.

    Article  MathSciNet  MATH  Google Scholar 

  7. Rabinowitz P: On a class of nonlinear Schrödinger equations. Z. Angew. Math. Phys. 1992, 43: 270-291. 10.1007/BF00946631

    Article  MathSciNet  MATH  Google Scholar 

  8. Wang X: On concentration of positive bound states of nonlinear Schrödinger equations. Commun. Math. Phys. 1993, 153: 229-244. 10.1007/BF02096642

    Article  MATH  Google Scholar 

  9. Wang X, Zeng B: On concentration of positive bound states of nonlinear Schrödinger equation with competing potential functions. SIAM J. Math. Anal. 1997, 28: 633-655. 10.1137/S0036141095290240

    Article  MathSciNet  MATH  Google Scholar 

  10. Del Pino M, Felmer P: Multi-peak bound states of nonlinear Schrödinger equations. Ann. Inst. Henri Poincaré, Anal. Non Linéaire 1998, 15(2):127-149. 10.1016/S0294-1449(97)89296-7

    Article  MathSciNet  MATH  Google Scholar 

  11. Gui C: Existence of multi-bump solutions for nonlinear Schrödinger equations. Commun. Partial Differ. Equ. 1996, 21: 787-820. 10.1080/03605309608821208

    Article  MathSciNet  MATH  Google Scholar 

  12. Lin C, Ni W-M, Takagi I: Large amplitude stationary solutions to a chemotaxis systems. J. Differ. Equ. 1988, 72: 1-27. 10.1016/0022-0396(88)90147-7

    Article  MathSciNet  MATH  Google Scholar 

  13. Ni WM, Takagi I: On the shape of least-energy solutions to a semilinear Neumann problem. Commun. Pure Appl. Math. 1990, 45: 819-851.

    MathSciNet  MATH  Google Scholar 

  14. Ni WM, Takagi I: Locating the peaks of least-energy solutions to a semilinear Neumann problem. Duke Math. J. 1993, 70: 247-281. 10.1215/S0012-7094-93-07004-4

    Article  MathSciNet  MATH  Google Scholar 

  15. Wang ZQ: On the existence of multiple, single peaked solutions for a semilinear Neumann problem. Arch. Ration. Mech. Anal. 1992, 120: 375-399. 10.1007/BF00380322

    Article  MATH  MathSciNet  Google Scholar 

  16. Gui C: Multi-peak solutions for a semilinear Neumann problem. Duke Math. J. 1996, 84: 739-769. 10.1215/S0012-7094-96-08423-9

    Article  MathSciNet  MATH  Google Scholar 

  17. Li YY: On a singularly perturbed equation with Neumann boundary condition. Commun. Partial Differ. Equ. 1998, 23: 487-545.

    MATH  MathSciNet  Google Scholar 

  18. Wei J, Winter M: Multiple boundary spike solutions for a wide class of singular perturbation problems. J. Lond. Math. Soc. 1999, 59(2):585-606. 10.1112/S002461079900719X

    Article  MathSciNet  Google Scholar 

  19. Cingolani S, Lazzo M: Multiple positive solutions to nonlinear Schrödinger equations with competing potential functions. J. Differ. Equ. 2000, 160: 118-138. 10.1006/jdeq.1999.3662

    Article  MathSciNet  MATH  Google Scholar 

  20. Ambrosetti A, Malchiodi A: Nonlinear Analysis and Semilinear Elliptic Problems. Cambridge University Press, Cambridge; 2007.

    Book  MATH  Google Scholar 

  21. Liu, W, Zhao, P: Critical semilinear Neumann problem with magnetic fields. Preprint

  22. Berestycki H, Lions PL: Nonlinear scalar field equations I. Existence of a ground state. Arch. Ration. Mech. Anal. 1983, 82: 313-345.

    MathSciNet  MATH  Google Scholar 

  23. Zhang J, Zou W: A Berestycki-Lions theorem revisited. Commun. Contemp. Math. 2012., 14(5): Article ID 1250033

    MATH  Google Scholar 

  24. Lions PL: The concentration-compactness principle in the calculus of variation. The locally compact case. II. Ann. Inst. Henri Poincaré, Anal. Non Linéaire 1984, 1: 223-283.

    MATH  Google Scholar 

  25. Alves CO, Figueiredo GM, Furtado MF: Multiple solutions for a nonlinear Schrödinger equation with magnetic fields. Commun. Partial Differ. Equ. 2011, 36: 1565-1586. 10.1080/03605302.2011.593013

    Article  MathSciNet  MATH  Google Scholar 

  26. Brezis H, Lieb E: A relation between pointwise convergence of functions and convergence of functional. Proc. Am. Math. Soc. 1983, 88: 486-490.

    Article  MathSciNet  MATH  Google Scholar 

  27. Lions PL: The concentration-compactness principle in the calculus of variations. The limit case. I. Rev. Mat. Iberoam. 1985, 1: 145-201.

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

This work was partially supported by the National Natural Science Foundation of China (11261052). The authors are grateful to Prof. Guowei Dai for pointing out several mistakes and valuable comments.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Wulong Liu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

WL carried out the genetic studies, participated in the sequence alignment and drafted the manuscript. PZ checked the references.

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

Liu, W., Zhao, P. Multiple solutions to nonlinear Schrödinger equations with critical growth. Bound Value Probl 2013, 199 (2013). https://doi.org/10.1186/1687-2770-2013-199

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords