Skip to main content

Multiple blowing-up and concentrating solutions for Liouville-type equations with singular sources under mixed boundary conditions

Abstract

In this article, we mainly construct multiple blowing-up and concentrating solutions for a class of Liouville-type equations under mixed boundary conditions:

- Δ v = ε 2 e v - 4 π i = 1 N α i δ p i , in Ω , ε ( 1 - t ) v ν + t b ( x ) v = 0 , on Ω ,

for ε small, where t ( 0 , 1 ] ,N { 0 } , { α 1 , α 2 , , α N } ( - 1 , + ) \ ( { 0 } ) , Ω is a bounded, smooth domain in 2 , Γ := {p1, ..., p N } Ω is the set of singular sources, δ p denotes the Dirac mass at p, ν denotes unit outward normal vector to Ω and b(x) > 0 is a smooth function on Ω.

2000 Mathematics Subject Classification: 35B25; 35J25; 35B38.

1 Introduction

In this article, we mainly investigate the mixed boundary value problem:

- Δ v = ε 2 e v - 4 π i = 1 N α i δ p i , in Ω , ε ( 1 - t ) v ν + t b ( x ) v = 0 , on Ω ,
(1)

for ε small, where t ( 0 , 1 ] ,N { 0 } , { α 1 , α 2 , , α N } ( - 1 , + ) \ ( { 0 } ) , Ω is a bounded, smooth domain in 2 , Γ:= {p1, ..., p N } Ω is the set of singular sources, δ p denotes the Dirac mass at p, ν denotes unit outward normal vector to Ω and b(x) > 0 is a smooth function on Ω.

Such problems occur in conformal geometry [1], statistical mechanics [24], Chern-Simons vortex theory [511] and several other fields of applied mathematics [1216]. In all these contexts, an interesting point is how to construct solutions which exactly "blow-up" and "concentrate" at some given points, whose location carries relevant information about the potentially geometrical or physical properties of the problem. However, the authors mainly consider the Dirichlet boundary value problem, and little is known for the problem with singular sources satisfying α i (-1, 0) for some i = 1, ..., N. The main purpose of this article is to study how to construct multiple blowing-up and concentrating solutions of the Equation (1) with the mixed boundary conditions and singular sources.

Let G t,ε denotes the Green's function of -∆ with mixed boundary conditions on Ω, namely for any y Ω,

- Δ x G t , ε ( x , y ) = 2 π δ y ( x ) , in Ω , ε ( 1 - t ) G t , ε ( x , y ) ν + t b ( x ) G t , ε ( x , y ) = 0 , on Ω ,
(2)

and let H t,ε (x, y) = G t,ε (x, y) + log |x - y| be its regular part. Set G1 = G1,εand H1 = H1,ε. Since ε exactly disappears in the Equation (2)|t = 1, G1 and H1 don't depend on ε. The Equation (1) is equivalent to solving for u=v+2 i = 1 N α i G t , ε ( x , p i ) , the regular part of v, the equation

- Δ u = ε 2 x - p 1 2 α 1 x - p N 2 α N e - 2 i = 1 N α i H t , ε ( x , p i ) e u , in Ω , ε ( 1 - t ) u ν + t b ( x ) u = 0 , on Ω .

Thus, we consider the more general model problem:

- Δ u = ε 2 x - p 1 2 α 1 x - p N 2 α N f ( x ) e u , in Ω , ε ( 1 - t ) u ν + t b ( x ) u = 0 , on Ω ,
(3)

where f:Ω is a smooth function such that f(pi) > 0 for any i = 1, ..., N. Set Ω' = {x Ω: f(x) > 0}, S ( x ) = x - p 1 2 α 1 x - p N 2 α N and Δ m = {p = (p1, ..., p m ) Ωm: p i = pj for some i ≠ j}.

It is known that for { α 1 , , α N } ( 0 , + ) \, or α i = 0 for any i = 1, ..., N, if u ε is a family of solutions of the Equation (3)|t = 1with inf Ω f > 0, which is not uniformly bounded from above for ε small, then u ε blows up at different points p k 1 ,, p k n + m with n + m ≥ 1, 0 ≤ nN, p= ( p k n + 1 , , p k n + m ) ( Ω \ Γ ) m \ Δ m and { p k 1 , , p k n } Γ, and satisfies the concentration property:

ε 2 x - p 1 2 α 1 x - p N 2 α N f ( x ) e u ε 8 π i = 1 n ( 1 + α k i ) δ p k i + 8 π i = n + 1 n + m δ p k i ,
(4)

in the sense of measures in Ω ̄ . Moreover, p= ( p k n + 1 , , p k n + m ) is a critical point of the function:

φ n , m ( p ) = i = n + 1 n + m H 1 ( p k i , p k i ) + 1 2 log f ( p k i ) S ( p k i ) + i , j = n + 1 , i j n + m G 1 ( p k j , p k i ) + 2 i = 1 n j = n + 1 n + m ( 1 + α k i ) G 1 ( p k j , p k i ) ,
(5)

(see [7, 1723]). An obvious problem for the Equation (3) is the reciprocal, namely the existence of multiple blowing-up solutions with concentration points near critical points of φ n,m .

The earlier result concerning the existence of multiple blowing-up and concentrating solutions of the Equation (3) is given by Baraket and Pacard in [24]. When t = 1 and α i = 0 for any i = 1,2, ..., N, they prove that any non-degenerate critical point p= ( p k n + 1 , , p k n + m ) of the function φ n,m with n = 0 generates a family of the solutions u ε which blow-up at p k n + 1 ,, p k n + m , and concentrate in the sense that (4) holds. Esposito [20] performs a similar asymptotic analysis and extends the previous result by allowing the presence of singular sources in the Equation (3)|t = 1, that is, { α 1 , , α N } ( 0 , + ) \. However, the asymptotic analysis method depends on the non-degenerate assumption of critical point of the function φ n,m so much that it pays in return at a price of the very complicated and accurate control on the asymptotics of the solutions.

In fact, the finite dimensional reduction method, used successfully in higher dimensional nonlinear elliptic equation involving critical Sobolev exponent (see [6, 25]), can avoid the technical difficulty in carrying out the asymptotic analysis method for the Equation (3). It is necessary to point out that the key step of the finite dimensional reduction is the analysis of the bounded invertibility of the corresponding linearized operator L of the Equation (3) at the suitable approximate solution. In [26, 27], the authors construct the approximate solution, carry out the finite dimensional reduction and use some stability assumptions of critical points of φ0,mto get the existence of multiple blowing-up and concentrating solutions for the Equation (3)|t = 1with Γ=, namely α i = 0 for any i = 1,2, ..., N. When { α 1 , , α N } ( 0 , + ) \, a similar result for the Equation (3)|t = 1under C0-stable assumption of critical point of φ n,m (see Definition 4.1) is also established in [28].

Here in the spirit of the finite dimensional reduction, we try to extend the result of the Equation (1) in [20, 28] by allowing the presence of singular sources 4 π i = 1 N α i δ p i with some α i (-1, 0) and Robin boundary conditions ε ( 1 - t ) v ν + t b ( x ) v = 0 with t (0,1). When we carry out the finite dimensional reduction, we need to get the invertibility of the desired linearized operator L for the Equation (3) under some α i (-1, 0). Obviously, the linearized operator L easily produces the singularities at some singular sources with α i (-1, 0), which makes trouble for the analysis of the bounded invertibility of L. But we can successfully get rid of it by introducing a suitable L-weighted norm (see (30) below) related with a "gap interval" (- 1, α0), where α0 = min{0, α1, ..., α N }. On the other hand, the presence of the term ε ( 1 - t ) u ν in the Equation (3)|0<t<1brings some new technical difficulties. A flexible approach exactly helps us overcome the difficulties by making use of the maximum principle. In addition, a weaker stable assumption of critical points of the function φ n,m also helps us construct multiple blowing-up and concentrating solutions of the Equation (3). As a consequence, we have the following result.

Theorem 1.1 Let 0 ≤ n ≤ N and m { 0 } such that n + m ≥ 1. Assume that N { 0 } , { α 1 , α 2 , , α N } ( - 1 , + ) \ ( { 0 } ) and p * = ( p n + 1 * , , p n + m * ) is a C0-stable critical point for φ n,m in (Ω' \ Γ)m \ Δ m with m ≥ 1 (see Definition 4.1). Then there exists a family of solutions u ε for the Equation (3) with the concentration property (4), which blow up at n-different points ( p k 1 , , p k n ) in Γ, and m-points p= ( p k n + 1 , , p k n + m ) in (Ω' \ Γ)m \ Δ m with φ n,m (p*) = φ n,m (p). Moreover, u ε remains uniformly bounded on Ω\ i = 1 n + m B λ ( p k i ) , and sup B λ ( p k i ) u ε + for any λ > 0.

Let us point out that from the proof of Theorem 1.1 Robin boundary condition can be considered as a perturbation of Dirichlet boundary condition for the problem (3) in using perturbation techniques to construct multiple blowing-up and concentrating solutions. Based on this point, we also consider the Dirichlet-Robin boundary value problem:

- Δ u = ε 2 x - p 1 2 α 1 x - p N 2 α N f ( x ) e u , in Ω , ε u ν + b ( x ) u = 0 , on T , u = 0 , on Ω \ T ,
(6)

