Equation of state
In physics and thermodynamics, an equation of state is a relation between state variables.^{[1]} More specifically, an equation of state is a thermodynamic equation describing the state of matter under a given set of physical conditions. It is a constitutive equation which provides a mathematical relationship between two or more state functions associated with the matter, such as its temperature, pressure, volume, or internal energy. Equations of state are useful in describing the properties of fluids, mixtures of fluids, solids, and even the interior of stars.
Contents
Overview
The most prominent use of an equation of state is to correlate densities of gases and liquids to temperatures and pressures. One of the simplest equations of state for this purpose is the ideal gas law, which is roughly accurate for weakly polar gases at low pressures and moderate temperatures. However, this equation becomes increasingly inaccurate at higher pressures and lower temperatures, and fails to predict condensation from a gas to a liquid. Therefore, a number of more accurate equations of state have been developed for gases and liquids. At present, there is no single equation of state that accurately predicts the properties of all substances under all conditions.
In addition, there are also equations of state describing solids, including the transition of solids from one crystalline state to another. There are equations that model the interior of stars, including neutron stars, dense matter (quark–gluon plasmas) and radiation fields. A related concept is the perfect fluid equation of state used in cosmology.
In practical context, the equations of state are instrumental for PVT calculation in process engineering problems and especially in petroleum gas/liquid equilibrium calculations. A successful PVT model based on a fitting equation of state can be helpful to determine the state of the flow regime, the parameters for handling the reservoir fluids, piping and sizing.
Historical
Boyle's law (1662)
Boyle's Law was perhaps the first expression of an equation of state. In 1662, the noted Irish physicist and chemist Robert Boyle performed a series of experiments employing a Jshaped glass tube, which was sealed on one end. Mercury was added to the tube, trapping a fixed quantity of air in the short, sealed end of the tube. Then the volume of gas was carefully measured as additional mercury was added to the tube. The pressure of the gas could be determined by the difference between the mercury level in the short end of the tube and that in the long, open end. Through these experiments, Boyle noted that the gas volume varied inversely with the pressure. In mathematical form, this can be stated as:
 pV = \mathrm{constant}.\,\!
The above relationship has also been attributed to Edme Mariotte and is sometimes referred to as Mariotte's law. However, Mariotte's work was not published until 1676.
Charles's law or Law of Charles and GayLussac (1787)
In 1787 the French physicist Jacques Charles found that oxygen, nitrogen, hydrogen, carbon dioxide, and air expand to the same extent over the same 80 kelvin interval. Later, in 1802, Joseph Louis GayLussac published results of similar experiments, indicating a linear relationship between volume and temperature:
 \frac{V_1}{T_1} = \frac{V_2}{T_2}.
Dalton's law of partial pressures (1801)
Dalton's Law of partial pressure states that the pressure of a mixture of gases is equal to the sum of the pressures of all of the constituent gases alone.
Mathematically, this can be represented for n species as:

 p_\text{total} = p_1+p_2+\cdots+p_n = p_\text{total} = \sum_{i=1}^n p_i.
The ideal gas law (1834)
In 1834 Émile Clapeyron combined Boyle's Law and Charles' law into the first statement of the ideal gas law. Initially the law was formulated as pV_{m} = R(T_{C} + 267) (with temperature expressed in degrees Celsius), where R is the gas constant. However, later work revealed that the number should actually be closer to 273.2, and then the Celsius scale was defined with 0 °C = 273.15 K, giving:
 \ pV_m = R (T_C+273.15).
Van der Waals equation of state (1873)
In 1873, J. D. van der Waals introduced the first equation of state derived by the assumption of a finite volume occupied by the constituent molecules.^{[2]} His new formula revolutionized the study of equations of state, and was most famously continued via the RedlichKwong equation of state and the Soave modification of RedlichKwong.
Major equations of state
For a given amount of substance contained in a system, the temperature, volume, and pressure are not independent quantities; they are connected by a relationship of the general form:
 {\ f(p,V,T) = 0}.
In the following equations the variables are defined as follows. Any consistent set of units may be used, although SI units are preferred. Absolute temperature refers to use of the Kelvin (K) or Rankine (°R) temperature scales, with zero being absolute zero.
 \ p = pressure (absolute)
 \ V = volume
 \ n = number of moles of a substance
 \ V_m = \frac{V}{n} = molar volume, the volume of 1 mole of gas or liquid
 \ T = absolute temperature
 \ R = ideal gas constant (8.3144621 J/(mol·K))
 \ p_c = pressure at the critical point
 \ V_c = molar volume at the critical point
 \ T_c = absolute temperature at the critical point
