Content deleted Content added
Quibbler II (talk | contribs) Signs of $\lambda$ where incorrectly flipped halfway through the calculation. Signs are fixed, results obviously didn't change. |
|||
Line 47:
0 & \equiv \delta \mathcal{L} \\
&= \delta \left( - \sum_i k_\text{B} \rho_i \ln \rho_i \right) + \delta \left( \lambda_1 - \sum_i \lambda_1 \rho_i \right) + \delta \left( \lambda_2 U - \sum_i \lambda_2 \rho_i E_i \right) \\
&= \sum_i \bigg[ \delta \Big( - k_\text{B} \rho_i \ln \rho_i \Big)
&= \sum_i \left[ \frac{\partial}{\partial \rho_i } \Big( - k_\text{B} \rho_i \ln \rho_i \Big) \, \delta ( \rho_i )
&= \sum_i \bigg[ -k_\text{B} \ln \rho_i - k_\text{B}
\end{align}</math>
Since this equation should hold for any variation <math> \delta ( \rho_i ) </math>, it implies that
<math display="block"> 0 \equiv - k_\text{B} \ln \rho_i - k_\text{B}
Isolating for <math> \rho_i </math> yields
<math display="block">\rho_i = \exp \left( \frac{-k_\text{B}
To obtain <math> \lambda_1 </math>, one substitutes the probability into the first constraint:
<math display="block">\begin{align}
1 &= \sum_i \rho_i \\
&= \exp \left( \frac{-k_\text{B}
\end{align}</math>
where '''<math> Z </math> is a constant number defined as the canonical ensemble partition function''':
<math display="block">Z \equiv \sum_i \exp \left( - \frac{\lambda_2}{k_\text{B}} E_i \right) .</math>
Isolating for <math> \lambda_1 </math> yields <math> \lambda_1 =
Rewriting <math> \rho_i </math> in terms of <math> Z </math> gives
<math display="block"> \rho_i = \frac{1}{Z} \exp \left( - \frac{\lambda_2}{k_\text{B}} E_i \right) .</math>
Rewriting <math> S </math> in terms of <math> Z </math> gives
<math display="block">\begin{align}
S &= - k_\text{B} \sum_i \rho_i \ln \rho_i \\
&= - k_\text{B} \sum_i \rho_i \left( - \frac{\lambda_2}{k_\text{B}} E_i - \ln(Z) \right) \\
&=
&=
\end{align}</math>
To obtain <math> \lambda_2 </math>, we differentiate <math> S </math> with respect to the average energy <math> U </math> and apply the [[first law of thermodynamics]], <math> dU = T dS - P dV </math>:
<math display="block">\frac{dS}{dU} =
Thus the canonical partition function <math> Z </math> becomes
|