Hostname: page-component-586b7cd67f-l7hp2 Total loading time: 0 Render date: 2024-11-21T22:39:29.164Z Has data issue: false hasContentIssue false

On large solutions for fractional Hamilton–Jacobi equations

Published online by Cambridge University Press:  11 July 2023

Gonzalo Dávila
Affiliation:
Departamento de Matemática, Universidad Técnica Federico Santa María, Casilla: v-110, Avda. España 1680, Valparaíso, Chile ([email protected]; [email protected])
Alexander Quaas
Affiliation:
Departamento de Matemática, Universidad Técnica Federico Santa María, Casilla: v-110, Avda. España 1680, Valparaíso, Chile ([email protected]; [email protected])
Erwin Topp
Affiliation:
Departamento de Matemática y C.C., Universidad de Santiago de Chile, Casilla 307 Santiago, 454003, Chile Instituto de Matemática, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ 21941-909, Brazil ([email protected])
Rights & Permissions [Opens in a new window]

Abstract

We study the existence of large solutions for nonlocal Dirichlet problems posed on a bounded, smooth domain, associated with fully nonlinear elliptic equations of order $2\,s$, with $s\in (1/2,\,1)$, and a coercive gradient term with subcritical power $0< p<2\,s$. Due to the nonlocal nature of the diffusion, new blow-up phenomena arise within the range $0< p<2\,s$, involving a continuum family of solutions and/or solutions blowing-up to $-\infty$ on the boundary. This is in striking difference with the local case studied by Lasry–Lions for the subquadratic case $1< p<2$.

Type
Research Article
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press on behalf of The Royal Society of Edinburgh

1. Introduction

In this paper, we study the boundary blow-up phenomena for solutions of fractional problems with coercive gradient with the form

(1.1)\begin{equation} (-\Delta)^s u + |Du|^p + \lambda u =f \quad \mbox{in} \ \Omega, \end{equation}

where $s \in (1/2,\, 1)$, $0 < p < 2\,s,\, \Omega \subset \mathbb {R}^N$ is a bounded, open set with smooth boundary, $f \in C(\Omega )$ and $\lambda \in \mathbb {R}$. Here $(-\Delta )^{s}$ denotes the fractional Laplacian of order $2\,s$, defined for smooth functions $u:\mathbb {R}^N\to \mathbb {R}$ as

(1.2)\begin{equation} (-\Delta)^s u(x) = {C_{N, s} \mathrm{P.V.} \int_{\mathbb{R}^N} \frac{u(x) - u(z)}{|x - z|^{N + 2s}}{\rm d}z}, \end{equation}

whenever the integral converges. Here, the principal value is understood in the Cauchy principal value sense. The constant $C_{N, s} > 0$ is a normalizing constant so that $(-\Delta )^{s}u\to -\Delta u$ as $s\to 1$ (in an adequate functional framework, see e.g. [Reference Di Neza, Palatucci and Valdinoci12]).

The question we address concerns the existence of solutions to (1.1) that become unbounded near the boundary. We prove the existence of multiple solutions to this problem under certain assumptions on the data.

Our method relies on the construction of solutions by Perron's method, as it is presented by Lasry and Lions in [Reference Lasry and Lions13], that is

(1.3)\begin{equation} -\Delta u + |Du|^p + \lambda u = f \quad \mbox{in} \ \Omega, \end{equation}

where $\Delta$ denotes the usual Laplacian, $1 < p \leq 2$, $\lambda > 0$. In this local setting, the Dirichlet problem with blow-up boundary condition associated with (1.3) is complemented by the expression

(1.4)\begin{equation} \lim_{x \in \Omega, \ x \to \partial \Omega} u(x) ={+}\infty. \end{equation}

It is proven in [Reference Lasry and Lions13] the existence of a unique large solution for the Dirichlet problem (1.3)(1.4) for $f \in L^\infty _{loc}(\Omega )$ with an appropriate growth on the boundary. Condition $p \leq 2$ is typically referred as the subcritical case.

Due to the nonlocal nature of the fractional diffusion (1.2), the Dirichlet problem associated with (1.1) requires that we impose a condition on $\Omega ^c$. We consider here the Dirichlet condition $u = \varphi$ in $\Omega ^c$, where $\varphi : \Omega ^c \to \mathbb {R}$ is a given function that satisfies the integrability condition $\varphi \in L^1_w(\Omega ^c)$, where for measurable set $E \subset \mathbb {R}^N$ we denote

\[ L^1_\omega(E) := \left\{ u \in L^1_{loc}(\mathbb{R}^N) : \int_{E} |u| \omega <{+}\infty\right\}, \quad \omega(y) := \frac{1}{(1 + |y|)^{N + 2s}}. \]

Summarizing, the first model for Dirichlet problem we consider here takes the form

