Consider the parabolic equation in $p: \mathbb R^2\to\mathbb R$


$$\partial\_t p + b(t)\partial\_x p + D(t,x)\partial^2\_{xx}p=0,$$


where $b$, $D$ are nice enough functions. I look for the continuity of the derivatives $\partial\_t p$, $\partial\_x p$ of the solution. It is known by Nash's paper ([Continuity of Solutions of Parabolic and Elliptic Equations](https://www.jstor.org/stable/2372841?seq=1#metadata_info_tab_contents)) that, under very reasonable conditions on $b$, $D$, we have the Hölder-type continuity of $p$. Is there any work concerning such continuity analysis of $\partial\_t p$, $\partial\_x p$?


PS: My idea is to consider $q\mathrel{:=}\partial\_x p$. Then


$$\partial\_t q + b(t)\partial\_xq + D(t,x)\partial^2\_{xx}q=-\partial\_x D(t,x)\partial\_{xx}p$$


is a similar parabolic equation for $q$ with an additional source. But the term $\partial\_x D(t,x)\partial\_{xx}p$ contains $\partial\_{xx}p$, which makes the estimation even harder….



Reference request: continuity of the derivatives of the (fundamental) solution to a parabolic equation
The equation for $q = p\_x$ can be written in divergence form as
$$q\_t + b(t)q\_x + (D(t,\,x)q\_x)\_x = 0,$$
so Nash's theorem (which applies to divergence-form equations) implies that $p\_x$ is Holder continuous under mild hypotheses on the coefficients (boundedness and measurability).


If the coefficients are more regular (e.g. in Holder classes) then parabolic Schauder estimates (e.g. in the book of Lieberman) give higher regularity of $p$, the general principle being that $p$ has twice as many spatial derivatives as temporal ones (by the scaling of the equation).


