Let $\rho:\mathbb R^d\times[0,\infty)\to(0,\infty)$ such that $\int \rho_t(x)\,dx=1$ for all $t\geq0\,$, $\rho$ is Holder-continuous (in both variables) and $\rho_t\in W^{1,1}(\mathbb R^d)$ for a.e. $t>0$.

Let $A$ be a positive definite symmetric matrix and $b\in C^2(\mathbb R^d;\mathbb R^d)$. Suppose that $\rho$ solves the following weak version of a parabolic PDE:
\begin{align} \tag{1} & \int_{\mathbb R^d} \varphi_t(x)\,\rho_t(x)\,dx \,-\, \int_{\mathbb R^d} \varphi_0(x)\,\rho_0(x)\,d x \;=\\[5pt] &=\, -\,\lim_{\epsilon\to0}\int_\epsilon^t\!\!\int_{\mathbb R^d} \Big(\partial_s\varphi_s(x)\,\rho_s(x) \,+\, A\,\nabla\!\varphi_s(x)\cdot\nabla\!\rho_s(x) \,+\, b(x)\cdot\nabla\!\varphi_s(x)\,\rho_s(x)\Big)\,d x\,d s \end{align}
for all $t>0$, for all $\varphi\in C(\mathbb R^d\!\times\![0,\infty))$ such that, for example, $\varphi_s(x)=\zeta(s)\,\psi(x)$ with $\zeta\in C^\infty(0,\infty)$ and $\psi\in C^\infty_c(\mathbb R^d)\,$.

**My question:** Is the identity (1) still true when $\varphi$ is replaced by $\log\rho$ ? Namely:
\begin{align} \tag{1'} & \int_{\mathbb R^d} \log\rho_t(x)\,\rho_t(x)\,dx \,-\, \int_{\mathbb R^d} \log\rho_0(x)\,\rho_0(x)\,d x \;=\\[5pt] &=\, -\,\int_0^t\!\!\int_{\mathbb R^d} \Big( A\,\frac{\nabla\!\rho_s(x)}{\rho_s(x)}\cdot\nabla\!\rho_s(x) \,+\, b(x)\cdot\nabla\!\rho_s(x) \Big)\,d x\,d s \quad ?\end{align}
Can we prove (1') by some approximation argument?

I know that:
$$ \int_0^t\!\!\int_{\mathbb R^d} |b(x)|^2\,\rho_s(x) \,d x\,d s \ <\infty \;;\qquad\qquad\int_0^t\!\!\int_{\mathbb R^d} \frac{|\nabla\rho_s(x)|^2}{\rho_s(x)} \,d x\, ds\ <\infty \;.$$

I also know from theory that there exists $\partial_s\rho$ bounded in the *dual space of* $\mathbb H_0^{1,1}(B\times I)$ (*see note at the end) for every $B\times I$ compactly contained in $\mathbb R^d\!\times\!(0,\infty)$ and $||\rho_\epsilon-\rho_0||_{L^1(\mathbb R^d)}\to0\ $ as $\epsilon\to0\,$.
_________________________________________________________________________________

Notice that the integrals in (1') are absolutely convergent.

Regarding the time-derivative term on the r.h.s. of (1), I expect that:
$$ \tag{2} \lim_{\epsilon\to0}\int_\epsilon^t\!\!\int_{\mathbb R^d} \partial_s\rho_s(x)\,d x\,d s \;=\, 0 \;$$

even if I'm not sure how to prove (2). I would take a ball $B_R$ of radius $R$ in $\mathbb R^d$ and say
$$ \int_\epsilon^t\!\!\int_{B_R} \partial_s\rho_s(x)\,d x\,d s \,=\, 
\int_{B_R} \int_\epsilon^t\partial_s\rho_s(x)\,d s\,d x \,=\, 
\int_{B_R} \big(\rho_t(x) - \rho_\epsilon(x)\big) \,d x $$
which in absolute value is bounded by
$$ \int_{B_R} \rho_t(x)\, dx \,-\, \int_{B_R} \rho_0(x) \,d x \,+\, ||\rho_\epsilon-\rho_0||_{L^1(\mathbb R^d)} \,\xrightarrow[\substack{\epsilon\to0,\,R\to\infty}]{}\, 1-1+0 = 0 $$
but I'm not sure this is equivalent to (2).
_____________________________________________________________________________________

(*) $\mathbb H_0^{1,1}(B\times I)$ denotes the space of functions $v:B\times I\to\mathbb R$ such that $v_s\in W^{1,1}_0(B)$ for a.e. $s\in I$ and $\int_I\int_B(|v_s(x)|^p+|\nabla v_s(x)|^p)\,d x\,d s<\infty$.