Isothermal–isobaric (NPT) ensemble

Equilibrium ensemble for fixed particle number, temperature, and pressure, with fluctuating volume.
Isothermal–isobaric (NPT) ensemble

The isothermal–isobaric ensemble (often called the NPT ensemble) describes a in contact with reservoirs that fix the temperature (see ) and the pressure (see ), while the particle number NN is held fixed. In this setting the volume VV is not fixed; it is a fluctuating quantity determined by mechanical equilibrium with the pressure reservoir.

Microscopically, one convenient classical description takes microstates (see ) as a pair consisting of a point xx in and a volume parameter V>0V>0. The energy of a microstate is given by a H(x;V)H(x;V), and integration over microstates uses the at fixed VV together with an outer integration over VV.

Probability weight

Let β\beta denote the inverse temperature (see ), with β=1/(kBT)\beta = 1/(k_B T) where kBk_B is the . In the NPT ensemble, the unnormalized weight of a microstate (x,V)(x,V) is proportional to

exp ⁣[β(H(x;V)+pV)], \exp\!\big[-\beta\big(H(x;V) + pV\big)\big],

so the probability density has the schematic form

P(dxdV)eβ(H(x;V)+pV)dxdV, \mathbb{P}(dx\,dV) \propto e^{-\beta(H(x;V)+pV)}\,dx\,dV,

where dxdx denotes the phase-space measure at fixed VV.

The combination H+pVH+pV is the enthalpy-like quantity that plays the role of “effective energy” when pressure is imposed rather than volume.

Normalization and thermodynamic potential

The normalization constant of the NPT distribution is the Δ(β,p,N)\Delta(\beta,p,N). It can be expressed as a Laplace transform of the Z(β,V,N)Z(\beta,V,N):

Δ(β,p,N)=0dVeβpVZ(β,V,N). \Delta(\beta,p,N) = \int_{0}^{\infty} dV\, e^{-\beta pV}\, Z(\beta,V,N).

In equilibrium statistical mechanics, Δ\Delta generates the through

G(T,p,N)=kBTlnΔ(β,p,N), G(T,p,N) = -k_B T \,\ln \Delta(\beta,p,N),

up to conventions that become irrelevant in the . This matches the thermodynamic Legendre structure: NPT is the natural ensemble for GG (compare with the which is natural for the ).

Averages and fluctuations

Ensemble averages are defined as in any : for an observable A(x,V)A(x,V),

ANPT=1Δ(β,p,N)0dVdx  A(x,V)eβ(H(x;V)+pV). \langle A\rangle_{NPT} = \frac{1}{\Delta(\beta,p,N)}\int_{0}^{\infty} dV \int dx\; A(x,V)\, e^{-\beta(H(x;V)+pV)}.

Derivatives of lnΔ\ln\Delta generate moments of VV:

V=lnΔ(βp),Var(V)=2lnΔ(βp)2. \langle V\rangle = -\frac{\partial \ln \Delta}{\partial(\beta p)}, \qquad \mathrm{Var}(V) = \frac{\partial^{2} \ln \Delta}{\partial(\beta p)^{2}}.

In particular, the isothermal compressibility can be read from volume fluctuations:

κT=1V(Vp)T=βVar(V)V. \kappa_T = -\frac{1}{\langle V\rangle}\left(\frac{\partial \langle V\rangle}{\partial p}\right)_T = \frac{\beta\,\mathrm{Var}(V)}{\langle V\rangle}.

These relations are instances of general fluctuation–response identities (see ).