Open access

The Use of Solvation Models in Gas Chromatography

Written By

Fabrice Mutelet

Submitted: 24 November 2011 Published: 24 October 2012

DOI: 10.5772/48381

From the Edited Volume

Chromatography - The Most Versatile Method of Chemical Analysis

Edited by Leonardo de Azevedo Calderon

Chapter metrics overview

3,673 Chapter Downloads

View Full Metrics

1. Introduction

Gas chromatography is widely used for determining thermodynamic properties of pure substances or solvent properties of binary mixtures. From retention data, the solute activity coefficient at infinite dilution, the gas–liquid partition coefficient and others thermodynamic properties of mixing can be easily obtained. Using these parameters and appropriate models allows understanding of the intermolecular interactions responsible for solvation in the stationary phase [1-2]. The solvation parameter model is now well established as a useful tool for obtaining quantitative structure–property relationships for chemical, biomedical and environmental processes. The model correlates a free-energy related property of a system to six free-energy derived descriptors describing molecular properties. The ultimate goal is to establish a suitable quantitative structure–property relationship (QSPR) to facilitate the prediction of further system properties for compounds lacking experimental values. Two broad strategies are generally employed in QSPR studies. The first approach is based on theoretical descriptors. The advantage of using the QSPR approach based on theoretical descriptors is that all of the necessary parameters for prediction can be calculated purely from the three-dimensional representation of the molecular structure of each of the compounds of the mixtures, including mixtures of chemically diverse compounds [3-4]. The main weakness of this approach is that the selected descriptors may be difficult to understand and the models may lack obvious chemical significance. The second approach is based on descriptors determined using experimental technique such as gas chromatography.

In this review, we will present the different possibility of using chromatographic methods to facilitate the rapid and convenient measurement of the LSER model presentsed in reference [5]. An application of the LSER model on ionic liquids will be presented. The last paragraph will be focused on the use of the LSER model coupled to a group contribution model for the estimation of the partition coefficient of organic compounds in ionic liquids.


2. Solvation models used in gas chromatography

2.1. Linear solvation energy relationship

In numerous solvation models, the partition of a solute between the gas phase and a solvent (or stationary phase) may be described by a cavity theory of solution [6-7] in which the solvation process is divided into three steps. In the step 1, a cavity of suitable size to accommodate the solute is created in the solvent. This process is endoergic because work is required to disrupt solvent–solvent interactions. In the second step, the solvent molecules around the cavity are reorganized from their original positions to the positions they will adopt when the solute is in equilibrium with the solvent. Of course, these positions are not fixed but are averages of solvent positions. The Gibbs energy change for such reorganization is assumed to be negligible, by analogy with the melting of a solid. However, enthalpy and entropy changes in reorganizations may be large, again by analogy with the melting of a solid. In the last step, the solute is introduced into the reorganized cavity, and various solute–solvent interactions are set up, all of which are exoergic.

Several thermodynamic properties are related to partitioning between water and other phases, for example octanol (Kow) or the pure compound itself (water solubility). These partitioning processes can be understood from thermodynamic concepts—like free energy, chemical potential and fugacity. The equilibrium partition constant between two phases, on a mole fraction basis, can be expressed as:

Ki 12=exp(Δ12Gi/RT)E1

where Δ12Gi is the Gibbs energy (or free energy) of transfer between the two phases, R is the general gas constant and T is the absolute temperature. The Gibbs energy sums up both the enthalpic (Δ12H) and entropic (Δ12S) effects resulting from changes in intermolecular interactions:


The calculation of partitioning from structural or other descriptors therefore requires the modelling of these effects. The Gibbs energy change can also be separated into the contributions of van der Waals and polar interactions, assuming that these are additive:


In the early 1980’s, Taft & Kamlet [8-12] have developed the basic concept of linear solvation energy relationships (LSERs). They have demonstrated for thousands of chemical systems that some property which is linearly related to either a free energy of reaction, a free energy of transfer, or an activation energy can be correlated with various fundamental molecular properties of the solvents or solutes involved. Chromatographic retention and in particular logarithmic retention factors (log k′), logarithmic partition coefficients (logKL) are linear free energy parameters and as such one can linearly correlate these data with the molecular properties of the solutes using the LSER model [13-16].

Figure 1.

Model of the solvation process.

The most recent representation of the LSER model proposed by Abraham and coworkers [5, 17-19] is given by equation (4)


Where SP is a solvation parameter related with the free energy change such as gas-liquid partition coefficient, specific retention volume or adjusted retention time at a given temperature. The capital letters represent the solutes properties and the lower case letters the complementary properties of the stationary phase. The solute descriptors are the excess molar refraction R2, dipolarity/ polarizabilityπ2H, hydrogen bond acidity basicity, α2Handβ2H, respectively, and the gas-liquid partition coefficient on n-hexadecane at 298 K, LogL16. The solute descriptors may be determined using inverse gas chromatography or estimated using a group contribution method. A databank of descriptors for about 3000 compounds may be found in the literature [2, 20, 21]. The coefficients c, e, s, a, b and l are not simply fitting coefficients, but they reflect complementary properties of the solvent phase. These coefficients are determined by mutiple linear regression of equation (4). This model was strongly applied to characterize chemicals products, petroleum fluids.

