Almost Lipschitz regularity for solutions of elliptic equations with discontinuous coefficients

Abstract

We are interested in the local higher integrability of solutions to elliptic equations with linear growth of the form

$$-\text{ div}A(x,Du)=f(x).$$

Under a Besov regularity assumption both on the partial map $x\mapsto A(x,\xi)$ and the datum $f$, we prove that the solutions are almost Lipschitz continuous, i.e. their gradients belong locally to $L^{q}$, for any finite exponent $q$. In turn, solutions are locally $\gamma$-Hölder continuous, for every $\gamma\in(0,1)$. The difficulty arising from the lack of an explicit second variation for the problem is overcome by testing the equation with a function proportional to a power of the finite difference quotient of the solution. To the best of our knowledge, this technique is used in this context for the first time. We also provide an example showing the sharpness of our result in the scale of Lebesgue spaces.

Keywords: Besov spaces, Gradient regularity, Hölder regularity.

MSC 2020: 49N60; 46E35; 35B65; 35J25.

1 Introduction

The aim of this paper is to prove the almost Lipschitz continuity of the local weak solutions of equations of the form

$$-\text{ div}A(x,Du)=f(x)\qquad\text{in }\Omega,$$ (1.1)

where $\Omega\subset\mathbb{R}^{n}$ is an open set, $n\geq 2$, $A:\Omega\times\mathbb{R}^{n}\to\mathbb{R}^{n}$ is a Carathéodory function such that $A(\cdot,0)\in L^{2}_{\textrm{loc}}(\Omega)$ and $f\in B^{\alpha}_{\frac{n}{2\alpha},p,\textrm{loc}}(\Omega)$, for $\alpha\in(0,1)$ and $p\in\left[1,\frac{2n}{n-2\alpha}\right]$. We assume that there exist positive constants $\nu,L$ such that

$$\langle A(x,\xi)-A(x,\eta),\xi-\eta\rangle\geq\nu|\xi-\eta|^{2}$$ (A1)

$$|A(x,\xi)-A(x,\eta)|\leq L|\xi-\eta|$$ (A2)

for a.e. $x\in\Omega$ and all $\xi,\eta\in\mathbb{R}^{n}$. Concerning the dependence on the $x$-variable, we assume that there exists a sequence of measurable non-negative functions $g_{k}\in L^{\frac{n}{\alpha}}_{\textrm{loc}}(\Omega)$ such that for every $\Omega^{\prime}\Subset\Omega$

$$\sum_{k=0}^{\infty}\|g_{k}\|^{p}_{L^{\frac{n}{\alpha}}(\Omega^{\prime})}<\infty,$$ (1.2)

for $1\leq p\leq\frac{2n}{n-2\alpha}$, and at the same time there exists $\mu\in[0,1]$ such that

$$|A(x,\xi)-A(y,\xi)|\leq|x-y|^{\alpha}(g_{k}(x)+g_{k}(y))(\mu^{2}+|\xi|^{2})^{\frac{1}{2}},$$ (A3)

for a.e. $x,y\in\Omega$ such that $2^{-k}\text{diam}(\Omega)\leq|x-y|<2^{-k+1}\text{diam}(\Omega)$ and for every $\xi\in\mathbb{R}^{n}$. We note that assumption (A3) tells us that the map $x\mapsto A(x,\xi)$ belongs locally to the Besov space $B^{\alpha}_{\frac{n}{\alpha},p}(\Omega)$ (see Theorem 4.9 below).

Definition 1.1.

A function $u\in W_{\textrm{loc}}^{1,2}(\Omega)$ is a local weak solution of equation (1.1) if and only if, for any test function $\varphi\in\mathcal{C}_{0}^{\infty}(\Omega)$, the following integral identity holds:

$$\int_{\Omega}\langle A(x,Du),D\varphi\rangle\,dx\,=\,\int_{\Omega}f\varphi\,dx.$$

A main topic in the analysis of equations as in (1.1) is to investigate how the regularity of the right-hand side $f$ transfers to the solutions under suitable assumptions on the coefficients of the operator $A$.
It is known that for solutions of the $p$-Laplacian type equations or systems with divergence form right-hand side such as

$$\mathrm{div}(|Du|^{p-2}Du)=\mathrm{div}(|F|^{p-2}F)\,,\quad p>1\,,$$

Calderón-Zygmund estimates of the type

$$F\in L^{q}_{\mathrm{loc}}(\Omega)\quad\Longrightarrow\quad Du\in L^{q}_{\mathrm{loc}}(\Omega)\,,\quad q>p\,,$$ (1.3)

hold. The scalar case was treated in [26], while the vectorial one in [13].
The result in (1.3) is still true for solutions of elliptic equations of the form

$$\mathrm{div}A(x,Du)=\mathrm{div}(|F|^{p-2}F)\,,$$ (1.4)

where the vector field $A$ is assumed to satisfy standard $p$-growth and ellipticity conditions and a continuity dependence on the $x$-variable (see [26]). In the case of systems of the type (1.4), Kristensen and Mingione [33] proved that (1.3) holds, provided $q$ is less than a threshold parameter if $n>2$, or for any $q$ if $n=2$.

Extensions of Calderón-Zygmund estimates to non-uniformly elliptic problems are addressed in [1], for systems with $p(x)$-growth conditions, and in [9], for double phase type equations. We also refer to [38] for a comprehensive Calderón-Zygmund theory in the context of measure data problems. In the setting of Morrey-Campanato spaces we recall the recent paper [40].

Iwaniec and Sbordone [27, 28] developed the theory of the integrability of the gradient of solutions for linear equations with coefficients of vanishing mean oscillation (VMO), which need not be continuous. A first non-linear counterpart for equations of $p$-Laplacian type can be found in [30]; see also [31] for global estimates.
Since then, the regularity theory for solutions to variational problems with VMO coefficients has attracted considerable attention in the last years. Baison et al. [2] established Calderón-Zygmund estimates (1.3) for solutions to equations (1.4) in the linear growth case. In the homogeneous case, that is when $F=0$, Cupini et al. [10] proved that the $W^{1,n}$-regularity of the coefficients implies that the gradient of solutions belongs to $L^{q}$, for any finite exponent $q$ (recall that $W^{1,n}$ embeds into the space of VMO functions). In turn, this gives that solutions are $\gamma$-Hölder continuous, for any exponent $\gamma\in(0,1)$. The Lipschitz continuity cannot be expected without assuming that the coefficients belong to a Sobolev space strictly contained in $W^{1,n}$ ([14, 15, 18, 22]) and the right-hand side $f$ belongs to a function space strictly contained in $L^{n}$ ([7, 35, 36, 41]) or to the Lorentz space $L^{n,\infty}$, if in addition it is radially decreasing ([39]).
Since the works [16, 37], it is known that if the right-hand side $f\in L^{\infty}$ and the operator $A$ is continuous in the $x$-variable, then solutions of (1.1) are locally Hölder continuous with exponent $\gamma$, for any $\gamma\in(0,1)$, but not locally Lipschitz continuous.

Our analysis naturally fits into the setting of problems with VMO coefficients; indeed, the Besov assumption (A3) implies that the map $x\mapsto A(x,\xi)$ is locally uniformly in VMO (see Section 2 for the precise definition).

Regarding variational problems with Besov coefficients, several results concerning the higher fractional differentiability of solutions have been established in recent years (see for example [2, 4, 3, 8, 11, 17, 6, 21, 23, 34]). Here, we are interested in the higher integrability properties of the gradient for solutions $u$ to equation (1.1). More precisely, we establish that $Du\in L^{q}_{\mathrm{loc}}(\Omega)$ for every finite $q>2$, under the hypothesis that the right-hand side belongs to the local Besov space $f\in B^{\alpha}_{\frac{n}{2\alpha},p,\mathrm{loc}}(\Omega)$. It is worth noting that, by virtue of Lemma 4.2 below, this fractional differentiability guarantees the summability condition $f\in L^{\frac{n}{\alpha}}_{\mathrm{loc}}(\Omega)$. Crucially, this setting falls outside the scope of classical Calderón-Zygmund estimates. Our result provides a novel extension of the standard Calderón-Zygmund theory: we demonstrate that arbitrary higher integrability of the gradient can be attained under a mild summability condition on $f$, provided this deficit is counterbalanced by a small fractional differentiability assumption.

The classical strategy for establishing Lipschitz regularity of solutions typically relies on differentiating the underlying equation and choosing a test function proportional to a suitable power of the gradient $Du$. However, in our setting, this standard procedure is precluded: under assumption (A3), the coefficients possess merely fractional differentiability and lack the smoothness required to justify weak differentiation. To overcome this obstruction, we build upon a discrete differentiation approach pioneered by Bögelein et al. [5] for the fractional $p$-Laplacian. The core strategy consists in testing equation (1.1) not with the gradient, but with a carefully chosen function proportional to a power of the finite difference quotient $\tau_{h}u$. By rigorously coupling this technique with a Moser iteration scheme , we establish the almost Lipschitz continuity of solutions. Specifically, we prove that $u\in\mathcal{C}^{0,\gamma}_{\mathrm{loc}}(\Omega)$ for every exponent $0<\gamma<1$. It is crucial to emphasize that the lack of $L^{\infty}_{\mathrm{loc}}$ bounds for the gradient is not a technical limitation of our proof, but rather an intrinsic feature of the problem. Because our coefficients belong to the Besov space $B^{\alpha}_{r,p}$ at the critical exponent $r=\frac{n}{\alpha}$, they may fail to be continuous. Consequently, the almost Lipschitz regularity achieved in Theorem 1.2 cannot be improved to Lipschitz continuity. We validate this sharpness by constructing a counterexample in Section 7.

We now state our main result.

Theorem 1.2.

Let $0<\alpha<1$ and $1\leq p\leq\frac{2n}{n-2\alpha}$. Let $u\in W^{1,2}_{\textrm{loc}}(\Omega)$ be a local weak solution of (1.1) and $f\in B^{\alpha}_{\frac{n}{2\alpha},p,\textrm{loc}}(\Omega)$. Then, $u\in W^{1,q}_{\textrm{loc}}(\Omega)$ for every $q\in[2,+\infty)$ and the following estimate

$$\left(\int_{B_{R/2}}|Du|^{q}dx\right)^{\frac{1}{q}} \leq\tilde{c}\left(\|f\|_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{2R})}+\|u\|_{W^{1,2}(B_{2R})}+1\right)$$

holds for every concentric balls $B_{R/2}\subset B_{2R}\Subset\Omega$ and for a positive constant $\tilde{c}=\tilde{c}(n,q,\alpha,\nu,L,p,R)$.

A straightforward consequence of the previous result is the following

Corollary 1.3.

Let $0<\alpha<1$ and $1\leq p\leq\frac{2n}{n-2\alpha}$. Let $u\in W^{1,2}_{\textrm{loc}}(\Omega)$ be a local weak solution of (1.1) and $f\in B^{\alpha}_{\frac{n}{2\alpha},p,\textrm{loc}}(\Omega)$. Then, $u\in\mathcal{C}^{0,\gamma}_{\textrm{loc}}(\Omega)$, for every $0<\gamma<1$, and we have the following estimate

$$[u]_{\mathcal{C}^{0,\gamma}(B_{R/2})} \leq\tilde{c}\left(\|f\|_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{2R})}+\|u\|_{W^{1,2}(B_{2R})}+1\right),$$

for every concentric balls $B_{R/2}\subset B_{2R}\Subset\Omega$ and for a positive constant $\tilde{c}=\tilde{c}(n,\gamma,\alpha,\nu,L,p,R)$.

A key role in the proof of Theorem 1.2 is to derive suitable estimates on second-order finite differences of the solution $u$. The first step is to establish some a priori estimates for sufficiently regular solutions. Inspired by the work [5], we test the equation (1.1) with

$$\varphi=\tau_{-h}(\eta^{2}|\tau_{h}u|^{q-2}\tau_{h}u),$$

where $q\geq 2$ and $\eta$ is a cut-off function. Thus, we obtain uniform estimates for the $B^{\frac{2\alpha}{q}+1}_{q,\frac{pq}{2}}$ norm of $\tau_{h}(\tau_{h}u)$, for every $q\geq 2$.

The next step is to improve the integrability of the gradient $Du$ from the natural energy space $L^{2}$ to $L^{q}$, for any arbitrarily large exponent $q>2$. To achieve this, we use a bootstrap argument inspired by Moser’s iteration, but suitable adapted to the fractional setting. Specifically, assuming $Du\in L^{q}_{\mathrm{loc}}(\Omega)$ for some $q\geq 2$, our second-order finite difference estimates, combined with Proposition 4.6, yield the refined Besov regularity $Du\in B^{\frac{2\alpha}{q}}_{q,\frac{pq}{2},\mathrm{loc}}(\Omega)$. By invoking sharp Sobolev-type embeddings for Besov spaces, this fractional differentiability translates directly into the higher integrability $Du\in L^{\frac{nq}{n-2\alpha}}_{\mathrm{loc}}(\Omega)$. A key difference from the classical Moser iteration for the Laplace equation is the finite nature of our scheme. Due to the structural properties of the fractional embeddings, the bootstrap cannot be carried out infinitely many times to yield a $L^{\infty}$ bound. Nevertheless, this finite sequence of iterations is perfectly calibrated to reach any targeted higher integrability exponent. Finally, Theorem 1.2 follows by proving that the a priori estimate is preserved in passing to the limit.
We conclude this introduction by describing the organization of the paper. After a brief presentation of the notation and preliminary results, in Section 4 we state some elementary properties and fundamental estimates for Besov spaces. Section 5 is devoted to the proof of Theorem 1.2. Eventually, in Section 6 we prove our main result by using an approximation argument. We conclude exhibiting an example of equation for which the result in Theorem 1.2 is sharp, in the sense that we cannot expect to obtain local Lipschitz continuity of solutions to (1.1) under the assumption (A3). This counterexample was constructed by Jin et al. [29] to show that solutions to elliptic equations with continuous coefficients may not be Lipschitz continuous.

