For $1 \leq p < \infty$ and initial-data $u _ 0 \in L ^p ( \mathbb{R} ^d)$ due to P. Li (Uniqueness of $L^1$ solutions for the Laplace equation and the heat equation on Riemannian manifolds) there exists a unique solution $u$ of the heat equation $u_t - \Delta u = 0$ in $(0, \infty) \times \mathbb{ R }^d$ which satisfies \begin{equation*} u \in C ( [0 , \infty ) ; L^p ( \mathbb{ R }^d) ) \quad \text{and} \quad u ( 0 , \cdot ) = u_0. \end{equation*} This solution is given by $u ( t ) = e ^{t \Delta} u_0$, the convolution with the fundamental solution.
My question is, whether there is known literature for the uniqueness in the case of non-negative $u_0 \in L^\infty ( \mathbb{ R }^d )$ for the heat equation in $(0, \infty) \times \mathbb{ R }^d$, where we replace $u \in C ( [0 , \infty ) ; L^p ( \mathbb{ R }^d) )$ by $u \in L^\infty ( (0, \infty) \times \mathbb{ R }^d)$ and $u(t,x) \to u_0(x)$ for almost every $x \in \mathbb{ R }^d$ as $t \to 0$. Is weak star convergency of $u(t, \cdot)$ to $u_0$ in $L^\infty(\mathbb{R}^d)$ as $t \to 0$ sufficient?