\rho(\vec{r},t)