2 Notation and preliminary results

Here, we introduce some notation. We will use the symbols $C$ or $c$ to denote general positive constants. Different occurrences from line to line will be still denoted using the same letters. Relevant dependencies on parameters will be emphasized using parentheses or subscripts. We denote by $B(x,r)=B_{r}(x)=\{y\in\mathbb{R}^{n}:|y-x|<r\}$ the ball centered at $x$ of radius $r$. We shall omit the dependence on the center and on the radius when no confusion arises. For a function $v\in L^{1}(B)$, the symbol

$v_{B}:=\displaystyle\fint_{B}v(x)dx=\dfrac{1}{|B|}\displaystyle\int_{B}v(x)dx$.

will denote the integral mean of the function $v$ over the ball $B$.

We recall the following estimate (see e.g. [19, Lemma 8.3]).

Lemma 2.1.

Let $1<p<+\infty$. There exists a constant $c=c(n,p)>0$ such that

$c^{-1}(|\xi|^{2}+|\eta|^{2})^{\frac{p-2}{2}}|\xi-\eta|^{2}\leq{\left||\xi|^{\frac{p-2}{2}}\xi-|\eta|^{\frac{p-2}{2}}\eta\right|^{2}}\leq c(|\xi|^{2}+|\eta|^{2})^{\frac{p-2}{2}}|\xi-\eta|^{2}$

for any $\xi,\eta\in\mathbb{R}^{n}$.

Now we state a well-known iteration lemma (see [19, Lemma 6.1] for the proof).

Lemma 2.2.

Let $\Phi:[\frac{R}{2},R]\rightarrow\mathbb{R}$ be a bounded nonnegative function, where $R>0$. Assume that for all $\frac{R}{2}\leq r<\rho\leq R$ it holds

$$\Phi(r)\leq\theta\Phi(\rho)+A+\dfrac{B}{(\rho-r)^{2}}+\dfrac{C}{(\rho-r)^{\gamma}},$$

where $\theta\in(0,1)$, $A$, $B$, $C\geq 0$ and $\gamma>0$ are constants. Then there exists a constant $c=c(\theta,\gamma)$ such that

$$\Phi\biggl(\dfrac{R}{2}\biggr)\leq c\biggl(A+\dfrac{B}{R^{2}}+\dfrac{C}{R^{\gamma}}\biggr).$$

As usual, we denote by $\mathcal{C}^{0,\gamma}(\Omega)$, with $0<\gamma<1$, the space of functions $v$ such that the following seminorm

$$[v]_{C^{0,\gamma}(\Omega)}=\sup_{x,y\in\Omega,x\neq y}\frac{|v(x)-v(y)|}{|x-y|^{\gamma}}$$

is finite.

Given a ball $B_{r}(x_{0})\subset\Omega$, let us define

$${A}_{B_{r}(x_{0})}(\xi)=\displaystyle\fint_{B_{r}(x_{0})}{A}(x,\xi)dx.$$

One can easily check that if ${A}(x,\xi)$ satisfies (A1) and (A2) they also hold true for the operator ${A}_{B}(\xi)$. Setting

$$\mathcal{H}(x,B_{r}(x_{0})):=\sup_{\xi\neq 0}\dfrac{|{A}(x,\xi)-{A}_{B_{r}(x_{0})}(\xi)|}{(\mu^{2}+|\xi|^{2})^{\frac{1}{2}}},$$ (2.1)

we will say that the map $x\mapsto{A}(x,\xi)$ is locally uniformly in VMO if for each compact set $K\subset\Omega$ we have that

$$\lim_{R\rightarrow 0}\sup_{r<R}\sup_{x_{0}\in K}\displaystyle\fint_{B_{r}(x_{0})}\mathcal{H}(x,B_{r}(x_{0}))dx=0.$$ (2.2)

We have that if ${A}(x,\xi)$ satisfies (A3) then it has the locally uniform VMO property (see [2] for the proof).

Lemma 2.3.

Let $A:\Omega\times\mathbb{R}^{n}\to\mathbb{R}$ be such that (A1), (A2) and (A3) hold. Then, $A$ is locally uniformly in VMO.

3 Difference quotient

We recall some properties of the finite difference quotient operator that will be needed in the sequel.

Definition 3.1.

Let $F$ be a function defined in an open set $\Omega\subset\mathbb{R}^{n}$ and let $h\in\mathbb{R}^{n}$. We call the difference quotient of $F$ with respect to $h$ the function

$$\tau_{h}^{1}F(x)=\tau_{h}F(x):=F(x+h)-F(x).$$

Moreover, for every $r\in\mathbb{N}$ we define the $(r+1)$-th finite difference quotient with respect to $h$ as

$$\tau^{r+1}_{h}F(x):=\tau_{h}(\tau^{r}_{h}F(x)).$$

The function $\tau_{h}F$ is defined in the set

$$\Omega_{|h|}:=\{x\in\Omega:\mathrm{dist}(x,\partial\Omega)>|h|\}.$$

We start with the description of some elementary properties that can be found, for example, in [19].

Proposition 3.2.

Let $F\in W^{1,p}(\Omega)$, with $p\geq 1$, and let $G:\Omega\rightarrow\mathbb{R}$ be a measurable function. Then
(i) $\tau_{h}F\in W^{1,p}(\Omega_{|h|})$ and

$$D_{i}(\tau_{h}F)=\tau_{h}(D_{i}F).$$

(ii) If at least one of the functions $F$ or $G$ has support contained in $\Omega_{|h|}$, then

$$\int_{\Omega}F\tau_{h}Gdx=\displaystyle\int_{\Omega}G\tau_{-h}Fdx.$$

(iii) We have

$$\tau_{h}(FG)(x)=F(x+h)\tau_{h}G(x)+G(x)\tau_{h}F(x).$$

The next result about the finite difference operator is a kind of integral version of Lagrange Theorem.

Lemma 3.3.

If $0<\rho<R,$ $|h|<\frac{R-\rho}{2},$ $1<p<+\infty$ and $F\in W^{1,p}(B_{R})$, then

$\displaystyle\int_{B_{\rho}}|\tau_{h}F(x)|^{p}dx\leq c(n,p)|h|^{p}\displaystyle\int_{B_{R}}|DF(x)|^{p}dx$.

Moreover, it holds

$\displaystyle\int_{B_{\rho}}|F(x+h)|^{p}dx\leq\displaystyle\int_{B_{R}}|F(x)|^{p}dx$.

4 Besov spaces

We give the definition of Besov spaces as done in [42, Section 2.5.12, Theorem 1].

Definition 4.1.

Let $1\leq p<+\infty$ and let $\alpha>0$ be a positive real number. Denote by $r$ the smallest integer larger than $\alpha$.
Let $1\leq q<+\infty$. We say that a function $v:\mathbb{R}^{n}\rightarrow\mathbb{R}$ belongs to the Besov space $B^{\alpha}_{p,q}(\mathbb{R}^{n})$ if, and only if, $v\in L^{p}(\mathbb{R}^{n})$ and