Classical ideal gas law
The classical ideal gas law may be written:
 {\ pV = nRT}.
In the form show above, the equation of state is thus
f(p, V, T) = pV  nRT = 0.
The ideal gas law may also be expressed as follows
 {\ p = \rho (\gamma1)e}
where \rho is the density, \gamma = C_p/C_v is the adiabatic index (ratio of specific heats), e=C_vT is the internal energy per unit mass (the "specific internal energy"), C_v is the specific heat at constant volume, and C_p is the specific heat at constant pressure.
Cubic equations of state
Cubic equations of state are called such because they can be rewritten as a cubic function of V_{m}.
Van der Waals equation of state
The Van der Waals equation of state may be written:
 {\left(p + \frac{a}{V_m^2}\right)\left(V_mb\right) = RT}
where V_m is molar volume. The substancespecific constants a and b can be calculated from the critical properties p_c, T_c and V_c (noting that V_c is the molar volume at the critical point) as:
 a = 3p_c \,V_c^2
 b = \frac{V_c}{3}.
Also written as
 a = \frac{27(R\,T_c)^2}{64p_c}
 b = \frac{R\,T_c}{8p_c}.
Proposed in 1873, the van der Waals equation of state was one of the first to perform markedly better than the ideal gas law. In this landmark equation a is called the attraction parameter and b the repulsion parameter or the effective molecular volume. While the equation is definitely superior to the ideal gas law and does predict the formation of a liquid phase, the agreement with experimental data is limited for conditions where the liquid forms. While the van der Waals equation is commonly referenced in textbooks and papers for historical reasons, it is now obsolete. Other modern equations of only slightly greater complexity are much more accurate.
The van der Waals equation may be considered as the ideal gas law, "improved" due to two independent reasons:
 Molecules are thought as particles with volume, not material points. Thus V_m cannot be too little, less than some constant. So we get (V_m  b) instead of V_m.
 While ideal gas molecules do not interact, we consider molecules attracting others within a distance of several molecules' radii. It makes no effect inside the material, but surface molecules are attracted into the material from the surface. We see this as diminishing of pressure on the outer shell (which is used in the ideal gas law), so we write (p + something) instead of p. To evaluate this ‘something’, let's examine an additional force acting on an element of gas surface. While the force acting on each surface molecule is ~\rho, the force acting on the whole element is ~\rho^2~\frac{1}{V_m^2}.
With the reduced state variables, i.e. V_{r}=V_{m}/V_{c}, P_{r}=P/P_{c} and T_{r}=T/T_{c}, the reduced form of the Van der Waals equation can be formulated:
 {\left(P_r + \frac{3}{V_r^2}\right)\left(3V_r1\right) = 8T_r}
The benefit of this form is that for given T_{r} and P_{r}, the reduced volume of the liquid and gas can be calculated directly using Cardano's method for the reduced cubic form:
 {V_r^3\left(\frac{1}{3}+ \frac{8T_r}{3P_r}\right)V_r^2+\frac{3V_r}{P_r}\frac{1}{P_r}=0}
For P_{r}<1 and T_{r}<1, the system is in a state of vapor–liquid equilibrium. The reduced cubic equation of state yields in that case 3 solutions. The largest and the lowest solution are the gas and liquid reduced volume.
RedlichKwong equation of state
 {p = \frac{R\,T}{V_mb}  \frac{a}{\sqrt{T}\,V_m\left(V_m+b\right)}}
 a = \frac{0.42748\,R^2\,T_c^{\,2.5}}{p_c}
 b = \frac{0.08662\,R\,T_c}{p_c}
Introduced in 1949, the RedlichKwong equation of state was a considerable improvement over other equations of the time. It is still of interest primarily due to its relatively simple form. While superior to the van der Waals equation of state, it performs poorly with respect to the liquid phase and thus cannot be used for accurately calculating vapor–liquid equilibria. However, it can be used in conjunction with separate liquidphase correlations for this purpose.
The RedlichKwong equation is adequate for calculation of gas phase properties when the ratio of the pressure to the critical pressure (reduced pressure) is less than about onehalf of the ratio of the temperature to the critical temperature (reduced temperature):
 \frac{p}{p_c} < \frac{T}{2T_c}.
