# Elementary Principles in Statistical Mechanics/Chapter I

ELEMENTARY PRINCIPLES IN STATISTICAL MECHANICS
CHAPTER I.
GENERAL NOTIONS. THE PRINCIPLE OF CONSERVATION OF EXTENSION-IN-PHASE.
We shall use Hamilton's form of the equations of motion for a system of $n$ degrees of freedom, writing $q_{1},\ldots q_{n}$ for the (generalized) coördinates, ${\dot {q}}_{1},\ldots {\dot {q}}_{n}$ for the (generalized) velocities, and
| $F_{1}\,dq_{1}+F_{2}\,dq_{2}\ldots +F_{n}\,dq_{n}$ | (1) |
$F_{1},\ldots F_{n}$
$p_{1}\ldots p_{n}$
| $p_{1}={\frac {d\epsilon _{p}}{d{\dot {q}}_{1}}},\quad p_{2}={\frac {d\epsilon _{p}}{d{\dot {q}}_{2}}},\quad \mathrm {etc.,}$ | (2) |
$\epsilon _{p}$
$\epsilon _{p}$
${\dot {q}}$
$q$
$d\epsilon _{p}/dq_{1}$
$p$
$q$
We have then
| ${\dot {q}}_{1}={\frac {d\epsilon _{p}}{dp_{1}}},\quad {\dot {p}}_{1}=-{\frac {d\epsilon _{p}}{dq_{1}}}+F_{1},\quad \mathrm {etc.}$ | (3) |
These equations will hold for any forces whatever. If the forces are conservative, in other words, if the expression (1) is an exact differential, we may set
| $F_{1}=-{\frac {d\epsilon _{q}}{dq_{1}}},\quad F_{2}=-{\frac {d\epsilon _{q}}{dq_{2}}},\quad \mathrm {etc.}$ | (4) |
$\epsilon _{q}$
$\epsilon$
| $\epsilon =\epsilon _{p}+\epsilon _{q}$ | (5) |
| ${\dot {q}}_{1}={\frac {d\epsilon }{dp_{1}}},\quad {\dot {p}}_{1}=-{\frac {d\epsilon }{dq_{1}}},\quad \mathrm {etc.}$ | (6) |
The potential energy ($\epsilon _{q}$) may depend on other variables beside the coördinates $q_{1},\ldots q_{n}$. We shall often suppose it to depend in part on coördinates of external bodies, which we shall denote by $a_{1}$, $a_{2}$, etc. We shall then have for the complete value of the differential of the potential energy
| $d\epsilon _{q}=-F_{1}\,dq_{1}\ldots -F_{n}\,dq_{n}-A_{1}\,da_{1}-A_{2}\,da_{2}-\mathrm {etc.} ,$ | (7) |
$A_{1}$
$A_{2}$
$\epsilon$
| ${\begin{aligned}d\epsilon ={\dot {q}}_{1}\,dp_{1}\ldots &+{\dot {q}}_{n}\,dp_{n}-{\dot {p}}_{1}\,dq_{1}\ldots \\&-{\dot {p}}_{n}\,dq_{n}-A_{1}\,da_{1}-A_{2}\,da_{2}-\mathrm {etc.} \end{aligned}}$ | (8) |
It will be observed that the kinetic energy ($\epsilon _{p}$) in the most general case is a quadratic function of the $p$'s (or ${\dot {q}}$'s) ​involving also the $q$'s but not the $a$'s ; that the potential energy, when it exists, is function of the $q$'s and $a$'s ; and that the total energy, when it exists, is function of the $p$'s (or ${\dot {q}}$s), the $q$'s, and the $a$'s. In expressions like $d\epsilon /dq_{1}$ the $p$'s, and not the ${\dot {q}}$'s, are to be taken as independent variables, as has already been stated with respect to the kinetic energy.
Let us imagine a great number of independent systems, identical in nature, but differing in phase, that is, in their condition with respect to configuration and velocity. The forces are supposed to be determined for every system by the same law, being functions of the coördinates of the system $q_{1},\ldots q_{n}$, either alone or with the coördinates $a_{1}$, $a_{2}$, etc. of certain external bodies. It is not necessary that they should be derivable from a force-function. The external coördinates $a_{1}$, $a_{2}$, etc. may vary with the time, but at any given time have fixed values. In this they differ from the internal coördinates $q_{1},\ldots q_{n}$, which at the same time have different values in the different systems considered.
Let us especially consider the number of systems which at a given instant fall within specified limits of phase, viz., those for which
| ${\begin{array}{ccc}p_{1}'<p_{1}<p_{1}'',&\qquad &q_{1}'<q_{1}<q_{1}''\\p_{2}'<p_{2}<p_{2}'',&\qquad &q_{2}'<q_{2}<q_{2}''\\\ldots &&\ldots \\p_{n}'<p_{n}<p_{n}'',&\qquad &q_{n}'<q_{n}<q_{n}''\end{array}}$ | (9) |
$p_{1}''-p_{1}'$
$q_{1}''-q_{1}'$
| $D(p_{1}''-p_{1}')\ldots (p_{n}''-p_{n}')(q_{1}''-q_{1}')\ldots (q_{n}''-q_{n}')$ | (10) |
| $D\,dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n},$ | (11) |
$D$
$p$
$q$
$t$
If we regard all possible phases as forming a sort of extension of $2n$ dimensions, we may regard the product of differentials in (11) as expressing an element of this extension, and $D$ as expressing the density of the systems in that element. We shall call the product
| $dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}$ | (12) |
$D$
It is evident that the changes which take place in the density of the systems in any given element of extension-in-phase will depend on the dynamical nature of the systems and their distribution in phase at the time considered.
In the case of conservative systems, with which we shall be principally concerned, their dynamical nature is completely determined by the function which expresses the energy ($\epsilon$) in terms of the $p$'s, $q$'s, and $a$'s (a function supposed identical for all the systems); in the more general case which we are considering, the dynamical nature of the systems is determined by the functions which express the kinetic energy ($\epsilon _{p}$) in terms of the $p$'s and $q$'s, and the forces in terms of the $q$'s and $a$'s. The distribution in phase is expressed for the time considered by $D$ as function of the $p$'s and $q$'s. To find the value of $dD/dt$ for the specified element of extension-in-phase, we observe that the number of systems within the limits can only be varied by systems passing the limits, which may take place in $4n$ different ways, viz., by the $p_{1}$ of a system passing the limit $p_{1}'$, or the limit $p_{1}''$, or by the $q_{1}$ of a system passing the limit $q_{1}'$ or the limit $q_{1}''$, etc. Let us consider these cases separately.
​In the first place, let us consider the number of systems which in the time $dt$ pass into or out of the specified element by $p_{1}$ passing the limit $p_{1}'$. It will be convenient, and it is evidently allowable, to suppose $dt$ so small that the quantities ${\dot {p}}_{1}\,dt$, ${\dot {q}}_{1}\,dt$, etc., which represent the increments of $p_{1}$, $q_{1}$, etc., in the time $dt$ shall be infinitely small in comparison with the infinitesimal differences $p_{1}''-p_{1}'$, $q_{1}''-q_{1}'$, etc., which determine the magnitude of the element of extension-in-phase. The systems for which $p_{1}$ passes the limit $p_{1}'$ in the interval $dt$ are those for which at the commencement of this interval the value of $p_{1}$ lies between $p_{1}'$ and $p_{1}'-{\dot {p}}_{1}\,dt$, as is evident if we consider separately the cases in which ${\dot {p}}_{1}$ is positive and negative. Those systems for which $p_{1}$ lies between these limits, and the other $p$'s and $q$'s between the limits specified in (9), will therefore pass into or out of the element considered according as ${\dot {p}}$ is positive or negative, unless indeed they also pass some other limit specified in (9) during the same interval of time. But the number which pass any two of these limits will be represented by an expression containing the square of $dt$ as a factor, and is evidently negligible, when $dt$ is sufficiently small, compared with the number which we are seeking to evaluate, and which (with neglect of terms containing $dt^{2}$) may be found by substituting ${\dot {p}}_{1}\,dt$ for $p_{1}''-p_{1}'$ in (10) or for $dp_{1}$ in (11).
The expression
| $D\,{\dot {p}}_{1}\,dt\,dp_{2}\ldots dp_{n}\,dq_{1}\ldots dq_{n}$ | (13) |
$p_{1}'$
$D$
${\dot {p}}$
$p_{1}''$
$p_{1}'$
$p_{1}''$
| ${\frac {d(D\,{\dot {p}}_{1})}{dp_{1}}}\,dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}\,dt$ | (14) |
$p_{1}'$
$p_{1}''$
The decrease in the number of systems within the limits due to systems passing the limits $q_{1}'$ and $q_{1}''$ may be found in the same way. This will give
| ${\bigg{\frac {d(D\,{\dot {p}}_{1})}{dp_{1}}}+{\frac {d(D\,{\dot {q}}_{1})}{dq_{1}}}{\bigg )}\,dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}\,dt$ | (15) |
$p_{1}'$
$p_{1}''$
$q_{1}'$
$q_{1}''$
| ${\frac {d{\dot {p}}_{1}}{dp_{1}}}+{\frac {d{\dot {q}}_{1}}{dq_{1}}}=0,$ | (16) |
| ${\bigg{\frac {dD}{dp_{1}}}{\dot {p}}_{1}+{\frac {dD}{dq_{1}}}{\dot {q}}_{1}{\bigg )}\,dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}\,dt$ | (17) |
If we prefix $\Sigma$ to denote summation relative to the suffixes $1\ldots n$, we get the total decrease in the number of systems within the limits in the time $dt$. That is,
| ${\begin{aligned}\Sigma {\Bigg{\frac {dD}{dp_{1}}}{\dot {p}}_{1}+{\frac {dD}{dq_{1}}}{\dot {q}}_{1}{\Bigg )}\,dp_{1}\ldots &dp_{n}\,dq_{1}\ldots dq_{n}\,dt=\\&-dD\,dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n},\end{aligned}}$ | (18) |
| ${\Bigg{\frac {dD}{dt}}{\Bigg )}_{p,q}=-\Sigma {\Bigg{\frac {dD}{dp_{1}}}{\dot {p}}_{1}+{\frac {dD}{dq_{1}}}{\dot {q}}_{1}{\Bigg )},$ | (19) |
$p$
$q$
| $\Sigma {\Bigg{\frac {dD}{dp_{1}}}{\dot {p}}_{1}+{\frac {dD}{dq_{1}}}{\dot {q}}_{1}{\Bigg )}=0.$ | (20) |
$p$
$q$
$(dD/dt)_{p,q}$
${\dot {p}}$
If we write equation (19) in the form
| ${\Bigg{\frac {dD}{dt}}{\Bigg )}_{p,q}dt+\Sigma {\Bigg{\frac {dD}{dp_{1}}}{\dot {p}}_{1}dt+{\frac {dD}{dq_{1}}}{\dot {q}}_{1}dt{\Bigg )}=0$ | (21) |
$D$
$t$
$p_{1},\ldots p_{n}$
$q_{1}\ldots q_{n}$
$D$
$t$
$p$
$q$
$D$
$p$
$q$
${\dot {p}}_{1}\,dt$
${\dot {q}}_{1}\,dt$
$p$
$q$
$dt$
$D$
In an ensemble of mechanical systems identical in nature and subject to forces determined by identical laws, but distributed in phase in any continuous manner, the density-in-phase is constant in time for the varying phases of a moving system; provided, that the forces of a system are functions of its coördinates, either alone or with the time.
This may be called the principle of conservation of density-in-phase. It may also be written
| ${\Bigg{\frac {dD}{dt}}{\Bigg )}_{a,\ldots h}=0,$ | (22) |
$a,\ldots h$
We may give to this principle a slightly different expression. Let us call the value of the integral
| $\int \ldots \int dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}$ | (23) |
When the phases bounding an extension-in-phase vary in the course of time according to the dynamical laws of a system subject to forces which are functions of the coördinates either alone or with the time, the value of the extension-in-phase thus bounded remains constant. In this form the principle may be called the principle of conservation of extension-in-phase. In some respects this may be regarded as the most simple statement of the principle, since it contains no explicit reference to an ensemble of systems.
Since any extension-in-phase may be divided into infinitesimal portions, it is only necessary to prove the principle for an infinitely small extension. The number of systems of an ensemble which fall within the extension will be represented by the integral
* $\int \ldots \int D\,dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}.$
$D$
* $D\int \ldots \int dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}$
$D$
* $\int \ldots \int dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n},$
Since the system of coördinates employed in the foregoing discussion is entirely arbitrary, the values of the coördinates relating to any configuration and its immediate vicinity do not impose any restriction upon the values relating to other configurations. The fact that the quantity which we have called density-in-phase is constant in time for any given system, implies therefore that its value is independent of the coördinates which are used in its evaluation. For let the density-in-phase as evaluated for the same time and phase by one system of coördinates be $D_{1}'$, and by another system $D_{2}'$. A system which at that time has that phase will at another time have another phase. Let the density as calculated for this second time and phase by a third system of coördinates be $D_{3}''$. Now we may imagine a system of coördinates which at and near the first configuration will coincide with the first system of coördinates, and at and near the second configuration will coincide with the third system of coördinates. This will give $D_{1}'=D_{3}''$. Again we may imagine a system of coördinates which at and near the first configuration will coincide with the second system of coördinates, and at and near the ​second configuration will coincide with the third system of coördinates. This will give $D_{2}'=D_{3}''$. We have therefore $D_{1}'=D_{2}'$.
It follows, or it may be proved in the same way, that the value of an extension-in-phase is independent of the system of coördinates which is used in its evaluation. This may easily be verified directly. If $q_{1},\ldots q_{n}$, $Q_{1},\ldots Q_{n}$ are two systems of coördinates, and $p_{1},\ldots p_{n}$, $P_{1},\ldots P_{n}$ the corresponding momenta, we have to prove that
| $\int \ldots \int dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}=\int \ldots \int dP_{1}\ldots dP_{n}\,dQ_{1}\ldots dQ_{n},$ | (24) |
| ${\frac {d(P_{1},\ldots P_{n},Q_{1},\ldots Q_{n})}{d(p_{1},\ldots p_{n},q_{1},\ldots q_{n})}}=1,$ | (25) |
$dQ/dp$
| ${\frac {d(P_{1},\ldots P_{n})}{d(p_{1},\ldots p_{n})}}{\frac {d(Q_{1},\ldots Q_{n})}{d(q_{1},\ldots q_{n})}}=1.$ | (26) |
${\dot {Q}}$
${\dot {q}}$
$p$
$q$
$d{\dot {Q}}_{r}/dp_{y}$
$q$
| ${\begin{aligned}{\frac {dP_{x}}{dp_{y}}}={\frac {d}{dp_{y}}}{\frac {d\epsilon _{p}}{d{\dot {Q}}_{x}}}&=\sum _{r=1}^{r=n}{\bigg{\frac {d^{2}\epsilon _{p}}{d{\dot {Q}}_{r}d{\dot {Q}}_{x}}}{\frac {d{\dot {Q}}_{r}}{dp_{y}}}{\bigg )}=\\&{\frac {d}{d{\dot {Q}}_{x}}}\sum _{r=1}^{r=n}{\bigg{\frac {d\epsilon _{p}}{d{\dot {Q}}_{r}}}{\frac {d{\dot {Q}}_{r}}{dp_{y}}}{\bigg )}={\frac {d}{d{\dot {Q}}_{x}}}{\frac {d\epsilon _{p}}{dp_{y}}}={\frac {d{\dot {q}}_{y}}{d{\dot {Q}}_{x}}}.\end{aligned}}$ | (27) |
* ${\dot {q}}_{y}=\sum _{r=1}^{r=n}{\bigg{\frac {dq_{y}}{dQ_{r}}}{\dot {Q}}_{r}{\bigg )}$,
| ${\frac {d{\dot {q}}_{y}}{d{\dot {Q}}_{x}}}={\frac {dq_{y}}{dQ_{x}}}.$ | (28) |
| ${\frac {d(P_{1},\ldots P_{n})}{d(p_{1},\ldots p_{n})}}={\frac {d({\dot {q}}_{1},\ldots {\dot {q}}_{n})}{d({\dot {Q}}_{1},\ldots {\dot {Q}}_{n})}}={\frac {d(q_{1},\ldots q_{n})}{d(Q_{1},\ldots Q_{n})}}.$ | (29) |
| ${\frac {d({\dot {q}}_{1},\ldots {\dot {q}}_{n})}{d(Q_{1},\ldots Q_{n})}}{\frac {d(Q_{1},\ldots Q_{n})}{d({\dot {q}}_{1},\ldots {\dot {q}}_{n})}}=1,$ | (30) |
The numerical value of an extension-in-phase will however depend on the units in which we measure energy and time. For a product of the form $dp\,dq$ has the dimensions of energy multiplied by time, as appears from equation (2), by which the momenta are defined. Hence an extension-in-phase has the dimensions of the $n$th power of the product of energy and time. In other words, it has the dimensions of the $n$th power of action, as the term is used in the `principle of Least Action.'
If we distinguish by accents the values of the momenta and coördinates which belong to a time $t'$, the unaccented letters relating to the time $t$, the principle of the conservation of extension-in-phase may be written
| $\int \ldots \int dp_{1}\ldots dp_{n}\,dq_{1}\ldots dq_{n}=\int \ldots \int dp_{1}'\ldots dp_{n}'\,dq_{1}'\ldots dq_{n}',$ | (31) |
| $\int \ldots \int dp_{1}\ldots dq_{n}=\int \ldots \int dp_{1}'\ldots dq_{n}',$ | (32) |
$t$
$t'$
* $\int \ldots \int dp_{1}\ldots dq_{n}=\int \ldots \int {\frac {d(p_{1}\ldots q_{n})}{d(p_{1}'\ldots q_{n}')}}dp_{1}'\ldots dq_{n}'$
| ${\frac {d(p_{1},\ldots q_{n})}{d(p_{1}'\ldots q_{n}')}}=1.$ | (33) |
* ${\frac {d(p_{1},\ldots q_{n})}{d(p_{1}',\ldots q_{n}')}}={\frac {d(p_{1},\ldots q_{n})}{d(p_{1}'',\ldots q_{n}'')}}{\frac {d(p_{1}'',\ldots q_{n}'')}{d(p_{1}',\ldots q_{n}')}},$
$t''$
$t$
$t'$
$t''$
| ${\frac {d}{dt}}{\frac {d(p_{1},\ldots q_{n})}{d(p_{1}',\ldots q_{n}')}}={\frac {d(p_{1}'',\ldots q_{n}'')}{d(p_{1}',\ldots q_{n}')}}{\frac {d}{dt}}{\frac {d(p_{1},\ldots q_{n})}{d(p_{1}'',\ldots q_{n}'')}}.$ | (34) |
$t''$
$t''$
$t$
* ${\frac {d(p_{1},\ldots q_{n})}{d(p_{1}'',\ldots q_{n}'')}}$
* ${\frac {d}{dt}}{\frac {d(p_{1},\ldots q_{n})}{d(p_{1}'',\ldots q_{n}'')}}={\frac {d{\dot {p}}_{1}}{dp_{1}''}}\ldots +{\frac {d{\dot {p}}_{n}}{dp_{n}''}}+{\frac {d{\dot {q}}_{1}}{dq_{1}''}}\ldots +{\frac {d{\dot {q}}_{n}}{dq_{n}''}}.$
$t=t''$
* ${\frac {d}{dt}}{\frac {d(p_{1},\ldots q_{n})}{d(p_{1}'',\ldots q_{n}'')}}=0,$
* ${\frac {d}{dt}}{\frac {d(p_{1},\ldots q_{n})}{d(p_{1}',\ldots q_{n}')}}=0.$
$t=t'$
Again, if we write $a,\ldots h$ for a system of $2n$ arbitrary constants of the integral equations of motion, $p_{1}$, $q_{1}$, etc. will be functions of $a,\ldots h$, and $t$, and we may express an extension-in-phase in the form
| $\int \ldots \int {\frac {d(p_{1},\ldots q_{n})}{d(a,\ldots h)}}da\ldots dh.$ | (35) |
$a,\ldots h$
| ${\frac {d}{dt}}{\frac {d(p_{1},\ldots q_{n})}{d(a,\ldots h)}}=0.$ | (36) |
* ${\frac {d(p_{1},\ldots q_{n})}{d(a,\ldots h)}}={\frac {d(p_{1},\ldots q_{n})}{d(p_{1}',\ldots q_{n}')}}{\frac {d(p_{1}',\ldots q_{n}')}{d(a,\ldots h)}}$
Since the coördinates and momenta are functions of $a,\ldots h$, and $t$, the determinant in (36) must be a function of the same variables, and since it does not vary with the time, it must be a function of $a,\ldots h$ alone. We have therefore
| ${\frac {d(p_{1},\ldots q_{n})}{d(a,\ldots h)}}=\operatorname {func.} (a,\ldots h).$ | (37) |
​It is the relative numbers of systems which fall within different limits, rather than the absolute numbers, with which we are most concerned. It is indeed only with regard to relative numbers that such discussions as the preceding will apply with literal precision, since the nature of our reasoning implies that the number of systems in the smallest element of space which we consider is very great. This is evidently inconsistent with a finite value of the total number of systems, or of the density-in-phase. Now if the value of $D$ is infinite, we cannot speak of any definite number of systems within any finite limits, since all such numbers are infinite. But the ratios of these infinite numbers may be perfectly definite. If we write $N$ for the total number of systems, and set
| $P={\frac {D}{N}},$ | (38) |
$P$
$N$
$D$
| $\int \ldots \int P\,dp_{1}\ldots dq_{n}$ | (39) |
| $P\,dp_{1}\ldots dq_{n}$ | (40) |
$dp_{1}\ldots dq_{n}$
$P$
$\eta$
If we substitute $NP$ and $Ne^{\eta }$ for $D$ in equation (19), we get
| ${\bigg{\frac {dP}{dt}}{\bigg )}_{p,q}=-\Sigma {\bigg{\frac {dP}{dp_{1}}}{\dot {p}}_{1}+{\frac {dP}{dq_{1}}}d{\dot {q}}_{1}{\bigg )},$ | (41) |
| ${\bigg{\frac {d\eta }{dt}}{\bigg )}_{p,q}=-\Sigma {\bigg{\frac {d\eta }{dp_{1}}}{\dot {p}}_{1}+{\frac {d\eta }{dq_{1}}}d{\dot {q}}_{1}{\bigg )}.$ | (42) |
The same substitutions in (22) give
| ${\bigg{\frac {dP}{dt}}{\bigg )}_{a,\ldots h}=0,$ | (43) |
| ${\bigg{\frac {d\eta }{dt}}{\bigg )}_{a,\ldots h}=0.$ | (44) |
$P$
$\eta$
$D$
The application of this principle is not limited to cases in which there is a formal and explicit reference to an ensemble of systems. Yet the conception of such an ensemble may serve to give precision to notions of probability. It is in fact customary in the discussion of probabilities to describe anything which is imperfectly known as something taken at random from a great number of things which are completely described. But if we prefer to avoid any reference to an ensemble of systems, we may observe that the probability that the phase of a system falls within certain limits at a certain time, is equal to the probability that at some other time the phase will fall within the limits formed by phases corresponding to the first. For either occurrence necessitates the other. That is, if we write $P'$ for the coefficient of probability of the phase $p_{1}',\ldots q_{n}'$ at the time $t'$, and $P''$ for that of the phase $p_{1}'',\ldots q_{n}''$ at the time $t''$, ​
| $\int \ldots \int P'\,dq_{1}'\ldots dq_{n}'=\int \ldots \int P''\,dp_{1}''\ldots dq_{n}'',$ | (45) |
$P'$
$P''$
* $P'\int \ldots \int \,dp_{1}'\ldots dq_{n}''=P''\int \ldots \int \,dp_{1}''\ldots dq_{n}''.$
* $P''=P'.$
With reference to an important class of cases this principle may be enunciated as follows.
When the differential equations of motion are exactly known, but the constants of the integral equations imperfectly determined, the coefficient of probability of any phase at any time is equal to the coefficient of probability of the corresponding phase at any other time. By corresponding phases are meant those which are calculated for different times from the same values of the arbitrary constants of the integral equations.
Since the sum of the probabilities of all possible cases is necessarily unity, it is evident that we must have
| $\mathop {\int \ldots \int } _{\rm {phases}}^{\rm {all}}P\,dp_{1}\ldots dq_{n}=1,$ | (46) |
* $N=\mathop {\int \ldots \int } _{\rm {phases}}^{\rm {all}}D\,dp_{1}\ldots dq_{n},$
$N$
​The values of the coefficient and index of probability of phase, like that of the density-in-phase, are independent of the system of coördinates which is employed to express the distribution in phase of a given ensemble.
In dimensions, the coefficient of probability is the reciprocal of an extension-in-phase, that is, the reciprocal of the $n$th power of the product of time and energy. The index of probability is therefore affected by an additive constant when we change our units of time and energy. If the unit of time is multiplied by $c_{t}$ and the unit of energy is multiplied by $c_{e}$, all indices of probability relating to systems of $n$ degrees of freedom will be increased by the addition of
| $n\log c_{t}+n\log c_{e}.$ | (47) |
* ↑ The use of the momenta instead of the velocities as independent variables is the characteristic of Hamilton's method which gives his equations of motion their remarkable degree of simplicity. We shall find that the fundamental notions of statistical mechanics are most easily defined, and are expressed in the most simple form, when the momenta with the coördinates are used to describe the state of a system.
* ↑ It will be observed, that although we call $\epsilon _{q}$ the potential energy of the system which we are considering, it is really so defined as to include that energy which might be described as mutual to that system and external bodies.
* ↑ In strictness, a finite number of systems cannot be distributed continuously in phase. But by increasing indefinitely the number of systems, we may approximate to a continuous law of distribution, such as is here described. To avoid tedious circumlocution, language like the above may be allowed, although wanting in precision of expression, when the sense in which it is to be taken appears sufficiently clear.
* ↑ The condition that the forces $F_{1},\ldots F_{n}$ are functions of $q_{1},\ldots q_{n}$ and $a_{1}$, $a_{2}$, etc., which last are functions of the time, is analytically equivalent to the condition that $F_{1},\ldots F_{n}$ are functions of $q_{1},\ldots q_{n}$ and the time. Explicit mention of the external coördinates, $a_{1}$, $a_{2}$, etc., has been made in the preceding pages, because our purpose will require us hereafter to consider these coördinates and the connected forces, $A_{1}$, $A_{2}$, etc., which represent the action of the systems on external bodies.
* ↑ If we regard a phase as represented by a point in space of $2n$ dimensions, the changes which take place in the course of time in our ensemble of systems will be represented by a current in such space. This current will be steady so long as the external coördinates are not varied. In any case the current will satisfy a law which in its various expressions is analogous to the hydrodynamic law which may be expressed by the phrases conservation of volumes or conservation of density about a moving point, or by the equation ${\frac {d{\dot {x}}}{dx}}+{\frac {d{\dot {y}}}{dy}}+{\frac {d{\dot {z}}}{dz}}=0.$ The analogue in statistical mechanics of this equation, viz., ${\frac {d{\dot {p}}_{1}}{dp_{1}}}+{\frac {d{\dot {q}}_{1}}{dq_{1}}}+{\frac {d{\dot {p}}_{2}}{dp_{2}}}+{\frac {d{\dot {q}}_{2}}{dq_{2}}}+\mathrm {etc.} =0$ may be derived directly from equations (3) or (6), and may suggest such theorems as have been enunciated, if indeed it is not regarded as making them intuitively evident. The somewhat lengthy demonstrations given above will at least serve to give precision to the notions involved, and familiarity with their use.
* ↑ The form of the equation ${\frac {d}{dp_{y}}}{\frac {d\epsilon _{p}}{d{\dot {Q}}_{x}}}={\frac {d}{d{\dot {Q}}_{x}}}{\frac {d\epsilon _{p}}{dp_{y}}}$ in (27) reminds us of the fundamental identity in the differential calculus relating to the order of differentiation with respect to independent variables. But it will be observed that here the variables ${\dot {Q}}_{x}$ and $p_{y}$ are not independent and that the proof depends on the linear relation between the ${\dot {Q}}$'s and the $p$'s.