Editing Prior probability (section)

==Prior probability in statistical mechanics==
While in Bayesian statistics the prior probability is used to represent initial beliefs about an uncertain parameter, in [[statistical mechanics]] the a priori probability is used to describe the initial state of a system.<ref>{{cite book |last=Iba |first=Y. |year=1989 |chapter=Bayesian Statistics and Statistical Mechanics |editor-last=Takayama |editor-first=H. |title=Cooperative Dynamics in Complex Physical Systems |series=Springer Series in Synergetics |volume=43 |pages=235–236 |publisher=Springer |location=Berlin |doi=10.1007/978-3-642-74554-6_60 |isbn=978-3-642-74556-0 }}</ref> The classical version is defined as the ratio of the number of [[elementary event]]s (e.g., the number of times a die is thrown) to the total number of events—and these considered purely deductively, i.e., without any experimenting. In the case of the die if we look at it on the table without throwing it, each elementary event is reasoned deductively to have the same probability—thus the probability of each outcome of an imaginary throwing of the (perfect) die or simply by counting the number of faces is 1/6. Each face of the die appears with equal probability—probability being a measure defined for each elementary event. The result is different if we throw the die twenty times and ask how many times (out of 20) the number 6 appears on the upper face. In this case time comes into play and we have a different type of probability depending on time or the number of times the die is thrown. On the other hand, the a priori probability is independent of time—you can look at the die on the table as long as you like without touching it and you deduce the probability for the number 6 to appear on the upper face is 1/6.  

