One can derive all the numerous thermodynamic potentials (Helmholtz, Gibbs, Grand, Enthalpy) by Legendre transformations, but I'm interested in seeing each from Stat Mech instead (ie taking the log of a partition function). I can do this easily for all but the Enthalpy, which has me stumped.
Easy example: Deriving Helmholtz
The Helmholtz energy can be derived by putting our system S in thermal contact with a heat bath B at temperature T and maximizing the total entropy:
Stot/kb=log∑ESΩB(Etot−ES)ΩS(ES)
If the heat bath is big enough to maintain constant temperature 1/T=(dSBdEB)VB, then we can say SB≈−ES/T+S0 for small ES. Then ΩB≈Ω0e−βES where S0=kblogΩ0, β=1/kbT, so
Stot/kb=log∑ESΩ0e−βESΩS(ES)=S0/kb+logZS where ZS=∑ESe−βESΩS(ES). So maximizing the total entropy is just maximizing ZS. If we define the Helmholtz Free Energy AS as −βAS=logZS, and use the Shannon entropy S/kb=−∑plogp, we see
SS/kb=−∑ΩS(ES)e−βESZlogΩS(ES)e−βESZ SS/kb=βZ∑ΩS(ES)ESe−βES+logZ SS/kb=−βZ∂βZ+logZ SS/kb=β⟨ES⟩−βAS AS=⟨ES⟩−TSS
The other thermodynamic potentials at fixed T are similarly easy to derive.
But now try deriving Enthalpy
The same procedure does not work for enthalpy, because,
Stot/kb=log∑VSΩB(Vtot−VS,EB0+pVS)ΩS(VS,ES0−pVS)
...if the bath is big enough to maintain a constant temperature, then its total entropy is constant as a function of VS. That is, if VS increases, the bath entropy decreases due to less volume, but the bath energy increases by the same amount due to increased energy. So, to first-order, ΩB is constant and the total entropy splits into a sum of both individual subsystem entropies.
Is there a way to derive the enthalpy from stat mech considerations, as there is for the other potentials, rather than by Legendre transforming the energy?
By "derive the enthalpy", I mean "derive that the quantity which should be minimized in equilibrium is given by H=⟨E⟩+p⟨V⟩."
No comments:
Post a Comment