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 of 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_b \, = \, 198.2 + \sum {G_i} }[/math]

This simple form assumes that the property (normal boiling point in the example) is strictly linear dependent from 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 pure 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 than 3/2.

[math]\displaystyle{ T_c \, = \, T_b \left[0.584 + 0.965 \sum {G_i} - {G_i}^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 ten parameters for ten groups a group interaction model needs already 45 parameters. Therefore, a group interaction model has normally not parameter for all possible combinations.

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 formation (ΔfG′°) and reaction (ΔrG′°) in biochemical systems: aqueous solution, temperature of 25℃ 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 the

  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=f(ΣGi2)
  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.

See also

References

  1. http://www.cas.org/newsevents/releases/research120810.html
  2. Constantinou L., Gani R., "New Group Contribution Method for Estimating Properties of Pure Compounds", AIChE J., 40(10), 1697-1710, 1994
  3. Nannoolal Y., Rarey J., Ramjugernath J., "Estimation of pure component properties Part 2. Estimation of critical property data by group contribution", Fluid Phase Equilib., 252(1-2), 1-27, 2007
  4. Jankowski, M. D., C. S. Henry,L. J. Broadbelt, and V. Hatzimanikatis. Group Contribution Method for Thermodynamic Analysis of Complex Metabolic Networks. Biophys. J. 95(3): 1487-1499, 2008
  5. Mavrovouniotis, M. L. Estimation of standard Gibbs energy changes of biotransformations. J. Biol. Chem. 266:14440-14445, 1991
  6. Mavrovouniotis, M. L. Group contributions for estimating standard Gibbs energies of formation of biochemical-compounds in aqueous-solution. Biotechnol. Bioeng. 36:1070-1082, 1990
  7. http://lcsbweb.epfl.ch/GCMWebSite/