01:24
@User198 You are so close that Tobias did not even want to stress the small difference. You have a weighted sum with classical probabilities, so the equations must be $$\begin{align}\tag1\rho(t)&=\sum_jp_j\left|\psi_j(t)\right>\!\left<\psi_j(t)\right|\\\tag2&=\sum_jp_jU(t-t_0)\left|\psi_j(t_0)\right>\!\left<\psi_j(t_0)\right|U^\dagger(t,t_0)\\\tag3&=U(t,t_0)\left(\sum_jp_j\left|\psi_j(t_0)\right>\!\left<\psi_j(t_0)\right|\right)U^\dagger(t-t_0)\end {align}$$
$$\tag4\therefore\qquad\rho(t)=U(t,t_0)\rho(t_0)U^\dagger(t,t_0)$$ as wanted
Hence, if you postulated quantum dynamics for wavefunctions, then you will automatically obtain the correct dynamics for density operators by this derivation.