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). Unfortunately, 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\; o\; 0\}frac\{Delta\; H\}\{Delta\; T\}\; ight\; \}=left\; (\; frac\{partial\; H\}\{partial\; T\}\; ight\; )\; \_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).

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").

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.

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).

The standard enthalpy 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).

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^circ\_\{form,298\}\; +\; [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^circ\_\{form\}\; =\; H(T)compound\; -\; sum\; left\; \{\; H(T)elements\; ight\; \}$

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

$G(T)\; =\; H(T)\; -\; T\; imes\; S(T)$

For a compound:

$Delta\; G^circ\_\{form\}\; =\; G(T)compound\; -\; sum\; left\; \{\; G(T)elements\; ight\; \}$

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^circ\_\{298\}\; -\; G^circ\_T)\; /\; T\; =\; S^circ\_T\; -\; (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\}\; ight\; )\; =\; -Delta\; G^circ\_\{form\}\; /(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\})\; +\; extstyle\; frac\; \{1\}\{2\}\; D(T^\{-0.5\})\; +\; 3E(T^2)\; ,$

The entropy equation is obtained by integrating the "C"_p/T equation:

$S^circ\_T\; =\; A(ln\; T)\; +\; 2B(T)\; +\; extstyle\; frac\; \{1\}\{2\}C(T^\{-2\})\; -\; D(T^\{\; extstyle\; -\; frac\; \{1\}\{2)\; +\; 1\; extstyle\; frac\; \{1\}\{2\}\; E(T^2)\; +\; F\text{'}$

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^circ\_\{form\}\; =\; (\; Delta\; A\; -\; Delta\; F\text{'}\; )T\; -\; Delta\; A\; (\; T\; ln\; T\; )\; -\; Delta\; B\; (\; T^2\; )\; +\; extstyle\; frac\; \{1\}\{2\}\; Delta\; C\; (\; T^\{-1\}\; )\; +\; 2\; Delta\; D\; (\; T^\{\; extstyle\; frac\; \{1\}\{2\}\; \}\; )$

$-\; extstyle\; frac\; \{1\}\{2\}\; Delta\; E\; (T^3)\; +\; Delta\; F\; +\; Delta\; H^circ\_\{form\; 298\}$

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.

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 [http://www.thermart.net] creates the following type of datafile, here for a standard pressure of 1 atm.

:* 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:

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.

ee also

* Chemical thermodynamics
* Physical chemistry
* Materials science
* Laws of thermodynamics
* Thermochemistry
* Standard temperature and pressure
* Dortmund Data Bank
* CALPHAD_(method)

References

*cite book|last = Barin|first=Ihsan|year=2004
title=Thermochemical Data of Pure Substances|publisher=Wiley-VCH
id=ISBN 3-527-30993-4
*cite book| last = Chase| first = M. W.
year = 1998| title = NIST - JANAF Thermochemical Tables
edition = Fourth edition| publisher = Journal of Physical and Chemical Reference Data
id =ISBN 1-56396-831-2| url = http://www.aip.org/pubs/books/jpcrd_books.html
*cite book| last = Cox| first=J. D.
coauthors=Wagman, D.D. and V. A. Medvedev, V. A.|year= 1989
title=CODATA Key Values for Thermodynamics|publisher=John Benjamins Publishing Co
id=ISBN 0-89116-758-7
*cite book| last = Hummel|first=Wolfgang
coauthors=Urs Berner, Enzo Curti, F. J. Pearson, and Tres Thoenen|year=2002
title=Nagra/Psi Chemical Thermodynamic Data Base|publisher=Universal Publishers
id=ISBN 1-58112-620-4
*cite book|last=Lide|first=David R.|coauthors=Henry V. Kehiaian|year=1994
title=CRC Handbook of Thermophysical and Thermochemical Data|edition=book & disk edition
id=ISBN 0-8493-0197-1
*cite journal|last=Pankratz|first=L. B.|year=1982
title=Thermodynamic Properties of Elements and Oxides
journal=U. S. Bureau of Mines Bulletin|volume=672
*cite journal|last=Pankratz|first=L. B.|year=1984
title=Thermodynamic Properties of Halides
journal=U. S. Bureau of Mines Bulletin|volume=674
*cite journal|last=Pankratz|first=L. B.|coauthors=A. D. Mah and S. W. Watson|year=1987
title=Thermodynamic Properties of Sulfides
journal=U. S. Bureau of Mines Bulletin|volume=689
*cite journal|last=Pankratz|first=L. B.|year=1994
title=Thermodynamic Properties of Carbides, Nitrides, and Other Selected Substances
journal=U. S. Bureau of Mines Bulletin|volume=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

* [http://webbook.nist.gov/ NIST WebBook] A gateway to the data collection of the National Institute of Standards and Technology.
* [http://www.thermodata.org THERMODATA] Thermochemical Databases and Softwares.
* [http://www.aiche.org/TechnicalSocieties/DIPPR/index.aspx DIPPR] The Design Institute for Physical Properties
* [http://dippr.byu.edu DIPPR 801] Critically evaluated thermophysical property database useful for chemical process design and equilibrium calculations.
* [http://www.mtdata-software.com MTDATA software and databases for calculation of thermodynamic properties and phase equilibria]
* [http://www.steamtablesonline.com Free Steam Tables Online] calculator based on IAPWS-IF97