I've been following Reif's Fundamentals of Statistical and Thermal Physics; there I came before the derivation of Liouville's theorem:
There I couldn't understood few things.
I could conceive the change in the number of systems in $\mathrm dt$ is given by $$\frac{\partial \rho}{\partial t}\; \mathrm dt\; (\mathrm d q_1,\mathrm dq_2, \ldots,\mathrm d q_f; \mathrm d p_1, \mathrm dp_2,\ldots, \mathrm d p_f)$$ where
- $$\rho(q_1,q_2,\ldots, q_f; p_1,p_2,\ldots, p_f ; t)\;\mathrm d q_1,\mathrm dq_2, \ldots,\mathrm d q_f; \mathrm d p_1, \mathrm dp_2,\ldots, \mathrm d p_f = \textrm{no of systems in the ensemble at $t$ in the phase-space volume}\;(\mathrm d q_1,\mathrm dq_2, \ldots,\mathrm d q_f; \mathrm d p_1, \mathrm dp_2,\ldots, \mathrm d p_f)$$
But then, I couldn't understand why the number of systems 'entering this volume in time $dt$ through the face $q_1$= constant' is given by the quantity $\rho(\dot{q_1}\mathrm dt, \mathrm dq_2, \ldots,\mathrm dp_f )\;.$
My questions are:
$\bullet$ How does $\rho(q_1,q_2,\ldots, q_f; p_1,p_2,\ldots, p_f ; t)((\dot{q_1}\mathrm dt)(\mathrm dq_2, \ldots,\mathrm dp_f ))$ represent the number of systems that would enter the volume in time-interval $\mathrm d t\;?$
$\bullet$ How does the evaluation of $\dot q_i$ at $q_1+\mathrm dq_1$ yield $\dot q_i +\dfrac{\partial \dot q_i}{\partial q_1}\mathrm dq_1 \;?$