2.2. Determination of LSER parameters of pure solutes

2.2.1. Determination of Log L16

To preserve the general character of equation (4), all characteristic parameters should be carefully determined and correlation between parameters should be avoided. Experimental procedures of successive determination of LSER parameters were described in the literature [20, 22-24]. Log L16 characterizes the most general interactions present in every physical system and should be determined before other parameters [24]. The original values of log L16 were determined from retention data of organic compounds on n-hexadecane coated packed columns at 298.2 K [5]. A number of papers proposed techniques based on the use of either capillary or packed columns in wide temperature ranges and replacing the n-hexadecane with other non-polar stationary phases. The gas liquid partition coefficient of a solute is directly related to the experimental capacity factor, k, by equation:


where VM and VS are the volumes of the mobile and stationary phases, respectively, and Φ is the phase ratio (VS/ VM). Experimental determination of log L16 is often very difficult. Adsorption phenomena introduce an important error in determination of the capacity factor. Zhang et al. [23] determined log L16 with capillary columns coated with n-hexadecane and concluded that results were not influenced by adsorption in this case. However, recent results presented in reference [25] showed that this improvement is not general and that an interfacial adsorption still exists with capillary columns. Moreover, it is very difficult to obtain absolute retention data using this technique. Li et al. [24] and Abraham et al. [25] studied influence of the solute support and of the stationary phase loading on adsorption phenomena. They concluded that the high loading ratio (up to 20 %) of the stationary phase and the high temperature of the column allow to reduce adsorption. In this case, the knowledge of R2, α2Hand π2H parameters is necessary. The number and the nature of parameters needed for calculation depend on the stationary phase used.

Serious difficulties arise when log L16 of non-volatile compounds is to be measured. This is due to the definition of log L16 itself. Indeed, direct experimental determination of log L16 of compounds less volatile than n-hexadecane is impossible. In the case of heavy compounds slightly more volatile than n-hexadecane the experiment is possible but difficult, especially at 298.2 K. Often it is recommended to measure retention times at higher temperatures and then extrapolate partition coefficients to the ambient temperature. In this case, the quality of results depends strongly on the extrapolation method used. The problem of the temperature dependence of retention times was often discussed in literature. A suitable extrapolation procedure was described in reference [26]. In the case of compounds less volatile than n-hexadecane several authors proposed to work with columns coated with long chain branched paraffins and to establish relationships between corresponding partition coefficients and log L16 [16, 26]. Defayes et al. [26] worked with apolane coated stationary phase (apolane is a C87H176 branched alkane). This column can be used at temperatures up to 550 K without a weight loss of the stationary phase. Moreover, it was shown that the effect of adsorption at the liquid-gas interface is negligible in this case [26]. However, this opinion is not generally accepted and Werckwerth et al. [16] found the influence of adsorption in the case of strongly polar compounds. The same authors observed that a linear relationship exists between gas-apolane partition coefficients, log L87 and log L16 and that the data obtained with apolane can be used to estimate the value of log L16. Moreover, they demonstrated that strong correlation between both partition coefficients exists also for log L87 determined at significantly higher temperature [16]. Recently, several authors investigated the use of predictive methods to estimate log L16 [27, 28]. This approach is particularly interesting to determine the log L16 of nonvolatile compounds.

2.2.2. Determination of LogL16 using capillary columns

A direct determination of the stationary phase mass is difficult in the case of capillary or Megabore columns and the use of relative methods to determine partition coefficient from retention data is often preferred. In this case, different approaches were proposed to calculate partition coefficient based on the well-established value of log Ln-hexane [23]. Corresponding equation is as follows:


where tR and tm are the solute retention time and the dead time of the column respectively. Retention data of the solute X and of n-hexane should be determined at the same temperature. Available data of log Ln-hexane on apolane were determined at 312.4 K [16]. On the other hand, it was shown that partition coefficients determined at two temperatures are linearly correlated [16].


New relationship can be established between the partition coefficient of the solute X at temperature T and the partition coefficient of the n-hexane at temperature T’. This relationship is based on the observation verified with several n-paraffin stationary phases that equation (7) for one stationary phase is reduced to the following form:


Equations (6) and (8) lead to relationship between the partition coefficient at temperature T of the solute X and the partition coefficient of n-hexane determined at temperature T’:


Apolane coated capillary columns are considered as an appropriate tool to determine log L16 of heavy compounds. Studies of Defayes et al. [26] and Werckwerth et al. [16] provided arguments supporting this opinion. As it can be seen in Table 1, results obtained with a similar but deactivated column are in good agreement with literature. Chromatographic peaks obtained with non-polar and polar compounds were symmetric. While this method gave good results at high temperatures, the column was deactivated irreversibly within few hours. A probable explanation of this phenomenon is that the adhesion between apolane and the deactivated silica does not assure the film stability at higher temperature [24, 27]. Our experience indicates that the use of commercially available apolane coated capillary columns should be limited to low temperature ranges. In the case of heavy compounds this implies very long retention times and imposes injection of samples of important volume, which induces adsorption effects. Consequently, to determine log L16 of heavy compounds we decided to use packed columns with long chain n-alkane stationary phase. Moreover, results obtained with a non-deactivated column indicate that retention times are influenced by adsorption phenomena. Indeed, polar solutes exhibited severely asymmetric peaks and their retention times strongly depended on the sample size. Retention times of alcohols are longer than the literature values that may indicate the presence of active sites.