(1.5)\begin{align} \left\{\begin{aligned} (-\Delta)^s u+|D u|^p + \lambda u & =f & & \mbox{in} \ \Omega,\\ u & =\varphi & & \mbox{in} \ \Omega^c,\\ \lim\limits_{x\in\Omega, \ x\to\partial\Omega}u & ={+}\infty. & & \end{aligned}\right. \end{align}

The study of large solutions has a long history starting with the work of Keller and Osserman (see [Reference Keller14] and [Reference Osserman15]), where they studied conditions on nonlinearities $f$ in order to find unbounded solutions to

\[ \Delta u=f(u) \quad \text{in }\Omega. \]

For recent developments and extensions see [Reference Alarcón, García-Melián and Quaas2], where a thorough list of references can be found. Existence of blow-up solutions for nonlocal Dirichlet problems has been studied in [Reference Abatangelo1, Reference Chen and Lv7Reference Chen and Veron10]. Of particular interest is the work of Abatangelo [Reference Abatangelo1], where the existence of an intriguing variety of blow-up solutions, for an ample class of reaction–diffusion problems with nonhomogeneous exterior data, is addressed. Roughly speaking, he constructs a fractional harmonic function that blows up near the boundary with a rate given by $d^{s-1}$ (here $d= d(x)$ denotes the distance function $\mathrm {dist}(x,\, \partial \Omega )$ for $x \in \Omega$). These results rely on Green functions and integral formulas for the fractional Laplacian. Using similar tools, existence of weak solutions for quasilinear equations with measure ingredients was already treated by Chen and Veron in [Reference Chen and Veron10]. We mention here that all the mentioned results deal with linear diffusion.

Here we understand the condition $p < 2\,s$ as strictly subcritical in the sense that the growth of the gradient is strictly less than the order of the diffusion. Our result shows a multiplicity phenomenon which is in high contrast with the second-order setting, and it is a consequence of the nonlocal nature of the problem. For instance, if we look on the existence result in [Reference Lasry and Lions13], the authors construct a Perron's solution through blow-up barriers. These are suitable powers of the distance function $d^\beta$. By a natural scaling property of the equation, the exponent $\beta = \frac {p - 2}{p - 1}$ ensures a good approximation for the problem. A logarithmic profile is found in the critical case $p=2$.

We follow the same program here to construct solutions. For introductory purposes, we concentrate on the case $u$ satisfies the homogeneous exterior condition $u = 0$ in $\Omega ^c$. Extending $d$ as zero outside $\Omega$ and for $p \in (1,\,2\,s)$, we have the corresponding exponent

(1.6)\begin{equation} \beta = \frac{2\,s - p}{1 - p}, \end{equation}

as the one making $d^\gamma$ a good ansatz for the problem (1.5). Since we require $d^\gamma$ to be integrable, this introduces new critical exponents of $p$, depending on the diffusive parameter $s$. By the method used to find it, we refer to this solution as a scale solution. Nevertheless, the existence of blow-up fractional harmonic functions involves the existence of a continuum of solutions which are not present in the local framework. Moreover, the nonlocal phenomenon also permits the existence of blow-up solutions to $-\infty$ in certain regimes of $p < 2\,s$. We will come back to more specific aspects of the problem later.

Next, we describe the general class of operators we will consider. For $s \in (1/2,\, 1)$ and constants $0 < \gamma \leq \Gamma < +\infty$, we consider the class $\mathcal {K}$ of measurable kernels $K : \mathbb {R}^N \setminus \{ 0\} \to \mathbb {R}$ such that $K(y) = K(-y)$ for all $y$ and satisfying the ellipticity condition

(1.7)\begin{equation} \gamma |y|^{-(N + 2s)} \leq K(y) \leq \Gamma |y|^{-(N + 2s)}, \quad y \neq 0. \end{equation}

For each $K \in \mathcal {K}$, we consider the linear operator

(1.8)\begin{equation} L_K u(x) := \mathrm{PV} \int_{\mathbb{R}^n} [u(x + y)-u(x)] K(y) {\rm d}y, \end{equation}

which is well defined for measurable $u: \mathbb {R}^N \to \mathbb {R}$ satisfying adequate regularity assumptions on $x$ and weighted integrability condition at infinity; typically, $u \in C^{1,1}$ in a neighbourhood of $x$ and $u \in L^1_\omega (\mathbb {R}^N)$. Thus, for a two-parameter family of kernels $\{ K_{ij}\}_{i \in I, j \in J} \subset \mathcal {K}$, and denoting $L_{ij} := L_{K_{ij}}$, we write

(1.9)\begin{equation} \mathcal{I} u(x):= \inf\limits_{i\in I} \sup\limits_{j \in J} L_{ij}u (x). \end{equation}

Associated with $\mathcal {I}$ of this type, we consider

(1.10)\begin{equation} \lambda_0(\mathcal{I}) = \inf_{x \in \Omega, i \in I, j \in J} \int_{\Omega^c} K_{ij}(x - y) {\rm d}y. \end{equation}

Notice that $0 < \lambda _0(\mathcal {I}) < +\infty$.

We focus on the fully homogeneous class of kernels with the form

(1.11)\begin{equation} K(z) = \frac{a(z/|z|)}{|z|^{N + 2s}}, \end{equation}

for some nonnegative, measurable function $a: S^{N - 1} \to \mathbb {R}$. In this setting, condition (1.7) turns out to be $\gamma \leq a \leq \Gamma$.

In order to describe our existence results, we require to introduce two exponents related to the gradient nonlinearity, depending on the order $s$. We write $p_i = p_i(s)$ for $i = 1,\, 2$ as

(1.12)\begin{equation} p_1= s + \frac{1}{2} \quad\mbox{and} \quad p_2=\frac{s+1}{2-s}, \end{equation}

A third exponent $p_0=\frac {2s}{2-s}$ also emerges in our analysis, but this plays a less important role. We notice that for $s \in (1/2,\, 1)$, we have $p_1 > 1$ and

\[ p_1(s)< p_2(s)< 2\,s, \quad p_2(1^-)=2, \quad p_1(1/2^+)=p_2(1/2^+)=1, \]

Thus, our main result is the existence of boundary blow-up solutions for the Dirichlet problem

(P)\begin{equation}\left\{ \begin{array}{@{}rl} {-\mathcal{I}(u) + |Du|^p + \lambda u = f} \quad & {{\rm in}\;\Omega ,} \\ {u = 0} \quad & {{\rm in}\;\Omega ^c,} \end{array}\right.\end{equation}

as it is described by the following.

Theorem 1.1 Let $s \in (1/2,\,1)$, $0 < p < 2\,s$, $\Omega \subset \mathbb {R}^N$ be a bounded domain with $C^2$ boundary, $f\in L^\infty (\Omega ) \cap C(\Omega )$. When $1 < p < 2\,s$, let $\beta$ as in (1.6).

Let $\mathcal {K}$ be a family of symmetric kernels satisfying (1.7) and (1.11), $\mathcal {I}$ a nonlinear operator with form (1.9), and $\lambda > -\lambda _0(\mathcal {I})$. Let $p_i$ be defined as in (1.12), $i=0,\,1,\,2$.

Then, we have the following existence results:

1. One parameter family of solutions (close to $s$-harmonic): If $0 < p < p_2$, there exist $\sigma > 0$ and a family of solutions $\{ u_t \}_{t \in \mathbb {R}, t \neq 0} \subset C^{\sigma }_{loc}(\Omega )$ to (P), such that for each $t$ we have

\[ d^{1 - s}(x) u_t(x) - t = O(d^{\theta}(x)), \]

for some $\theta > 0$ depending on $p$. In particular, if $t_1 < t_2$, then

\[ u_{t_1} < u_{t_2} \quad \mbox{in} \ \Omega. \]

Moreover, if $p$ additionally satisfies $p < p_0$, then we can take $\theta >1 - s$ .

2. Positive scale solution: If $p_1 < p < p_2$, then there exist $\sigma > 0$ and a constant $T > 0$ and a function $u \in C^{\sigma }(\Omega )$ solving (P) such that

\[ d(x)^{-\beta}u(x) - T = O(d(x)^\theta), \]

for some $\theta > 0$.

3. Negative scale solution: For $p_2 < p < 2\,s$, then there exist $\sigma > 0$, $T > 0$ and a solution $u \in C^{\sigma }(\Omega )$ of (P) such that

\[ d^{-\beta}(x)u(x) + T = O(d(x)^\theta), \]

for some $\theta > 0$.

For the reader who is familiarized with nonlocal problems, condition $\lambda > -\lambda _0(\mathcal {I})$ is a condition ensuring comparison principle, see for instance [Reference Topp19]. Since we employ an approximation procedure to construct the solutions (c.f. proposition 2.2), comparison principle plays an important role.

As we previously mentioned, we construct barriers that blow-up (to $+\infty$) as powers of the distance function. This allows us to construct scale solutions behaving like $d^\beta$ with $\beta$ as in (1.6). This imposes the condition $p < p_2$ in order for the power function to be integrable. In the regime $p > p_2$, the same argument allows us to get a negative blow-up solution.

Other types of large solutions, not present in the local case, emerge here. For instance in case 1, the family $\{ u_t \}$ is constructed ‘around’ a blow-up, fractional harmonic function. Our key technical result (c.f. proposition 3.1) indicates that the function $d^{s-1}$ is close to being harmonic for $\mathcal {I}$. In case 3, we construct barriers perturbing $d^{s-1}$ with other power-type functions of lower order, leading to the existence of solutions that diverge to $-\infty$.

There are plenty of open questions after this work that we believe deserve to be investigated. For instance, we could not cover the natural critical exponents as $p=p_2$ and $p = 2\,s$. None of the solutions found in theorem 1.1 converge to a blow-up (or blow-down) solution of (P) when $p \to p_2$. In fact, they surprisingly converge to the unique bounded viscosity solution to (P), see remark 4.1. On the other hand, the critical case $p = 2\,s$ resembles the case $p=2$ in (1.3), for which a logarithmic blow-up profile for the solution is obtained, see [Reference Lasry and Lions13]. In our context, the construction of barriers with a log profile is hard to handle at a technical level, and therefore we did not pursue it in this work.

A question about applications that emerges here has to do with the connection of problem (P) with stochastic optimal control problems. In the local setting presented in [Reference Lasry and Lions13], the solution to (1.3) turns out to be the value function of an infinite horizon stochastic optimal control problem with a cost involving $f$ and a ‘feedback’ term depending on $p$. The admissible drifts are those preventing the trajectories of the stochastic process to exit the domain, leading to the blow-up of the associated value function. Here, for each $K$ like in (1.11), its associated linear operator $L_K$ is the infinitesimal generator of a $s$-stable pure jump Lévy process, see for instance [Reference Bertoin5, Reference Ros-Oton and Serra17, Reference Sato18]. The connection of fractional Hamilton–Jacobi equations and stochastic optimal control problems with jumps has been shown to hold in some cases, for instance in the unrestricted state case $(\Omega = \mathbb {R}^N)$. The fully nonlinear structure of (P) involves stochastic differential equations with controlled random parts, see for instance [Reference Pham16, Reference Sato18]. The paper is organized as follows: In § 2, we provide the notion of solution we use here and a general Perron's method. In § 3, we provide precise estimates of the nonlocal operator applied to powers of the distance function. In § 4, we provide the proof of theorem 1.1, which is accomplished by constructing sub and supersolutions based on the results of § 3. Finally, in § 5, we provide some extensions that include unbounded right-hand side and nonhomogeneous exterior Dirichlet data.

2. Perron's solutions

Given $\varphi \in L^1_\omega (\Omega ^c)$, $\lambda > -\lambda _0(\mathcal {I})$ and $f \in C(\Omega )$, we are interested in viscosity solutions to the problem

(2.1)\begin{align} \left\{\begin{aligned} -\mathcal{I} u+|D u|^p + \lambda u & =f & & \mbox{in} \ \Omega,\\ u & =\varphi & & \mbox{in} \ \Omega^c,\\ \lim\limits_{x\in\Omega, \ x\to\partial\Omega}u & ={+}\infty, & & \end{aligned}\right. \end{align}

and its blow-up version to minus infinity, that is, $\lim \limits _{x\in \Omega, \ x\to \partial \Omega }u = -\infty$.

We start with the notion of viscosity solution, see [Reference Barles, Chasseigne and Imbert3].

Definition 2.1 A function $u \in L^1_\omega (\mathbb {R}^N)$, upper semicontinuous in $\Omega$, is a viscosity subsolution to the Dirichlet problem

(2.2)\begin{align} \left\{\begin{aligned} -\mathcal{I} u+|D u|^p + \lambda u & =f & & \mbox{in} \ \Omega,\\ u & =\varphi & & \mbox{in} \ \Omega^c, \end{aligned}\right. \end{align}

if $u \leq \varphi$ in $\Omega ^c$, and, for every $x_0 \in \Omega$ and every function $\phi \in L^1_\omega (\mathbb {R}^N) \cap C^2(\Omega )$ such that $u(x_0) = \phi (x_0)$, $u \leq \phi$ in $B_\delta (x_0)$ for some $\delta > 0$, we have the inequality

\[ -\mathcal{I} u^\phi_{\delta, x_0}(x_0) + |D \phi(x_0)|^p + \lambda u(x_0) \leq f(x_0), \]

where $u^\phi _{\delta, x_0}: \mathbb {R}^N \to \mathbb {R}$ is the function defined as $u^\phi _{\delta, x_0}(x) = \phi (x)$ in $B_\delta (x_0)$, $u^\phi _{\delta, x_0}(x) = u(x)$ in $B_\delta ^c(x_0)$.

A function $u \in L^1_\omega (\mathbb {R}^N)$, lower semicontinuous in $\Omega$, is a viscosity supersolution to the Dirichlet problem (2.2) if $u \geq \varphi$ in $\Omega ^c$ and for every $x_0 \in \Omega$ and every function $\phi \in L^1_\omega (\mathbb {R}^N) \cap C^2(\Omega )$ such that $u(x_0) = \phi (x_0)$, $u \geq \phi$ in $B_\delta (x_0)$ for some $\delta > 0$, we have the inequality

\[ -\mathcal{I} u^\phi_{\delta, x_0}(x_0) + |D \phi(x_0)|^p + \lambda v(x_0) \geq f(x_0), \]

where $u^\phi _{\delta, x_0}$ is defined as before.

A function $u \in L^1_\omega (\mathbb {R}^N) \cap C(\Omega )$ is a solution to (2.2) if $u = \varphi$ in $\Omega ^c$ and is simultaneously a viscosity sub and supersolution to the problem.

Finally, we say that $u$ is a strict subsolution (resp. supersolution) to (2.2) if there exists $\epsilon > 0$ such that $u$ satisfies the viscosity inequality with $f(x_0)-\epsilon$ (resp. $f(x_0) + \epsilon$) instead of $f(x_0)$, for all $x_0 \in \Omega$.

Existence and uniqueness for solutions $u \in C(\mathbb {R}^N)$ can be found in [Reference Barles, Chasseigne and Imbert3], in particular it attains the boundary data imposed by $\varphi$. However, definition 2.1 allows the possibility to have solutions which are unbounded in $\Omega$. We use an approximation procedure based on Perron's method.

Proposition 2.2 Let $\Omega$ be a bounded, smooth domain, $f \in C(\Omega ),$ $0< p\leq 2\,s$ and $\lambda >-\lambda _0(\mathcal {I})$. Suppose there exist a supersolution $\bar U$ and a subsolution ${\underline U}$ of (2.1) with $\bar U,\, {\underline U}\in C(\Omega )\cap L^1_\omega (\mathbb {R}^N),$ and ${\underline U} = \bar U = \varphi$ in $\Omega ^c,$ and such that one of them is strict. Furthermore, assume that

(2.3)\begin{equation} \bar U\geq{\underline U} \quad \mbox{in} \ \mathbb{R}^N, \quad \lim\limits_{x\in\Omega,\ x\to\partial\Omega}{\underline U}={+}\infty. \end{equation}

Then there exists a solution $u \in C^\alpha (\Omega ) \cap L^1_w(\mathbb {R}^N)$ of (2.1) satisfying ${\underline U}\leq u\leq \bar U$.

An analogous result can be stated for sub and supersolutions ${\underline U} \leq \bar U$ with $\bar U(x) \to -\infty$ as $x \to \partial \Omega$.

Proof. We assume ${\underline U}$ is a strict subsolution, the other case follows the same lines. Let $\Omega _n = \{x\in \Omega : \text {dist}(x,\,\partial \Omega )>1/n\}$. For $k,\, n \in \mathbb {N}$, let $w_{n, k}$ be a continuous function in $\bar \Omega \setminus \Omega _{n + k + 1}$ such that $w_{n,k} = \varphi$ on $\partial \Omega$, $w_{n,k} = {\underline U}$ on $\partial \Omega _{n + k + 1}$ (say, the harmonic function in $\Omega \setminus \bar \Omega _{n + k+ 1}$ satisfying the mentioned boundary conditions).

Now, let $U_{n,k} : \Omega _{n}^c \to \mathbb {R}$ given by

\[ U_{n,k}(x) = \left \{ \begin{array}{@{}ll} {\underline U}(x) \quad & \mbox{if} \ x \in \Omega_{n + k + 1} \setminus \Omega_n, \\ \min\{ w_{n,k}(x), {\underline U}(x) \} \quad & \mbox{if} \ x \in \Omega \setminus \Omega_{n + k + 1}, \\ \varphi(x) & \mbox{if} \ x \in \Omega^c. \end{array} \right . \]

Notice that $U_{n,k}$ is continuous, and $\min _{\partial \Omega } \{ \varphi \} \leq U_{n, k} \leq {\underline U}$ in $\Omega \setminus \Omega _{n + k + 1}$, from which, by dominated convergence theorem, we have that

\[ \int_{\Omega \setminus \Omega_{n + k + 1}} |U_{n,k}(y) - {\underline U}(y) |K(x - y) {\rm d}y \to 0, \]

as $k \to \infty$, uniformly in $x \in \Omega _n,\, K \in \mathcal {K}$, for $n$ fixed. Then, since ${\underline U}$ is a strict subsolution, the above estimate implies that for each $n$, there exists $k(n)$ such that, for each $k \geq k(n)$, the function $U_{n,k}$ is a viscosity subsolution to the Dirichlet problem

(2.4)\begin{align} \left \{ \begin{aligned} -\mathcal{I} u+|Du|^p+\lambda u & =f - \dfrac{1}{n} \quad & & \text{in }\Omega_n, \\ u & = U_{n,k} \quad & & \text{in }\Omega_n^c. \end{aligned} \right . \end{align}

By a similar argument, using that $\bar U$ is a viscosity supersolution for the problem in $\Omega$, we can construct $\bar U_{n,k} \in C(\mathbb {R}^N)$ a supersolution to (2.4) with $\bar U_{n,k} = \bar U$ in $\Omega _{n + k + 1}$, $U_{n,k} \leq \bar U_{n,k}$ in $\Omega$, and such that $\bar U_{n,k} = \varphi$ in $\Omega ^c$, for all $k \geq k(n)$ (relabelling $k(n)$ if necessary). Thus, by Theorem 1 in [Reference Barles, Chasseigne and Imbert3], there exists a unique viscosity solution $u_{n,k} \in C(\mathbb {R}^N)$ for (2.4). Moreover, it satisfies ${\underline U} \leq u_{n,k} \leq \bar U$ in $\Omega _{n + k + 1}$, for all $n$ and $k \geq k(n)$, and by construction we have $u_{n,k} \in L^1_\omega (\mathbb {R}^N)$ uniformly in $n$ and $k \geq k(n)$. Comparison principles are available by the assumption $\lambda > -\lambda _0 (\mathcal {I})$ for all $n,\, k$ large.

Thus, the family $\{ u_{n,k(n)} \}_{n \in \mathbb {N}}$ has uniform interior $C^{\alpha }$ estimates by the results of [Reference Guy, Emmanuel and Cyril4]. Using this and that the family is uniformly bounded in compact sets of $\Omega$, we can use stability results of viscosity solutions to conclude the result, taking $n \to +\infty$.

Remark 2.3 We notice that the above result holds if we assume that $\bar U$ (resp. $\bar U$) is a viscosity subsolution (resp. supersolution) to (2.1) which is strict in each compact subset of $\Omega$.

3. Technical lemmas

We use the notation $d: \mathbb {R}^N \to \mathbb {R}$ such that $d(x) = \mathrm {dist}(x,\, \partial \Omega )$ for $x \in \Omega$, and $d(x) = 0$ for $x \in \Omega ^c$. Since the domain is smooth, we have the existence $\delta _0 > 0$ such that $d$ is a $C^2$ function on the set $\Omega _{\delta } = \{ x \in \Omega : d(x) < \delta \}$ for all $\delta < \delta _0$. Given $\tau \in (-1,\, 2\,s)$, we denote $d^\tau : \mathbb {R}^N \to \mathbb {R}$ such that $d^\tau (x) = (d(x))^\tau$ for $x \in \Omega$ and zero in $\Omega ^c$, with the convention and $d^0 = \chi _\Omega$.

For a function $u: \mathbb {R}^N \to \mathbb {R}$ measurable, $A \subset \mathbb {R}^N$ measurable set and $x \in \mathbb {R}^N$, we denote

\[ L_K[A]u(x) = \mathrm{P.V.} \int_A[u(x + z) - u(x)]K(z){\rm d}z. \]

Here, we only assume (1.7). For $K$ in this class and $\rho > 0$, we introduce the notation

(3.1)\begin{equation} K^\rho(z) = \rho^{N + 2s} K(\rho z), \quad z \neq 0. \end{equation}

Notice that $K^\rho$ satisfies (1.7) with the same ellipticity constants as $K$, and that if $K$ satisfies (1.11), then $K^\rho = K$.

Finally, we consider the extremal Pucci operators associated with the family of kernels $\mathcal {K}$ as

\[ \mathcal{M}^+ u(x) = \sup_{K \in \mathcal{K}} L_K u(x), \quad \mathcal{M}^- u(x) = \inf_{K \in \mathcal{K}} L_K u(x). \]

The main result of this section is the following

Proposition 3.1 Let $\Omega \subset \mathbb {R}^N$ be a bounded domain with $C^2$ boundary, $s\in (0,\,1),$ and let $\mathcal {K}$ be a family of kernels satisfying (1.7).

Then, for each $\tau \in (-1,\, 2s),$ there exists $\delta > 0$ such that

\[ \mathcal{I} d^{\tau}(x) = d^{\tau - 2s}(x) (c(d(x), \tau) + O(d(x)^s)), \quad x \in \Omega_\delta, \]

where

\[ c(d(x), \tau) = \inf_{i \in I} \sup_{j \in J} \mathrm{P.V.} \int_{\mathbb{R}^N} [(1 + z_N)_+^\tau - 1]K_{ij}^{d(x)}(z){\rm d}z. \]

Finally, if we additionally assume (1.11) , then $c( \mathcal {I},\, d(x),\, \tau )= c(\tau )$ and this constant satisfies $c(-1^+) =+\infty,$ $c(2s^-) = +\infty,$ $c(s - 1) = c(s) = 0,$ $c(\tau ) > 0$ if $\tau \in (-1, s-1) \cup (s,\, 2\,s)$ and $c(\tau ) < 0$ for $\tau \in (s-1,\, s)$.

Before we continue with the proof of the proposition, we will introduce some notation. When (1.11) holds, then we denote

\[ c(\mathcal{I}, \tau)=c( \mathcal{I}, d(x), \tau) \]

and we will omit the dependence on $\mathcal {I}$ whenever the context is clear. In particular, if we define $\tilde{\mathcal {I}}$ as the operator $- \mathcal {I} (- \cdot )$ then it satisfies (1.11) so we will denote $\tilde c(\tau )=c( \tilde{\mathcal {I}},\, d(x),\, \tau )$. Moreover, since $\mathcal {M}^\pm$ also satisfies (1.11) we will denote $c^\pm (\tau )=c( \mathcal {M}^\pm,\, d(x),\, \tau )$.

We recall that for each $x \in \partial \Omega$, there exists an open set $\mathcal {U} \subset \mathbb {R}^{N - 1}$ containing the origin, $r > 0$ and a $C^2$ function $\psi _x: \mathcal {U} \to \mathbb {R}$ such that $\partial \Omega \cap B_{\mathbb {R}^N}(x,\,r) \subset \mathcal {R}_x \{ x + (z',\, \psi _z(z')) : z' \in \mathcal {U} \}$, for some rotation matrix $\mathcal {R}_x$. By compactness and regularity of $\partial \Omega$, we have a finite number of charts covering $\partial \Omega$ with uniform $C^2$ bounds.

Now, let $x \in \Omega$, and denote $\rho = d(x)$. After rotation, we assume $x = \rho e_N$, the projection of $x$ to $\partial \Omega$ is the origin, and therefore, that the local chart $\psi$ corresponding to this point satisfies $\psi (0) = 0,\, D\psi (0) = 0$. Thus, we have the existence of $C_\Omega > 0$ such that

(3.2)\begin{equation} |\psi(z')| \leq C_\Omega |z'|^2. \end{equation}

The key technical step to prove proposition 3.1 is the following

Lemma 3.2 Let $K \in \mathcal {K}$ satisfying (1.7), and $\tau \in (-1,\, 2\,s)$. Let $\eta \in (0,\,1)$, $Q_\eta = B_\eta ' \times (-\eta,\, \eta ) \subset \mathbb {R}^N$. Let $x = (0',\, \rho ) \in \Omega$ such that the projection of $x$ to $\partial \Omega$ is the origin, and that $x \in Q_{\eta }$. Denote

\[ I := \lim_{\epsilon \to 0} \int_{Q_\eta \setminus B_\epsilon} [(\rho + y_N - \psi(y'))_+^\tau - \rho^\tau] K(y) {\rm d}y. \]

Then, the limit exists. For each $\eta > 0$ small enough, and all $\rho > 0$ small enough in terms of $\eta,$ we have the expansion

(3.3)\begin{equation} I = \rho^{\tau - 2s} (c_K(\rho, \tau) + O(\rho^s) + O(\rho^{\tau + 1})), \end{equation}

where

\[ c_K(\rho, \tau) = \mathrm{P.V.} \int_{\mathbb{R}^N} [(1 + z_N)_+^\tau - 1] K^\rho(z) {\rm d}z, \]

and the $O$ terms depend only on $N,\,s,\, \Omega,\, 1 + \tau,$ $\eta$ and the ellipticity constants.

We use this estimate to prove our main result of this section.

Proof of proposition 3.1. Proof of proposition 3.1

For an arbitrary linear operator $L = L_K$ in the family, and for $\eta$ as in lemma 3.2, we write

(3.4)\begin{equation} L d^\tau (x) = L[Q_\eta] d^\tau(x) + L[Q_\eta^c] d^\tau(x). \end{equation}

It is easy to see that if $\rho < \eta /4$, we have

\[ |L[Q_\eta^c] d^\tau(x)| \leq \int_{Q_\eta^c} |d^\tau(x + y) - \rho^\tau|K(y){\rm d}y \leq C\Lambda (c_\tau \eta^{-(N + 2s)} + \rho^\tau \eta^{{-}2s}), \]

for some $C > 0$ just depending on $N,\, s,\, \Omega$.

From now on, we concentrate on $L[Q_\eta ] d^\tau (x)$ in (3.4). For each $z = (z',\, z_N) \in \Omega \cap Q_\eta$, we have

\[ d(z) \leq z_N - \psi(z'), \]

from which we directly have

(3.5)\begin{equation} L d^\tau (x) \geq I, \end{equation}

with $I$ as in lemma 3.2.

On the other hand, by the smoothness of the domain, we use Lemma 3.1 in [Reference Chen, Felmer and Quaas9], from which we get the existence of a constant $C_\Omega > 0$ just depending on $\Omega$ and $N$ such that

\[ d(y) \geq (y_N - \psi(y'))(1 - C_\Omega|y'|^2), \]

for $y$ close to the boundary, near $x_0$. Thus, taking $\eta$ small enough in terms of $C_\Omega$, we also have

\[ L d^\tau (x) \leq \lim_{\epsilon \to 0^+} \int_{Q_\eta \setminus B_\epsilon} [(\rho + y_N - \psi(y'))_+^\tau (1 + C_\Omega |y'|^2) - \rho^\tau] K(y) {\rm d}y, \]

and from here, it is easy to see that

(3.6)\begin{equation} L d^\tau (x) \leq I + C_\Omega \int_{Q_\eta} (\rho + y_N - \psi(y'))^\tau_+ |y'|^2 K(y){\rm d}y. \end{equation}

Now, for the second term in the last expression, we can write

\begin{align*} & \int_{Q_\eta} (\rho + y_N - \psi(y'))^\tau_+ |y'|^2 K(y){\rm d}y \\ & \quad\leq \Lambda \int_{B_\eta'} |y'|^{2 - N - 2s} \int_{\psi(y') - \rho}^{\eta} (\rho + y_N - \psi(y'))^\tau {\rm d}y_N {\rm d}y' \\ & \quad\leq \Lambda \frac{1}{1 + \tau} (\rho + \eta/4)^{1 + \tau}\int_{B_\eta'} |y'|^{2 - N - 2s} {\rm d}y' \\ & \quad\leq C \frac{\eta^{1 + \tau}}{1 + \tau}, \end{align*}

where we have used the fact that $\psi (z') \leq \eta /4$ for $|z'| \leq \eta$. Hence, replacing in (3.6) and using (3.5), we conclude that

(3.7)\begin{equation} I \leq L d^\tau (x) \leq I + C C_\Omega \frac{1}{1 + \tau}, \end{equation}

where the constant $C > 0$ depends on $N,\, s,\, \Lambda$

Then, by lemma 3.2, we get

\[ L d^\tau(x) = \rho^{\tau - 2s} (c_K(\rho, \tau) + O(\rho^s) + O(\rho^{1 + \tau}) + O(\rho^{2\,s - \tau})). \]

We get from here that

\[ \mathcal{I} d^\tau(x) = \rho^{\tau - 2s} \Bigg{(} \inf_{i \in I} \sup_{j \in J} c_{K_{ij}}(\rho, \tau) + O(\rho^s) + O(\rho^{1 + \tau}) \Bigg{)}, \]

from which the first result follows.

For the last part of the proposition, since the operator is of form (1.11) we have $K_{ij}^\rho = K_{ij}$, then by Lemma 2.1 in [Reference Ros-Oton and Serra17] we arrive at

(3.8)\begin{equation} c(\tau):=\inf_{i \in I} \sup_{j \in J} c_{K_{ij}}(\tau) = \inf_i \sup_j \{ -\tilde a_{ij} \ (-\Delta)_{\mathbb{R}}^s w_\tau (1) \} \end{equation}

where $(\Delta _{\mathbb {R}})^s$ denotes the fractional Laplacian in dimension one, $w_\gamma (t) = t_+^{\tau }$ for $t \in \mathbb {R}$, and

\[ \tilde a_{ij} := \int_{S^{N - 1}} |\theta_N|^{2s} a_{ij}(\theta) d\sigma(\theta), \]

where $\sigma$ denotes the $N-1$ dimensional Hausdorff measure in the unit sphere. Now if we define $c_1(\tau )= (-\Delta _{\mathbb {R}})^s w_\tau (1)$, the qualitative properties follow since the function $c_1$ is strictly concave in (-1,2 s) by proposition 3.1 of [Reference Chen, Felmer and Quaas9] and $c_1(-1^+) =+\infty$, $c_1(2s^-) = +\infty$. Moreover, $c_1(s)= c_1(s-1) = 0$ by Lemma 6.2 of [Reference Ros-Oton and Serra17].

The rest of the section is devoted to the

Proof of lemma 3.2. Proof of lemma 3.2

We concentrate on the case $\tau < 0$ since it is the most difficult due to the unboundedness of $d^\tau$.

We write

\[ I = I_0 + \tilde I, \]

with

\begin{align*} I_0& = \mathrm{P.V.} \int_{Q_\eta} [(\rho + y_N)_+^\tau - \rho^{\tau}] K(y){\rm d}y, \\ \tilde I & = \int_{Q_\eta} [(\rho + y_N - \psi(y'))_+^\tau - (\rho + y_N)_+^\tau] K(y){\rm d}y, \end{align*}

where the last integral is well defined since, using (3.2), we have

\[ |(\rho + y_N - \psi(y'))_+^\tau - (\rho + y_N)^\tau| \leq C\rho^{\tau - 2} |y|^2, \]

for all $|y|$ small enough in terms of $\rho$. This is enough to compensate the singularity of the kernel $K$ and pass to the limit as $\epsilon \to 0$ using dominated convergence theorem.

1. Estimate for $\tilde I$. Rescaling, we have

\[ \tilde I = \rho^{\tau - 2s} \int_{Q_{\eta/\rho}} [(1 + z_N - \tilde \psi(z'))_+^\tau - (1 + z_N)_+^\tau] K^\rho(z){\rm d}z, \]

where $\tilde \psi (z') := \rho ^{-1} \psi (\rho z')$ and $K^\rho (z) = \rho ^{N + 2s} K(\rho z)$. Notice that by the ellipticity condition, we have

\[ \gamma |z|^{-(N + 2s)} \leq K^\rho(z) \leq \Lambda |z|^{-(N + 2s)}, \quad z \in \mathbb{R}^N \setminus \{ 0 \}, \ \rho > 0. \]

For this, we divide the integral in several parts. Namely, we consider the splitting

(3.9)\begin{equation} d^{2\,s - \tau} \tilde I = I_1 + I_2 + I_3, \end{equation}

where for $i = 1,\,2,\,3$, we denote

\begin{align*} & I_i = \int_{A_i} [(1 + z_N - \tilde \psi(z'))_+^\tau - (1 + z_N)_+^\tau] K^\rho(z){\rm d}z, \ \mbox{with} \\ & A_1 = B_1' \times \left(-\frac{\eta}{\rho}, \frac{\eta}{\rho}\right), \ A_2 = B_{\frac{\eta}{\sqrt{\rho}}}' \setminus B_1' \times \left(-\frac{\eta}{\rho}, \frac{\eta}{\rho}\right), A_3 = Q_{\frac{\eta}{\rho}} \setminus (A_1 \cup A_2). \end{align*}

For $I_1$, we notice that by the assumptions on the chart $\psi$ (c.f. (3.2)), we have $|\tilde \psi (z')| \leq \rho |z'|^2$ if $|z'| \leq 1$. Then, we perform the subdivision

\begin{align*} & I_1 = I_{11} + I_{12} + I_{13}, \\ & A_{11} = B_1' \times \left(-\frac{\eta}{\rho}, -\frac{1}{2}\right), \ A_{12} = B_1' \times \left(-\frac{1}{2}, \frac{1}{2}\right), \ A_{13} = B_1' \times \left(\frac{1}{2}, \frac{\eta}{\rho}\right), \end{align*}

where we have adopted a similar notation as in (3.9).

For $I_{11}$, we make a subdivision with the form

\begin{align*} I_{11} & \leq \int_{B_1'} \int_{\tilde \psi_-(z') - 1}^{{-}1} (1 + z_N - \tilde \psi(z'))^\tau K^\rho(z){\rm d}z \\ & \quad+ \int_{B_1'} \int_{\tilde \psi_+(z') - 1}^{{-}1} [(1 + z_N - \tilde \psi_+(z'))^\tau - (1 + z_N)^\tau]K^\rho(z){\rm d}z \\ & = : I_{111} + I_{112}, \end{align*}

where, for $a \in \mathbb {R}$, we have written $a = a_+ + a_-$.

Using the ellipticity condition and integrating by parts, we can write

\begin{align*} I_{112} &\le \Gamma \int_{B_{{1}}^{\prime}} {\int_{{\tilde{\psi }}_ + ({z}^{\prime})-1}^{-1/2} {\displaystyle{{{(1 + z_N-{\tilde{\psi }}_ + ({z}^{\prime}))^\tau} -{(1 + z_N)^\tau} } \over {|z|^{N + 2s}}}} } {\rm d}z_N{\rm d}{z}^{\prime} \\ & = \displaystyle{\Gamma \over {1 + \tau }}\int_{B_{{1}}^{\prime}} {\left\{ {\displaystyle{{{(1 + z_N-{\tilde{\psi }}_ + ({z}^{\prime}))}^{\tau + 1}-{(1 + z_N)}^{\tau + 1}} \over {|z|^{N + 2s}}}\Bigg{|}_{z_N = {\tilde{\psi }}_ + ({z}^{\prime})-1}^{z_N = -1/2} } \right.} {\rm } \\ & \quad + (N + 2{\mkern 1mu} s)\left. {\int_{{\tilde{\psi }}_ + -1}^{-1/2} {\displaystyle{{{(1 + z_N-{\tilde{\psi }}_ + ({z}^{\prime}))}^{\tau + 1}-{(1 + z_N)}^{\tau + 1}} \over {|z|^{N + 2{\kern 1pt} s + 1}}}} z_N{\rm d}z_N} \right\}{\rm d}{z}^{\prime} \\ & \le \displaystyle{{C\Gamma } \over {1 + \tau }}\int_{B_{{1}}^{\prime}} {\left\{ {\displaystyle{{(\tau + 1)\rho |{z}^{\prime}|^2 + \rho ^{\tau + 1}|{z}^{\prime}|^{2(\tau + 1)}} \over {{(1 + |{z}^{\prime}|^2)}^{(N + 2s)/2}}}} \right.} {\rm } \\ & \quad + \left. {\int_{{\tilde{\psi }}_ + ({z}^{\prime})-1}^{-1/2} {\displaystyle{{{(1 + z_N-{\tilde{\psi }}_ + ({z}^{\prime}))}^{\tau + 1}-{(1 + z_N)}^{\tau + 1}} \over {|z|^{N + 2{\kern 1pt} s + 1}}}} z_N{\rm d}z_N} \right\}{\rm d}{z}^{\prime}{\rm } \\ & \le \displaystyle{{C\Gamma } \over {1 + \tau }}\left\{ {\rho ^{\tau + 1} + \int_{B_{{1}}^{\prime}} {\int_{{\tilde{\psi }}_ + ({z}^{\prime})-1}^{-1/2} {} } } \right. \\ & \qquad \qquad \quad \left. { \times \displaystyle{{{(1 + z_N-{\tilde{\psi }}_ + ({z}^{\prime}))^{\tau + 1}}-{(1 + z_N)^{\tau + 1}}} \over {|z|^{N + 2{\kern 1pt} s + 1}}}z_N{\rm d}z_N{\rm d}{z}^{\prime}} \right\},\end{align*}

for some universal constant $C > 0$. Here, we have used that $|\tilde \psi (z')| \leq C \rho |z'|^2$ for $|z'| \leq 1$. Since $1 + \tau > 0$, we see that

\[I_{112} \le \displaystyle{{C\Gamma } \over {1 + \tau }}{\rm }\left\{ {\rho ^{\tau + 1} + \int_{B_{{1}}^{\prime}} {\displaystyle{{\rho |{z}^{\prime}|^2} \over {{(1 + |{z}^{\prime}|^2)^{{{N + 2{\kern 1pt} s + 1} \over 2}}}}}} \int_{{\tilde{\psi }}_ + -1}^{-1/2} {{(1 + z_N)}^\tau } {\rm d}z_N{\rm d}{z}^{\prime}} \right\}, \]

and from this, we conclude that

\[ I_{112} \leq \frac{C \Gamma}{1 + \tau} \rho^{\tau + 1}. \]

For $I_{111}$, by direct integration and the estimates for $\tilde \psi$, we see that

\[ I_{111} \leq C\Gamma \int_{B_1'} |\tilde \psi(z')|^{1 + \tau} {\rm d}z \leq C \Gamma \rho^{1 + \tau}. \]

For the lower bound, we proceed similarly, noticing that this time we can write

(3.10)\begin{equation} \begin{aligned} I_{11} & \geq{-}\int_{B_1'} \int_{{-}1}^{- 1 + \tilde \psi_+(z')} (1 + z_N)^\tau K^\rho(z){\rm d}z \\ & \quad+ \int_{B_1'} \int_{- 1 + \tilde \psi_+(z')}^{{-}1/2} [(1 + z_N - \tilde \psi_+(z'))^\tau - (1 + z_N)^\tau] K^\rho(z){\rm d}z. \end{aligned} \end{equation}

From here, by direct integration in the first term in the last inequality, and a similar procedure leading to the estimate concerning $I_{112}$ above for the second term, we conclude that

(3.11)\begin{equation} -\frac{C\Gamma}{1 + \tau} \rho^{\tau + 1} \leq I_{11} \leq \frac{C \Gamma}{1 + \tau} \rho^{\tau + 1}, \end{equation}

For $I_{12}$, we use that $|\tilde \psi (z')| \leq C \rho |z'|^2$ to perform a first-order Taylor expansion to get

(3.12)\begin{equation} -C \Gamma \rho \leq I_{12} \leq C \Gamma \int_{B_1'} \int_{{-}1/2}^{1/2} \frac{\rho |z'|^2}{|z|^{N + 2s}} {\rm d}z_N {\rm d}z' \leq C \Gamma \rho. \end{equation}

For $I_{13}$, we perform a Taylor expansion again, from which we can write

\begin{align*} I_{13} & \leq C\Gamma \rho \int_{B_1'} |z'|^2 \int_{1/2}^{\eta/\rho} \frac{{\rm d}z_N}{(|z_N|^2 + |z'|^2)^{\frac{N + 2s} {2}}}{\rm d}z' \\ & \leq C\Gamma \rho \int_{B_1'} |z'|^2 |z'|^{-(N + 2\,s) + 1} \int_{0}^{+\infty} \frac{{\rm d}t}{(t^2 + 1)^{\frac{N + 2s} {2}}}{\rm d}z'. \end{align*}

A similar lower bound can be easily obtained, from which we conclude that $-C \Gamma \rho \leq I_{13} \leq C\Gamma \rho.$ Gathering this estimate together with (3.12) and (3.11) lead us to

(3.13)\begin{equation} I_1 = O(\rho^{\tau + 1}). \end{equation}

Now we proceed with $I_{2}$ in (3.9). We write

\begin{align*} & I_2 = I_{21} + I_{22}, \ \mbox{with} \\ & A_{21} = B_{\eta/\sqrt{\rho}}' \setminus B_1' \times \left(-\frac{\eta}{\rho}, -\frac{1}{2}\right), \ A_{22} = B_{\eta/\sqrt{\rho}}' \setminus B_1' \times \left(-\frac{1}{2}, \frac{\eta}{\rho}\right), \end{align*}

where we have adopted a similar notation as in (3.9).

We start by noticing that $|\tilde \psi (z')| \leq C \eta ^2$ when $|z'| \leq \eta /\sqrt {\rho }$, and from here, fixing $\eta > 0$ universally small, we can write

\begin{align*} I_{21} & \leq \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} \int_{\tilde \psi_-(z') - 1}^{{-}1} (1 + z_N - \tilde \psi(z'))^\tau K^\rho(z){\rm d}z \\ & \qquad+ \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} \int_{\tilde \psi_+(z') - 1}^{{-}1/2} [(1 + z_N - \tilde \psi_+(z'))^\tau - (1 + z_N)^\tau]K^\rho(z){\rm d}z \\ & \quad= : I_{211} + I_{212}. \end{align*}

For $I_{211}$, we have

\begin{align*} I_{211} \leq & C \Gamma \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} |z'|^{-(N + 2s)} \int_{\tilde \psi_-(z') - 1}^{{-}1} (1 + z_N - \tilde \psi(z'))^\tau {\rm d}z_N {\rm d}z' \\ \leq & C \Gamma \rho^{\tau + 1} \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} |z'|^{-(N + 2\,s) + 2(\tau + 1)} {\rm d}z' \\ \leq & \frac{C}{1 - 2s}\Gamma \rho^{1 + \tau}. \end{align*}

For $I_{212}$, by taking $\eta$ small enough but independent of $\rho$, we have

\begin{align*} I_{21} & \leq C \Lambda \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} \frac{1}{(1 + |z'|)^{N+2s}} \int_{\tilde \psi_+(z') - 1}^{{-}1/2} [(1 + z_N - \tilde \psi_+(z'))^\tau\\ & \quad - (1 + z_N)^\tau] {\rm d}z_N {\rm d}z' \\ & \leq C \frac{\Lambda}{1 + \tau} \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} \frac{1}{(1 + |z'|)^{N+2s}} \left(\left(\frac{1}{2} - \tilde \psi_+(z')\right)^{\tau + 1}\right.\\ & \quad \left.- \,\left(\frac{1}{2}\right)^{\tau + 1} + \tilde \psi_+(z')^{\tau + 1}\right) {\rm d}z' \\ & \leq C \frac{\Lambda}{1 + \tau} \Bigg{(} \rho \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} \frac{|z'|^2}{(1 + |z'|)^{N+2s}} {\rm d}z'\\& \quad + \rho^{\tau + 1} \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} \frac{|z'|^{2(\tau + 1)}}{(1 + |z'|)^{N+2s}} {\rm d}z' \Bigg{)} \\ & \leq C \rho^{1 + \tau}, \end{align*}

and from here, collecting the above estimates, we conclude that $I_{21} \leq C \Gamma \rho ^{1 + \tau }$. In a similar fashion as in (3.10) but applied to $I_{21}$, we arrive at the estimate

(3.14)\begin{equation} I_{21} = O(\rho^{1 + \tau}). \end{equation}

For $I_{22}$, taking $\eta$ universally small, there exist $0 < c,\, C < +\infty$ such that

\[ c (1 + z_N)^{\tau - 1} |\tilde \psi_+(z')| \leq (1 + z_N - \tilde \psi_+(z'))^\tau - (1 + z_N)^\tau = C (1 + z_N)^{\tau - 1} |\tilde \psi_+(z')|, \]

from which we can write

\[ I_{22}= O(1) \int_{B_{\eta/\sqrt{\rho}}' \setminus B_1'} \frac{\rho |z'|^2}{(1 + |z'|)^{N+2s}} \int_{{-}1/2}^{+\infty} (1 + z_N)^{\tau - 1} {\rm d}z_N {\rm d}z', \]

and from here we conclude that $I_{22} = O(\rho )$. This together with (3.14) lead us to

(3.15)\begin{equation} I_{2} = O(\rho^{1 + \tau}). \end{equation}

Now we deal with $I_3$. This time we consider the splitting

\begin{align*} & I_3 = I_{31} + I_{32}, \ \mbox{with} \\ & A_{31} = \{ (z', z_N) : z' \in B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}', \ z_N \in (-\eta/\rho, \tilde \psi_+(z') + 1) \}, \\ & A_{32} = \{ (z', z_N) : z' \in B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}', \ z_N \in (\tilde \psi_+(z') + 1, \eta/\rho) \}, \end{align*}

where we have adopted the notation in (3.9).

For $I_{31}$, we see that

\begin{align*} I_{31} & \leq \Lambda \int_{B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}'} \int_{\tilde \psi_+(z') - 1}^{\tilde \psi_+(z') + 1} \frac{(1 + z_N - \tilde \psi_+(z'))^\tau}{|z'|^{N + 2s}}{\rm d}z_N {\rm d}z' \\ & \leq C\Lambda \int_{B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}'} \frac{1}{|z'|^{N + 2s}} \int_{0}^{2} (1 + t)^\tau {\rm d}t {\rm d}z' \\ & \leq C\Lambda \int_{\eta/\sqrt{\rho}}^{\eta/\rho} r^{-(2 + 2s)}{\rm d}r, \end{align*}

from which we conclude that $I_{31} \leq C\Lambda \rho ^{s + 1/2}$.

On the other hand, since $\tau < 0$, for $I_{32}$, we can write

\begin{align*} I_{32} & \leq \Lambda \int_{B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}'} \int_{\tilde \psi_+(z') + 1}^{+\infty} \frac{(1 + z_N - \tilde \psi_+(z'))^\tau}{(z_N^2 + |z'|^2)^{(N + 2s)/2}}{\rm d}z_N {\rm d}z' \\ & \leq \Lambda \int_{B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}'} \int_{1}^{+\infty} \frac{1}{(z_N^2 + |z'|^2)^{(N + 2s)/2}}{\rm d}z_N {\rm d}z' \\ & \leq \Lambda \int_{B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}'} \frac{1}{|z'|^{N + 2\,s - 1}} \int_{0}^{+\infty} \frac{1}{(t^2 + 1)^{(N + 2s)/2}}{\rm d}t {\rm d}z', \end{align*}

and from here, we conclude that

\[ I_{32} \leq C\Lambda \int_{\eta/\sqrt{\rho}}^{+\infty} r^{-(1 + 2s)} {\rm d}r \leq C \Lambda \rho^s. \]

Then, collecting the previous estimates, we conclude that

\[ I_3 \leq C \rho^s. \]

For the lower bound, we see that

\[ I_3 \geq{-} \int_{B_{\eta/\rho}' \setminus B_{\eta/\sqrt{\rho}}'} \int_{-\frac{\eta}{\rho}}^{{-}1} (1 + z_N)^\tau K^\rho(z){\rm d}z, \]

and arguing similarly as before, we arrive at $I_3 \geq -C \rho ^s$. Hence, we conclude that $I_3 = O(\rho ^s)$. Using this estimate, together with (3.13) and (3.15) and replacing them into (3.9), we conclude that

(3.16)\begin{equation} \tilde I = \rho^{\tau - 2s} (O(\rho^{1 + \tau}) + O(\rho^{s})), \end{equation}

where the $O$-term depends on $N,\, s,\, \Omega,\, \eta,\, \frac {1}{1 + \tau }$.

2. Estimate for $I_0$. Rescaling, we have

\begin{align*} I_0 & = \rho^{\tau - 2s} \mathrm{P.V.} \int_{Q_{\eta/\rho}} [(1 + z_N)_+^\tau - 1] K^\rho(z){\rm d}z \\ & = \rho^{\tau - 2s} \mathrm{P.V.} \int_{\mathbb{R}^N} [(1 + z_N)_+^\tau - 1] K^\rho (z){\rm d}z + \rho^{\tau - 2s} \int_{Q_{\eta/\rho}^c} [(1 + z_N)_+^\tau - 1] K^\rho(z){\rm d}z\\ & = : \rho^{\tau - 2s} c_K(\rho, \tau) + I_{01}. \end{align*}

Notice that

\[{-}C \Gamma \rho^{\tau - 2s} \rho^{2s} \leq I_{01}. \]

On the other hand

\begin{align*} I_{01} & \leq \Gamma \rho^{\tau - 2s} \int_{Q_{\eta/\rho}^c \cap \{ z_N >{-}1 \} } (1 + z_N)^{\tau} |z|^{-(N + 2s)} {\rm d}z \\ & \leq C\Gamma \rho^{- 2s} \int_{B_{\eta/\rho}'} \int_{\eta/\rho}^{+\infty} |z|^{-(N + 2s)}{\rm d}z + C\Gamma \rho^{\tau - 2s} \frac{1}{1 + \tau}\int_{B_{\eta/\rho}'^c} |z'|^{-(N + 2s)}{\rm d}z' \\ & \leq C \Gamma (1 + \rho^{1 + \tau}). \end{align*}

Summarizing, we have

\[ I_{01} = O(\rho^\tau). \]

Joining this together with (3.16), we conclude the result.

4. Proof of theorem 1.1

This section is entirely devoted to the proof of theorem 1.1. We recall that for a given family of kernels $\mathcal {K}$ satisfying (1.7), we consider the extremal operators associated with this family as

\[ \mathcal{M}^+ u(x) = \sup_{K \in \mathcal{K}} L_K u(x), \quad \mathcal{M}^{-} u(x) = \inf_{K \in \mathcal{K}} L_K u(x). \]

For every operator with form (1.9), each admissible function $u$ and $x \in \mathbb {R}^N$, we have

(4.1)\begin{equation} \mathcal{M}^{-} u(x) \leq\mathcal{I}(v+u)(x)-\mathcal{I} v(x) \leq \mathcal{M}^+ u(x), \end{equation}

see [Reference Caffarelli and Silvestre6].

Proof of theorem 1.1. Proof of theorem 1.1

We prove each case separately. We give the general remark that both sub and supersolutions we construct here are strict.

Case 1: For $t > 0$ and $s-1 < \gamma < 2\,s - 1< s$ to be fixed, we consider

\[ U^-_t=td^{s-1} - C_1d^{\gamma}, \]

for some $C_1>0$ to be chosen.

Then, using proposition 3.1 together with assumption (1.11) , we have $c(s\!-\!1) = 0$ and $c^+(\gamma )<0$. Then using (4.1), we find

\begin{align*} -\mathcal{I} U^-_t +|DU^-_t|^p & \leq{-}\mathcal{I}(td^{s-1})+C_1 \mathcal{M}^+(d^\gamma)+|DU^-_t|^p \\ & \leq tO(d^{{-}1})+C_1c^+(\gamma) d^{\gamma-2s}+|(t (s-1) d^{s - 2} - C_1 \gamma d^{\gamma - 1}) Dd|^p \\ & = tO(d^{{-}1})+C_1c^+(\gamma) d^{\gamma-2s} \\ & \quad+|t (s-1)|^p d^{(s - 2)p}|1 - C_1 t^{{-}1} (s-1)^{{-}1}\gamma d^{\gamma + 1 - s}|^p. \end{align*}

Since $p < p_2$, we can take $\gamma > s - 1$ such that $(s - 2)p > \gamma -2\,s$ (notice that if $p < p_0$, then $\gamma$ can be taken positive).

For such a $\gamma$, we take

\[ \bar C_1 = |t (s-1)|^p/|c^+(\gamma)|, \]

to conclude that there exists $\bar c_1 > 0$ such that, for every $\epsilon > 0$, we take $C_1 = \bar C_1 - \epsilon$ in the expression above to obtain that

\[ -\mathcal{I} U^-_t +|DU^-_t|^p \leq tO(d^{{-}1}) - \bar c_1 \epsilon d^{\gamma - 2s}, \]

for each $d \leq d_\epsilon$ for some $d$ small enough in terms of $\epsilon,\, t,\, C_1,\, \gamma.$ Since $s > 1/2$, and $\gamma < 2s-1$, we can take $d_\epsilon$ smaller to conclude that

\[ -\mathcal{I} U^-_t +|DU^-_t|^p +{\leq}{-} \frac{\bar c_1 \epsilon}{2} d^{\gamma - 2s} \leq{-} \| f \|_\infty , \quad \mbox{for} \ d(x) \leq d_\epsilon, \]

and therefore we have constructed a subsolution near the boundary. A straightforward computation tells us that the function

\[ U_t^-{-} C\chi_\Omega \]

is a viscosity subsolution to the problem in $\Omega$, when $C = C_\epsilon$ is taken large . Thus, we have constructed the subsolution.

In a similar way, we can construct a supersolution in $\Omega$ with the form

\[ U^+_t=td^{s-1} + (\bar C_1+ \epsilon) d^{\gamma} + \mathcal{C}_\epsilon \chi_\Omega. \]

We can apply the Perron method from § 2 to conclude the existence of a solution to problem (P) satisfying the desired rate near the boundary, where $\theta = \gamma + 1 - s$. If $p < p_0$, then we can take $\gamma > 0$ in the above analysis.

Similar arguments hold for $t<0$.

Case 2: Let $\beta$ as in the statement of the theorem. Hence, $-1< \beta < s-1$ and for some $s-1 < \gamma < 0$ and $T,\, C_1 > 0$ to be fixed, denote

\[ U^-{=} Td^\beta - d^\gamma. \]

We invoke proposition 3.1 again, noticing that if (1.11) holds, and by choice of $\beta$ and range of $\gamma$, we have $c(\beta )<0$ and $c^+(\gamma )<0$. Then using again (4.1), we find

\begin{align*} -\mathcal{I} U^-{+}|DU^-|^p & \leq{-}c(\beta) T d^{\beta - 2s} + O(d^{\beta - s}) + c^+(\gamma) d^{\gamma - 2s} \\ & \quad + (T |\beta|)^p d^{(\beta - 1)p}|1 - T^{{-}1} \beta^{{-}1} \gamma d^{\gamma - \beta}|^p \\ & \leq ({-}c(\beta) T + (T |\beta|)^p) d^{\beta - 2s} + O(d^{\beta - s}) + c^+(\gamma) d^{\gamma - 2s}. \end{align*}

Then, taking $T = \bar T$ such that

\[ c(\beta) \bar T = \bar T^p |\beta|^p \]

concludes that

\[ -\mathcal{I} U^-{+} |DU^-|^p \leq O(d^{\beta - s})+ c^+(\gamma) d^{\gamma - 2s}, \]

in a neighbourhood of $\partial \Omega$. Note that it is possible to fix $\gamma$ close enough to $s - 1$ in order to have $\gamma - 2\,s < \beta - s$, and then, for each point close to the boundary, we have

\[ -\mathcal{I} U^-{+} |DU^-|^p \leq \frac{c^+(\gamma)}{2} d^{\gamma - 2s}. \]

We conclude the existence of a subsolution in $\Omega$ in the same way as before.

For the supersolution, we consider

\[ U^+{=} \bar T d^{\beta} + d^\gamma, \]

and proceed as before.

Case 3: The proof is similar to the previous case, but we provide the details for completeness. Since $p>p_2$, we have $s-1<\beta <0$. For $T,\, C>0$ and $\gamma \in (0,\,2s-1)$ to be fixed, we define

\[ U ={-}Td^{\beta} - Cd^\gamma. \]

Notice that in this case $c(\gamma ) < 0$.

As before, writing $\tilde{\mathcal {I}}$ as the operator $- \mathcal {I} (- \cdot )$, using proposition 3.1 and (1.11), together with (4.1), we have that

\begin{align*} -\mathcal{I} U + |DU|^p & \leq T \tilde {\mathcal{I}} (d^\beta ) + C \mathcal{M}^+ (d^\gamma) + |T \beta d^{\beta - 1} + C\gamma d^{\gamma - 1}|^p \\ & \leq d^{\beta - 2s} (T \tilde c(\beta) + O(d^s)) + C c^+(\gamma) d^{\gamma - 2s} + |T \beta d^{\beta - 1} + C\gamma d^{\gamma - 1}|^p, \end{align*}

where $\tilde c(\beta ) < 0$ since $s-1<\beta <0$. At this point, we notice that since $\beta < 0$ and $C,\, \gamma > 0$, we have

\[ |T \beta d^{\beta - 1} + C\gamma d^{\gamma - 1}| \leq T |\beta| d^{\beta - 1}, \]

for all $d = d(x)$ small enough. Using this, we get that

(4.2)\begin{equation} -\mathcal{I} U + |DU|^p \leq d^{\beta - 2s} (T \tilde c(\beta) + O(d^s)) + C c^+(\gamma) d^{\gamma - 2s} + T^p |\beta|^p d^{(\beta - 1)p}, \end{equation}

and fixing $T = T^* := (-\tilde c(\beta ) |\beta |^p)^{\frac {1}{p - 1}} > 0$, since $\beta - 2\,s = (\beta - 1)p$, we conclude that

\[ -\mathcal{I} U + |DU|^p \leq O(d^{\beta - s}) + Cc^+(\gamma) d^{\gamma - 2s} . \]

Since we also have $\beta - s>-1>\gamma -2\,s$, by similar arguments used before, we can fix $C > 0$ large enough to find that $U$ is a subsolution.

Now, for $C > 0$ and $\gamma > 0$, consider the function $V = -T^* d^\beta + C d^\gamma$. Notice that $U \leq V$. By a similar computation as above, we get that $V$ is a supersolution for the problem.

Remark 4.1 Some remarks concerning the asymptotics as $s \to 1$ for $p \in (1,\,2)$ fixed. By the construction of the barriers in theorem 1.1, we see that the one-parameter solutions, and the negative scale solutions (cases 1 and 3) shall always converge to the unique, bounded solution to problem (1.3) with $u = 0$ on $\partial \Omega$. On the other hand, positive scale solutions (case 2) converge to the unique large solution to (1.3). This is a consequence of well-known stability results of viscosity solutions, the estimates for each of the solutions in cases 1, 2 and 3 and the uniqueness of the limit equation.

Concerning ‘critical’ cases, we do not know if there exist blow-up solutions for the case $p=p_2$ (and $s$ fixed). In fact, the barriers constructed provide estimates that imply that the solutions found for $p < p_2$ and $p > p_2$ tend to the bounded solution of the limit problem as $p \to p_2^\pm$. For the case $p = 2\,s$, it was neither possible to construct large solutions by approximation as $p \to 2s^-$ ($s$ fixed). The local setting suggests to search for large solutions with logarithmic profile, but we leave that analysis for a future work.

5. Extensions

In this section, we provide a discussion about possible extensions of theorem 1.1 for more general operators and data.

First recall that the estimates found in theorem 1.1 are crucial in order to prove theorem 1.1. In order to obtain a similar result than theorem 1.1 for general kernels, it suffices that the family of kernels $K$ satisfies the requirement

(5.1)\begin{equation} c(d(x), \tau) = c(\tau) + O(d(x)^\alpha), \end{equation}

for some $\alpha > 0$ and some $c: (-1,\, 2\,s) \to \mathbb {R}$ that is independent of $x$.

For example, this condition holds when $K = K_a$ has the form

(5.2)\begin{equation} K_a(z) = \frac{a(z)}{|z|^{N + 2s}}, \quad z \neq 0, \end{equation}

where $a: \mathbb {R}^N \to \mathbb {R}$ is a nonnegative, measurable function with a uniform modulus of continuity at $z = 0$. More specifically, if we consider the class $a \in \mathcal {A}$ such that there exists $C,\, r > 0$ and $\alpha ' > 0$ such that

(5.3)\begin{equation} \sup_{a \in \mathcal{A} } \{ |a(z) - a(0) | \} \leq C|z|^{\alpha'}, \quad \mbox{for all} \ z \in B_r, \end{equation}

and use the class of kernels $\{ K_a \}_{a \in \mathcal {A}}$, then expansion (5.1) holds for some $0 < \alpha < \alpha '$.

In fact, the key step to arrive at (5.1) comes by the expansion, for $K = K_a$, given by

\begin{align*} _K(d(x), \tau) & ={-}C_{N, s}^{{-}1} a(0) (-\Delta)^s (x_N)_+^\tau (e_N) \\ & \quad+ \mathrm{P.V.} \int_{\mathbb{R}^N} [(1 + z_N)_+^\tau - 1] \frac{(a(d(x) z) - a(0)){\rm d}z}{|z|^{N + 2s}}. \end{align*}

The first term in the right-hand side is independent of $\rho$. The regularity assumption on $a$ allows us to control the second term in the righ-hand side as an error term of order $O(\rho ^{\alpha })$ for some $\alpha < \alpha '$. For this, we divide the integrand as $\mathbb {R}^N = B_{R_\rho } \cup B_{R_\rho }^c$ for some $R_\rho$ large depending on $\rho$. The exterior part is controlled by the tails of the kernel and the boundedness of $a$. The inner part requires the modulus of continuity $a$, for which we have the restriction that $R_\rho << \rho ^{-1}$. See Lemma 3.1 in [Reference Dávila, Quaas and Topp11] for details.

Using this, for the class of kernels satisfying (5.1)(5.3), we have

\[ \mathcal{I} d^\tau(x) = d^{\tau - 2s}(x) \inf_{i \in I} \sup_{j \in J} \Bigg{(} c_{ij}(\tau, N, s) + O(d(x)^\alpha) + O(d(x)^s)\Bigg{)}, \]

where

\[ c_{ij}(\tau, N, s) ={-}\tilde C_{N,s} a_{ij}(0)(-\Delta)^s_{\mathbb{R}} (x)_+^\tau (1). \]

This time, the Isaacs form of $\mathcal {I}$ is indexed by $\{ a_{ij} \}_{i \in I, j \in J}$ through the kernels $K_{a_{ij}}$ defined in (5.2). As in (3.8), we have denoted as $(-\Delta )_{\mathbb {R}}^s$ the one-dimensional fractional Laplacian of order $2\,s$ and since $a_{ij}(0) \geq a_0 > 0$, we can use the estimate regarding $\mathcal {I}$ as in the proof of theorem 1.1.

Concerning the exterior data, we reduce the problem to the homogeneous case by inserting it in the source term. Namely, for $\varphi \in L^1_\omega (\Omega ^c)$, denote

\[ \tilde\varphi(x)=\left\{\begin{array}{@{}ll} 0 & \text{for }x\in\Omega,\\ \varphi(x) & \text{for }x\in\Omega^c. \end{array}\right. \]

Now, if $u$ is a solution of

(5.4)\begin{equation} \left \{ \begin{array}{@{}rl} - \mathcal{I} u + |Du|^p + \lambda u = f \quad & \mbox{in} \ \Omega, \\ u = \varphi \quad & \mbox{in} \ \Omega^c, \end{array} \right . \end{equation}

then $\tilde u=u-\tilde \varphi$ is a solution of

(5.5)\begin{equation} \left\{\begin{array}{@{}rl} - \mathcal{I}_\varphi \tilde u+|D \tilde u|^p+\lambda \tilde u=f, & \text{in }\Omega,\\ \tilde u=0, & \text{in }\Omega^c \end{array}\right. \end{equation}

where $\mathcal {I}_\varphi v = \mathcal {I}(v+\tilde \varphi )$ and satisfies the extremal inequalities

\[ \mathcal{I} v+\mathcal{M}^-\tilde\varphi \leq \mathcal{I}_\varphi v\leq \mathcal{I} v+\mathcal{M}^+\tilde\varphi. \]

Observe that $\mathcal {M}^\pm \varphi (x)$ is well defined for $x \in \Omega$. Moreover, when $\varphi$ is bounded in a neighbourhood of $\partial \Omega$, then $x \mapsto \mathcal {M}^\pm \tilde \varphi (x)$ would be of order $d^{-2s}(x)$.

Problem (5.5) allows us to consider zero exterior condition at the expense of imposing the extra assumptions on the ‘new’ right-hand side involving $f$ and the extremal operators evaluated at $\tilde \varphi$. More specifically, we will have to deal with a function $\tilde u$ satisfying

(5.6)\begin{equation} f - \mathcal{M}^+ \tilde \varphi \leq{-}\mathcal{I} \tilde u+|D \tilde u|^p+\lambda \tilde u \leq f - \mathcal{M}^- \tilde \varphi \quad \mbox{in} \ \Omega, \end{equation}

in the viscosity sense, and vanishing in $\Omega ^c$.

By the method used here, we can consider an unbounded right-hand side $f$ satisfying the following blow-up conditions at the boundary. $(H1)$ There exist $\eta >0$ and $C_f \geq 0$ such that

\[{-}C_f \leq \liminf\limits_{x\to \partial \Omega} d(x)^{s+1-\eta} f(x) \leq \limsup\limits_{x\to \partial \Omega} d(x)^{s+1-\eta} f(x) \leq C_f. \]

$(H2)$ There exist $C_f > 0$ such that

\[{-}C_f \leq \liminf\limits_{x\to \partial \Omega} d(x)^{2s} f(x) \leq \limsup\limits_{x\to \partial \Omega} d(x)^{2s} f(x) \leq C_f. \]

Next, we present our version of theorem 1.1 in the current generality.

Theorem 5.1 Let $s \in (1/2,\,1),$ $0 < p < 2\,s,$ $\Omega \subset \mathbb {R}^N$ be a bounded domain with $C^2$ boundary, $f\in C(\Omega )$ and $\varphi \in C(\Omega ^c) \cap L_\omega ^1(\Omega ^c)$.

Assume $\mathcal {K}$ as in theorem 1.1, or with form (5.2)(5.3). Let $\mathcal {I}$ be a nonlinear operator with form (1.9), and let $p_i$ be defined as in (1.12), $i=0,\,1,\,2$. Assume that $\lambda >-\lambda _0(\mathcal {I})$.

Then, we have the following existence results:

1. One parameter family of solutions (close to $s$-harmonic): If $0 < p < p_2$ and $f$ satisfies $(H1)$ , there exists $\sigma > 0$ and a family of solutions $\{ u_t \}_{t \in \mathbb {R}, t \neq 0} \subset C^{\sigma }(\Omega )$ to (5.4), such that for each $t$ we have

\[ d^{1 - s}(x) u_t(x) - t = O(d^{\theta}(x)), \]

for some $\theta > 0$ depending on $p$. In particular, if $t_1 < t_2$, then

\[ u_{t_1} < u_{t_2} \quad \mbox{in} \ \Omega. \]

2. Positive scale solution: If $p_1 < p < p_2$ and $f$ satisfies $(H1)$, then there exist $\sigma > 0$ and a constant $T > 0$ and a function $u \in C^{\sigma }(\Omega )$ solving (5.4) such that

\[ d(x)^{-\beta}u(x) - T = O(d(x)^\theta), \]

for some $\theta > 0$.

3. Negative scale solution: For $p_2 < p < 2\,s$ and $f$ satisfies $(H2),$ then there exist $\sigma > 0,$ $T > 0$ and a solution $u \in C^{\sigma }(\Omega )$ of (5.4) such that

\[ d^{-\beta}(x)u(x) + T = O(d(x)^\theta), \]

for some $\theta > 0$.

The proof of theorem 5.1 follows the same lines as theorem 1.1. We first observe that there exists $C_\varphi > 0$ such that, for all $x \in \Omega$ close to the boundary, we have

\[ |\mathcal{M}^+ \tilde \varphi(x)|, |\mathcal{M}^- \tilde \varphi(x)| \leq C_\varphi d(x)^{{-}2s}. \]

In case 1, we consider a lower barrier with the form

\[ U_t(x) = td^{s - 1}(x) - C_1 d^\gamma(x) - C_2 \chi_\Omega(x), \]

where in view of inequalities (5.6), we properly choose $\gamma \in (s-1,\,0)$ close to $s-1$ in terms of $\eta$, $C_1 > 0$ large enough in terms of $C_f$, and $C_2 > 0$ large in terms of $C_\varphi$.

For case 2, we consider a lower barrier with the form

\[ U(x) = T d^\beta(x) - C_1 d^\gamma(x) - C_2 \chi_\Omega(x), \]

with $T > 0$ as in the proof of theorem 1.1, and a choice for $\gamma \in (s-1,\, 0)$, $C_1 > 0$ and $C_2 > 0$ similar to the mentioned in the paragraph above.

Finally, for case 3., we consider

\[ U(x) ={-}Td^\beta(x) - C_1 d^\gamma(x) - C_2 \chi_\Omega(x), \]

with $T > 0$ as in the proof of theorem 1.1, $\gamma \in (0,\, 2\,s - 1)$, and $C_1,\, C_2 > 0$ similarly as above.

Upper bounds for all the cases can be constructed in a similar fashion.

Acknowledgements

G. D. was partially supported by Fondecyt Grant No. 190209, and Fondecyt Grant No. 1230635. A. Q. was partially supported by Fondecyt Grant No. 1190282. E. T. was partially supported by Fondecyt Grant No. 1201897.

References

Abatangelo, N.. Large $s$-harmonic functions and boundary blow-up solutions for the fractional laplacian. Discrete Contin. Dyn. Syst. 35 (2015), 55555607.10.3934/dcds.2015.35.5555CrossRefGoogle Scholar
Alarcón, S., García-Melián, J. and Quaas, A.. Keller–Osserman type conditions for some elliptic problems with gradient terms. J. Differ. Equ. 252 (2012), 886914.CrossRefGoogle Scholar
Barles, G., Chasseigne, E. and Imbert, C.. On the Dirichlet problem for second order elliptic integro-differential equations. Indiana Univ. Math. J. 57 (2008), 213246.CrossRefGoogle Scholar
Guy, B., Emmanuel, C. and Cyril, I.. Hölder continuity of solutions of second-order non-linear elliptic integro-differential equations. J. Eur. Math. Soc. (JEMS) 13 (2011), 126.Google Scholar
Bertoin, J.. Levy processes, Cambridge Tracts in Mathematics, Vol. 121 (Cambridge University Press, Cambridge, 1996).Google Scholar
Caffarelli, L. and Silvestre, L.. Regularity theory for fully nonlinear integro-differential equations. Communicat. Pure Appl. Math. 62 (2009), 597638.CrossRefGoogle Scholar
Chen, H. and Lv, G.. Boundary blow-up solutions to nonlocal elliptic equations with gradient nonlinearity. Communicat. Contemporary Math. 19 (2017), 1650051.10.1142/S0219199716500516CrossRefGoogle Scholar
Chen, H., Hajaiej, Y. and Wang, H.. Boundary blow-up solutions to fractional elliptic equations in a measure framework. Discrete Continuous Dynam. Syst. - A 36 (2016), 18811903.Google Scholar
Chen, H., Felmer, P. and Quaas, A.. Large solutions to elliptic equations involving fractional Laplacian. Ann. Inst. Henri Poincare Analyse non lineaire 32 (2015), 11991228.10.1016/j.anihpc.2014.08.001CrossRefGoogle Scholar
Chen, H. and Veron, L.. Semilinear fractional elliptic equations with gradient nonlinearity involving measures. J. Funct. Anal. 266 (2014), 54675492.10.1016/j.jfa.2013.11.009CrossRefGoogle Scholar
Dávila, G., Quaas, A. and Topp, E.. Continuous viscosity solutions for nonlocal Dirichlet problems with coercive gradient terms. Math. Annalen 369 (2017), 12111236.CrossRefGoogle Scholar
Di Neza, E., Palatucci, G. and Valdinoci, E.. Hitchhiker's guide to the fractional Sobolev spaces. Bull. Sci. Math. 136 (2012), 521573.CrossRefGoogle Scholar
Lasry, J. M. and Lions, P. L.. Nonlinear elliptic equations with singular boundary conditions and stochastic control with state constraints. Math. Ann. 283 (1989), 583630.CrossRefGoogle Scholar
Keller, J. B.. On solutions of $\Delta u = f (u)$. Comm. Pure Appl. Math. 10 (1957), 503510.10.1002/cpa.3160100402CrossRefGoogle Scholar
Osserman, R.. On the inequality $\Delta u = f (u)$. Pacific J. Math. 7 (1957), 16411647.CrossRefGoogle Scholar
Pham, H.. Optimal stopping of controlled jump diffusion processes: a viscosity solution approach. J. Math. Systems Estim. Control 8 (1998), 27.Google Scholar
Ros-Oton, X. and Serra, J.. Boundary regularity for fully nonlinear integro-differential equations. Duke Math. J. 165 (2016), 20792154.10.1215/00127094-3476700CrossRefGoogle Scholar
Sato, K.-I.. Levy Processes and Infinitely Divisible Distributions (Cambridge University Press, Cambridge, 1999).Google Scholar
Topp, E.. Existence and uniqueness for integro-differential equations with dominating drift terms. Comm. Partial Differ. Equ. 39 (2014), 15231554.10.1080/03605302.2014.900567CrossRefGoogle Scholar