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?

NVE→allow E to fluctuateNVT→allow V to fluctuateNPT\text{NVE} \xrightarrow{\text{allow } E \text{ to fluctuate}} \text{NVT} \xrightarrow{\text{allow } V \text{ to fluctuate}} \text{NPT}
NVT→allow 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])=−kB∑ipilog⁡piS([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
Er≫E,Vr≫V,Nr≫NE_r\gg E,\quad V_r \gg V, \quad N_r \gg N
log⁡Ωr(Et−E, Vr−V, Nr−N)≈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)∼Ω⋅Ωr∼eS(E,N,V)/kB⋅e−βE⋅eβμN⋅e−β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=−kBTln⁡ZF = -k_BT \ln Z
Isothermal-Isobaric (NPT)N,p,TN, p, TE,VE, VΔ=∫Z e−βpVdV\Delta = \int Z\, e^{-\beta pV} dVG=−kBTln⁡ΔG = -k_BT \ln \Delta
Grand Canonical (μ\muVT)μ,V,T\mu, V, TE,NE, NΞ=∑ZN eβμ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)Pi∼e−β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)Pi∼e−β(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)Pi∼eβ(μN−Ei) P_i \sim e^{\beta (\mu N - E_i)} P(N,E)∼eS(N,E)/kB+β(μN−E) P(N, E) \sim e^{S(N,E)/k_B + \beta (\mu N - E)}

Extensive vs Intensive Variables

dU=TdS−pdV+μ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 energyUU—S,V,NS, V, NdU=TdS−pdV+μdNdU = TdS - pdV + \mu dN
HelmholtzFFU−TSU - TST,V,NT, V, NdF=−SdT−pdV+μdNdF = -SdT - pdV + \mu dN
EnthalpyH HU+PVU + PVS,P,NS, P, NdH=TdS+VdP+μdNdH = TdS + VdP + \mu dN
GibbsG GU−TS+PVU - TS + PVT,P,NT, P, NdG=−SdT+VdP+μdNdG = -SdT + VdP + \mu dN
Grand potentialΨ \PsiU−TS−μNU - TS - \mu NT,V,μT, V, \mudΨ=−SdT−pdV−Ndμ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=kBT V κ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