Soave modification of RedlichKwong
 p = \frac{R\,T}{V_mb}  \frac{a\,\alpha}{V_m\left(V_m+b\right)}
 a = \frac{0.427\,R^2\,T_c^2}{P_c}
 b = \frac{0.08664\,R\,T_c}{P_c}
 \alpha = \left(1 + \left(0.48508 + 1.55171\,\omega  0.15613\,\omega^2\right) \left(1T_r^{\,0.5}\right)\right)^2
 T_r = \frac{T}{T_c}
Where ω is the acentric factor for the species.
This formulation for \alpha is due to Graboski and Daubert. The original formulation from Soave is:
 \alpha = \left(1 + \left(0.48 + 1.574\,\omega  0.176\,\omega^2\right) \left(1T_r^{\,0.5}\right)\right)^2
for hydrogen:
 \alpha = 1.202 \exp\left(0.30288\,T_r\right).
We can also write it in the polynomial form, with:
 A = \frac{a\,\alpha\,P}{R^2\,T^2}
 B = \frac{b\,P}{R\,T}
then we have:
 0 = Z^3Z^2+Z\left(ABB^2\right)AB \;,
where R is the universal gas constant and Z=PV/(RT) is the compressibility factor.
In 1972 G. Soave^{[3]} replaced the 1/√(T) term of the RedlichKwong equation with a function α(T,ω) involving the temperature and the acentric factor (the resulting equation is also known as the SoaveRedlichKwong equation). The α function was devised to fit the vapor pressure data of hydrocarbons and the equation does fairly well for these materials.
Note especially that this replacement changes the definition of a slightly, as the T_c is now to the second power.
PengRobinson equation of state
 p=\frac{R\,T}{V_mb}  \frac{a\,\alpha}{V_m^2+2bV_mb^2}
 a = \frac{0.457235\,R^2\,T_c^2}{p_c}
 b = \frac{0.077796\,R\,T_c}{p_c}
 \alpha = \left(1 + \kappa \left(1T_r^{\,0.5}\right)\right)^2
 \kappa = 0.37464 + 1.54226\,\omega  0.26992\,\omega^2
 T_r = \frac{T}{T_c}
In polynomial form:
 A = \frac{a\alpha p}{ R^2\,T^2}
 B = \frac{b p}{RT}
 Z^3  (1B)\ Z^2 + (A2B3B^2)\ Z (ABB^2B^3) = 0 \;,
where \omega is the acentric factor of the species, R is the universal gas constant and Z=PV/(RT) is compressibility factor.
The PengRobinson equation was developed in 1976 at The University of Alberta in order to satisfy the following goals:^{[4]}
 The parameters should be expressible in terms of the critical properties and the acentric factor.
 The model should provide reasonable accuracy near the critical point, particularly for calculations of the compressibility factor and liquid density.
 The mixing rules should not employ more than a single binary interaction parameter, which should be independent of temperature pressure and composition.
 The equation should be applicable to all calculations of all fluid properties in natural gas processes.
For the most part the PengRobinson equation exhibits performance similar to the Soave equation, although it is generally superior in predicting the liquid densities of many materials, especially nonpolar ones. The departure functions of the Peng–Robinson equation are given on a separate article.
PengRobinsonStryjekVera equations of state
PRSV1
A modification to the attraction term in the PengRobinson equation of state published by Stryjek and Vera in 1986 (PRSV) significantly improved the model's accuracy by introducing an adjustable pure component parameter and by modifying the polynomial fit of the acentric factor.^{[5]}
The modification is:
 \kappa = \kappa_0 + \kappa_1 \left(1+T_r^{0.5}\right) \left(0.7T_r\right)
 \kappa_0 = 0.378893+1.4897153\,\omega0.17131848\,\omega^2+0.0196554\,\omega^3
where \,\kappa_1 is an adjustable pure component parameter. Stryjek and Vera published pure component parameters for many compounds of industrial interest in their original journal article. At reduced temperatures above 0.7, they recommend to set \,\kappa_1 = 0 and simply use \kappa = \kappa_0 . For alcohols and water the value of \,\kappa_1 may be used up to the critical temperature and set to zero at higher temperatures.^{[5]}
PRSV2
A subsequent modification published in 1986 (PRSV2) further improved the model's accuracy by introducing two additional pure component parameters to the previous attraction term modification.^{[6]}
The modification is:
 \kappa = \kappa_0 + \left[\kappa_1+\kappa_2\left(\kappa_3T_r\right)\left(1T_r^{0.5}\right)\right]\left(1+T_r^{0.5}\right) \left(0.7T_r\right)
 \kappa_0 = 0.378893+1.4897153\,\omega0.17131848\,\omega^2+0.0196554\,\omega^3