where T Ω is a relatively closed subset and Ω\T. This together with other similar mixed boundary value problems can be founded in [29, 30]. For the problem (6), we obtain the following result.

Theorem 1.2 Under the assumption of Theorem 1.1, then there exists a family of solutions u ε for the Equation (6) with the concentration property (4), which blow up at n-different points ( p k 1 , , p k n ) in Γ, and m-points p= ( p k n + 1 , , p k n + m ) in (Ω' \ Γ)m\ Δ m with φn,m(p*) = φn,m(p). Moreover, u ε remains uniformly bounded on Ω\ i = 1 n + m B λ ( p k i ) , and sup B λ ( p k i ) u ε + for any λ > 0.

Finally, it is very interesting to mention that to prove the above results we need to choose the classification solutions of the following Liouville-type equation to construct concentrating solutions of the Equation (1) or (3):

- Δ u = z 2 γ e u , in 2 , 2 z 2 γ e u < + , γ > - 1 ,
(7)

given by

u ( z ) = log 8 ( 1 + γ ) 2 μ 2 ( μ 2 + z γ + 1 - c 2 ) 2 ,
(8)

with μ>0,c if γ { 0 } , c = 0 if γ ( - 1 , + ) \ ( { 0 } ) (see [5, 11, 31, 32]). Using these classification solutions scaled up and projected to satisfy the mixed boundary conditions up to a right order, the initial approximate solutions can be built up. Then through the finite dimensional reduction procedure and the notions of stability of critical points of the asymptotic reductional functional φ n,m , multiple blowing-up and concentrating solutions can be constructed as a small additive perturbation of the initial approximations.

This article is organized as follows. In Section 2, we will construct the approximate solution and rewrite the Equation (3) in terms of a linearized operator L. In Section 3, we give the invertibility of the linearized operator L, carry out the finite dimensional reduction and get the asymptotical expansion of the functional of the Equation (3) with respect to the suitable approximate solution. In Section 4, we give the proofs of Theorems 1.1 and 1.2.

2 Construction of the approximate solution

In this section, we will construct the approximate solution for the Equation (3). Let μ i , i = 1, ..., N + m, be positive numbers and set

α i = 0 , i = N + 1 , , N + m ,

and

Q i ( x ) = S ( x ) x - p i 2 α i .

Obviously, Q i (x) = S(x) for any i = N + 1, ..., N + m. Then the function

u i ( x ) = log 8 μ i 2 ( 1 + α i ) 2 ( μ i 2 ε 2 + x - p i 2 ( 1 + α i ) ) 2 f ( p i ) Q i ( p i ) ,
(9)

satisfies

- Δ u i = ε 2 x - p i 2 α i f ( p i ) Q i ( p i ) e u i , in 2 .
(10)

Set {k1, ..., k n } {1, ..., N} and k n+i = N + i for any i = 1, ..., m.

We hope to take i = 1 n + m u k i as an initial approximate solution of the problem (3). So we modify it to be

U ( x ) : = i = 1 n + m U k i = i = 1 n + m ( u k i + H k i t ) ,
(11)

where H k i t ( x ) (x) is the solution of

Δ H k i t = 0 in Ω , ε ( 1 - t ) ν H k i t + t b ( x ) H k i t = - ε ( 1 - t ) u k i ν + t b ( x ) u k i , on Ω .
(12)

Then U k i := u k i + H k i t satisfies

- Δ U k i = ε 2 x - p k i 2 α k i f ( p k i ) Q k i ( p k i ) e u k i , in Ω , ε ( 1 - t ) ν U k i + t b ( x ) U k i = 0 , on Ω .
(13)

Lemma 2.1 For t (0, 1] and p k i Ω ,

H k i t ( x ) = 4 ( 1 + α k i ) H t , ε ( x , p k i ) - log 8 μ k i 2 ( 1 + α k i ) 2 f ( p k i ) Q k i ( p k i ) + O ( ε 2 ) ,
(14)

uniformly in C ( Ω ̄ ) and in C loc 2 ( Ω ) for ε small.

Proof. Set z t ( x ) = H k i t ( x ) - 4 ( 1 + α k i ) H t , ε ( x , p k i ) + log 8 μ k i 2 ( 1 + α k i ) 2 f ( p k i ) Q k i ( p k i ) . Then z t (x) satisfies

Δ z t ( x ) = 0 , in Ω , ε ( 1 - t ) z t ( x ) ν + t b ( x ) z t ( x ) = F t ( x ) , on Ω ,

where

F t ( x ) = 4 ε ( 1 - t ) ( 1 + α k i ) x - p k i 2 α k i ν ( x ) ( x - p k i ) μ k i 2 ε 2 + x - p k i 2 ( 1 + α k i ) - ν ( x ) ( x - p k i ) x - p k i 2 + t b ( x ) 2 log μ k i 2 ε 2 + x - p k i 2 ( 1 + α k i ) - 4 ( 1 + α k i ) log x - p k i .

For any t (0, 1], it is easy to check F t x L Ω = O ε 2 . If t = 1, from the maximum principle and smooth function b(x) > 0, it follows

max Ω ̄ z 1 ( x ) = max Ω z 1 ( x ) C ( b ) F 1 ( x ) L ( Ω ) = O ( ε 2 ) .

If 0 < t < 1, from the maximum principle with the Robin boundary condition (see [[33], Lemma 2.6]), it also follows

max Ω ̄ z t ( x ) 1 t C ( b ) F t ( x ) L ( Ω ) = O ( ε 2 ) .

Thus using the interior estimate of derivative of harmonic function (see [[34], Theorem 2.10]), there holds

max K D α z t ( x ) 2 α dist ( K , Ω ) α max Ω ̄ z t ( x ) = O ( ε 2 ) ,

for any compact subset K of Ω, any t (0, 1] and any multi-index α with |α| ≤ 2, which derives (14) uniformly in C ( Ω ̄ ) and in C loc 2 ( Ω ) for ε small. □

From this lemma we can construct the approximate solution U ( x ) = i = 1 n + m ( u k i + H k i t ) , which satisfies the mixed boundary conditions. On the other hand, we hope that the error U ( x ) - u k i is smaller near every p k i . In fact, we can realize this point by further choosing positive number μ k i such that

log 8 μ k i 2 ( 1 + α k i ) 2 f ( p k i ) Q k i ( p k i ) = 4 ( 1 + α k i ) H t , ε ( p k i , p k i ) + 4 j = 1 , j i j = n + m ( 1 + α k j ) G t , ε ( p k i , p k j ) .
(15)

Consider the scaling of the solution of the Equation (3)

v ( y ) = u ( ε y ) + 4 log ε ,

then v(y) satisfies

- Δ v = S ( ε y ) f ( ε y ) e v , in Ω ε , ( 1 - t ) v ν + t b ( ε y ) v = 4 t b ( ε y ) log ε , on Ω ε ,
(16)

where Ω ε = 1 ε Ω.. We also set p k i = 1 ε p k i and define the new approximation in expanded variables as V(y) = U(εy) + 4 log ε. Furthermore, set

ρ k i = ε 1 1 + α k i ,
(17)

and

W ( y ) = S ( ε y ) f ( ε y ) e V ( y ) .
(18)

Obviously, ρ k n + i =ε for all i = 1, ..., m.

Here, we want to see how well -∆V(y) match with W(y) through V(y). A simple computation shows

- Δ V ( y ) = - ε 2 Δ x U ( x ) = - ε 2 i = 1 n + m Δ x u k i ( x ) + H k i t ( x ) = i = 1 n + m ε ρ k i 2 8 μ k i 2 ( 1 + α k i ) 2 ε y - p k i ρ k i 2 α k i μ k i 2 + ε y - p k i ρ k i 2 ( 1 + α k i ) 2 .

Then given a small number δ > 0, if ε y - p k i ρ k i > δ ρ k i for all i = 1, ..., n + m,

- Δ V ( y ) = O ( ε 4 ) ,
(19)

and if ε y - p k i ρ k i δ ρ k i for some i,

- Δ V ( y ) = ε ρ k i 2 8 μ k i 2 ( 1 + α k i ) 2 ε y - p k i ρ k i 2 α k i μ k i 2 + ε y - p k i ρ k i 2 ( 1 + α k i ) 2 + O ( ε 4 ) .
(20)

On the other hand, if ε y - p k i ρ k i > δ ρ k i for all i = 1, ..., n + m, obviously,

W ( y ) = O ( ε 4 ) ,
(21)

and if ε y - p k i ρ k i δ ρ k i for some i,

W ( y ) = ε y - p k i 2 α k i Q k i ( ε y ) f ( ε y ) e V ( y ) = ε 4 8 μ k i 2 ( 1 + α k i ) 2 ε y - p k i 2 α k i ε 2 μ k i 2 + ε y - p k i 2 ( 1 + α k i ) 2 f ( ε y ) Q k i ( ε y ) f ( p k i ) Q k i ( p k i ) × exp H k i t ( ε y ) + j = 1 , j i n + m log 8 μ k j 2 ( 1 + α k j ) 2 ε 2 μ k j 2 + ε y - p k j 2 ( 1 + α k j ) 2 f ( p k j ) Q k j ( p k j ) + H k j t ( ε y ) .

Now from (14), (15) and (17), we have

W ( y ) = ε ρ k i 2 8 μ k i 2 ( 1 + α k i ) 2 ε y - p k i ρ k i 2 α k i μ k i 2 + ε y - p k i ρ k i 2 ( 1 + α k i ) 2 1 + O ρ k i ε y - p k i ρ k i + O ( ε 2 ) .
(22)

In summary, we set

R ( y ) = Δ V ( y ) + W ( y ) ,
(23)

and if ε y - p k i ρ k i > δ ρ k i for all i = 1, ..., n + m,

R ( y ) = O ( ε 4 ) ,
(24)

while ε y - p k i ρ k i δ ρ k i for some i,

R ( y ) = ε ρ k i 2 8 μ k i 2 ( 1 + α k i ) 2 ε y - p k i ρ k i 2 α k i μ k i 2 + ε y - p k i ρ k i 2 ( 1 + α k i ) 2 O ρ k i ε y - p k i ρ k i + O ( ε 2 ) .
(25)

In the rest of this article, we try to find a solution v of the form v = V + ϕ of the Equation (16). In terms of ϕ, the problem (3) becomes

L ϕ = Δ ϕ + W ϕ = - [ R + N ( ϕ ) ] , in Ω ε , ( 1 - t ) ϕ ν + t b ( ε y ) ϕ = 0 , on Ω ε ,
(26)

where

N ( ϕ ) = W [ e ϕ - 1 - ϕ ] .
(27)

3 The finite dimensional reduction

In this section, we will carry out the finite dimensional reduction to solve the Equation (26). First of all, we need to get the desired invertibility of linearized operation L. Set

z i 0 ( z ) = z 2 ( 1 + α k i ) - μ k i 2 z 2 ( 1 + α k i ) + μ k i 2 , for i = 1 , 2 , , n + m , z i j ( z ) = 4 z j z 2 + μ k i 2 , for i = n + 1 , , n + m , j = 1 , 2 , L i ϕ = Δ ϕ + 8 μ k i 2 ( 1 + α k i ) 2 z 2 α k i μ k i 2 + z 2 ( 1 + α k i ) 2 ϕ , for i = 1 , , n + m .

A basic fact to get the needed invertibility is that the linearized operator L formally approaches to the operator L i under suitable dilations and translations, which have some well-known properties that any bounded solution of L t ϕ = 0 is

  • a linear combination of zi 0and z ij for i = n + 1, ..., n + m, j = 1, 2 (see [24, 35]);

  • proportional to z i0 for 0< α k i and i = 1, 2, ..., n (see [20, 28, 36]).

Remark 3.1 These properties of the operator L i have been discussed in the above papers only if 0< α k i for i = 1, ..., n, or α k i =0 for i = n + 1, ..., n + m. In fact, if -1< α k i <0 for some i = 1, ..., n, the operator L i has also the corresponding properties.

Lemma 3.2 For -1<α { 0 } , any bounded solution ϕ of

Δ ϕ + 8 μ 2 ( 1 + α ) 2 z 2 α μ 2 + z 2 ( 1 + α ) 2 ϕ = 0 , z \ { 0 } ,
(28)

is proportional to z 2 ( 1 + α ) - μ 2 z 2 ( 1 + α ) + μ 2 .

Proof. If we express the bounded solution ϕ of the Equation (28) in Fourier expansion form as follow

ϕ ( z ) = n = - + u n ( r ) e - i n θ , z = r e i θ ,

u n (r) is a bounded nontrivial solution of the equation

u ( r ) + 1 r u ( r ) - n 2 r 2 u + 8 μ 2 ( 1 + α ) 2 r 2 α ( μ 2 + r 2 ( 1 + α ) ) 2 u = 0 .
(29)

Since any solution of -∆u = eu in is given by the Liouville formula

ln 8 F ( z ) 2 1 + F ( z ) 2 2 ,

for any meromorphic function F defined on { z : F ( z ) 0 } , the function

ln 8 μ 2 ( 1 + α ) 2 1 + n + α + 1 α + 1 a z n 2 ( μ 2 + z 2 ( 1 + α ) 1 + a z n 2 ) 2 ,

with any n and a < α + 1 n + α + 1 , is the solution of -∆u = |z|2αeu in \ { 0 } . Moreover, its derivative with respect to a at a = 0

ϕ n ( z ) = 1 α + 1 ( n + α + 1 ) μ 2 + ( n - α - 1 ) z 2 ( 1 + α ) μ 2 + z 2 ( 1 + α ) z n , n ,

is a solution of the Equation (29) with r = |z|.

For |n| ≥ 1, since {ϕ n (r), ϕ -n (r)} is a set of linearly independent solutions of the second order linear homogeneous ODE (29), any bounded solution is a linear combination of ϕ n (r) and ϕ -n (r). However, ϕ|n|(r) ( resp. ϕ -n (r) ) tends to 0 ( resp. ∞ ) as r 0 and ϕ|n|(r) ( resp. ϕ -|n| (r) ) tends to ∞ ( resp. 0 ) as r + ∞, which implies that the Equation (29) ||n|≥1 has no bounded nontrivial solution.

For n = 0, ϕ 0 ( z ) =- z 2 ( 1 + α ) - μ 2 z 2 ( 1 + α ) + μ 2 is a bounded solution of the Equation (29)|n = 0, that is, of the Equation (28). We claim that there does not exist the second linearly independent bounded solution of the Equation (29)|n = 0. Otherwise, let ω be another linearly independent bounded solution of (29)|n = 0. Writing ω(r) = c(r)ϕ0(r), we get that

c ( r ) ϕ 0 + c ( r ) 2 ϕ 0 + 1 r ϕ 0 = 0 .

Then there exists a constant C > 0 such that

c ( r ) = C r ϕ 0 2 ( r ) = C ( r 2 ( 1 + α ) + μ 2 ) 2 r ( r 2 ( 1 + α ) - μ 2 ) 2 ~ C r for r small, c ( r ) ~ C log r for r small .

Hence, ω(r) ~ C log r for r small, which implies ω(r) is unbounded on (0, + ∞). It contradicts the assumption that ω is bounded. □

Let us denote

Z i 0 ( y ) = z i 0 ε y - p k i ρ k i , for i = 1 , 2 , , n + m , Z i j ( y ) = z i j ε y - p k i ρ k i , for i = n + 1 , , n + m , j = 1 , 2 , χ ι ( y ) = χ ε y - p k i ρ k i , for i = 1 , 2 , , n + m ,

where χ(r) is a smooth, non-increasing cut-off function such that for a large but fixed number R0 > 0, χ(r) = 1 if rR0, and χ(r) = 0 if rR0 + 1. Additionally, set α0 = min{0, α1, ..., α N }. For any α (-1, α0), we introduce the Banach space

C n , m : = { ψ L ( Ω ε ) : ψ n , m < + } ,

with the norm

ψ n , m = sup y Ω ε ψ ( y ) ε 2 + i = 1 , α i < 0 n ε ρ i 2 ε y - p i ρ i 2 α i 1 + ε y - p i ρ i 4 + 2 α + 2 α i + i = 1 , α i 0 n + m ε ρ i 2 1 + ε y - p i ρ i 4 + 2 α .
(30)

Now to get the invertibility of the linearized operator L, we only need to solve the following linear problems: given h of class C n , m C 0 , β Ω ε with β (0,1), for m ≥ 1 and 0 ≤ nN, we find a function ϕ and scalars c ij , i = n + 1, ..., n + m, j = 1, 2, such that

L ϕ = Δ ϕ + W ϕ = h + j = 1 2 i = n + 1 n + m c i j χ i Z i j , in Ω ε , ( 1 - t ) ν ϕ + t b ( ε y ) ϕ = 0 , on Ω ε , Ω ε χ i Z i j ϕ d y = 0 , i = n + 1 , , n + m , j = 1 , 2 ,
(31)

and for m = 0 and 1 ≤ n ≤ N, we find a function ϕ such that

L ϕ = Δ ϕ + W ϕ = h , in Ω ε , ( 1 - t ) ν ϕ + t b ( ε y ) ϕ = 0 , on Ω ε .
(32)

Proposition 3.1 (i) If m ≥ 1 and 0 ≤ nN, given a fixed number δ > 0, there exist positive numbers ε0 and C such that for any points p k l , l = n+ 1, ..., n + m, in Ω', with

dist( p k l , Ω ) δ , p k l - p k i δ , f o r l i a n d i = 1 , , n + m ,
(33)

there is a unique solution ϕ L ( Ω ε ) , c n + 1 ,, c n + m , of the Equation (31), which satisfies

ϕ C log 1 ε h n , m ,
(34)

for all ε < ε0 and t (0, 1]. Moreover, the map p' ϕ is C1 and

D p ϕ C log 1 ε 2 h n , m ,
(35)

where p : = 1 ε p k n + 1 , , 1 ε p k n + 1 .

(ii) If m = 0 and 1 ≤ n ≤ N, there exist positive numbers ε0 and C such that there is a unique solution ϕ L ε ) of the Equation (32), which satisfies

