---
abstract: 'We consider uniformly parabolic equations and inequalities of second order in the non-divergence form with drift $$-u_{t}+Lu=-u_{t}+\sum_{ij}a_{ij}D_{ij}u+\sum b_{i}D_{i}u=0\,(\geq0,\,\leq0)$$ in some domain $\Omega\subset \mathbb{R}^{n+1}$. We prove a variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate with $L^{p}$ norm of the inhomogeneous term for some number $p<n+1$. Based on it, we derive the growth theorems and the interior Harnack inequality. In this paper, we will only assume the drift $b$ is in certain Morrey spaces defined below which are critical under the parabolic scaling but not necessarily to be bounded. This is a continuation of the work in [@GC].'
address: 'Department of Mathematics, The University of Chicago, 5734 South University Avenue, Chicago, IL 60615, U.S.A'
author:
- Gong Chen
date: '05/27/15'
title: 'Non-divergence Parabolic Equations of Second Order with Critical Drift in Morrey Spaces'
---

Introduction
============

General Introduction
--------------------

The qualitative properties of solutions to partial differential equations have been intensively studied for a long time. Following [@GC], in this note, we continue our discussion on the qualitative properties of solutions to the uniform parabolic equation of non-divergence form with drift, $$-u_{t}+Lu:=-u_{t}+\sum_{ij}a_{ij}D_{ij}u+\sum_{i}b_{i}D_{i}u=0\label{eq:01}$$ and the associated inequalities: $-u_{t}+Lu\geq 0$ and $-u_{t}+Lu\leq 0$. Throughout the paper, we use the notations $D_{i}:=\frac{\partial}{\partial x_{i}},\, D_{ij}:=\frac{\partial^{2}}{\partial x_{i}\partial x_{j}}$ and $u_{t}:=\frac{\partial u}{\partial t}$. We assume $b=\left(b_{1},\ldots,b_{n}\right)$ and $a_{ij}$’s are real measurable, $a_{ij}$’s also satisfy the [*uniform parabolicity condition*]{} $$\forall\xi\in\mathbb{R}^{n},\,\nu^{-1}\left|\xi\right|^{2}\leq\sum_{i,j=1}^{n}a_{ij}(X)\xi_{i}\xi_{j},\,\,\,\,\,\,\,\sum_{i,j=1}^{n}a^{2}_{ij}\leq \nu^{2{}}\label{eq:02}$$ with some constant $\nu\geq1$, $\forall X=(x,t)$ in the domain of definition $\Omega \subset \mathbb{R}^{n+1}$.

For the drift $b$, we will only require it is in certain Morrey spaces which are critical under the parabolic scaling. To formulate our setting more precisely, we define Morrey spaces as following: given some constants $p,q\geq1$ and $\alpha\geq0$ satisfying $$\frac{n}{p}+\frac{2}{q}-\alpha=1,\label{eq:04}$$ on the domain of definition $\Omega$, we define $$M_{p,q}^{\alpha}(\Omega):=\left\{ f\in L_{x}^{p}L_{t}^{q}(\Omega);\,||f||_{M_{p,q}^{\alpha}(\Omega)}:=\sup_{Q_{r}\subset\Omega,\,r>0}r^{-\alpha}||f||_{L_{x}^{p}L_{t}^{q}(Q_{r})}<\infty\right\} \label{eq:24}$$ where $Q_{r}$ is the standard parabolic cylinder defined in Definition \[cylinder\] and $$\left\Vert f\right\Vert _{L_{x}^{p}L_{t}^{q}}:=\left(\intop\left[\int\left|f(x,t)\right|^{q}dt\right]^{\frac{p}{q}}dx\right)^{\frac{1}{p}}.$$ We will focus on a particular case with $p=q=\frac{1}{\alpha}=n+1$, $$b\in M_{n+1,n+1}^{\frac{1}{n+1}}(\Omega)$$ with $$\sup_{Q_{r}\subset\Omega,\,r>0}\left[\frac{1}{r}\intop_{Q_{r}}|b|^{n+1}dxdt\right]=: S(\Omega)<\infty.\label{eq:23}$$ By “critical”, we mean that with the $M_{n+1}^{\frac{1}{n+1}}(\Omega)$ norm, the drift is scaling invariant under the parabolic scaling: for $k>0$, $$x\rightarrow k^{-1}x,\,\, t\rightarrow k^{-2}t.$$ Indeed, suppose $u$ satisfies $$-u_{t}+\sum_{ij}a_{ij}D_{ij}u+\sum_{i}b_{i}D_{i}u=0.$$ in a domain $Q\in \mathbb{R}^{n+1}$. Then for any constant $k>0$, let $$\tilde{x}=k^{-1}x,\,\tilde{t}=k^{-2}t.$$ Then $\widetilde{u}\left(\tilde{x},\tilde{t}\right)=u\left(r\tilde{x},r^{2}\tilde{t}\right)$ satisfies the equation $$-\widetilde{u}_{\tilde{t}}+\sum_{ij}\widetilde{a_{ij}}D_{ij}\widetilde{u}+\sum_{i}\tilde{b_{i}}D_{i}\widetilde{u}=0,$$ in $Q_{k}:=\{(x,t),\,(kx,k^{2}t)\in Q \}$. Note that $\tilde{b}=kb$, so $$S^{\frac{1}{n+1}}(\Omega_{k})=\left\Vert \tilde{b}\right\Vert _{M_{n+1,n+1}^{\frac{1}{n+1}}(\Omega_{k})}=\left\Vert b\right\Vert _{M_{n+1,n+1}^{\frac{1}{n+1}}(\Omega)}=S^{\frac{1}{n+1}}(\Omega).$$

In general, regarding the scaling, intuitively, there is a competition between the transport term and the diffusion part. One might expect that for the supercritical scaling case, $\frac{n}{p}+\frac{2}{q}-\alpha>1$: the solutions of the equations have discontinuities [@SVZ]. For the critical situation we are considering here, we have Hölder continuous solutions, see Theorem \[hoelder\]. Finally, if the drift is subcritical with respect to the scaling, i.e. $\frac{n}{p}+\frac{2}{q}-\alpha<1$, we expect the solutions will be smooth. We should notice $L_{x,t}^{n+1}$ ($p=q=n+1,\,\alpha=0$) is supercritical with respect to the parabolic scaling. We will discuss a concrete example in the appendix.

We will concentrate on the growth theorems and the interior Harnack inequality for parabolic equations in non-divergence form with critical drift. In order to derive them, we prove a variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate, Theorem \[variant\]. This variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate enable us to estimate the supremum of a solution to $-u_{t}+Lu=f$ in a bounded Lipschitz domain in $\mathbb{R}^{n+1}$ in terms of the Dirichlet data on the boundary and the $L^{p}$ norm of $f$ with some constant $p<n+1$ depending on $n,\,\nu,\,S$.

With the assumptions and preparations above, the main results in this paper are then expressed by Theorems \[variant\] and \[harnack\].

\[space\] Given $p\geq1$, for any open set $\Omega\subset\mathbb{R}^{n+1}$, we define the space $$W_{p}(\Omega):= C(\Omega)\cap W_{p,p}^{2,1}(\Omega),\label{eq:W}$$ where $f\in W_{p,p}^{2,1}(\Omega)$ means $f_{t},\, D_{i}f,\, D_{ij}f\in\left(L_{x}^{p}L_{t}^{p}\right)_{loc}$.

\[variant\] Under the assumptions above, there are constants $p:=p(\nu,\, n,\, S)<n+1$ and $N$ depending on such that if $u\in C(\overline{\Omega})\cap W_{p}(\Omega)$ satisfies $-u_{t}+Lu\geq f$ in $\Omega$ with $f\in L^{p}(\Omega)$, and $u\leq0$ on $\partial_{p}\Omega$, then $$\sup_{\Omega}u\leq Nr^{2-\frac{n+2}{p}}||f||_{L^{p}}\label{abpkt}$$ where $r$ is the diameter of $\Omega$.

With the help of the variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate , we can obtain:

\[harnack\] Suppose $u\in C(\overline{Q_{2r}(Y)})\cap W_{p}(Q_{2r}(Y))$ and $-u_{t}+Lu=0$ in $Q_{2r}(Y)$, $Y=(y,s)\in\mathbb{R}^{n+1}$ and $r>0$. If $u\geq0$, then $$\sup_{Q^{0}}u\leq N\inf_{Q_{r}}u,\label{eq:harc}$$ where $N=N(n,\nu,S)$ and $Q^{0}=B_{r}(y)\times(s-3r^{2},s-2r^{2})$.

Harnack inequalities have many important applications, not only in differential equations, but also in other areas, such as diffusion processes, geometry, etc. Unlike the classical maximum principle, the interior Harnack inequality is far from obvious. For elliptic and parabolic equations with measurable coefficients in the divergence form, it was proved by Moser in the papers [@M61],[@M64]. However, a similar result for non-divergence equations was obtained 15 years later after Moser’s papers by Krylov and Safonov [@KS], [@S80] in 1978-80. Their proofs relied on some improved versions of growth theorems from the book by Landis [@EML]. These growth theorems control the behavior of (sub-, super-) solutions of second order elliptic and parabolic equations in terms of the Lebesgue measure of areas in which solutions are positive or negative. So certainly, if some estimate can directly import the information about measure, it should be useful. In [@FS], Ferretti and Safonov used growth theorems as a common background for both divergence and non-divergence equations and used these three growth theorems to derive the interior Harnack inequality. Even in the one-dimensional case, the Harnack inequality fails for equations of a “joint” structure, which combine both divergence and non-divergence parts. One can find detailed discussion in [@CS13].