where \,\kappa_1, \,\kappa_2, and \,\kappa_3 are adjustable pure component parameters.
PRSV2 is particularly advantageous for VLE calculations. While PRSV1 does offer an advantage over the Peng–Robinson model for describing thermodynamic behavior, it is still not accurate enough, in general, for phase equilibrium calculations.^{[5]} The highly nonlinear behavior of phaseequilibrium calculation methods tends to amplify what would otherwise be acceptably small errors. It is therefore recommended that PRSV2 be used for equilibrium calculations when applying these models to a design. However, once the equilibrium state has been determined, the phase specific thermodynamic values at equilibrium may be determined by one of several simpler models with a reasonable degree of accuracy.^{[6]}
One thing to note is that in the PSRV equation, the parameter fit is done in a particular temperature range which is usually below the critical temperature. Above the critical temperature, the PRSV alpha function tends to diverge and become arbitrarily large instead of tending towards 0. Because of this, alternate equations for alpha should be employed above the critical point. This is especially important for systems containing hydrogen which is often found at temperatures far above its critical point. Several alternate formulations have been proposed. Some well known ones are by Twu et all or by Mathias and Copeman.
Elliott, Suresh, Donohue equation of state
The Elliott, Suresh, and Donohue (ESD) equation of state was proposed in 1990.^{[7]} The equation seeks to correct a shortcoming in the PengRobinson EOS in that there was an inaccuracy in the van der Waals repulsive term. The EOS accounts for the effect of the shape of a nonpolar molecule and can be extended to polymers with the addition of an extra term (not shown). The EOS itself was developed through modeling computer simulations and should capture the essential physics of the size, shape, and hydrogen bonding.
 \frac{p V_m}{RT}=Z=1 + Z^{\rm{rep}} + Z^{\rm{att}}
where:
 Z^{\rm{rep}} = \frac{4 c \eta}{11.9 \eta}
 Z^{\rm{att}} = \frac{z_m q \eta Y}{1+ k_1 \eta Y}
and
 c is a "shape factor", with c=1 for spherical molecules
 For nonspherical molecules, the following relation is suggested:
 c=1+3.535\omega+0.533\omega^2 where \omega is the acentric factor
 The reduced number density \eta is defined as \eta=\frac{v^* n}{V}
where
 v^* is the characteristic size parameter
 n is the number of molecules
 V is the volume of the container
The characteristic size parameter is related to the shape parameter c through
 v^*=\frac{kT_c}{P_c}\Phi
where
 \Phi=\frac{0.0312+0.087(c1)+0.008(c1)^2}{1.000+2.455(c1)+0.732(c1)^2} and k is Boltzmann's constant.
Noting the relationships between Boltzmann's constant and the Universal gas constant, and observing that the number of molecules can be expressed in terms of Avogadro's number and the molar mass, the reduced number density \eta can be expressed in terms of the molar volume as
 \eta=\frac{R T_c}{P_c}\Phi\frac{1}{V_m}.
The shape parameter q appearing in the Attraction term and the term Y are given by
 q=1+k_3(c1) (and is hence also equal to 1 for spherical molecules).
 Y=\exp\left(\frac{\epsilon}{kT}\right)  k_2
where \epsilon is the depth of the squarewell potential and is given by
 \frac{\epsilon}{k} =\frac{1.000+0.945(c1)+0.134(c1)^2}{1.023+2.225(c1)+0.478(c1)^2}
 z_m, k_1, k_2 and k_3 are constants in the equation of state:
 z_m = 9.49 for spherical molecules (c=1)
 k_1 = 1.7745 for spherical molecules (c=1)
 k_2 = 1.0617 for spherical molecules (c=1)
 k_3 = 1.90476.
The model can be extended to associating components and mixtures of nonassociating components. Details are in the paper by J.R. Elliott, Jr. et al. (1990).^{[7]}
Noncubic equations of state
Dieterici equation of state
 \ p(Vb)=RTe^{a/RTV}
where a is associated with the interaction between molecules and b takes into account the finite size of the molecules, similar to the Van der Waals equation.
The reduced coordinates are:
 \ T_{c} = \frac{a}{4Rb}, \ p_{c} = \frac{a}{4b^{2}e^{2}}, \ V_{c} = 2b.
