[guided]The comparison argument we need is not the pure heat-equation maximum principle, because after subtracting a constant barrier the nonlinear term produces a bounded coefficient multiplying the error. We therefore prove the exact linear estimate needed in the sequel.
Let $Q_T := (0,T) \times \Omega$ and define its parabolic boundary by
\begin{align*}
\partial_p Q_T := \bigl(\{0\} \times \overline{\Omega}\bigr) \cup \bigl([0,T] \times \partial\Omega\bigr).
\end{align*}
Suppose
\begin{align*}
w \in C^{1,2}(Q_T) \cap C([0,T] \times \overline{\Omega})
\end{align*}
and $c: Q_T \to \mathbb{R}$ is a bounded real-valued function such that
\begin{align*}
\partial_t w - \Delta w + c w \leq 0 \quad \text{in } Q_T
\end{align*}
and $w \leq 0$ on $\partial_p Q_T$. We want to show that $w$ cannot become positive inside the cylinder.
The coefficient $c$ may have either sign, so we remove it by exponential damping. Choose $M > \|c\|_{L^\infty(Q_T)}$ and define
\begin{align*}
v: [0,T] \times \overline{\Omega} \to \mathbb{R}, \qquad v(t,x) := e^{-Mt}w(t,x).
\end{align*}
The factor $e^{-Mt}$ is positive, so $v$ and $w$ have the same sign. Also $v \leq 0$ on $\partial_p Q_T$.
For $(t,x)\in Q_T$, differentiating the product gives
\begin{align*}
\partial_t v(t,x)-\Delta v(t,x)
=
e^{-Mt}\bigl(\partial_t w(t,x)-\Delta w(t,x)-Mw(t,x)\bigr).
\end{align*}
Using $\partial_t w-\Delta w\leq -cw$ and $w=e^{Mt}v$, this becomes
\begin{align*}
\partial_t v(t,x)-\Delta v(t,x)
\leq
-(M+c(t,x))v(t,x).
\end{align*}
The choice of $M$ ensures $M+c(t,x)>0$ for every $(t,x)\in Q_T$.
Now assume $v$ is positive somewhere. Since $[0,T]\times\overline{\Omega}$ is compact and $v$ is continuous, the maximum
\begin{align*}
A := \max_{[0,T]\times\overline{\Omega}} v
\end{align*}
is positive. The possible difficulty is that the maximum of $v$ might occur at the final time $T$, where the hypotheses do not give a two-sided time derivative. To avoid using any derivative at $T$, choose $\alpha$ with $0<\alpha<A$ and define the first time at which the spatial maximum reaches $\alpha$:
\begin{align*}
t_\alpha := \inf\{t\in[0,T] : \max_{x\in\overline{\Omega}} v(t,x)\geq\alpha\}.
\end{align*}
Because $v\leq0$ at $t=0$ and because $\alpha$ is strictly below the global maximum $A$, continuity gives $t_\alpha\in(0,T)$ and
\begin{align*}
\max_{x\in\overline{\Omega}}v(t_\alpha,x)=\alpha.
\end{align*}
Choose $x_\alpha\in\overline{\Omega}$ such that $v(t_\alpha,x_\alpha)=\alpha$. The lateral boundary condition excludes $x_\alpha\in\partial\Omega$, since $v\leq0$ there and $\alpha>0$; hence $x_\alpha\in\Omega$.
At the point $(t_\alpha,x_\alpha)$ we may use the classical derivatives of $v$, because $t_\alpha<T$ and $x_\alpha\in\Omega$. For every $s<t_\alpha$ and every $x\in\overline{\Omega}$, the definition of the first hitting time gives $v(s,x)<\alpha=v(t_\alpha,x_\alpha)$. Therefore the time derivative satisfies $\partial_t v(t_\alpha,x_\alpha)\geq0$. Also $x_\alpha$ is a spatial maximum point of $v(t_\alpha,\cdot)$ in the [open set](/page/Open%20Set) $\Omega$, so the second-derivative test gives $\Delta v(t_\alpha,x_\alpha)\leq0$. Combining these two facts,
\begin{align*}
\partial_t v(t_\alpha,x_\alpha)-\Delta v(t_\alpha,x_\alpha)\geq0.
\end{align*}
But the differential inequality gives
\begin{align*}
\partial_t v(t_\alpha,x_\alpha)-\Delta v(t_\alpha,x_\alpha)
\leq
-(M+c(t_\alpha,x_\alpha))v(t_\alpha,x_\alpha)<0,
\end{align*}
because $M+c(t_\alpha,x_\alpha)>0$ and $v(t_\alpha,x_\alpha)=\alpha>0$. This contradiction proves $v\leq0$, and the positivity of $e^{-Mt}$ gives $w\leq0$.[/guided]