Dear all,
I would like to prove the exponential decay of the *derivatives* of a solution to the following equation in $\mathbb{R}^N$:
$$
\sqrt{-\Delta+m^2} u +u= f(u),
$$
where I can assume that $m \neq 0$, that $f$ is smooth and that the solution $u \in H^{1/2}(\mathbb{R}^N)$ is at least Hölder-continuous. I also know that $u(x) = O(e^{-C|x|})$ at infinity. If the equation were local (like $-\Delta u + u = f(u)$), the usual approach would consist in using interior Schauder estimates for $\nabla u$ (or some Harnack-like estimate); I did not find any precise reference for Schauder theory of non-local equations, except in the case $m=0$ (the fractional laplace equation).

Any help is welcome.