Next Article in Journal
Two New Modified Regularized Methods for Solving the Variational Inclusion and Null Point Problems
Previous Article in Journal
Hybrid Newton–Sperm Swarm Optimization Algorithm for Nonlinear Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interior Bubbling Solutions for an Elliptic Equation with Slightly Subcritical Nonlinearity

by
Khalil El Mehdi
1,2,*,† and
Fatimetou Mohamed Salem
2,†
1
Department of Mathematics, College of Science, Qassim University, P.O. Box 6644, Buraydah 51452, Saudi Arabia
2
Faculté des Sciences et Techniques, Université de Nouakchott, Nouakchott 2373, Mauritania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Accepted: 15 March 2023 / Published: 17 March 2023
(This article belongs to the Section Difference and Differential Equations)

Abstract

:
In this paper, we considered the Neumann elliptic equation ( P ε ) : Δ u + K ( x ) u = u ( n + 2 ) / ( n 2 ) ε , u > 0 in Ω , u / ν = 0 on Ω , where Ω is a smooth bounded domain in R n , n 6 , ε is a small positive real and K is a smooth positive function on Ω ¯ . Using refined asymptotic estimates of the gradient of the associated Euler–Lagrange functional, we constructed simple and non-simple interior bubbling solutions of ( P ε ) which allowed us to prove multiplicity results for ( P ε ) provided that ε is small. The existence of non-simple interior bubbling solutions is a new phenomenon for the positive solutions of subcritical problems.

1. Introduction

In this paper, we consider the following nonlinear elliptic equation
( P q ) : Δ u + μ u = u q , u > 0 in Ω u ν = 0 on Ω ,
where 1 < q < , Ω is a smooth bounded domain in R n , n 3 , and μ is a positive real.
Equation ( P q ) arises in many areas of applied sciences, including the Gierer–Meinhardt model in biological pattern formation [1] and the Keller– Segel model in chemotaxis [2].
The problem ( P q ) has been studied extensively in the last two decades, and it is well known that the situation depends on the exponent q, the parameter μ , and the dimension n.
For a subcritical case, that is q < n + 2 n 2 ; Lin, Ni, and Takagi [3] showed that ( P q ) does not have any non constant solution for small μ , whereas nonconstant solutions appear for large μ which blow up, as μ goes to infinity, at one or several points. The least energy solution blows up, for large μ , at the maximum point of the mean curvature of the boundary [4,5]. This concentration phenomenon was shown by several authors (see [6] for good review).
Due to the lack of compactness of the associated variational Euler–Lagrange functional, the critical case, that is q = n + 2 n 2 , is much more difficult to handle. For n = 3 , Zhu [7] showed that when Ω is convex ( P q ) has only constant solution for small μ (see also [8]). However, for μ , small and n = 4 , 5 , 6 , nonconstant solutions exist (see [9] when Ω is a ball and [8,10] for general domains). In [11,12], the authors proved the existence of solutions for large μ . As in the subcritical case, solutions blow up, for large μ , at one or several boundary points as μ goes to infinity (see [13,14,15,16,17,18,19,20,21,22]). In [23], Rey proved, in contrast with the subcritical case, that at least one blow-up point must be on the boundary. In [24,25], the authors studied the problem for fixed μ , when the exponent q is close to the critical one, i.e., q = n + 2 n 2 ε and ε is a small nonzero number. When the dimension n = 3 , they proved the existence of single interior peak solution and for n 4 , they also constructed a single boundary peak solution.
In this work, we consider the case where the constant μ in ( P q ) is replaced by a function K, and we are going to prove that multiple interior bubbling solutions exist provided that q is close enough to the critical exponent. More precisely, we considered the following nonlinear elliptic equation with subcritical nonlinearity
( P ε ) : Δ u + K ( x ) u = u p ε , u > 0 in Ω u n = 0 on Ω ,
where Ω is a smooth bounded domain in R n , n 6 , ε is a small positive real, p + 1 = ( 2 n ) / ( n 2 ) is the critical Sobolev exponent for the embedding H 1 ( Ω ) L q ( Ω ) , and K is a C 3 positive function on Ω ¯ .
As we mentioned above, there have been many works on ( P q ) (i.e., when μ is constant). On the contrary, in the case of problem ( P ε ) there is a very poor literature. To our knowledge, the only results are made on the whole space [26,27,28].
In the first part of this work, our aim is to construct simple interior bubbling solutions for problem ( P ε ) . This construction allows us to prove a multiplicity result for problem ( P ε ) in connection with the number of critical points of K.
In the second part, we constructed interior bubbling solutions with clustered bubbles at a critical point of K. Such solutions give a new phenomenon for the positive solutions, that is, the existence of non simple blow up points in the interior for the subcritical problem. This phenomenon is known for changing-sign solution with Dirichlet boundary conditions on some symmetric domains (see [29]).
The remainder of the paper is organized as follows: in Section 2, we state our main results. In Section 3, we set up the analytical framework of the problems ( P ε ) , introduce the neighbourhood at infinity and its parametrization. We also give a precise estimates of the infinite dimensional part. Section 4 is devoted to the expansion of the gradient of the associated functional. In Section 5, we prove Theorems 1 and 2, while Section 6 is devoted to the proof Theorems 3 and 4.

2. Main Results

To state ours results, we need to introduce some notation. Problem ( P ε ) has a variational structure. Indeed, its solutions correspond to the positive critical points of the functional
I ε ( u ) : = 1 2 u 2 1 p + 1 ε Ω | u | p + 1 ε , with   p : = n + 2 n 2 .
defined on H 1 ( Ω ) equipped with the norm . and its corresponding inner product defined by:
u 2 = Ω u 2 + K u 2 , ( u , v ) = Ω u v + K u v .
For ε = 0 , the functional I ε fails to satisfy the Palais Smale condition, and the reason for such a lack of compactness is the existence of almost solution of the equation ( P ε ) . These almost solutions, called bubbles, are defined as follows:
δ a , λ ( x ) = c 0 λ ( n 2 ) / 2 ( 1 + λ 2 x a 2 ) ( n 2 ) / 2 , λ > 0 , a , x R n , c 0 = ( n ( n 2 ) ) n 2 4 ,
which are the only solutions to the problem [30]
Δ u = u n + 2 n 2 , u > 0 in R n .
The aim of our first result is to construct simple interior bubbling solutions. More precisely, we prove:
Theorem 1. 
Let n 6 and let K : Ω ¯ R be a smooth positive function having m non- degenerate critical points y 1 , ..., y m . Then, for any N m , there exists ε 0 > 0 small such that for any ε ( 0 , ε 0 ] , problem ( P ε ) admits a solution u ε , y i 1 , , y i N satisfying: u ε , y i 1 , , y i N develops exactly one bubble at each point y i j and converges weakly to zero in H 1 ( Ω ) as ε 0 . More precisely there exist λ 1 , ε , ..., λ N , ε having the same order as ε 1 / 2 and points a i j , ε y i j for all j such that
| | u ε , y i 1 , , y i N j = 1 N δ a i j , ε , λ i j , ε | | 0 , a s ε 0 .
Theorem 1 allows us to obtain the following multiplicity result in connection with the number of critical points of K.
Theorem 2. 
Let n 6 and let K : Ω ¯ R be a smooth positive function having m non-degenerate critical points. Then, for ε > 0 small, ( P ε ) admits at least 2 m 1 solutions.
Our aim in the next result was to construct interior bubbling solutions with clustered bubbles at a critical point of K. To this aim, we introduced some notation. For N N and y a critical point of K, we defined the following function
F N , y ( z 1 , , z N ) = j = 1 N D 2 K ( y ) ( z j , z j ) k r 1 | z k z r | n 2 ,
where ( z 1 , , z N ) ( R n ) N such that z i z j if i j .
Our result reads as follows:
Theorem 3. 
Let n 6 and let y be a non-degenerate critical point of K. Let N N * and assume that the function F N , y has a non-degenerate critical point ( z ¯ 1 , , z ¯ N ) . Then, for any integer, there exists ε 0 > 0 small such that for any ε ( 0 , ε 0 ] , problem ( P ε ) admits a solution u ε , y with the following properties
u ε , y = i = 1 N α i , ε δ a i , ε , λ i , ε + v ε ,
with
v ε 0 a s ε 0 ,
| α i , ε 1 | c 1 ε ln 2 ε i
ε c 2 1 λ i , ε c 2 ε i
| a i , ε y + σ ε n 4 2 n z ¯ i | ε n 4 2 n η 0 i ,
where c 1 , c 2 are positive constants, η 0 is a small positive real and σ is the constant defined by (135). Moreover, if for each N, F N , y has a non-degenerate critical point, then, problem ( P ε ) has an arbitrary number of non-constant distinct solutions provided that ε is small.
Remark 1. 
If we assume that D 2 K ( y ) ( · , · ) is negative definite, we see that for each N, F N , y tends to near the boundary of its definition domain and therefore F N , y achieves its maximum.
Our last result deals with the case of many interior blow up points. Namely, we prove:
Theorem 4. 
Let n 6 and let y 1 , , y m be non-degenerate critical points of K. For each k m and N 1 , , N m N * , we assume that the function F N k , y k has a non-degenerate critical point ( z ¯ k , 1 , , z ¯ k , N k ) . Then, there exists ε 0 > 0 small such that for any ε ( 0 , ε 0 ] , problem ( P ε ) admits a solution u ε , y 1 , , y m satisfying
u ε , y 1 , , y m = i = 1 N 1 α 1 , i , ε δ a 1 , i , ε , λ 1 , i , ε + + i = 1 N m α m , i , ε δ a m , i , ε , λ m , i , ε + v ε ,
with v ε 0 as ε 0 and for each k m the coefficients α k , i , ε , the speeds λ k , i , ε and the points a k , i , ε satisfy the properties (5), (6), and (7), respectively.
Remark 2. 
We believe that our results should also be true in small dimensions, but we need more careful computations. We will come back to these dimensions in future work.
The strategy of the proof of our results is based on refined asymptotic estimates of the gradient of the associated Euler–Lagrange functional in the so-called neighborhood of critical points at infinity [31]. The aim is to find the equilibrium conditions satisfied by the concentration parameters. These balancing conditions are obtained by testing the equation by vector fields which are the dominant term of the gradient with respect to the parameters of concentration. The analysis of these balancing conditions allows us to obtain all the information we need to prove our results.