$$[v]_{B^{\alpha}_{p,q}(\mathbb{R}^{n})}:=\biggl(\displaystyle\int_{\mathbb{R}^{n}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{|\tau^{r}_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}<+\infty.$$ (4.1)

Equivalently, we could simply say that $v\in L^{p}(\mathbb{R}^{n})$ and $\frac{\tau^{r}_{h}{v}}{|h|^{\alpha}}\in L^{q}\bigl(\frac{dh}{|h|^{n}};L^{p}(\mathbb{R}^{n})\bigr)$.
When $q=+\infty$, the Besov space $B^{\alpha}_{p,\infty}(\mathbb{R}^{n})$ consists of the functions $v\in L^{p}(\mathbb{R}^{n})$ such that

$$[v]_{B^{\alpha}_{p,\infty}(\mathbb{R}^{n})}:=\displaystyle\sup_{h\in\mathbb{R}^{n}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{|\tau^{r}_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{1}{p}}<+\infty.$$ (4.2)

Accordingly, for $1\leq q\leq+\infty$, the Besov space $B^{\alpha}_{p,q}(\mathbb{R}^{n})$ is normed with

$$\|v\|_{B^{\alpha}_{p,q}(\mathbb{R}^{n})}:=\|v\|_{L^{p}(\mathbb{R}^{n})}+[v]_{B^{\alpha}_{p,q}(\mathbb{R}^{n})}.$$ (4.3)

Observe that, if $1\leq q<+\infty$, integrating for $h\in B(0,\delta)$ for a fixed $\delta>0,$ an equivalent norm is obtained, because

$\biggl(\displaystyle\int_{\{|h|\geq\delta\}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{|\tau^{r}_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}\leq c(n,\alpha,p,q,\delta)\|v\|_{L^{p}(\mathbb{R}^{n})}$.

In the case $q=+\infty$, one can simply take the supremum over $|h|\leq\delta$ and obtain an equivalent norm. By construction, one has $B^{\alpha}_{p,q}(\mathbb{R}^{n})\subset L^{p}(\mathbb{R}^{n})$. One also has the following version of the Sobolev embeddings (a proof can be found at [25, Proposition 7.12]).

Lemma 4.2.

Suppose that $0<\alpha<1$.
(a) If $1<p<\frac{n}{\alpha}$ and $1\leq q\leq p^{*}_{\alpha}=\frac{np}{n-\alpha p}$, then there is a continuous embedding $B^{\alpha}_{p,q}(\mathbb{R}^{n})\subset L^{p^{*}_{\alpha}}(\mathbb{R}^{n})$.
(b) If $p=\frac{n}{\alpha}$ and $1\leq q\leq+\infty$, then there is a continuous embedding $B^{\alpha}_{p,q}(\mathbb{R}^{n})\subset BMO(\mathbb{R}^{n})$,
where $BMO$ denotes the space of functions with bounded mean oscillations [19, Chapter 2].

We recall the following inclusions between Besov spaces ([25, Proposition 7.10 and Formula (7.35)]).

Lemma 4.3.

Suppose that $0<\beta<\alpha<1$.
(a) If $1<p<+\infty$ and $1\leq q\leq t\leq+\infty$, then $B^{\alpha}_{p,q}(\mathbb{R}^{n})\subset B^{\alpha}_{p,t}(\mathbb{R}^{n})$.
(b) If $1<p<+\infty$ and $1\leq q,t\leq+\infty$, then $B^{\alpha}_{p,q}(\mathbb{R}^{n})\subset B^{\beta}_{p,t}(\mathbb{R}^{n})$.
(c) If $1\leq q\leq+\infty$, then $B^{\alpha}_{\frac{n}{\alpha},q}(\mathbb{R}^{n})\subset B^{\beta}_{\frac{n}{\beta},q}(\mathbb{R}^{n})$.

Differences in (4.1) and (4.2) can be replaced by derivatives (see [42, Remark 2, page 113]). Let $1\leq p,q<+\infty$ and

$$\alpha=\sigma+m,\quad 0<\sigma<1\quad\text{and}\quad m\in\mathbb{N}_{0}.$$

Then, the norm

$$\sum_{|k|\leq m}\|D^{k}v\|_{L^{p}(\mathbb{R}^{n})}+\sum_{|k|\leq m}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}D^{k}v(x)|^{p}}{|h|^{\sigma p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$ (4.4)

is equivalent to (4.3) (with the usual modification if $q=+\infty$).

Given a domain $\Omega\subset\mathbb{R}^{n}$, we say that a function $v:\mathbb{R}^{n}\rightarrow\mathbb{R}$ belongs to the local Besov space $B^{\alpha}_{p,q,\text{loc}}$ if $\varphi v\in B^{\alpha}_{p,q}(\mathbb{R}^{n})$ whenever $\varphi\in\mathcal{C}^{\infty}_{0}(\Omega)$. It is worth noticing that one can prove suitable version of Lemmas 4.2 and 4.3, by using local Besov spaces.

The following lemma provides a characterization of local Besov spaces $B^{\alpha}_{p,q,\text{loc}}$, which will be fundamental for our subsequent analysis (for the proof see [2, Lemma 7], in the case $0<\alpha<1$, and [20, Lemma 2.9], in the case $1\leq\alpha<2$).

Lemma 4.4.

Let $1\leq p<+\infty$, $1\leq q\leq+\infty$ and $0<\alpha<2$. A function $v\in L^{p}_{\text{loc}}(\Omega)$ belongs to the local Besov space $B^{\alpha}_{p,q,\text{loc}}$, if, and only if,

$\biggl\|\dfrac{\tau^{r}_{h}v}{|h|^{\alpha}}\biggr\|_{L^{q}\bigl(\frac{dh}{|h|^{n}};L^{p}(B)\bigr)}<+\infty,$

for any ball $B\subset 2B\subset\Omega$ with radius $r_{B}$. Here the measure $\frac{dh}{|h|^{n}}$ is restricted to the ball $B(0,r_{B})$ on the h-space.

Now, we prove the following local version of the Sobolev embedding stated at Lemma 4.2 .

Proposition 4.5.

Suppose that $0<\alpha<1$, $1<p<\frac{n}{\alpha}$ and $1\leq q\leq p^{*}_{\alpha}$. Then, for every function $v\in B^{\alpha}_{p,q,\text{loc}}(\Omega)$ the following local estimate

$$\|v\|_{L^{\frac{np}{n-\alpha p}}(B_{\rho})} \leq c\|v\|_{L^{p}\left(B_{\frac{R+\rho}{2}}\right)}+\dfrac{c}{(R-\rho)^{\alpha}}\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)}$$

$$\qquad\qquad+c\biggl(\displaystyle\int_{|h|\leq\frac{R-\rho}{4}}\biggl(\displaystyle\int_{B_{\frac{3R+\rho}{4}}}\dfrac{|\tau_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$ (4.5)

holds for every balls $B_{\rho}\subset B_{R}\Subset\Omega$ and for a positive constant $c=c(n,p,q,\alpha)$.

Proof.

Let us consider radii $0<\rho<R$ and balls $B_{\rho}\subset B_{R}\Subset\Omega$ and a cut-off function $\eta\in\mathcal{C}_{0}^{\infty}(B_{\frac{R+\rho}{2}})$ such that $0\leq\eta\leq 1$, $\eta=1$ on $B_{\rho}$, $|D\eta|\leq\frac{C}{R-\rho}$. By virtue of Lemma 4.2 $(a)$, we have

$$\|v\|_{L^{\frac{np}{n-\alpha p}}(B_{\rho})}\leq\|\eta v\|_{L^{\frac{np}{n-\alpha p}}(\mathbb{R}^{n})}\leq c\|\eta v\|_{B^{\alpha}_{p,q}(\mathbb{R}^{n})},$$ (4.6)

where $c$ is a positive constant depending on $n$, $p$, $q$ and $\alpha$.

We take care of the norm $\|\eta v\|_{B^{\alpha}_{p,q}(\mathbb{R}^{n})}$ in (4.6). It holds

$$\|\eta v\|_{B^{\alpha}_{p,q}(\mathbb{R}^{n})} =\|\eta v\|_{L^{p}(\mathbb{R}^{n})}+[\eta v]_{B^{\alpha}_{p,q}(\mathbb{R}^{n})}$$

$$\leq\|v\|_{L^{p}\left(B_{\frac{R+\rho}{2}}\right)}+\biggl(\displaystyle\int_{|h|\leq\frac{R-\rho}{4}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{\eta(x+h)^{p}|\tau_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$

$$\qquad\qquad+\biggl(\displaystyle\int_{|h|>\frac{R-\rho}{4}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{\eta(x+h)^{p}|\tau_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$

$$\qquad\qquad+\biggl(\displaystyle\int_{|h|\leq\frac{R-\rho}{4}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}\eta|^{p}|v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$

$$\qquad\qquad+\biggl(\displaystyle\int_{|h|>\frac{R-\rho}{4}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}\eta|^{p}|v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$

$$\leq\|v\|_{L^{p}\left(B_{\frac{R+\rho}{2}}\right)}+A_{1}+A_{2}+A_{3}+A_{4},$$ (4.7)

where we used $(iii)$ of Proposition 3.2.
Using that $\eta\in\mathcal{C}_{0}^{\infty}(B_{\frac{R+\rho}{2}})$ is such that $0\leq\eta\leq 1$ and the fact that $\alpha>0$, we get the following estimates

$$A_{1}\leq\biggl(\displaystyle\int_{|h|\leq\frac{R-\rho}{4}}\biggl(\displaystyle\int_{B_{\frac{3R+\rho}{4}}}\dfrac{|\tau_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}},$$ (4.8)

$$A_{2}\leq 2\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)}\biggl(\underbrace{\displaystyle\int_{|h|>\frac{R-\rho}{4}}\dfrac{dh}{|h|^{n+\alpha q}}}_{I}\biggr)^{\frac{1}{q}}=\dfrac{c(n,q,\alpha)}{(R-\rho)^{\alpha}}\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)},$$ (4.9)

$$A_{4} \leq\biggl(\displaystyle\int_{|h|>\frac{R-\rho}{4}}\biggl(\displaystyle\int_{B_{\frac{3R+\rho}{4}}}\dfrac{(\eta(x+h)+\eta(x))^{p}|v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$

$$\leq\biggl(\displaystyle\int_{|h|>\frac{R-\rho}{4}}\biggl(\displaystyle\int_{B_{\frac{3R+\rho}{4}}}\dfrac{2^{p}|v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$

$$\leq 2\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)}\biggl(\underbrace{\displaystyle\int_{|h|>\frac{R-\rho}{4}}\dfrac{dh}{|h|^{n+\alpha q}}}_{II}\biggr)^{\frac{1}{q}}=\dfrac{c(n,q,\alpha)}{(R-\rho)^{\alpha}}\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)}.$$ (4.10)

The integral $A_{3}$ can be estimated using the property $|D\eta|\leq\frac{C}{R-\rho}$ and the fact that $\alpha<1$ in the following way

$$A_{3} \leq\biggl(\displaystyle\int_{|h|\leq\frac{R-\rho}{4}}\biggl(\displaystyle\int_{\mathbb{R}^{n}}\|D\eta\|_{\infty}^{p}|v(x)|^{p}|h|^{(1-\alpha)p}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}}$$

$$\leq\dfrac{C}{R-\rho}\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)}\biggl(\underbrace{\displaystyle\int_{|h|\leq\frac{R-\rho}{4}}|h|^{(1-\alpha)q}\dfrac{dh}{|h|^{n}}}_{III}\biggr)^{\frac{1}{q}}=\dfrac{c(n,q,\alpha)}{(R-\rho)^{\alpha}}\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)}.$$ (4.11)

Here, the integrals $I$, $II$ and $III$ were computed using polar coordinates.
Now, combining estimates (4.6), (4.7), (4.8), (4.9), (4.10) and (4.11), we infer

$$\|v\|_{L^{\frac{np}{n-\alpha p}}(B_{\rho})} \leq c\|v\|_{L^{p}\left(B_{\frac{R+\rho}{2}}\right)}+\dfrac{c}{(R-\rho)^{\alpha}}\|v\|_{L^{p}\left(B_{\frac{3R+\rho}{4}}\right)}$$

$$\qquad\qquad+c\biggl(\displaystyle\int_{|h|\leq\frac{R-\rho}{4}}\biggl(\displaystyle\int_{B_{\frac{3R+\rho}{4}}}\dfrac{|\tau_{h}v(x)|^{p}}{|h|^{\alpha p}}dx\biggr)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\biggr)^{\frac{1}{q}},$$

for a positive constant $c=c(n,p,q,\alpha)$. ∎

The next result will help us to deal with second order differences in the fractional context.

Proposition 4.6.

Let $1\leq p,q<+\infty$, $0<\alpha<1$, $M>0$ and $0<r<\rho<R$. There exists a positive constant $c=c(n,p,q,\alpha)$ such that whenever $v\in L^{p}(B_{R})$ satisfies

$$\int_{B_{\frac{R-\rho}{8}}}\left(\int_{B_{R}}\dfrac{|\tau_{h}(\tau_{h}v)|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\leq M^{q},$$ (4.12)

then $v\in B^{1+\alpha}_{p,q,\textrm{loc}}(B_{R})$ and

$$\int_{B_{\rho-r}}\left(\int_{B_{r}}\dfrac{|\tau_{h}Dv|^{p}}{|h|^{p\alpha}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\leq \,cM^{q}+c\|v\|^{q}_{L^{p}(B_{R})}+\frac{c}{(R-\rho)^{q(1+\alpha)}}\|v\|_{L^{p}(B_{R})}^{q}$$

$$\qquad+\frac{c}{(R-\rho)^{q\alpha}}\|Dv\|_{L^{p}(B_{R})}^{q}.$$ (4.13)

Proof.

For $r<\rho$ given in the assumptions, let us consider a cut-off function $\eta\in\mathcal{C}_{0}^{\infty}(B_{\frac{R+\rho}{2}})$ such that $0\leq\eta\leq 1$, $\eta=1$ on $B_{\rho}$, $|D\eta|\leq\frac{C}{R-\rho}$ and $|D^{2}\eta|\leq\frac{C}{(R-\rho)^{2}}$.

The equivalence between the norms (4.3) and (4.4) and the properties of $\eta$ give that

$$\int_{B_{\rho-r}}\left(\int_{B_{r}}\dfrac{|\tau_{h}Dv|^{p}}{|h|^{p\alpha}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\leq\int_{\mathbb{R}^{n}}\left(\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}D(\eta v)|^{p}}{|h|^{p\alpha}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\leq c\|v\|^{q}_{L^{p}(B_{R})}+c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}(\tau_{h}(\eta v))|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\qquad+c\int_{|h|>\frac{R-\rho}{8}}\left(\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}(\tau_{h}(\eta v))|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}.$$ (4.14)

The last integral on the right-hand side of (4.14) can be estimated as follows

$$\int_{|h|>\frac{R-\rho}{8}}\left(\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}(\tau_{h}(\eta v))|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\leq 4^{q}\|v\|^{q}_{L^{p}(B_{R})}\displaystyle\int_{|h|>\frac{R-\rho}{8}}\dfrac{dh}{|h|^{n+(1+\alpha)q}}=\dfrac{c(n,q,\alpha)}{(R-\rho)^{(1+\alpha)q}}\|v\|_{L^{p}(B_{R})}^{q}.$$ (4.15)

We observe that, by (iii) of Proposition 3.2, we have that

$$\tau_{h}(\tau_{h}(\eta v))(x)=\eta(x+2h)\tau^{2}_{h}v(x)+2\tau_{h}\eta(x+h)\tau_{h}v(x)+v(x)\tau^{2}_{h}\eta(x).$$ (4.16)

Therefore, by (4.16), the properties of $\eta$ and (4.12), the second term on the right-hand side of (4.14) can be estimated as follows

$$\int_{B_{\frac{R-\rho}{8}}}\left(\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}(\tau_{h}(\eta v))|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\leq c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{\mathbb{R}^{n}}|\eta(x+2h)|^{p}\frac{|\tau^{2}_{h}v|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\qquad+c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{\mathbb{R}^{n}}|\tau_{h}\eta(x+h)|^{p}\frac{|\tau_{h}v|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\qquad+c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{\mathbb{R}^{n}}|v|^{p}\frac{|\tau^{2}_{h}\eta(x)|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\leq c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{B_{\frac{3R+\rho}{4}}}\frac{|\tau^{2}_{h}v|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\qquad+c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{B_{\frac{3R+\rho}{4}}}\|D\eta\|^{p}_{\infty}\frac{|\tau_{h}v|^{p}}{|h|^{p\alpha}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\qquad+c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{B_{\frac{3R+\rho}{4}}}|h|^{p(1-\alpha)}|v|^{p}\|D^{2}\eta\|^{p}_{\infty}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\leq cM^{q}+cI+cII.$$ (4.17)

Exploiting Lemma 3.3 and the fact $|D\eta|\leq\frac{c}{R-\rho}$, we estimate the integral $I$ as follows

$$I\leq \,c\int_{B_{\frac{R-\rho}{8}}}\left(\int_{B_{\frac{3R+\rho}{4}}}\|D\eta\|^{p}_{\infty}|Dv|^{p}|h|^{p(1-\alpha)}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\leq \,\dfrac{c}{(R-\rho)^{q}}\|Dv\|_{L^{p}(B_{R})}^{q}\int_{B_{\frac{R-\rho}{8}}}|h|^{q(1-\alpha)}\dfrac{dh}{|h|^{n}}=\dfrac{c}{(R-\rho)^{q\alpha}}\|Dv\|_{L^{p}(B_{R})}^{q},$$ (4.18)

where the last integral was computed using polar coordinates.
By using the fact $|D^{2}\eta|\leq\frac{c}{(R-\rho)^{2}}$, we get the following estimate for $II$

$$II\leq \,\dfrac{c}{(R-\rho)^{2q}}\|v\|_{L^{p}(B_{R})}^{q}\int_{B_{\frac{R-\rho}{8}}}|h|^{q(1-\alpha)}\dfrac{dh}{|h|^{n}}=\dfrac{c}{(R-\rho)^{q(1+\alpha)}}\|v\|_{L^{p}(B_{R})}^{q}.$$ (4.19)

Putting (4.18) and (4.19) in (4.17), we find

$$\int_{B_{\frac{R-\rho}{8}}}\left(\int_{\mathbb{R}^{n}}\dfrac{|\tau_{h}(\tau_{h}(\eta v))|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\qquad\leq cM^{q}+\dfrac{c}{(R-\rho)^{q\alpha}}\|Dv\|_{L^{p}(B_{R})}^{q}+\dfrac{c}{(R-\rho)^{q(1+\alpha)}}\|v\|_{L^{p}(B_{R})}^{q}.$$ (4.20)

Eventually, inserting (4.15), (4.20) in (4.14), we obtain the desired estimate (4.13).

∎

Combining Propositions 4.5 and 4.6, we obtain the following

Corollary 4.7.

Suppose that $0<\alpha<1$, $1<p<\frac{n}{\alpha}$ and $1\leq q\leq p^{*}_{\alpha}$. Let $M>0$ and $0<s<t<\overline{R}$. There exists a positive constant $c=c(n,p,q,\alpha)$ such that whenever $u\in L^{p}(B_{\overline{R}})$ satisfies

$$\int_{B_{\frac{\overline{R}-t}{8}}}\left(\int_{B_{\overline{R}}}\dfrac{|\tau_{h}(\tau_{h}u)|^{p}}{|h|^{p(1+\alpha)}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\leq M^{q},$$

then $Du\in L^{\frac{np}{n-\alpha p}}_{\textrm{loc}}(B_{\overline{R}})$ and

$$\|Du\|_{L^{\frac{np}{n-\alpha p}}(B_{s})} \leq cM+c\|u\|_{L^{p}(B_{\overline{R}})}+\frac{c}{(\overline{R}-t)^{1+\alpha}}\|u\|_{L^{p}(B_{\overline{R}})}+\frac{c}{(\overline{R}-t)^{\alpha}}\|Du\|_{L^{p}(B_{\overline{R}})}$$

$$\qquad+c\|Du\|_{L^{p}(B_{\overline{R}})}+\dfrac{c}{(t-s)^{\alpha}}\|Du\|_{L^{p}(B_{\overline{R}})}.$$

Proof.

Let us consider radii $0<s<t<\overline{R}$ and balls $B_{s}\subset B_{t}\subset B_{\overline{R}}$. We can use Proposition 4.6 with $R=\overline{R}$, $\rho=t$ and $r=\frac{3t+s}{4}$, thus getting by (4.13) the following inequality

$$\int_{B_{\frac{t-s}{4}}}\left(\int_{B_{\frac{3t+s}{4}}}\dfrac{|\tau_{h}Du|^{p}}{|h|^{p\alpha}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}} \leq\int_{B_{\frac{t-s}{4}}}\left(\int_{B_{t}}\dfrac{|\tau_{h}Du|^{p}}{|h|^{p\alpha}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}$$

$$\leq cM^{q}+c\|u\|^{q}_{L^{p}(B_{\overline{R}})}+\frac{c}{(\overline{R}-t)^{q(1+\alpha)}}\|u\|_{L^{p}(B_{\overline{R}})}^{q}$$

$$\qquad+\frac{c}{(\overline{R}-t)^{q\alpha}}\|Du\|_{L^{p}(B_{\overline{R}})}^{q}.$$

On the other hand, estimate (4.5) with $v=Du$, $R=t$ and $\rho=s$ reads as

$$\|Du\|_{L^{\frac{np}{n-\alpha p}}(B_{s})} \leq c\|Du\|_{L^{p}\left(B_{\frac{t+s}{2}}\right)}+\dfrac{c}{(t-s)^{\alpha}}\|Du\|_{L^{p}\left(B_{\frac{3t+s}{4}}\right)}$$

$$\qquad\qquad+c\left(\displaystyle\int_{B_{\frac{t-s}{4}}}\left(\displaystyle\int_{B_{\frac{3t+s}{4}}}\dfrac{|\tau_{h}Du(x)|^{p}}{|h|^{\alpha p}}dx\right)^{\frac{q}{p}}\dfrac{dh}{|h|^{n}}\right)^{\frac{1}{q}}.$$

The conclusion follows combining these two inequalities. ∎

We conclude recalling that Besov spaces of fractional order $\alpha\in(0,1)$ can be characterized in pointwise terms. We give the following definition according to [32].

Definition 4.8.

Given a measurable function $v:\mathbb{R}^{n}\rightarrow\mathbb{R}$, a fractional $\alpha$-Hajlasz gradient for $v$ is a sequence $\{g_{k}\}_{k}$ of measurable, non-negative functions $g_{k}:\mathbb{R}^{n}\rightarrow\mathbb{R}$, together with a null set $N\subset\mathbb{R}^{n}$, such that the inequality

$|v(x)-v(y)|\leq(g_{k}(x)+g_{k}(y))|x-y|^{\alpha}$

holds whenever $k\in\mathbb{Z}$ and $x,y\in\mathbb{R}^{n}\setminus N$ are such that $2^{-k}\leq|x-y|<2^{-k+1}$. We say that $\{g_{k}\}_{k}\in l^{q}(\mathbb{Z};L^{p}(\mathbb{R}^{n}))$ if

$\|\{g_{k}\}_{k}\|_{l^{q}(L^{p})}=\biggl(\displaystyle\sum_{k\in\mathbb{Z}}\|g_{k}\|^{q}_{L^{p}(\mathbb{R}^{n})}\biggr)^{\frac{1}{q}}<+\infty.$

The following result was proved in [32].

Theorem 4.9.

Let $0<\alpha<1,$ $1\leq p<+\infty$ and $1\leq q\leq+\infty$. Let $v\in L^{p}(\mathbb{R}^{n})$. Then $v\in B^{\alpha}_{p,q}(\mathbb{R}^{n})$ if, and only if, there exists a fractional $\alpha$-Hajlasz gradient $\{g_{k}\}_{k}\in l^{q}(\mathbb{Z};L^{p}(\mathbb{R}^{n}))$ for $v$. Moreover,

$\|v\|_{B^{\alpha}_{p,q}(\mathbb{R}^{n})}\simeq\inf\|\{g_{k}\}_{k}\|_{l^{q}(L^{p})},$

where the infimum runs over all possible fractional $\alpha$-Hajlasz gradients for $v$.

5 A priori estimate

We recall that by assumption $f\in B^{\alpha}_{\frac{n}{2\alpha},p,\textrm{loc}}(\Omega)$, by the very definition of Besov space, we know that $f\in L^{\frac{n}{2\alpha}}_{\textrm{loc}}(\Omega)$. By standard elliptic regularity (see for example [12, Chapter 10]), under assumptions (A1) and (A2) and $f\in L^{\frac{n}{2\alpha}}_{\textrm{loc}}(\Omega)$, local weak solutions $u$ of (1.1) are locally bounded in $\Omega$ and satisfy the following $L^{\infty}$-bound

$$\|u\|_{L^{\infty}(B_{\rho}(x_{0}))}\leq\,cR^{1-\alpha}\|f\|_{L^{\frac{n}{2\alpha}}(B_{R}(x_{0}))}+\dfrac{c}{(R-\rho)^{n/(1-\alpha)}}\left(\fint_{B_{R}(x_{0})}u^{2}dx\right)^{\frac{1}{2}}$$ (5.1)

for every radii $0<\rho<R$, where $B_{R}(x_{0})\Subset\Omega$ and $c$ is a positive constant depending on $n$, $\alpha$, $\nu$ and $L$.

The next Theorem is devoted to the proof of an a priori estimate, which is the main tool in the proof of Theorem 1.2.

Theorem 5.1.

Let $\delta\geq 1$, $0<\alpha<1$ and $1\leq p\leq\frac{2n}{n-2\alpha}$. Let $u\in W^{1,2}_{\textrm{loc}}(\Omega)$ be a local weak solution of (1.1) such that $Du\in L^{\frac{n(\delta+1)}{n-2\alpha}}_{\textrm{loc}}(\Omega)$. Then, for every $x_{0}\in\Omega$ there exists a radius $R_{\delta}=R_{\delta}(x_{0},n,\delta,\alpha,\nu,L,p)\leq 1$ and a ball $B_{R_{\delta}}(x_{0})\Subset\Omega$ such that we have the following estimate

$$\left(\int_{B_{r}(x_{0})}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{n-2\alpha}{n}} \leq\frac{\tilde{c}}{(R-r)^{2(1+\alpha)}}\Biggl[\int_{B_{R}(x_{0})}|u|^{\delta+1}dx$$

$$\qquad+\int_{B_{R}(x_{0})}|Du|^{\delta+1}dx+[f]^{\delta+1}_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+1\Biggr],$$ (5.2)

for every radii $0<r<R\leq R_{\delta}$ and for a positive constant $\tilde{c}=\tilde{c}(n,\delta,\alpha,\nu,L,p)$.

Proof.

Fix a ball $B_{R_{\delta}}(x_{0})\Subset\Omega$ and consider radii $0<r<\sigma<t^{\prime}<R\leq R_{\delta}$ (the radius $R_{\delta}$ will be choosen later). We divide the interval $(\sigma,t^{\prime})$ into five equal parts by means of $\rho^{\prime}$, $\rho$, $s$ and $t$, so that $r<\sigma<\rho^{\prime}<\rho<s<t<t^{\prime}<R\leq R_{\delta}$, with

$$\rho^{\prime}-\sigma=\rho-\rho^{\prime}=s-\rho=t-s=t^{\prime}-t=\frac{t^{\prime}-\sigma}{5}.$$ (5.3)

For the sake of notation, from now on we shall omit the dependence of the balls on the center $x_{0}$.

Let us consider a cut-off function $\eta\in C_{0}^{\infty}(B_{t})$, with $\eta=1$ on $B_{s}$, $0\leq\eta\leq 1$, $|D\eta|\leq\frac{c}{t-s}$ and $|h|\leq\frac{t^{\prime}-\sigma}{40}$. We test the equation (1.1) with the function

$$\varphi=\tau_{-h}(\eta^{2}|\tau_{h}u|^{\delta-1}\tau_{h}u)\qquad(\delta\geq 1)$$

and then we use (ii) of Proposition 3.2, thus getting

$$I=\int_{\Omega}\langle\,\tau_{h}\left(A(x,Du)\right),D(\eta^{2}|\tau_{h}u|^{\delta-1}\tau_{h}u)\,\rangle=\int_{\Omega}\eta^{2}\,\tau_{h}f\,|\tau_{h}u|^{\delta-1}\tau_{h}u=I_{5}.$$ (5.4)

The integral $I$ can be written as follows

$$I =\int_{\Omega}\langle A(x+h,Du(x+h))-A(x+h,Du(x)),D(\eta^{2}|\tau_{h}u|^{\delta-1}\tau_{h}u)\rangle dx$$

$$\qquad+\int_{\Omega}\langle A(x+h,Du(x))-A(x,Du(x)),D(\eta^{2}|\tau_{h}u|^{\delta-1}\tau_{h}u)\rangle dx$$

$$=\delta\int_{\Omega}\langle A(x+h,Du(x+h))-A(x+h,Du(x)),\eta^{2}|\tau_{h}u|^{\delta-1}\tau_{h}Du\,\rangle dx$$

$$\qquad+2\int_{\Omega}\langle A(x+h,Du(x+h))-A(x+h,Du(x)),\eta D\eta|\tau_{h}u|^{\delta-1}\tau_{h}u\,\rangle dx$$

$$\qquad+\delta\int_{\Omega}\langle A(x+h,Du(x))-A(x,Du(x)),\eta^{2}|\tau_{h}u|^{\delta-1}\tau_{h}Du\,\rangle dx$$

$$\qquad+2\int_{\Omega}\langle A(x+h,Du(x))-A(x,Du(x)),\eta D\eta|\tau_{h}u|^{\delta-1}\tau_{h}u\,\rangle dx$$

$$=I_{1}+I_{2}+I_{3}+I_{4}.$$ (5.5)

From equalities (5.4) and (5.5), we deduce that

$$I_{1}\leq|I_{2}|+|I_{3}|+|I_{4}|+|I_{5}|.$$ (5.6)

By virtue of assumption (A1), we infer

$$I_{1}\geq\delta\nu\int_{\Omega}\eta^{2}|\tau_{h}Du|^{2}\,|\tau_{h}u|^{\delta-1}dx.$$ (5.7)

Now we consider the term $|I_{2}|$. By using hypothesis (A2), Young’s inequality, the properties of $\eta$ we have

$$I_{2} \leq 2L\int_{\Omega}\eta\,|\tau_{h}Du||D\eta||\tau_{h}u|^{\delta}$$

$$\leq\delta\frac{\nu}{4}\int_{\Omega}\eta^{2}|\tau_{h}Du|^{2}|\tau_{h}u|^{\delta-1}dx+\frac{c(\nu,L,n,\delta)}{(t-s)^{2}}\int_{B_{t}}|\tau_{h}u|^{\delta+1}dx$$

$$\leq\delta\frac{\nu}{4}\int_{\Omega}\eta^{2}|\tau_{h}Du|^{2}|\tau_{h}u|^{\delta-1}dx+\frac{c(\nu,L,n,\delta)}{(t-s)^{2}}|h|^{\delta+1}\int_{B_{t^{\prime}}}|Du|^{\delta+1}dx,$$ (5.8)

where in the last inequality we used Lemma 3.3. From condition (A3), Young’s inequality and the properties of $\eta$, we derive for $2^{-k-1}\frac{t^{\prime}-\sigma}{40}\leq|h|<2^{-k}\frac{t^{\prime}-\sigma}{40}$, $k\in\mathbb{N},$

$$I_{3} \leq\delta|h|^{\alpha}\int_{\Omega}\eta^{2}\,\left(g_{k}(x+h)+g_{k}(x)\right)|\tau_{h}u|^{\delta-1}\,|\tau_{h}Du|\,(\mu^{2}+|Du|^{2})^{\frac{1}{2}}\,dx$$

$$\leq\delta\frac{\nu}{4}\int_{\Omega}\eta^{2}\,|\tau_{h}u|^{\delta-1}\,|\tau_{h}Du|^{2}\,dx$$

$$\qquad+c(\nu)\delta|h|^{2\alpha}\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{2}|\tau_{h}u|^{\delta-1}\,(1+|Du|^{2})\,dx$$

$$=:\delta\frac{\nu}{4}\int_{\Omega}\eta^{2}\,|\tau_{h}u|^{\delta-1}\,|\tau_{h}Du|^{2}\,dx+c(\nu)\delta|h|^{2\alpha}I^{\prime}_{3}.$$ (5.9)

Using Hölder’s Inequality and Lemma 3.3, we get

$$I^{\prime}_{3} \leq\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}\,\left(\int_{B_{t}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}\frac{2}{\delta+1}}\,\left(\int_{B_{t}}|\tau_{h}u|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{n-2\alpha}{n}\frac{\delta-1}{\delta+1}}$$

$$\leq c(n,\delta,\alpha)|h|^{\delta-1}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}\,\left(\int_{B_{t}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}\frac{2}{\delta+1}}\,$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{n-2\alpha}{n}\frac{\delta-1}{\delta+1}}$$

$$\leq c(n,\delta,\alpha)|h|^{\delta-1}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}\,\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}}.$$ (5.10)

Hence, it follows that

$$I_{3} \leq\delta\frac{\nu}{4}\int_{\Omega}\eta^{2}\,|\tau_{h}u|^{\delta-1}\,|\tau_{h}Du|^{2}\,dx$$

$$\qquad+c(n,\nu,\delta,\alpha)|h|^{\delta-1+2\alpha}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}}.$$ (5.11)

From (A3), Young’s inequality, the properties of $\eta$ and Lemma 3.3, we get for $2^{-k-1}\frac{t^{\prime}-\sigma}{40}\leq|h|<2^{-k}\frac{t^{\prime}-\sigma}{40}$, $k\in\mathbb{N}$

$$I_{4} \leq 2|h|^{\alpha}\int_{\Omega}\eta\,\left(g_{k}(x+h)+g_{k}(x)\right)|D\eta|\,|\tau_{h}u|^{\delta}\,(\mu^{2}+|Du|^{2})^{\frac{1}{2}}\,dx$$

$$\leq c\int_{B_{t}}|D\eta|^{2}\,|\tau_{h}u|^{\delta+1}\,dx$$

$$\qquad+c\int_{B_{t}}|h|^{2\alpha}\left(g_{k}(x+h)+g_{k}(x)\right)^{2}|\tau_{h}u|^{\delta-1}\,(1+|Du|^{2})\,dx$$

$$\leq\frac{c}{(t-s)^{2}}\int_{B_{t}}\,|\tau_{h}u|^{\delta+1}\,dx+c|h|^{2\alpha}I^{\prime}_{3}$$

$$\leq\frac{c}{(t-s)^{2}}|h|^{\delta+1}\int_{B_{t}}\,|Du|^{\delta+1}\,dx+c|h|^{2\alpha}I^{\prime}_{3},$$ (5.12)

where $I^{\prime}_{3}$ has been defined at (5).
Inserting (5) in (5), we deduce

$$I_{4} \leq\frac{c}{(t-s)^{2}}|h|^{\delta+1}\int_{B_{t}}\,|Du|^{\delta+1}\,dx$$

$$\qquad+c(n,\delta,\alpha)|h|^{\delta-1+2\alpha}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}\,\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}}.$$ (5.13)

Using Hölder’s inequality with exponents $\left(\frac{n(\delta+1)}{\delta(n-2\alpha)},\frac{n(\delta+1)}{n+2\alpha\delta}\right)$, the properties of $\eta$ and Lemma 3.3, we have

$$I_{5} \leq\int_{\Omega}\eta^{2}\,|\tau_{h}f|\,|\tau_{h}u|^{\delta}dx$$

$$\leq\left(\int_{B_{t}}|\tau_{h}f|^{\frac{n(\delta+1)}{n+2\alpha\delta}}dx\right)^{\frac{n+2\alpha\delta}{n(\delta+1)}}\,\left(\int_{B_{t}}|\tau_{h}u|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}}$$

$$\leq|h|^{\delta}c(n,\delta,\alpha)R^{\frac{n-2\alpha}{\delta+1}}\left(\int_{B_{t}}|\tau_{h}f|^{\frac{n}{2\alpha}}dx\right)^{\frac{2\alpha}{n}}\,\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}}$$

$$\leq c(n,\delta,\alpha)R^{\frac{n-2\alpha}{\delta+1}}|h|^{\delta+\alpha}\left(\int_{B_{t}}\frac{|\tau_{h}f|^{\frac{n}{2\alpha}}}{|h|^{\frac{n}{2}}}dx\right)^{\frac{2\alpha}{n}}\,\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}},$$ (5.14)

where we also used that $\frac{n(\delta+1)}{n+2\alpha\delta}\leq\frac{n}{2\alpha}$ for every $\delta\geq 1$.
Inserting (5.7), (5), (5), (5) and (5) in (5.6), we get

$$\delta\nu\int_{\Omega}\eta^{2}|\tau_{h}Du|^{2}\,|\tau_{h}u|^{\delta-1}dx \leq\delta\frac{\nu}{2}\int_{\Omega}\eta^{2}|\tau_{h}Du|^{2}|\tau_{h}u|^{\delta-1}dx+\frac{c(\nu,L,n,\delta)}{(t-s)^{2}}|h|^{\delta+1}\int_{B_{t^{\prime}}}|Du|^{\delta+1}dx$$

$$\qquad+c(n,\nu,\delta,\alpha)|h|^{\delta-1+2\alpha}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}}$$

$$\qquad+c(n,\delta,\alpha)R^{\frac{n-2\alpha}{\delta+1}}|h|^{\delta+\alpha}\left(\int_{B_{t}}\frac{|\tau_{h}f|^{\frac{n}{2\alpha}}}{|h|^{\frac{n}{2}}}dx\right)^{\frac{2\alpha}{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}}.$$ (5.15)

Reabsorbing the first integral on the right-hand side of (5) by the left-hand side and using that $\eta=1$ on $B_{s}$, we obtain

$$\int_{B_{s}}|\tau_{h}Du|^{2}\,|\tau_{h}u|^{\delta-1}dx \leq\frac{c(\nu,L,n,\delta)}{(t-s)^{2}}|h|^{\delta+1}\int_{B_{t^{\prime}}}|Du|^{\delta+1}dx$$

$$\qquad+c(n,\nu,\delta,\alpha)|h|^{\delta-1+2\alpha}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}}$$

$$\qquad+c(n,\nu,\delta,\alpha)R^{\frac{n-2\alpha}{\delta+1}}|h|^{\delta+\alpha}\left(\int_{B_{t}}\frac{|\tau_{h}f|^{\frac{n}{2\alpha}}}{|h|^{\frac{n}{2}}}dx\right)^{\frac{2\alpha}{n}}\,\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}}.$$ (5.16)

Exploiting the identity

$$\left|D\left(|\tau_{h}u|^{\frac{\delta-1}{2}}\,\tau_{h}u\right)\right|^{2}=\left(\frac{\delta+1}{2}\right)^{2}|\tau_{h}Du|^{2}\,|\tau_{h}u|^{\delta-1}$$

and Lemma 3.3, it follows that

$$\int_{B_{\rho}}\left|\tau_{\lambda}\left(|\tau_{h}u|^{\frac{\delta-1}{2}}\,\tau_{h}u\right)\right|^{2}dx\leq \,c\,|\lambda|^{2}\int_{B_{s}}\left|D\left(|\tau_{h}u|^{\frac{\delta-1}{2}}\,\tau_{h}u\right)\right|^{2}dx$$

$$= \,c\,|\lambda|^{2}(\delta+1)^{2}\int_{B_{s}}|\tau_{h}Du|^{2}\,|\tau_{h}u|^{\delta-1}.$$ (5.17)

Choosing $\lambda=h$, which is legitimate by our choice of the step size $|h|<\frac{t^{\prime}-\sigma}{40}$ and putting (5.16) in (5.17), we derive

$$\int_{B_{\rho}}\left|\tau_{h}\left(|\tau_{h}u|^{\frac{\delta-1}{2}}\,\tau_{h}u\right)\right|^{2}dx \leq\dfrac{c(\nu,L,n,\delta)}{(t-s)^{2}}|h|^{\delta+3}\int_{B_{t^{\prime}}}|Du|^{\delta+1}dx$$

$$\qquad+{c(n,\nu,\delta,\alpha)}|h|^{\delta+1+2\alpha}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}}$$

$$\qquad+{c(n,\nu,\delta,\alpha)}R^{\frac{n-2\alpha}{\delta+1}}|h|^{\delta+\alpha+2}\left(\int_{B_{t}}\frac{|\tau_{h}f|^{\frac{n}{2\alpha}}}{|h|^{\frac{n}{2}}}dx\right)^{\frac{2\alpha}{n}}\,$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}}.$$ (5.18)

An application of Lemma 2.1 with $p=\delta+1$ and the fact that $\delta+1\geq 2$ give that

$$|\tau_{h}(\tau_{h}u)|^{\delta+1}\leq\left|\tau_{h}\left(|\tau_{h}u|^{\frac{\delta-1}{2}}\,\tau_{h}u\right)\right|^{2}.$$ (5.19)

Combining (5.18) and (5.19) and then dividing both sides by $|h|^{\delta+1+2\alpha}$, we get

$$\int_{B_{\rho}}\frac{|\tau_{h}(\tau_{h}u)|^{\delta+1}}{|h|^{\delta+1+2\alpha}}dx \leq\dfrac{c(\nu,L,n,\delta)}{(t-s)^{2}}|h|^{2(1-\alpha)}\int_{B_{t^{\prime}}}|Du|^{\delta+1}dx$$

$$\qquad+{c(n,\nu,\delta,\alpha)}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{2\alpha}{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{n-2\alpha}{n}}$$

$$\qquad+{c(n,\nu,\delta,\alpha)}R^{\frac{n-2\alpha}{\delta+1}}|h|^{1-\alpha}\left(\int_{B_{t}}\frac{|\tau_{h}f|^{\frac{n}{2\alpha}}}{|h|^{\frac{n}{2}}}dx\right)^{\frac{2\alpha}{n}}\,$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}}.$$ (5.20)

Now, raise to the power $\frac{p}{2}$, with $p$ appearing in (1.2), integrate over the ball $B(0,(t^{\prime}-\sigma)/40)$ w.r.t. the measure $\frac{dh}{|h|^{n}}.$ Since the functions $g_{k}$ are defined for $2^{-k-1}\frac{t^{\prime}-\sigma}{40}\leq|h|<2^{-k}\frac{t^{\prime}-\sigma}{40}$, we interpret the ball $B(0,(t^{\prime}-\sigma)/40)$ as

$$B(0,(t^{\prime}-\sigma)/40)=\bigcup_{k=0}^{\infty}B(0,2^{-k}(t^{\prime}-\sigma)/40)\setminus B(0,2^{-k-1}(t^{\prime}-\sigma)/40)=:\bigcup_{k=0}^{\infty}E_{k}.$$

Then, by (5.20), we infer the following estimate

$$\int_{B_{\frac{\rho-\rho^{\prime}}{8}}(0)}\left(\int_{B_{\rho}}\frac{|\tau_{h}(\tau_{h}u)|^{\delta+1}}{|h|^{\delta+2\alpha+1}}dx\right)^{\frac{p}{2}}\frac{dh}{|h|^{n}} \leq\frac{c}{(t-s)^{p}}\int_{B_{\frac{t^{\prime}-\sigma}{40}}(0)}\left(\int_{B_{t^{\prime}}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}|h|^{{(1-\alpha)p}}\frac{dh}{|h|^{n}}$$

$$\qquad+c\sum_{k=0}^{\infty}\int_{E_{k}}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{p\alpha}{n}}\frac{dh}{|h|^{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{(n-2\alpha)p}{2n}}$$

$$\qquad+{c}R^{\frac{p(n-2\alpha)}{2(\delta+1)}}\int_{B_{\frac{t^{\prime}-\sigma}{40}}(0)}\left(\int_{B_{t}}\frac{|\tau_{h}f|^{\frac{n}{2\alpha}}}{|h|^{\frac{n}{2}}}dx\right)^{\frac{\alpha p}{n}}|h|^{\frac{(1-\alpha)p}{2}}\frac{dh}{|h|^{n}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}{\frac{p}{2}}},$$ (5.21)

with $c=c(n,\delta,\alpha,\nu,L,p)$ positive constant. Here, in the integral on the left-hand side, we used the fact that $B(0,(t^{\prime}-\sigma)/40)=B(0,(\rho-\rho^{\prime})/8)$, which follows from the definition of radii given in (5.3).

Now, we take care of the term

$$J=\sum_{k=0}^{\infty}\int_{E_{k}}\left(\int_{B_{t}}\left(g_{k}(x+h)+g_{k}(x)\right)^{\frac{n}{\alpha}}dx\right)^{\frac{p\alpha}{n}}\frac{dh}{|h|^{n}}.$$

Following the arguments in [2], we write the right-hand side of the previous estimate in polar coordinates, so $h\in E_{k}$ if, and only if, $h=m\xi$ for some $2^{-k-1}\frac{t^{\prime}-\sigma}{40}\leq m<2^{-k}\frac{t^{\prime}-\sigma}{40}$ and some $\xi$ in the unit sphere $\mathbb{S}^{n-1}$ on $\mathbb{R}^{n}$. We denote by $dS(\xi)$ the surface measure on $\mathbb{S}^{n-1}$. Setting $\Delta_{m\xi}g_{k}(x)=g_{k}(x+m\xi)$ and $m_{k}=2^{-k}(t^{\prime}-\sigma)/40$, we infer

$$J\leq \displaystyle\sum_{k=0}^{\infty}\displaystyle\int_{m_{k+1}}^{m_{k}}\displaystyle\int_{\mathbb{S}^{n-1}}\biggl(\displaystyle\int_{B_{t}}(g_{k}(x+m\xi)+g_{k}(x))^{\frac{n}{\alpha}}dx\biggr)^{\frac{p\alpha}{n}}dS(\xi)\dfrac{dm}{m}$$

$$= \displaystyle\sum_{k=0}^{\infty}\displaystyle\int_{m_{k+1}}^{m_{k}}\displaystyle\int_{\mathbb{S}^{n-1}}\|\Delta_{m\xi}g_{k}+g_{k}\|_{L^{\frac{n}{\alpha}}(B_{t})}^{p}dS(\xi)\dfrac{dm}{m}.$$

We note that for each $\xi\in\mathbb{S}^{n-1}$ and $m_{k-1}\leq m\leq m_{k}$

$$\|\Delta_{m\xi}g_{k}+g_{k}\|_{L^{\frac{n}{\alpha}}(B_{t})}\leq\, \|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{t-m_{k}\xi})}+\|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{t})}$$

$$\leq \,2\|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{R})},$$

hence

$$J\leq c(n)\,2^{p}\,\text{log }2\,\sum_{k=0}^{\infty}\|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{R})}^{p},$$ (5.22)

which is finite by assumption (1.2).

Inserting (5.22) in (5.21), we find that

$$\int_{B_{\frac{\rho-\rho^{\prime}}{8}}(0)}\left(\int_{B_{\rho}}\frac{|\tau_{h}(\tau_{h}u)|^{\delta+1}}{|h|^{\delta+2\alpha+1}}dx\right)^{\frac{p}{2}}\frac{dh}{|h|^{n}} \leq\frac{c}{(t-s)^{p}}\left(\int_{B_{t^{\prime}}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}\int_{B_{R}(0)}|h|^{{(1-\alpha)p}}\frac{dh}{|h|^{n}}$$

$$\qquad+c\sum_{k=0}^{\infty}\|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{R})}^{p}\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{(n-2\alpha)p}{2n}}$$

$$\qquad+cR^{\frac{p(n-2\alpha)}{2(\delta+1)}}\,[f]^{\frac{p}{2}}_{B^{\alpha}_{\frac{n}{2\alpha},p}(B_{R})}\left(\int_{B_{R}(0)}|h|^{{(1-\alpha)p}}\frac{dh}{|h|^{n}}\right)^{\frac{1}{2}}$$

$$\qquad\cdot\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{(n-2\alpha)\delta}{n(\delta+1)}{\frac{p}{2}}}.$$ (5.23)

We observe that, since $\alpha\in(0,1)$, the integral

$$\int_{B_{R}(0)}|h|^{{(1-\alpha)p}}\frac{dh}{|h|^{n}}=c(n,\alpha,p,R)$$

is finite.

Now, since $p\in[1,\frac{2n}{n-2\alpha}]$, we can use Corollary 4.7 with $\delta+1$, $\frac{p(\delta+1)}{2}$ and $\frac{2\alpha}{\delta+1}$ in place of $p$, $q$ and $\alpha$ respectively. This, together with (5.23) and (5.3), gives for every $r<\sigma<t^{\prime}<R$ the following estimate

$$\left(\int_{B_{\sigma}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{p(n-2\alpha)}{2n}} \leq c\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}+\frac{c}{(t^{\prime}-\sigma)^{p(1+\alpha)}}\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+c\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{c}{(t^{\prime}-\sigma)^{{\alpha p}}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{c}{(t^{\prime}-\sigma)^{p}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\gamma\sum_{k=0}^{\infty}\|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{R})}^{p}\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{(n-2\alpha)p}{2n}}$$

$$\qquad+cR^{\frac{p(n-2\alpha)}{2(\delta+1)}}\,[f]^{\frac{p}{2}}_{B^{\alpha}_{\frac{n}{2\alpha},p}(B_{R})}\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}\,dx\right)^{\frac{(n-2\alpha)\delta p}{2n(\delta+1)}},$$ (5.24)

for a positive constants $c=c(n,\delta,\alpha,\nu,L,p)$ and $\gamma=\gamma(n,\delta,\alpha,\nu,L,p)$.

Then, applying Young’s inequality in the last term on the right-hand side of (5.24), we derive

$$\left(\int_{B_{\sigma}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{p(n-2\alpha)}{2n}} \leq c\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}+\frac{c}{(t^{\prime}-\sigma)^{p(1+\alpha)}}\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+c\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{{c}}{(t^{\prime}-\sigma)^{\alpha p}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{{c}}{(t^{\prime}-\sigma)^{p}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\gamma\sum_{k=0}^{\infty}\|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{R})}^{p}\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}\,dx\right)^{\frac{(n-2\alpha)p}{2n}}$$

$$\qquad+c\gamma R^{\frac{p(n-2\alpha)}{2}}\sum_{k=0}^{\infty}\|g_{k}\|_{L^{\frac{n}{\alpha}}(B_{R})}^{p}$$

$$\qquad+\frac{1}{4}\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}\,dx\right)^{\frac{(n-2\alpha)p}{2n}}+{c}R^{\frac{p(n-2\alpha)}{2}}\,[f]^{\frac{p(\delta+1)}{2}}_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R})}.$$

where $c$ and $\gamma$ are positive constants depending at most on $(n,\delta,\alpha,\nu,L,p).$

We choose a radius $0<R_{\delta}\leq 1$, depending on $(x_{0},n,\delta,\alpha,\nu,L,p),$ such that for $0<R\leq R_{\delta}$ it holds

$$\gamma\sum_{k=0}^{\infty}\|g_{k}\|^{p}_{L^{\frac{n}{\alpha}}(B_{R})}<\dfrac{1}{4}.$$ (5.25)

We show that such a radius $R_{\delta}$ does exist. Indeed, by assumption (1.2), we have that there exists $m\in\mathbb{N}$ such that

$$\sum_{k=m+1}^{\infty}\|g_{k}\|^{p}_{L^{\frac{n}{\alpha}}(B_{R})}<\dfrac{1}{8\gamma}.$$

Moreover, by the absolute continuity of the integral, for every $k\in\{0,\dots,m\}$ there exists $R_{\delta,k}$ such that if $R\leq R_{\delta,k}$, then

$$\|g_{k}\|^{p}_{L^{\frac{n}{\alpha}}(B_{R})}<\dfrac{1}{8(m+1)\gamma}.$$

Hence, (5.25) follows by choosing $R_{\delta}=\min\{R_{\delta_{k}},\,k=0,\dots,m\}$.

Therefore, if $0<R\leq R_{\delta}$, we get for every $r<\sigma<t^{\prime}<R$

$$\left(\int_{B_{\sigma}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{p(n-2\alpha)}{2n}} \leq c\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}+\frac{c}{(t^{\prime}-\sigma)^{p(1+\alpha)}}\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+c\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{{c}}{(t^{\prime}-\sigma)^{\alpha p}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{{c}}{(t^{\prime}-\sigma)^{p}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{1}{4}\left(\int_{B_{t^{\prime}}}(1+|Du|^{\frac{n(\delta+1)}{n-2\alpha}})\,dx\right)^{\frac{(n-2\alpha)p}{2n}}$$

$$\qquad+cR^{\frac{p(n-2\alpha)}{2}}$$

$$\qquad+\frac{1}{4}\left(\int_{B_{t^{\prime}}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}\,dx\right)^{\frac{(n-2\alpha)p}{2n}}+cR^{\frac{p(n-2\alpha)}{2}}\,[f]^{\frac{p(\delta+1)}{2}}_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R})}.$$ (5.26)

Thus applying the iteration Lemma 2.2, we infer

$$\left(\int_{B_{r}}|Du|^{\frac{n(\delta+1)}{n-2\alpha}}dx\right)^{\frac{p(n-2\alpha)}{2n}} \leq\tilde{c}\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}+\frac{\tilde{c}}{(R-r)^{p(1+\alpha)}}\left(\int_{B_{R}}|u|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\tilde{c}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{\tilde{c}}{(R-r)^{\alpha p}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\frac{\tilde{c}}{(R-r)^{p}}\left(\int_{B_{R}}|Du|^{\delta+1}dx\right)^{{\frac{p}{2}}}$$

$$\qquad+\tilde{c}R^{\frac{p(n-2\alpha)}{2}}+{\tilde{c}}R^{\frac{p(n-2\alpha)}{2}}\,[f]^{\frac{p(\delta+1)}{2}}_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R})},$$

for every concentric balls $B_{r}\subset B_{R}\subseteq B_{R_{\delta}}\Subset\Omega$, where $\tilde{c}=\tilde{c}(n,\delta,\alpha,\nu,L,p)$ is a positive constant. Eventually, raising to the power $2/p$, we get the desired estimate. ∎

Theorem 5.1 allows us to implement a Moser-type iteration scheme that improves the regularity of a local weak solution to (1.1) from $W^{1,2}_{\textrm{loc}}(\Omega)$ to $W^{1,q}_{\textrm{loc}}(\Omega)$, for every $q\in[2,+\infty)$.

Theorem 5.2.

Let $0<\alpha<1$ and $1\leq p\leq\frac{2n}{n-2\alpha}$. Let $u\in W^{1,2}_{\textrm{loc}}(\Omega)$ be a local weak solution of (1.1) such that $Du\in L^{\frac{2n}{n-2\alpha}}_{\textrm{loc}}(\Omega)$. Then, $u\in W^{1,q}_{\textrm{loc}}(\Omega)$ for every $q\in[2,+\infty)$ and for every $x_{0}\in\Omega$ there exists a radius $\tilde{R}=\tilde{R}(x_{0},n,q,\alpha,\nu,L,p)\leq 1$ and a ball $B_{\tilde{R}}(x_{0})\Subset\Omega$ such that we have the following estimate

$$\left(\int_{B_{R/2}(x_{0})}|Du|^{q}dx\right)^{\frac{1}{q}} \leq\tilde{c}\Biggl\{\|u\|_{L^{\infty}(B_{R}(x_{0}))}+[f]_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+\left(\int_{B_{R}(x_{0})}|Du|^{2}dx\right)^{\frac{1}{2}}+1\Biggr\}.$$ (5.27)

for every radius $0<R\leq\tilde{R}$ and for a positive constant $\tilde{c}=\tilde{c}(n,q,\alpha,\nu,L,p,R)$.

Proof.

We define a sequence of exponents $(\delta_{i})_{i\in\mathbb{N}_{0}}$ by

$$\begin{cases}\delta_{0}=1,\\ \delta_{i+1}=\dfrac{n(\delta_{i}+1)}{n-2\alpha}-1=2\left(\dfrac{n}{n-2\alpha}\right)^{i+1}-1,\qquad i\geq 0.\end{cases}$$

Clearly $\delta_{i}\to+\infty$ as $i\to+\infty$. Hence, for all $q\in[2,+\infty)$ there exists $i_{*}\in\mathbb{N}$ such that $\delta_{i_{*}-1}+1<q\leq\delta_{i_{*}}+1$.

For any $i\in\{0,1,\dots,i_{*}\}$, let $0<R_{\delta_{i}}\leq 1$ be the radius such that (5.25) holds. Now, fix a radius $R$ such that $R<\min\{R_{\delta_{i}}:i=0,1,\dots,i_{*}\}$ and define a sequence of radii $(\rho_{i})_{i\in\mathbb{N}_{0}}$ by

$$\rho_{i}=\dfrac{R}{2}\left(1+\dfrac{1}{2^{i}}\right).$$

By the a priori assumption $Du\in L^{\frac{2n}{n-2\alpha}}_{\mathrm{loc}}(\Omega)$, we are legitimate to use estimate (5.2) with $\delta=\delta_{i-1}$, $r=\rho_{i}$ and $R=\rho_{i-1}$, thus getting

$$\left(\int_{B_{\rho_{i}}(x_{0})}|Du|^{\delta_{i}+1}dx\right)^{\frac{\delta_{i-1}+1}{\delta_{i}+1}} \leq\frac{\tilde{c_{i}}}{(\rho_{i-1}-\rho_{i})^{2(1+\alpha)}}\Biggl[\int_{B_{\rho_{i-1}}(x_{0})}|u|^{\delta_{i-1}+1}dx$$

$$\qquad+\int_{B_{\rho_{i-1}}(x_{0})}|Du|^{\delta_{i-1}+1}dx+[f]^{\delta_{i-1}+1}_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+1\Biggr].$$ (5.28)

Since $\rho_{i-1}-\rho_{i}=\frac{R}{2^{i+1}}$, we get

$$\left(\int_{B_{\rho_{i}}(x_{0})}|Du|^{\delta_{i}+1}dx\right)^{\frac{\delta_{i-1}+1}{\delta_{i}+1}} \leq\frac{\tilde{c_{i}}2^{(i+1)(1+\alpha)}}{R^{2(1+\alpha)}}\Biggl[\int_{B_{\rho_{i-1}}(x_{0})}|u|^{\delta_{i-1}+1}dx$$

$$\qquad+\int_{B_{\rho_{i-1}}(x_{0})}|Du|^{\delta_{i-1}+1}dx+[f]^{\delta_{i-1}+1}_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+1\Biggr].$$ (5.29)

Define

$$K_{i}=\max\{1,\tilde{c_{i}}2^{2(i+1)(1+\alpha)}\}.$$

Then, we have