At the beginning, the interior Harnack inequality was proved with bounded drift. Later on, this condition was relaxed to subcritical drift $b$. For the subcritical case, we can always rescale the problem. In small scale, the drift will work like a perturbation from the case without drift. But for the critical situation, our common tricks do not work. One can find a historical overview of this progress in [@NU]. For non-divergence elliptic equations of second order, in [@S10], Safonov shown the interior Harnack inequality for the scaling critical case $b\in L^{n}$. In [@GC], the author proved the interior Harnack inequality for parabolic equations of second order in non-divergence form with the drift $b$ in certain Lebesgue spaces which are scaling invariant. In this note which appears as a companion of the earlier paper [@GC], we consider the case when the drift $b\in M_{n+1,n+1}^{\frac{1}{n+1}}(\Omega)$ which is again scaling invariant. Similar results for both divergence form elliptic and parabolic equations are presented in [@NU].

We will follow the unified approach to growth theorems and the interior Harnack inequality developed in [@FS]. For this purpose, we need to prove three growth theorems and derive the interior Harnack inequality as a consequence for parabolic equations with critical drift formulated as above. In order to take the measure conditions into account and help us carry out growth theorems, we discuss the variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate . Although we only consider the case $b\in M_{n+1,n+1}^{\frac{1}{n+1}}(\Omega)$, one can see from the proofs, our approach works well for other pairs $(p,q,\alpha)$ satisfying condition with $\alpha>0$ provided the associated standard Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate holds, see Sections 2 and 3. For the sake of simplicity, we assume that all functions (coefficients and solutions) are smooth enough. It is easy to get rid of extra smoothness assumptions by means of standard approximation procedures, see Section 7. We should notice that it is important to have appropriate estimates for solutions with constants depending only on the prescribed quantities, such as the dimension $n$, the parabolicity constant, etc., but not depending on “additional” smoothness.

This paper is organized as follows: In Section 1, we introduce our basic assumptions and notations. In Section 2, we formulate a weak version of the classical maximum principle, the Alexandrov-Bakelman-Pucci-Krylov-Tso estimate, and some consequences of it. In Section 3, we establish a variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso based on some estimates of Green’s function. In Sections 4, 5, 6, we formulate and prove three growth theorems and prove the interior Harnack inequality. Finally, in Section 7, we use approximation to show all results are valid without smoothness assumption. In the appendix, an example of loss of continuity of the solution to a parabolic equation with drift $b\in L_{x,t}^{n+1}$ will be presented.

Notations:
----------

In this paper, we use summation convention.

“$A:=B$” or “$B=:A$” is the definition of $A$ by means of the expression $B$.

$\mathbb{R}^{n}$ is the n-dimensional Euclidean space, $n\geq1$, with points $x=(x_{1},\ldots,x_{n})^{t}$, where $x_{i}$’s are real numbers. Here the symbol $t$ stands for the transposition of vectors which indicates that vectors in $\mathbb{R}^{n}$ are treated as column vectors. For $x=(x_{1},\ldots,x_{n})^{t}$ and $y=(y_{1},\ldots,y_{n})^{t}$ in $\mathbb{R}^{n}$, the scalar product $(x,y):=\Sigma x_{i}y_{i}$, the length of $x$ is $|x|:=(x,x)^{\frac{1}{2}}$.

For a Borel set $\Gamma\subset\mathbb{R}^{n}$, $\bar{\Gamma}:=\Gamma\cup\partial\Gamma$ is the closure of $\Gamma$, $|\Gamma|$ is the n-dimensional Lebesgue measure of $\Gamma$. Sometimes we use the same notation for the surface measure of a subset $\Gamma$ of a smooth surface $S$.

For real numbers $c$, we denote $c_{+}:=\max(c,0)$, $c_{-}:=\max(-c,0)$.

In order to formulate our results, we need some standard definitions and notations for the setting of parabolic equations.

\[cylinder\] Let $Q$ be an open connected set in $\mathbb{R}^{n+1}$, $n\geq1$. The parabolic boundary $\partial_{p}Q$ of $Q$ is the set of all points $X_{0}=(x_{0},t_{0})\in\partial Q$, such that there exists a continuous function $x=x(t)$ on the interval $[t_{0},t_{0}+\delta)$ with values in $\mathbb{R}^{n}$, such that $x(t_{0})=x_{0}$ and $(x(t),t)\in Q$ for all $t\in(t_{0},t_{0}+\delta)$. Here $x=x(t)$ and $\delta>0$ depend on $X_{0}$. In particular, for cylinders $Q_{U}=U\times(0,T)$ with $U\subset \mathbb{R}^n$, the parabolic boundary $\partial_{p}Q_{U}:=\left(\partial_{x}Q_{U}\right)\cup\left(\partial_{t}Q_{U}\right)$, where $\partial_{x}Q_{U}:=(\partial U)\times(0,T)$, $\partial_{t}Q_{U}:=U\times\{0\}$.

We will use the following notation for the “standard” parabolic cylinder. For $Y=(y,s)$ and $r>0$, we define $Q_{r}(Y):=B_{r}(y)\times(s-r^{2},s),$ where $B_{r}(y):=\{x\in\mathbb{R}^{n}:|x-y|<r\}$.

Preliminaries
=============

In this section, we briefly discuss some well-known theorems and results which are crucial for us to carry out the discussion in the later parts of this paper. We use the notation $u\in W_{n+1}$ in the sense of Definition \[space\].