ϕ C log 1 ε h n , 0 ,
(36)

for all ε < ε0 and t (0, 1].

These results can be established through some technical lemmas. First for the linear Equation (32) under the additional orthogonality conditions with respect to Zi 0, i = 1, ..., n + m, and Z ij , i = n + 1, ..., n + m, j = 1, 2, we prove the following priori estimates.

Lemma 3.3 (i) If m ≥ 1 and 0 ≤ n ≤ N, given a fixed number δ > 0, there exist positive numbers ε0 and C such that for any points p k l , l = n+ 1, ..., n + m, in Ω', which satisfy the relation (33), and any solution ϕ of the Equation (32) with t (0, 1] under the orthogonality conditions

Ω ε χ i Z i 0 ϕ d y = 0 , i = 1 , , n + m , Ω ε χ i Z i j ϕ d y = 0 , i = n + 1 , , n + m , j = 1 , 2 ,
(37)

one has

ϕ C h n , m ,
(38)

for all ε < ε0.

(ii) If m = 0 and 1 ≤ n ≤ N, there exist positive numbers ε 0 and C such that for any solution ϕ of the Equation (32) with t (0, 1] under the orthogonality conditions

Ω ε χ i Z i 0 ϕ d y = 0 , i = 1 , , n ,
(39)