Table 1.

Partition coefficient log L16 exp. at 298.2K obtained with a C50H102 packed column at 373.2K. log L16 a : data taken from [5,14,21,31].

2.2.3. Determination of LogL16 using packed columns

Problems of the capillary column stability encouraged us to review the possibility of application of packed columns for determining log L16 of non-volatile compounds. Stationary phases used were long chain n-alkanes, n-hexatriacontane and n-pentacontane. They were used at temperatures up to 320 K without significant loss of weight. The essential problem encountered with packed columns concerned adsorption effects [30, 23]. Mutelet & Rogalski [2] used teflon columns, inert and stable up to 330 K. Selecting an appropriate support material can reduce the adsorption on the surface of the support. Preliminary tests showed that the best results were obtained with the Chromosorb PAW DMCS and the Chromosorb WHP. Both supports were loaded with 25% of n-pentacontane. The fact that with Chromosorb PAW DMCS retention times depend on the sample size and chromatographic peaks are asymmetric indicates the presence of adsorption. Moreover, retention times of alcohols are longer than expected, indicating the presence of active adsorption sites on the support surface. The Chromosorb WHP support has a lower specific area and a smaller concentration of hydroxyl groups which reduces the adsorption. Results obtained are in good agreements with literature data for most of the compounds studied and retention times depend only slightly on the sample size. However, retention times observed with polyaromatic hydrocarbons are still longer than expected. Good results were obtained by deactivating the column (Table 1) with Silyl 8, as recommended by [31]. The use of packed columns with Chromosorb WHP coated with n-alkane and deactivated with Silyl 8 made it possible to obtain a homogenous set of log L16 in good agreement with literature data.

2.2.4. Determination of LogL16 using temperature gradient method

The packed column technique can be used to measure log L16 data of volatile organic compounds. The reasonable limit of application of this method is the retention time of n-eicosane. Experimental difficulties make hazardous quantitative determination log L16 of heavier compounds. To enlarge the applicability of chromatographic methods to organic compounds less volatile than n-eicosane, a method based on the temperature gradient chromatography can be used. Recently, Donovan [32] showed that retention times of heavy organic compounds obtained in a gradient mode are linearly related to the logarithm of the vapor pressure at 298.2 K. The authors used a DB-1 megabore column at high flow rates of the gas phase. This method making it possible to reduce considerably retention time was applied to determine vapor pressures of pesticides and polyaromatic hydrocarbons. Nevertheless, the stationary phase DB-1 is slightly polar [24]. Corresponding system parameters of the poly(dimethylsiloxane) immobilized in DB-1 column were published in the reference [24]. Values determined at t= 60°C are as follows: r =0, s = 0.211, a = 0.308 and b = 0. Therefore, experimental results obtained with a DB-1 column can be used to determine log L16 only if LSER parameters expressing solute polarity are known. No general relationship between the reduced retention time and log L16 valid with all organic compounds can be obtained without the knowledge of above parameters.

However, this approach can be used to establish relationship between the reduced retention time and log L16 within a series of compounds. Indeed, polar parameters vary only slightly and in a regular way within a series. Moreover, certain parameters decrease strongly with rising temperature [24]. Therefore, it can be supposed that the effect of the stationary phase polarity is nearly constant within a homologous series of moderately polar compounds. Measurements performed in a gradient mode with several homologous series confirmed this hypothesis. However, linear relationship does not afford the precision required for the log L16 determination. It was noticed that not only the reduced retention time tR but also the corresponding temperature T is needed to establish the appropriate function. Function log L16 = f(tR,T) is linear with R=0.996 that is not enough to represent the log L16 with the precision required. We found that the suitable function is as follows:


In the case of n-alkanes, function f(tR) was obtained with log L16 literature data of n-alkanes from n-dodecane up to n-docosane. The plot of f(tR) and values of parameters determining this function are given in Figure 1. The log L16 of n-alkanes up to n=38 calculated with equation (10) using gradient mode results are presented in Table 2. It is reasonable to suppose that partition coefficients of heavy n-alkanes up to approximately C45H92 can be obtained with gradient method. It should be pointed out that the present approach based on the gradient mode chromatography can be used only to determine log L16 within a homologous series of moderately polar compounds. The use of the present method with less polar stationary phase (recently, Li et al. [24] shown that in the case of SPB columns r = 0, a = 0, b = 0) can facilitate the study of polar compounds and perhaps obtain more general results.

Table 2.

Retention times tR and partition coefficients log L16 exp of n-alkanes determined at 298.2°C using the gradient method on DB-1 column. Log L16 a,b a: -literature data [5,14,21,31]. b: -estimates calculated using the group contribution method proposed by reference [2].

Figure 2.

Plot of f(tR) function established with temperature gradient method; aromatics and polyaromatics (ο), n-alkanes ().

