Higher Sobolev regularity on the mixed local and nonlocal p-Laplace equations

We develop a systematic study of the interior Sobolev regularity of weak solutions to the mixed local and nonlocal $p$-Laplace equations. To be precise, we show that the weak solution $u$ belongs to $W^{2, p}_\mathrm{loc}$ and even $W^{2, 2}_{\rm loc}$ Sobolev spaces in the subquadratic case, while $|\nabla u|^{\frac{p-2}{2}}\nabla u$ is of the class $W^{1, 2}_\mathrm{loc}$ in the superquadratic scenario, both of which coincide with that of the classical $p$-Laplace equations. Moreover, an improved higher fractional differentiability and integrability result $u\in W^{1+\beta, q}_\mathrm{loc}$ is proved in the full range $p\in (1, \infty)$ for any $q\in [\max\{p, 2\}, \infty)$ and $\beta\in(0, \frac 2q)$. The main analytical tools are the finite difference quotient technique, suitable energy method and tail estimates. As far as we know, our results are new within the context of such mixed problems.