one has

ϕ C h n , 0 ,
(40)

for all ε < ε0.

Remark 3.4 The idea behind these estimates partly comes from observing the linear Equation (32) with h = 0 on bounded set B i , R := y Ω ε : ε y - p k i ρ k i < R for ε small. After a translation and a rotation so that Ω ε converges to the whole plan 2 , the Equation (32) approaches L i ϕ = 0 in 2 . As a result, the solution of the Equation (32) under the additional orthogonality conditions (37) should be zero.

Proof. Case (i): First consider the "inner norm" ϕ l = sup i = 1 n + m B i , R ̄ ϕ and the "boundary norm" ϕ o = sup Ω ε ϕ , we claim that there is a constant C > 0 such that if Lϕ = h in Ω ε , then

ϕ C ϕ l + ϕ o + h n , m .
(41)

We will establish it with the help of suitable barrier.

Consider that the function g ( z ) = z 2 ( 1 + α ) - 1 z 2 ( 1 + α ) + 1 is a radial solution in 2 of

Δ g ( z ) + 8 ( 1 + α ) 2 z 2 α 1 + z 2 ( 1 + α ) 2 g ( z ) = 0 ,

we define a bounded comparison function

Z ( y ) = i = 1 n + m g a ε y - p k i ρ k i , y Ω ε ,

with a > 0. Set R a = 1 a 3 1 2 ( 1 + α ) . While ε y - p k i ρ k i R a for all i = 1, ..., n + m,

- Δ Z ( y ) = i = 1 n + m ε ρ k i 2 a 2 8 ( 1 + α ) 2 a ε y - p k i ρ k i 2 α 1 + a ε y - p k i ρ k i 2 ( 1 + α ) 2 g a ε y - p k i ρ k i i = 1 n + m ε ρ k i 2 a 2 4 ( 1 + α ) 2 a ε y - p k i ρ k i 2 α 1 + a ε y - p k i ρ k i 2 ( 1 + α ) 2 > i = 1 n + m ε ρ k i 2 a - 2 ( 1 + α ) ( 1 + α ) 2 ε y - p k i ρ k i 2 α + 4 .
(42)

Moreover, according to (21) and (22), on the same region,

W ( y ) Z ( y ) C i = 1 n + m ε ρ k i 2 1 ε y - p k i ρ k i 2 α k i + 4 .
(43)

So if a is small enough to satisfy (1 + α)2a-2(1+α) > C + 1, R a is sufficiently large. As a result, by (42) and (43), for any R ≥ R a , we have Z(y) > 0 and L(Z) < 0 in Ω R , ε c : = Ω ε \ i = 1 n + m B i , R .

Let M be a large number such that for all i = 1, ..., n + m, ΩB ( p k i , M ) . Consider now the solution of the problem

- Δ ψ i = ε ρ k i 2 4 ε y - p k i ρ k i 2 α + 4 + 4 ε 2 , R < ε y - p k i ρ k i < M ρ k i , ψ i ( y ) = 0 , for ε y - p k i ρ k i = R and ε y - p k i ρ k i = M ρ k i .

A direct computation shows

ψ i ( y ) = φ i ( r i ) - φ i M ρ k i log r i R log M R ρ k i ,

where

r i = ε y - p k i ρ k i ,

and

φ i ( t ) = 1 ( 1 + α ) 2 1 R 2 ( 1 + α ) - 1 t 2 ( 1 + α ) + ρ k i 2 ( R 2 - t 2 ) .

For the sake of the convenience, we choose R larger if necessary. Then it easily see that these functions ψ i , i = 1, ..., n + m, have a uniform bound independent of ε.

Now we can construct the needed barrier:

ϕ ̃ ( y ) = 2 ( ϕ l + ϕ o ) Z ( y ) + h n , m i = 1 n + m ψ i ( y ) .

It is easy to check that L ϕ ̃ <h=Lϕ in Ω R , ε c , and ϕ ̃ ϕ on Ω R , ε c . Since Z(y) > 0 and LZ(y) < 0 in Ω R , ε c , from the maximum principle (see [[37], Theorem 10, Chap. 2 ]), it follows that ϕ ̃ ϕ in Ω R , ε c . Similarly, - ϕ ̃ ϕ in Ω R , ε c , which derives the estimate (41).

We prove the priori estimate (38) by contradiction. Assume that there exist a sequence ε k → 0, points p k l k , l = n + 1, ..., n + m, in Ω' which satisfy relation (33), functions h k with ║h k n,m → 0, solutions ϕ k with ║ϕ k = 1, such that

L ϕ k = Δ ϕ k + W ϕ k = h k , in Ω ε , ( 1 - t ) ν ϕ k + t b ( ε y ) ϕ k = 0 , on Ω ε , Ω ε χ i Z i 0 ϕ k d y = 0 , i = 1 , , n + m , Ω ε χ i Z i 0 ϕ k d y = 0 , i = n + 1 , , n + m , j = 1 , 2 .

Then from the estimate (41), ║ϕ k l κ or ║ϕ k o κ for some κ > 0. Briefly set ε:= ε k , p k i : p k i k . If ║ϕ k l κ, with no loss of generality, we assume that sup B i , R ϕ k k for some i Then if we set ϕ ^ k ( z ) = ϕ k ρ k i ε z + p k i ε and ĥ k ( z ) = h k ρ k i ε z + p k i ε , ϕ ^ k satisfies

Δ ϕ ^ k + 8 μ k i 2 ( 1 + α k i ) 2 z 2 α k i [ μ k i 2 + z 2 ( 1 + α k i ) ] 2 [ 1 + O ( ρ k i z ) + O ( ε 2 ) ] ϕ ^ k = ρ k i ε 2 ĥ k ,

for z B R (0). Obviously, for any q 1 , - 1 α we easily get ρ k i ε 2 ĥ k 0 in Lq(B R (0)). Since 8 μ k i 2 ( 1 + α k i ) 2 z 2 α k i [ μ k i 2 + z 2 ( 1 + α k i ) ] 2 is bounded in Lq(B R (0)) and ϕ ^ k =1, elliptic regularity theory readily implies that ϕ ^ k converges uniformly over compact subsets near the origin to a bounded nontrivial solution ϕ ^ of the equation

L i ϕ ^ = Δ ϕ ^ + 8 μ k i 2 ( 1 + α k i ) 2 z 2 α k i [ μ k i 2 + z 2 ( 1 + α k i ) ] 2 ϕ ^ = 0 , in 2 .

From Lemma 3.2, this equation implies that ϕ ^ is proportional to zi 0for i = 1, ..., n, or a linear combination of zi 0and z ij for i = n + 1, ..., n + m, j = 1, 2. However, our assumed orthogonality conditions (37) on ϕ k pass to limit and yield the corresponding conditions (37) on ϕ ^ , which means ϕ ^ 0. Hence, it is absurd because ϕ ^ is nontrivial.

If ║ϕ k o ≥ κ and ║ϕ k l → 0, there exists a point q Ω and a number R1 > 0 such that sup Ω ε B R 1 ( q ) ϕ k ( y ) k>0 with q = 1 ε q. Consider ϕ ^ k ( y ) = ϕ k ( y - q ) and let us translate and rotate Ω ε so that q' = 0 and Ω ε approaches the upper half-plan + 2 . Since ε q - p k i ρ k i > δ ρ k i for all i = 1, ..., n + m, ϕ ^ k ( z ) satisfies

Δ ϕ ^ k ( y ) + O ( ε 4 ) ϕ ^ k ( y ) = h k ( y ) , in Ω ε \ i = 1 n + m B i , δ ρ k i - 1 , ( 1 - t ) ν ϕ ^ k + t b ( ε y ) ϕ ^ k = 0 , on Ω ε ,