2.2.5. Determination of polar LSER parameters

The excess molar refraction E is defined by the difference between the value for the solute molar refraction and the molar refraction for an alkane of the same characteristic volume:

R2= MRX(solute)  MRX(alkane of the same VX)E11

The solute molar refraction is calculated from the following equation:


Where VX is the specific volume in cm3.mol-1/100 and n the refractive index of the solute.

Abraham et al. [18] set out to construct scales of solute hydrogen bond acidity and hydrogen bond basicity using logK values for reaction 1 in tetrachloromethane. The authors set out logK values for a series of acids against 45 given bases. If log K values for acids in a given reference base is plotted against log K values for acids in another reference base, a series of straight lines is observed with an intersection at a magic point of -1.1 log units. It is found that:


Where LB and DB characterize the base and logKAH values characterize the series of acids. The α2H and β2Hparameters are then defined by:

The α2H and β2Hparameters can be determined at the same time with π2Hparameter. Retention data (such as retention volume, partition coefficients) of the solute in three stationary phases at least of different polarity can be used to determine theα2H, β2Handπ2H. Each stationary phase can be described by equation (16):

At first, the LSER parameters c, s, a, b and l of each stationary phase are determined by multiple linear regression using solutes for which R2, α2H, β2H, π2Hand logL16 are known. Then, α2H, β2H, π2Hcan be determined by multiple linear regression.

2.2.6. Group contribution Method for calculation of LSER parameters of organic compounds

Predictive methods allow to calculate these physico-chemical parameters which are inacessible via direct experiment. This alternative is particularly interesting in the case of log L16 of nonvolatile compounds. We consider that experimental methods described in the preceeding paragraph are useful for determination of LSER parameters of volatile and moderatly non-volatile compounds. Therefore, the large data bank of log L16 values already available in the literature can be used to establish group additivity rules and to predict log L16 of less volatile compounds. Havelec & Sevcik [27,28] presented a general group contribution method making it possible to calculate accurate estimates of log L16 of about 2000 organic compounds. The number of groups necessary to obtain good estimates of log L16 depends on the complexity of the molecular structure and rises in the case of polyfunctional molecules. This explains a high number of adjustable parameters used in the model [27,28]. For instance, log L16 of non-aromatic hydrocarbons is established with 33 parameters and 9 structural contributions. The total number of all group parameters, interactional parameters and structural contributions is of 131. The contribution of a given group is represented in the reference [27,28] with three parameters related to the structure of the molecule and to its interactions with the stationary phase. As log L16 is dependent on the solute vapor pressure and on the infinite dilution activity coefficient this approach is basically correct. However, molecular interactions are always related to n-hexadecane and certain parameters can be correlated. Platts et al. [20] recently proposed a new predictive method based on a careful analysis of contributions of various functional groups to establish log L16 and other LSER parameters. Therefore, molecular segments were defined in view to obtain good estimates of each. The log L16 of hydrocarbons is calculated with 9 parameters only. This method was established with 81 parameters, using a databank of 1908 compounds. A new model was proposed to calculating log L16 for nonvolatile organic compounds with special attention paid to heavy hydrocarbons. Data for 550 organic compounds containing mainly hydrocarbons and members of homologous series were used in regression. Basic heteroatom segments were taken into account but the polyfunctional organic compounds were not dealt with. Values of log L16 were taken from literature [5, 21]. To elaborate the group contribution method a simple and efficient approach was used. Accordingly, log L16 of the compounds X was calculated with the following expression:


where ci is the contribution of the group “i” and ni is the number of groups “i” in the compound X.

Platts et al. [20] have developed and tested additive models for six important molecular LFER descriptors, namely, R2, α2H, β2H, π2Hand logL16. Five of these six, all bar α2H are calculated from a single set of 81 atom and group fragments, while β2H is calculated from a separate set of 51 fragments. In general, the linear fit obtained with these additive models is good, with R2 and log L16 in particular giving excellent correlation. Splitting the data into training and test sets has also tested the predictive ability of such models, and is found to be almost as accurate as the full regressions. The performance of the method in calculating descriptors for “difficult” structures, ones containing intramolecular interactions such as hydrogen bonds, has been analyzed. Variations in descriptors due to such interactions are generally found to be reproduced, though inevitably some small discrepancies are found. In conclusion, this model is particurlarly powerful and useful for the prediction of LSER parameters of heavy and complicated molecules.


3. Application of linear solvation energy relationship on ionic liquids

