arxiv: v1 [math.ap] 25 Jul 2012

Transcription

1 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN: REGULARITY UP TO THE BOUNDARY XAVIER ROS-OTON AND JOAQUIM SERRA arxiv: v1 [math.ap] 25 Jul 2012 Abstract. We study the regularity up to the boundary of solutions to the Dirichlet problem for the fractional Laplacian. We prove that if u is a solution of s u = g in Ω, u 0 in R n \Ω, for some s 0, 1 and g L Ω, then u is C s R n and u/δ s Ω is C α up to the boundary Ω for some α 0, 1, where δx = distx, Ω. For this, we develop a fractional analog of the Krylov boundary Harnack method. Moreover, under further regularity assumptions on g we obtain higher order Hölder estimates for u and u/δ s. Namely, the C β norms of u and u/δ s in the sets {x Ω : δx ρ} are controlled by Cρ s β and Cρ α β, respectively. These regularity results are crucial tools in our proof of the Pohozaev identity for the fractional Laplacian [19, 20]. 1. Introduction and results Let s 0, 1 and g L Ω, and consider the fractional elliptic problem { 1.1 s u = g in Ω u = 0 in R n \Ω, in a bounded domain Ω R n, where 1.2 s ux uy ux = c n,s PV dy R x y n n+2s and c n,s is a normalization constant. Problem 1.1 is the Dirichlet problem for the fractional Laplacian. There are classical results in the literature dealing with the interior regularity of s-harmonic functions, or more generally for equations of the type 1.1. However, there are few results on regularity up to the boundary. This is the topic of study of the paper. Our main result establishes the Hölder regularity up to the boundary Ω of the function u/δ s Ω, where δx = distx, Ω. For this, we develop an analog of the Krylov [17] boundary Harnack method for problem 1.1. As in Krylov s work, our proof applies also to operators with bounded Key words and phrases. Fractional Laplacian, Dirichlet problem, regularity, boundary Harnack inequality. The authors were supported by grants MTM C03-01, MTM C04-01 Spain, and 2009SGR345 Catalunya. 1

2 2 XAVIER ROS-OTON AND JOAQUIM SERRA measurable coefficients, more precisely those of the type 1.5. This will be treated in a future work [21]. In this paper we only consider the constant coefficient operator s, since in this case we can establish more precise regularity results. Most of them will be needed in our subsequent work [20], where we find and prove the Pohozaev identity for the fractional Laplacian, announced in [19]. For 1.1, in addition to the Hölder regularity up to the boundary for u/δ s, we prove that any solution u is C s R n. Moreover, when g is not only bounded but Hölder continuous, we obtain better interior Hölder estimates for u and u/δ s. The Dirichlet problem for the fractional Laplacian 1.1 has been studied from the point of view of probability, potential theory, and PDEs. The closest result to the one in our paper is that of Bogdan [2], establishing a boundary Harnack inequality for nonnegative s-harmonic functions. It will be described in more detail later on in the Introduction in relation with Theorem 1.2. Related regularity results up to the boundary have been proved in [16] and [7]. In [16] it is proved that u/δ s has a limit at every boundary point when u solves the homogeneous fractional heat equation. The same is proven in [7] for a free boundary problem for the fractional Laplacian. Some other results dealing with various aspects concerning the Dirichlet problem are the following: estimates for the heat kernel of the parabolic version of this problem and for the Green function, e.g., [3, 10]; an explicit expression of the Poisson kernel for a ball [18]; and the explicit solution to problem 1.1 in a ball for g 1 [13]. In addition, the interior regularity theory for viscosity solutions to nonlocal equations with bounded measurable coefficients is developed in [9]. The first result of this paper gives the optimal Hölder regularity for a solution u of 1.1. The proof, which is given in Section 2, is based on two ingredients: a suitable upper barrier, and the interior regularity results for the fractional Laplacian. Given g L Ω, we say that u is a solution of 1.1 when u H s R n is a weak solution see Definition 2.1. When g is continuous, the notions of weak solution and of viscosity solution agree; see Remark We recall that a domain Ω satisfies the exterior ball condition if there exists a positive radius ρ 0 such that all the points on Ω can be touched by some exterior ball of radius ρ 0. Proposition 1.1. Let Ω be a bounded Lipschitz domain satisfying the exterior ball condition, g L Ω, and u be a solution of 1.1. Then, u C s R n and u C s R n C g L Ω, where C is a constant depending only on Ω and s. This C s regularity is optimal, in the sense that a solution to problem 1.1 is not in general C α for any α > s. This can be seen by looking at the problem { 1.3 s u = 1 in B r x 0 u = 0 in R n \B r x 0,

3 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 3 for which its solution is explicit. For any r > 0 and x 0 R n, it is given by [13, 3] 2 2s Γn/2 1.4 ux = Γ n+2s r 2 x x 0 2 s in B r x 0. 2 Γ1 + s It is clear that this solution is C s up to the boundary but it is not C α for any α > s. Since solutions u of 1.1 are C s up to the boundary, and not better, it is of importance to study the regularity of u/δ s up to Ω. For instance, our recent proof [20, 19] of the Pohozaev identity for the fractional Laplacian uses in a crucial way that u/δ s is Hölder continuous up to Ω. This is the main result of the present paper and it is stated next. For local equations of second order with bounded measurable coefficients and in non-divergence form, the analog result is given by a theorem of N. Krylov [17], which states that u/δ is C α up to the boundary for some α 0, 1. This result is the key ingredient in the proof of the C 2,α boundary regularity of solutions to fully nonlinear elliptic equations F D 2 u = 0 see [15, 6]. For our nonlocal equation 1.1, the corresponding result is the following. Theorem 1.2. Let Ω be a bounded C 1,1 domain, g L Ω, u be a solution of 1.1, and δx = distx, Ω. Then, u/δ s Ω can be continuously extended to Ω. Moreover, we have u/δ s C α Ω and u/δ s C α Ω C g L Ω for some α > 0 satisfying α < min{s, 1 s}. The constants α and C depend only on Ω and s. To prove this result we use the method of Krylov see [15]. It consists of trapping the solution between two multiples of δ s in order to control the oscillation of the quotient u/δ s near the boundary. For this, we need to prove, among other things, that s δ0 s is bounded in Ω, where δ 0 x = distx, R n \ Ω is the distance function in Ω extended by zero outside. This will be guaranteed by the assumption that Ω is C 1,1. To our knowledge, the only previous results dealing with the regularity up to the boundary for solutions to 1.1 or its parabolic version were the ones by K. Bogdan [2] and S. Kim and K. Lee [16]. The first one [2] is the boundary Harnack principle for nonnegative s-harmonic functions, which reads as follows: assume that u and v are two nonnegative functions in a Lipschitz domain Ω, which satisfy s u 0 and s v 0 in Ω B r x 0 for some ball B r x 0 centered at x 0 Ω. Assume also that u v 0 in B r x 0 \ Ω. Then, the quotient u/v is C α B r/2 x 0 for some α 0, 1. In [4] the same result is proven in open domains Ω, without any regularity assumption. While the result in [4] assumes no regularity on the domain, we need to assume Ω to be C 1,1. This assumption is needed to compare the solutions with the function δ s. As a counterpart, we allow nonzero right hand sides g L Ω and also changingsign solutions. In C 1,1 domains, our results in Section 3 which are local near any

4 4 XAVIER ROS-OTON AND JOAQUIM SERRA boundary point extend Bogdan s result. For instance, assume that u and v satisfy s u = g and s v = h in Ω, u v 0 in R n \ Ω, and that h is positive in Ω. Then, by Theorem 1.2 we have that u/δ s and v/δ s are C α Ω functions. In addition, by the Hopf lemma for the fractional Laplacian we find that v/δ s c > 0 in Ω. Hence, we obtain that the quotient u/v is C α up to the boundary, as in Bogdan s result for s-harmonic functions. As in Krylov s result, our method can be adapted to the case of nonlocal elliptic equations with bounded measurable coefficients. Namely, in another paper [21] we will prove the boundary Harnack principle for solutions to Lu = g in Ω, u 0 in R n \ Ω, where g L Ω, 2ux ux + y ux y 1.5 Lux = dy, R n y T Axy n+2s 2 and Ax is a symmetric matrix, measurable in x, and with 0 < λid Ax ΛId. A second result for the parabolic problem related to ours is contained in [16]. The authors show that any solution of t u + s u = 0 in Ω, u 0 in R n \ Ω, satisfies the following property: for any t > 0 the function u/δ s is continuous up to the boundary Ω. Our results were motivated by the study of nonlocal semilinear problems s u = fu in Ω, u 0 in R n \ Ω, more specifically, by the Pohozaev identity that we establish in [20]. Its proof requires the precise regularity theory up to the boundary developed in the present paper see Corollary 1.6 below. Other works treating the fractional Dirichlet semilinear problem, which deal mainly with existence of solutions and symmetry properties, are [22, 23, 12, 1]. In the semilinear case, g = fu and therefore g automatically becomes more regular than just bounded. When g has better regularity, the next two results improve the preceding ones. The proofs of these results require the use of the following weighted Hölder norms, a slight modification of the ones in Gilbarg-Trudinger [14, Section 6.1]. Throughout the paper, and when no confusion is possible, we use the notation C β U with β > 0 to refer to the space C k,β U, where k is the is greatest integer such that k < β and where β = β k. This notation is specially appropriate when we work with s in order to avoid the splitting of different cases in the statements of regularity results. According to this, [ ] C β U denotes the C k,β U seminorm [u] C β U = [u] C k,β U = sup D k ux D k uy. x,y U, x y x y β Moreover, given an open set U R n with U, we will also denote d x = distx, U and d x,y = min{d x, d y }.

5 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 5 Definition 1.3. Let β > 0 and σ β. Let β = k + β, with k integer and β 0, 1]. For w C β U = C k,β U, define the seminorm [w] σ = sup d β+σ D k wx D k wy x,y. x,y U x y β For σ > 1, we also define the norm σ k while for 1 < σ < 0, w σ = l=0 sup x U w σ = w C σ U + as follows: in case that σ 0, + [w] σ, d l+σ x D l wx k l=1 sup x U d l+σ x D l wx + [w] σ. Note that σ is the rescale order of the seminorm [ ] σ, in the sense that [wλ ]σ /λ = λ σ [w] σ. When g is Hölder continuous, the next result provides optimal estimates for higher order Hölder norms of u up to the boundary. Proposition 1.4. Let Ω be a bounded domain, and β > 0 be such that neither β nor β + 2s is an integer. Let g C β Ω be such that g s β;ω <, and u Cs R n be a solution of 1.1. Then, u C β+2s Ω and u s β+2s;ω C u C s R n + g s β;ω, where C is a constant depending only on Ω, s, and β. Next, the Hölder regularity up to the boundary of u/δ s in Theorem 1.2 can be improved when g is Hölder continuous. This is stated in the following theorem, whose proof uses a nonlocal equation satisfied by the quotient u/δ s in Ω see 4.2 and the fact that this quotient is C α Ω. Theorem 1.5. Let Ω be a bounded C 1,1 domain, and let α 0, 1 be given by Theorem 1.2. Let g L Ω be such that g s α α;ω <, and u be a solution of 1.1. Then, u/δ s C α Ω C γ Ω and u/δ s α γ;ω C g L Ω + g s α α;ω, where γ = min{1, α + 2s} and C is a constant depending only on Ω and s. Finally, we apply the previous results to the semilinear problem { 1.6 s u = fx, u in Ω u = 0 on R n \Ω, where Ω is a bounded C 1,1 domain and f is a Lipschitz nonlinearity.

6 6 XAVIER ROS-OTON AND JOAQUIM SERRA In the following result, the meaning of bounded solution is that of bounded weak solution see definition 2.1 or that of viscosity solution. By Remark 2.11, these two notions coincide. Also, by f C 0,1 loc Ω R we mean that f is Lipschitz in every compact subset of Ω R. Corollary 1.6. Let Ω be a bounded and C 1,1 domain, f C 0,1 loc Ω R, u be a bounded solution of 1.6, and δx = distx, Ω. Then, a u C s R n and, for every β [s, 1 + 2s, u is of class C β Ω and [u] C β {x Ω : δx ρ} Cρ s β for all ρ 0, 1. b The function u/δ s Ω can be continuously extended to Ω. Moreover, there exists α 0, 1 such that u/δ s C α Ω. In addition, for all β [α, s + α], it holds the estimate [u/δ s ] C β {x Ω : δx ρ} Cρ α β for all ρ 0, 1. The constants α and C depend only on Ω, s, f, u L R n, and β. The paper is organized as follows. In Section 2 we prove Propositions 1.1 and 1.4. In Section 3 we prove Theorem 1.2 using the Krylov method. In Section 4 we prove Theorem 1.5 and Corollary 1.6. Finally, the Appendix deals with some basic tools and barriers which are used throughout the paper. 2. Optimal Hölder regularity for u In this section we prove that, assuming Ω to be a bounded Lipschitz domain satisfying the exterior ball condition, every solution u of 1.1 belongs to C s R n. For this, we first establish that u is C β in Ω, for all β 0, 2s, and sharp bounds for the corresponding seminorms near Ω. These bounds yield u C s R n as a corollary. First, we make precise the notion of weak solution to problem 1.1. Definition 2.1. We say that u is a weak solution of 1.1 if u H s R n, u 0 a.e. in R n \ Ω, and s/2 u s/2 v dx = gv dx R n Ω for all v H s R n such that v 0 in R n \ Ω. We recall first some well known interior regularity results for linear equations involving the operator s, defined by 1.2. The first one states that w C β+2s B 1/2 whenever w C β R n and s w C β B 1. Recall that, throughout this section and in all the paper, we denote by C β, with β > 0, the space C k,β, where k is an integer, β 0, 1], and β = k + β. Proposition 2.2. Assume that w C R n solves s w = h in B 1 and that neither β nor β + 2s is an integer. Then, w C β+2s B 1/2 C w C β R n + h C β B 1,

7 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 7 where C is a constant depending only on n, s, and β. Proof. Follow the proof of Proposition in [24], where the same result is proved with B 1 and B 1/2 replaced by the whole R n. The second result states that w C β B 1/2 for each β 0, 2s whenever w L R n and s w L B 1. Proposition 2.3. Assume that w C R n solves s w = h in B 1. Then, for every β 0, 2s, w C β B 1/2 C w L R n + h L B 1, where C is a constant depending only on n, s, and β. Proof. Follow the proof of Proposition in [24], where the same result is proved in the whole R n. The third result is the analog of the first, with the difference that it does not need to assume w C β R n, but only w C β B 2 and 1 + x n 2s wx L 1 R n. Corollary 2.4. Assume that w C R n is a solution of s w = h in B 2, and that neither β nor β + 2s is an integer. Then, w C β+2s B 1/2 1 C + x n 2s wx L 1 R n + w C β B 2 + h C β B 2 where the constant C depends only on n, s, and β. Proof. Let η C R n be such that η 0 outside B 2 and η 1 in B 3/2. Then w := wη C R n and s w = h := h s w1 η. Note that for x B 3/2 we have s w1 η y w1 η x = c n,s dy. R n \B 3/2 x y n+2s From this expression we obtain that and for all γ 0, β], s w1 η L B 1 C 1 + y n 2s wy L 1 R n [ s w1 η] C γ B 1 C 1 + y n 2s γ wy L 1 R n C 1 + y n 2s wy L 1 R n for some constant C that depends only on n, s, β, and η. Therefore while we also clearly have h C β B 1 C h C β B x n 2s wx L 1 R n, w C β R n C w C β B 2. The constants C depend only on n, s, β and η. Now, we finish the proof by applying Proposition 2.2 with w replaced by w.

8 8 XAVIER ROS-OTON AND JOAQUIM SERRA Finally, the fourth result is the analog of the second one, but instead of assuming w L R n, it only assumes w L B 2 and 1 + x n 2s wx L 1 R n. Corollary 2.5. Assume that w C R n is a solution of s w = h in B 2. Then, for every β 0, 2s, w C β B 1/2 1 C + x n 2s wx L 1 R n + w L B 2 + h L B 2 where the constant C depends only on n, s, and β. Proof. Analog to the proof of Corollary 2.4. As a consequence of the previous results we next prove that every solution u of 1.1 is C s R n. First let us find an explicit upper barrier for u to prove that u Cδ s in Ω. This is the first step to obtain the C s regularity. To construct this we will need the following result, which is proved in the Appendix. Lemma 2.6 Supersolution. There exist C 1 > 0 and a radial continuous function ϕ 1 Hloc s Rn satisfying s ϕ 1 1 in B 4 \ B 1 ϕ 1 0 in B ϕ 1 C 1 x 1 s in B 4 \ B 1 1 ϕ 1 C 1 in R n \ B 4. The upper barrier for u will be constructed by scaling and translating the supersolution from Lemma 2.6. The conclusion of this barrier argument is the following. Lemma 2.7. Let Ω be a bounded domain satisfying the exterior ball condition and let g L Ω. Let u be the solution of 1.1. Then, ux C g L Ωδ s x for all x Ω, where C is a constant depending only on Ω and s. In the proof of Lemma 2.7 it will be useful the following Claim 2.8. Let Ω be a bounded domain and let g L Ω. Let u be the solution of 1.1. Then, u L R n Cdiam Ω 2s g L Ω where C is a constant depending only on n and s. Proof. The domain Ω is contained in a large ball of radius diam Ω. Then, by scaling the explicit supersolution for the ball given by 1.4 we obtain the desired bound. We next give the

9 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 9 Proof of Lemma 2.7. Since Ω satisfies the exterior ball condition, there exists ρ 0 > 0 such that every point of Ω can be touched from outside by a ball of radius ρ 0. Then, by scaling and translating the supersolution ϕ 1 from Lemma 2.6, for each of this exterior tangent balls B ρ0 we find an upper barrier in B 2ρ0 \ B ρ0 vanishing in B ρ0. This yields the bound u Cδ s in a ρ 0 -neighborhood of Ω. By using Claim 2.8 we have the same bound in all of Ω. Repeating the same argument with u we find u Cδ s, as wanted. The following lemma gives interior estimates for u and yields, as a corollary, that every bounded weak solution u of 1.1 in a C 1,1 domain is C s R n. Lemma 2.9. Let Ω be a bounded domain satisfying the exterior ball condition, g L Ω, and u be the solution of 1.1. Then, u C β Ω for all β 0, 2s and for all x 0 Ω we have the following seminorm estimate in B R x 0 = B δx0 /2x 0 : 2.2 [u] C β B R x 0 CRs β g L Ω, where C is a constant depending only on Ω, s, and β. Proof. Recall that if u solves 1.1 in the weak sense and η ɛ is the standard mollifier then s u η ɛ = g η ɛ in B R for ɛ small enough. Hence, we can regularize u, obtain the estimates, and then pass to the limit. In this way we may assume that u is smooth. Note that B R x 0 B 2R x 0 Ω. Let ũy = ux 0 + Ry. We have that 2.3 s ũy = R 2s gx 0 + Ry in B 1. Furthermore, using that u C u L R n + g L Ω δ s in Ω by Lemma 2.7 we obtain 2.4 ũ L B 1 C u L R n + g L Ω R s and, observing that ũy C u L R n + g L Ω R s 1 + y s in all of R n, y n 2s ũy L 1 R n C u L R n + g L Ω R s, with C depending only on Ω and s. Next we use Corollary 2.5, which taking into account 2.3, 2.4, and 2.5, yields ũ Cβ B 1/4 C u L R n + g L Ω R s for all β 0, 2s, where C = CΩ, s, β. Finally, we observe that [u] Cβ B R/4 x 0 = R β [ũ] Cβ B 1/4. Hence, by an standard covering argument, we find the estimate 2.2 for the C β seminorm of u in B R x 0. We now prove the C s regularity of u.

10 10 XAVIER ROS-OTON AND JOAQUIM SERRA Proof of Proposition 1.1. By Lemma 2.9, taking β = s we obtain 2.6 ux uy C u x y s L R n + g L Ω for all x, y such that y B R x with R = δx/2. We want to show that 2.6 holds, perhaps with a bigger constant C = CΩ, s, for all x, y Ω, and hence for all x, y R n since u 0 outside Ω. Indeed, observe that after a Lipschitz change of coordinates, the bound 2.6 remains the same except for the value of the constant C. Hence, we can flatten the boundary near x 0 Ω to assume that Ω B ρ0 x 0 = {x n > 0} B 1 0. Now, 2.6 holds for all x, y satisfying x y γx n for some γ = γω 0, 1 depending on the Lipschitz map. Next, let z = z, z n and w = w, w n be two points in {x n > 0} B 1/4 0, and r = z w. Let us define z = z, z n + r, z = z, z n + r and z k = 1 γ k z + γ k z and w k = γ k w + 1 γ k w, k 0. Then, using that bound 2.6 holds whenever x y γx n, we have uz k+1 uz k C z k+1 z k s = C γ k z zγ 1 s Cγ k z z. Moreover, since x n > r in all the segment joining z and w, splitting this segment into a bounded number of segments of length less than γr, we obtain Therefore, uz uw k 0 as wanted. C k 0 u z u w C z w s Cr s. uz k+1 uz k + u z u w + k 0 γ k r s + Cr u L Rn s + g L Ω C u L R n + g L Ω z w s, uw k+1 uw k The following lemma is similar to Proposition 2.2 but it involves the weighted norms introduced above. It will be used to prove Proposition 1.4 and Theorem 1.5. Lemma Let s and α belong to 0, 1, and β > 0. Let U be an open set with nonempty boundary. Assume that neither β nor β + 2s is an integer, and α < 2s. Then, 2.7 w α β+2s;u C w C α R n + s w 2s α for all w with finite right hand side. The constant C depends only on n, s, α, and β.

11 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 11 Proof. Step 1. We first control the C β+2s norm of w in balls B R x 0 with R = d x0 /2. Let x 0 U and R = d x0 /2. Define wy = wx 0 + Ry wx 0 and note that and This is because w C α B 1 R α [w] C α R n 1 + y n 2s wy L 1 R n Cn, sr α [w] C α R n. and α < 2s. Note also that wy = wx 0 + Ry wx 0 R α y α [w] C α R n s w C β B 1 = R2s+β s w C β B R x 0 Rα s w 2s α. Therefore, using Corollary 2.4 we obtain that w C β+2s B 1/2 CRα [w] C α R n + s w 2s α where the constant C depends only on n, s, α, and β. Scaling back we obtain k R l α D l w L B R/2 x 0 + R 2s+β α [w] C β+2s B R/2 x l=1 C w C α R n + s w 2s α α;u, where k denotes the greatest integer less that β + 2s and C = Cn, s. This bound holds, with the same constant C, for each ball B R x 0, x 0 U, where R = d x0 /2. Step 2. Next we claim that if 2.8 holds for each ball B dx/2x, x U, then 2.7 holds. It is clear that this already yields 2.9 k l=1 d k α x sup D k ux C x U, w C α R n + s w 2s α where k is the greatest integer less than β + 2s. To prove this claim we only have to control [w] α β+2s;u γ 0, 1 be such that β + 2s = k + γ. We next bound D k wx D k wy x y γ see Definition 1.3. Let when d x d y and x y d x /2. This will yield the bound for [w] α β+2s;u, because if x y < d x /2 then y B dx/2x, and that case is done in Step 1. We proceed differently in the cases k = 0 and k 1. If k = 0, then If k 1, then d β+2s α wx wy x = x y 2s+β d β+2s α D k wx D k wy x x y γ dx x y β+2s α wx wy x y α C w C α R n. γ dx d β+2s α γ x D k wx D k wy C w α k;u x y,

12 12 XAVIER ROS-OTON AND JOAQUIM SERRA where we have used that β + 2s α γ = k α. Finally, noting that for x B R x 0 we have R d x0 2.8, 2.9 and the definition of w α α+2s;u in 1.3. Finally, to end this section, we prove Proposition R, 2.7 follows from Proof of Proposition 1.4. Set α = s in Lemma Remark When g is continuous, the notions of bounded weak solution and viscosity solution of 1.1 and hence of 1.6 coincide. Indeed, let u H s R n be a weak solution of 1.1. Then, from Proposition 1.1 it follows that u is continuous up to the boundary. Let u ε and g ε be the standard regularizations of u and g by convolution with a mollifier. It is immediate to verify that, for ε small enough, we have s u ε = g ε in every subdomain U Ω in the classical sense. Then, noting that u ε u and g ε g locally uniformly in Ω, and applying the stability property for viscosity solutions [9, Lemma 4.5], we find that u is a viscosity solution of 1.1. Conversely, every viscosity solution of 1.1 is a weak solution. This follows from three facts: the existence of weak solution, that this solution is a viscosity solution as shown before, and the uniqueness of viscosity solutions [9, Theorem 5.2]. As a consequence of this, if g is continuous, any viscosity solution of 1.1 belongs to H s R n since it is a weak solution. This fact, which is not obvious, can also be proved without using the result on uniqueness of viscosity solutions. Indeed, it follows from Proposition 1.4 and Lemma 4.4, which yield a stronger fact: that s/2 u L p R n for all p <. Note that although we have proved Proposition 1.4 for weak solutions, its proof is also valid with almost no changes for viscosity solutions. 3. Boundary regularity In this section we study the precise behavior near the boundary of the solution u to problem 1.1, where g L Ω. More precisely, we prove that the function u/δ s Ω has a C α Ω extension. This is stated in Theorem 1.2. This result will be a consequence of the interior regularity results of Section 2 and an oscillation lemma near the boundary, which can be seen as the nonlocal analog of Krylov s boundary Harnack principle; see Theorem 4.28 in [15]. The following proposition and lemma will be used to establish Theorem 1.2. They are proved in the Appendix. Proposition D solution in half space, [7]. The function ϕ 0, defined by { 0 if x ϕ 0 x = x s if x 0, satisfies s ϕ 0 = 0 in R +.

13 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 13 The lemma below gives a subsolution in B 1 \ B 1/4 whose support is B 1 R n and such that it is comparable to 1 x s in B 1. Lemma 3.2 Subsolution. There exist C 2 > 0 and a radial function ϕ 2 = ϕ 2 x satisfying s ϕ 2 0 in B 1 \ B 1/4 ϕ 2 = 1 in B 1/4 3.2 ϕ 2 x C 2 1 x s in B 1 ϕ 2 = 0 in R n \ B 1. To prove Hölder regularity of u/δ s Ω up to the boundary, we will control the oscillation of this function in sets near Ω whose diameter goes to zero. To do it, we will set up an iterative argument as it is done for second order equations. Let us define the sets in which we want to control the oscillation and also auxiliary sets that are involved in the iteration. Definition 3.3. Let κ > 0 be a fixed small constant and let κ = 1/2 + 2κ. We may take, for instance κ = 1/16, κ = 5/8. Given a point x 0 in Ω and R > 0 let us define D R = D R x 0 = B R x 0 Ω and D + κ R = D+ κ R x 0 = B κ Rx 0 {x Ω : x νx 0 2κR}, where νx 0 is the unit outward normal at x 0 ; see Figure 3.1. By C 1,1 regularity of the domain, there exists ρ 0 > 0, depending on Ω, such that the following inclusions hold for each x 0 Ω and R ρ 0 : 3.3 B κr y D R x 0 for all y D + κ R x 0, and 3.4 B 4κR y 4κRνy D R x 0 and B κr y 4κRνy D + κ R x 0 for all y D R/2, where y Ω is the unique boundary point satisfying y y = disty, Ω. Note that, since R ρ 0, y D R/2 is close enough to Ω and hence the point y 4κRνy lays on the line joining y and y ; see Remark 3.4 below. Remark 3.4. Throughout the paper, ρ 0 > 0 is a small constant depending only on Ω, which we assume to be a bounded C 1,1 domain. Namely, we assume that 3.3 and 3.4 hold whenever R ρ 0, for each x 0 Ω, and also that every point on Ω can be touched from both inside and outside Ω by balls of radius ρ 0. In other words, given x 0 Ω, there are balls of radius ρ 0, B ρ0 x 1 Ω and B ρ0 x 2 R n \ Ω, such that B ρ0 x 1 B ρ0 x 2 = {x 0 }. A useful observation is that all points y in the segment that joins x 1 and x 2 through x 0 satisfy δy = y x 0. Recall that δ = dist, Ω.

14 14 XAVIER ROS-OTON AND JOAQUIM SERRA D R D + κ R x 0 B R/2 B κ R y Ω B R Figure 3.1. The sets D R and D + κ R In the rest of this section, by s u K we mean that either s u = g in the weak sense for some g L satisfying g L K or that u satisfies K s u K in the viscosity sense. The first and main step towards Theorem 1.2 is the following. Proposition 3.5. Let Ω be a bounded C 1,1 domain, and u be such that s u K in Ω and u 0 in R n \ Ω, for some constant K. Given any x 0 Ω, let D R be as in Definition 3.3. Then, there exist α 0, 1 and C depending only on Ω and s but not on x 0 such that 3.5 sup u/δ s inf u/δ s CKR α D R D R for all R ρ 0, where ρ 0 > 0 is a constant depending only on Ω. To prove Proposition 3.5 we need three preliminary lemmas. We start with the first one, which might be seen as the fractional version of Lemma 4.31 in [15]. Recall that κ 1/2, 1 is a fixed constant throughout the section. It may be useful to regard the following lemma as a bound by below for inf DR/2 u/δ s, rather than an upper bound for inf D + u/δ s. κ R Lemma 3.6. Let Ω be a bounded C 1,1 domain, and u be such that u 0 in all of R n and s u K in D R, for some constant K. Then, there exists a positive constant C, depending only on Ω and s, such that 3.6 inf u/δ s C inf u/δ s + KR s D + κ D R/2 R for all R ρ 0, where ρ 0 > 0 is a constant depending only on Ω. Proof. Step 1. We do first the case K = 0. Let R ρ 0, and let us call m = inf D + u/δ s 0. We have u mδ s mκr s on D + κ κ R R. The second inequality is a consequence of 3.3.

15 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 15 We scale the subsolution ϕ 2 in Lemma 3.2 as follows, to use it as lower barrier: By 3.2 we have ψ R x := κr s ϕ 2 x 4κR. s ψ R 0 in B 4κR \ B κr ψ R = κr s in B κr ψ R 4 s C 2 4κR x s in B 4κR \ B κr ψ R 0 in R n \ B 4κR. Given y D R/2, we have either y D + κ R or δy < 4κR, by 3.4. If y D+ κ R it follows from the definition of m that m uy/δy s. If δy < 4κR, let y be the closest point to y on Ω and ỹ = y +4κνy. Again by 3.4, we have B 4κR ỹ D R and B κr ỹ D + κ R. But recall that u mκrs in D + κ R, s u = 0 in Ω, and u 0 in R n. Hence, ux mψ R x ỹ in all R n and in particular u/δ s 4 s C 2 m on the segment joining y and ỹ, that contains y. Therefore, 3.7 inf u/δ s C inf u/δ s. D + κ D R/2 R Step 2. If K > 0 we consider ũ to be the solution of { s ũ = 0 in D R ũ = u in R n \ D R. By Step 1, 3.7 holds with u replaced by ũ. On the other hand, w = ũ u satisfies s w K and w 0 outside D R. Recall that points of Ω can be touched by exterior balls of radius less than ρ 0. Hence, using the rescaled supersolution KR 2s ϕ 1 x/r from Lemma 2.6 as upper barrier and we readily prove, as in the proof of Lemma 2.7, that w C 1 KR s δ s in D R. Thus, 3.6 follows. The second lemma towards Proposition 3.5, which might be seen as the fractional version of Lemma 4.35 in [15], is the following. Lemma 3.7. Let Ω be a bounded C 1,1 domain, and u be such that u 0 in all of R n and s u K in D R, for some constant K. Then, there exists a positive constant C, depending on Ω and s, such that 3.8 sup u/δ s C inf u/δ s + KR s D + D + κ R κ R for all R ρ 0, where ρ 0 > 0 is a constant depending only on Ω.

16 16 XAVIER ROS-OTON AND JOAQUIM SERRA Proof. Step 1. Consider first the case K = 0. In this case 3.8 follows from the Harnack inequality for the fractional Laplacian [18] note that we assume u 0 in all R n. Indeed, by 3.3, for each y D + κ R we have B κry D R and hence s u = 0 in B κr y. Then we may cover D + κ R by a finite number of balls B κr/2 y i, using the same scaled covering for all R ρ 0, to obtain sup u C inf u. B κr/2 y i B κr/2 y i Then, 3.8 follows since κr/2 s δ s 3κR/2 s in B κr/2 y i by 3.3. Step 2. When K > 0, we prove 3.8 by using a similar argument as in Step 2 in the proof of Proposition 3.6. Before proving Lemma 3.9 we give an extension lemma see [11, Theorem 1, Section 3.1] where the case α = 1 is proven in full detail. Lemma 3.8. Let α 0, 1] and V R n a bounded domain. There exists a nonlinear map E : C 0,α V C 0,α R n satisfying Ew w in V, [Ew] C 0,α R n [w] C 0,α V, and Ew L R n w L V for all w C 0,α V. Proof. It is immediate to check that { Ewx = min min z V satisfies the conditions since, for all x, y, z in R n, {wz + [w] C αv z x α }, w L V z x α z y α + y x α. We can now give the third lemma towards Proposition 3.5. This lemma, which is related to Proposition 3.1, is crucial. It states that δ s Ω, extended by zero outside Ω, is an approximate solution in a neighborhood of Ω inside Ω. Lemma 3.9. Let Ω be a bounded C 1,1 domain, and δ 0 = δχ Ω be the distance function in Ω extended by zero outside Ω. Let α = min{s, 1 s}, and ρ 0 be given by Remark 3.4. Then, s δ s 0 belongs to C α Ω ρ0, where Ω ρ0 = Ω {δ < ρ 0 }. In particular, s δ s 0 C Ω in Ω ρ0, where C Ω is a constant depending only on Ω and s. Proof. Fix a point x 0 on Ω and denote, for ρ > 0, B ρ = B ρ x 0. Instead of proving that Rn s δ0 s δ 0 x s δ 0 y s = c n,s PV dy x y n+2s }

17 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 17 is C α Ω B ρ0 as a function of x, we may equivalently prove that 3.9 PV B2ρ0 δ 0 x s δ 0 y s This is because the difference x y n+2s dy belongs to C α Ω B ρ0. 1 s δ s δ 0 x 0 PV c n,s B2ρ0 s δ 0 y s δ 0 x dy = x y Rn\B2ρ0 s δ 0 y s dy n+2s x y n+2s belongs to C s B ρ0, since δ0 s is C s R n and x n 2s is integrable and smooth outside a neighborhood of 0. To see 3.9, we flatten the boundary. Namely, consider a C 1,1 change of variables X = Ψx, where Ψ : B 3ρ0 V R n is a C 1,1 diffeomorphism, satisfying that Ω is mapped onto {X n = 0}, Ω B 3ρ0 is mapped into R n +, and δ 0 x = X n +. Such diffeomorphism exists because we assume Ω to be C 1,1. Let us respectively call V 1 and V 2 the images of B ρ0 and B 2ρ0 under Ψ. Let us denote the points of V V by X, Y. We consider the functions x and y, defined in V, by x = Ψ 1 X and y = Ψ 1 Y. With these notations, we have and therefore x y = DΨ 1 XX Y + O X Y 2, 3.10 x y 2 = X Y T AXX Y + O X Y 3, where AX = DΨ 1 X T DΨ 1 X is a symmetric matrix, uniformly positive definite in V 2. Hence, PV δ 0 x B2ρ0 s δ 0 y s X n dy = PV x y V2 s + Y n s + gx, Y dy, n+2s X Y T AXX Y n+2s 2 where we have denoted X Y T AXX Y gx, Y = x y 2 n+2s 2 JY and J = det DΨ 1. Note that we have g C 0,1 V 2 V 2, since Ψ is C 1,1 and we have Now we are reduced to proving that 3.11 ψ 1 X := PV X n V2 s + Y n s + X Y T AXX Y n+2s 2 gx, Y dy, belongs to C α V + 1 as a function of X, where V + 1 = V 1 {X n > 0}.

18 18 XAVIER ROS-OTON AND JOAQUIM SERRA To prove this, we extend the Lipschitz function g C 0,1 V 2 V 2 to all R n. Namely, consider the function g = Eg C 0,1 R n R n provided by Proposition 3.8, which satisfies g g in V 2 V 2 and g C 0,1 R n R n g C 0,1 V 2 V 2. By the same argument as above, using that V 1 V 2, we have that ψ 1 C α V + 1 if and only if so is the function ψx = PV Rn X n s + Y n s + X Y T AXX Y n+2s 2 g X, Y dy. Furthermore, from g define g C 0,1 V 2 R n by gx, Z = g X, X+MZ det M, where M = MX = DΨX. Then, using the change of variables Y = X + MZ we deduce ψx = PV Rn X n s + e n X + MZ s + gx, ZdZ. Z n+2s Next, we prove that ψ C α R n, which concludes the proof. Indeed, taking into account that the function X n s + is s-harmonic in R n + by Proposition 3.1 we obtain Rn e X s + e X + Z s + PV dz = 0 Z n+2s for every e R n and for every X such that e X > 0. Thus, letting e = e T nm and X = M 1 X we deduce PV Rn X n s + e n X + MZ s + Z n+2s dz = 0 for every X such that e T nm M 1 X > 0, that is, for every X R n +. Therefore, it holds φx, 0 φx, Z ψx = gx, Z gx, 0 dz, R Z n n+2s where φx, Z = e n X + MZ s + satisfies [φ] C s V 2 R n C, and g C 0,1 V 2 R n C. Let us finally prove that ψ belongs to C α V 1 +. To do it, let X and X be in V 1 +. Then, we have ψx ψ X = ΘX, X, Z R Z n n+2s dz,

19 where 3.12 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 19 ΘX, X, Z = φx, 0 φx, Z gx, Z gx, 0 φ X, 0 φ X, Z g X, Z g X, 0 = φx, 0 φx, Z φ X, 0 + φ X, Z gx, Z gx, 0 φ X, 0 φ X, Z gx, Z gx, 0 g X, Z + g X, 0. Now, on the one hand, it holds 3.13 ΘX, X, Z C Z 1+s, since [φ] C s V 2 R n C and g C 0,1 V 2 R n C. On the other hand, it also holds 3.14 ΘX, X, Z C X X s min{ Z, Z s }. Indeed, we only need to observe that gx, Z gx, 0 g X, Z + g X, 0 C min { min{ Z, 1}, X X } C min{ Z 1 s, 1} X X s. Thus, letting r = X X and using 3.13 and 3.14, we obtain ψx ψ X ΘX, X, Z dz R Z n n+2s C Z Br 1+s Cr s min{ Z, Z s } dz + dz Z n+2s R n \B r Z n+2s as desired. Next we prove Proposition 3.5. Cr 1 s + C max{r 1 s, r s }, Proof of Proposition 3.5. By considering u/k instead of u we may assume that K = 1, that is, that s u 1 in Ω. Then, by Claim 2.8 we have u L R n C for some constant C depending only on Ω and s. Let ρ 0 > 0 be given by Remark 3.4. Fix x 0 Ω. We will prove that there exist constants C 0 > 0, ρ 1 0, ρ 0, and α 0, 1, depending only on Ω and s, and monotone sequences m k and M k such that, for all k 0, 3.15 M k m k = 4 αk, 1 m k m k+1 < M k+1 M k 1, and 3.16 m k C 1 0 u/δ s M k in D Rk = D Rk x 0, where R k = ρ 1 4 k. Note that 3.16 is equivalent to the following inequality in B Rk instead of D Rk recall that D Rk = B Rk Ω m k δ s 0 C 1 0 u M k δ s 0 in B Rk = B Rk x 0, where R k = ρ 1 4 k.

20 20 XAVIER ROS-OTON AND JOAQUIM SERRA If there exist such sequences, then 3.5 holds for all R ρ 1 with C = 4 α C 0 /ρ α 1. Then, by increasing the constant C if necessary, 3.5 holds also for every R ρ 0. Next we construct {M k } and {m k } by induction. By Lemma 2.7, we find that there exist m 0 and M 0 such that 3.15 and 3.16 hold for k = 0 provided we pick C 0 large enough depending on Ω and s. Assume that we have sequences up to m k and M k. We want to prove that there exist m k+1 and M k+1 which fulfill the requirements. Let 3.18 u k = C 1 0 u m k δ s 0. We will consider the positive part u + k of u k in order to have a nonnegative function in all of R n to which we can apply Lemmas 3.6 and 3.7. Let u k = u + k u k. Observe that, by induction hypothesis, 3.19 u + k = u k and u k = 0 in B R k. Moreover, C 1 0 u m j δ s 0 in B Rj for each j k. Therefore, by 3.18 we have u k m j m k δ s 0 m j M j + M k m k δ s 0 4 αj + 4 αk δ s 0 in B Rj. But clearly 0 δ s 0 R s j = ρ s 14 js in B Rj, and therefore using R j = ρ 1 4 j u k ρ α 1 R s jr α j R α k in B Rj for each j k. Thus, since for every x B R0 \ B Rk we find 3.20 u k x ρ α there is j < k such that x x 0 < R j = ρ 1 4 j 4 x x 0, 1 R α+s k 4x x 0 R k s 4x x 0 R k By 3.20 and 3.19, at x B Rk /2x 0 we have 0 s u k x = c x + y n,s dy c n,s ρ α 1 u k x+y / B Rk y n+2s = Cρ α 1 R α s k y R k /2 ε 0 ρ α 1 R α s k, z 1/2 R α+s k α 1 8y s 8y R k R k 8z s 8z α 1 z n+2s α dz outside B Rk. 1 y n 2s dy where ε 0 = ε 0 α 0 as α 0 since 8z α 1. Therefore, writing u + k = C 1 0 u m k δ0 s + u k and using Lemma 3.9, we have s u + k C 1 0 s u + m k s δ s 0 + s u k C0 1 + C Ω + ε 0 ρ α C 1 ρ s α 1 + ε 0 ρ α 1 1 R α s k R α s k in D Rk /2.

21 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 21 In the last inequality we have just used R k ρ 1 and α s. Now we can apply Lemmas 3.6 and 3.7 with u in its statements replaced by u + k, recalling that = u k = C0 1 u m k δ s in D Rk to obtain sup D + κ R k / u + k C0 1 u/δ s m k C inf D + C κ 0 1 u/δ s m k + C 1 ρ s α 1 + ε 0 ρ α 1 R α R k /2 k C inf DRk /4 C 1 0 u/δ s m k + C 1 ρ s α 1 + ε 0 ρ α 1 R α k. Next we can repeat all the argument upside down, that is, with the functions u k = M k δ s u instead of u k. In this way we obtain, instead of 3.21, the following: 3.22 sup M k C0 1 u/δ s C inf M k C0 1 u/δ s + C 1 ρ s α 1 + ε 0 ρ α 1 R α k. D + D Rk κ /4 R k /2 Adding 3.21 and 3.22 we obtain 3.23 M k m k C inf = C C 1 D Rk /4 0 u/δ s m k + inf D Rk /4 M k C0 1 u/δ s + C 1 ρ s α 1 + ε 0 ρ α 1 R α k R α k, inf C0 1 u/δ s sup C0 1 u/δ s + M k m k + C 1 ρ s α 1 + ε 0 ρ α 1 D Rk+1 D Rk+1 and thus, using that M k m k = 4 αk and R k = ρ 1 4 k, sup C0 1 u/δ s inf C0 1 u/δ s C 1 + C D Rk+1 D C 1ρ s 1 + ε 0 4 αk. Rk+1 Now we choose α and ρ 1 small enough so that C 1 C + C 1ρ s 1 + ε 0 α 4 α. This is possible since ε 0 α 0 as α 0 and the constants C and C 1 do not depend on α nor ρ 1 they depend only on Ω and s. Then, we find sup C0 1 u/δ s inf C0 1 u/δ s 4 αk+1, D Rk+1 D Rk+1 and thus we are able to choose m k+1 and M k+1 satisfying 3.15 and Finally, we give the: Proof of Theorem 1.2. Define v = u/δ s Ω and K = g L Ω. As in the proof of Proposition 3.5, by considering u/k instead of u we may assume that s u 1 in Ω and that u L Ω C for some constant C depending only on Ω and s. First we claim that there exist constants C, M > 0, α 0, 1 and β 0, 1, depending only on Ω and s, such that

22 22 XAVIER ROS-OTON AND JOAQUIM SERRA i v L Ω C. ii For all x Ω, it holds the seminorm bound [v] C β B R/2 x C 1 + R M, where R = distx, R n \ Ω. iii For each x 0 Ω and for all ρ > 0 it holds sup v inf v Cρ α. B ρx 0 Ω B ρx 0 Ω Indeed, it follows from Lemma 2.7 that v L Ω C for some C depending only on Ω and s. Hence, i is satisfied. Moreover, if β 0, 2s, it follows from Lemma 2.9 that for every x Ω, [u] C β B R/2 x CR β, β 0, 2s, where R = δx. But since Ω is C 1,1, then provided δx < ρ 0 we will have δ s L B R/2 x CR s and [δ s ] C 0,1 B R/2 x CR s 1 and hence, by interpolation, [δ s ] C β B R/2 x CR s β for each β 0, 1. Thus, since v = uδ s, we find [v] C β B R/2 x C 1 + R s β for all x Ω and β < min{1, 2s}. Therefore hypothesis ii is satisfied. The constants C depend only on Ω and s. In addition, using Proposition 3.5 and that v L Ω C, we deduce that hypothesis iii is satisfied. Now, we claim that i-ii-iii lead to [v] C α Ω C, for some α 0, 1 depending only on Ω and s. Indeed, let x, y Ω, R = distx, R n \ Ω disty, R n \ Ω, and r = x y. Let us see that vx vy Cr α for some α > 0. If r 1 then it follows from i. Assume r < 1, and let p 1 to be chosen later. Then, we have the following dichotomy: Case 1. Assume r R p /2. Let x 0, y 0 Ω be such that x x 0 = distx, R n \ Ω and y y 0 = disty, R n \ Ω. Then, using iii and the definition of R we deduce vx vy vx vx 0 + vx 0 vy 0 + vy 0 vy CR α Cr α/p. Case 2. Assume r R p /2. Hence, since p 1, we have y B R/2 x. Then, using ii we obtain vx vy C1 + R M r β C 1 + r M/p r β Cr β M/p.

23 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 23 To finish the proof we only need to choose p > M/β and take α = min{ α/p, β M/p}. 4. Interior estimates for u/δ s The main goal of this section is to prove the C γ bounds in Ω for the function u/δ s in Theorem 1.5. To prove this result we find an equation for the function v = u/δ s Ω, that is derived below. This equation is nonlocal, and thus, we need to give values to v in R n \ Ω, although we want an equation only in Ω. It might seem natural to consider u/δ s, which vanishes outside Ω since u 0 there, as an extension of u/δ s Ω. However, such extension is discontinuous through Ω, and it would lead to some difficulties. Instead, we consider a C α R n extension of the function u/δ s Ω, which is C α Ω by Theorem 1.2. Namely, throughout this section, let v be the C α R n extension of u/δ s Ω given by Lemma 3.8. Let δ 0 = δχ Ω, and note that u = vδ s 0 in R n. Then, using 1.1 we have gx = s vδ s 0 = v s δ s 0 + δ s 0 s v I s v, δ s 0 in Ω ρ0 = {x Ω : δx < ρ 0 }, where w1 x w 1 y w 2 x w 2 y 4.1 I s w 1, w 2 x = c n,s dy x y n+2s R n and ρ 0 is a small constant depending on the domain; see Remark 3.4. Here, we have used that s w 1 w 2 = w 1 s w 2 +w 2 s w 1 I s w 1, w 2, which follows easily from 1.2. This equation is satisfied pointwise in Ω ρ0, since g is C α in Ω. We have to consider Ω ρ0 instead of Ω because the distance function is C 1,1 there and thus we can compute s δ0. s In all Ω the distance function δ is only Lipschitz and hence s δ0 s is singular for s 1. 2 Thus, the following is the equation for v: 4.2 s v = 1 gx v s δ s δ0 s 0 + I s v, δ0 s in Ω ρ0. From this equation we will obtain the interior estimates for v. More precisely, we will obtain a priori bounds for the interior Hölder norms of v, treating δ0 s I s v, δ0 s as a lower order term. For this, we consider the weighted Hölder norms given by Definition 1.3. Recall that, in all the paper, we denote C β the space C k,β, where β = k + β with k integer and β 0, 1]. In Theorem 1.2 we have proved that u/δ s Ω is C α Ω for some α 0, 1, with an estimate. From this C α estimate and from the equation for v 4.2, we will find next the estimate for u/δ s α γ;ω stated in Theorem 1.5. The proof of this result relies on some preliminary results below. Next lemma is used to control the lower order term δ0 s I s v, δ0 s in the equation 4.2 for v.

24 24 XAVIER ROS-OTON AND JOAQUIM SERRA Lemma 4.1. Let Ω be a bounded C 1,1 domain, and U Ω ρ0 be an open set. Let s and α belong to 0, 1 and satisfy α + s 1 and α < s. Then, 4.3 I s w, δ0 s s α α;u C [w] C α R n + [w] α α+s;u, for all w with finite right hand side. The constant C depends only on Ω, s, and α. To prove Lemma 4.1 we need the next Lemma 4.2. Let U R n be a bounded open set. Let α 1, α 2, 0, 1 and β 0, 1] satisfy α i < β for i = 1, 2, α 1 + α 2 < 2s, and s < β < 2s. Assume that w 1, w 2 C β U. Then, 4.4 I s w 1, w 2 2s α 1 α 2 2β 2s;U C [w 1 ] C α 1 R n + [w 1 ] α 1 [w 2 ] C α 2 R n + [w 2 ] α 2 for all functions w 1, w 2 with finite right hand side. The constant C depends only on α 1, α 2, n, β, and s. Proof. Let x 0 U and R = d x0 /2, and denote B ρ = B ρ x 0. Let K = [w 1 ] C α 1 R n + [w 1 ] α 1 [w 2 ] C α 2 R n + [w 2 ] α 2., First we bound I s w 1, w 2 x 0. w1 x 0 w 1 y w2 x 0 w 2 y I s w 1, w 2 x 0 C dy R x n 0 y n+2s R α 1+α 2 2β [w 1 ] α 1 C [w 2] α 2 z 2β dz + B R 0 z n+2s [w 1 ] C α 1 R + C n [w 2 ] C α 2 R n z α 1+α 2 CR α 1+α 2 2s K. R n \B R 0 z n+2s dz Let x 1, x 2 B R/2 x 0 B 2R x 0. Next, we bound I s w 1, w 2 x 1 I s w 1, w 2 x 2. Let η be a smooth cutoff function such that η 1 on B 1 0 and η 0 outside B 3/2 0. Define x η R x0 x = η and w i = w i w i x 0 η R, i = 1, 2. R Note that we have w i L R n = w i L B 3R/2 αi 3R [w i ] 2 CαiRn

25 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 25 and [ w i ] C β R n C [w i ] C β B 3R/2 η L B 3R/2 + w i w i 0 L B 3R/2 [w i ] C β B 3R/2 CR α i β [w i ] C α ir n + [w i ] α i. Let ϕ i = w i w i x 0 w i and observe that ϕ i vanishes in B R. Hence, ϕ i x 1 = ϕ i x 2 = 0, i = 1, 2. Next, let us write I s w 1, w 2 x 1 I s w 1, w 2 x 2 = c n,s J 11 + J 12 + J 21 + J 22, where w1 x 1 w 1 y w 2 x 1 w 2 y J 11 = dy R x n 1 y n+2s w1 x 2 w 1 y w 2 x 2 w 2 y dy, R x n 2 y n+2s w1 x 1 w 1 y ϕ 2 y w1 x 2 w 1 y ϕ 2 y J 12 = + dy, R n \B R x 1 y n+2s x 2 y n+2s w2 x 1 w 2 y ϕ 1 y w2 x 2 w 2 y ϕ 1 y J 21 = + dy, R n \B R x 1 y n+2s x 2 y n+2s and ϕ 1 yϕ 2 y J 22 = R n \B R x 1 y ϕ 1yϕ 2 y dy. n+2s x 2 y n+2s We now bound separately each of these terms. Bound of J 11. We write J 11 = J J11 2 where J11 1 w1 x 1 w 1 x 1 + z w 1 x 2 + w 1 x 2 + z w 2 x 1 w 2 x 1 + z = dz, R n z n+2s J11 2 w1 x 2 w 1 x 2 + z w 2 x 1 w 2 x 1 + z w 2 x 2 + w 2 x 2 + z = dz. R z n n+2s To bound J11 1 we proceed as follows R α1 β [w J ] α 1 z β R α2 β [w 2 ] α 2 z β dz + B r0 z n+2s R α1 β [w 1 ] α 1 + rβ R α2 β [w 2 ] α 2 z β dz R n \B r0 CR α 1+α 2 2β r 2β 2s K. Similarly, J 2 11 CR α 1+α 2 2β r 2β 2s K. z n+2s

26 26 XAVIER ROS-OTON AND JOAQUIM SERRA and Bound of J 12 and J 21. We write J 12 = J J12 2 where J12 1 = ϕ 2 y w 1x 1 w 1 x 2 dy R n \B R x 1 y n+2s J12 2 = ϕ 2 y w 1 x 2 w 1 y { } 1 R n \B R x 1 y 1 dy. n+2s x 2 y n+2s To bound J12 1 we recall that ϕ 2 x 1 = 0 and proceed as follows J12 1 C x 1 y α 2 R α1 β [w 1 ] α 1 [ϕ 2 ] rβ C 0,α 2 R n dy R n \B R x 1 y n+2s CR α 1+α 2 β 2s r β K CR α 1+α 2 2β r 2β 2s K. We have used that [ϕ 2 ] C α 2 R n = [w w] C α 2 R n 2[w] C α 2 R n, r R, and β < 2s. To bound J 2 12, let Φz = z n 2s. Note that, for each γ 0, 1], we have 4.5 Φz 1 z Φz 2 z C z 1 z 2 γ z n 2s γ for all z 1, z 2 in B R/2 0 and z R n \ B R 0. Then, using that ϕ 2 x 2 = 0, J12 2 C x 2 y α 1+α 2 x 1 x 2 2β 2s [ϕ 2 ] C α 2 R n [ϕ 2 ] C α 2 R n dy R n \B R x 2 y n+2β CR α 1+α 2 2β r 2β 2s K. This proves that J 12 CR α 1+α 2 2β r 2β 2s K. Changing the roles of α 1 and α 2 we obtain the same bound for J 21. Bound of J 22. Using again ϕ i x i = 0, i = 1, 2, we write J 22 = ϕ1 x 1 ϕ 1 y ϕ 2 x 1 ϕ 2 y R n \B R 1 x 1 y n+2s 1 x 2 y n+2s Hence, using again 4.5, J 22 C x 1 y α 1+α 2 x 1 x 2 2β 2s [ϕ 2 ] C 0,α 2 R [ϕ n 2 ] C 0,α 2 R n dy R n \B R x 1 y n+2β CR α 1+α 2 2β r 2β 2s K. dy. Summarizing, we have proven that for all x 0 such that d x = 2R and for all x 1, x 2 B R/2 x 0 it holds and I s δ s 0, wx 0 CR α 1 α 2 2s K I s δ s 0, wx 1 I s δ s 0, wx 2 x 1 x 2 2β 2s CR α 1+α 2 2β [w] α α+s;u + [w] C α R n. This yields 4.4, as shown in Step 2 in the proof of Lemma 2.10.

27 THE DIRICHLET PROBLEM FOR THE FRACTIONAL LAPLACIAN 27 Next we prove Lemma 4.1. Proof of Lemma 4.1. The distance function δ 0 is C 1,1 in Ω ρ0 and since U Ω ρ0 we have d x δ 0 x for all x U. Hence, it follows that [δ s 0] C s R n + [δ s 0] s CΩ, β for all β [s, 2]. Then, applying Lemma 4.2 with w 1 = w, w 2 = δ0, s α 1 = α, α 2 = s, and β = s + α, we obtain I s w, δ0 s s α 2α;U C [w] C α R n + [w] α α+s;u, and hence 4.3 follows. Using Lemma 4.1 we can now prove Theorem 1.5 and Corollary 1.6. Proof of Theorem 1.5. Let U Ω ρ0. We prove first that there exist α 0, 1 and C, depending only on s and Ω and not on U, such that u/δ s α α+2s;u C g L Ω + g s α α;ω. Then, letting U Ω ρ0 we will find that this estimate holds in Ω ρ0 with the same constant. To prove this, note that by Theorem 1.2 we have u/δ s C α Ω C s, Ω g L Ω. Recall that v denotes the C α R n extension of u/δ s Ω given by Lemma 3.8, which satisfies v C α R n = u/δ s C α Ω. Since u Cα+2s Ω and δ C 1,1 Ω ρ0, it is clear that v α α+2s;u < it is here where we use that we are in a subdomain U and not in Ω ρ0. Next we obtain an a priori bound for this seminorm in U. To do it, we use the equation 4.2 for v: s v = 1 gx v s δ s δ s 0 + Iδ0, s v in Ω ρ0 = {x Ω : δx < ρ 0 }. Now we will se that this equation and Lemma 2.10 lead to an a priori bound for v α α+2s;u. To apply Lemma 2.10, we need to bound s v 2s α α;u. Let us examine the three terms on the right hand side of the equation. First term. Using that d x = distx, U < distx, Ω = δx for all x U we obtain that, for all α s, δ s g 2s α α;u C s, Ω g s α α;ω. Second term. We know from Lemma 3.9 that, for α min{s, 1 s}, s δ s 0 C α Ω ρ0 C s, Ω.

28 28 XAVIER ROS-OTON AND JOAQUIM SERRA Hence, δ s v s δ s 0 2s α α;u Third term. From Lemma 4.1 we know that Iv, δ0 s s α α;u Cn, s, α diamω s δ s v s δ s 0 s α α;u C s, Ω v C α R n C s, Ω g L Ω. v C α R n + [v] α α+s;u, and hence δ s Iv, δ0 s 2s α α;u Cn, s, Ω, α v C α R n + [v] α α+s;u Cn, s, Ω, α, ε 0 v C α R n + ε 0 v α α+2s;u for each ɛ 0 > 0. The last inequality is by standard interpolation. Now, using Lemma 2.10 we deduce v α α+2s;u C v C α R n + s v 2s α α;u C v C α R n + δ s g 2s α α;u + δ s v s δ0 s 2s α α;u + Iv, δ0 s s α α;u Cs, Ω, α, ε 0 g L Ω + g s α α;ω + Cε 0 v α α+2s;u, and choosing ε 0 small enough we obtain v α α+2s;u C g L Ω + g s α α;ω. Furthermore, letting U Ω ρ0 we obtain that the same estimate holds with U replaced by Ω ρ0. Finally, in Ω \ Ω ρ0 we have that u is C α+2s and δ s is uniformly positive and C 0,1. Thus, we have u/δ s C γ Ω \ Ω ρ0, where γ = min{1, α + 2s}, and the theorem follows. Next we give the Proof of Corollary 1.6. a It follows from Proposition 1.1 that u C s R n. The interior estimate follow by applying repeatedly Proposition 1.4. b It follows from Theorem 1.2 that u/δ s Ω C α Ω. The interior estimate follows from Theorem 1.5. The following two lemmas are closely related to Lemma 4.2 and are needed in [20] and in Remark 2.11 of this paper. Lemma 4.3. Let U be an open domain and α and β be such that α s < β and β s is not an integer. Let k be an integer such that β = k + β with β 0, 1]. Then, 4.6 [ s/2 w] s α β s;u C w C α R n + w α,