Research

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

Yibin Chang* and Haitao Yang

### Author affiliations

Department of Mathematics, Zhejiang University, Hangzhou 310027, China

For all author emails, please log on.

Boundary Value Problems 2012, 2012:33  doi:10.1186/1687-2770-2012-33

 Received: 17 October 2011 Accepted: 23 March 2012 Published: 23 March 2012

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

### 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, ..., pN} ⊂ Ω 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.

##### Keywords:
multiple blowing-up and concentrating solution; Liouville-type equation; singular source; mixed boundary conditions; finite dimensional reduction

### 1 Introduction

- Δ 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, ..., pN} ⊂ 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 [2-4], Chern-Simons vortex theory [5-11] and several other fields of applied mathematics [12-16]. 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 Gt,ε 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 Ht,ε(x, y) = Gt,ε(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, ..., pm) ∈ m: pi = 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 = 1 with inff > 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 ≤ n N, 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,17-23]). 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,m to get the existence of multiple blowing-up and concentrating solutions for the Equation (3)|t = 1 with Γ = , namely αi = 0 for any i = 1,2, ..., N. When { α 1 , , α N } ( 0 , + ) \ , a similar result for the Equation (3)|t = 1 under 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<1 brings 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, Qi(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, ..., kn} ⊂ {1, ..., N} and kn+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 zt(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 Li under suitable dilations and translations, which have some well-known properties that any bounded solution of Ltϕ = 0 is

- a linear combination of zi0 and zij for i = n + 1, ..., n + m, j = 1, 2 (see [24,35]);

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

Remark 3.1 These properties of the operator Li 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 Li 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 θ ,

un(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 r R0, and χ(r) = 0 if r R0 + 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 ≤ n N, we find a function ϕ and scalars cij, 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 ≤ 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 ', 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 Zi0, i = 1, ..., n + m, and Zij, 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 Liϕ = 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, Ra is sufficiently large. As a result, by (42) and (43), for any R ≥ Ra, 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 hk with ║hkn,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), ║ϕkl κ or ║ϕko κ for some κ > 0. Briefly set ε:= εk, p k i : p k i k . If ║ϕkl κ, 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 BR(0). Obviously, for any q 1 , - 1 α we easily get ρ k i ε 2 ĥ k 0 in Lq(BR(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(BR(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 zi0 for i = 1, ..., n, or a linear combination of zi0 and zij 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 ║ϕko ≥ κ and ║ϕkl → 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

Δ