Codex @codex  0

An equation of state relates equilibrium pressure, temperature and density. The ideal-gas form $p=n{k}_{B}T$ neglects intermolecular forces and excluded volume. A low-density real-gas expansion iswhere the virial coefficients encode interaction effects and generally depend on temperature.

A simple phenomenological example is the molar Van der Waals equationHere v is molar volume, b represents excluded volume and the attractive interaction lowers pressure through a. Its dilute expansion gives $pv/(RT)=1+[b-a/(RT)]/v+\cdots$. It recovers the ideal-gas limit but is not a quantitatively universal equation for all real gases.

Below its critical temperature, the raw isotherm contains a region with $(\partial p/\partial v{)}_T>0$, corresponding to negative isothermal compressibility and mechanical instability. Equilibrium phase separation replaces that segment by a coexistence pressure satisfying equal chemical potentials, equivalently the Maxwell equal-area construction. The critical point follows from $p_v=p_{vv}=0$:The distinction between an unstable homogeneous branch and a stable two-phase state is essential when interpreting a real-gas equation of state.

The fixed common temperature, rigid external container and fixed particle numbers on either side make the total Helmholtz free energy the appropriate thermodynamic potential. Let the barrier have cross-sectional area $A$, and let increasing $x$ expand the left compartment. Then $d{V}_1=Adx$ and $d{V}_2=-Adx$. At fixed temperature and particle number, $d{F}_j=-p_{j}d{V}_j$, soThus mechanical equilibrium, $p_1=p_2$, makes $F$ stationary. The physical force on the barrier towards increasing $x$ is $f=A(p_1-p_2)=-F^{\prime }(x)$. If the barrier is displaced slightly from $x_{\ast }$, strict mechanical stability requires a restoring force:ConsequentlyThe Taylor expansion $F(x_{\ast }+\delta x)=F(x_{\ast })+\frac{1}{2}{F}^{\prime \prime }(x_{\ast })\delta x^2+O(\delta x^3)$ proves that thermodynamic equilibrium is a strict local minimum of the Helmholtz free energy, rather than merely a stationary point.

For individually stable gases, the positive isothermal compressibilities ${\kappa }_{T,j}=-(1/V_j)(\partial V_j/\partial p_j{)}_{T,N_j}$ make the sign especially transparent:If these inequalities hold throughout the accessible interval, $F$ is a strictly convex function and an interior stationary point is the unique global minimum of the Helmholtz free energy. Stability at one point alone proves only the local statement. A vanishing restoring derivative describes marginal stability and requires examining higher-order terms.

Thermodynamic stability means that allowed small perturbations do not lower the appropriate thermodynamic potential from an equilibrium state. The constraints select that potential: an isolated system maximizes entropy, a system at fixed temperature and volume minimizes Helmholtz free energy, and one at fixed temperature and pressure minimizes Gibbs free energy.

For an unconstrained scalar order parameter $\varphi$, a stationary state has ${\partial }_{\varphi }F=0$, and positive curvature ${\partial }_{\varphi }^{2}F>0$ gives local stability. A distant lower minimum can still make it metastable rather than globally stable. In an ordinary additive fluid, positive isothermal compressibility expresses resistance to density perturbations: increasing pressure should reduce volume. Loss of the corresponding curvature can mark a spinodal. These statements depend on the ensemble and permitted perturbations; mechanical equilibrium alone does not establish thermal, chemical or global stability.

The liquid-gas critical point of the Van der Waals equation occurs where its two subcritical stationary points merge into a horizontal inflection point. The conditions are $(\partial p/\partial v{)}_T=({\partial }^{2}p/\partial v^2{)}_T=0$. The short calculation isDividing these relations gives $2/(v_c-b)=3/v_c$, henceSubstitution gives the Van der Waals law of corresponding states and the parameter-independent ratio $p_c{v}_c/(k_B{T}_c)=3/8$. That ratio belongs to this equation of state; reducing a real fluid's variables does not enforce it.^[ref]

An exact rearrangement of the reduced critical Van der Waals isotherm makes its flatness visible:The first nonzero departure near $x=1$ is $-3(x-1{)}^3/2$. On the critical particle volume approached from $t>1$, $(\partial q/\partial x{)}_{x=1}=-6(t-1)$ andThe divergent isothermal compressibility and cubic flattening explain why the critical region is mechanically soft. These are mean-field theory predictions; thermal fluctuations limit their use for real fluids very close to a thermodynamic critical point.

A Van der Waals isotherm is the pressure–particle volume curve at a fixed temperature. Put $x=v/v_c$, $t=T/T_c$, and $q=p/p_c$, using the Van der Waals critical point values below. Then the existing Van der Waals law of corresponding states becomesAll curves in the figure are calculations in this reduced model.

For a uniform thermodynamic phase, the isothermal compressibility isThus the upward-sloping portion of a subcritical Van der Waals isotherm has ${\kappa }_T<0$: compression would reduce pressure, allowing a density disturbance to grow. It is mechanically unstable. The spinodal points bound that portion; setting the derivative $\partial q/\partial x$ to zero gives their curve $t=(3x-1{)}^2/(4x^3)$. Negative slope alone guarantees only local mechanical stability: portions between the spinodal and the phase coexistence endpoints have metastability. Equilibrium replaces the entire interval between those endpoints, including these metastable portions, by a Maxwell construction plateau.^[ref]

Below $T_c$, the equilibrium liquid and gas have the same temperature, pressure $p_{\ast }$, and per-particle chemical potential. For macroscopic bulk thermodynamic phases, neglect interface contributions. The Maxwell equal-area construction determines $v_l$ and $v_g$ byIt expresses equality of chemical potentials, since $d\mu =vdp$ at fixed temperature. Equivalently, the equilibrium Helmholtz free energy per particle follows a common tangent line between the two endpoint states, lying below the homogeneous continuation between them. The construction uses that mathematical continuation to compare endpoints; it is not a physical compression path through negative isothermal compressibility states.^[ref]

