Thermodynamic databases for pure substances

Thermodynamic databases contain information about thermodynamic properties for substances, the most important being enthalpy, entropy, and Gibbs free energy. Numerical values of these thermodynamic properties are collected as tables or are calculated from thermodynamic datafiles. Data is expressed as temperature-dependent values for one mole of substance at the standard pressure of 101.325 kPa (1 atm), or 100 kPa (1 bar). Both of these definitions for the standard condition for pressure are in use.

Thermodynamic data

Thermodynamic data is usually presented as a table or chart of function values for one mole of a substance (or in the case of the steam tables, one kg). A thermodynamic datafile is a set of equation parameters from which the numerical data values can be calculated. Tables and datafiles are usually presented at a standard pressure of 1 bar or 1 atm, but in the case of steam and other industrially important gases, pressure may be included as a variable. Function values depend on the state of aggregation of the substance, which must be defined for the value to have any meaning. The state of aggregation for thermodynamic purposes is the standard state, sometimes called the reference state, and defined by specifying certain conditions. The normal standard state is commonly defined as the most stable physical form of the substance at the specified temperature and a pressure of 1 bar or 1 atm. However, since any non-normal condition could be chosen as a standard state, it must be defined in the context of use. A physical standard state is one that exists for a time sufficient to allow measurements of its properties. The most common physical standard state is one that is stable thermodynamically (i.e., the normal one). It has no tendency to transform into any other physical state. If a substance can exist but is not thermodynamically stable (for example, a supercooled liquid), it is called a metastable state. A non-physical standard state is one whose properties are obtained by extrapolation from a physical state (for example, a solid superheated above the normal melting point, or an ideal gas at a condition where the real gas is non-ideal). Metastable liquids and solids are important because some substances can persist and be used in that state indefinitely. Thermodynamic functions that refer to conditions in the normal standard state are designated with a small superscript °. The relationship between certain physical and thermodynamic properties may be described by an equation of state.

Enthalpy, heat content and heat capacity

It is very difficult to measure the absolute amount of any thermodynamic quantity involving the internal energy (e.g. enthalpy), since the internal energy of a substance can take many forms, each of which has its own typical temperature at which it begins to become important in thermodynamic reactions. It is therefore the change in these functions that is of most interest. The isobaric change in enthalpy H above the common reference temperature of 298.15 K (25 °C) is called the high temperature heat content, the sensible heat, or the relative high-temperature enthalpy, and called henceforth the heat content. Different databases designate this term in different ways; for example H_T-H₂₉₈, H°-H°₂₉₈, H°_T-H°₂₉₈ or H°-H°(T_r), where T_r means the reference temperature (usually 298.15 K, but abbreviated in heat content symbols as 298). All of these terms mean the molar heat content for a substance in its normal standard state above a reference temperature of 298.15 K. Data for gases is for the hypothetical ideal gas at the designated standard pressure. The SI unit for enthalpy is J/mol, and is a positive number above the reference temperature. The heat content has been measured and tabulated for virtually all known substances, and is commonly expressed as a polynomial function of temperature. The heat content of an ideal gas is independent of pressure (or volume), but the heat content of real gases varies with pressure, hence the need to define the state for the gas (real or ideal) and the pressure. Note that for some thermodynamic databases such as for steam, the reference temperature is 273.15 K (0 °C).

The heat capacity C is the ratio of heat added to the temperature increase. For an incremental isobaric addition of heat:

C_{P}(T)=\left\{\lim _{\Delta T\to 0}{\frac {\Delta H}{\Delta T}}\right\}=\left({\frac {\partial H}{\partial T}}\right)_{p}

C_p is therefore the slope of a plot of temperature vs. isobaric heat content (or the derivative of a temperature/heat content equation). The SI units for heat capacity are J/(mol·K).

Molar heat content of four substances in their designated states above 298.15 K and at 1 atm pressure. CaO(c) and Rh(c) are in their normal standard state of crystalline solid at all temperatures. S₂(g) is a non-physical state below about 882 K and NiO(g) is a non-physical state at all temperatures.