$$\left(\int_{B_{\rho_{i}}(x_{0})}|Du|^{\delta_{i}+1}dx\right)^{\frac{\delta_{i-1}+1}{\delta_{i}+1}} \leq\frac{K_{i}}{R^{2(1+\alpha)}}\Biggl\{\int_{B_{\rho_{i-1}}(x_{0})}|u|^{\delta_{i-1}+1}dx$$

$$\qquad+\int_{B_{\rho_{i-1}}(x_{0})}|Du|^{\delta_{i-1}+1}dx+\,[f]^{\delta_{i-1}+1}_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+1\Biggr\}.$$ (5.30)

Raising to the power $\frac{1}{\delta_{i-1}+1}$, we find that

$$\left(\int_{B_{\rho_{i}}(x_{0})}|Du|^{\delta_{i}+1}dx\right)^{\frac{1}{\delta_{i}+1}} \leq\frac{K_{i}^{\frac{1}{\delta_{i-1}+1}}2^{\frac{2}{\delta_{i-1}+1}}}{R^{\frac{2(1+\alpha)}{\delta_{i-1}+1}}}\Biggl\{\left(\int_{B_{\rho_{i-1}}(x_{0})}|u|^{\delta_{i-1}+1}dx\right)^{\frac{1}{\delta_{i-1}+1}}$$

$$\qquad+\left(\int_{B_{\rho_{i-1}}(x_{0})}|Du|^{\delta_{i-1}+1}dx\right)^{\frac{1}{\delta_{i-1}+1}}+\,[f]_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+1\Biggr\}.$$ (5.31)

Using the boundedness of $u$ and the fact $\rho_{i-1}\leq R\leq 1,$ we infer

$$\left(\int_{B_{\rho_{i}}(x_{0})}|Du|^{\delta_{i}+1}dx\right)^{\frac{1}{\delta_{i}+1}} \leq\frac{K_{i}^{\frac{1}{\delta_{i-1}+1}}2^{\frac{2}{\delta_{i-1}+1}}}{R^{\frac{2(1+\alpha)}{\delta_{i-1}+1}}}\Biggl\{c(n)^{\frac{1}{\delta_{i-1}+1}}\|u\|_{L^{\infty}(B_{R}(x_{0}))}$$

$$\qquad+\left(\int_{B_{\rho_{i-1}}(x_{0})}|Du|^{\delta_{i-1}+1}dx\right)^{\frac{1}{\delta_{i-1}+1}}+\,[f]_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+1\Biggr\}.$$ (5.32)

Setting

$$\overline{K_{i}}=K_{i}^{\frac{1}{\delta_{i-1}+1}}2^{\frac{2}{\delta_{i-1}+1}}c(n)^{\frac{2}{\delta_{i-1}+1}},\qquad\Gamma=\|u\|_{L^{\infty}(B_{R}(x_{0}))}+[f]_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R}(x_{0}))}+1,$$

the previous estimate reads as

$$\left(\int_{B_{\rho_{i}}(x_{0})}|Du|^{\delta_{i}+1}dx\right)^{\frac{1}{\delta_{i}+1}} \leq\frac{\overline{K}_{i}}{R^{\frac{2(1+\alpha)}{\delta_{i-1}+1}}}\Biggl\{\Gamma+\left(\int_{B_{\rho_{i-1}}(x_{0})}|Du|^{\delta_{i-1}+1}dx\right)^{\frac{1}{\delta_{i-1}+1}}\Biggr\}.$$ (5.33)

Iterating this estimate for $i\in\{1,\dots,i_{*}\}$, where $i_{*}$ has been chosen at the beginning of the proof, we get

$$\left(\int_{B_{\rho_{i_{*}}}(x_{0})}|Du|^{\delta_{i_{*}}+1}dx\right)^{\frac{1}{\delta_{i_{*}}+1}} \leq\tilde{K}_{R}\Biggl\{\Gamma+\left(\int_{B_{R}(x_{0})}|Du|^{2}dx\right)^{\frac{1}{2}}\Biggr\},$$ (5.34)

where $\tilde{K_{R}}=i_{*}\prod_{i=1}^{i_{*}}\frac{\overline{K_{i}}}{R^{\frac{2(1+\alpha)}{\delta_{i-1}+1}}}$. Thanks to Hölder’s inequality, we have

$$\left(\int_{B_{R/2}(x_{0})}|Du|^{q}dx\right)^{\frac{1}{q}} \leq c(n,q,\delta_{i_{*}},R)\tilde{K}_{R}\Biggl\{\Gamma+\left(\int_{B_{R}(x_{0})}|Du|^{2}dx\right)^{\frac{1}{2}}\Biggr\},$$ (5.35)

that is the conclusion. ∎

6 Proof of Theorem 1.2

In this section, we give the proof of Theorem 1.2. First, we define a sequence of regular problems with a unique smooth solution $u_{\varepsilon}$, satisfying the a priori estimates. Then, we show that the solutions $u_{\varepsilon}$ converge to the original solution of our problem, which subsequently inherits the same regularity properties.

Proof of Theorem 1.2.

Let us fix a non-negative smooth kernel $\phi\in\mathcal{C}^{\infty}_{0}(B_{1}(0))$ such that $\int_{B_{1}(0)}\phi=1$ and consider the corresponding family of mollifiers $(\phi_{\varepsilon})_{\varepsilon>0}$. For a given ball $B_{2R}\Subset\Omega$ and every $0<\varepsilon<\text{dist}(B_{R},\partial\Omega)$, we set

$$A_{\varepsilon}(x,\xi)=\int_{B_{\varepsilon}(0)}\phi_{\varepsilon}(y)A(x-y,\xi)dy,\quad f_{\varepsilon}(x)=f*\phi_{\varepsilon}(x),$$

for a.e. $x\in B_{R}$ and every $\xi\in\mathbb{R}^{n}$. One can easily check that $A_{\varepsilon}(x,\xi)$ is such that $A_{\varepsilon}(\cdot,0)\in L^{2}(B_{R})$ and satisfies assumptions (A1) and (A2). Moreover, denoting by $g^{\varepsilon}_{k}=g_{k}*\phi_{\varepsilon}$, from (1.2) and (A3) it follows that

$$\sum_{k=0}^{\infty}\|g^{\varepsilon}_{k}\|^{p}_{L^{\frac{n}{\alpha}}(B_{R})}<\infty$$ (6.1)

and

$$|A_{\varepsilon}(x,\xi)-A_{\varepsilon}(y,\xi)|\leq|x-y|^{\alpha}(g^{\varepsilon}_{k}(x)+g^{\varepsilon}_{k}(y))(\mu^{2}+|\xi|^{2})^{\frac{1}{2}},$$ (A4)

for a.e. $x,y\in B_{R}$ such that $2^{-k}R\leq|x-y|<2^{-k+1}R$ and for every $\xi\in\mathbb{R}^{n}$.

Let $u\in W^{1,2}_{\textrm{loc}}(\Omega)$ be a local weak solution to (1.1) and let $u_{\varepsilon}\in u+W^{1,2}_{0}(B_{R})$ the solution to the Dirichlet problem

$$\begin{cases}-\text{ div}A_{\varepsilon}(x,Du_{\varepsilon})=f_{\varepsilon}(x)&\qquad\text{in }B_{R},\\ u_{\varepsilon}=u&\qquad\text{on }\partial B_{R}.\end{cases}$$ (6.2)

By standard elliptic regularity (see for example [19, Chapter 8]), we know that $u_{\varepsilon}\in\mathcal{C}^{1,\beta}_{\textrm{loc}}(B_{R})$, for some $\beta\in(0,1)$.

Since $u_{\varepsilon}\in u+W_{0}^{1,2}(B_{R})$ is a solution to (6.2) and $u\in W^{1,2}(\Omega)$ is a solution to (1.1), choosing $\varphi=u_{\varepsilon}-u$ as test function for both equations, we obtain

$$\int_{B_{R}}\langle A_{\varepsilon}(x,Du_{\varepsilon}),D(u_{\varepsilon}-u)\rangle dx=\int_{B_{R}}f_{\varepsilon}(u_{\varepsilon}-u)dx,$$ (6.3)

and

$$\int_{B_{R}}\langle A(x,Du),D(u_{\varepsilon}-u)\rangle dx=\int_{B_{R}}f(u_{\varepsilon}-u)dx.$$ (6.4)

Subtracting the previous two equalities, we derive

$$I=\int_{B_{R}}\langle A_{\varepsilon}(x,Du_{\varepsilon})-A(x,Du)\,,\,D(u_{\varepsilon}-u)\rangle dx=\int_{B_{R}}\left(f_{\varepsilon}-f\right)(u_{\varepsilon}-u)dx.$$ (6.5)

We can write the integral on the left-hand side of (6.5) as follows

$$I=\int_{B_{R}}\langle A_{\varepsilon}(x,Du_{\varepsilon})-A_{\varepsilon}(x,Du)\,,\,Du_{\varepsilon}-Du\rangle dx+\int_{B_{R}}\langle A_{\varepsilon}(x,Du)-A(x,Du)\,,\,Du_{\varepsilon}-Du\rangle dx,$$ (6.6)

thus getting

$$\int_{B_{R}}\langle A_{\varepsilon}(x,Du_{\varepsilon})-A_{\varepsilon}(x,Du)\,,\,Du_{\varepsilon}-Du\rangle dx= \int_{B_{R}}\left(f_{\varepsilon}-f\right)(u_{\varepsilon}-u)dx$$

$$- \int_{B_{R}}\langle A_{\varepsilon}(x,Du)-A(x,Du)\,,\,Du_{\varepsilon}-Du\rangle dx.$$ (6.7)

Then, by assumption (A1) and Hölder’s inequality, we infer

$$\nu\int_{B_{R}}\lvert Du_{\varepsilon}-Du\rvert^{2}dx \leq\int_{B_{R}}\langle A_{\varepsilon}(x,Du_{\varepsilon})-A_{\varepsilon}(x,Du)\,,\,Du_{\varepsilon}-Du\rangle dx$$

$$\leq\int_{B_{R}}\lvert f_{\varepsilon}-f\rvert\,\lvert u_{\varepsilon}-u\rvert dx$$

$$\qquad+\int_{B_{R}}\lvert A_{\varepsilon}(x,Du)-A(x,Du)\rvert\,\lvert Du_{\varepsilon}-Du\rvert dx$$

$$\leq\left(\int_{B_{R}}\lvert f_{\varepsilon}-f\rvert^{\frac{2n}{n+2}}dx\right)^{\frac{n+2}{2n}}\left(\int_{B_{R}}\lvert u_{\varepsilon}-u\rvert^{\frac{2n}{n-2}}dx\right)^{\frac{n-2}{2n}}$$

$$\qquad+\left(\int_{B_{R}}\lvert A_{\varepsilon}(x,Du)-A(x,Du)\rvert^{2}dx\right)^{\frac{1}{2}}\left(\int_{B_{R}}\lvert Du_{\varepsilon}-Du\rvert^{2}dx\right)^{\frac{1}{2}}$$

$$\leq c\left(\int_{B_{R}}\lvert Du_{\varepsilon}-Du\rvert^{2}dx\right)^{\frac{1}{2}}\Biggl\{\left(\int_{B_{R}}\lvert f_{\varepsilon}-f\rvert^{\frac{2n}{n+2}}dx\right)^{\frac{n+2}{2n}}$$

$$\qquad+\left(\int_{B_{R}}\lvert A_{\varepsilon}(x,Du)-A(x,Du)\rvert^{2}dx\right)^{\frac{1}{2}}\Biggr\},$$

where in the last line we used Sobolev-Poincaré inequality. Next, we divide both sides by

$$\left(\int_{B_{R}}\lvert Du_{\varepsilon}-Du\rvert^{2}dx\right)^{\frac{1}{2}}$$

and square the resulting inequality. This yields

$$\int_{B_{R}}\lvert Du_{\varepsilon}-Du\rvert^{2}dx \leq c\left(\int_{B_{R}}\lvert f_{\varepsilon}-f\rvert^{\frac{2n}{n+2}}dx\right)^{\frac{n+2}{n}}+c\int_{B_{R}}\lvert A_{\varepsilon}(x,Du)-A(x,Du)\rvert^{2}dx$$

$$\leq c\left(\int_{B_{R}}\lvert f_{\varepsilon}-f\rvert^{\frac{n}{2\alpha}}dx\right)^{\frac{\alpha}{n}}+c\int_{B_{R}}\lvert A_{\varepsilon}(x,Du)-A(x,Du)\rvert^{2}dx,$$ (6.8)

where we used that $\frac{2n}{n+2}<\frac{n}{\alpha}$. By assumption (A2) and the fact that $A_{\varepsilon}(x,0)\in L^{2}(B_{R})$, we have

$$\lvert A_{\varepsilon}(x,Du)\rvert \leq\lvert A_{\varepsilon}(x,Du)-A_{\varepsilon}(x,0)\rvert+\lvert A_{\varepsilon}(x,0)\rvert$$

$$\leq L|Du|+|A_{\varepsilon}(x,0)|\in L^{2}(B_{R}).$$

Since

$$A_{\varepsilon}(x,0)\rightarrow A(x,0)\quad\text{in }L^{2}(B_{R})\text{ as }\varepsilon\to 0$$

and

$$A_{\varepsilon}(x,Du)\rightarrow A(x,Du)\quad\text{a.e.\ in }B_{R}\text{ as }\varepsilon\to 0,$$

by Lebesgue dominated convergence Theorem, we get