The LSER model may be used to characterize the stationary phases in chromatography. In this case, a large number of solutes (between 20 and 50) for which LSER parameters R2, α2H, β2H, π2Hand logL16 are known have to be injected. The LSER parameters characterizing the stationary phases (c, s, r, a, b and l) are determined by multiple linear regression. In the literature, there is a large amount of data of partition coefficients or activity coefficients measured by gas–liquid chromatography or by dilutor technique. Some system constants for various ionic liquids and classical solvents at 25 °C are summarized in Table 3. The data for the 1-ethanol-3-methylimidazolium tetrafluoroborate, 1-ethanol-3-methylimidazolium hexafluorophosphate, 1,3-dimethylimidazolium dimethylphosphate and 1-ethyl-3-methylimidazolium diethylphosphate [33], 1-Butyl-3-methylimidazolium tetrafluoroborate [34], n-Acryloyloxypropyl-N-methylimidazolium bromide and n-Methacryloyloxyhexyl-N-methylimidazolium bromide [35], 1-Propenyl-3-alkyl-imidazolium bromide [36], 1-butyl-3-methylimidazolium octyl sulfate and 1-ethyl-3-methylimidazolium tosylate [37], Triethylsulphonium bis(trifluoromethylsulfonyl)imide [38], 1-Methyl-3-ethylimidazolium bis(trifluorosulfonyl)-amide and 1.2-Dimethyl-3-ethylimidazolium bis(trifluorosulfonyl)-amide [39] were taken from the sources indicated. Poole & Poole [40] found that the system constants of LSER model for the room temperature ionic liquids fall into the range e = −0.62 to 0.86, s = 1.7–2.8, a = 2.1–7.3, b = 0–1.07, and l = 0.35–0.96. Compared with the scale of the polar organic solvents e = −0.60 to 0.82, s = 0.54–2.8, a = 0.28–5.50, b = 0–4.8, and l = −0.21 to 0.98, we can see that both scales are similar indicating that the solvation properties for the room temperature ionic liquids are classical and fit quite well into the scales developed for polar molecular solvents.

The (c + lL) term gives information on the effect of cohesion of the ionic liquids on solute transfer from the gas phase. In general, the ionic liquids are cohesive solvents; they interact weakly via nonbonding and π-electrons (r system constant is zero) and are not much different to other polar non-ionic liquids. The ionic liquids are roughly as dipolar/polarizable as classical solvents. The hydrogen-bond basicity of the ionic liquid (a system constants) are considerably larger than values obtained for non phases (0-2.1) [1]. The hydrogen-bond basicity of Ils depends on the anion grafted on the cation. ionic liquids can be slightly more hydrogen-bond basic than dimethyl sulfoxide and N-methylpyrrolidinone, and are weak to moderate hydrogen-bond acids, similar to the aliphatic alcohols. From Table 3 and data collected in the reference [40], we can see that the hydrogen-bond acidity of the ionic liquids depends largely on the cation and is lower for the 1,3-dialkylimidazolium salts with an alkyl group at C-2 position than 1,3-dialkylimidazolium salts.


4. Predictive models based on LSER model coupled to a group contribution method

Solvation model may be also used to set up correlation between thermodynamic properties and LSER parameters. Abraham et al. [41,42] reported mathematical correlations based on the general Abraham solvation parameter model for the gas-to-solvent, KL, and water-to-solvent, P, partition coefficients. Recently, [43-46] modified the Abraham solvation parameter model:

LogKL= ccation+ canion+ (ecation+ eanion) ·E+ (scation+ sanion) ·S+                 + (acation+ aanion) ·A+ (bcation+ banion) ·B+ (lcation+ lanion) ·LE18
LogP= ccation+ canion+ (ecation+ eanion) ·E+ (scation+ sanion) ·S+ (acation+ aanion) ·A + (bcation+ banion) ·B+ (vcation+ vanion) ·VE19

by rewriting each of the six solvent equation coefficients as a summation of their respective cation and anion contribution. The dependent variables in equations (18) and (19) are solutes descriptors as follows: E and S refer to the excess molar refraction in units of (cm3.mol-1)/10 and dipolarity/polarizability descriptors of the solute, respectively, A and B are measures of the solute hydrogen-bond acidity and basicity, V is the McGowan volume in units of (cm3.mol-1)/100 and L is the logarithm of the gas-to-hexadecane partition coefficient at 298 K. Sprunger et al. calculated equation coefficients for 8 cations and 4 anions using a database that contained 584 experimental log KL and 571 experimental log P values. No loss in predictive accuracy was observed by separating the equation coefficients into individual cation-specific and anion-specific values. The major advantage of splitting the equation coefficients into individual cation-specific and anion-specific contributions is that one can make predictions for more ILs. In Sprunger’s approach, the major advantage of splitting the equation coefficients into individual cation-specific and anion-specific contributions is that one can make predictions for more Ils. Most of the cations are alkylimidazolium. The use of this model is somewhat limited since it can not be extrapolated to alkylimidazolium ionic liquids not initially defined by the method (e.g. with long alkyl chains).

Table 3.

LSER descriptors of ionic liquids determined at 313.15 K.

In the development of Mutelet et al. [47], the cation with its alkyl chains is splitted in different contributions: (CH3, CH2, N, CHcyclic…). The approach allows to have a predictive model. The aim of this work was to develop a group contribution method allowing to estimate the log KL and Log P of organic compounds in ionic liquids at 298 K. Using the LSER model proposed by Abraham, the group contribution method expresses LSER coefficients ci, ei, si, ai, bi and li of equation (19). orci', ei', si', ai', bi'and vi of equation (20) by:

LogKL=i21ni×ci+i21ni×ei·E+i21ni×si·S+i21ni×ai·A+i21ni×bi·B+i21ni×li ·LLogP=i21ni×ci'+i21ni×ei'·E+i21ni×si'·S+i21ni×ai'·A+i21ni×bi'·B+i21ni×vi·VE20