with ( 1 - t ) Ω ε ϕ ^ k 2 +t Ω ε b ( ε y ) ϕ ^ k 2 <C. Moreover, we easily get h k (y) → 0 in Ω ε \ i = 1 n + m B i , δ ρ k i - 1 . While t = 1, it is obvious to see that ϕ ^ k ( y ) =0 on Ω ε . So it is absurd because of sup Ω ε B R 1 ( q ) ϕ k ( y ) k>0. On the other hand, for any t (0,1), elliptic regularity theory with the Robin boundary condition (see [30, 34, 38] and the references therein) implies that ϕ ^ k converges uniformly on compact subsets near the origin to a bounded nontrivial solution ϕ ^ of the equation

Δ ϕ ^ = 0 , in + 2 , ( 1 - t ) ν ϕ ^ + t b ( 0 ) ϕ ^ = 0 , on + 2 ,

with ( 1 - t ) + 2 ϕ ^ 2 +tb ( 0 ) + 2 ϕ ^ 2 <C. It follows that its bounded solution ϕ ^ is zero. Hence, it is also absurd because ϕ ^ is nontrivial, which derives the priori estimate (38) of the case (i). Since the proof of the case (ii) is similar to that of the case (i), we omit it.□

We will give next the priori estimate for the solution of the Equation (32) that satisfies orthogonality conditions with respect to Z ij , i = n + 1, ..., n + m, j = 1, 2, only.

Lemma 3.5. (i) If m ≥ 1 and 0 ≤ n ≤ N, given a fixed number δ > 0, there exist positive numbers ε0 and C such that for any points p k l , l = n + 1, ..., n + m, in Ω', which satisfy the relation (33), and any solution ϕ of the Equation (32) with t (0, 1] under the orthogonality conditions

Ω ε χ i Z i j ϕ d y = 0 , i = n + 1 , , n + m , j = 1 , 2 ,
(44)

one has

ϕ C log 1 ε h n , m ,
(45)

for all ε < ε0.

(ii) If m = 0 and 1 ≤ n ≤ N, there exist positive numbers ε0 and C such that for any solution ϕ of the Equation (32) with t (0, 1], one has

ϕ C log 1 ε h n , 0 ,
(46)

for all ε < ε0.

Proof. Case (i): Let ϕ satisfy the Equation (32) under the orthogonality conditions (44). We will modify ϕ to satisfy the orthogonality conditions (37). To realize this point, we consider some related modifications with compact support of the functions Zi 0, i = 1, ..., n + m.

Let R > R0 + 1 be large and fixed, and let i 0 be the solution of the equation

Δ i 0 + 8 μ k i 2 ( 1 + α k i ) 2 z 2 α k i ( μ k i 2 + z 2 ( 1 + α k i ) ) 2 i 0 = 0 , for R < z < δ 3 ρ k i , i 0 ( z ) = z i 0 ( R ) on z = R , i 0 ( z ) = 0 on z = δ 3 ρ k i .

A simple computation shows that this solution is explicitly given by

i 0 ( z ) = z i 0 ( r ) 1 - R r d s s z i 0 2 ( s ) R δ 3 ρ k i d s s z i 0 2 ( s ) , r = z .

Set

i 0 ( y ) = i 0 ε y - p k i ρ k i , η 1 i ( y ) = η 1 ε y - p k i ρ k i , η 2 i ( y ) = η 2 4 ρ k i ε y - p k i ρ k i ,

where η1(r) and η2(r) are smooth cut-off functions with the properties: η1(r) = 1 for r < R, η1(r) = 0 for r>R+1, η 1 ( r ) 2; η2(r) = 1 for r < δ, η2(r) = 0 for r> 4 δ 3 , η 2 ( r ) C η 2 ( r ) C. We define a test function

Z ̃ i 0 = η 1 i Z i 0 + ( 1 - η 1 i ) η 2 i i 0 .

Obviously, Z ̃ i 0 ( y ) = Z i 0 ( y ) if ε y - p k i ρ k i <R, and Z ̃ i 0 ( y ) =0 if ε y - p k i ρ k i > δ 3 ρ k i , in particular, Z ̃ i 0 ( y ) =0 near ∂Ωɛ.

Now we modify ϕ to satisfy the orthogonality conditions with respect to Zi 0, i = 1, ..., n+ m, and set

ϕ ̃ = ϕ + i = 1 m + n d i Z ̃ i 0 ,

where the numbers d i are chosen such that

d i Ω χ i Z i 0 2 + Ω χ i Z i 0 ϕ = 0 , i = 1 , , m + n .

Thus

L ϕ ̃ = h + i = 1 m + n d i L Z ̃ i 0 , in Ω ε , ( 1 - t ) ν ϕ ̃ + t b ( ε y ) ϕ ̃ = 0 , on Ω ε ,
(47)

and ϕ ̃ satisfies all the orthogonality conditions in (37). From Lemma 3.3 (i), we have

ϕ ̃ C h n , m + i = 1 n + m d i L Z ̃ i 0 n , m .
(48)

In order to get the estimate (45) of ϕ, we need to give the sizes of d i and L Z ̃ i 0 n , m for any t (0, 1]. Multiplying the first equation of (47) by Z ̃ j 0 , integrating by parts and using the mixed boundary conditions of (47), we get

L Z ̃ j 0 , ϕ ̃ = Z ̃ j 0 , h + d j L Z ̃ j 0 , Z ̃ j 0 ,
(49)

where f , g = Ω ε fg. A simple computation shows that Z ̃ j 0 , h C h n , m , which in combination with (48) and (49) yields

d j L Z ̃ j 0 , Z ̃ j 0 C h n , m 1 + L Z ̃ j 0 n , m + C L Z ̃ j 0 n , m i = 1 m + n d i L Z ̃ i 0 n , m .
(50)

From some similar computations (see [[26], Lemma 3.2] and [[28], Lemma 3.3]), there exists a constant C > 0 independent of ε such that

L Z ̃ j 0 n , m C 1 log ρ k j ,
(51)

and

L Z ̃ j 0 , Z ̃ j 0 - C log ρ k j 1 + O 1 log ρ k j ,
(52)

which combined with (50) yields

d j C log 1 ρ k j h n , m .
(53)

Furthermore, from (48), (51), (53), and the definitions of ϕ ̃ and ρ k j , we have

ϕ C log 1 ε h n , m .
(54)

Similarly, the proof of the case (ii) can also be done, we omit it. □

Proof of Proposition 3.1. Case (i): From Lemma 3.5 (i), and the Fredholm's alternative theory with Robin boundary condition instead of Dirichlet boundary condition if necessary (see [34, 38] and the references therein), the proof can be similarly given through those in [[26], pp. 61-63].

Case (ii): Since the priori estimate (36) of the solution of the Equation (32) has been established in Lemma 3.5 (ii), we can use the Fredholm's alternative and obtain the unique solution of the Equation (32). □

Let us now introduce the auxiliary nonlinear problems: for m ≥ 1 and 0 ≤ n < N, we find the function ϕ and scalars c ij , i = n + 1, ..., n + m, j = 1, 2, such that

L ϕ = - R + N ( ϕ ) + j = 1 2 i = n + 1 n + m c i j χ i Z i j , in Ω ε , ( 1 - t ) ν ϕ + t b ( ε y ) ϕ = 0 , on Ω ε , Ω ε χ i Z i j ϕ d y = 0 , i = n + 1 , , n + m , j = 1 , 2 ,
(55)

and for m = 0 and 1 ≤ n ≤ N, we find the solution ϕ of the nonlinear Equation (26).

The following result can be proved using standard arguments as in [26, 28].

Proposition 3.2 (i) If m ≥ 1 and 0 n ≤ N, given a fixed number δ > 0, there exist positive numbers ε0 and C such that for any points p k l , l= n + 1, ..., n + m, in Ω' satisfying the relation (33), there is a unique solution for the Equation (55) which satisfies

ϕ C ρ log ε ,
(56)

for all ε < ε0 and t (0, 1]. Moreover, the map p'ϕ is C1 and

D p ϕ C ρ log ε 2 ,
(57)

where ρ:= max 1 i n + m ρ k i and p := 1 ε p k n + 1 , , 1 ε p k n + m .

(ii) If m = 0 and 1≤n≤N, there exist positive numbers ε0 and C such that there is a unique solution for the Equation (26) which also satisfies the estimate (56) for all ε < ε0 and t (0, 1].

Now we only need to find a solution to the Equation (26) with m ≥ 1 and 0 ≤ n ≤ N, and hence to the Equation (55) if p = 1 ε p k n + 1 , , 1 ε p k n + m is such that

c i j ( p ) = 0 , i = n + 1 , , n + m , j = 1 , 2 .
(58)

Let us introduce the energy functional of the Equation (3), namely for t = 1,

J ε , 1 ( u ) = 1 2 Ω u 2 - ε 2 Ω S ( x ) f ( x ) e u ,

and for t (0,1),

J ε , t ( u ) = ( 1 - t ) 1 2 Ω u 2 - ε 2 Ω S ( x ) f ( x ) e u + t 2 ε Ω b ( x ) u 2 .

Furthermore, we define

F ε , t ( p ) = J ε , t ( U ( p ) + ϕ ̃ ( p ) ) ,
(59)

where p= ( p k n + 1 , , p k n + m ) ( Ω ) m and ϕ ̃ ( p ) ( x ) =ϕ x ε , p x with the solution ϕ of the Equation (55).

