One fundamental function can generate a whole property table
A pressure equation alone can fit p–V–T data while leaving energy and heat-capacity behaviour disconnected. A fundamental equation written as Helmholtz energy a(T,ρ) contains more structure. Its derivatives generate pressure, internal energy, entropy and response functions that are linked by thermodynamics.
Modern reference equations often use the dimensionless form α=a/(RT), split into ideal-gas and residual contributions and expressed in reduced temperature and density. The power of the formulation comes with a responsibility: derivatives, reference states, validity ranges and stability checks must all remain consistent.
This guide is a specialist extension of the Gold Standard of Thermodynamics. The sequence is deliberate: explain the model, calculate one case, change the context, and check what remains valid.
Find your route through this guide
| If you want to… | Start | Continue |
|---|---|---|
| understand the mechanism | Chapter 1 | Chapters 2–3 |
| work through the calculation | Chapter 4 | Chapter 5 |
| interpret evidence and diagnostics | Chapter 6 | Chapter 7 |
| transfer the model to a changed case | Chapter 8 | Chapter 9 |
| check limits, FAQs and practice | Chapter 10 | Chapters 11–12 |
Chapter contents
Foundations and setup
Calculation and evidence
1. Why start from Helmholtz energy
For a simple compressible pure fluid, molar Helmholtz energy a=u−Ts naturally uses temperature T and molar density ρ as independent variables. Its differential da=−s dT+(p/ρ²)dρ shows that entropy and pressure are first derivatives of one function.
Because the properties share a parent function, Maxwell relations are built into a smooth formulation. Fitting unrelated correlations for each property can violate those cross-derivative links.
| Derivative of a(T,ρ) | Property |
|---|---|
| −(∂a/∂T)ρ | Molar entropy s |
| ρ²(∂a/∂ρ)T | Pressure p |
| a+Ts | Molar internal energy u |
2. Reduced variables and the ideal–residual split
Reference equations commonly define α=a/(RT), inverse reduced temperature τ=Tᵣ/T and reduced density δ=ρ/ρᵣ. Then α(τ,δ)=α⁰(τ,δ)+αʳ(τ,δ).
The ideal part carries ideal-gas behaviour and reference-state terms. The residual part represents intermolecular departure from ideality and vanishes as density approaches zero for a well-formed model.
Reducing parameters may be critical constants or mixture reducing functions. They are model definitions, not interchangeable scaling choices.
3. Pressure and caloric properties from derivatives
For a pure fluid in the common reduced formulation, the compressibility factor is Z=p/(ρRT)=1+δ αʳ_δ. Residual internal energy satisfies uʳ/(RT)=τ αʳ_τ, and residual enthalpy satisfies hʳ/(RT)=τ αʳ_τ+δ αʳ_δ.
Heat capacities require second derivatives. This is why an equation that looks smooth in α can still misbehave if its curvature is poorly fitted or evaluated outside the valid range.
Derivative notation must state what is held constant. Confusing τ and T derivatives reverses signs and scaling factors.
4. Worked example: a toy residual term
Take a deliberately simple residual model αʳ=−Aδτ with A=0.10 at δ=0.50 and τ=1.20. Then αʳ=−0.060, αʳ_δ=−Aτ=−0.12 and αʳ_τ=−Aδ=−0.05.
| Quantity | Formula | Value |
|---|---|---|
| Z | 1+δαʳ_δ | 0.94 |
| uʳ/(RT) | ταʳ_τ | −0.06 |
| hʳ/(RT) | ταʳ_τ+δαʳ_δ | −0.12 |
One term changes several properties coherently. That coherence is the point; the numerical model is too simple for engineering use.
5. Reference states and the ideal contribution
Absolute internal energy and entropy require conventions. Helmholtz equations encode reference-state constants in the ideal part so property differences and tabulated values follow the chosen convention.
Pressure and many derivative response properties are unaffected by additive constants, while enthalpy and entropy values can shift. Comparing two software packages without aligning reference states can create a false discrepancy.
Always record molar versus mass basis, gas constant, reference state and composition basis.
6. Heat capacities, compressibility and sound speed
Second derivatives of α generate cᵥ, cₚ, isothermal compressibility, expansivity and speed of sound. These properties are sensitive tests because they probe curvature rather than only the fitted surface value.
A model that matches pressure data may still miss heat capacity or sound speed if the derivatives are wrong. Reference-quality fitting therefore uses multiple property types.
Near the critical region, response functions can change sharply and numerical evaluation becomes especially demanding.
7. Thermodynamic consistency checks
Mixed derivatives of a must commute where the function is smooth, producing Maxwell relations. Calculated cₚ−cᵥ must agree with the expansion/compressibility identity. Pressure recovered from α must match pressure used in any inverse state calculation.
| Check | Warning sign |
|---|---|
| Low-density limit | Residual terms do not vanish |
| Mixed derivatives | Numerical order changes the result |
| Heat capacities | Negative values in a stable ordinary state |
| Mechanical stability | (∂p/∂ρ)T≤0 outside intended metastable handling |
| Sound speed | Negative squared value |
8. Transfer check: halve the density
Keep A=0.10 and τ=1.20 but reduce δ from 0.50 to 0.25. For the toy term, αʳ becomes −0.030, Z becomes 0.97, uʳ/(RT)=−0.03 and hʳ/(RT)=−0.06.
The departures move halfway toward zero because this particular residual term is linear in density. A real multiparameter equation contains several powers and damping terms, so not every residual property scales linearly.
What transfers is the derivative procedure and the ideal-gas limit. The simple scaling does not transfer to an arbitrary equation.
9. Phase boundaries and density roots
At a given subcritical temperature, a pressure equation can admit several density roots. Stable vapour and liquid states are selected through phase-equilibrium conditions such as equal chemical potential and pressure, not by taking any numerical root.
Inside unstable regions, mechanical stability fails. Metastable extensions may be useful for research, but they should not be confused with equilibrium property values.
Two-phase states require quality or composition information in addition to T and p because the phase densities differ.
10. Validity ranges and misconceptions
- A Helmholtz formulation does not make extrapolation safe; every fitted equation has a validity range and uncertainty.
- Residual does not mean measurement error; it means departure from the ideal contribution.
- Additive reference constants affect some tabulated values even when pressure is unchanged.
- Transport properties such as viscosity need separate models; they are not ordinary derivatives of equilibrium Helmholtz energy.
- Mixtures require composition derivatives and mixing/reducing rules beyond the pure-fluid equations.
- A numerical density root is not automatically the stable phase.
Use the fluid-specific equation and its cited uncertainty, not a generic ‘REFPROP value’ detached from the underlying correlation.
11. Frequently asked questions
Why use density rather than pressure as an independent variable?
Helmholtz energy naturally uses T and density, and density distinguishes coexisting phase states that can share T and p.
What does the residual term approach at zero density?
It should approach zero so the model recovers ideal-gas behaviour.
Does one equation provide viscosity?
No. Equilibrium thermodynamic properties come from the fundamental equation; viscosity and conductivity use separate transport correlations.
12. A reproducible property checklist
| Stage | Record |
|---|---|
| Equation | Fluid, version, source and validity range |
| Variables | T, density, composition and reducing parameters |
| Basis | Molar or mass-specific units and R |
| Derivatives | Analytic or numerical method and tolerances |
| Phase | Root selection and equilibrium test |
| Verification | Limits, Maxwell identities and independent checkpoints |
The gold standard is not a long table of digits. It is a traceable fundamental equation whose derivatives, phase choice, reference convention and uncertainty all agree.
Reading routes · Previous chapter · Contents · Continue reading
Continue through the Gold Standard ecosystem
Return to the Gold Standard of Thermodynamics, the Gold Standard of Education parent, or the How X Works hub.
Read next in this subject
- Same Endpoints, Different Pressure–Volume Paths
- Heating Through a Table of Changing Heat Capacities
- Finite Bodies, Entropy Generation and the Direction of Heat Flow
- Joule–Thomson Inversion, Throttling and Temperature Change