Consider the following evolution equation
$$u_t=\Delta u$$ in a bounded and regular open subset $\Omega$ of $\mathbb{R}^N$, with smooth initial conditions $u_0\geq 0$ and homogeneous Dirichlet boundary conditions.
It is known that this equation has a smooth global solution $u$. By differentiating $E(t):=\int_\Omega u^2 dx$, using integration parts and Poincaré's inequality we can prove that
$$|u(t)|_{L^2(\Omega)}\leq e^{ -\frac{1}{c^2}t}|u(0)|_{L^2(\Omega)}$$ where $c$ is the Poincaré constant. This means that the solution decays exponentially in $L^2(\Omega)$.
My question is: Does this decay result holds in $L^p(\Omega)$ for $p>2$ or maybe even $L^\infty(\Omega)$? The fact that $p=2$ is very important here because after the integration by parts we get the term $\int_\Omega |\nabla u|^2dx$ which we know how to estimate in terms of $\int_\Omega u^2 dx$ via Poincaré's inequality.