Molar heat capacity of four substances in their designated states at 1 atm pressure. CaO(c) and Rh(c) are in their normal standard state of crystalline solid at all temperatures. S₂(g) is a non-physical state below about 882 K and NiO(g) is a non-physical state at all temperatures.

Enthalpy change of phase transitions

When heat is added to a condensed-phase substance, its temperature increases until a phase change temperature is reached. With further addition of heat, the temperature remains constant while the phase transition takes place. The amount of substance that transforms is a function of the amount of heat added. After the transition is complete, adding more heat increases the temperature. In other words, the enthalpy of a substance changes isothermally as it undergoes a physical change. The enthalpy change resulting from a phase transition is designated ΔH. There are four types of enthalpy changes resulting from a phase transition. To wit:

Enthalpy of transformation. This applies to the transformations from one solid phase to another, such as the transformation from α-Fe (bcc ferrite) to $\gamma$ -Fe (fcc austenite). The transformation is designated ΔH_tr.
Enthalpy of fusion or melting. This applies to the transition of a solid to a liquid and is designated ΔH_m.
Enthalpy of vaporization. This applies to the transition of a liquid to a vapor and is designated ΔH_v.
Enthalpy of sublimation. This applies to the transition of a solid to a vapor and is designated ΔH_s.

C_p is infinite at phase transition temperatures because the enthalpy changes isothermally. At the Curie temperature, C_p shows a sharp discontinuity while the enthalpy has a change in slope.

Values of ΔH are usually given for the transition at the normal standard state temperature for the two states, and if so, are designated with a superscript °. ΔH for a phase transition is a weak function of temperature. In some texts, the heats of phase transitions are called latent heats (for example, latent heat of fusion).

Molar enthalpy of zinc above 298.15 K and at 1 atm pressure, showing discontinuities at the melting and boiling points. The ΔH°m of zinc is 7323 J/mol, and the ΔH°v is 115 330 J/mol.

Enthalpy change for a chemical reaction

An enthalpy change occurs during a chemical reaction. For the special case of the formation of a compound from the elements, the change is designated ΔH_form and is a weak function of temperature. Values of ΔH_form are usually given where the elements and compound are in their normal standard states, and as such are designated standard heats of formation, as designated by a superscript °. The ΔH°_form undergoes discontinuities at a phase transition temperatures of the constituent element(s) and the compound. The enthalpy change for any standard reaction is designated ΔH°_rx.

Standard molar heat of formation of ZnBr₂(c,l) from the elements, showing discontinuities at transition temperatures of the elements and the compound.

Entropy and Gibbs energy

The entropy of a system is another thermodynamic quantity that is not easily measured. However, using a combination of theoretical and experimental techniques, entropy can in fact be accurately estimated. At low temperatures, the Debye model leads to the result that the atomic heat capacity C_v for solids should be proportional to T³, and that for perfect crystalline solids it should become zero at absolute zero. Experimentally, the heat capacity is measured at temperature intervals to as low a temperature as possible. Values of C_p/T are plotted against T for the whole range of temperatures where the substance exists in the same physical state. The data are extrapolated from the lowest experimental temperature to 0 K using the Debye model. The third law of thermodynamics states that the entropy of a perfect crystalline substance becomes zero at 0 K. When S₀ is zero, the area under the curve from 0 K to any temperature gives the entropy at that temperature. Even though the Debye model contains C_v instead of C_p, the difference between the two at temperatures near 0 K is so small as to be negligible.

The absolute value of entropy for a substance in its standard state at the reference temperature of 298.15 K is designated S°₂₉₈. Entropy increases with temperature, and is discontinuous at phase transition temperatures. The change in entropy (ΔS°) at the normal phase transition temperature is equal to the heat of transition divided by the transition temperature. The SI units for entropy are J/(mol·K).

Absolute entropy of strontium. The solid line refers to the entropy of strontium in its normal standard state at 1 atm pressure. The dashed line refers to the entropy of strontium vapor in a non-physical state.