3. Analytical Framework

Let n 6 and y 1 , , y m be distinct critical points of K. The analytical framework is similar to the one considered in [32]. We are going to proceed to the suitable parametrization of the associated variational problem. For N m and ν 0 positive small, we introduce the following set
O ( N , ν 0 ) = { ( α , λ , a ) ( R + ) N × ( R + ) N × Ω N : | α i 1 | < ν 0 , λ i > ν 0 1 , ε ln λ i < ν 0 , d ( a i , ) > c , and ε i j < ν 0 } ,
where c is a positive constant and ε i j : = λ i λ j + λ j λ i + λ i λ j | a i a j | 2 ( 2 n ) / 2 .
  • For each ( α , λ , a ) O ( N , ν 0 ) , we consider u = i = 1 N α i δ a i , λ i + v , with v E a , λ , where
    E a , λ = { v H 1 ( Ω ) / Ω v δ i = Ω v δ i λ i = Ω v δ i a i j = 0 , 1 i N , 1 j n } ,
    with δ i = δ a i , λ i .
As usual in these type of equations, we first studied the infinite dimensional variable v. To this aim, we need to perform and expansion of I ε with respect to v E a , λ . Writing U = i = 1 N α i δ i , we see that
I ε U + v = 1 2 Ω ( U + v ) 2 + 1 2 Ω K ( U + v ) 2 1 p + 1 ε Ω U + v p + 1 ε = I ε ( U ) + f ε , v + 1 2 Q ε ( v ) + R ε ( v )
with
f ε , v = Ω K U v Ω U p ε v ,
Q ε ( v ) = Ω v 2 + Ω K v 2 ( p ε ) Ω U p 1 ε v 2 ,
R ε ( v ) = O ( v m i n ( 3 , p + 1 ε ) ) .
In addition the derivatives of R ε satisfy
R ε ( v ) = O ( v m i n ( 2 , p ε ) ) and R ε ( v ) = O ( v m i n ( 1 , p 1 ε ) ) .
Now, to simplify the expression of Q ε , we prove the following lemma
Lemma 1. 
Let ( α , λ , a ) O ( N , ν 0 ) and v E a , λ . Then, for ε small, the following statement holds
Q ε ( v ) = Ω v 2 + Ω K v 2 p i = 1 N Ω δ i p 1 v 2 + o ( v 2 ) ,
where Q ε is defined in (11).
Proof. 
First, since ε ln λ i 0 , we have
δ i ε = c 0 ε λ i ε ( n 2 ) 2 1 + O ( ε ln ( 1 + λ i 2 x a i 2 = 1 + o ( 1 ) .
Second, we observe that
ε Ω U p 1 ε v 2 c ε i = 1 N Ω δ i p 1 ε v 2 c i = 1 N Ω δ i p + 1 2 n v 2 c ε v 2
and
Ω U p 1 ε v 2 = i = 1 N Ω ( α i δ i ) p 1 ε v 2 + i j O Ω ( δ i δ j ) p 1 2 v 2 .
However, using estimate E 2 of [31], we see that
Ω ( δ i δ j ) p 1 2 v 2 c Ω ( δ i δ j ) n n 2 2 n v 2 c ( ε i j n n 2 ln ε i j 1 ) 2 n v 2 c ε 2 n 2 ( ln ε i j 1 ) 2 n v 2 .
Using (14) and the fact that | α i 1 | is small, we obtained
Ω ( α i δ i ) p 1 ε v 2 = ( 1 + o ( 1 ) ) Ω δ i p 1 v 2 = Ω δ i p 1 v 2 + o ( v 2 ) .
Combining (15)–(18), our lemma follows. □
To proceed further, we need to prove the uniform coercivity of the quadratic form Q ε . Notice that such a kind of coercivity has been proved by Bahri [31] in the case of Dirichlet boundary conditions. Such a result has been adapted in [33] for our case when the concentration points do not approach each other. Hpwever, here we need a more general result which holds even if the points are close from each other. More precisely, without any assumption on the distance between the points, we have:
Proposition 1. 
Let ( α , λ , a ) O ( N , ν 0 ) and v E a , λ . Then, for ε small, there exists a constant C > 0 such that
Q ε ( v ) C v 2 v E a , λ .
Proof. 
By Lemma 1, we have
Q ε ( v ) = Q ( v ) + o ( v 2 ) ,
where
Q ( v ) = Ω v 2 + Ω K v 2 p i = 1 N Ω δ i p 1 v 2 .
Now, let P v be the projection of v onto H 0 1 ( Ω ) defined by
Δ P v = Δ v in Ω , P v = 0 on Ω
and define w = v P v . Clearly, we have
Ω P v w = Ω P v Δ w + Ω P v w ν = 0 .
Thus
Ω | v | 2 = Ω | P v | 2 + Ω | w | 2 .
Now, let Ω 0 be a compact subset of Ω such that B ( a i , r 0 ) Ω 0 for all i with a fixed positive real r 0 . For y Ω 0 , since w is a harmonic function in Ω, we see that
| w ( y ) | c Ω | G 0 ν ( x , y ) | | w ( x ) | d x C ( Ω 0 ) Ω | w ( x ) | d x C w H 1 ,
where G 0 is the Green’s function of the Laplace operator with Dirichlet boundary conditions.
  • In the same way, we have
| w ( y ) | C ( Ω 0 ) w H 1 y Ω 0 .
Next, for 1 i N , taking
ψ i δ i , λ i δ i λ i , 1 λ i δ i ( a i ) j , 1 j n ,
we observe that
Ω P v ψ i = Ω v ψ i Ω w ψ i = Ω 0 w ψ i + O 1 λ i ( n 2 ) / 2 Ω \ Ω 0 | w | = O w H 1 Ω 0 | ψ i | + w H 1 λ i ( n 2 ) / 2 = O w H 1 λ i ( n 2 ) / 2 ,
where we have used (23) and the fact that Ω 0 | ψ i | = O 1 λ i ( n 2 ) / 2 . Now, we write
Q ( v ) = Q 0 ( P v ) + Ω | w | 2 + Ω K v 2 p i = 1 N Ω δ i p 1 ( 2 P v + w ) w ,
where
Q 0 ( P v ) = Ω | P v | 2 p i = 1 N Ω δ i p 1 ( P v ) 2 .
However, using (22), we obtain
Ω δ i p 1 ( 2 | P v | + | w | ) | w | = Ω 0 δ i p 1 ( 2 | P v | + | w | ) | w | + Ω \ Ω 0 δ i p 1 ( 2 | P v | + | w | ) | w | C w H 1 ( P v H 1 + w H 1 ) Ω 0 δ i 8 n n 2 4 n + 2 2 n + Ω \ Ω 0 δ i p + 1 2 n .
By easy computations, we have
Ω 0 δ i 8 n n 2 4 n + 2 2 n + Ω \ Ω 0 δ i p + 1 2 n = ln 2 / 3 λ i λ i 2 if n = 6 1 λ i 2 if n 6 .
Thus
Ω δ i p 1 ( 2 | P v | + | w | ) | w | = o ( v 1 H 1 2 + v 2 H 1 2 ) .
To proceed further, let
v 1 = P v in Ω , v 1 = 0 in R n \ Ω .
Clearly v 1 D 1 , 2 ( R n ) . Now, we split
v 1 = i = 1 N ( n + 2 ) γ i ψ i + v ¯ 1 .
where v ¯ 1 satisfies R n v ¯ 1 ψ i = 0 for all i and ψ i goes through all the N ( n + 2 ) functions defined by (24).
Note that, by easy computations, the following estimates can be verified
Ω | δ i | 2 = c + O 1 λ i n 2 ; Ω δ i δ k = O ε i k for k i ,
Ω δ i ( λ i δ i λ i ) = O 1 λ i n 2 ; Ω δ i ( λ k δ k λ k ) = O ε i k for k i ,
Ω δ i ( 1 λ i δ i a i j ) = O 1 λ i n 1 ; Ω δ i ( 1 λ k δ k a k j ) = O ε i k for k i ,
Ω δ i ( λ i δ i λ i ) = O 1 λ i n 2 ; Ω ( λ i δ i λ i ) δ k = O ( ε i k ) for k i ,
Ω ( λ i δ i λ i ) ( λ i δ i λ i ) = c + O 1 λ i n 2 ; Ω ( λ i δ i λ i ) ( λ k δ k λ k ) = O ε i k for k i ,
Ω ( λ i δ i λ i ) ( 1 λ i δ i a i j ) = O 1 λ i n 1 ; Ω ( λ i δ i λ i ) ( 1 λ k δ k a k j ) = O ( ε i k ) for k i ,
Ω δ i ( 1 λ i δ i a i j ) = O 1 λ i n 1 ; Ω ( 1 λ i δ i a i j ) δ k = O ( ε i k ) for k i ,
Ω ( 1 λ i δ i a i j ) ( λ i δ i λ i ) = O 1 λ i n 1 ; Ω ( 1 λ i δ i a i j ) ( λ k δ k λ k ) = O ε i k for k i ,
Ω ( 1 λ i δ i a i j ) ( 1 λ i δ i a i l ) = c δ j l + O 1 λ i n 1 ; Ω ( 1 λ i δ i a i j ) ( 1 λ k δ k a k j ) = O ( ε i k ) for k i ,
where δ j l denotes the Kronecker delta, i.e., δ j l = 1 if j = l and δ j l = 0 if j l .
Now, using (25), (29)–(37), we obtain
R n v 1 ψ i = c γ i + o ( γ j ) = Ω P v ψ i = O w H 1 λ l 2 n 2 = o ( w H 1 ) .
This implies that
R n | v ¯ 1 | 2 = R n | v 1 | 2 + o ( w H 1 2 ) .
Using (39), we get
Q 0 ( P v ) = R n | v 1 | 2 p R n δ i p 1 v 1 2 = R n | v ¯ 1 | 2 p R n δ i p 1 v ¯ 1 2 + o ( v ¯ 1 L 2 2 + w H 1 2 ) .
Combining (38) and Proposition 3.1 of [31], we obtain
Q 0 ( P v ) ρ R n | v ¯ 1 | 2 + o ( v ¯ 1 L 2 2 + w H 1 2 ) ρ 2 R n | v 1 | 2 + o ( w H 1 2 ) = ρ 2 Ω | P v | 2 + o ( w H 1 2 ) .
Therefore (41) and (26) imply that
Q ( v ) ρ 2 Ω | P v | 2 + Ω | w | 2 + μ Ω v 2 + o ( P v 2 + w H 1 2 ) C v 2
which completes the proof of the proposition. □
Now, we are able to deal with the infinite dimensional variable v. More precisely, we prove the following:
Proposition 2. 
Let ( α , λ , a ) O ( N , ν 0 ) . Then, for ε small, there exists a unique v ¯ E a , λ satisfying
I ε ( i = 1 N α i δ a i , λ i + v ¯ ) , h = 0 h E a , λ .
Moreover, v ¯ satisfies v ¯ c T ( ε , a λ ) with
T ( ε , a , λ ) : = ε + i = 1 N 1 λ i n 2 2 + T 1 ( a , λ ) + i = 1 N T 2 ( λ i ) ,
where
T 1 ( a , λ ) = i j ε i j n + 2 2 ( n 2 ) ln n + 2 2 n ε i j 1 a n d T 2 ( λ i ) = ln 2 3 λ i / λ i 2 i f n = 6 c λ i 2 i f n 7 .
Proof. 
Applying Proposition 1 and estimate (13), we see that the implicit function theorem implies, for ε small, the existence of v ¯ E a , λ satisfying the estimate v ¯ = O ( f ε ) where f ε is defined by (10). However, we have
Ω K U v c Ω δ i v c v i = 1 N Ω δ i 2 n n + 2 n + 2 2 n .
By easy computations, we obtain
Ω δ i 2 n n + 2 n + 2 2 n = T 2 ( λ i ) .
Using the fact that, for a , b > 0 and β > 0 , we have
( a + b ) β = a β + b β + O ( ( a b ) β 2 ) ,
we see that
Ω U p ε v = i = 1 N α i p ε Ω δ i p ε v + i j O Ω ( δ i δ j ) p ε 2 v .
Using (14) and estimate E 2 of [31], we obtain
Ω ( δ i δ j ) p ε 2 v c Ω ( δ i δ j ) p 2 v c v Ω ( δ i δ j ) n n 2 n + 2 2 n c v ε i j n n 2 ln ε i j 1 n + 2 2 n .
For the first term in right hand side of (44), we use (14) to derive
Ω δ i p ε v = c 0 ε λ i ε ( n 2 2 Ω δ i p v + O ε Ω δ i p | v | ln ( 1 + λ i 2 x a i 2 ) = c 0 ε λ i ε ( n 2 2 Ω δ i p v + O ε v R n δ i p + 1 ln 2 n n + 2 ( 1 + λ i 2 x a i 2 ) ( n + 2 ) / ( 2 n ) = c 0 ε λ i ε ( n 2 2 Ω δ i p v + O ( ε v ) .
However, since v E a , λ , we have
Ω δ i p v = Ω Δ δ i v = Ω δ i ν v c λ i n 2 2 Ω v c λ i n 2 2 v .
Combining (42)–(47), we easily derive our proposition. □

4. Expansion of the Gradient in the Neighbourhood at Infinity

In this section, we are going to give asymptotic expansions of the gradient of the Euler-Lagrange function I ε . We start by the expansion with respect to the gluing parameter α i s . Namely, we prove:
Proposition 3. 
Let ( α , λ , a ) O ( N , ν 0 ) , v E a , λ and u = i = 1 N α i δ a i , λ i + v . Then, for ε small, the following statement holds
I ε ( u ) , δ i = α i S n 1 α i p 1 ε λ i ε ( n 2 ) 2 + O R 1 ,
where
S n = c 0 2 n n 2 R n d x ( 1 + x 2 ) n a n d R 1 = ε + 1 λ i 2 + j i 1 ( λ i λ j ) n 2 2 + j i ε i j + v 2 .
Proof. 
For h H 1 ( Ω ) , we notice that
I ε ( u ) , h = Ω u h + Ω K u h Ω u p 1 ε u h ,
where
u = i = 1 N α i δ i + v : = U + v .
Taking h = δ i in (48), we see that
Ω u δ i = α i Ω δ i 2 + j i α j Ω δ j δ i .
However, we have
Ω δ j δ i = Ω Δ δ j δ i + Ω δ j ν δ i = Ω δ j p δ i + Ω δ j ν δ i = O ε i j + 1 ( λ i λ j ) n 2 2
and
α i Ω δ i 2 = α i S n + O 1 λ i n 2 .
We also have
Ω K u δ i = α i Ω K δ i 2 + j i α j Ω K δ j δ i + Ω K v δ i .
For the last term in the right hand side of (53), we have
Ω K v δ i c v δ i c v T 2 ( λ i ) ,
where T 2 is defined in Proposition 2.
  • For the other term, using Lemma 2.2 of [34], we observe that
Ω K δ j δ i c Ω δ j δ i c ε i j for i j
Ω K δ i 2 c Ω δ i 2 = O 1 λ i 2 .
It remains to deal with the last term in the right hand side of (48). To this aim, we write
Ω 1 = { x : v 1 2 U } , Ω 2 = Ω \ Ω 1 ,
where U is defined in (49).
Observe that
Ω 2 u p 1 ε u δ i Ω 2 v p + 1 ε c v p + 1 ε ,
Ω 1 u p 1 ε u δ i = Ω 1 U p ε δ i + ( p ε ) Ω 1 U p 1 ε v δ i + O Ω 1 U p 2 ε v 2 δ i ,
Ω 1 U p 2 ε v 2 δ i Ω 1 U p 1 ε v 2 = O ( v 2 ) ,
Ω 1 U p ε δ i = α i p ε Ω δ i p + 1 ε + O j i Ω δ j p ε δ i + j i Ω δ i p ε δ j + O Ω 2 v p + 1 ε = α i p ε Ω δ i p + 1 ε + j i O ( ε i j ) + O ( v 2 ) ,
Ω 1 U p 1 ε v δ i = Ω U p 1 ε v δ i + O Ω 2 v p + 1 ε = α i p 1 ε Ω δ i p ε v + O j i Ω i δ i p 1 ε | v | δ j + Ω \ Ω i δ j p 1 | v | δ i + O ( v 2 ) ,
where
Ω i = { x Ω : j i α j δ j ( x ) α i δ i ( x ) } and ( Ω i ) c = Ω \ Ω i .
Observe that
j i Ω i δ i p 1 ε | v | δ j + Ω \ Ω i δ j p 1 | v | δ i c v j i Ω ( δ i δ j ) p / 2 | v | v Ω ( δ i δ j ) n n 2 n + 2 2 n = O v ε i j n + 2 2 ( n 2 ) ln n + 2 2 n ε i j 1 .
Lastly, we write
α i p ε Ω δ i p + 1 ε = α i p ε c 0 ε λ i ε ( n 2 ) 2 Ω δ i p + 1 + O ε Ω δ i p + 1 ln ( 1 + λ i 2 | x a i | 2 ) = α i p ε c 0 ε λ i ε ( n 2 ) 2 S n + O ε + 1 λ i n
and, since v E a , λ ,
Ω δ i p ε v = c 0 ε λ i ε ( n 2 ) 2 Ω δ i p v + O ε Ω δ i p | v | ln ( 1 + λ i 2 x a i 2 ) = c 0 ε λ i ε ( n 2 ) 2 Ω Δ δ i v + O ε v = O v λ i ( n 2 ) / 2 + ε v .
Combining (48)–(65), we easily derive our proposition. □
Now we are going to provide a balancing formula involving the rate of the concentration λ i and the self-interaction of bubbles ε i j . Namely, we prove:
Proposition 4. 
Let ( α , λ , a ) O ( N , ν 0 ) , v E a , λ and u = i = 1 N α i δ a i , λ i + v . Then, for ε small, the following statement holds for all 1 i N
I ε ( u ) , λ i δ i λ i = α i p c 2 ε α i c 1 K ( a i ) λ i 2 + O ( R 2 ) ,
where
R 2 = ε 2 + v 2 + j = 1 N 1 ( λ j λ i ) n 2 2 + j i ε i j + v T 2 ( λ i ) + T 3 ( λ i ) ,
c 1 = n 2 2 c 0 2 R ( x 2 1 ) d x ( 1 + x 2 ) n 1 > 0 , c 2 = ( n 2 2 ) 2 R n c 0 2 n n 2 ( x 2 1 ) ( 1 + x 2 ) n + 1 ln ( 1 + x 2 ) d x
and
T 3 ( λ i ) λ i 4 i f n 7 λ i 4 ln λ i i f n = 6 .
Proof. 
We will take h = λ i δ i λ in (48) and estimate each term.
  • Observe that:
Ω K u λ i δ i λ i = j = 1 N α j Ω K δ j λ i δ i λ i + Ω K v λ i δ i λ i ,
Ω K v λ i δ i λ i c Ω v δ i c v T 2 ( λ i ) ,
| Ω K δ j λ i δ i λ i | c Ω δ j δ i c ε i j j i ,
where we have used in the last formula Lemma 2.2 of [34].
For r > 0 small, we have
Ω K δ i λ i δ i λ i = B a i , r K δ i λ i δ i λ i + Ω \ B a i , r K δ i λ i δ i λ i
Ω \ B a i , r K δ i λ i δ i λ i C Ω \ B a i , r δ i 2 c Ω \ B a i , r d x λ i n 2 | x a i | 2 n 4 C λ i n 2 .
On B a i , r : = B i , we write
B i K δ i λ i δ i λ i = K a i B i δ i λ i δ i λ i + O B i x a i 2 δ i 2 .
Notice that
B i | x a i | 2 δ i 2 c B i λ i n 2 | x a i | 2 1 + λ i 2 | x a i | 2 n 2 d x c λ i 4 0 λ i r t n + 1 1 + t 2 n 2 d t = O ( T 3 ( λ i ) ) ,
B i δ i λ i δ i λ i = ( n 2 ) 2 c 0 2 B i λ i n 2 1 λ i 2 x a i 2 1 + λ i 2 x a i 2 n 1 = ( n 2 ) 2 c 0 2 λ i 2 R n 1 | x | 2 1 + | x | 2 n 1 d x + O 1 λ i n 4 = c 1 λ i 2 + O 1 λ i n 2 .
Combining (66)–(73), we obtain
Ω K u λ i δ λ i = K ( a i ) c 1 α i λ i 2 + O 1 λ i n 2 + j i ε i j + v T 2 ( λ i ) + T 3 ( λ i ) .
Next, we are going to estimate the first term in the right hand side of (48) with h = λ i δ i λ i . We observe that
Ω u λ i δ i λ i = j = 1 N α j Ω δ j λ i δ i λ i = j = 1 N α j Ω Δ δ j λ i δ i λ i + Ω δ j ν λ i δ i λ i = j = 1 N α j Ω δ j p λ i δ i λ i + Ω δ j ν λ i δ i λ i .
However, using estimate E 1 of [31], we have
Ω δ j p λ i δ i λ i c Ω δ j p δ i c ε i j j i .
We also have
Ω δ i p λ i δ i λ i = R n δ i p λ i δ i λ i R n \ Ω δ i p λ i δ i λ i = R n \ Ω δ i p λ i δ i λ i = O 1 λ i n ,
Ω δ j ν λ i δ i λ i c λ i λ j n 2 2 i , j .
Combining (75)–(78), we see that
Ω u λ i δ i λ i = O j = 1 N 1 ( λ i λ j ) ( n 2 ) / 2 + j i ε i j .
Now, we deal with the last term in the right hand side of (48) with h = λ i δ i λ i . To this aim, we write
Ω | u | p 1 ε u λ i δ i λ i = Ω 2 | u | p 1 ε u λ i δ i λ i + Ω 1 | u | p 1 ε u λ i δ i λ i ,
where Ω 1 , Ω 2 are defined in (57) and u = i = 1 N α i δ i + v : = U + v .
Observe that
Ω 2 | u | p 1 ε u λ i δ i λ i = O Ω 2 | v | p + 1 ε = O v p + 1 ε = o v 2 .
In addition, we write
Ω 1 | u | p 1 ε u λ i δ i λ i = Ω 1 U p ε λ i δ i λ i + p ε Ω 1 U p 1 ε v λ i δ i λ i + O Ω 1 U p 2 ε v 2 δ i = Ω U p ε λ i δ i λ i + p ε Ω U p 1 ε v λ i δ i λ i + O Ω U p 1 ε v 2 + Ω 2 v p + 1 ε .
For the remaining term in (82), we have
Ω U p 1 ε v 2 + Ω 2 v p + 1 ε c Ω U p 1 v 2 + c Ω 2 v p + 1 c v 2 .
For the other terms, from one hand, we write
Ω U p 1 ε v λ i δ i λ i = Ω ( α i δ i ) p 1 ε v λ i δ i λ i + O Ω i ( α i δ i ) p 2 ( j i α j δ j ) | v | | λ i δ i λ i | + O ( Ω i ) c ( j i α j δ j ) p 1 | v | | λ i δ i λ i | ,
where Ω i is defined in (63).
Note that, since n 6 , we have p 1 : = 4 / ( n 2 ) 1 . Thus, it follows that
Ω i ( α i δ i ) p 2 ( j i α j δ j ) | v | | λ i δ i λ i | + ( Ω i ) c ( j i α j δ j ) p 1 | v | | λ i δ i λ i | C Ω i ( α i δ i ) p 1 ( j i α j δ j ) | v | + C ( Ω i ) c ( j i α j δ j ) p 1 | v | δ i C j i Ω ( δ i δ j ) p / 2 | v | c v Ω ( δ i δ j ) n / ( n 2 ) ( n + 2 ) / ( 2 n ) C v ε i j ( n + 2 ) / ( 2 ( n 2 ) ) ln ( ε i j 1 ) ( n + 2 ) / ( 2 n ) .
Now, as the computation of (65), we obtain
Ω δ i p 1 ε λ i δ i λ i v = O ε v + v λ i n 2 2 .
On the other hand, we write
Ω U p ε λ i δ i λ i = Ω ( α i δ i ) p ε λ i δ i λ i + O Ω δ i p ε ( j i α j δ j ) + j i Ω δ j p ε δ i .
For the remaining term in (87), using the fact that ε ln λ i 0 , we see that
Ω δ i p ε ( j i α j δ j ) + j i Ω δ i p ε δ i c j i ε i j .
Now, since ε ln λ i 0 , we notice that
δ i ε ( x ) = c 0 ε λ i ε ( n 2 ) 2 ( 1 + ε n 2 2 ln ( 1 + λ i 2 | x a i | 2 ) ) + O ( ε 2 ln ( 1 + λ i 2 | x a i | 2 ) .
Thus
Ω δ i p ε λ i δ i λ i = c 0 ε λ i ε ( n 2 ) 2 Ω δ i p λ i δ i λ i + ε ( n 2 ) 2 c 0 ε λ i ε ( n 2 ) 2 Ω δ i p λ i δ i λ i ln 1 + λ i 2 | x a i | 2 + O ε 2 δ i p + 1 ln 2 1 + λ i 2 x a i 2 .
However, we have
Ω δ i p δ i λ i = R n δ i p λ i δ i λ i R n \ Ω δ i p λ i δ i λ i = O R n \ Ω δ i p + 1 = O 1 λ i n
and
Ω δ i p λ i δ i λ i ln ( 1 + λ i 2 x a i 2 ) = n 2 2 R n c 0 2 n n 2 ( 1 x 2 ) ( 1 + x 2 ) n + 1 ln ( 1 + x 2 ) d x + O ln λ i λ i n .
This implies that
Ω δ i p ε λ i δ i λ i = c 2 ε c 0 ε λ i ε ( n 2 ) 2 + O ε 2 + ln λ i λ i n .
Combining (80)–(91), we obtain
Ω u p 1 ε u λ i δ i λ i = α i p ε c 2 ε c 0 ε λ i ε ( n 2 ) 2 + O ε 2 + v 2 + j i ε i j + j = 1 N 1 ( λ i λ j ) n 2 2 .
Clearly (74), (79), and (82) give the desired result. □
Next, we provide a balancing condition involving the point of concentration.
Proposition 5. 
Let ( α , λ , a ) O ( N , ν 0 ) , v E a , λ and u = i = 1 N α i δ a i , λ i + v . Then for ε small and for each i N , the following statement holds
I ε ( u ) , 1 λ i δ i a i = α i c 4 λ i 3 K a i + c 3 j i α j λ i ε i j a i 1 c 0 ε α j p 1 ε λ j ε ( n 2 ) 2 c 0 ε α i p 1 ε λ i ε ( n 2 ) 2 + R 3 ,
where
R 3 = O ε 2 + v 2 + k r ε k r n n 2 ln ε k r 1 + j = 1 N 1 λ j n 2 2 λ i n 2 + j i λ j | a i a j | ε i j n + 1 n 2 + T 2 2 ( λ i ) + 1 λ i j i ε i j ln n 2 n ,
c 3 = c 0 2 n n 2 R n d x ( 1 + x 2 ) n + 2 2 , c 4 = ( n 2 ) n c 0 2 R n | x | 2 ( 1 + | x | 2 ) n 1 d x .
Proof. 
As in the proof of the previous proposition, we will take h = 1 λ i δ i a i in (48) and estimate each term.
First, using estimate F 11 of [31], we observe that for i j , we have
Ω δ j 1 λ i δ i a i = Ω Δ δ j 1 λ i δ i a i + Ω δ j ν 1 λ i δ i a i = Ω δ j p 1 λ i δ i a i + O 1 λ j n 2 2 1 λ i n 2 = c 3 1 λ i ε i j a i + O λ j | a i a j | ε i j n + 1 n 2 + O 1 λ j n 2 2 1 λ i n 2
and
Ω δ i 1 λ i δ i a i = R n \ Ω δ i 1 λ i δ i a i = O 1 λ i n 1 .
Combining (93) and (94), we get
Ω u 1 λ i δ i a i = c 3 j i α j 1 λ i ε i j a i + O j i λ j | a i a j | ε i j n + 1 n 2 + j = 1 N 1 λ j n 2 2 1 λ i n 2 .
Second, we notice that
Ω K v 1 λ i δ i a i c Ω | v | δ i c v T 2 λ i , Ω K δ j 1 λ i δ i a i c Ω δ j δ i λ i | x a i | c λ i Ω ( δ i δ j ) n n 2 n 2 n Ω d x | x a | n / 2 2 n
c λ i ε i j ln ( n 2 ) / n ε i j 1 for i j ,
Ω \ B i K δ i 1 λ i δ i a i c λ i n 1
and
B i K δ i 1 λ i δ i a i = B i K ( a i ) ( x a i ) δ i 1 λ i δ i a i + O B i δ i 1 λ i δ i a i ( x a i ) 3 ,
where B i = B ( a i , r ) with r > 0 small.
Note that
B i K x a i a i δ i 1 λ i δ i a i j = l = 1 n K x l a i B i δ i 1 λ i δ i a i j x a i l = K x j a i B i δ i 1 λ i δ i a i j x a i j = K x j a i ( n 2 ) c 0 2 B i λ i n 1 x a i j 2 1 + λ i 2 x a i 2 n 1 d x = K x j a i ( n 2 ) c 0 2 n B i λ i n 1 x a i 2 1 + λ i 2 x a i 2 n 1 d x = K x j a i ( n 2 ) n c 0 2 λ i 3 B ( 0 , λ i r ) | x | 2 1 + x 2 n 1 d x = K x j a i c 4 λ i 3 + O 1 λ i n 1 .
This implies that
B i K a i x a i δ i 1 λ i δ i a i = c 4 λ i 3 K a i + O 1 λ i n 1 .
We also note that, since 1 λ i δ i a i c λ i | x a i | δ i , we have
B i δ i 1 λ i δ i a i | x a i | 3 d x c λ i B i δ i 2 | x a i | 2 d x c λ i B i λ i n 2 x a i 2 1 + λ i 2 x a i 2 n 2 d x c λ i 5 if n 7 λ i 5 ln λ i i f n 6 .
Combining (96)–(101), we obtain
Ω K u 1 λ i δ i a i = α i c 4 λ i 3 K a i + O 1 λ i n 1 + 1 λ i j i ε i j ln n 2 n ε i j 1 + v 2 + T 2 2 ( λ i ) .
Next, we are going to estimate the last term of the right hand side of (48) with h = 1 λ i δ i a i .
  • To this aim, we observe that
Ω 2 | u | p 1 ε u 1 λ i δ i a i Ω 2 | v | p + 1 ε = O v p + 1 ε
and
Ω 1 u p ε 1 λ i δ i a i = Ω 1 U p ε 1 λ i δ i a i + ( p ε ) Ω 1 U p 1 ε v 1 λ i δ i a i + O Ω 1 U p 2 ε v 2 1 λ i δ i a i = Ω U p ε 1 λ i δ i a i + ( p ε ) Ω U p 1 ε v 1 λ i δ i a i + O Ω U p 1 ε v 2 + Ω 2 | v | p + 1 ε ,
where Ω 1 and Ω 2 are defined in (57).
However, we have
Ω U p 1 ε v 2 + Ω 2 | v | p + 1 ε = O ( v 2 )
and
Ω U p ε 1 λ i δ i a i = Ω ( α i δ i ) p ε 1 λ i δ i a i + ( p ε ) Ω ( α i δ i ) p ε 1 ( j i α j δ j ) 1 λ i δ i a i + Ω ( j i α j δ j ) p ε 1 λ i δ i a i + O Ω i ( α i δ i ) p 2 ε ( j i α j δ j ) 2 + ( j i α j δ j ) p ε 1 λ i δ i a i + O Ω \ Ω i α i δ i ( j i α j δ j ) p ε 1 + ( α i δ i ) p 1 ε ( j i α j δ j ) | 1 λ i δ i a i | ,
where Ω i is defined in (63).
We need to estimate each term of the right hand side of (105). First, we observe that
Ω i [ ] 1 λ i δ i a i + Ω \ Ω i [ ] 1 λ i δ i a i C j i Ω ( δ i δ j ) n n 2 C j i ε i j n n 2 ln ε i j 1 .
Second, by oddness we have
Ω ( α i δ i ) p ε 1 λ i δ i a i = O Ω \ B ( a i , r ) δ i p | 1 λ i δ i a i | = O Ω \ B ( a i , r ) δ i p + 1 = O 1 λ i n ,
where r is a fixed small positive constant.
Third, we write
Ω ( α i δ i ) p ε 1 1 λ i δ i a i δ j = c 0 ε λ i ε ( 2 n ) 2 Ω δ i p 1 1 λ i δ i a i δ j + O ε Ω δ i p ε δ j ln ( 1 + λ i 2 | x a i | 2 = c 0 ε λ i ε ( 2 n ) 2 R n δ i p 1 1 λ i δ i a i δ j + O R n \ Ω δ i p δ j + ε Ω ( δ i δ j ) δ i p 1 ln ( 1 + λ i 2 | x a i | 2 ) = c 0 ε λ i ε ( 2 n ) 2 c 3 p λ i ε i j a i + O λ j | a i a j | ε i j n + 1 n 2 + λ i ( n + 2 ) 2 λ j 2 n 2 + ε ε i j ln n 2 n ε i j 1 .
Lastly, we write
Ω ( j i α j δ j ) p ε 1 λ i δ i a i = j i Ω ( α j δ j ) p ε 1 λ i δ i a i + O k { i , j } Ω δ k δ j p 1 | 1 λ i δ i a i | = j i α j p ε c 0 ε λ i ε ( 2 n ) 2 Ω δ j p 1 λ i δ i a i + O ε Ω δ j p | 1 λ i δ i a i | ln ( 1 + λ j 2 | x a j | 2 ) + O k { i , j } Ω δ k δ j p 1 | 1 λ i δ i a i | .
However, we have
k { i , j } Ω δ k δ j p 1 | 1 λ i δ i a i | l r Ω ( δ l δ r ) n / ( n 2 ) C l r ε l r n / ( n 2 ) ln ε l r 1 , Ω δ j p | 1 λ i δ i a i | ln ( 1 + λ j 2 | x a j | 2 ) Ω δ j p 1 ( δ i δ j ) ln ( 1 + λ j 2 | x a j | 2 )
C Ω ( δ i δ j ) n / ( n 2 ) n 2 n C ε i j ln n 2 n ε i j 1 ,
Ω δ j p 1 λ i δ i a i = R n δ j p 1 λ i δ i a i + O R n \ Ω δ j p δ i = c 3 λ i ε i j a i + O λ j | a i a j | ε i j n + 1 n 2 + λ j ( n + 2 ) 2 λ i 2 n 2 .
To deal with the second integral in the right hand side of (103), we write
Ω U p 1 ε v 1 λ i δ i a i = Ω ( α i δ i ) p 1 ε v 1 λ i δ i a i + O Ω i ( α i δ i ) p 2 ε ( j i α j δ j ) | v | | 1 λ i δ i a i | + O Ω \ Ω i ( j i α j δ j ) p 1 ε | v | | 1 λ i δ i a i | ,
where Ω i is defined in (63).
However, since n 6 , we have
Ω i ( α i δ i ) p 2 + ε ( j i α j δ j ) | v | | 1 λ i δ i a i | + Ω i c ( j i α j δ j ) p 1 + ε | v | | 1 λ i δ i a i | C Ω i ( α i δ i ) p 1 ( j i α j δ j ) | v | | + C Ω i c ( j i α j δ j ) p 1 | v | δ i C j i Ω ( δ i δ j ) p 2 | v | c v Ω ( δ i δ j ) n n 2 n + 2 2 n C v ε i j n + 2 2 ( n 2 ) ln n + 2 2 n ε i j 1 .
To complete the estimate (113), we write
Ω δ i p 1 ε v 1 λ i δ i a i = c 0 ε λ i ε ( 2 n ) 2 Ω δ i p 1 v 1 λ i δ i a i + O ε Ω δ i p 1 | v | | 1 λ i δ i a i | ln ( 1 + λ i 2 | x a i | 2 ) = c 0 ε λ i ε ( 2 n ) 2 1 p < 1 λ i δ i a i , v > L 2 Ω ν 1 λ i δ i a i v + O ( ε v ) = O v [ ε + λ i n / 2 ] .
Combining the above estimates, we obtain
Ω | u | p ε u 1 λ i δ i a i = c 3 c 0 ε j i α j 1 λ i ε i j a i k = i , j α k p 1 ε λ k ε n 2 2 + O k r ε k r n n 2 ln ε k r 1 + O j i 1 λ i n + 2 2 λ j n 2 2 + 1 λ i n 2 2 λ j n + 2 2 + v λ i n 2 + ε v + j i λ j | a i a j | ε i j n + 1 n 2 + v 2 + O v T 1 ( a , λ ) + 1 λ i n + ε j i ε i j ln ( n 2 ) / n ε i j 1 .
Combining (95), (102) and (116), we easily obtain the desired result. □

5. Construction of Simple Interior Bubbling Solutions

In this section, we assume that n 6 and we take N m , where m is the number of critical points of K, and let y 1 , ..., y N be non-degenerate distinct critical points of K. As in [32], the strategy of the proof of Theorem 1 is the following: let
M ( N , ε ) = { ( α , λ , a , v ) ( R + ) N × ( R + ) N × Ω N × H 1 ( Ω ) : | α i 1 | < c ε ln 2 ε , 1 c < λ i 2 ε < c , | a i y i | < c ε 1 / 5 1 i N , v E a , λ and v < c ε } ,
where c is a positive constant and E a , λ is defined by (9).
In addition, we consider the following function
ψ ε : M ( N , ε ) R , ( α , λ , a , v ) ψ ε ( α , λ , a , v ) = I ε i = 1 N α i δ a i , λ i + v .
We notice that ( α , λ , a , v ) is a critical point of ψ ε if and only if u = i = 1 N α i δ i + v is a critical point of I ε . Thus, we need to look for critical point of ψ ε . Since the variable v belongs to E a , λ , the lagrange multiplier theorem allows us to get the following proposition.
Proposition 6. 
( α , λ , a , v ) M ( N , ε ) is a critical point of ψ ε if and only if there exists ( A , B , C ) R N × R N × ( R n ) N such that the following holds:
ψ ε = k = 1 N ( A k φ k + B k ϕ k + j = 1 n c k j ξ k j ) ,
where
φ k ( α , λ , a , v ) = Ω v δ k , ϕ k ( α , λ , a , v ) = λ k Ω v δ k λ k , ξ k j ( α , λ , a , v ) = 1 λ k Ω v δ k a k j .
In other words, ( α , λ , a , v ) is a critical point of ψ ε is equivalent to the following system:
ψ ε α i ( α , λ , a , v ) = 0 i 1 , , N
ψ ε λ i ( α , λ , a , v ) = B i Ω v λ i 2 δ i λ i 2 + j = 1 n C i j Ω v 1 λ i 2 δ i λ i a i j i 1 , , N
ψ ε a i ( α , λ , a , v ) = B i Ω v λ i 2 δ i λ i a i + j = 1 n C i j Ω v 1 λ i 2 δ i a i j a i i 1 , , N
ψ ε v ( α , λ , a , v ) = k = 1 N A k φ k v + B k ψ k v + j = 1 n C k j ξ k j v .
  • To prove Theorem 1, we will make a careful study of the previous equations. Notice that
ψ ε α i ( α , λ , a , v ) = I ε i = 1 N α i δ i + v , δ i
ψ ε λ i ( α , λ , a , v ) = I ε i = 1 N α i δ i + v , α i δ i λ i
ψ ε a i ( α , λ , a , v ) = I ε i = 1 N α i δ i + v , α i δ i a i
ψ ε v ( α , λ , a , v ) = I ε i = 1 N α i δ i + v .
Let ( α , λ , a , 0 ) M ( N , ε ) , where M ( N , ε ) is defined by (117). We are going to solve the system defined by Equations (118)–(121). Clearly v ¯ , obtained in Proposition 2, satisfies Equation (121). To simplify the notation, we will write in the sequel v instead of v ¯ . Combining Equations (118)–(125), we see that u = i = 1 N α i δ i + v is a critical point of I ε if and only if ( α , λ , a ) solves the following system for each 1 i N
E α i I ε ( u ) , δ i )   =   0
E λ i I ε ( u ) , α i δ i λ i = B i Ω v λ i 2 δ i λ i 2 + j = 1 n C i j Ω v 1 λ i 2 δ i λ i a i j
E a i I ε ( u ) , α i δ i a i = B i Ω v λ i 2 δ i λ i a i + j = 1 n C i j Ω v 1 λ i 2 δ i a i j a i .
Recall that ( α , λ , a , 0 ) M ( N , ε ) , a i is close to the critical point y i of K and 0 < 1 / c λ i / λ j c for some fixed positive constant. Therefore | a i a j | c > 0 and ε i j = O ( ε ( n 2 ) / 2 ) .
Next, we state the following crucial estimates, which are a direct consequence of Propositions 2–5.
Lemma 2. 
For ε small, the following statements hold:
v c ε 1 i f n 7 , ε | ln ε | 2 3 i f n = 6 , R 1 = O ( ε ) , R 2 ; R 3 c ε 2 i f n 6 ε 2 | ln ε | 4 / 3 i f n = 6 ,
where R 1 , R 2 , and R 3 are defined in Propositions 3, 4, and 5, respectively.
To study the system ( E α i ) , ( E λ i ) , ( E a i ) , we need to estimate the constants A i s , B i s , and C i j s , which appear in equations ( E α i ) , ( E λ i ) , and ( E a i ) . This is the goal of the following lemma:
Lemma 3. 
Let ( α , λ , a , 0 ) M ( N , ε ) . Then, for ε small, the following statements hold:
A i = O ( ε ln 2 ε ) , B i = O ( ε ) a n d C i j = O ( ε 3 / 2 ) 1 i N , 1 j n .
Proof. 
Applying ψ ε v (see (121)) to the functions δ i , λ i δ i λ i and 1 λ i δ i a i j , we obtained the following quasi-diagonal system
c A i + O ε n 2 2 k = 1 N ( | A k | + | B k | + j = 1 n | C k j ) = ψ ε v , δ i c B i + O ε n 2 2 k = 1 N ( | A k | + | B k | + j = 1 n | C k j ) = ψ ε v , λ i δ i λ i c C i l + O ε n 2 2 k = 1 N ( | A k | + | B k | + j = 1 n | C k j ) = ψ ε v , 1 λ i δ i a i l ,
where
c = R n | δ i | 2 , c = R n λ i δ i λ i 2 and c = R n 1 λ i δ i a i j 2 .
Combining Propositions 3–5, estimates (29)–(37) and the fact that ( α , λ , a , 0 ) M ( N , ε ) , we see that for all i { 1 , , N } we have
ψ ε v , δ i = O ( ε ln 2 ε ) , ψ ε v , λ i δ i λ i = O ( ε ) and ψ ε v , 1 λ i δ i a i = O ( ε 3 2 ) .
This implies that
M ( A , B , C ) T = ( O ( ε ln 2 ε ) , O ( ε ) , O ( ε 3 / 2 ) T ,
where M is the matrix defined by
m i j = O ( ε n 2 2 ) i j and m i i = c + O ( ε n 2 2 ) for 1 i N c + O ( ε n 2 2 ) for N + 1 i 2 N c + O ( ε n 2 2 ) for 2 N + 1 i N ( n + 2 ) ,
where c, c , and c are defined in (129).
Clearly, M is an invertible matrix. We see that M and M 1 have the same form. Thus, Lemma 3 follows. □
Next, we are going to analyse equations ( E α i ) , ( E λ i ) , ( E a i ) . To obtain an easy system to solve, we considered the following change of variables
β i = 1 α i p 1 , 1 λ i 2 = c 2 c 1 1 K ( y i ) ε ( 1 + i ) , z i = a i y i 1 i N ,
where c 1 and c 2 are defined in Proposition 4.
This change of variables allows us to rewrite the system in the following simple form:
Lemma 4. 
For ε small, equations ( E α i ) , ( E λ i ) , ( E a i ) are equivalent to the following system
( S ) β i = O ( ε | ln ε | ) 1 i N i = O ( | β i | + | z i | 2 + ε 1 / 2 ) 1 i N D 2 K ( y i ) ( z i , . ) = O ε 1 / 2 + | z i | 2 1 i N .
Proof. 
Using the fact that
α i p 1 ε = α i p 1 + O ( ε ) and λ i ε n 2 2 = 1 + O ( ε ln λ i ) ,
we see that ( E α i ) is equivalent to
( E α i ) β i = O ( ε | ln ε | ) 1 i N .
For the second equation ( E λ i ) , Proposition 4, Lemmas 2 and 3 imply that
( 1 β i ) c 2 ε c 1 K ( a i ) λ i 2 c ε min ( 2 , ( n 2 ) / 2 ) if n 6 ε 2 | ln ε | 4 / 3 if n = 6 .
Writing
K ( a i ) = K ( y i ) + O ( | z i | 2 ) ,
we obtain the second equation in the system ( S ) .
Lastly, writing
K ( a i ) = D 2 K ( y i ) ( z i , . ) + O ( | z i | 2 )
and using Proposition 5, Lemmas 2 and 3, we see that equation ( E a i ) is equivalent to the third equation in the system ( S ) , which completes the proof of Lemma 4. □
Now, we are ready to prove our results related to the construction of simple interior bubbling solutions. We note that, since Theorem 2 is a straightforward consequence of Theorem 1, we only need to prove this result.
Proof of Theorem 1 
The system ( S ) , given in Lemma 4, can be rewritten
β i = U 1 , i ( ε , β , , z ) = O ( ε | ln ε | ) 1 i N i = U 2 , i ( ε , β , , z ) = O ( | β i | + | z i | 2 + ε 1 / 2 ) 1 i N D 2 K ( y i ) ( z i , . ) = U 3 , i ( ε , β , , z ) = O ε 1 / 2 + | z i | 2 1 i N ,
where β = ( β 1 , , β N ) , = ( 1 , , N ) and z = ( z 1 , , z N ) .
Defining the following linear map
l : R N × R N × ( R n ) N R N × R N × ( R n ) N ( β , , z ) ( β , , D 2 K ( y 1 ) ( z 1 , . ) , , D 2 K ( y N ) ( z N , . ) ) ,
we see that the system is equivalent to
l ( β , , z ) = ( U 1 ( ε , β , , z ) , U 2 ( ε , β , , z ) , U 3 ( ε , β , , z ) ) ,
where
U 1 ( ε , β , , z ) = ( U 1 , 1 , , U 1 , N ) , U 2 ( ε , β , , z ) = ( U 2 , 1 , , U 2 , N ) , U 3 ( ε , β , , z ) = ( U 3 , 1 , , U 3 , N ) .
In addition, since y i is a non-degenerate critical point of K, we deduced that l is invertible. This implies that (131) is equivalent to
( β , , z ) = l 1 ( U 1 ( ε , β , , z ) , U 2 ( ε , β , , z ) , U 3 ( ε , β , , z ) ) .
Choosing r positive small and ( β , , z ) B ( 0 , r ) ¯ , we obtain
l 1 ( U 1 ( ε , β , , z ) ,   U 2 ( ε , β , , z ) , U 3 ( ε , β , , z ) ) l 1 ( ε , β , , z ) , U 2 ( ε , β , , z ) , U 3 ( ε , β , , z ) ) C l 1 ( | z | 2 + | e 1 / 4 ) C l 1 2 r 2
and hence if we choose r 1 / ( 2 C l 1 , we see that the function
f : B ( 0 , r ) ¯ B ( 0 , r ) ¯ , f ( β , , z ) = l 1 ( U 1 ( ε , β , , z ) , U 2 ( ε , β , , z ) , U 3 ( ε , β , , z ) )
is well-defined and continuous. Thus, applying Brouwer’s fixed point theorem, we derived that f has a fixed point. Therefore, the system ( S ) has at least one solution ( β ε , ε , z ε ) for ε small. To complete the proof of the theorem, it remains to be proven that the constructed function u ε = i = 1 N α i δ i + v is positive. To this aim, we first remark, since v 0 , that u ε 0 for ε small. By construction, u ε satisfies
( P ε ) : Δ u ε + K ( x ) u ε = | u ε | p 1 ε u ε in Ω u ε n = 0 on Ω .
Multiplying ( P ε ) by u ε = max ( 0 , u ε ) and integrating on Ω, we obtain
Ω | u ε | 2 + Ω K ( u ε ) 2 = Ω ( u ε ) p + 1 ε .
However, we have
Ω ( u ε ) p + 1 ε c Ω ( u ε ) p + 1 p + 1 ε p + 1 C Ω | u ε | 2 + Ω K ( u ε ) 2 p + 1 ε 2
which implies that
either u ε 0 or Ω ( u ε ) p + 1 C 2 ( n 2 ) / ( 4 ε ( n 2 ) ) .
However, since u ε | v ε | and v ε L p + 1 0 as ε 0 , we derive that u ε 0 and u ε 0 . Thus, using the maximum principal, u ε has to be positive. This completes the proof of the theorem. □

6. Construction of Clustered Bubbling Solutions

This section is devoted to the proof of Theorems 3 and 4. We start by proving Theorem 3. Let n 6 , y be a non-degenerate critical point of K and ( z ¯ 1 , , z ¯ N ) be a non-degenerate critical point of F N , y , where F N , y is defined by (3). The strategy of the proof of Theorem 3 is the same as that of Theorem 1. We start by introducing a neighborhood of the desired constructed solutions. Let
V ( N , y , ε ) = { ( α , λ , a , v ) ( R + ) N × ( R + ) N × Ω N × H 1 ( Ω ) : | α i 1 | < c ε ln 2 ε , 1 c < λ i 2 ε < c , | a i y ε γ σ z ¯ i | ε γ η 0 1 i N , v E a , λ and v < c ε } ,
where E a , λ is defined by (9),
γ = n 4 2 n and σ = c 3 c 4 1 / n c 2 c 1 K ( y ) γ
with c 3 , c 4 , c 2 and c 1 are the constants defined in Propositions 4 and 5. As in the proof of Theorem 1, we reduced the problem to a finite dimensional system. Proposition 2 allows us to obtain such a reduction by finding v ¯ verifying the Equation (121). Thus, we are looking for ( α , λ , a ) V N , y , ε solution of the system defined by Equations (126)–(128), where V N , y , ε is defined by
V N , y , ε = { ( α , λ , a ) ( R + ) N × ( R + ) N × Ω N : ( α , λ , a , 0 ) V ( N , y , ε ) 1 i N } .
As in the proof of Theorem 1 and in order to work with a simpler system, we made the following change of variables:
β i = 1 α i p 1 , 1 λ i 2 = c 2 c 1 1 K ( y ) ε ( 1 + i ) , ( a i y ) = ε γ σ ( ζ i + z ¯ i ) 1 i N ,
where c 1 and c 2 are defined in Proposition 4.
Using these changes of variables, it is easy to see that
ε i j = 1 ( λ i λ j | a i a j | 2 ) n 2 2 ( 1 + O ( ε 4 / n ) ) = O ε 2 γ + 1
ε i j a i = ( n 2 ) ε γ + 1 c 2 n 2 2 σ n 1 ( c 1 K ( y ) ) n 2 2 ( ζ j + z ¯ j ζ i z ¯ i ) ζ i + z ¯ i ζ j z ¯ j n 1 + R ( i , j ( 1 + O ( ε 4 n ) ) = O ( ε γ + 1 ) ,
where
R ( i , j ) = n 2 4 i + n 2 4 j + O i 2 + O j 2 .
Next, using Propositions 2–5, we deduced that the following estimates hold:
Lemma 5. 
For ε small, the following statements hold:
v c ε i f n 7 , ε | ln ε | 2 3 i f n = 6 , R 1 = O ( ε ) , R 3 c ε 2 | ln ε | i f n 7 ε 2 | ln ε | 4 / 3 i f n = 6 , R 2 = O ε 2 γ + 1 ,
where R 1 , R 2 , and R 3 are defined in Propositions 3, 4, and 5, respectively.
Now, arguing as in the proof of Lemma 3, we derive that the constants A i s , B i s , and C i j s , which appear in equations ( E α i ) , ( E λ i ) and ( E a i ) , satisfy the following estimates:
Lemma 6. 
Let ( α , λ , a , 0 ) V ( N , y , ε ) . Then, for ε small, the following statements hold:
A i = O ( ε ln 2 ε ) , B i = O ( ε ) a n d C i j = O ( ε γ + 3 / 2 ) 1 i N , 1 j n .
Next, our aim is to rewrite equations ( E α i ) , ( E λ i ) , ( E a i ) in a simple form.
Lemma 7. 
For ε small, equations ( E α i ) , ( E λ i ) , ( E a i ) are equivalent to the following system
( S ) β i = O ( ε | ln ε | ) 1 i N i = O | β i | + ε 2 γ 1 i N D 2 K ( y ) ( ζ i , . ) ( n 2 ) j i ζ j ζ i | z ¯ j z ¯ i | n + n ( n 2 ) j i z ¯ j z ¯ i | z ¯ j z ¯ i | 2 , ζ j ζ i ( n 2 ) j i z ¯ j z ¯ i | z ¯ j z ¯ i | n n 6 4 i + n 2 4 j = O ( R 4 ) 1 i N ,
where
R 4 = j = 1 N j 2 + j = 1 N | ζ j | 2 + ε γ + ε 1 2 γ | ln ε | 4 / 3 .
Proof. 
First, using Proposition 3, Lemma 5 and the fact that ( α , λ , a ) V N , y , ε , we see that ( E α i ) is equivalent to
( F α i ) β i = O ( ε | ln ε | ) 1 i N .
Second, using Lemma 6, we write
I ε ( u ) , λ i δ i λ i = O | B i | v L 2 λ i 2 δ i λ i 2 L 2 + j = 1 N | C i j | v L 2 2 δ i λ i a i j L 2 = O ε v λ i + ε γ + 3 / 2 v = O ε 3 / 2 v .
Using Proposition 4 and Lemma 5, we obtain
c 2 ε c 1 K ( a i ) λ i 2 = O ε | β i | + ε 2 γ + 1 .
However, we have
K ( a i ) = K ( y ) + O | a i y | 2 = K ( y ) + O ε 2 γ .
This implies that ( E λ i ) is equivalent to
( F λ i ) i = O | β i | + ε 1 4 / n 1 i N .
To deal with the third equation ( E a i ) , we write
I ε ( u ) , 1 λ i δ i a i = O | B i | v L 2 2 δ i λ i a i L 2 + j = 1 N | C i j | v L 2 1 λ i 2 2 δ i a i a i j L 2 = O ε v + ε γ + 3 / 2 v = O ε v .
However, combining (139), Proposition 5 and Lemma 5, (141) becomes
c 4 λ i 3 K ( a i ) c 5 ε 1 + γ λ i j i ( ζ j + z ¯ j ζ i z ¯ i ) | ζ j + z ¯ j ζ i z ¯ i | n 1 + n 2 4 i + n 2 4 j + O i 2 + j 2 = O ε 2 | ln ε | 4 / 3 ,
where
c 5 = c 3 ( n 2 ) σ n 1 c 2 c 1 K ( y ) n 2 2 .
Observe that
K ( a i ) = D 2 K ( y ) ( a i y , . ) + O | a i y | 2 = ε γ σ D 2 K ( y ) ( ζ i + z ¯ i , . ) + O ε 2 γ .
Combining (142) and (143), we obtain
D 2 K ( y ) ( ζ i + z ¯ i , . ) ( n 2 ) j i ( ζ j + z ¯ j ζ i z ¯ i ) | ζ j + z ¯ j ζ i z ¯ i | n 1 + n 6 4 i + n 2 4 j + O i 2 + j 2 = O ε γ + ε 1 2 γ | ln ε | 4 3 = O ε 1 2 γ | ln ε | 4 3 .
However, we have
D 2 K ( y ) ( ζ i + z ¯ i , . ) = D 2 K ( y ) ( ζ i , . ) + D 2 K ( y ) ( z ¯ i , . )
and
( ζ j + z ¯ j ζ i z ¯ i ) | ζ j + z ¯ j ζ i z ¯ i | n = ( z ¯ j z ¯ i ) | z ¯ j z ¯ i | n + ( ζ j ζ i ) | z ¯ j z ¯ i | n 1 n ( z ¯ j z ¯ i ) | z ¯ j z ¯ i | 2 , ζ j ζ i + O j = 1 N | ζ j | 2 = ( z ¯ j z ¯ i ) | z ¯ j z ¯ i | n n ( z ¯ j z ¯ i ) | z ¯ j z ¯ i | 2 , ζ j ζ i ( z ¯ j z ¯ i ) | z ¯ j z ¯ i | n + ( ζ j ζ i ) | z ¯ j z ¯ i | n + O j = 1 N | ζ j | 2 .
Combining (144)–(146) and the fact that ( z ¯ 1 , , z ¯ N ) is a critical point of F N , y , we see that equation ( E a i ) is equivalent to
( F a i ) D 2 K ( y ) ( ζ i , . ) ( n 2 ) j i ζ j ζ i | z ¯ j z ¯ i | n + n ( n 2 ) j i z ¯ j z ¯ i | z ¯ j z ¯ i | 2 , ζ j ζ i ( n 2 ) j i z ¯ j z ¯ i | z ¯ j z ¯ i | n n 6 4 i + n 2 4 j = O ( R 4 ) ,
which completes the proof of Lemma 7. □
Now, we are ready to prove our results related to the construction of clustered bubbling solutions.
Proof of Theorem 3 
Note that the system ( F a 1 ) , ..., ( F a N ) is equivalent to
1 2 D 2 F N , y ( z ¯ 1 , , z ¯ N ) ( ζ 1 , , ζ N ) ( n 2 ) j i z ¯ j z ¯ i | z ¯ j z ¯ i | n n 6 4 i + n 2 4 j = O j = 1 N j 2 + j = 1 N | ζ j | 2 + ε γ + ε 1 2 γ | ln ε | 4 / 3 .
As in the proof of Theorem 1, we defined a linear map by taking the left hand side of the system defined by ( F α i ) , ( F λ i ) , and ( F a i ) . Since ( z ¯ 1 , , z ¯ N ) is a non-degenerate critical point of F N , y , we deduced that such a linear map is invertible and arguing as in the proof of Theorem 1, we derive that the system ( S ) has a solution ( β ε , ε , ζ ε ) for ε small. This implies that ( P ε ) admits a solution u ε , y = i = 1 N α i , ε δ a i , ε , λ i , ε + v ε . and by construction, properties (4)–(7) are satisfied. The proof of Theorem 3 is thereby completed. □
Proof of Theorem 4 
For y 1 , y 2 two non-degenerate critical points of K, we remark that, for a 1 , ε , a 2 , ε such that | a 1 , ε y 1 | < < 1 and | a 2 , ε y 2 | < < 1 , we see that the interaction between the bubbles δ a 1 , ε , λ 1 , ε and δ a 2 , ε , λ 1 , ε is of the order of 1 ( λ 1 , ε λ 2 , ε ) ( n 2 ) / 2 = O ( ε ( n 2 ) ) / 2 ) which is a negligible term in front of ε 3 / 2 for n 6 . This implies that the interaction goes into the remainder term and thus, we may separate each pack alone. Hence, arguing as in the proof of Theorem 3 and taking a new system ( ( S 1 ) , ..., ( S m ) ) with each ( S i ) represents the system studied in the proof of Theorem 3, and the proof of the theorem follows. □

Author Contributions

K.E.M. and F.M.S.: conceptualization, methodology, investigation, writingoriginal draft, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

Researchers would like to thank the Deanship of Scientific Research, Qassim University for funding publication of this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gierer, A.; Meinhardt, H. A theory of biological pattern formation. Kybernetik 1972, 12, 30–39. [Google Scholar] [CrossRef] [Green Version]
  2. Keller, E.F.; Segel, L.A. Initiation of slime mold aggregation viewed as an instability. J. Theor. Biol. 1970, 26, 399–415. [Google Scholar] [CrossRef]
  3. Lin, C.S.; Ni, W.N.; Takagi, I. Large amplitude stationary solutions to a chemotaxis system. J. Differ. Equ. 1988, 72, 1–27. [Google Scholar] [CrossRef] [Green Version]
  4. Ni, W.N.; Takagi, I. On the shape of least-energy solutions to a semilinear Neumann problem. Commun. Pure Appl. Math. 1991, 44, 819–851. [Google Scholar] [CrossRef]
  5. Ni, W.-M.; Takagi, I. Locating the peaks of least-energy solutions to a semi-linear Neumann problem. Duke Math. J. 1993, 70, 247–281. [Google Scholar] [CrossRef]
  6. Ni, W.-M. Diffusion, cross-diffusion, and their spike-layer steady states. Not. Am. Math. Soc. 1998, 45, 9–18. [Google Scholar]
  7. Zhu, M. Uniqueness results through a priori estimates, I. A three dimensional Neumann problem. J. Differ. Equ. 1999, 154, 284–317. [Google Scholar] [CrossRef] [Green Version]
  8. Wei, J.; Xu, B.; Yang, W. On Lin-Ni’s Conjecture in dimensions four and six. Sci. China Math. 2019, 49, 281–306. [Google Scholar] [CrossRef]
  9. Adimurthi; Yadava, S.L. Existence and nonexistence of positive radial solutions of Neumann problems with critical Sobolev exponents. Arch. Ration. Mech. Anal. 1991, 115, 275–296. [Google Scholar] [CrossRef]
  10. Rey, O.; Wei, J. Arbitrary number of positive solutions for elliptic problem with critical nonlinearity. J. Eur. Math. Soc. 2005, 7, 449–476. [Google Scholar] [CrossRef] [Green Version]
  11. Adimurthi; Mancini, G. The Neumann problem for elliptic equations with critical nonlinearity, “A tribute in honour of G. Prodi”. Sc. Norm. Super. Pisa 1991, 9–25. [Google Scholar]
  12. Wang, X.J. Neumann problem of semilinear elliptic equations involving critical Sobolev exponents. J. Differ. Equ. 1991, 93, 283–310. [Google Scholar] [CrossRef] [Green Version]
  13. Adimurthi; Mancini, G. Geometry and topology of the boundary in the critical Neumann problem. J. Reine Angew. Math. 1994, 456, 1–18. [Google Scholar]
  14. Adimurthi; Pacella, F.; Yadava, S.L. Interaction between the geometry of the boundary and positive solutions of a semilinear Neumann problem with critical nonlinearity. J. Funct. Anal. 1993, 113, 318–350. [Google Scholar] [CrossRef] [Green Version]
  15. Gui, C.; Lin, C.S. Estimates for boundary-bubbling solutions to an elliptic Neumann problem. J. Reine Angew. Math. 2002, 546, 201–235. [Google Scholar] [CrossRef]
  16. Lin, C.S.; Wang, L.; Wei, J. Bubble Accumulations in an Elliptic Neumann Problem with Critical Sobolev Exponent. Calc. Var. Partial. Differ. Equ. 2007, 30, 153–182. [Google Scholar] [CrossRef]
  17. Maier-Paape, S.; Schmitt, K.; Wang, Z.Q. On Neumann problems for semilinear elliptic equations with critical nonlinearity: Existence and symmetry of multi-peaked solutions. Commun. Partial. Differ. Equ. 1997, 22, 1493–1527. [Google Scholar]
  18. Ni, W.-M.; Pan, X.B.; Takagi, I. Singular behavior of least-energy solutions of a semi-linear Neumann problem involving critical Sobolev exponents. Duke Math. J. 1992, 67, 1–20. [Google Scholar] [CrossRef]
  19. Rey, O. An elliptic Neumann problem with critical nonlinearity in three dimensional domains. Commun. Contemp. Math. 1999, 1, 405–449. [Google Scholar] [CrossRef]
  20. Wang, Z.Q. The effect of domain geometry on the number of positive solutions of Neumann problems with critical exponents. Differ. Integral Equ. 1995, 8, 1533–1554. [Google Scholar] [CrossRef]
  21. Wang, Z.Q. High energy and multi-peaked solutions for a nonlinear Neumann problem with critical exponent. Proc. R. Soc. Edinb. Sect. A 1995, 125, 1003–1029. [Google Scholar] [CrossRef]
  22. Wang, Z.Q. Construction of multi-peaked solutions for a nonlinear Neumann problem with critical exponent in symmetric domains. Nonlinear Anal. Theory Methods Appl. 1996, 27, 1281–1306. [Google Scholar] [CrossRef]
  23. Rey, O. The question of interior blow-up points for an elliptic Neumann problem: The critical case. J. Math. Pures Appl. 2002, 81, 655–696. [Google Scholar] [CrossRef] [Green Version]
  24. Rey, O.; Wei, J. Blow-up solutions for an elliptic Neumann problem with sub- or supercritical nonlinearity, I: N = 3. J. Funct. Anal. 2004, 212, 472–499. [Google Scholar] [CrossRef] [Green Version]
  25. Rey, O.; Wei, J. Blow-up solutions for an elliptic Neumann problem with sub- or supercritical nonlinearity, II: N ≥ 4. Ann. L’i.H.P. Anal. Non Linéaire 2005, 22, 459–484. [Google Scholar] [CrossRef]
  26. Wang, X.; Wei, J. On the equation Δu + K(x)u = u n + 2 n 2 ± ϵ 2 in Rn. Rend. Circ. Mat. Palermo 1995, 2, 365–400. [Google Scholar] [CrossRef]
  27. Wang, X. On location of blow-up of ground states of semilinear elliptic equations in Rn involving critical Sobolev exponents. J. Differ. Equ. 1996, 127, 148–173. [Google Scholar] [CrossRef] [Green Version]
  28. Pan, X. Blow-up behavior of ground states of semilinear elliptic equations in R n involving critical Sobolev exponents. J. Differ. Equ. 1992, 99, 78–107. [Google Scholar] [CrossRef] [Green Version]
  29. Pistoia, A.; Weith, T. Sign-changing bubble-tower solutions in a slightly subcritical semi-linear Dirichlet problem. Ann. L’i.H.P. Anal. Non Linéaire 2007, 24, 325–340. [Google Scholar]
  30. Caffarelli, L.; Gidas, B.; Spruck, J. Asymptotic symmetry and local behavior of semilinear elliptic equations with critical Sobolev growth. Commun. Pure Appl. Math. 1989, 42, 271–297. [Google Scholar] [CrossRef]
  31. Bahri, A. Critical Points at Infinity in Some Variational Problems; Research Notes in Mathematics, 182; Longman-Pitman: London, UK, 1989. [Google Scholar]
  32. Bahri, A.; Li, Y.; Rey, O. On a variational problem with lack of compactness: The topological effect of the critical points at infinity. Calc. Var. Partial. Differ. Equ. 1995, 3, 67–93. [Google Scholar] [CrossRef]
  33. Rey, O. Boundary effect for an elliptic Neumann problem with critical nonlinearity. Commun. Partial. Differ. Equ. 1997, 22, 1055–1139. [Google Scholar] [CrossRef]
  34. Dammak, Y. A non-existence result for low energy sign-changing solutions of the Brezis–Nirenberg problem in dimensions 4, 5 and 6. J. Differ. Equ. 2017, 263, 7559–7600. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

El Mehdi, K.; Mohamed Salem, F. Interior Bubbling Solutions for an Elliptic Equation with Slightly Subcritical Nonlinearity. Mathematics 2023, 11, 1471. https://0-doi-org.brum.beds.ac.uk/10.3390/math11061471

AMA Style

El Mehdi K, Mohamed Salem F. Interior Bubbling Solutions for an Elliptic Equation with Slightly Subcritical Nonlinearity. Mathematics. 2023; 11(6):1471. https://0-doi-org.brum.beds.ac.uk/10.3390/math11061471

Chicago/Turabian Style

El Mehdi, Khalil, and Fatimetou Mohamed Salem. 2023. "Interior Bubbling Solutions for an Elliptic Equation with Slightly Subcritical Nonlinearity" Mathematics 11, no. 6: 1471. https://0-doi-org.brum.beds.ac.uk/10.3390/math11061471

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop