The Neumann problem for the fractional Laplacian: regularity up to the boundary

We study the regularity up to the boundary of solutions to the Neumann problem for the fractional Laplacian. We prove that if u is a weak solution of . Å/su D f in and Nsu D 0 in c, then u is C˛ up to the boundary for some ˛>0. Moreover, in case s > 1 2 , we show that u 2 C2s 1C˛. /. To prove these results we need, among other things, a delicate Moser iteration on the boundary with some logarithmic corrections. Our methods allow us to treat as well the Neumann problem for the regional fractional Laplacian, and we establish the same boundary regularity result. Prior to our results, the interior regularity for these Neumann problems was well understood, but near the boundary even the continuity of solutions was open.