In statistical mechanics, e.g., that of a gas contained in a finite volume <math> V </math>, both the spatial coordinates <math> q_i </math> and the momentum coordinates <math> p_i </math> of the individual gas elements (atoms or molecules) are finite in the phase space spanned by these coordinates. In analogy to the case of the die, the a priori probability is here (in the case of a continuum) proportional to the phase space volume element <math> \Delta q\Delta p</math> divided by <math>h</math>, and is the number of standing waves (i.e., states) therein, where <math> \Delta q</math> is the range of the variable <math> q </math>  and <math> \Delta p</math> is the range of the variable <math> p</math> (here for simplicity considered in one dimension). In 1 dimension (length <math>L</math>) this number or statistical weight or a priori weighting  is <math> L \Delta p/h </math>. In customary 3 dimensions (volume <math>V</math>) the corresponding number can be calculated to be <math>V 4\pi p^2\Delta p/h^3</math>.<ref>{{cite book |first=H. J. W. |last=Müller-Kirsten |title=Basics of Statistical Physics |edition=2nd |publisher=World Scientific |location=Singapore |year=2013 |at=Chapter 6 }}</ref> In order to understand this quantity as giving a number of states in quantum (i.e., wave) mechanics, recall that in quantum mechanics every particle is associated with a matter wave which is the solution of a [[Schrödinger equation]]. In the case of free particles (of energy <math>\epsilon = {\bf p}^2/2m</math>) like those of a gas in a box of volume <math>V = L^3</math> such a matter wave is explicitly
<math display="block"> \psi \propto \sin(l\pi x/L)\sin(m\pi y/L)\sin(n\pi z/L),</math>
where <math> l, m, n </math> are integers. The number of different  <math>(l,m,n)</math> values and hence states in the region between <math>p, p+dp, p^2 = {\bf p}^2,</math> is then found to be the above expression <math> V4\pi p^2dp/h^3</math> by considering the area covered by these points. 
Moreover, in view of the [[uncertainty relation]], which in 1 spatial dimension is
<math display="block">\Delta q\Delta p \geq h ,</math>
these states are indistinguishable (i.e., these states do not carry labels). An important consequence is a result known as [[Liouville's theorem (Hamiltonian)|Liouville's theorem]], i.e., the time independence of this phase space volume element and thus of the a priori probability. A time dependence of this quantity would imply known information about the dynamics of the system, and hence would not be an a priori probability.<ref>{{cite book |first=A. |last=Ben-Naim |title=Entropy Demystified |publisher=World Scientific |location=Singapore |year=2007 }}</ref> Thus the region
<math display="block"> \Omega := \frac{\Delta q\Delta p}{\int \Delta q\Delta p},\;\;\;  \int \Delta q\Delta p = \mathrm{const.}, </math>
when differentiated with respect to time <math> t</math> yields zero (with the help of Hamilton's equations): The volume at time <math> t</math> is the same as at time zero. One describes this also as conservation of information. 

In the full quantum theory one has an analogous conservation law. In this case, the phase space region is replaced by a subspace of the space of states expressed in terms of a projection operator <math> P</math>, and instead of the probability in phase space, one has the probability density
<math display="block"> \Sigma: = \frac{P}{\text{Tr}(P)},\;\;\; N = \text{Tr}(P) = \mathrm{const.},</math>
where <math>N</math> is the dimensionality of the subspace. The conservation law in this case is expressed by the unitarity of the [[S-matrix]]. In either case, the considerations assume a closed isolated system. This closed isolated system is a system with (1) a fixed energy <math>E</math> and (2) a fixed number of particles <math>N</math> in (c) a state of equilibrium. If one considers a huge number of replicas of this system, one obtains what is called a ''[[microcanonical ensemble]]''. It is for this system that one postulates in quantum statistics the "fundamental postulate of equal a priori probabilities of an isolated system." This says that the isolated system in equilibrium occupies each of its accessible states with the same probability. This fundamental postulate therefore allows us to equate the a priori probability to the degeneracy of a system, i.e., to the number of different states with the same energy.

===Example===
The following example illustrates the a priori probability (or a priori weighting) in (a) classical and (b) quantal contexts.
{{ordered list
| list-style-type = lower-alpha
| 1 = '''Classical a priori probability'''

Consider the rotational  energy E of a diatomic molecule with moment of inertia I in spherical polar coordinates <math> \theta, \phi</math> (this means <math>q</math> above is here <math>\theta, \phi</math>), i.e.
<math display="block"> E = \frac{1}{2I}\left(p^2_{\theta} + \frac{p^2_{\phi}}{\sin^2\theta}\right).</math>
The <math>(p_{\theta}, p_{\phi})</math>-curve for constant E and <math>\theta </math> is an ellipse of area 
<math display="block">\oint dp_{\theta}dp_{\phi} = \pi \sqrt{2IE}\sqrt{2IE}\sin\theta = 2\pi IE\sin\theta .</math>
By integrating over <math>\theta </math> and <math>\phi </math> the total volume of phase space  covered for constant energy E is 
<math display="block">\int^{\phi=2\pi}_{0}\int^{\theta=\pi}_0 2I\pi E\sin\theta d\theta d\phi = 8\pi^2 IE = \oint dp_{\theta}dp_{\phi}d\theta d\phi,</math>
and hence the classical a priori weighting in the energy range <math>dE </math> is
:<math>\Omega \propto</math>  (phase space volume at <math>E+dE</math>) minus  (phase space volume at <math>E</math>) is given by <math> 8{\pi}^2 I dE.</math>
| 2 = '''Quantum a priori probability'''

Assuming that the number of quantum states in a range <math>\Delta q \Delta p </math> for each direction of motion is given, per element, by a factor <math>\Delta q\Delta p/h</math>, the number of states in the energy range dE is, as seen under (a) <math>8\pi^2I dE/h^2 </math> for the rotating diatomic molecule. From wave mechanics it is known that the energy levels of a
rotating diatomic molecule are given by
<math display="block"> E_n = \frac{n(n+1)h^2}{8\pi^2 I}, </math>
each such level being (2n+1)-fold degenerate. By evaluating <math>  dn/dE_n = 1/(dE_n/dn)</math>
one obtains
<math display="block">\frac{dn}{dE_n} = \frac{8\pi^2 I}{(2n+1)h^2}, \;\;\; (2n+1) dn = \frac{8\pi^2 I}{h^2} dE_n.</math>
Thus by comparison with <math>\Omega </math> above, one finds that the approximate number of states in the range dE is given by the degeneracy, i.e.
<math display="block"> \Sigma \propto (2n+1)dn.</math>
Thus the a priori weighting in the classical context (a) corresponds to the a priori weighting here in the quantal context (b).
In the case of the one-dimensional simple harmonic oscillator of natural frequency <math> \nu </math> one finds correspondingly: (a) <math> \Omega \propto dE/\nu </math>, and (b) <math> \Sigma \propto dn </math> (no degeneracy).
Thus in quantum mechanics the a priori probability is effectively a measure of the [[Degeneracy (mathematics)|degeneracy]], i.e. the number of states having the same energy. 

In the case of the hydrogen atom or Coulomb potential (where the evaluation of the phase space volume for constant energy is more complicated) one knows that the quantum mechanical degeneracy is <math> n^2 </math> with <math> E\propto 1/n^2 </math>. Thus in this case <math> \Sigma \propto n^2 dn </math>.
}}
===Priori probability and distribution functions===

In statistical mechanics (see any book) one derives the so-called [[distribution function (physics)|distribution functions]] <math>f</math> for various statistics. In the case of [[Fermi–Dirac statistics]] and [[Bose–Einstein statistics]] these functions are 
respectively
<math display="block"> f^{FD}_i = \frac{1}{e^{(\epsilon_i - \epsilon_0)/kT}+1}, \quad f^{BE}_i = \frac{1}{e^{(\epsilon_i-\epsilon_0)/kT}-1}.</math>
These functions are derived for (1) a system in dynamic equilibrium (i.e., under steady, uniform conditions) with (2) total (and huge) number of particles <math>N = \Sigma_in_i</math> (this condition determines the constant <math>\epsilon_0</math>), and (3) total energy <math>E = \Sigma_in_i\epsilon_i</math>, i.e., with each of the <math>n_i</math> particles having the energy <math>\epsilon_i</math>. An important aspect in the derivation is the taking into account of the indistinguishability of particles and states in quantum statistics, i.e., there particles and states do not have labels. In the case of fermions, like electrons, obeying the [[Pauli principle]] (only one particle per state or none allowed), one has therefore
<math display="block"> 0 \leq f^{FD}_i \leq 1, \quad  \text{whereas} \quad 0 \leq f^{BE}_i \leq \infty.</math>
Thus <math>f^{FD}_i</math> is a measure of the fraction of states actually occupied by electrons at energy <math>\epsilon_i</math> and temperature <math>T</math>. On the other hand, the a priori probability <math>g_i</math> is a measure of the number of wave mechanical states available. Hence
<math display="block"> n_i = f_ig_i.</math>
Since <math>n_i</math> is constant under uniform conditions (as many particles as flow out of a volume element also flow in steadily, so that the situation in the element appears static), i.e., independent of time <math>t</math>, and <math>g_i</math> is also independent of time <math>t</math> as shown earlier, we obtain
<math display="block"> \frac{df_i}{dt} = 0, \quad f_i = f_i(t, {\bf v}_i, {\bf r}_i).</math>
Expressing this equation in terms of its partial derivatives, one obtains the [[Boltzmann transport equation]]. How do coordinates <math>{\bf r}</math> etc. appear here suddenly? Above no mention was made of electric or other fields. Thus with no such fields present we have the Fermi-Dirac distribution as above. But with such fields present we have this additional dependence of <math>f</math>.