Si todavía es real para usted. La densidad de probabilidad
$$
e^{-\beta f\left( \mu\right) }
$$
es la llamada canónica (Gibbs) de distribución.Hay un montón de métodos de cómo
para derivar de ella. Puedo reproducir la más simple.
Imaginemos que el sistema tiene la Hamiltoniana $H\left( \mu\right) $ y te gustaría estudiar para una temperatura determinada. Con el fin de ajustar la temperatura de poner su sistema en el interior de un termostato, de manera que su sistema de intercambios de energía sólo con la
termostato, pero el volumen y el número de partículas es constante. Vamos
supongamos que el termostato es un gran tanque lleno de un gas ideal, de modo que su
energía:
$$
h=\sum_{i=1}^{N}\frac{P_{i}^{2}}{2m}.
$$
El total del sistema (el sistema+termostato) es aislado por lo tanto la energía total es
fijo. Por lo tanto, la distribución con respecto al total de la energía es un
delta-función:
$$
\rho\left( E\right) =\Lambda\delta\left( h+Que\right) ,
$$
donde $\Lambda$ es una cierta normalización factor que
$$
\int \rho\left( E\ \ derecho)\,d\Gamma =1,\qquad\left( 1\right)
$$
donde $d\Gamma$ es un elemento de la totalidad del espacio de fase:
$$
d\Gamma=d\mu\prod_{i=1}^{N}d^{3}P_{i}\,d^{3}Q_{i}.
$$
Vamos a integrar a cabo todos los grados de libertad del termostato:
$$
\rho\left( H\right) =\int\prod_{i=1}^{N}d^{3}P_{i}\,d^{3}Q_{i}
\,\Lambda\delta\left( H+\sum_{i=1}^{N}\frac{P_{i}^{2}}{2m}-E\right) =\Lambda
V^{N}\int\prod_{i=1}^{N}d^{3}P_{i}\,\,\delta\left( H+\sum_{i=1}^{N}
\frac{P_{i}^{2}}{2m}-E\ \ derecho) ,
$$
donde $V$ es el volumen de termostato. La integración puede ser una medida de
simplificado de la siguiente manera:
$$
\int\left[\prod_{i=1}^{N}d^{3}P_{i}\right] f(\epsilon)=\frac{2\pi^{3N/2}}{\Gamma\left( 3N/2\right)
}\int\left[\epsilon^{3N-1}d\epsilon\right] f(\epsilon),
$$
donde
$$
\epsilon^{2}=\sum_{i=1}^{N}P_{i}^{2}.
$$
Así, la integración puede procederse de la siguiente manera:
$$
\rho\left( H\right) =\Lambda V^{N}\frac{2\pi^{3N/2}}{\Gamma\left(
3N/2\right) }\int\,d\epsilon\,\epsilon^{3N-1}\,\delta\left( H+\frac
{\epsilon^{2}}{2m}-E\right) =\Lambda V^{N}\frac{2m\pi^{3N/2}}{\Gamma\left(
3N/2\right) }\,\a la izquierda( E-H\right) ^{\frac{3N}{2}-1}.
$$
Ahora vamos a considerar el $N\rightarrow\infty$ límite, así que
$$
\frac{E}{N}\approx\frac{h}{N}=\frac{3T}{2}.
$$
La distribución toma la forma:
$$
\rho\left( H\right) \sim\left( 1-\frac{H}{E}\right) ^{3N-2}\approx\left(
1-\frac{H}{\frac{3N}{2}T}\right) ^{\frac{3N}{2}-1}\approx\exp\left(
-\frac{H}{T}\right) .
$$
La normalización factor puede ser encontrado a partir de la normalización de la condición (1).
Finalmente, la densidad de probabilidad para la energía de tu sistema tiene la forma:
$$
\rho\left( H\right) =\frac{e^{-\beta f\left( \mu\right) }}{Z},\quad Z=\int
d\mu\,e^{-\beta f\left( \mu\right) }.
$$
De hecho, el resultado es independiente de la naturaleza del termostato, ver, por ejemplo, L. D. Landau, E. M. Lifshitz, Volumen 5 del Curso de Física Teórica, física Estadística Parte 1 Ch. III, La distribución de Gibbs.