Una solución para su ejercicio se obtienen mediante el método de los multiplicadores de lagrange.
Las limitaciones que tenemos para satisfacer son tres:
1) $max_{P_j} S[p]$
2) $\sum_{j} P_j E_j = E$
3) $\sum_{j} P_j =1$
Eligiendo como parámetros de variación $\lambda$$\gamma$, los cálculos siga como:
$S_{\lambda,\gamma}=-\sum_{j} P_j ln{P_j} + \lambda(\sum_{j} P_j E_j - E)+\gamma(\sum_{j} P_j -1)$
$\frac{ \partial S_{\lambda,\gamma} }{ \partial{P_j} }=0\Longrightarrow -ln{P_j}-1+\lambda E_j+\gamma=0 \Longrightarrow P_j=\frac{ e^{- \lambda E_j} }{ e^{1-\gamma} }$
Llamar a $Z$ la normalización constante, tenemos:
$P_j= \frac{ e^{-\lambda E_j} }{Z}$
$\sum_{j} P_j=1 \rightarrow e^{1-\gamma}=\sum_{j}e^{-\lambda E_j}:=Z $
Los valores de la $P_j$ maximizar la entropía desde:
$\frac{ \partial^2 S_{\lambda , \gamma}}{ \partial P^2_j}=\frac{-1}{P_j}$
El parámetro $\lambda$ puede determinarse a partir de:
$\frac{ \sum_{j} E_j e^{-\lambda E_j} }{ Z } = E$
Ahora voy a tratar de mostrar cómo el parámetro de $\lambda$ puede ser interpretado como la inversa de la temperatura. Vamos a introducir la función de $F(\lambda , E_j ):=ln{Z}$. Tenemos:
$dF=\frac{\partial F}{\partial \lambda} d\lambda + \sum_{j} \frac{\partial F}{\partial E_j} dE_j=-Ed\lambda-\lambda \sum_{j} \frac{N_j}{N} dE_j $ donde $P_j=\frac{N_j}{N}$
La última fórmula puede escribirse como:
$d(F+E\lambda)=\lambda(dE-\sum_{j} \frac{N_j}{N} dE_j):=\lambda dQ$
La última asociación se sigue del hecho físico de que, si miramos desde el punto de vista de la mecánica cuántica si te gusta) $-\sum_{j} \frac{N_j}{N} dE_j$ puede ser interpretado como el trabajo realizado en el sistema para variar los niveles de energía de $E_j$ $E_j+dE_j$ $dE$es la variación de la energía interna. Por lo $dE-\sum_{j} \frac{N_j}{N} dE_j$ es la cantidad de calor $dQ$ intercambiado por el ensemble con el exterior.
Dado que sólo tenemos una diferencial exacta que involucren calor de la termodinámica se puede concluir que el $\lambda=1/T$
Ahora podemos definir:
$f:=\frac{-1}{\beta} ln{Z}=E-TS$
Por lo que el máximo de los valores de la entropía están incrustados en la definición de energía libre y determinar el valor mínimo de la misma.
El mismo procedimiento se aplica si partimos de la energía libre:
Restricciones: $E=\sum_{j} P_j E_j$ $\sum_{j} P_j = 1$
Parámetro del multiplicador de lagrange: $\mu$
$f_\mu=\sum_{j}P_jE_j+T\sum_{j}P_j ln(P_j)+\mu(\sum_{j} P_j -1)$
$\frac{\partial f}{\partial P_j}=E_j + T ln{P_j} + T + \mu =0 \Longrightarrow P_j=e^{\frac{-E_j}{T}} e^{\frac{\mu}{T}-1} $
$\sum_{i}P_i=1 \rightarrow e^{\frac{-\mu}{T}-1}=\frac{1}{\sum_{i} e^{\frac{-E_i}{T}} }:=\frac{1}{Z} \Longrightarrow P_i= \frac{e^{\frac{-E_i}{T}}}{Z}$