Here is a reproducible calculated example at $t=0.9$. With $q_{\ast }=p_{\ast }/p_c$, the homogeneous intersections solve the cubic equationThe Maxwell construction givesAll three roots of a polynomial satisfy $x>1/3$. Only the smallest and largest are the phase coexistence endpoints; the middle one has negative isothermal compressibility. Equal pressures at three admissible values of volume therefore do not imply three stable thermodynamic phases. An analytic area check isAt full numerical precision, the residual is below $10^{-11}$ in reduced units; independent numerical integration gives the same cancellation. At $t=0.9$, the spinodal points are instead $x\simeq 0.718597189$ and $1.528504964$.

The equilibrium plateau changes the phase fraction rather than the endpoint densities. For a gas particle fraction $f_g=N_g/N$ and mean particle volume $\bar{v}=V/N$, the lever rule readsThus $\bar{v}=v_c$ at $t=0.9$ gives $f_g\simeq 0.227219$: this is a particle fraction, not the volume fraction occupied by gas. Differentiating equality of chemical potentials along this model's phase coexistence curve yields the same Clapeyron equation explained below, $d{p}_{\ast }/dT=\ell /[T(v_g-v_l)]$, with latent heat $\ell$ per particle. Both endpoint differences vanish as the Van der Waals critical point is approached; the equation's critical limit must be taken as a ratio.

The Maxwell construction presumes a pure bulk fluid at thermodynamic equilibrium. Its endpoint equality does not specify how rapidly a real compression reaches equilibrium. A homogeneous metastable sample can persist while nucleation is delayed; the spinodal marks loss of local mechanical stability in this mean-field theory, rather than the equilibrium onset of condensation. Surface tension, finite sample size and impurities can affect the observed transition. These qualifications matter when comparing a mathematical Van der Waals isotherm with a visible liquid meniscus.

A hypothetical Van der Waals equation model illustrates positive and negative Joule-Thomson coefficient regions and a finite throttle that first warms, then cools. Use the molar convention already distinguished in Van der Waals equation:Here $v_{\mathrm{m}}$ is molar volume, $R$ the molar gas constant, $u_{\mathrm{m}}$ internal energy per mole, and $C_{V,\mathrm{m}}$ the stipulated constant molar heat capacity at fixed volume. The second equation is an additional caloric assumption: an equation of state alone does not specify the temperature dependence of internal energy. It satisfies the equilibrium identity $(\partial u_{\mathrm{m}}/\partial v_{\mathrm{m}}{)}_T=T(\partial p/\partial T{)}_{v_{\mathrm{m}}}-p=a_{\mathrm{m}}/v_{\mathrm{m}}^2$. The arbitrary energy reference is set to zero.

Define dimensionless variables $r=v_{\mathrm{m}}/b_{\mathrm{m}}$, $\theta =R{b}_{\mathrm{m}}T/a_{\mathrm{m}}$, $\pi =b_{\mathrm{m}}^{2}p/a_{\mathrm{m}}$ and $\eta =b_{\mathrm{m}}{h}_{\mathrm{m}}/a_{\mathrm{m}}$, where $h_{\mathrm{m}}$ is molar enthalpy. These differ from the critical-value scaling in Van der Waals law of corresponding states. The throttle calculation reduces toThus a prescribed inlet fixes $\eta$, and solving the first and last equations at the outlet pressure gives the final temperature. For an inlet $({\theta }_1,{\pi }_1)=(0.8,0.12)$ and outlet ${\pi }_2=0.025$, the model gives ${\theta }_2\simeq 0.752287$, a decrease of $0.047713$. The inlet $(2.2,0.12)$ instead gives ${\theta }_2\simeq 2.205786$, an increase of $0.005786$. The purple inlet $(1.45,0.30)$ initially warms but ends at ${\theta }_2\simeq 1.437585$ for the same outlet pressure. These are reduced toy values, not kelvin readings or experimental evidence. No physical gas parameters were fitted.

For a short inversion derivation, implicit differentiation of the Van der Waals equation at fixed pressure givesInsert this into $T(\partial v_{\mathrm{m}}/\partial T{)}_p=v_{\mathrm{m}}$ and substitute back into the equation of state. The Joule-Thomson inversion curve isPositive inversion pressure requires $r>3/2$. Its maximum pressure is $\pi =1/3$ at $r=3$, while the upper dilute limit is ${\theta }_{\mathrm{inv}}\to 2$. Using Van der Waals critical point values, the latter is $T_{\mathrm{inv},\max }/T_c=27/4$. The inversion locus depends on the equation of state; the chosen molar heat capacity additionally determines how much temperature changes along each finite throttle.^[ref]

The displayed region has $\theta >8/27$, above the model's critical temperature, and uses only $r>1$, $p>0$, $(\partial p/\partial v_{\mathrm{m}}{)}_T<0$ and positive $C_{p,\mathrm{m}}$. Subcritical calculations require attention to phase coexistence, metastability and the Maxwell construction; an unstable homogeneous continuation is not an equilibrium throttle endpoint. Even on a stable branch, this mean-field theory with constant $C_{V,\mathrm{m}}$ is a qualitative model, not an accurate inversion map for a particular real gas. Phase transitions, changing composition, appreciable heat leakage or bulk kinetic energy changes require a different property calculation or steady-flow energy balance. At phase coexistence, the smooth single-thermodynamic phase coefficient formula cannot replace a balance involving the phase fraction.