\[thm:(Alexandrov-Bakelman-Pucci-Krylo\]\[thm:ABPKT\] Suppose $u\in W_{n+1}(\Omega)$, $\Omega\subset Q_{r}$ and $-u_{t}+Lu\geq f$. If $\sup_{\partial_{p}\Omega}u\leq0$, then $$\sup_{\Omega}u\leq N\left(r^{\frac{n}{n+1}}+\left\Vert b\right\Vert _{L^{n+1}}^{n}\right)||f||_{L^{n+1}}\label{eq:AK-1}$$ where $N=N(n,\nu)$.

One can find the detailed proof of the above standard version of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate in [@AIN] and [@GL2]. From , one can see if we take $r^{\frac{n}{n+1}}$ out of the bracket, we will have $$\sup_{\Omega}u\leq\sup_{\partial\Omega}u+Nr^{\frac{n}{n+1}}\left(1+\left\Vert b\right\Vert^{n} _{M_{n+1}^{\frac{1}{n+1}}(\Omega)}\right)||f||_{L^{n+1}}.$$ So it is natural to consider the case $b\in M_{n+1}^{\frac{1}{n+1}}(\Omega)$.

In [@AIN], Nazarov shown the Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate holds for the drift $b\in L_{x}^{p}L_{t}^{q}$, i.e., $$\left\Vert b\right\Vert _{L_{x}^{p}L_{t}^{q}}=\left(\intop\left[\int\left|b(x,t)\right|^{q}dt\right]^{\frac{p}{q}}dx\right)^{\frac{1}{p}}<\infty,$$ for $$\frac{n}{p}+\frac{2}{q}\leq1,\,\,\,\,p,q\geq1.$$ The proof was based on Krylovs ideas and methods [@NVK]. We believe if $b$ is in other scaling invariant Morrey space, and the Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate holds for the corresponding case, then the proofs in this note also hold.

\[thm:MP\] Let $Q$ be a bounded open set in $\mathbb{R}^{n+1}$, and let a function $u\in C^{2,1}\left(\bar{Q}\backslash\partial_{p}Q\right)\cap C(\bar{Q})$ satisfy the inequality $-u_{t}+Lu\geq0$ in $Q$. Then $$\sup_{Q}u=\sup_{\partial_{p}Q}u\label{eq:-1}$$

As an easy consequence of the maximal principle and the Alexandrov-Bakelman-Pucci-Krylov-Tso estimate, we have the well-known comparison principle.

\[thm:Comparison\]Let $Q$ be a bounded domain in $\mathbb{R}^{n+1}$, $u,\, v\in C^{2,1}\left(\bar{Q}\backslash\partial_{p}Q\right)\cap C(\bar{Q})$, $-u_{t}+Lu\leq-v_{t}+Lv$ in $Q$, and $u\geq v$ on $\partial_{p}Q$, then $u\geq v$ on $\bar{Q}$.

A variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso Estimate
===========================================================

Estimates of Green’s function: 
-------------------------------

In this subsection, we show some estimates for Green’s function following [@FSt; @GL] in order to show the variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate .

We consider the non-divergence second order parabolic equations of the form

$$-u_{t}+Lu=-u_{t}+\sum_{i,j=1}^{n}a_{ij}D_{ij}u+\sum_{i=1}^{n}b_{i}D_{i}u\label{eq:21}$$

defined on some Lipschitz domain $\Omega\subset\mathbb{R}^{n+1}$.

We define Green’s function $G$ : $\Omega\times\Omega\rightarrow\mathbb{R}$ satisfies following properties: if $u(x,t)$ has the form

$$u(x,t)=\intop_{\Omega}G(x,t,y,s)f(y,s)\,dyds\label{eq:22}$$

then it solves the Dirichlet problem: $$-u_{t}+Lu=-f,$$ and $$u=0$$ on $\partial_{p}\Omega$. Throughout this subsection, we will assume all coefficients are smooth, then the existence of Green’s function is guaranteed. We will adapt the ideas in [@GL] to our parabolic setting. The key step is to verify Lemma 2.1 in [@GL] holds in parabolic case with $b\in M_{n+1}^{\frac{1}{n+1}}(\Omega)$. We formulate the following lemma which is similar to Lemma 2.1 in [@GL] under condition .

\[lem:Green1\]Let $(x_{1},t_{1})\in\Omega$, then we can choose $1>\rho>0$ depending on $\nu$ ,$n$ and $S$ such that $Q_{2\rho}(x_{1},t_{1})\subset\Omega$ and there is a constant $C$ depends on $\nu$, $n$ and $S$, we have

$$\left(\intop_{Q_{\rho}(x_{1},t_{1})}G(x,t,y,s)^{\frac{n+1}{n}}\,dyds\right)^{\frac{n}{n+1}}\leq\frac{C}{\rho^{\frac{n+2}{n+1}}}\intop_{Q_{\frac{\rho}{2}}(x_{1},t_{1})}G(x,t,y,s)\,dyds\label{eq:25}$$

for any $(x,t)\in\Omega$ and $t\leq t_{1}$ where $G$ is Green’s function defined as .

Without loss of generality, we may assume $(x_{1},t_{1})=(0,0)$. We will also use $Q_{r}$ to denote $Q_{r}(0,0)$. Clearly, it will be sufficient to show that

$$\intop_{Q_{\rho}}G(x,t,y,s)f(y,s)\,dyds\leq\frac{C}{\rho^{\frac{n+2}{n+1}}}\intop_{Q_{\frac{\rho}{2}}}G(x,t,y,s)\,dyds\label{eq:26}$$

for any non-negative function $f\in L^{n+1}(Q_{\rho})$ with $\intop_{Q_{\rho}}f^{n+1}=1$. We fix such a $f$, and then construct $u_{1}$ and $u_{2}$ as following:

$$u_{1}(x,t)=\intop_{Q_{\rho}}G(x,t,y,s)f(y,s)\,dyds,\label{eq:27}$$

$$u_{2}(x,t)=\intop_{Q_{\frac{\rho}{2}}}G(x,t,y,s)\,dyds.\label{eq:28}$$

We define $$\eta(x,t):=1-\frac{1}{4\rho^{2}}\left\Vert x\right\Vert ^{2}-\frac{1}{4\rho^{2}}t.\label{eq:29}$$ There are positive constants $N_{1}$, $\delta$ and $q$ (determined by the parabolicity $\nu$ and the dimension $n$), such that

$$-D_{t}\eta^{q}+\sum_{ij}a_{ij}D_{ij}\eta^{q}\geq\begin{cases}
0 & Q_{2\rho}\backslash Q_{\frac{\rho}{2}}\\
-N_{1}\rho^{-2} & Q_{\frac{\rho}{2}}
\end{cases},\label{eq:210}$$

and $$\eta^{q}>\delta\label{eq:211}$$ in $Q_{\rho}$.\
Let $u$ be the solution of the Dirichlet problem $$-u_{t}+Lu=-f$$ in $Q_{2\rho}$ and $$u=0$$ on $\partial_{p}Q_{2\rho}$.

By the Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate we have $$u\leq NN_{2}(\rho)^{\frac{n}{n+1}}=k_{1}\rho^{\frac{n}{n+1}},\label{eq:212}$$ where $k_{1}$ depends on $\nu$, $S$ and $n$.

We set $C=\frac{2N_{1}NN_{2}}{\delta}$, $$w:=u_{1}-\frac{C}{\rho^{\frac{n+2}{n+1}}}u_{2,}\label{eq:213}$$

$$\bar{w}:=w-u+\frac{2NN_{2}\rho^{\frac{n}{n+1}}}{\delta}\eta^{q}\label{eq:214}$$

and

$$M:=\max\left\{ 0,\sup_{\partial_{p}Q_{2\rho}}w\right\} .\label{eq:215}$$

Then by some computations, we know that $$-\overline{w}_{t}+L\overline{w}\geq N_{3}|b|\label{eq:216}$$ in $Q_{2\rho}$ where $N_{3}$ depends on $\nu$, $S$ and $n$. And $\overline{w}\leq M$ on the parabolic boundary of $Q_{2\rho}$. By the Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate again, we have $$\overline{w}\leq NN_{2}\rho S{}^{\frac{1}{n+1}}+M\label{eq:217}$$ in $Q_{2\rho}$. In $Q_{\rho}$ with $\eta^{q}>\delta$, we obtain

$$w\leq NN_{2}\rho S{}^{\frac{1}{n+1}}+M+NN_{2}(\rho)^{\frac{n}{n+1}}-\frac{2NN_{2}\rho^{\frac{n}{n+1}}}{\delta}\eta^{q}\leq M\label{eq:218}$$

when we take $\rho$ small enough since $\rho<\rho{}^{\frac{n}{n+1}}$ when $\rho<1$. Here the smallness condition only depends on prescribed constants.

By our construction, it is clear that $$-w_{t}+Lw=0\label{eq:219}$$ in $\Omega\backslash Q_{\rho}$. $$w\leq M\label{eq:220}$$ in $\overline{Q_{\rho}}$ and the parabolic boundary of $\Omega$. By the maximal principle, we have $M=0$. So $w\leq0$ in all $(x,t)\in\Omega$ with $t\leq0$. Therefore, after we decipher $w$, we get for $\rho$ small enough (here the smallness condition only depends on $\nu$, $S$ and $n$), we conclude that $$\intop_{Q_{\rho}}G(x,t,y,s)f(y,s)\,dyds\leq\frac{C}{\rho^{\frac{n+2}{n+1}}}\intop_{Q_{\frac{\rho}{2}}}G(x,t,y,s)\,dyds,$$ which implies $$\left(\intop_{Q_{\rho}(x_{1},t_{1})}G(x,t,y,s)^{\frac{n+1}{n}}\,dyds\right)^{\frac{n}{n+1}}\leq\frac{C}{\rho^{\frac{n+2}{n+1}}}\intop_{Q_{\frac{\rho}{2}}(x_{1},t_{1})}G(x,t,y,s)\,dyds.$$

With Lemma \[lem:Green1\], we can proceed to the prove of the variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate similar to results in [@GL] and [@FSt]. Since every quantity we are considering here is scaling invariant, we may rescale our setting to a domain with diameter $1$. So it suffices estimate the integrability of Green’s function in a domain with diameter $1$.

\[thm:Green2\]Under the same assumptions as above, then there are constants $q:=q((\nu,\, n,\, S)>\frac{n+1}{n}$ and $C$ only depending on $\nu$, $n$ and $S$ such that

$$\left(\intop_{\Omega}G(x,t,y,s)^{q}\,dyds\right)^{\frac{1}{q}}\leq C\label{eq:221}$$

where $G$ is Green’s function and the diameter of $\Omega$ is $1$.

By our assumptions $\Omega$ is bounded. Suppose $$\Omega'=\left\{ X:=(x,t)\in\Omega,\, d_{\Omega}(X)<\frac{1}{2}\right\} \label{eq:222}$$ where

$$d_{\Omega}(X):=\sup\left\{ \rho>0:\, Q_{\rho}(X)\subset\Omega\right\} .$$ For arbitrary $(x',t')\in\Omega'$ and let $r<\frac{1}{2}$, then $Q_{r}(x',t')\subset\Omega$. For each $Q_{r}(x',t')$, we can use finite many $Q_{\rho}(x_{i},t_{i})$ to cover it, where $\rho$ is chosen small enough so that the conditions of Lemma \[lem:Green1\] are satisfied. It is clear that the number of $Q_{\rho}(x_{i},t_{i})$ can be bounded by a positive constant $C_{1}$ depending on $\nu$, $n$ and $S$ and the diameter of $\Omega$. Let $\rho=C_{2}r$, then we have

$$\left(\intop_{Q_{r}(x',t')}G(x,t,y,s)^{\frac{n+1}{n}}\,dyds\right)^{\frac{n}{n+1}}\leq\frac{C_{3}}{\rho^{\frac{n+2}{n+1}}}\intop_{Q_{r}(x',t')}G(x,t,y,s)\,dyds\label{eq:223}$$

where $C_{3}$ depends on $C_{1}$, $C_{2}$ and $C$ in Lemma \[lem:Green1\].

We rewrite then inequality as $$\left(\fint_{Q_{r}(x',t')}G(x,t,y,s)^{\frac{n+1}{n}}\,dyds\right)^{\frac{n}{n+1}}\leq\ C_{3}\fint_{Q_{r}(x',t')}G(x,t,y,s)\,dyds\label{eq:2233}$$

By Gehring’s lemma [@MM], we get

$$\left(\intop_{\Omega'}G(x,t,y,s)^{q}\,dyds\right)^{\frac{1}{q}}\leq C_{4}\intop_{\Omega'}G(x,t,y,s)\,dyds\label{eq:224}$$

for some $$q>\frac{n+1}{n},$$ where $C_{4}$ depends on $C_{1}$, $C_{2}$ and $C$ in Lemma \[lem:Green1\]. $$\intop_{\Omega'}G(x,t,y,s)\,dyds< C_{5}$$ where the constant $C_{5}$ depends on $\nu$, $n$ and $S$ by the standard Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate, Theorem \[thm:(Alexandrov-Bakelman-Pucci-Krylo\]. Hence there is a constant $C'$ such that

$$\left(\intop_{\Omega'}G(x,t,y,s)^{q}\,dyds\right)^{\frac{1}{q}}\leq C'\label{eq:225}$$

where $C'$ only depends on $\nu$, $n$ and $S$. In order to get the same result for $\Omega$, we can extend the all coefficients to a domain $\widetilde{\Omega}$ with comparable quantities, such that $\Omega\subset\widetilde{\Omega}$ and $\forall X:=(x,t)\in\Omega$, $d_{\widetilde{\Omega}}(X)<1$. Then Green’s function $\widetilde{G}$ for $\widetilde{\Omega}$ will satisfy the same result as . And we know $\widetilde{G}\geq G$ by the maximal principle and the comparison principle. Hence we have

$$\left(\intop_{\Omega}G(x,t,y,s)^{q}\,dyds\right){}^{\frac{1}{q}}\leq C\label{eq:226}$$

where $C$ depends on $\nu$, $n$ and $S$.

We can see from the proof, we just need $b\in M_{n+1}^{\alpha}(\Omega)$ for $\alpha$ positive, the above arguments work too provided the associated Aleksandrov-Bakelman-Pucci-Krylov-Tso Estimate holds. The point we choose $\alpha=\frac{1}{n+1}$ is that the space is scaling invariant under the parabolic scaling.

A variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso Estimate:
------------------------------------------------------------

In this subsection, we still keep all the assumptions above. We will obtain a variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate for a solution of the second order parabolic equation differential equation $-u_{t}+Lu=f$ in a bounded Lipschitz domain in $\mathbb{R}^{n+1}$ with Dirichlet problem condition. We estimate the maximal of $u$ in terms of the Dirichlet data (boundary value) on the boundary and the $L^{p}$ norm of $f$ with $p<n+1$. For elliptic and parabolic equations without drift, one can find some references in [@CS; @EML; @GL; @FSt].

Under the same assumptions above, define $p=s \frac{q}{q-1}$, where $q>\frac{n+1}{n}$ is the constant from Theorem \[thm:Green2\]. Suppose $u\in C(\overline{\Omega})\cap W_{p}(\Omega)$ satisfies $-u_{t}+Lu\geq f$ in $\Omega$ and $u\leq0$ on $\partial_{p}\Omega$. Then there is a constant $N$ depends on such that

$$\sup_{\Omega}u\leq Nr^{2-\frac{n+2}{p}}||f||_{L^{p}}\label{eq:31}$$

where $r$ is the diameter of $\Omega$.

Without loss of generality, we may assume the diameter of $\Omega$ is 1. We first prove the result holds in case 1: when $u$ is smooth and the coefficients of the equations are smooth (Hölder). Then we use approximation to show in case 2: under general condition, the estimate holds.\
*Case 1[<span style="font-variant:small-caps;">:</span>]{}* we can represent $u$ using:

$$u(x,t)=\intop_{\Omega}G(x,t,y,s)\left[u_{t}-Lu\right]\,dyds+\intop_{\partial_{p}\Omega}\tilde{G}(x,t,y,s)u\,dyds\label{eq:32}$$

where $G$ is Green’s function and $\tilde{G}$ is from the Riesz representation theorem, since $\Omega$ is a Lipchitz domain. It is clear both of $G$ and $\tilde{G}$ are non-negative, hence under our assumptions we have

$$u(x,t)\leq\intop_{\Omega}G(x,t,y,s)\left|f\right|dyds\label{eq:33}$$

then by Theorem \[thm:Green2\] and Hölder inequality, it is clear that $$\sup_{\Omega}u\leq N\left\Vert f\right\Vert {}_{L^{p}},\label{eq:34}$$ where $N$ depends on $\nu$, $n$ and $S$.\
*Case 2[<span style="font-variant:small-caps;">:</span>]{}* First we assume $u\in C(\overline{\Omega})\cap W_{p}^{2,1}(\Omega)$ and other functions are still smooth. Under these conditions, we pick up a sequence $u^{i}$ which is smooth such that $u^{i}\rightarrow u$. It is clear there is a sequence $\{f_{i}\}_{i}^{\infty}$ such that $-(u_{i})_{t}+Lu_{i}\geq f_{i}$ and $f_{i}\rightarrow f$ in $L^{p}(\Omega)$. So it is clear that $$\sup_{\Omega}u\leq N||f||_{L^{p}}.\label{eq:35}$$ Finally, in general situation. Suppose $a_{ij}^{k}\rightarrow a_{ij}$ and $b_{i}^{k}\rightarrow b_{i}$ almost everywhere as $k\rightarrow\infty$. Define $$L_{k}u=\sum_{ij}a_{ij}^{k}D_{ij}u+\sum b_{i}^{k}D_{i}u,$$ and $$f_{k}=f+(L-L_{k})u.\label{eq:36}$$ It is clear that $f_{k}\rightarrow f$ in $L^{p}(\Omega)$. Therefore $$\sup_{\Omega}u\leq N||f||_{L^{p}}\label{eq:37}$$ holds.

Finally, after rescaling, we obtain in the most general setting, $$\sup_{\Omega}u\leq Nr^{2-\frac{n+2}{p}}||f||_{L^{p}}.$$

First Growth Theorem
====================

Suppose $R$ is the region in a cylinder where a subsolution $u$ of our equation is positive. The first growth theorem, Theorem \[thm:(First-Growth-Theorem)\], basically tells us if the measure of $R$ is small, then the maximal value of $u$ over half of the cylinder is strictly less than the maximal value over the whole cylinder. In other words, it gives us some quantitative decay properties. The variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate enables us to import information about the measure into our estimates.

\[thm:(First-Growth-Theorem)\] Let a function $u\in C^{2,1}(\overline{Q_{r}})$ where $r>0$ and $Q_{r}=Q_{r}(Y)$, in $\mathbb{R}^{n+1}$ containing $Y:=(y,s)$. Suppose $-u_{t}+Lu\geq0$ in $Q_{r}$, then $\forall\beta_{1}\in(0,1)$, there exists $0<\mu<1$ such that if we know $$\left|\left\{ u>0\right\} \cap Q_{r}(Y)\right|\leq\mu|Q_{r}(Y)|,\label{eq:129}$$ then

$$\mathcal{M}_{\frac{r}{2}}(Y)\leq\beta_{1}\mathcal{M}_{r}(Y),\label{eq:130}$$

where $$\mathcal{M}_{r}(Y):=\max_{Q_{r}(Y)}u_{+}$$ We also notice that $\beta_{1}\rightarrow0^{+}$ as $\mu\rightarrow0^{+}$.

First of all, we make some reductions. In our problem, we want to show under some conditions, given $-u_{t}+Lu\geq0$ in a cylinder $Q_{r}(Y)$, and some information about the set $\{u\leq0\}$, we want to show that $$\mathcal{M}_{\frac{r}{2}}(Y)\leq\beta_{1}\mathcal{M}_{r}(Y).$$ Clearly, in order to derive the above estimate, we only need to consider positive part of $u$. We observe that to obtain the above estimate, it actually suffices to show $$u(Y)\leq\beta_{1}\mathcal{M}_{r}(Y),\label{eq:extra}$$ for some $\beta_{1}\in(0,1)$. Indeed, for an arbitrary point $Z\in Q_{\frac{r}{2}}(Y)$, we notice $Q_{\frac{r}{2}}(Z)\subset Q_{r}(Y)$, we can apply the above estimate to $Q_{\frac{r}{2}}(Z)$ with $Y$ replaced by $Z$ and $r$ replaced by $\frac{r}{2}$ with some measure condition $\mu'$. In consistent with the measure condition in the first growth theorem, we also observe that $$\left|\left\{ u>0\right\} \cap Q_{\frac{r}{2}}(Z)\right|\leq\left|\left\{ u>0\right\} \cap Q_{r}(Y)\right|\leq\mu\left|Q_{r}\right|=2^{n+2}\mu\left|Q_{\frac{r}{2}}(Z)\right|.$$ So we just need to take $\mu=2^{-n-2}\mu'$ for the measure condition in the first growth theorem.

Since every quantity is scaling invariant, we might assume $r=1$. And we can multiply $u$ by a constant, so without loss of generality, we can also assume $\mathcal{M}_{1}(Y)=1$. Also we assume $u(Y)>0$, otherwise the result is trivial. $$v(X)=v(x,t)=u(x,t)+t-s-|x-y|^{2}\label{eq:43}$$ in $Q:=\left\{ v>0\right\} \cap Q_{1}(Y)$. Clearly, $Q\neq\emptyset$ since $v(Y)=u(Y)>0$ and $Y\in\partial Q_{1}(Y)$. It is easy to see that $v\leq u$ in $Q$. By the measure condition, we have $$\left|Q\right|\leq\left|\left\{ u>0\right\} \cap Q_{1}(Y)\right|\leq\mu_{1}|Q_{1}(Y)|\leq\mu_{1}.$$ Note that $v\leq0$ on $\partial_{p}Q_{1}(Y)$, so $v=0$ on $\partial_{p}Q$. Since $-u_{t}+Lu\geq0$, we know that $$(-\partial_{t}+L)v\geq0-1-2trace(a_{ij})-2|b|\geq-1-2nv^{-1}-2|b|.\label{eq:44}$$ By the variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate with some constant $p<n+1$ and Hölder inequality, $$u(Y)\leq N(\nu,n,S)\left\Vert -1-2nv^{-1}-2|b|\right\Vert _{L^{p}(Q)}\leq N_{1}(\nu,n,S)\left(\mu^{\frac{1}{p}}+S^{\frac{1}{n+1}}\mu^{\frac{1}{p}-\frac{1}{n+1}}\right).\label{eq:45}$$ Now we can pick $\mu$ small enough so that for fixed $\beta_{1}$, then we have $$u(Y)<\beta_{1}.\label{eq:46}$$ It is also clear from the construction, $\beta_{1}\rightarrow0^{+}$ as $\mu\rightarrow0^{+}$.

With the first growth theorem, we can do the following useful argument which is helpful for us to find a non-degenerate point to build a bridge between two regions we are interested in. Without loss of generality, we still assume $r=1$, for $X\in Q_{1}(Y)$, we define $$d(X):=\sup\left\{ \rho>0:\, Q_{\rho}(X)\subset Q_{1}(Y)\right\} .\label{eq:47}$$ Roughly here $d$ plays roles of weights with which we can make sure the point we are interested in is not degenerate, i.e., it is in the interior of the cylinder. For $\gamma>0$, we consider $d^{\gamma}u(x)$ instead of $u(x)$. $d^{\gamma}u(x)$ is a continuous function in $\overline{Q_{1}(Y)}$. Clearly, $d(Y)=1$, we obtain $$u(Y)=d^{\gamma}u(Y)\leq M:=\sup_{Q_{1}(Y)}d^{\gamma}u.\label{eq:48}$$ By our construction, $d^{\gamma}u$ vanishes on $\partial_{p}Q_{1}$, so $\exists X_{0}\in\overline{Q_{1}(Y)}\backslash\partial_{p}Q_{1}$ such that $$M=d^{\gamma}u(X_{0}).\label{eq:49}$$ Let $r_{0}:=\frac{1}{2}d(X_{0})$, we consider the intermediate region $Q_{r_{0}}(X_{0})$, In this region, we have $$\forall X\in Q_{r_{0}}(X_{0}),\,\, d(X)\geq r_{0}.$$ Therefore, we conclude that $$\sup_{Q_{r_{0}}(X_{0})}u\leq r_{0}^{-\gamma}\sup_{Q_{r_{0}}(X_{0})}d^{\gamma}u\leq r_{0}^{-\gamma}M\leq2^{\gamma}u(X_{0}).\label{eq:410}$$ Now, we define $v=u-\frac{1}{2}u(X_{0})$, then $$v(X_{0})=\frac{1}{2}u(X_{0})\geq2^{-1-\gamma}\sup_{Q_{r_{0}}(X_{0})}u>2^{-1-\gamma}\sup_{Q_{r_{0}}(X_{0})}v.$$ From the first growth theorem, Theorem \[thm:(First-Growth-Theorem)\], we know $\exists\mu(n,\nu,\gamma,S)\in(0,1]$ such that Theorem \[thm:(First-Growth-Theorem)\] holds with $\beta_{1}=2^{-1-\gamma}$. Now the above inequality tells us that $v$ does not satisfy the measure condition in the first growth theorem. So $$\left|\left\{ v>0\right\} \cap Q_{r_{0}}(X_{0})\right|=\left|\left\{ u>\frac{1}{2}u(X_{0})\right\} \cap Q_{r_{0}}(X_{0})\right|>\mu\left|Q_{r_{0}}(X_{0})\right|.\label{eq:411}$$ Now, we can show an integral estimate which is equivalent to the first growth theorem.

\[thm:Lepsilon\]Let a function $u\in C^{2,1}(Q_{r})$, where $Q_{r}:=Q_{r}(Y)$, $Y=(y,s)\in\mathbb{R}^{n+1}$, $r>0$. If $-u_{t}+Lu\geq0$ in $Q_{r}$, then for arbitrary $p>0$, we obtain $$u_{+}^{p}(Y)\leq\frac{N}{\left|Q_{r}\right|}\intop_{Q_{r}}u_{+}^{p}\,dX,\label{eq:412}$$ where $N$ only depends on $n,\,\nu,\, S,\, p$.

Since the quantities we are considering are scaling invariant, we might assume $r=1$. By the similar argument as above, we choose $$\gamma=\frac{n+2}{p}.\label{eq:413}$$ We get $$\left|\left\{ u>\frac{1}{2}u(X_{0})\right\} \cap Q_{r_{0}}(X_{0})\right|>\mu\left|Q_{r_{0}}(X_{0})\right|,\label{eq:414}$$ where $\mu=\mu(n,\nu,\gamma,S)\in(0,1]$. With the same notations as above we have $$\begin{aligned}
u_{+}^{p}(Y) & \leq & M^{p}=\left(u\right)^{p}(X_{0})=\left(2r_{0}\right)^{\gamma p}u^{p}(X_{0})\nonumber \\
 & \leq & \frac{\left(2r_{0}\right)^{\gamma p}}{\left|\left\{ u>\frac{1}{2}u(X_{0})\right\} \cap Q_{r_{0}}(X_{0})\right|}\intop_{\left\{ u>\frac{1}{2}u(X_{0})\right\} \cap Q_{r_{0}}(X_{0})}(2u)^{p}\,dX\nonumber \\
 & \leq & \frac{2^{\gamma p+p}r_{0}^{\gamma p}}{\mu\left|Q_{r_{0}}(X_{0})\right|}\intop_{Q_{1}}u_{+}^{p}\,dX.\label{eq:415}\end{aligned}$$ By our construction, it gives $r_{0}^{\gamma p}=r_{0}^{n+2}$, so the above estimate implies $$u_{+}^{p}(Y)\leq\frac{N}{\left|Q_{1}\right|}\intop_{Q_{1}}u_{+}^{p}dX.\label{eq:416}$$

We have seen the first growth theorem implies Theorem \[thm:Lepsilon\]. Actually, we can also obtain the first growth theorem from Theorem \[thm:Lepsilon\]. Indeed, from Theorem \[thm:Lepsilon\], we have $$u_{+}^{p}(Y)\leq\frac{N}{\left|Q_{r}\right|}\intop_{Q_{r}}u_{+}^{p}\,dX\leq N\frac{\left|\left\{ u>0\right\} \cap Q_{r}(Y)\right|}{\left|Q_{r}\right|}\sup_{Q_{r}}u_{+}^{p}\leq\mu N\sup_{Q_{r}}u_{+}^{p}.\label{eq:417}$$ Now it is trivial to see $\beta_{1}\rightarrow0^{+}$ as $\mu\rightarrow0^{+}$ which is in consistent with the conditions in the first growth theorem.

The idea we used to find a non-degenerate point above will be also helpful when we prove the interior Harnack inequality.

Second Growth Theorem
=====================

Before we establish the second growth theorem, Theorem \[thm:(Second-Growth-Theorem).\], we need to prove some intermediate results based on the comparison principle and the Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate. Let us first do some preliminary calculations in order to carry out some comparison arguments.

For fixed number $\alpha>0$ and $0<\epsilon<1$, in the cylinder $Q=B_{r}(0)\times(-r^{2},(\alpha-1)r^{2})$, we can define $$\psi_{0}=\frac{(1-\epsilon^{2})(t+r^{2})}{\alpha}+\epsilon^{2}r^{2}\label{eq:51}$$ and $$\psi_{1}=\left(\psi_{0}-|x|^{2}\right)_{+}\label{eq:52}$$ where $\left(\cdot\right)_{+}$ means positive part of the function. And we also define $$\psi=\psi_{1}^{2}\psi_{0}^{-q}\label{eq:53}$$ for some number $q\geq2$ to be determined later. First of all, we notice $\psi$ is $C^{2,1}$ in $\widetilde{Q}:=\left\{ (x,t)|\,|x|^{2}<\psi_{0},\,-r^{2}<t<(\alpha-1)r^{2}\right\} $. It is clear that $-\psi_{t}+L\psi=0$ if $\psi_{0}\leq|x|^{2}$. Now if $\psi_{0}>|x|^{2}$, by some computations, we obtain $$-\psi_{t}+a_{ij}D_{ij}\psi=\psi_{0}^{-q}\left[8a_{ij}x_{i}x_{j}-4\psi_{1}trace(a_{ij})+\frac{(1-\epsilon^{2})q}{\alpha\psi_{0}}\psi_{1}^{2}-2\frac{(1-\epsilon^{2})}{\alpha}\psi_{1}\right].$$ Set $F_{1}=\frac{2}{\alpha}+8n\nu^{-1}$, and $\xi=\frac{\psi_{1}}{\psi_{0}}$ then

$$-\psi_{t}+a_{ij}D_{ij}\psi\geq\psi_{0}^{1-q}\left[\frac{(1-\epsilon^{2})q}{\alpha}\xi^{2}-F_{1}\xi+8\lambda\right].\label{eq:54}$$

Pick $$q=2+\frac{\alpha}{32(1-\epsilon^{2})},\label{eq:55}$$ so that the quadratic form in is non-negative. Then we can conclude that $$-\psi_{t}+\sum_{ij}a_{ij}D_{ij}\psi\geq0,\,\,\,\,\forall(x,t)\in Q$$ We also notice that $$\psi(x,-r^{2})\leq(\epsilon r)^{-2q+4},\,\,\,\,\forall|x|\leq r,$$ and $$\psi\left(x,(\alpha-1)r^{2}\right)\geq\frac{9}{16}r^{-2q+4},\,\,\,\,\,\forall|x|\leq\frac{r}{2}.\label{eq:56}$$ Finally, we notice that by the monotonicity of $\psi$ with respect to $t\in\left[-r^{2},(\alpha-1)r^{2}\right]$ for $x=0$, we obtain $$\psi\left(0,t\right)\geq\frac{9}{16}r^{-2q+4}.\label{eq:57}$$ Now consider $-u_{t}+Lu\leq0$ , $u>0$ in $\Omega$. In Lemma \[lem:samegrowth\], fist of all, we establish that at least for a cylinder short enough, if we have a lower bound on some interior portion of the bottom, then it has a quantitative lower bound for the same portion on the top of the short cylinder. We also notice that the shortness only depends on the prescribed constants but not $u$. Then we can iterate this process to get a quantitative lower bound for an arbitrary time.

\[lem:samegrowth\]Let $\alpha$ be a positive constant and $-u_{t}+Lu\leq0$, $u>0$ in $\Omega$. Suppose $Q:=B_{r}(0)\times(-r^{2},(\alpha-1)r^{2})\subset\Omega$ and $u>0$ in $B_{r}(0)\times(-r^{2},(\alpha-1)r^{2})$. Then for $\epsilon<\frac{1}{2}$, there are positive constants $C_{1}=C_{1}(n,\nu)$ and $m=m(n,\nu,\alpha)$ such that if $$u\geq\ell\label{eq:58}$$ on $B_{\epsilon r}(0)\times\{-r^{2}\}$, then $$u\geq C_{1}\epsilon^{m}\ell\label{eq:59}$$ on $B_{\epsilon r}(0)\times\{(\alpha-1)r^{2}\}$.

*[<span style="font-variant:small-caps;">Step 1:</span>]{}*[ <span style="font-variant:small-caps;"></span> ]{} $$-\psi_{t}+L\psi\geq b_{i}D_{i}\psi=-4\psi_{1}\psi_{0}^{-q}(b,x)\geq-4|b|r(\frac{1}{2})^{-2q}r^{2-2q}.\label{eq:510}$$ Consider $$v=u-\ell(\epsilon r)^{2q-4}\psi\label{eq:511}$$ It is clear $\psi=0$ for $|x|=r$. So we can conclude that $v\geq0$ on $\partial_{p}\widetilde{Q}$ by the above calculations. Finally, we apply the Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate to $-v$, we get $$v\geq-N(n,\nu,S)\ell\epsilon{}^{-4}S\label{eq:512}$$ in $\widetilde{Q}$. In other words, we get $$u\geq\ell(\epsilon r)^{2q-4}\psi-N(n,\nu,S)\ell\epsilon^{-4}S.\label{eq:513}$$ So we have $$u(x,(\alpha-1)r^{2})\geq\ell\epsilon{}^{-4}\left[\frac{9}{16}(\epsilon)^{2q}-N(n,\nu,S)S\right].\label{eq:514}$$ In particular, with the monotonicity of $\psi$ with respect to $t$ when $x=0$, we obtain $$u(0,t)\geq\ell\epsilon{}^{-4}\left[\frac{9}{16}(\epsilon)^{2q}-N(n,\nu,S)S\right],\,\forall t\in\left[-r^{2},(\alpha-1)r^{2}\right].\label{eq:515}$$ By the similar calculations as above but with the variant of Aleksandrov-Bakelman-Pucci-Krylov-Tso applied to the region $B_{r}(0)\times(-r^{2},(h\alpha-1)r^{2})$, we get $$u(0,(h\alpha-1)r^{2})\geq\ell\epsilon{}^{-4}\left[\frac{9}{16}(\epsilon)^{2q}-N(n,\nu,S)S^{\frac{1}{n+1}}h^{\frac{1}{p}-\frac{1}{n+1}}\right].\label{eq:516}$$ Pick $h=h(n,\nu,S,\epsilon)$ small, we know that $$\left[\frac{9}{16}(\epsilon)^{2q}-N(n,\nu,S)S^{\frac{1}{n+1}}(\alpha h)^{\frac{1}{p}-\frac{1}{n+1}}\right]\geq\frac{1}{2}\epsilon^{2q}.\label{eq:517}$$ So we can conclude $$u(0,t)\geq C_{1}\epsilon^{k}\ell\label{eq:518}$$ for $t\leq(h-1)r^{2}$ where $k=2q-4$ and $C_{1}$ does not depend on $u$.\
*[<span style="font-variant:small-caps;">Step 2:</span>]{}*[ <span style="font-variant:small-caps;"></span> ]{}For arbitrary $t=(\alpha h-1)r^{2}$, and $x\in B_{\epsilon}(0)$, we can use a slanted cylinder with radius $\epsilon$ to connect $B_{\epsilon}(x)\times\{(\alpha h-1)r^{2}\}$ and $B_{\epsilon}(0)\times\{-r^{2}\}$. We can use a change of coordinate to reduce the slanted cylinder to a regular cylinder. We notice that with $k_{i}:=\frac{y_{i}}{s}$ and $\frac{|y_{i}|}{s}=|k_{i}|\leq\frac{|y|}{s}\leq K$, then define $w_{i}=x_{i}-k_{i}t$ and $z=t$. In this coordinate, the slanted cylinder is transformed to a standard cylinder. The equation with respect to the new coordinate is $$-u_{z}+\sum_{ij}a_{ij}D_{w_{i}w_{j}}u+\sum_{i}(b_{i}+k_{i})D_{w_{i}}u\leq0.\label{eq:519}$$ Then we apply the standard cylinder results to the equation with respect to coordinate $(w,z)$. We can do the same argument for all $x\in B_{\epsilon}(0)$. We have find $h_{x}$ such that $$u(x,(\alpha h_{x}-1)r^{2})\geq C_{2}\epsilon^{k}\ell.\label{eq:520}$$ Since $K$ is uniformly bounded above, so indeed, $h_{x}$ and $C_{2}$ only depend on $\epsilon$, $n,\nu,S$. In particular, $h_{x}$ can be uniformly bounded from below. Finally, we take $$h_{0}=\inf_{x\in B_{\epsilon}(0)}h_{x}>0,\label{eq:521}$$ then for $x\in B_{\epsilon}(0)$, $t=(\alpha h_{0}-1)r^{2}$, we obtain $$u(x,t)\geq C_{2}\epsilon^{k}\ell.\label{eq:522}$$

*[<span style="font-variant:small-caps;">Step 3:</span>]{}*[ <span style="font-variant:small-caps;"></span> ]{}Now for the general case, let $\alpha$ be positive constant as above, we can pick $h_{0}$ based on our discussion above. Finally by a simple iteration argument, we get the above result. Therefore $\exists m=m(n,S,\nu,\alpha)$, such that at least we conclude $$u\geq C_{1}\epsilon^{m}\ell\label{eq:523}$$ on $B_{\epsilon}(0)\times\{(\alpha-1)r^{2}\}$.

In fact, in consistent with the Lemma 7.39 in [@GL2], we can show that $$u(x,(\alpha-1)r^{2})\geq C_{3}\epsilon^{m}\ell$$ on $B_{\frac{1}{2}r}(0)\times\{(\alpha-1)r^{2}\}$.

For a fixed point $Y=(y,s)\in\mathbb{R}^{n+1}$ with $s>0$, and $r>0$, we define the slanted cylinder $$V_{r}=V_{r}(Y):=\left\{ X=(x,t)\in\mathbb{R}^{n+1};\,\left|x-\frac{t}{s}y\right|<r,0<t<s\right\} .\label{eq:122}$$

Now the useful slanted cylinder lemma [@FS] follows easily from Lemma \[lem:samegrowth\] after we apply Lemma \[lem:samegrowth\] to $1-u$ after we multiply $u$ by a constant to reduce our problem to the case $1=\sup_{V_{r}(Y)}u$.

\[lem:(Slant-Cylinder-Lemma)\] Let a function $u\in C^{2,1}(\overline{V_{r}})$ satisfy $-u_{t}+Lu\geq0$ in a slanted cylinder $V_{r}$, which is defined in with $Y=(y,s)\in\mathbb{R}^{n+1}$, $s>0$, $r>0$, such that

$$K^{-1}r|y|\leq s\leq Kr^{2}\label{eq:525}$$

where $K>1$ is a constant. In addition, suppose $u\leq0$ on $D_{r}:=B_{r}(0)\times\{0\}$. Then

$$u(Y)\leq\beta_{2}\sup_{V_{r}(Y)}u_{+}\label{eq:526}$$

with a constant $\beta_{2}=\beta_{2}(\nu,n,K,S)<1$.

With the slanted cylinder lemma, we are ready to prove the second growth theorem. The slanted cylinder lemma, i.e., Lemma \[lem:(Slant-Cylinder-Lemma)\] above plays a crucial role in this section to build a connection between different time slides. The second growth theorem helps us construct some control of the oscillation between different time slides. We follow the arguments in [@FS].

After rescaling and translation in $\mathbb{R}^{n+1}$, we reduce our problem to $r=1$, and $(z,\tau)=(0,0)\in\mathbb{R}^{n+1}$. For an arbitrary point $Y'\in Q_{\frac{1}{2}}(Y)$, we can apply the slanted cylinder lemma to the slanted cylinder $V_{\rho}(Y')\subset Q_{1}(Y)$. Note that in this situation, the constant $K$ in slanted cylinder lemma only depends on $\rho$. Therefore, with the parameter $\beta_{2}$ from the slanted cylinder lemma, we have $$u(Y')\leq\beta_{2}\sup_{V_{\rho}(Y')}u_{+}\leq\beta_{2}\sup_{Q_{1}(Y)}u_{+}.$$ The above estimate holds for all $Y'\in Q_{\frac{1}{2}}(Y)$. Therefore we get $$\sup_{Q_{\frac{1}{2}}(Y)}u_{+}\leq\beta_{2}\sup_{Q_{1}(Y)}u_{+}.$$

Now we establish an estimate similar to above with more explicit dependence of the constant on the ratio $\rho/r$.

\[lem:quotient\]Let a function $v\in C^{2,1}\left(\overline{Q_{r}}\right)$ satisfy $v\geq0$, $-v_{t}+Lv\leq0$ in $Q_{r}:=Q_{r}(Y),Y=(y,s)\in\mathbb{R}^{n+1},r>0$. For arbitrary disks $D_{\rho}:=B_{\rho}(z)\times\{\tau\}$ and $D^{0}:=B_{\frac{r}{2}}(y)\times\{\sigma\}$, such that $B_{\rho}(z)\subset B_{r}(y)$ and $$s-r^{2}\leq\tau<\tau+h^{2}r^{2}\leq\sigma\leq s,\label{eq:529}$$ where $h\in(0,1)$ is a constant. Then $$\inf_{D_{\rho}}v\leq\left(\frac{2r}{\rho}\right)^{\gamma}\inf_{D^{0}}v\label{eq:530}$$ where $\gamma=\gamma(n,\nu,h,S)$.

Without loss of generality, we may assume $m:=\inf_{D_{\rho}}v>0$, $r=1$, $z=0$, $\tau=0$, $\sigma=s=h^{2}$. So $D_{\rho}=B_{\rho}(0)\times\{0\}$. We can apply an additional linear transformation along $t$-axis, we can also reduce the proof to the case $h=1$. Now fix the integer $k$ such that $2^{-k-1}<\rho\leq2^{-k}$, and for $j=0,1,\ldots$, and define $y^{j}:=y^{*}+2^{-j}(y-y^{*})$, $B^{j}:=B_{2^{-j}}(y^{j})$, where $y^{*}:=\frac{\rho}{1-\rho}y$, $Y^{j}:=\left(y^{j},4^{-j}\right)$, $Q^{j}:=Q_{2^{-j}}(Y^{j})$, $D^{j}:=B_{2^{-j-1}}(y^{j})\times\{4^{-j}\}$. By construction, $0=y^{*}+\rho(y-y^{*})$, so that $$B_{\rho}(0),\, B^{j}\in\{B_{\theta}\left(y^{*}+\theta(y-y^{*})\right);\,0\leq\theta\leq1\}.$$ Then by the assumption, $B_{\rho}(0)\subset B_{1}(y)$ it follows $\left|y\right|\leq1-\rho$, $\left|y-y^{*}\right|\leq1$, and $$B^{k+1}\subset B_{\rho}(0)\subset B^{k}\subset B^{k-1}\subset\ldots\subset B^{1}\subset B^{0}=B_{1}(y).$$ Now apply Theorem \[thm:(Second-Growth-Theorem).\] to the function $u=1-\frac{1}{m}v$ in $Q^{k}$ with $$r=2^{-k},\,\rho=2^{-k-1},\, Y=Y^{k},\, z=0,\,\tau=0.$$ Then we conclude that $$\sup_{D_{k}}u\leq\sup_{Q_{2^{-k-1}(Y^{k})}}u\leq\beta_{3}\sup_{Q^{k}}u\leq\beta_{3}=\beta_{2}(n,\nu,S,\frac{1}{2})<1,$$ which is equivalent to $$\inf_{D_{\rho}}v=m\leq(1-\beta_{3})^{-1}\inf_{D^{k}}v=2^{\gamma}\inf_{D^{k}}v,$$ where $\gamma:=-\log_{2}(1-\beta_{3})>0$. Similarly, if $k\geq1$, we also have $$\inf_{D^{j}}v\leq2^{\gamma}\inf_{D^{j-1}}v,$$ for $j=1,2,\ldots,k$. Finally we have $$\inf_{D_{\rho}}v\leq2^{\gamma}\inf_{D^{k}}v\leq2^{2\gamma}\inf_{D^{k-1}}v\leq\ldots\leq2^{(k+1)\gamma}\inf_{D^{0}}v\leq\left(\frac{2r}{\rho}\right)^{\gamma}\inf_{D^{0}}v.$$

Interior Harnack Inequality
===========================

We also need the third growth theorem in order to establish the interior Harnack inequality. The first growth theorem tells us if $\mu\rightarrow0^{+}$ then $\beta_{1}\rightarrow0^{+}$. The third growth theorem tells us if we have a nice control of the measure of the set $\left\{ u>0\right\} $ near the bottom, then we can have a more precise estimate. In other words, if we have the similar measure condition for $$Q^{0}:=Q_{\frac{r}{2}}(Y^{0}), Y^{0}=\left(y,s-\frac{3}{4}r^{2}\right).$$ Then if $\mu<1$, then $\beta_{1}<1$. The proof of it is long and technical but independent of the specific structure of the equations. One can find a detailed proof in, for example, [@FS; @KS; @GC]. We just formulate the results here.

\[thm:(Third-Growth-Theorem).\] Let a function $u\in C^{2,1}\left(\overline{Q_{r}}\right)$, where $Q_{r}=Q_{r}(Y)$, $Y=(y,s)\in\mathbb{R}^{n+1}$, $r>0$, and let $-u_{t}+Lu\geq0$ in $Q_{r}$. In addition, we assume $$\left|\left\{ u>0\right\} \cap Q^{0}\right|\leq\mu\left|Q^{0}\right|,\label{eq:71}$$ where $$Q^{0}:=Q_{\frac{r}{2}}(Y^{0}),\,\, Y^{0}=\left(y,s-\frac{3}{4}r^{2}\right)\label{eq:72}$$ and $\mu<1$ is a constant. Then we have $$\mathcal{M}_{\frac{r}{2}}(Y)\leq\beta\mathcal{M}_{r}(Y)\label{eq:73}$$ with a constant $$\beta:=\beta(n,\nu,S,\mu)<1.$$

\[lower\]Let a function $v\in C^{2,1}\left(\overline{Q_{r}}\right)$ be such that $v\geq0$, $-v_{t}+Lv\leq0$ in $Q_{r}$, and $$\left|\left\{ v\geq1\right\} \cap Q^{0}\right|>(1-\mu)\left|Q^{0}\right|.\label{eq:74}$$ Then $$v\geq1-\beta>0\label{eq:75}$$ on $Q_{\frac{r}{2}}$ where $\beta=\beta(n,\nu,\mu,S)<1$ for $\mu<1$.

With the above preparation, we now are ready to establish the interior Harnack inequality.

Suppose $u\in C^{2,1}(Q_{2r}(Y)) \cap C(\overline{Q_{2r}(Y)})$ and $-u_{t}+Lu=0$ in $Q_{2r}(Y)$, $Y=(y,s)\in\mathbb{R}^{n+1}$ and $r>0$. If $u\geq0$, then $$\sup_{Q^{0}}u\leq N\inf_{Q_{r}}u,\label{eq:76}$$ where $N=N(n,\nu,S)$ and $Q^{0}=B_{r}(y)\times(s-3r^{2},s-2r^{2})$.

We will build a non-degenerate intermediate region to get a quantitative relation between two regions we are interested in with the help of three growth theorems.

After rescaling and translating as necessary, we can assume $Y=0$ and $r=1$. Now $Q_{1}=B_{1}(0)\times(-1,0)$, $Q^{0}=B_{1}(0)\times(-3,-2)$. It is easy to see that if we define $d(X):=\sup\left\{ \rho>0:\, Q_{\rho}(X)\subset Q_{2}(0)\right\} $, then $d(X)\geq1$ in $Q^{0}$. Hence, if we consider $Q^{1}:=B_{2}(0)\times(-3,-2)$ we conclude that $$\sup_{Q^{0}}u\leq M:=\sup_{Q^{1}}d^{\gamma}u,\label{eq:77}$$ where $\gamma$ is chosen at the same as the $\gamma$ in Lemma \[lem:quotient\] with $h=\frac{1}{2}$. From the discussion before Theorem \[thm:Lepsilon\], we can find $\exists X_{0}\in\overline{Q^{1}}\backslash\left[\partial_{p}Q^{1}\cap\partial_{p}Q_{2}\right]$ such that $$d^{\gamma}u(X_{0})=M.\label{eq:78}$$ Similarly as above, we define $$\rho=\frac{1}{4}d(X_{0})\in(0,\frac{1}{2}],\label{eq:79}$$ and $$Q_{0}=Q_{\rho}(X_{0})\cap\left\{ u>\frac{1}{2}u(X_{0})\right\} .\label{eq:710}$$ By the above discussion, we obtain $$\left|Q_{0}\right|>\mu_{1}\left|Q_{\rho}(X_{0})\right|$$ for some constant $\mu_{1}=\mu_{1}(n,\nu,S,\gamma)>0$. Now we apply the Corollary \[lower\] with $$v=\frac{2}{u(X_{0})}u,\,\, Q_{r}=Q_{2\rho}(Y_{0}),\,\, Y_{0}=(x_{0},t_{0}+3\rho^{2}),\,\, Q^{0}=Q_{\rho}(X_{0}),\,\,1-\mu=\mu_{1}.$$ Then we have $$u\geq\beta u(X_{0})\label{eq:711}$$ on $Q_{\rho}(Y_{0})$ with $\beta=\beta(n,\nu,S)>0$. Next we apply Lemma \[lem:quotient\] with $$v=u,\,\, r=2,\,\, D_{\rho}=B_{\rho}(x_{0})\times\{t_{0}+2\rho^{2}\}\subset\overline{Q_{\rho}(Y_{0})},$$ and $$D^{0}=B_{1}(0)\times\{\tau\},\,\,\forall\tau\in(-1,0).$$ So we have $$\beta u(X_{0})\leq\inf_{D_{\rho}}u\leq\left(\frac{4}{\rho}\right)^{\gamma}\inf_{Q_{1}(0)}u.\label{eq:712}$$ Finally, with the help of the intermediate region, we conclude that $$\sup_{Q^{0}}u\leq M=d^{\gamma}u(X_{0})=\left(4\rho\right)^{\gamma}u(X_{0})\leq\beta^{-1}4^{2\gamma}\inf_{Q_{1}(0)}u.\label{eq:713}$$ Taking $N=N(n,\nu,S)=\beta^{-1}4^{2\gamma}$ gives the desired result.

It is well-known that it is easy to derive the Hölder continuity of solutions from the Harnack inequality by standard oscillation and iteration arguments.

\[hoelder\] Suppose $u\in W_{p}$ where $p<n+1$ is from Theorem \[variant\], and $u$ is a solution of $-u_{t}+Lu=0$ in $Q_{r}$. Then u is Hölder continuous in $Q_{\frac{r}{2}}$.

Approximation
=============

In all the proofs from above sections, we always assume $u$ is $C^{2,1}$ in stead of $W_{p}$ where $p<n+1$ is from Theorem \[variant\]. In this section, we briefly show we can use an approximation argument to show that all results hold for $u\in W_{p}(Q_{2r})$ in the sense of Definition \[space\] and $p$ is from the variant Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate . Throughout, we assume $$u\geq0,\,\,-u_{t}+Lu=-u_{t}+\sum_{ij}a_{ij}D_{ij}u+\sum_{i}b_{i}D_{i}u=0\label{eq:81}$$ in $Q_{2r}$. We can approximate $a_{ij}$, $b_{i}$ and $u$ by smooth functions $a_{ij}^{\epsilon}\rightarrow a_{ij}$, $b_{i}^{\epsilon}\rightarrow b_{i}$ a.e. as $\epsilon\rightarrow0^{+}$. And $u^{\epsilon}\rightarrow u$ in $W_{p}^{2,1}$ as $\epsilon\rightarrow0^{+}$. Then $$f^{\epsilon}=-u_{t}^{\epsilon}+L^{\epsilon}u^{\epsilon}=-u_{t}^{\epsilon}+\sum_{ij}a_{ij}^{\epsilon}D_{ij}u^{\epsilon}+\sum_{i}b_{i}^{\epsilon}D_{i}u^{\epsilon}\rightarrow0\label{eq:82}$$ in $L_{loc}^{p}(Q_{2r})$ as $\epsilon\rightarrow0^{+}$. We know the existence of solutions to equations with smooth coefficients, therefore we can write $$u^{\epsilon}=v^{\epsilon}+w^{\epsilon},$$ where $$-v_{t}^{\epsilon}+L^{\epsilon}v^{\epsilon}=0$$ in $Q_{2r}$ and $$v^{\epsilon}=u^{\epsilon}$$ on $\partial_{p}Q_{2r}$: $$-w_{t}^{\epsilon}+L^{\epsilon}w^{\epsilon}=f^{\epsilon},$$ $$w^{\epsilon}=0$$ on $\partial_{p}Q_{2r}$. By the variant Aleksandrov-Bakelman-Pucci-Krylov-Tso estimate , we know $w^{\epsilon}\rightarrow0$ in $L^{\infty}$ and $v^{\epsilon}$ satisfies the Harnack inequality. Finally, by an easy limiting argument, $u$ also satisfies the Harnack inequality.

Appendix
========

As we mentioned in the introduction, when drift $b\in L_{x,t}^{n+1}$, we do not expect the solution to have Hölder continuity since $ L_{x,t}^{n+1}$ is supercritical with respect to the parabolic scaling. In this appendix, we present a concrete example. We consider the parabolic equation in $1+1$ dimensions with drift $b\in L_{x,t}^{2}(\mathbb{R}^{2})$, $$u_{t}+b\nabla u-\Delta u=0.\label{eq:81}$$ We define for $t\in[0,1]$ $$b(x,t)=a(t)\begin{cases}
1 & -r(t)\leq x<0\\
-1 & 0<x\leq r(t)\\
0 & x\notin[-r(t),\,0)\cup(0,\, r(t)]
\end{cases}\label{eq:82}$$ and if $t\notin[0,1]$, $b=0$. We set $a(t)=(1-t)^{-\beta}$ and $r(t)=(1-t)^{\alpha}$ where $\beta$ and $\alpha$ to be determined later. First of all, by the integrability condition of $b$, we see $\int_{0}^{1}(1-t)^{\alpha-2\beta}<\infty$, we get $\alpha-2\beta>-1$.

We try to construct an odd function $\phi$ so that we can do a comparison argument. For $0\leq x\leq1$, we define $\ensuremath{\phi(x)=\sin(\pi x/2)}$ and $\phi=1$ for $x>1$. Notice that for $x\in[0,1)$, we have $-\Delta\phi\leq C\phi$ for some constant $C$. In particular, based on our specific choice, we take the constant $C=\left(\frac{2}{\pi}\right)^{2}$ . Finally, we extend this $\phi$ oddly to the whole line. We consider $$v(x,t)=\exp\left[-C\int_{0}^{t}(1-s)^{-2\alpha}\right]\phi\left(x/r(t)\right)\label{eq:83}$$ which requires $-2\alpha>-1$. We try to verify on $\left(0,\, r(t)\right)$, $v$ is a subsolution for $t\in[0,1]$ $$v_{t}=\exp\left[-C\int_{0}^{t}(1-s)^{-2\alpha}\right]\left(x\alpha(1-t)^{-a-1}\phi'-C(1-t)^{-2\alpha}\phi\right)\label{eq:84}$$

$$b\nabla v=-(1-t)^{-\beta}(1-t)^{-\alpha}\phi'\exp\left[-C\int_{0}^{t}(1-s)^{-2\alpha}\right]\label{eq:85}$$

$$-\Delta v=-(1-t)^{-2\alpha}\phi''\exp\left[-C\int_{0}^{t}(1-s)^{-2\alpha}\right]\label{eq:86}$$

By construction, $$-(1-t)^{-2\alpha}\phi''\exp\left[-C\int_{0}^{t}(1-s)^{-2\alpha}\right]-C(1-t)^{-2\alpha}\phi\exp\left[-C\int_{0}^{t}(1-s)^{-2\alpha}\right]\leq0.\label{eq:87}$$ We only need to verify that $$C\left(x\alpha(1-t)^{-a-1}\phi'\right)-(1-t)^{-\beta}(1-t)^{-\alpha}\phi'\exp\left[-C\int_{0}^{t}(1-s)^{-2\alpha}\right]\leq0\label{eq:88}$$ for $x\in\left(0,r(t)\right)$ and $t\in[0,1)$. Since $\phi'$ is nonnegative, it suffices to check $$x\alpha(1-t)^{-a-1}-(1-t)^{-\beta}(1-t)^{-\alpha}\leq0,\,\, x\in[0,r(t)),\,\, t\in[0,1).\label{eq:89}$$ $$x\alpha(1-t)^{-a-1}\leq\alpha(1-t)^{-\alpha-1}(1-t)^{\alpha}=\alpha(1-t)^{-1}.\label{eq:810}$$ We pick $-\beta-\alpha+1<0$. Also $-2\alpha>-1$ and $\alpha-2\beta>-1$. So $\alpha\in(\frac{1}{3},\frac{1}{2})$, we can pick $\alpha=\frac{5}{12}$. Then we pick $\beta=\frac{2}{3}$. This pair satisfies all of our conditions. So with our choice of $\alpha$ and $\beta,$ we can see our $v(x,t)$ is a subsolution to our equation in $[0,r(t))\times[0,1)$.

By the symmetry of our equation and the oddness of $\phi$, $v$ is a subsolution of our equation $[0,r(t))\times[0,1)$ and is a supersolution on $(-r(t),0]\times[0,1)$. If we take $u$ to be the solution of our equation with initial data $v(x,0)$. Then when $t$ approach $1$, we look at the oscillation on the ball $B_{r(t)}(0)$. We have $$2\leq Osc_{B_{r(t)}(0)}v\leq Osc_{B_{r(t)}(0)}u.\label{eq:811}$$ Since $r(t)\rightarrow0$ as $t\rightarrow1$, we conclude that $u$ will have a discontinuity at the origin when $t\rightarrow1$.

Acknowledgment
==============

This work was initiated when I was still at the University of Minnesota Twin Cities. I would like to thank Professor Mikhail Safonov for suggesting this interesting problem to me, and for many motivating discussions. I also want to thank Professor Luis Silvestre for discussions on the supercritical scaling case and some related comments.

[CS13]{} Chen, G. Non-divergence parabolic equations of second order with critical drift in Lebesgue spaces. arXiv:1511.01215

Chen, G.; Safonov, M. V. On second order elliptic and parabolic equations of mixed type. arXiv:1401.0351

Crandall, M. G.; Swiech, A. A note on generalized maximum principles for elliptic and parabolic PDE. In Evolution Equations: Proceedings in Honor of J.A. Goldstein’s 60th Birthday, vol. 234 of Lecture Notes in Pure and Appl. Math., pp. 121–127, Dekker, New York, NY, USA, 2003

Escauriaza, L. $W^{2,n}$ a priori estimates for solutions to fully non-linear equations. Indiana Univ. Math. J. 42 (1993), 413423.

Fabes, E. B.; Safonov, M. V.; Yuan, Y. Behavior near the boundary of positive solutions of second order parabolic equations. II. Trans. Amer. Math. Soc. 351, no. 12 (1999), 4947–4961.

Fabes, E.; Stroock, D. The $L^{p}$ integrability of Green’s functions and fundamental solutions for elliptic and parabolic equations. Duke Math. J. 51 (1984), 997-1016.

Ferretti, E.; Safonov, M. V. Growth theorems and Harnack inequality for second order parabolic equations. In Harmonic Analysis and Boundary Value Problems (Fayetteville, AR, 2000), vol. 277 of Contemp. Math., American Mathematical Society, Providence, RI, USA, 2001, 87–112.

Krylov, N. V. Estimates of the maximum of the solution of a parabolic equation and estimates of the distribution of a semimartingale. (Russian) Mat. Sb. (N.S.) 130(172) (1986), no. 2, 207–221.

Krylov, N. V.; Safonov, M. V. A certain property of solutions of parabolic equations with measurable coefficients. Izvestia Akad. Nauk SSSR, ser. Matem. 44, No. 1 (1980), 161–175 in Russian; English translation in Math. USSR Izvestija, 16, no. 1 (1981), 151–164.

Landis, E. M. Second Order Equations of Elliptic and Parabolic Type. Nauka, Moscow, 1971 in Russian; English transl.: Amer. Math. Soc., Providence, RI, 1997.

Ladyzhenskaya, O. A.; V. A. Solonnikov, V. A.; Ural’tseva, N. N. Linear and Quasi-linear Equations of Parabolic Type. Nauka, Moscow, 1967 in Russian; English transl.: Amer. Math. Soc., Providence, RI, 1968

Lieberman, G. M. Maximum estimates for oblique derivative problems with right-hand side in $L^{p}$, $p<n$. Manuscr. Math., 112, 459–472 (2003).

Lieberman, G. M. Second order parabolic differential equations. World scientific. 2005 [\[]{}Rev. ed.[\]]{}

Moser, J. On Harnack’s theorem for elliptic differential equation. Comm. Pure Appl. Math. 14 (1961), 577–591.

Moser, J. A Harnack inequality for parabolic differential equations. Comm. Pure and Appl. Math. 17 (1964), 101–134; and correction in: Comm. Pure and Appl. Math. 20 (1967), 231–236.

Nash, J. Continuity of solutions of parabolic and elliptic equations. Amer. J. Math. 80 (1958), 931–954

Nazarov, A. I. Interpolation of linear spaces and maximum estimates for solutions of parabolic equations.(Russian) Part. Diff. Eq., Acad. Nauk SSSR, Siberian Dep. Math. Inst., Novosibirsk, 1987, pp. 50–72. MR 0994027.

Nazarov, A. I.; and Ural’tseva, N. N. The Harnack inequality and related properties of solutions of elliptic and parabolic equations with divergence-free lower-order coefficients. Algebra i Analiz 23 (2011), no. 1, 136–168. Milman, M. A note on Gehring?s lemma, Ann. Acad. Sci. Fenn. Math. 21 (1996), 389-398. Safonov, M. V. Harnack inequality for elliptic equations and the Holder property of their solutions. Zap. Nauchn. Sem. Leningrad. Otdel. Mat. Inst. Steklov (LOMI) 96 (1980), 272–287 in Russian; English translation in J. Soviet Math. 21, no. 5 (1983), 851–863.

Safonov, M. V. Estimates near the boundary for solutions to second order parabolic equations. Invited Lecture at the ICM in Berlin, August 18-27, 1998. In: Documenta Matematica, Extra Volume ICM 1998, vol. 1, pp. 637-647.

Safonov, M. V. Non-divergence elliptic equations of second order with unbounded drift. AMS Transl. Series 2. 229 (2010), 211–232.

Silvestre, L.; Vicol, V.; Zlato¨, A. On the loss of continuity for super-critical drift-diffusion equations. Arch. Rational Meeh. Anal. 207(2013), 845–877.

Tso, K. On an Aleksandrov-Bakelman maximum principle for second-order parabolic equations. Comm. Partial Differential Equations 10 (1985), 543–553.