The standard entropy change for the formation of a compound from the elements, or for any standard reaction is designated ΔS°_form or ΔS°_rx. The entropy change is obtained by summing the absolute entropies of the products minus the sum of the absolute entropies of the reactants. Like enthalpy, the Gibbs energy G has no intrinsic value, so it is the change in G that is of interest. Furthermore, there is no change in G at phase transitions between substances in their standard states. Hence, the main functional application of Gibbs energy from a thermodynamic database is its change in value during the formation of a compound from the standard-state elements, or for any standard chemical reaction (ΔG°_form or ΔG°_rx). The SI units of Gibbs energy are the same as for enthalpy (J/mol).

Standard heat and Gibbs energy change for the reaction:

Pb(c,l)+2HCl(g)\Rightarrow PbCl_{2}+H_{2}(g)

The ΔH°_rx shows discontinuities at the melting points of Pb (600.65 K) and PbCl₂ (771 K). ΔG°_rx is not discontinuous at these phase transition temperatures, but does undergo a change in slope, which is almost imperceptible on the chart.

Additional functions

Compilers of thermochemical databases may contain some additional thermodynamic functions. For example, the absolute enthalpy of a substance H(T) is defined in terms of its formation enthalpy and its heat content as follows:

H(T)=\Delta H_{form,298}^{\circ }+[H_{T}-H_{298}]

For an element, H(T) and [H_T - H₂₉₈] are identical at all temperatures because ΔH°_form is zero, and of course at 298.15 K, H(T) = 0. For a compound:

\Delta H_{form}^{\circ }=H(T)compound-\sum \left\{H(T)elements\right\}

Similarly, the absolute Gibbs energy G(T) is defined by the absolute enthalpy and entropy of a substance:

G(T)=H(T)-T\times S(T)

For a compound:

\Delta G_{form}^{\circ }=G(T)compound-\sum \left\{G(T)elements\right\}

Some tables may also contain the Gibbs energy function (H°_298.15 – G°_T)/T which is defined in terms of the entropy and heat content.

(H_{298}^{\circ }-G_{T}^{\circ })/T=S_{T}^{\circ }-(H_{T}-H_{298})/T

The Gibbs energy function has the same units as entropy, but unlike entropy, exhibits no discontinuity at normal phase transition temperatures.

The log₁₀ of the equilibrium constant K_eq is often listed, which is calculated from the defining thermodynamic equation.

\log _{10}\left(K_{eq}\right)=-\Delta G_{form}^{\circ }/(19.1448T)

Thermodynamic databases

A thermodynamic database consists of sets of critically evaluated values for the major thermodynamic functions. Originally, data was presented as printed tables at 1 atm and at certain temperatures, usually 100° intervals and at phase transition temperatures. Some compilations included polynomial equations that could be used to reproduce the tabular values. More recently, computerized databases are used which consist of the equation parameters and subroutines to calculate specific values at any temperature and prepare tables for printing. Computerized databases often include subroutines for calculating reaction properties and displaying the data as charts.

Thermodynamic data comes from many types of experiments, such as calorimetry, phase equilibria, spectroscopy, composition measurements of chemical equilibrium mixtures, and emf measurements of reversible reactions. A proper database takes all available information about the elements and compounds in the database, and assures that the presented results are internally consistent. Internal consistency requires that all values of the thermodynamic functions are correctly calculated by application of the appropriate thermodynamic equations. For example, values of the Gibbs energy obtained from high-temperature equilibrium emf methods must be identical to those calculated from calorimetric measurements of the enthalpy and entropy values. The database provider must use recognized data analysis procedures to resolve differences between data obtained by different types of experiments.

All thermodynamic data is a non-linear function of temperature (and pressure), but there is no universal equation format for expressing the various functions. Here we describe a commonly used polynomial equation to express the temperature dependence of the heat content. A common six-term equation for the isobaric heat content is:

H_{T}-H_{298}=A(T)+B(T^{2})+C(T^{-1})+D(T^{0.5})+E(T^{3})+F\,

