Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Summary of Ensembles

The Big Picture: Why Multiple Ensembles?

NVEallow E to fluctuateNVTallow V to fluctuateNPT\text{NVE} \xrightarrow{\text{allow } E \text{ to fluctuate}} \text{NVT} \xrightarrow{\text{allow } V \text{ to fluctuate}} \text{NPT}
NVTallow N to fluctuateμVT\text{NVT} \xrightarrow{\text{allow } N \text{ to fluctuate}} \mu\text{VT}

Entropy and the Origin of Ensemble Distributions

Entropy is given by the Shannon-Gibbs entropy formula:

S([p])=kBipilogpiS([p]) = -k_B \sum_{i} p_i \log p_i

where pip_i is the probability of the iith microstate.

Physical Constraints for Equilibrium

For a system to maintain equilibrium values, the microstate probabilities must satisfy the following conditions:

  1. Normalization:

    ipi=1\sum_{i} p_i = 1
  2. Constraint to Maintain the Expectation Value of an Observable XX (e.g., Energy or Volume):

    ipiXi=X,ipiYi=Y\sum_{i} p_i X_i = \langle X \rangle, \quad \sum_{i} p_i Y_i = \langle Y \rangle
ErE,VrV,NrNE_r\gg E,\quad V_r \gg V, \quad N_r \gg N
logΩr(EtE,VrV,NrN)constβEβPV+βμN\log \Omega_r(E_t-E,\, V_r-V,\, N_r-N)\approx \text{const} - \beta E -\beta PV + \beta \mu N
P(E,N,V)ΩΩreS(E,N,V)/kBeβEeβμNeβPVP(E, N, V) \sim \Omega \cdot \Omega_r \sim e^{S(E, N, V)/k_B} \cdot e^{-\beta E} \cdot e^{\beta \mu N} \cdot e^{-\beta PV}

Comparison of Ensembles


EnsembleFixedFluctuatingPartition FunctionThermodynamic Potential
Microcanonical (NVE)N,V,EN, V, EΩ(E)\Omega(E)S=kBlnΩS = k_B \ln \Omega
Canonical (NVT)N,V,TN, V, TEEZ=eβEZ = \sum e^{-\beta E}F=kBTlnZF = -k_BT \ln Z
Isothermal-Isobaric (NPT)N,p,TN, p, TE,VE, VΔ=ZeβpVdV\Delta = \int Z\, e^{-\beta pV} dVG=kBTlnΔG = -k_BT \ln \Delta
Grand Canonical (μ\muVT)μ,V,T\mu, V, TE,NE, NΞ=ZNeβμN\Xi = \sum Z_N\, e^{\beta \mu N}Ψ=kBTlnΞ=PV\Psi = -k_BT \ln \Xi = -PV

EnsembleP(microstate) P(\text{microstate}) P(macrostate) P(\text{macrostate})
Microcanonical (NVE)Pi=1Ω(E) P_i = \frac{1}{\Omega(E)} P(E)eS(E)/kB P(E) \sim e^{S(E)/k_B} (entropy-dominated)
Canonical (NVT)PieβEi P_i \sim e^{-\beta E_i} P(E)eS(E)/kBβE P(E) \sim e^{S(E)/k_B - \beta E} (entropy-weighted by energy)
Isothermal-Isobaric (NPT)Pieβ(Ei+PV) P_i \sim e^{-\beta (E_i + PV)} P(E,V)eS(E,V)/kBβ(E+PV) P(E, V) \sim e^{S(E, V)/k_B - \beta (E + PV)}
Grand Canonical (μ\muVT)Pieβ(μNEi) P_i \sim e^{\beta (\mu N - E_i)} P(N,E)eS(N,E)/kB+β(μNE) P(N, E) \sim e^{S(N,E)/k_B + \beta (\mu N - E)}

Extensive vs Intensive Variables

dU=TdSpdV+μdN+BdM+=ifidXidU = TdS-pdV + \mu dN + BdM + \dots = \sum_i f_i dX_i

Laplace Transform and Ensemble Connections

Legendre Transform and Thermodynamic Potentials

PotentialLegendre transformNatural variablesTotal differential
Internal energyUUS,V,NS, V, NdU=TdSpdV+μdNdU = TdS - pdV + \mu dN
HelmholtzFFUTSU - TST,V,NT, V, NdF=SdTpdV+μdNdF = -SdT - pdV + \mu dN
EnthalpyH HU+PVU + PVS,P,NS, P, NdH=TdS+VdP+μdNdH = TdS + VdP + \mu dN
GibbsG GUTS+PVU - TS + PVT,P,NT, P, NdG=SdT+VdP+μdNdG = -SdT + VdP + \mu dN
Grand potentialΨ \PsiUTSμNU - TS - \mu NT,V,μT, V, \mudΨ=SdTpdVNdμd\Psi = -SdT - pdV - Nd\mu

Partition Functions and Legendre Transforms

Ψ(f1,,fn,Xn+1,,XN)=U(X1,,XN)(f1X1++fnXn)\Psi(f_1, \dots, f_{n}, X_{n+1}, \dots, X_{N}) = U(X_1, \ldots, X_N) - (f_1 X_1+\cdots+f_nX_n)
Z(f1,,fn,Xn+1,,XN)=eβΨ(f1,,fn,Xn+1,,XN)Z(f_1, \dots, f_n, X_{n+1}, \dots, X_N) = e^{-\beta \Psi(f_1, \dots, f_{n}, X_{n+1}, \dots, X_{N})}

Fluctuation-Response Relations

Energy fluctuations (Canonical Ensemble):

$$

$$

\sigma^2_E = k_B T^2 C_V

$$

$$

Volume fluctuations (Isothermal-Isobaric Ensemble):

$$

$$

\sigma^2_V = k_B T, V, \kappa_T

$$

$$

Particle number fluctuations (Grand Canonical Ensemble):

$$

$$

\sigma^2_N = \frac{\langle N \rangle^2 k_B T \kappa_T}{V}

$$

$$

Unified pattern

EnsembleFluctuating variableResponse functionFluctuation formula
NVTEnergyEEHeat capacityCVC_VσE2=kBT2CV\sigma_E^2 = k_BT^2 C_V
NPTVolumeVVCompressibilityκT\kappa_TσV2=kBTVκT\sigma_V^2 = k_BT\, V\, \kappa_T
μ\muVTParticle numberNNCompressibilityκT\kappa_TσN2=N2kBTκTV\sigma_N^2 = \frac{N^2 k_BT \kappa_T}{V}

Key Insights