Where ni is the number of group i present in the ionic liquid.

Mutelet et al. [47] proposed to extend the temperature dependent GC-LSER in view of determining the partition coefficient of organic compounds in ionic liquids. The GC-LSER can be rewriting as followed:


The experimental data used to calculate Abraham’s model ion-specific equation coefficients were taken from the collection of [43-46] and were updated with recent data. A total of 1450 gas-liquid partition coefficients and 1410 water-to-liquid partition coefficients were used for the calculation. Solutes were mainly n-alkanes, cycloalkanes, alkenes, alkynes, aromatics, alcohols, ethers, aldehydes, ketones, chloroalkanes. The E-scale varies from 0 to 1.5, the S-scale from 0 to 1.72, the A-scale from 0 to 1.04, the B-scale from 0 to 1.28, the L-scale from –1.200 to 7.833 and the V-scale from 0.109 to 1.799. The dataset is composed of 27 imidazolium based ionic liquids, 3 ammonium, 3 pyridinium and 4 pyrolidinium based ionic liquids. The authors also add sulphonium and phosphonium ionic liquids although only one set of KL (or P) data may be found for these families. The twenty one groups which are defined in this method are listed in Table 4. The decomposition into groups of the ionic liquids is very easy, that is as simple as possible. No substitution effects are considered. No exceptions are defined. In Figure 3 are represented all ionic liquids studied in this work. Five groups are defined to describe the chains R1, R2, R3 and R4 grafted on the cation: CH3, CH2, -O-, -O-Ncycl and –OH. These groups allow the calculation of partition coefficients of alkyl based ionic liquids but also functionalized ionic liquids such as ether, alcohols. The remaining seven groups are: CH2cyclic, CHcyclic, Ccyclic, Ncyclic, N+ (ammonium cation), P+ (phosphonium cation) and S+ (sulphonium cation).

More precisely, Ncyclic represents two structures: N+=andN. Nine groups are used for anions: bis(trifluoromethylsulfonyl)imide :(TF)2N, hexafluorophosphate:PF6, tetrafluoroborate:BF4, ethylsulfate:EtSO4, octylsulfate:OcSO4, thiocyanate:SCN, trifluoromethylsulfonate :CF3SO3, trifluoroacetate : ACF3and dicyanamide:(CN)2N. As an example, let’s have a look at the decomposition of 1-butyl-3-methylimidazolium hexafluorophosphate. In this case, the decomposition of the molecule into elementary groups is: 2 group 1 (-CH3) + 3 group 2 (-CH2) + 3 group 7 ( Ccyclic) + 2 group 9 (Ncyclic) + 1 group 14 (PF6-).

Figure 3.

Cation of six families of ionic liquids.

Group contribution model coupled to LSER (GC-LSER) for estimating the gas-to-ionic liquids partition coefficients and water-to-ionic liquids partition coefficients allows to predict with good accuracy Log KL and Log P at 298 K of not only alkyl based ionic liquids but also functionalized ionic liquids. The parameters of the group contribution methods were determined for imidazolium, pyridinium, pyrrolidinium, phosphonium, ammonium and sulphonium based ionic liquids containing several different anions. A comparison between the experimental and calculated values showed that the proposed models describe the experimental data available with a mean absolute error of about 0.15 log unit. While the model is probably somewhat limited in prediction for pyridinium and pyrrolidinium based ionic liquids because of the poor dataset for these cations, results obtained are satisfactory.

Table 4.

Description of the 21 groups used for the estimation of LogKL and LogP.


5. Conclusion

