Fractional semilinear Neumann problems arising from a fractional Keller-Segel model

We consider the following fractional semilinear Neumann problem on a smooth bounded domain Omega subset of R-n, n >= 2, {(-epsilon Delta)(1/2)u + u = u(p), in Omega, partial derivative(nu)u = 0, on partial derivative Omega, u > 0, in Omega, where epsilon > 0 and 1 < p < (n + )/(n - 1). This is the fractional version of the semilinear Neumann problem studied by Lin-Ni-Takagi in the late 1980's. The problem arises by considering steady states of the Keller-Segel model with nonlocal chemical concentration diffusion. Using the semigroup language for the extension method and variational techniques, we prove existence of nonconstant smooth solutions for small epsilon, which are obtained by minimizing a suitable energy functional. In the case of large epsilon we obtain nonexistence of nonconstant solutions. It is also shown that as epsilon -> 0 the solutions u(epsilon) tend to zero in measure on Omega, while they form spikes in (Omega) over bar. The regularity estimates of the fractional Neumann Laplacian that we develop here are essential for the analysis. The latter results are of independent interest.