The finite dimensional variational reduction is meaningful in view of the following property.

Proposition 3.3 If p= ( p k n + 1 , , p k n + m ) ( Ω ) m satisfying the relation(33) is a critical point of F ε,t with t (0, 1], then U ( p ) + ϕ ̃ ( p ) is a critical point of J ε,t , namely a solution of the Equation (3). Besides, on any compact subsets S of (Ω' \ Γ)m \ Δ m the following expansion holds

F ε , t ( p ) = J ε , t ( U ( p ) ) + θ ε , t ( p ) ,
(60)

where

θ ε , t ( p ) 0 , u n i f o r m l y o n S ,
(61)

for ε small.

Proof. Step 1: Let us define for t = 1,

I ε , 1 ( v ) = 1 2 Ω ε v 2 - Ω ε S ( ε y ) f ( ε y ) e v ,

and for t (0,1),

I ε , t ( v ) = ( 1 - t ) 1 2 Ω ε v 2 - Ω ε S ( ε y ) f ( ε y ) e v + t 2 Ω ε b ( ε y ) ( v - 4 ε log ε ) 2 .

Then F ε , t ( p ) = J ε , t ( U ( p ) + ϕ ̃ ( p ) ) = I ε , t ( V ( p ) + ϕ ( p ) ) , where p = 1 ε p. Moreover, for k = n + 1, ..., n + m, l = 1, 2, it holds

p k l F ε , t ( p ) = ε - 1 D I ε , t ( V + ϕ ) p k l V + p k l ϕ .
(62)

Since ϕ(p') is a solution of the Equation (55), v = V(p') + ϕ(p') satisfies

v + S ( ε y ) f ( ε y ) e v = j = 1 2 i = n + 1 n + m c i j χ i Z i j , in Ω ε , ( 1 - t ) v ν + t b ( x y ) v = 4 t b ( ε y ) log ε , on Ω ε .
(63)

By (62), (63), D p F ε,t (p) = 0 implies for t (0, 1],

j = 1 2 i = n + 1 n + m c i j Ω ε χ i Z i j ( p k l V + p k l ϕ ) = 0 .

From the definition of V, it can be directly checked p k l V=- Z k l +o ( 1 ) , where o(1) is in the sense of the L-norm for ε small. Since D p ϕ Cρ log ε 2 , p k l V+ p k l ϕ=- Z k l +o ( 1 ) . Hence it follows

j = 1 2 i = n + 1 n + m c i j Ω ε χ i Z i j ( Z k l + o ( 1 ) ) = 0 , k = n + 1 , , n + m , l = 1 , 2 ,

which is a strictly diagonal dominant system. This implies that c ij = 0, i = n + 1, ..., n + m, j = 1, 2. By (63), U ( p ) + ϕ ̃ ( p ) is a critical point of J ε,t , that is, a solution of the Equation (3).

Step 2: Set θ ̃ ε , t ( p ) = I ε , t ( V ( p ) + ϕ ( p ) ) - I ε , t ( V ( p ) ) . Using DI ε,t (V + ϕ)[ϕ] = 0, a Taylor expansion and an integration by parts, it follows that for t (0, 1),

θ ̃ ε , t ( p ) = I ε , t ( V + ϕ ) - I ε , t ( V ) = 0 1 D 2 I ε , t ( V + s ϕ ) [ ϕ ] 2 ( 1 - s ) d s = 0 1 ( 1 - t ) Ω ε ϕ 2 - S ( ε y ) f ( ε y ) e V + s ϕ ϕ 2 + t Ω ε b ( ε y ) ϕ 2 ( 1 - s ) d s = ( 1 - t ) 0 1 Ω ε N ( ϕ ) + R ϕ + W 1 - e s ϕ ϕ 2 ( 1 - s ) d s ,
(64)

and similarly, for t = 1,

θ ̃ ε , 1 ( p ) = 0 1 Ω ε N ( ϕ ) + R ϕ + W 1 - e s ϕ ϕ 2 ( 1 - s ) d s .
(65)

Note that ║ϕ = O(ρ|log ε|), ║N (ϕ) n,m = O(ρ2| log ε|2), ║R║ n,m = O(ρ), and ║W║ n,m = O(1). Then from (64), (65), it is easy to deduce for t (0, 1],

θ ̃ ε , t ( p ) = O ( ρ 2 log ε ) .
(66)

Hence, from (66), the expansion (60) satisfies the property (61).□

Finally, we need to write the precisely asymptotical expansion of J ε,t (U). To realize it, we first establish the following result:

Lemma 3.6 Assume that points p k n + 1 , , p k n + m , satisfy the relation (33), then fort (0,1) and i = 1, ..., n + m, there hold

G t , ε ( x , p k i ) = G 1 ( x , p k i ) + O ( ε ) ,
(67)

and

H t , ε ( x , p k i ) = H 1 ( x , p k i ) + O ( ε ) ,
(68)

uniformly in C ( Ω ̄ ) and in C loc 2 ( Ω ) for ε small.

Proof. Set z i t ( x ) = G t , ε ( x , p k i ) - G 1 ( x , p k i ) . Then z i t ( x ) satisfies

z i t ( x ) = 0 , in Ω , ε ( 1 - t ) z i t ( x ) ν + t b ( x ) z i t ( x ) = F i t ( x ) , on Ω ,

where F i t ( x ) =-ε ( 1 - t ) ν G 1 ( x , p k i ) . Since H 1 ( x , p k i ) has the gradient estimate on Ω (see [[26], p. 76])

x H 1 ( x , p k i ) C 1 min 1 x - p k i , 1 dist ( p k i , Ω ) + C 2 C 1 1 δ + C 2 ,

we can easily get F i t ( x ) L ( Ω ) =O ( ε ) . Using the same technique with the proof of Lemma 2.1, we can also get z i t ( x ) =O ( ε ) uniformly in C ( Ω ̄ ) and in C loc 2 ( Ω ) for ε small, which means (67). From the definition of the regular part of Green function, we can also derive (68). □

Proposition 3.4 The following asymptotical expansions hold for t = 1,

J ε , 1 ( U ) = 8 π m log 8 + 8 π i = 1 n ( 1 + α k i ) log 8 ( 1 + α k i ) 2 f ( p k i ) Q k i ( p k i ) - 16 π i , j = 1 , i j n ( 1 + α k i ) ( 1 + α k j ) G 1 ( p k j , p k i ) - 16 π log ε i = 1 n + m ( 1 + α k i ) - 16 π i = 1 n ( 1 + α k i ) 2 H 1 ( p k i , p k i ) - 16 π i = 1 n + m ( 1 + α k i ) + ρ Θ ε ( p ) - 16 π φ n , m ( p ) ,
(69)

and for t (0,1),

J ε , t ( U ) = 8 π ( 1 - t ) m log 8 - 16 π ( 1 - t ) i = 1 n + m ( 1 + α k i ) - 16 π ( 1 - t ) log ε i = 1 n + m ( 1 + α k i ) + 8 π ( 1 - t ) i = 1 n ( 1 + α k i ) log 8 ( 1 + α k i ) 2 f ( p k i ) Q i ( p k i ) - 16 π ( 1 - t ) i = 1 n ( 1 + α k i ) 2 H 1 ( p k i , p k i ) - 16 π ( 1 - t ) i , j = 1 , i j n ( 1 + α k i ) ( 1 + α k j ) G 1 ( p k i , p k j ) + ρ Θ ε ( p ) - 16 π ( 1 - t ) φ n , m ( p ) ,
(70)

where φ n,m (p) is defined by (5), and for ε small, Θ ε is a bounded, smooth function of p = ( p k n + 1 , , p k n + m ) , uniformly on points p k n + 1 , , p k n + m in Ω' satisfying the relation (33).

Proof. According to [[26, 28], Lemma 6.1], it only remains to discuss the asymptotical expansion of the energy J ε,t (U) with respect to t (0,1). By (11), it follows

Ω U 2 = i , j = 1 n + m Ω U i U j = i , j = 1 n + m Ω U i U j ν - Ω U i Δ U j ,

and

Ω b ( x ) U 2 = i , j = 1 n + m Ω b ( x ) U i U j ,

which together with the Equation (13) yields

1 - t 2 Ω U 2 + t 2 ε Ω b ( x ) U 2 = - 1 - t 2 i , j = 1 n + m Ω U i Δ U j .

Furthermore

J ε , t ( U ) = ( 1 - t ) 1 2 Ω U 2 - ε 2 Ω S ( x ) f ( x ) e U + t 2 ε Ω b ( x ) U 2 = - ( 1 - t ) 1 2 i , j = 1 n + m Ω U i Δ U j + ε 2 Ω S ( x ) f ( x ) e U .
(71)

From (9), (13), and (14), it implies

- Ω U i Δ U j = Ω ε 2 | x - p k j | 2 α k j f ( p k j ) Q k j ( p k j ) ( u k i + H k i t ) e u k j = Ω 8 μ k j 2 ε 2 ( 1 + α k j ) 2 | x - p k j | 2 α k j ( μ k j 2 ε 2 + | x - p k j | 2 ( 1 + α k j ) ) 2 . [ log 1 ( μ k i 2 ε 2 + | x - p k i | 2 ( 1 + α k i ) ) 2 + 4 ( 1 + α k i ) H t , ε ( x , p k i ) + O ( ε 2 ) ] = Ω ρ k j 16 μ k j 2 ( 1 + α k j ) 2 α k j ( μ k j 2 + | y | 2 ( 1 + α k j ) ) 2 log 1 μ k i 2 ε 2 + | ρ k j y + p k j - p k i | 2 ( 1 + α k i ) + Ω ρ k j 32 μ k j 2 ( 1 + α k i ) ( 1 + α k j ) 2 | y | 2 α k j ( μ k j 2 + | y | 2 ( 1 + α k j ) ) 2 H t , ε ( ρ k j y + p k j , p k i , ) + O ( ε 2 ) .

Then if i ≠ j for any i and j,

- Ω U i Δ U j = 32 π ( 1 + α k i ) ( 1 + α k j ) G t , ε ( p k i , p k j ) + O ( ρ ) ,
(72)

and if i = j for any i,

- Ω U i Δ U i = 32 π ( 1 + α k i ) 2 H t , ε ( p k i , p k i ) - 32 π ( 1 + α k i ) log ( μ k i ε ) - 16 π ( 1 + α k i ) + O ( ρ ) .
(73)

On the other hand, from (21) and (22), it follows

Ω ε W = i = 1 n + m B i , δ ρ k i - 1 W + Ω ε \ i = 1 n + m B i , δ ρ k i - 1 W = i = 1 n + m B i , δ ρ k i - 1 ε ρ k i 2 8 μ k i 2 ( 1 + α k i ) 2 ε y - p k i p k i 2 α k i [ μ k i 2 + ε y - p k i p k i 2 ( 1 + α k i ) ] 2 1 + O ρ k i ε y - p k i p k i + O ( ε 2 ) + O ( ε 2 ) = i = 1 n + m z δ ρ k i 8 μ k i 2 ( 1 + α k i ) 2 z 2 α k i μ k i 2 + z 2 ( 1 + α k i ) 2 1 + O ( ρ k i z ) + O ( ε 2 ) d z + O ( ε 2 ) = i = 1 n + m 8 π ( 1 + α k i ) + O ( ρ ) .

As a result, it derives

ε 2 Ω S ( x ) f ( x ) e U = Ω ε W ( y ) = i = 1 n + m 8 π ( 1 + α k i ) + O ( ρ ) .
(74)

Now using the choice for μ k i 's by (15), together with (71)-(74), it holds

J ε , t ( U ) = 8 π ( 1 - t ) i = 1 n + m ( 1 + α k i ) log 8 ( 1 + α k i ) 2 f ( p k i ) Q k i ( p k i ) - 16 π ( 1 - t ) i = 1 n + m ( 1 + α k i ) 2 H t , ε ( p k i , p k i ) - 16 π ( 1 - t ) i = 1 n + m ( 1 + α k i ) - 16 π ( 1 - t ) i , j = 1 , i j n + m ( 1 + α k i ) ( 1 + α k j ) G t , ε ( p k j , p k i ) - 16 π ( 1 - t ) log ε i = 1 n + m ( 1 + α k i ) + O ( ρ ) ,

which derives the asymptotical expansion (70) by (5), (67), and (68). □

4 Proofs of theorems

In this section, we carry out the proofs of Theorems 1.1 and 1.2 basing on the finite dimensional reduction. Now we introduce the definition of C0-stable critical point of the function φ n,m just like in [28, 36, 39].

Definition 4.1. We say that p is a C0-stable critical point of φ n,m in (Ω' \ Γ)m \ Δ m , which says that if for any sequence of the functions ψ j such that ψ j φ n,m uniformly on the compact subsets of (Ω' \ Γ)m \ Δ m , ψ j has a critical point ξ j such that ψ j (ξ j ) → + φ n,m .

In particular, if p is a strict local maximum or minimum point of φ n,m , p is a C0-stable critical point of φ n,m .

Proof of Theorem 1.1. Case (i): m ≥ 1 and 0 ≤ n ≤ N. Let

v y = V p y + ϕ p y , y Ω ̄ ε ,

where ϕ is the unique solution of the problem (55), which is established in Proposition 3.2. From Proposition 3.3, v(y) is a solution of the Equation (16), namely ϕ ̃ p x is a solution of the Equation (3) if p ε = ( p k n + 1 , ε , , p k n + m , ε ) satisfying the relation (33) is a critical point of the function F ε,t (p) with t (0, 1]. This implies that we only need to find a critical point p ε of the following function in (Ω' \ Γ)m \ Δ m

F ̃ ε , t p = F ε , 1 p - β , for t = 1 , F ε , t p - 1 - t β , for t 0 , 1 ,
(75)

where

β = 8 π m log 8 + 8 π i = 1 n ( 1 + α k i ) log 8 ( 1 + α k i ) 2 f ( p k i ) Q k i ( p k i ) - 16 π i , j = 1 , i j n ( 1 + α k i ) ( 1 + α k j ) G 1 ( p k j p k i ) - 16 π log ε i = 1 n + m ( 1 + α k i ) - 16 π i = 1 n ( 1 + α k i ) 2 H 1 ( p k i , p k i ) - 16 π i = 1 n + m ( 1 + α k i ) .

From Propositions 3.3 and 3.4, it follows that

F ̃ ε , t ( p ) = θ ε , 1 ( p ) + ρ Θ ε ( p ) - 16 π φ n , m ( p ) , for t = 1 , θ ε , t ( p ) + ρ Θ ε ( p ) - 16 π ( 1 - t ) φ n , m ( p ) , for t ( 0 , 1 ) ,
(76)

where θ ε,t (p) = O(ρ2| log ε|) for any t (0, 1], and Θ ε (p) is uniformly bounded on any compact subset S of (Ω' \ Γ)m\ Δ m for ε small. Then F ̃ ε , 1 ( p ) -16π φ n , m ( p ) , and F ̃ ε , t ( p ) -16π ( 1 - t ) φ n , m ( p ) for any t (0,1), uniformly on for ε small. By Definition 4.1, there exists a critical point p ε = ( p k n + 1 , ε , , p k n + m , ε ) of the function F ̃ ε , t such that F ̃ ε , 1 ( p ε ) - 16 π φ n , m ( p * ) , and F ̃ ε , t ( p ε ) -16π ( 1 - t ) φ n , m ( p * ) for any t (0,1). Moreover, up to a subsequence, there exists p= ( p k n + 1 , , p k n + m ) ( Ω \ Γ ) m \ Δ m such that p ε p for ε small, and φ n,m (p*) = φ n,m (p). Hence, u ε =U ( p ε ) + ϕ ̃ ( p ε ) is a family of solutions of the Equation (3). As a consequence, from the related properties of U(p ε ) and ϕ ̃ ( p ε ) , we easily know that for any λ > 0, u ε is uniformly bounded on Ω\ i = 1 n + m B λ ( p k i ) , and sup B λ ( p k i ) u ε + for ε small.

Finally, we show that u ε satisfies the concentration property:

I = ε 2 Ω x - p 1 2 α 1 x - p N 2 α N f ( x ) e u ε Ψ 8 π i = 1 n + m ( 1 + α k i ) Ψ ( p k i ) , Ψ C ( Ω ̄ ) ,
(77)

for ε small. In fact, using the inequality |es - 1|e|s||s| for any s and the estimate (56), we obtain

I = ε 2 Ω x - p 1 2 α 1 x - p N 2 α N f ( x ) e U ( p ε ) ( x ) Ψ + o ( 1 ) = Ω ε W ( p ε ) ( y ) Ψ ( ε y ) + o ( 1 ) .

Then from the asymptotical expression (21) and (22) of W(p ε )(y), we can easily get

Ω ε W ( p ε ) ( y ) Ψ ( ε y ) = 8 π i = 1 n + m ( 1 + α k i ) Ψ ( p k i ) + o ( 1 ) ,

which implies (77).

Case (ii): m = 0 and 1 ≤ n ≤ N. From Proposition 3.2 (ii), we find that u ε ( x ) =U ( p ) ( x ) +ϕ ( 1 ε x ) is a family of solutions of the Equation (3) with p= ( p k 1 , , p k n ) . Then we can get the needed multiple blowing-up and concentrating properties of u ε through the similar proof of Case (i). □

In order to give the proof of Theorem 1.2, we need a version of the maximum principle under Dirichlet-Robin boundary conditions, which is the extension of the corresponding one with respect to Dirichlet or Robin boundary condition only.

Lemma 4.2 Assume that T Ω is a relatively closed subset, b > 0 is a smooth function on T,F:T is a smooth function. If u is a solution of the equation

Δ u = 0 , i n Ω , u ν + λ b ( x ) u = F , o n T , u = 0 , o n Ω \ T ,

where λ > 0, there exists a constant C(b) > 0 only depending on b(x) such that

u L ( Ω ) + dist ( x , Ω ) u L ( Ω ) C ( b ) λ F L ( T ) .

Proof. The proof is similar to that of Lemma 2.6 in [33]. □

Proof of Theorem 1.2. Using the maximum principle with Dirichlet-Robin boundary conditions instead of Robin boundary condition if necessary (see Lemma 4.2), the proof can be similarly given out through that of Theorem 1.1. □

References

  1. Aubin T: Some Nonlinear Problems in Riemannian Geometry. Springer-Verlag, Berlin; 1998.

    Book  Google Scholar 

  2. Caglioti E, Lions PL, Marchioro C, Pulvirenti M: A special class of stationary flows for two-dimensional Euler equations: a statistical mechanics description. Commun Math Phys 1992, 143: 501-525. 10.1007/BF02099262

    Article  MathSciNet  Google Scholar 

  3. Caglioti E, Lions PL, Marchioro C, Pulvirenti M: A special class of stationary flows for two-dimensional Euler equations: a statistical mechanics description II. Commun Math Phys 1995, 174: 229-260. 10.1007/BF02099602

    Article  MathSciNet  Google Scholar 

  4. Kiessling MKH: Statistical mechanics of classical particles with logarithmic interactions. Commun Pure Appl Math 1993, 46: 27-56. 10.1002/cpa.3160460103

    Article  MathSciNet  Google Scholar 

  5. Chen WX, Li CM: Classification of solutions of some nonlinear elliptic equations. Duke Math J 1991, 63: 615-622. 10.1215/S0012-7094-91-06325-8

    Article  MathSciNet  Google Scholar 

  6. Rey O: The role of the Green's function in a nonlinear elliptic equation involving the critical Sobolev exponent. J Funct Anal 1990, 89: 1-52. 10.1016/0022-1236(90)90002-3

    Article  MathSciNet  Google Scholar 

  7. Suzuki T: Two-dimensional Emden-Fowler equation with exponential nonlinearity. In Nonlinear Diffusion Equations and Their Equilibrium States. Volume 3. Birkhäuser Boston, Boston, MA; 1992:493-512. (Gregynog, 1989), Progr. Nonlinear Differential Equations Appl., 7

    Chapter  Google Scholar 

  8. Tarantello G: Multiple condensate solutions for the Chern-Simons-Higgs theory. J Math Phys 1996, 37: 3769-3796. 10.1063/1.531601

    Article  MathSciNet  Google Scholar 

  9. Tarantello G: Selfdual Gauge Field Vortices, An Analytical Approach. In Progress in Nonlinear Differential Equations and Their Applications. Volume 72. Birkhäuser Boston, Inc., Boston, MA; 2008.

    Google Scholar 

  10. Tarantello G: On some elliptic problems in the study of selfdual Chern-Simons vortices. In Geometric analysis and PDEs. Springer, Dordrecht; 2009:117-175. Lecture Notes in Math., 1977

    Chapter  Google Scholar 

  11. Tarantello G: Analytical, geometrical and topological aspects of a class of mean field equations on surfaces. Discrete Continuous Dyn Syst A 2010, 28: 931-973.

    Article  MathSciNet  Google Scholar 

  12. Bebernes J, Eberly D: Mathematical problems from combustion theory. Springer-Verlag, New York; 1989.

    Book  Google Scholar 

  13. Chandrasekhar S: An Introduction to The Study of Stellar Structure. Dover Publications, Inc., New York; 1957.

    Google Scholar 

  14. Gel'fand IM: Some problems in the theory of quasilinear equations. Am Math Soc Transl 1963, 29: 295-381.

    Google Scholar 

  15. Murray JD: Mathematical biology. Springer-Verlag, Berlin; 1989.

    Book  Google Scholar 

  16. Mignot F, Murat F, Puel JP: Variation d'un point de retournement par rapport au domaine. Commun Partial Diff Equ 1979, 4: 1263-1297. 10.1080/03605307908820128

    Article  MathSciNet  Google Scholar 

  17. Bartolucci D, Tarantello G: Liouville type equations with singular data and their applications to periodic multivortices for the electroweak theory. Commun Math Phys 2002, 229: 3-47. 10.1007/s002200200664

    Article  MathSciNet  Google Scholar 

  18. Bartolucci D, Montefusco E: On the shape of blow-up solutions to a mean field equation. Nonlinearity 2006, 19: 611-631. 10.1088/0951-7715/19/3/005

    Article  MathSciNet  Google Scholar 

  19. Brezis H, Merle F: Uniform estimates and blow-up behavior for solutions of - u = V ( x ) eu in two dimensions. Comm Partial Diff Equ 1991, 16: 1223-1253. 10.1080/03605309108820797

    Article  MathSciNet  Google Scholar 

  20. Esposito P: Blowup solutions for a Liouville equation with singular data. SIAM J Math Anal 2005, 36: 1310-1345. 10.1137/S0036141003430548

    Article  MathSciNet  Google Scholar 

  21. Li YY, Shafrir I: Blow-up analysis for solutions of - u = Veuin dimension two. Indiana Univ Math J 1994, 43: 1255-1270. 10.1512/iumj.1994.43.43054

    Article  MathSciNet  Google Scholar 

  22. Ma L, Wei JC: Convergence for a Liouville equation. Comment Math Helv 2001, 76: 506-514. 10.1007/PL00013216

    Article  MathSciNet  Google Scholar 

  23. Nagasaki K, Suzuki T: Asymptotic analysis for two-dimensional elliptic eigenvalue problems with exponentially dominated nonlinearities. Asymptotic Anal 1990, 3: 173-188.

    MathSciNet  Google Scholar 

  24. Baraket S, Parcard F: Construction of singular limits for a semilinear elliptic equation in dimension 2. Calc var Partial Diff Equ 1998, 6: 1-38.

    Article  Google Scholar 

  25. Ambrosetti A, Malchiodi A: Perturbation methods and semilinear elliptic problems on Rn. Birkhäuser Verlag, Basel; 2006.

    Google Scholar 

  26. del Pino M, Kowalczyk M, Musso M: Singular limits in Liouville-type equations. Calc var Partial Diff Equ 2005, 24: 47-81. 10.1007/s00526-004-0314-5

    Article  MathSciNet  Google Scholar 

  27. Esposito P, Grossi M, Pistoria A: On the existence of blowing-up solutions for a mean field equation. Ann Inst H Poincaré Anal Non Linéaire 2005, 22: 227-257. 10.1016/j.anihpc.2004.12.001

    Article  Google Scholar 

  28. Zhao CY: Singular limits in a Liouville-type equation with singular sources. Houston J Math 2008, 34: 601-621.

    MathSciNet  Google Scholar 

  29. Costabel M, Dauge M: A singularly perturbed mixed boundary value problem. Com-mun Partial Diff Equ 1996, 21: 1919-1949. 10.1080/03605309608821249

    Article  MathSciNet  Google Scholar 

  30. Pinchover Y, Saadon T: On positivity of solutions of degenerate boundary value problems for second-order elliptic equations. Israel J Math 2002, 132: 125-168. 10.1007/BF02784508

    Article  MathSciNet  Google Scholar 

  31. Chen WX, Li CM: Qualitative properties of solutions to some nonlinear elliptic equations in R2. Duke Math J 1993, 71: 427-439. 10.1215/S0012-7094-93-07117-7

    Article  MathSciNet  Google Scholar 

  32. Prajapat J, Tarantello G:On a class of elliptic problems in 2 : symmetry and uniqueness results. Proc Roy Soc Edinburgh Sect A 2001, 131: 967-985. 10.1017/S0308210500001219

    Article  MathSciNet  Google Scholar 

  33. Dávila J, Kowalczyk M, Montenegro M: Critical points of the regular part of the harmonic Green function with Robin boundary condition. J Funct Anal 2008, 255: 1057-1101. 10.1016/j.jfa.2007.11.023

    Article  MathSciNet  Google Scholar 

  34. Gilbarg D, Trudinger NS: Elliptic Partial Differential Equations of Second Order. Springer-Verlag, Berlin; 2001. Reprint of the 1998 edition

    Google Scholar 

  35. Chen CC, Lin CS: Sharp Estimates for Solutions of Multi-Bubbles in Compact Riemann Surfaces. Commun Pure Appl Math 2002, 55: 728-771. 10.1002/cpa.3014

    Article  Google Scholar 

  36. Esposito P, Pistoia A, Wei J:Concentrationg solutions for the Hénon equation in 2 . J Anal Math 2006, 100: 249-280. 10.1007/BF02916763

    Article  MathSciNet  Google Scholar 

  37. Protter MH, Weinberger HF: Maximum Principles in Differential Equations. Springer-Verlag, New York; 1984.

    Book  Google Scholar 

  38. Daners D: Robin boundary value problems on arbitrary domains. Trans Am Math Soc 2000, 352: 4207-4236. 10.1090/S0002-9947-00-02444-2

    Article  MathSciNet  Google Scholar 

  39. Wei JC, Wei L, Zhou F: Mixed interior and boundary nodal bubbling solutions for a sinh-Poisson equation. Pacific J Math 2011, 250: 225-256. 10.2140/pjm.2011.250.225

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

This study was supported by N. N. S. F. C. (11171214).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yibin Chang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All authors typed, read and approved the final manuscript.

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Chang, Y., Yang, H. Multiple blowing-up and concentrating solutions for Liouville-type equations with singular sources under mixed boundary conditions. Bound Value Probl 2012, 33 (2012). https://doi.org/10.1186/1687-2770-2012-33

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2012-33

Keywords