Physics:Group-contribution method

From HandWiki

A group-contribution method in chemistry is a technique to estimate and predict thermodynamic and other properties from molecular structures.

Introduction

In today's chemical processes hundreds of thousands of components are used. The Chemical Abstracts Service registry lists 56 million substances,[1] but many of these are only of scientific interest.

Process designers need to know some basic chemical properties of the components and their mixtures. Experimental measurement is often too expensive.

Predictive methods can replace measurements if they provide sufficiently good estimations. The estimated properties cannot be as precise as well-made measurements, but for many purposes the quality of estimated properties is sufficient. Predictive methods can also be used to check the results of experimental work.

Principles

Principle of a group-contribution method

A group-contribution method uses the principle that some simple aspects of the structures of chemical components are always the same in many different molecules. The smallest common constituents are the atoms and the bonds. The vast majority of organic components, for example, are built of carbon, hydrogen, oxygen, nitrogen, halogens, and maybe sulfur or phosphorus. Together with a single, a double, and a triple bond there are only ten atom types (not including astatine) and three bond types to build thousands of components. The next slightly more complex building blocks of components are functional groups, which are themselves built from few atoms and bonds.

A group-contribution method is used to predict properties of pure components and mixtures by using group or atom properties. This reduces the number of needed data dramatically. Instead of needing to know the properties of thousands or millions of compounds, only data for a few dozens or hundreds of groups have to be known.

Additive group-contribution method

The simplest form of a group-contribution method is the determination of a component property by summing up the group contribution:

[math]\displaystyle{ T_\text{b}[\text{K}] = 198.2022567824111 + \sum G_i. }[/math]

This simple form assumes that the property (normal boiling point in the example) is strictly linearly dependent on the number of groups, and additionally no interaction between groups and molecules are assumed. This simple approach is used, for example, in the Joback method for some properties, and it works well in a limited range of components and property ranges, but leads to quite large errors if used outside the applicable ranges.

Additive group contributions and correlations

This technique uses the purely additive group contributions to correlate the wanted property with an easy accessible property. This is often done for the critical temperature, where the Guldberg rule implies that Tc is 3/2 of the normal boiling point, and the group contributions are used to give a more precise value:

[math]\displaystyle{ T_\text{c} = T_\text{b} \left[0.584 + 0.965 \sum G_i - \left(\sum G_i\right)^2 \right]^{-1}. }[/math]

This approach often gives better results than pure additive equations because the relation with a known property introduces some knowledge about the molecule. Commonly used additional properties are the molecular weight, the number of atoms, chain length, and ring sizes and counts.

Group interactions

For the prediction of mixture properties it is in most cases not sufficient to use a purely additive method. Instead the property is determined from group-interaction parameters:

[math]\displaystyle{ P = f(G_{ij}), }[/math]

where P stands for property, and Gij for group-interaction value.

A typical group-contribution method using group-interaction values is the UNIFAC method, which estimates activity coefficients. A big disadvantage of the group-interaction model is the need for many more model parameters. Where a simple additive model only needs 10 parameters for 10 groups, a group-interaction model needs already 45 parameters. Therefore, a group-interaction model has normally not parameter for all possible combinations[clarify].

Group contributions of higher orders

Some newer methods[2] introduce second-order groups. These can be super-groups containing several first-order (standard) groups. This allows the introduction of new parameters for the position of groups. Another possibility is to modify first-order group contributions if specific other groups are also present.[3]

If the majority of group-contribution methods give results in gas phase, recently, a new such method[4] was created for estimating the standard Gibbs free energy of formationfG′°) and reaction (ΔrG′°) in biochemical systems: aqueous solution, temperature of 25 °C and pH = 7 (biochemical conditions). This new aqueous-system method is based on the group-contribution method of Mavrovouniotis.[5][6]

A free-access tool of this new method in aqueous condition is available on the web.[7]

Determination of group contributions

Group contributions are obtained from known experimental data of well defined pure components and mixtures. Common sources are thermophysical data banks like the Dortmund Data Bank, Beilstein database, or the DIPPR data bank (from AIChE). The given pure component and mixture properties are then assigned to the groups by statistical correlations like e. g. (multi-)linear regression.

