We study the asymptotic behaviour of the least energy solutions to the following class of nonlocal Neumann problems:

$$ \begin{align*} \begin{cases} { d(-\Delta)^{s}u+ u= \vert u\vert^{p-1}u } & \text{in } \Omega, \\ {u>0} & \text{in } \Omega, \\ { \mathcal{N}_{s}u=0 } & \text{in } \mathbb{R}^{n}\setminus \overline{\Omega}, \end{cases} \end{align*} $$

where $\Omega \subset \mathbb {R}^{n}$ is a bounded domain of class $C^{1,1}$, $1<p<({n+s})/({n-s}),\,n>\max \{1, 2s \}, 0<s<1, d>0$ and $\mathcal {N}_{s}u$ is the nonlocal Neumann derivative. We show that for small $d,$ the least energy solutions $u_d$ of the above problem achieve an $L^{\infty }$-bound independent of $d.$ Using this together with suitable $L^{r}$-estimates on $u_d,$ we show that the least energy solution $u_d$ achieves a maximum on the boundary of $\Omega $ for d sufficiently small.

In this paper, we study the global regularity for regular Monge-Ampère type equations associated with semilinear Neumann boundary conditions. By establishing a priori estimates for second order derivatives, the classical solvability of the Neumann boundary value problem is proved under natural conditions. The techniques build upon the delicate and intricate treatment of the standard Monge-Ampère case by Lions, Trudinger, and Urbas in 1986 and the recent barrier constructions and second derivative bounds by Jiang, Trudinger, and Yang for the Dirichlet problem. We also consider more general oblique boundary value problems in the strictly regular case.