Regardless of the equation format, the heat of formation of a compound at any temperature is ΔH°_form at 298.15 K, plus the sum of the heat content parameters of the products minus the sum of the heat content parameters of the reactants. The C_p equation is obtained by taking the derivative of the heat content equation.

C_{P}=A+2B(T)-C(T^{-2})+\textstyle {\frac {1}{2}}D(T^{-0.5})+3E(T^{2})\,

The entropy equation is obtained by integrating the C_p/T equation:

S_{T}^{\circ }=A(\ln T)+2B(T)+\textstyle {\frac {1}{2}}C(T^{-2})-D(T^{\textstyle -{\frac {1}{2}}})+1\textstyle {\frac {1}{2}}E(T^{2})+F'

F' is a constant of integration obtained by inserting S° at any temperature T. The Gibbs energy of formation of a compound is obtained from the defining equation ΔG°_form = ΔH°_form – T(ΔS°_form), and is expressed as

\Delta G_{form}^{\circ }=(\Delta A-\Delta F')T-\Delta A(T\ln T)-\Delta B(T^{2})+\textstyle {\frac {1}{2}}\Delta C(T^{-1})+2\Delta D(T^{\textstyle {\frac {1}{2}}})

-\textstyle {\frac {1}{2}}\Delta E(T^{3})+\Delta F+\Delta H_{form298}^{\circ }

For most substances, ΔG°_form deviates only slightly from linearity with temperature, so over a short temperature span, the seven-term equation can be replaced by a three-term equation, whose parameter values are obtained by regression of tabular values.

\Delta G_{form}^{\circ }=\alpha T+\beta (T\ln T)+\chi

Depending on the accuracy of the data and the length of the temperature span, the heat content equation may require more or fewer terms. Over a very long temperature span, two equations may be used instead of one. It is unwise to extrapolate the equations to obtain values outside the range of experimental data used to derive the equation parameters.

Thermodynamic datafiles

The equation parameters and all other information required to calculate values of the important thermodynamic functions are stored in a thermodynamic datafile. The values are organized in a format that makes them readable by a thermodynamic calculation program or for use in a spreadsheet. For example, the Excel-based thermodynamic database FREED [1] creates the following type of datafile, here for a standard pressure of 1 atm.

Thermodynamic datafile for MgCl₂(c,l,g) from FREED. Some values have truncated significant figures for display purposes. The explanation for the values is shown below.

Row 1. Molar mass of species, density at 298.15 K, ΔH°_{form 298.15}, S°_298.15. and the upper temperature limit for the file.
Row 2. Number of C_p equations required. Here, three because of three species phases.
Row 3. Values of the five parameters for the first C_p equation; temperature limit for the equation.
Row 4. Values of the five parameters for the second C_p equation; temperature limit for the equation.
Row 5. Values of the five parameters for the third C_p equation; temperature limit for the equation.
Row 6. Number of H_T - H₂₉₈ equations required.
Row 7. Values of the six parameters for the first H_T - H₂₉₈ equation; temperature limit for the equation, and ΔH°_trans for the first phase change.
Row 8. Values of the six parameters for the second H_T - H₂₉₈ equation; temperature limit for the equation, and ΔH°_trans for the second phase change.
Row 9. Values of the six parameters for the third H_T - H₂₉₈ equation; temperature limit for the equation, and ΔH°_trans for the third phase change.
Row 10. Number of ΔH°_form equations required. Here five; three for species phases and two because one of the elements has a phase change.
Row 11. Values of the six parameters for the first ΔH°_form equation; temperature limit for the equation.
Row 12. Values of the six parameters for the second ΔH°_form equation; temperature limit for the equation.
Row 13. Values of the six parameters for the third ΔH°_form equation; temperature limit for the equation.
Row 14. Values of the six parameters for the fourth ΔH°_form equation; temperature limit for the equation.
Row 15. Values of the six parameters for the fifth ΔH°_form equation; temperature limit for the equation.
Row 16. Number of ΔG°_form equations required.
Row 17. Values of the seven parameters for the first ΔG°_form equation; temperature limit for the equation.
Row 18. Values of the seven parameters for the second ΔG°_form equation; temperature limit for the equation.
Row 19. Values of the seven parameters for the third ΔG°_form equation; temperature limit for the equation.
Row 20. Values of the seven parameters for the fourth ΔG°_form equation; temperature limit for the equation.
Row 21. Values of the seven parameters for the fifth ΔG°_form equation; temperature limit for the equation.

Most computerized databases will create a table of thermodynamic values using the values from the datafile. For MgCl₂(c,l,g) at 1 atm pressure:

Thermodynamic properties table for MgCl₂(c,l,g), from the FREED datafile. Some values have truncated significant figures for display purposes.

The table format is a common way to display thermodynamic data. The FREED table gives additional information in the top rows, such as the mass and amount composition and transition temperatures of the constituent elements. Transition temperatures for the constituent elements have dashes ------- in the first column in a blank row, such as at 922 K, the melting point of Mg. Transition temperatures for the substance have two blank rows with dashes, and a center row with the defined transition and the enthalpy change, such as the melting point of MgCl₂ at 980 K. The datafile equations are at the bottom of the table, and the entire table is in an Excel worksheet. This is particularly useful when the data is intended for making specific calculations.

References

Barin, Ihsan (2004). Thermochemical Data of Pure Substances. Wiley-VCH. ISBN 3-527-30993-4.
Chase, M. W. (1998). NIST - JANAF Thermochemical Tables (Fourth ed.). Journal of Physical and Chemical Reference Data. ISBN 1-56396-831-2.
Cox, J. D.; Wagman, Donald D.; Medvedev, Vadim A. (1989). CODATA Key Values for Thermodynamics. John Benjamins Publishing Co. ISBN 0-89116-758-7.
Hummel, Wolfgang; Urs Berner; Enzo Curti; F. J. Pearson; Tres Thoenen (2002). Nagra/Psi Chemical Thermodynamic Data Base. Universal Publishers. ISBN 1-58112-620-4.
Lide, David R.; Henry V. Kehiaian (1994). CRC Handbook of Thermophysical and Thermochemical Data (book & disk ed.). Boca Raton: CRC Press. ISBN 0-8493-0197-1.
Pankratz, L. B. (1982). "Thermodynamic Properties of Elements and Oxides". U. S. Bureau of Mines Bulletin. 672.
Pankratz, L. B. (1984). "Thermodynamic Properties of Halides". U. S. Bureau of Mines Bulletin. 674.
Pankratz, L. B.; A. D. Mah; S. W. Watson (1987). "Thermodynamic Properties of Sulfides". U. S. Bureau of Mines Bulletin. 689.
Pankratz, L. B. (1994). "Thermodynamic Properties of Carbides, Nitrides, and Other Selected Substances". U. S. Bureau of Mines Bulletin. 696.
Robie, Richard A., and Bruce S. Hemingway (1995). Thermodynamic Properties of Minerals . . . at Higher Temperatures, U. S. Geological Survey Bulletin 2131.
Yaws, Carl L. (2007). Yaws Handbook of Thermodynamic Properties for Hydrocarbons & Chemicals, Gulf Publishing Company. ISBN 1-933762-07-1.
Gurvich, L.V., Veitz, I.V., et al. (1989) Thermodynamic Properties of Individual Substances. Fourth edition, Hemisphere Pub Co. NY, L., Vol.1 in 2 parts.

External links

NIST WebBook A gateway to the data collection of the National Institute of Standards and Technology.
NASA Glenn ThermoBuild A web interface to generate tabulated thermodynamic data.
Burcat's Thermodynamic Database Database for more than 3,000 chemical species.
DIPPR The Design Institute for Physical Properties
DIPPR 801 Critically evaluated thermophysical property database useful for chemical process design and equilibrium calculations.
Free Steam Tables Online calculator based on IAPWS-IF97
FACT-Web programs Various on-line tools for obtaining thermodynamic data and making equilibrium calculations.
Mol-Instincts A chemical database based on Quantum Mechanics and QSPR, providing thermodynamic properties for millions of compounds.