Important steps during the development of a new method are:

  1. Evaluation of the quality of available experimental data, elimination of wrong data, finding of outliers.
  2. Construction of groups.
  3. Searching additional simple and easily accessible properties that can be used to correlate the sum of group contributions with the examined property.
  4. Finding a good but simple mathematical equation for the relation of the group contribution sum with the wanted property. The critical pressures, for example, is often determined as Pc = fGi2).
  5. Fitting the group contribution.

The reliability of a method mainly relies on a comprehensive data bank where sufficient source data have been available for all groups. A small data base may lead to a precise reproduction of the used data but will lead to significant errors when the model is used for the prediction of other systems.

Group contribution methods

Joback method

The Joback method was published in 1984 by Kevin G. Joback. It can be used to estimate critical temperature, critical pressure, critical volume, standard ideal gas enthalpy of formation, standard ideal gas Gibbs energy of formation, ideal gas heat capacity, enthalpy of vaporization, enthalpy of fusion, normal boiling point, freezing point, and liquid viscosity.[8] The Joback method is a first-order method, and does not account for molecular interactions.

Ambrose method

The Ambrose method was published by Douglas Ambrose in 1978 and 1979. It can be used to estimate critical temperature, critical pressure, and critical volume. In addition to the molecular structure, it requires normal boiling point for estimating critical temperature and molecular weight for estimating critical pressure.[9][10]

Nannoolal method

The Nannoolal method was published by Yash Nannoolal et al in 2004. It can be used to estimate the normal boiling point. It includes first-order and second-order contributions.[11]

See also

References

  1. "CAS Chemistry Research Report". http://www.cas.org/newsevents/releases/research120810.html. 
  2. Constantinou, Leonidas; Gani, Rafiqul (1994). "New group contribution method for estimating properties of pure compounds". AIChE Journal 40 (10): 1697–1710. doi:10.1002/aic.690401011. 
  3. Nannoolal, Yash; Rarey, Jürgen; Ramjugernath, Deresh (2007). "Estimation of pure component properties". Fluid Phase Equilibria 252 (1–2): 1–27. doi:10.1016/j.fluid.2006.11.014. 
  4. Jankowski, Matthew D.; Henry, Christopher S.; Broadbelt, Linda J.; Hatzimanikatis, Vassily (2008). "Group Contribution Method for Thermodynamic Analysis of Complex Metabolic Networks". Biophysical Journal 95 (3): 1487–1499. doi:10.1529/biophysj.107.124784. PMID 18645197. Bibcode2008BpJ....95.1487J. 
  5. Mavrovouniotis, M. L. (1991). "Estimation of standard Gibbs energy changes of biotransformations". The Journal of Biological Chemistry 266 (22): 14440–5. doi:10.1016/S0021-9258(18)98705-3. PMID 1860851. 
  6. Mavrovouniotis, Michael L. (1990). "Group contributions for estimating standard gibbs energies of formation of biochemical compounds in aqueous solution". Biotechnology and Bioengineering 36 (10): 1070–1082. doi:10.1002/bit.260361013. PMID 18595046. 
  7. "GCM Software". http://lcsbweb.epfl.ch/GCMWebSite/. 
  8. Joback, K. G. (1984). A Unified Approach to Physical Property Estimation Using Multivariate Statistical Techniques (PDF) (MS). Massachusetts Institute of Technology.
  9. Ambrose, D. (1978). "Correlation and Estimation of Vapour-Liquid Critical Properties. I. Critical Temperatures of Organic Compounds". National Physical Laboratory Report. 
  10. Ambrose, D. (1979). "Correlation and Estimation of Vapour-Liquid Critical Properties. II. Critical Pressures and Volumes of Organic Compounds". National Physical Laboratory Report. 
  11. Nannoolal, Yash; Rarey, Jürgen; Ramjugernath, Deresh; Cordes, Wilfried (10 December 2004). "Estimation of pure component properties: Part 1. Estimation of the normal boiling point of non-electrolyte organic compounds via group contributions and group interactions". Fluid Phase Equilibria 226: 45–63. doi:10.1016/j.fluid.2004.09.001.