Sigma-Profiles for COSMO-RS Type Models

The DDB σ Profile Data Bank DDB-Sigma

For the use with COSMO-RS type models like e.g. COSMO-RS, COSMO-SAC and COSMO-RS(Ol), a large data bank with the required σ profiles has been prepared. As the quantumchemical calculations become very demanding for large molecules, group contribution methods have been developed to predict σ profiles. These complimentary methods are delivered with Predict Pure and lead to reliable results for not too complex or multifunctional structures.

After careful analysis of the results from different modern DFT functionals (BP, B3LYP) with different basis sets and of different commercial products (DMOL3, Turbomole, Gaussian) it was decided to use the software package Gaussian 03 with the basis set 6-311G(d,p) as this provided the lowest overall deviation and the most reliable results (less large deviations from experimental data). Unluckily calculation times are much longer for this hybrid functional compared to the BP functional that can be accelerated by the RI-method by a factor of 10-20.

σ profiles from DDB-Sigma can be used within the DDB software and in most third party software (in COSMOTherm® a specially optimized parametrization for these profiles is available).

A brief description of the theoretical basis of COSMO-RS type models is given below. For more detailed information, please consult the literature given at the bottom of the page.


Short introduction to COSMO-RS type models

As a possible alternative to claasical predictive methods based on the group contribution approach, COSMO-RS type models originally developed by A. Klamt have attracted some interest in the recent years.

The model uses the results of a so called COSMO (COnductor-like Screening MOdel)-calculation, in which a molecular shaped cavity is constructed around the molecule and a large number of Coulomb charges are placed on this surface. The individual charges and the structure and charge distribution of the molecule are now optimized to find the minimum energy of the system. This can be achieved by different empirical, semiempirical or ab-initio methods. COSMO therefore allows to approximately calculate the electrostatic interaction of a molecule with a solvent.

In most cases, quantumchemical DFT methods are employed as these currently represent a good compromise between computational effort and reliability of the result.

In case of flexible molecules or components with different tautomeric forms, the calculation needs to be repeated for different conformers and/or tautomers.

In order to be used for COSMO-RS type models, an infinite dielectric constant is used to approximate the solvation of the molecule in an ideal conductor.

From the optimized charges on the molecular cavity, a list of charged surface segments can be generated. As the shielding charge density  (SCD) on the molecular surface generated by this procedure is rather detailed but contacting molecules are only able to "feel" the mean SCD difference over the whole contact area, some averaging of the SCD is required.

Based on this information, the electrostatic misfit energy can now be calculated for each possible contact between two molecules in a COSMO-RS (COSMO for Real Solvents) calculation by a simple formula. An empirical terms has been added to account for hydrogen bonding. As the shielding charges are now replaced by the space-distributed charges of the two molecules, both molecules should in principal change their conformation and charge distribution but this effect is neglected to simplify the procedure.

Instead of calculating all different mutual contacts of the three-dimensional representations of the molecules, a further simplification is introduced. From the SCD values on the molecular surface, a SCD distribution is calculated that represents the probability to find a certain SCD on the cavity surface. In practical implementations, the SCD distribution or σ profile is mostly represented by the probability values or surface areas for a number of discrete SCD intervals (e.g. 50 or 80 intervals). In this step, all structural information (relative position of the segments towards each other) is lost.

It is now assumed, that each segment can get into contact with each other segment and the probabilty and energy of all segment contacts in equilibrium is calculated using a self consistent set of equations. This leads to segment activity coefficients in the pure component and mixture. In the same way as in UNIFAC or mod. UNIFAC, the residual activity coefficients can now be calculated from the segment (or group in UNIFAC) activity coefficients. A combinatorial contribution is calculated from the cavity surfaces and volumes in a similar manner as in UNIFAC. An overview on the calculation procedure is given in the following figure:

In UNIFAC or mod. UNIFAC, the groups represent unique entities with specific chemical behavior while in COSMO-RS, the molecule is broken down into segments which only differ by their SCD. UNIFAC and mod. UNIFAC require regressed interaction parameters between the groups while in case of COSMO-RS, empirically modified electrostatic interactions and a hydrogen bonding contribution are used.

While this allows to calculate the real mixture behavior based on the sigma profiles of the components only, there is also a serious disadvantage. In contrast to the group contribution methods it is not easily possible to "train" the method by adjusting a larger number of parameters to the large amount of readily available information. The model thus has to rely on the asumption that the complex behavior of liquid mixtures can be calculated from only a few simple and universal equations and empirical constants.

Very detailed evaluations of the COSMO-RS flavors COSMO-SAC and COSMO-RS(Ol) have been presented by Mu, Rarey and Gmehling (see literature below).


Scientific Papers

2009 Group Contribution Prediction of Surface Charge Density Distribution of Molecules for COSMO-SAC    Mu T., Rarey J., Gmehling J.    Journal AIChE J., 55, 12, 3298 3300 (2009)
2009 Application of COSMO-RS Type Models to the Prediction of Excess Enthalpies    Constantinescu D., Rarey J., Gmehling J.    Journal Ind.Eng.Chem.Res., 277, 1,8710-8725 (2009)
2008 Reply to: Comments on Performance of COSMO-RS with Sigma Profiles from Different Model Chemistries    Mu T., Rarey J., Gmehling    Journal Ind.Eng.Chem.Res., 47, 3, 989 989 (2008)
2008 Conductor-Like Screening Model for Real Solvents ( COSMO-RS)    Mu T., Gmehling J.    Journal Huaxue-jinzhan, 20, 10, 1487 1494 (2008)
2007 Performance of COSMO-RS with Sigma Profiles from Different Model Chemistries    Mu T., Rarey J., Gmehling J.    Journal Ind.Eng.Chem.Res., 46, 20, 6612 6629 (2007)
2007 Group Contribution Prediction of Surface Charge Density Profiles for COSMO-RS(Ol)    Mu T., Rarey J., Gmehling J.    Journal AIChE J., 53, 12, 3231 3240 (2007)
2005 Weiterentwicklung thermodynamischer Vorhersagemodelle durch Einsatz von quantenchemischen Methoden    Grensemann H.    PhD thesis    Thesis, Carl-von-Ossietzky-Universität Oldenburg, 1-235, 2005
2005 Rebuttal to the Comments of Andreas Klamt on "Performance of a COSMO-RS Model in Comparison to Classical Group Contribution Methods"    Grensemann H., Gmehling J.    Journal Ind.Eng.Chem.Res., 44, 17, 7043 7044 (2005)
2005 Systems with Ionic Liquids: Measurement of VLE and gamma-infinite Data and Prediction of their Thermodynamic Behavior using Original UNIFAC, mod. UNIFAC(Do) and COSMO-RS(Ol)    Kato R., Gmehling J.    Journal J.Chem.Thermodyn., 37, 6, 603 619 (2005)
2005 Performance of a Conductor-Like Screening Model for Real Solvents Model in Comparison to Classical Group Contribution Methods    Grensemann H., Gmehling J.    Journal Ind.Eng.Chem.Res., 44, 5, 1610 1624 (2005)