$$A_{\varepsilon}(x,Du)\rightarrow A(x,Du)\quad\text{in }L^{2}(B_{R})\text{ as }\varepsilon\to 0.$$

Therefore, letting $\varepsilon\to 0$ in (6.8), we have that

$$u_{\varepsilon}\to u\quad\text{strongly in }W^{1,2}_{0}(B_{R}).$$

Now, $u_{\varepsilon}$ satisfies inequalities (5.1) and (5.27), i.e.

$$\|u_{\varepsilon}\|_{L^{\infty}(B_{R})}\leq\,cR^{1-\alpha}\|f_{\varepsilon}\|_{L^{\frac{n}{2\alpha}}(B_{2R})}+\dfrac{c}{R^{n/(1-\alpha)}}\left(\fint_{B_{2R}}u_{\varepsilon}^{2}dx\right)^{\frac{1}{2}}$$

and

$$\left(\int_{B_{R/2}}|Du_{\varepsilon}|^{q}dx\right)^{\frac{1}{q}} \leq\tilde{c}\Biggl\{\|u_{\varepsilon}\|_{L^{\infty}(B_{R})}+[f_{\varepsilon}]_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{R})}+\left(\int_{B_{R}}|Du_{\varepsilon}|^{2}dx\right)^{\frac{1}{2}}+1\Biggr\}.$$

Thus, we have

$$\left(\int_{B_{R/2}}|Du_{\varepsilon}|^{q}dx\right)^{\frac{1}{q}} \leq\tilde{c}\left(\|f_{\varepsilon}\|_{B_{\frac{n}{2\alpha},p}^{\alpha}(B_{2R})}+\|u_{\varepsilon}\|_{W^{1,2}(B_{2R})}+1\right),$$

where $\tilde{c}$ is a positive constant independent of $\varepsilon$. Eventually, letting $\varepsilon\to 0$ and using a covering argument, we conclude that $u\in W^{1,q}_{\textrm{loc}}(\Omega)$. ∎

7 Example

In this section, we exhibit an example showing that Theorem 1.2 is sharp, in the sense that under assumption (A3), one cannot, in general, expect solutions to (1.1) to be locally Lipschitz continuous. More precisely, we provide an equation satisfying assumption (A3) whose solution satisfies

$$Du\in L^{q}_{\mathrm{loc}}\quad\text{for every }q<\infty,\qquad\text{but}\qquad Du\notin L^{\infty}_{\mathrm{loc}}.$$

We recall the following

Proposition 7.1.

Setting for every $r\in(0,1)$ and a fixed $r_{0}>1$

$$v(r)=\log\frac{r_{0}}{r},$$

the function $u(x)=x_{1}v(|x|)$ solves the equation

$$-\ \mathrm{div}(A(x)Du)=0\qquad\text{in }B_{1}(0),$$

where the matrix $A(x)=(a_{ij}(x))_{\begin{subarray}{c}1\leq i\leq n\\ 1\leq j\leq n\end{subarray}}$ is defined by

$$a_{ij}(x)=\delta_{ij}+\gamma(|x|)\left(\delta_{ij}-\frac{x_{i}x_{j}}{|x|^{2}}\right),\qquad\gamma(r)=\frac{-n}{(n-1)\log\frac{r_{0}}{r}}.$$ (7.1)

Moreover, $u$ is such that $Du\in L^{q}_{\mathrm{loc}}(B_{1}(0))$ for every $q>1$, but $Du\notin L^{\infty}_{\mathrm{loc}}(B_{1}(0))$.

See [29, Proposition 1.5] for the proof.

In order to prove that the operator $A(x,\xi)=A(x)\xi$, with $A(x)$ defined in (7.1), satisfies the assumption (A3), it is sufficient to show that

$$A(x)\in W^{1,n}_{\textrm{loc}}(B_{1}(0)),$$ (7.2)

since by [42, Remark 3, Section 2.7.1] and Lemma 4.3, it holds

$$W^{1,n}_{\textrm{loc}}(B_{1}(0))\subset B^{\beta}_{\frac{n}{\beta},n,\textrm{loc}}(B_{1}(0))\subset B^{\alpha}_{\frac{n}{\alpha},p,\textrm{loc}}(B_{1}(0)),$$

for every $0<\alpha<\beta<1$.
By virtue of the pointwise characterization of Sobolev spaces in [24], (7.2) is equivalent to saying that

$$\lvert A(x)-A(y)\rvert\leq|x-y|\left(h(x)+h(y)\right),$$

for some non-negative function $h\in L^{n}_{\mathrm{loc}}(B_{1}(0))$ and for all $x,y\in B_{1}(0)$.

From the definition of $a_{ij}(x)$, we have

$$\lvert A(x)-A(y)\rvert =\sum_{i,j}\lvert a_{ij}(x)-a_{ij}(y)\rvert$$

$$=\sum_{i,j}\left\lvert\left(\gamma(|x|)-\gamma(|y|)\right)\delta_{ij}-\gamma(|x|)\frac{x_{i}x_{j}}{|x|^{2}}+\gamma(|y|)\frac{y_{i}y_{j}}{|y|^{2}}\right\rvert$$

$$\leq n\left\lvert\gamma(|x|)-\gamma(|y|)\right\rvert+\sum_{i,j}\left\lvert\gamma(|x|)\frac{x_{i}x_{j}}{|x|^{2}}-\gamma(|y|)\frac{y_{i}y_{j}}{|y|^{2}}\right\rvert$$

$$=:nI_{1}+I_{2}.$$ (7.3)

Without loss of generality, we assume $|x|>|y|$. From (7.1), we obtain

$$I_{1} =\frac{n}{n-1}\left\lvert-\frac{1}{\log\frac{r_{0}}{|x|}}+\frac{1}{\log\frac{r_{0}}{|y|}}\right\rvert$$

$$=\frac{n}{n-1}\left\lvert\frac{\log|y|-\log|x|}{\log\frac{r_{0}}{|x|}\log\frac{r_{0}}{|y|}}\right\rvert.$$

Since

$$\lvert\log|y|-\log|x|\rvert =\int_{|y|}^{|x|}\frac{1}{t}\,dt\leq\int_{|y|}^{|x|}\frac{1}{|y|}\,dt=\frac{1}{|y|}(|x|-|y|)$$

$$\leq(|x|-|y|)\left(\frac{1}{|y|}+\frac{1}{|x|}\right)$$

$$\leq|x-y|\left(\frac{1}{|y|}+\frac{1}{|x|}\right),$$

where we used the fact that $\frac{1}{t}\leq\frac{1}{|y|}$. So, we have

$$I_{1}\leq \,\frac{n}{n-1}|x-y|\left(\frac{1}{|x|}+\frac{1}{|y|}\right)\frac{1}{\log\frac{r_{0}}{|x|}\log\frac{r_{0}}{|y|}}$$

$$\leq \,\frac{n}{n-1}|x-y|\left(\frac{1}{|x|\,\log\frac{r_{0}}{|x|}\log\frac{r_{0}}{|y|}}+\frac{1}{|y|\,\log\frac{r_{0}}{|x|}\log\frac{r_{0}}{|y|}}\right)$$

$$\leq \,c(n,r_{0})|x-y|\left(\frac{1}{|x|\,\log\frac{r_{0}}{|x|}}+\frac{1}{|y|\,\log\frac{r_{0}}{|y|}}\right),$$ (7.4)

where we used that $\frac{1}{\log\frac{r_{0}}{|x|}}$ is bounded in $B_{1}(0)$.
We take care of $I_{2}$

$$I_{2} =\left\lvert\gamma(|x|)\frac{x_{i}x_{j}}{|x|^{2}}-\gamma(|y|)\frac{x_{i}x_{j}}{|x|^{2}}+\gamma(|y|)\frac{x_{i}x_{j}}{|x|^{2}}-\gamma(|y|)\frac{y_{i}y_{j}}{|y|^{2}}\right\rvert$$

$$\leq\left\lvert\gamma(|x|)-\gamma(|y|)\right\rvert\frac{|x_{i}x_{j}|}{|x|^{2}}+|\gamma(|y|)|\left\lvert\frac{x_{i}x_{j}}{|x|^{2}}-\frac{y_{i}y_{j}}{|y|^{2}}\right\rvert$$

$$=:I^{\prime}_{2}+I^{\prime\prime}_{2}.$$ (7.5)

The term $I^{\prime}_{2}$ can be estimated as follows

$$I_{2}^{\prime}\leq\left\lvert\gamma(|x|)-\gamma(|y|)\right\rvert\left(\frac{x_{i}^{2}}{2}+\frac{x_{j}^{2}}{2}\right)\frac{1}{|x|^{2}}\leq\left\lvert\gamma(|x|)-\gamma(|y|)\right\rvert.$$ (7.6)

On the other hand, for $I_{2}^{{}^{\prime\prime}}$ we have

$$I_{2}^{\prime\prime} =|\gamma(|y|)|\,\frac{\left|x_{i}x_{j}|y|^{2}-y_{i}y_{j}|y|^{2}+y_{i}y_{j}|y|^{2}-y_{i}y_{j}|x|^{2}\right|}{|x|^{2}|y|^{2}}$$

$$\leq|\gamma(|y|)|\,\frac{\left\lvert x_{i}x_{j}-y_{i}y_{j}\right\rvert|y|^{2}+|y_{i}y_{j}|\left(|x|^{2}-|y|^{2}\right)}{|x|^{2}|y|^{2}}$$

$$\leq|\gamma(|y|)|\,\left(\frac{\left\lvert x_{i}x_{j}-y_{i}y_{j}\right\rvert}{|x|^{2}}+\frac{|y_{i}y_{j}|\left(|x|^{2}-|y|^{2}\right)}{|x|^{2}|y|^{2}}\right)$$

$$\leq|\gamma(|y|)|\,\left(\frac{\left\lvert x_{i}x_{j}-x_{i}y_{j}+x_{i}y_{j}-y_{i}y_{j}\right\rvert}{|x|^{2}}+\left(\frac{y_{i}^{2}}{2}+\frac{y_{j}^{2}}{2}\right)\frac{\left(|y|-|x|\right)\left(|y|+|x|\right)}{|x|^{2}|y|^{2}}\right)$$

$$\leq|\gamma(|y|)|\,\left(\frac{|x_{j}||x_{i}-y_{j}|+|y_{j}||x_{i}-y_{i}|}{|x|^{2}}+\frac{1}{2}\frac{\left(|y|-|x|\right)\left(|y|+|x|\right)}{|x|^{2}}\right)$$

$$\leq|\gamma(|y|)|\,\left(\frac{|x-y|(|x|+|y|)}{|x|^{2}}+\frac{1}{2}|y-x|\left(\frac{|y|}{|x|^{2}}+\frac{1}{|x|}\right)\right)$$

$$\leq\frac{c(n)}{\log\frac{r_{0}}{|y|}}\,|y-x|\left(\frac{|y|}{|x|^{2}}+\frac{1}{|x|}\right)$$

$$\leq c(n)|y-x|\left(\frac{|y|}{|x|^{2}\,\log\frac{r_{0}}{|y|}}+\frac{1}{|x|\,\log\frac{r_{0}}{|y|}}\right)$$

$$\leq c(n)|y-x|\left(\frac{1}{|x|\,\log\frac{r_{0}}{|x|}}+\frac{1}{|y|\,\log\frac{r_{0}}{|y|}}\right),$$ (7.7)

where we used that $|x|>|y|$. Inserting (7.6) and (7) in (7), we infer

$$I_{2} \leq\left\lvert\gamma(|x|)-\gamma(|y|)\right\rvert+c(n)|y-x|\left(\frac{1}{|x|\,\log\frac{r_{0}}{|x|}}+\frac{1}{|y|\,\log\frac{r_{0}}{|y|}}\right)$$

$$\leq c(n)|y-x|\left(\frac{1}{|x|\,\log\frac{r_{0}}{|x|}}+\frac{1}{|y|\,\log\frac{r_{0}}{|y|}}\right),$$ (7.8)

where the first term of r.h.s. has been estimated as (7). Combining (7) and (7), we find

$$\lvert A(x)-A(y)\rvert \leq c(n,r_{0})|x-y|\left(\frac{1}{|x|\,\log\frac{r_{0}}{|x|}}+\frac{1}{|y|\,\log\frac{r_{0}}{|y|}}\right).$$ (7.9)

In order to obtain assumption (A3), it is sufficient to choose

$$h(x):=\frac{1}{|x|\log\frac{r_{0}}{|x|}}.$$

It remains to check whether $\frac{1}{|x|\log\frac{r_{0}}{|x|}}\in L^{n}(B_{1}(0))$. We have that

$$\int_{B_{1}(0)}\left(\frac{1}{|x|\log\frac{r_{0}}{|x|}}\right)^{n}dx =c(n)\int_{0}^{1}\frac{1}{\rho\log^{n}\frac{r_{0}}{\rho}}\,d\rho<+\infty.$$

Acknowledgments. The authors are members of the Gruppo Nazionale per l’Analisi Matematica, la Probabilità e le loro Applicazioni (GNAMPA) of the Istituto Nazionale di Alta Matematica (INdAM). The authors have been partially supported through the INdAM - GNAMPA 2026 Project “Esistenza e regolarità per soluzioni di equazioni ellittiche e paraboliche anisotrope” (CUP: E53C25002010001). In addition S. Russo has also been supported through the project: Sustainable Mobility Center (Centro Nazionale per la Mobilità Sostenibile – CNMS) - SPOKE 10, grant number E63C22000930007.