Virial equations of state
Virial equation of state
 \frac{pV_m}{RT} = 1 + \frac{B}{V_m} + \frac{C}{V_m^2} + \frac{D}{V_m^3} + \dots
 B = V_c \,
 C = \frac{V_c^2}{9}
Although usually not the most convenient equation of state, the virial equation is important because it can be derived directly from statistical mechanics. This equation is also called the Kamerlingh Onnes equation. If appropriate assumptions are made about the mathematical form of intermolecular forces, theoretical expressions can be developed for each of the coefficients. In this case B corresponds to interactions between pairs of molecules, C to triplets, and so on. Accuracy can be increased indefinitely by considering higher order terms. The coefficients B, C, D, etc. are functions of temperature only.
It can also be used to work out the Boyle Temperature (the temperature at which B = 0 and ideal gas laws apply) from a and b from the Van der Waals equation of state, if you use the value for B shown below:
 B = b  \frac{a}{RT}.
The BWR equation of state
 p=\rho RT + \left(B_0 RTA_0  \frac{C_0}{T^2} + \frac{D_0}{T^3}  \frac{E_0}{T^4}\right) \rho^2 + \left(bRTa\frac{d}{T}\right) \rho^3 + \alpha\left(a+\frac{d}{T}\right) \rho^6 + \frac{c\rho^3}{T^2}\left(1 + \gamma\rho^2\right)\exp\left(\gamma\rho^2\right)
where
 p = pressure
 ρ = the molar density
Values of the various parameters for 15 substances can be found in K.E. Starling (1973). Fluid Properties for Light Petroleum Systems.
Multiparameter equations of state
Helmholtz Function form
Multiparameter equations of state (MEOS) can be used to represent pure fluids with high accuracy, in both the liquid and gaseous states. MEOS's represent the Helmholtz function of the fluid as the sum of ideal gas and residual terms. Both terms are explicit in reduced temperature and reduced density  thus:
\frac{a(T,\rho)}{RT} = \frac{a^o(T,\rho)+a^r(T,\rho)}{RT}= \alpha^o(\tau,\delta)+ \alpha^r(\tau,\delta)
Where:
\tau = \frac{T_r}{T}, \delta = \frac{\rho}{\rho_r}
The reduced density and temperature are typically, though not always, the critical values for the pure fluid. Other thermodynamic functions can be derived from the MEOS by using appropriate derivatives of the Helmholtz function; hence, because integration of the MEOS is not required, there are few restrictions as to the functional form of the ideal or residual terms. Typical MEOS use upwards of 50 fluid specific parameters, but are able to represent the fluid's properties with high accuracy. MEOS are available currently for about 50 of the most common industrial fluids including refrigerants. Mixture models also exist.
Other equations of state of interest
Stiffened equation of state
When considering water under very high pressures (typical applications are underwater nuclear explosions, sonic shock lithotripsy, and sonoluminescence) the stiffened equation of state is often used:
 p=\rho(\gamma1)e\gamma p^0 \,
where e is the internal energy per unit mass, \gamma is an empirically determined constant typically taken to be about 6.1, and p^0 is another constant, representing the molecular attraction between water molecules. The magnitude of the correction is about 2 gigapascals (20,000 atmospheres).
The equation is stated in this form because the speed of sound in water is given by c^2=\gamma(p+p^0)/\rho.
Thus water behaves as though it is an ideal gas that is already under about 20,000 atmospheres (2 GPa) pressure, and explains why water is commonly assumed to be incompressible: when the external pressure changes from 1 atmosphere to 2 atmospheres (100 kPa to 200 kPa), the water behaves as an ideal gas would when changing from 20,001 to 20,002 atmospheres (2000.1 MPa to 2000.2 MPa).
This equation mispredicts the specific heat capacity of water but few simple alternatives are available for severely nonisentropic processes such as strong shocks.
Ultrarelativistic equation of state
An ultrarelativistic fluid has equation of state
 p=\rho_m c_s^2
where p is the pressure, \rho_m is the mass density, and c_s is the speed of sound.
Ideal Bose equation of state
The equation of state for an ideal Bose gas is
 pV_m=RT~\frac{\textrm{Li}_{\alpha+1}(z)}{\zeta(\alpha)} \left(\frac{T}{T_c}\right)^\alpha
