Gradient Hölder regularity in mixed local and nonlocal linear parabolic problem

We prove the local Hölder regularity of weak solutions to the mixed local nonlocal parabolic equation of the form \begin{equation*} u_t-Δu+$\text{P.V.}$\int_{$\mathbb{R}$^{n}} {\frac{u(x,t)-u(y,t)}{{\left|x-y\right|}^{n+2s}}}dy=0, \end{equation*} where $0<s<1$; for every exponent ${\alpha }_0\in (0,1)$. Here, $\mathrm{\Delta }$ is the usual Laplace operator. Next, we show that the gradients of weak solutions are also $\alpha$-Hölder continuous for some $\alpha \in (0,1)$. Our approach is purely analytic and it is based on perturbation techniques.