The solvation parameter model is suitable for describing the retention properties of molecules in chromatographic systems. To establish the system properties requires identification of a group of compounds with well known descriptor values. We have shown that all LSER parameters of solutes may be determined using gas chromatography or experimental techniques. The solvation model may be used either for the physico-chemical characterization of the stationary phases or for the establishment of a suitable quantitative structure–property relationship to facilitate the prediction of further system properties for compounds lacking experimental values.


  1. 1. PooleC. F.Chromatographic and spectroscopic methods for the determination of solvent properties of room temperature ionic liquidsJournal of ChromatographyA 2004103749
  2. 2. MuteletF.RogalskiM.Experimental determination and prediction of the gas-liquid n-hexadecane partition coefficients J. Chromatogr. A 2001923153
  3. 3. RussomC. L.BretonR. L.WalkerJ. D.BraduryS. P.An overview of the use of quantitative structure-activity relationships for ranking and prioritizing large chemical inventories for environmental risk assessments Environ. Toxicol. Chem. 2003221810
  4. 4. KatritzkyA. R.MaranU.LobonovV. S.KarelsonM.Structurally Diverse Quantitative Structure-Property Relationship Correlations of Technologically Relevant Physical Properties J. Chem. Inf. Comput. Sci. 2000401
  5. 5. AbrahamM. H.GrellierP. L.Mc GillR. A.Determination of Olive Oil-Gas and Hexadecane-Gas Partition Coefficients, and calculation of the corresponding Olive Oil-Water and Hexadecane-Water Partition Coefficients. J. Chem. Soc. PERKIN TRANS II 1987797803
  6. 6. VithaM.CarrP. W.The chemical interpretation and practice of linear solvation energy relationships in chromatographyJournal of Chromatography A20061126143
  7. 7. AbrahamM. H.PooleC. F.PooleS. K.Classification of stationary phases and other materials by gas chromatography Journal of Chromatography A, 199984279
  8. 8. KamletM. J.AbboudJ. L.TaftR. W.1977The solvatochromic comparison method. 6. The π scale of solvent polarities. Journal of the American Chemical Society 1977; 189960276038
  9. 9. KamletM. J.AbboudJ.L. M.AbrahamM. H.TaftR. W.Linear solvation energy relationships. 23. A comprehensive collection of the solvatochromic parameters, π*, α, and β, and some methods for simplifying the generalized solvatochromic equation. Journal of Organic Chemistry 1983174828772887
  10. 10. KamletM. J.TaftR. W.Linear Solvation Energy Relationships. Local Empirical Rules-- or Fundamental Laws of Chemistry? A Reply to the Chemometricians., Acta Chem. Scand., B 198539611
  11. 11. KamletM. J.DohertyR. M.AbrahamM. H.MarcusY.TaftR. W.Linear solvation energy relationships. 46. An improved equation for correlation and prediction ofoctanol/water partition coefficients of organic nonelectrolytes (Including strong hydrogen bond donor solutes). Journal of Physical Chemistry, (1998189252445255
  12. 12. TaftR. W.KamletM. J.The solvatochromic comparison method. 2. The α-scale of solvent hydrogen-bond donor (HBD) acidities. Journal of the American Chemical Society 1976109828862894
  13. 13. AbrahamM. H.WhitingG. S.Hydrogen-bondingPart. 2.Characterization of soybean oil and prediction of activity coefficients in soybean oil from inverse gas chromatographic data. J. Am. Oil Chem. Soc., 1992691236
  14. 14. AbrahamM. H.Andonian-HanftvanJ.Osei-OwusuJ. P.SakellariouP.UrietaJ. S.LópezM. C.FuchsR.Hydrogenbonding.Part2.The solvation properties of methylene iodide. J. Chem. Soc., Perkin Trans. II 19933299
  15. 15. ParkJ. H.CarrP. W.Interpretation of normal-phase solvent strength scales based on linear solvation energy relationships using the solvatochromic parameters π*, α and β. J. Chromatogr. 19894653123136
  16. 16. WeckwerthJ. D.CarrP. W.VithaM. F.NasehzadehA. A.comparisonof.gas-hexadecanegas-apolanepartition.coefficientsAnal.Chem. 1998703712
  17. 17. AbrahamM. H.WhitingG. S.DohertyR. M.HydrogenBonding.Part1.NewA.Methodfor.theCharacterization.ofG. L. C.Stationary-ThePhases.LaffordData.SetJ. Chem. Soc. PERKIN TRANS II 199014511460
  18. 18. AbrahamM. H.WhitingG. S.DohertyR. M.ShuelyW. J.Hydrogenbonding. X. V. I. A.newsolute.solvationparameter. π. H.fromgas.chromatographicdata.J. Hydrogen bonding XVI. A new solute solvation parameter, π2H, from gas chromatographic data. J. Chromatogr. 1991587213
  19. 19. AbrahamM. H.Scales of Solute Hydrogen-bonding: Their construction and Application to Physicochemical and Biochemical ProcessesChem. Soc. Rev. 19932273
  20. 20. PlattsJ. A.ButinaD.AbrahamM. H.HerseyA. Estimation of molecular linear free energy relation descriptors using a group contribution approach. Journal of Chemical Information and Computer Sciences 1999; 39835 .
  21. 21. AbrahamM. H.PlattsJ. A.Hydrogen bond structural group constants.Journal of Organic Chemistry 2001663484
  22. 22. WeckwerthJ. D.CarrP. W.Study of Interactions in Supercritical Fluids and Supercritical Fluid Chromatography by Solvatochromic Linear Solvation Energy Relationships.Anal. Chem. 199870714041410
  23. 23. ZhangY.DallasA. J.CarrP. W.Critical comparison of gas hexadecane partitioncoefficients as measured with packed and open-tubular capillary columns J. Chromatogr. 19936384356
  24. 24. LiQ.PooleC. F.KiridenaW.KoziolW. W.Chromatographic methods for the determination of the logL16 solute descriptorAnalyst, 20001252180
  25. 25. AbrahamM. H.Andonian-HaftenJ.MyDu. C.Osei-OwusuJ. P.SakellariouP.ShuelyW. J.PooleC. F.PooleS. K.Comparison of uncorrected retention data on a capillary and a packed hexadecane column with corrected retention data on a packed squalane column J. Chromatogr. A 1994688125
  26. 26. DefayesG.FritzD. F.GörnerT.HuberG.De ReyffC.KovatsE.szJ.ChromatogrOrganicsolutes.inparaffin.solventsInfluence.ofthe.sizeof.thesolvent.moleculeon.solutiondata.1990500139
  27. 27. HavelecP.SevcikJ. G. K.Concept of additivity for a nonpolar solute-solvent criterion Log L(16)- Nonaromatic compounds J. Chromatogr. A 1994677319
  28. 28. HavelecP.SevcikJ. G. K.1996Extended additivity model of parameter log(L(16)) ,J. Phys. Chem. Ref. Data 1996; 251483
  29. 29. TranchantJ.Manuel pratique de chromatographie en phase gazeuseMasson, Paris, 1995
  30. 30. DallasA. J.CarrP. W.Direct chromatographic comparison of the relative adsorption activity of various types of capillary transfer tubing Anal. Chim. Acta 199125183
  31. 31. KerstenB. R.PooleC. F.Influence of concurrent retention mechanisms on the determination of stationary phase selectivity in gas chromatography J. Chromatogr. 19873991
  32. 32. DonovanS. F.1996New method for estimating vapor pressure by the use of gas chromatography J. Chromatogr. A 1996; 749123
  33. 33. RevelliA.L.MuteletF.JaubertJ.N.Partition coefficients of organic compounds in new imidazolium based ionic liquids using inverse gas chromatographyJournal of ChromatographyA 200912162347754786
  34. 34. RevelliA.L.MutetetF.TurmineM.SolimandoR.JaubertJ.N.Activity coefficients at infinite dilution of organic compounds in 1 -butyl-3-methylimidazolium tetrafluoroborate using inverse gas chromatographyJournal of Chemical and Engineering Data2009
  35. 35. MuteletF.JaubertJ.N.RogalskiM.HarmandJ.SindtM.MieloszynskiJ.L.Activity coefficients at infinite dilution of organic compounds in 1-(meth)acryloyloxyalkyl-3-methylimidazolium bromide using inverse gas chromatography. Journal of Physical Chemistry B 20081121237733785
  36. 36. MuteletF.JaubertJ.N.RogalskiM.BoukherissaM.DickoA.Thermodynamic properties of mixtures containing ionic liquids: Activity coefficients at infinite dilution of organic compounds in 1-propyl boronic acid-3-alkylimidazolium bromide and 1-propenyl-3-alkylimidazolium bromide using inverse gas chromatographyJournal of Chemical and Engineering Data200651412741279
  37. 37. MuteletF.JaubertJ.N.Accurate measurements of thermodynamic properties of solutes in ionic liquids using inverse gas chromatographyJournal of Chromatography A2006256 EOF267 EOF
  38. 38. DomańskaU.MarciniakA.Activity coefficients at infinite dilution measurements for organic solutes and water in the ionic liquid triethylsulphonium bis(trifluoromethylsulfonyl)imide. Journal of Chemical Thermodynamics 2009416754758
  39. 39. KrummenM.WasserscheidP.GmehlingJ.Measurement of activity coefficients at infinite dilution in ionic liquids using the dilutor techniqueJ. Chem. Eng. Data 200247614111417
  40. 40. PooleC. F.PooleS. K.Extraction of organic compounds with room temperature ionic liquidsJournal of Chromatography A201012171622682286
  41. 41. AbrahamM. H.ZissimosA. M.HuddlestonJ. G.WillauerH. D.RogersR. D.AcreeW. E.Jr Some novel liquid partitioning systems: Water-ionic liquids and aqueous biphasic systemsInd. Eng. Chem. Res. 200342413
  42. 42. AcreeW. E.Jr AbrahamM. H.The analysis of solvation in ionic liquids and organic solvents using the Abraham model linear free energy relationship. J. Chem. Technol. Biotechnol. 2006811441
  43. 43. SprungerL.ClarkM.AcreeW. E.Jr AbrahamM. H.Characterization of room-temperature ionic liquids by the Abraham model with cation-specific and anion-specific equation coefficients.J. Chem. Inf. Model. 20074711231129
  44. 44. Sprunger, L.M.; Proctor, A.; Acree, W.E., Jr. & Abraham, M.H. LFER correlations for room temperature ionic liquids: Separation of equation coefficients into individual cation-specific and anion-specific contributions.Fluid Phase Equilib. 2008265104
  45. 45. SprungerL. M.AchiS. S.Acree JrW. E.AbrahamM. H.LeoA. J.HoekmanD.Correlation and prediction of solute transfer to chloroalkanes from both water and the gas phaseFluid Phase Equilib. 2009281144
  46. 46. SprungerL. M.GibbsJ.ProctorA.Acree JrW. E.AbrahamM. H.MengY.YaoC.AndersonJ. L.Linear free energy relationship correlations for room temperature ionic liquids: revised cation-specific and anion-specific equation coefficients for predictive applications covering a much larger area of chemical spaceInd. Eng. Chem. Res. 2009484145
  47. 47. MuteletF.Ortega-VillaV.MoiseJ.AcreeC.Jr W. E.BakerG. A.Prediction of Partition Coefficients of Organic Compounds in Ionic Liquids Using a Temperature-Dependent Linear Solvation Energy Relationship with Parameters Calculated through a Group Contribution Method Journal of Chemical and Engineering Data 2011563598

Written By

Fabrice Mutelet

Submitted: 24 November 2011 Published: 24 October 2012