where α is an exponent specific to the system (e.g. in the absence of a potential field, α=3/2), z is exp(μ/kT) where μ is the chemical potential, Li is the polylogarithm, ζ is the Riemann zeta function, and T_{c} is the critical temperature at which a Bose–Einstein condensate begins to form.
JonesWilkinsLee equation of state for explosives (JWLequation)
The equation of state from JonesWilkinsLee is used to describe the detonation products of explosives.
 p=A \cdot \left( 1  \frac{\omega}{R_1 \cdot V} \right) \cdot \exp (R_1 \cdot V) + B \cdot \left( 1  \frac{\omega}{R_2 \cdot V} \right) \cdot \exp (R_2 \cdot V) + \frac{\omega \cdot e_0}{V}
Material  \rho_0\, (g/cm^{3})  v_{D}\, (m/s)  p_{CJ}\, (GPa)  A\, (GPa)  B\, (GPa)  R_1\,  R_2\,  \omega\,  e_0\, (GPa) 

TNT  1.630  6930  21.0  373.8  3.747  4.15  0.90  0.35  6.00 
Composition B  1.717  7980  29.5  524.2  7.678  4.20  1.10  0.35  8.50 
PBX 9501^{[9]}  1.844  36.3  852.4  18.02  4.55  1.3  0.38  10.2 
Equations of state for solids
 p(V)=3K_0\left(\frac{1\eta}{\eta^2}\right)\exp\left[\tfrac{3}{2}\left(K_0^\prime1\right)(1\eta)\right]~,~~ \eta:=(V/V_0)^{1/3}~,~~ K_0^\prime := \frac{dK_0}{dp}
 p(V) = \frac{K_0}{K_0'} \left[\left(\frac{V}{V_0}\right)^{K_0'}  1\right] \,.
 p(V)=\frac{3K_0}{2} \left[\left(\frac{V_0}{V}\right)^\frac{7}{3}  \left(\frac{V_0}{V}\right)^\frac{5}{3}\right] \left\{1+\frac{3}{4}\left(K_0^\prime4\right) \left[\left(\frac{V_0}{V}\right)^\frac{2}{3}  1\right]\right\}.
 p(V) = \begin{cases} k_1~\xi + k_2~\xi^2 + k_3~\xi^3 + \Delta p & \qquad \text{Compression} \\ k_1~\xi & \qquad \text{Tension} \end{cases} ~;~~ \xi := \cfrac{V_0}{V}1
 p(V)  p_0 = \frac{\Gamma}{V} (e  e_0)
 p(V) =  \beta \left(\frac{V}{V_0}\right)^n \ln\left(\frac{V}{V_0}\right)
where \beta = K_0 is the bulk modulus at equilibrium volume V_0 and n is about 2 and related to the Grüneisen parameter n = \frac{1}{6}  \gamma_G
See also
References
 ^ Perrot, Pierre (1998). A to Z of Thermodynamics. Oxford University Press.
 ^ van der Waals, J. D. (1873). On the Continuity of the Gaseous and Liquid States (doctoral dissertation). Universiteit Leiden.
 ^ Soave, G. Equilibrium Constants from a Modified RedlichKwong Equation of State, Chem. Eng. Sci.,1 9 7 2, 27, 11971203
 ^ Peng, D. Y., and Robinson, D. B. (1976). "A New TwoConstant Equation of State". Industrial and Engineering Chemistry: Fundamentals 15: 59–64.
 ^ ^{}a ^{b} ^{c} Stryjek, R. and Vera, J. H. (1986). "PRSV: An improved PengRobinson equation of state for pure compounds and mixtures". The Canadian Journal of Chemical Engineering 64: 323–333.
 ^ ^{}a ^{b} Stryjek, R. and Vera, J. H. (1986). "PRSV2: A cubic equation of state for accurate vapor—liquid equilibria calculations". The Canadian Journal of Chemical Engineering 64: 820–826.
 ^ ^{}a ^{b} J. Richard, Jr. Elliott, S. Jayaraman Suresh, Marc D. Donohue (1990). "A Simple Equation of State for Nonspherical and Associating Molecules". Ind. Eng. Chem. Res. 29 (7): 1476–1485.
 ^ B. M. Dobratz, P. C. Crawford (1985). LLNL Explosives Handbook: Properties of Chemical Explosives and Explosive Simulants. University of California; Lawrence Livermore National Laboratory; Report UCRL5299; Rev.2;.
 ^ Wilkins, Mark L. (1999), Computer Simulation of Dynamic Phenomena, Springer, p. 80
Bibliography
 Elliot & Lira, (1999). Introductory Chemical Engineering Thermodynamics, Prentice Hall.

