Solvation parameter model

In a series of papers published throughout the 1980s, Colin Poole and his co-workers investigated the solvation properties of a wide range of alkylammonium and, to a lesser extent, phosphonium salts. Parameters such as McReynolds phase constants were calculated by using the ionic liquids as stationary phases for gas chromatography and analysis of the retention of a variety of probe compounds. However, these analyses were found to be unsatisfactory and were abandoned in favour of an analysis that used Abraham s solvation parameter model [5]. 

In initial work, a total of 17 different ILs were evaluated by the solvation parameter model [8]. Ten of these ILs were comprised of imidazolium or pyrolidinium cations paired with different anions. Many of these compounds represent the traditional class of IL solvents that have been used extensively in organic synthesis reactions or in other analytical uses. The remaining seven ILs consisted of substituted ammonium cations that have proven to be successful analyte matrices in matrix-assisted laser desorption ionization (MALDI) mass spectrometry [11]. 

Retention of Rohrschneider-McReynolds standards of selected chiral alcohols and ketones was measured to determine the thermodynamic selectivity parameters of stationary phases containing (- -)-61 (M = Pr, Eu, Dy, Er, Yb, n = 3, R = Mef) dissolved in poly(dimethylsiloxane) . Separation of selected racemic alcohols and ketones was achieved and the determined values of thermodynamic enantioselectivity were correlated with the molecular structure of the solutes studied. The decrease of the ionic radius of lanthanides induces greater increase of complexation efficiency for the alcohols than for the ketone coordination complexes. The selectivity of the studied stationary phases follows a common trend which is rationalized in terms of opposing electronic and steric effects of the Lewis acid-base interactions between the selected alcohols, ketones and lanthanide chelates. The retention of over fifty solutes on five stationary phases containing 61 (M = Pr, Eu, Dy, Er, Yb, n = 3, R = Mef) dissolved in polydimethylsiloxane were later measured ". The initial motivation for this work was to explore the utility of a solvation parameter model proposed and developed by Abraham and coworkers for complexing stationary phases containing metal coordination centers. Linear solvation... 

Larrivee, M. L. and Poole, C. F. 1994. Solvation parameter model for the prediction of breakthrough volumes in solid-phase extraction with particle-loaded membranes. Anal. Chem., 66. 139-146. 

Hoover, K.R., Acree, W.E., Jr and Abraham, M.H. (2005) Chemical toxicity correlations for several fish species based on the Abraham solvation parameter model. Chem. Res. Toxicol., 18, 1497-1505. 

Calculation of solute descriptor values for the solvation parameter model... 

An example of part of the output for fitting the solvation parameter model to a reversed-phase chromatographic system by multiple linear regression analysis... 

Figure 1.3. Variation of the system constants of the solvation parameter model (section 1.4.3) with temperature for 37 % (v/v) propan-2-oI in water on the porous polymer PLRP-S stationary phase. The m constant reflects the difference in cohesion and dispersive interactions, r constant loan-pair electron interactions, s constant dipole-type interactions, a constant hydrogen-bond basicity and b constant hydrogen-bond acidity between the mobile and stationary phases. (From ref. [89] The Royal Society of Chemistry).

The master retention equation of the solvation parameter model relating the above processes to experimentally quantifiable contributions from all possible intermolecular interactions was presented in section 1.4.3. The system constants in the model (see Eq. 1.7 or 1.7a) convey all information of the ability of the stationary phase to participate in solute-solvent intermolecular interactions. The r constant refers to the ability of the stationary phase to interact with solute n- or jr-electron pairs. The s constant establishes the ability of the stationary phase to take part in dipole-type interactions. The a constant is a measure of stationary phase hydrogen-bond basicity and the b constant stationary phase hydrogen-bond acidity. The / constant incorporates contributions from stationary phase cavity formation and solute-solvent dispersion interactions. The system constants for some common packed column stationary phases are summarized in Table 2.6 [68,81,103,104,113]. Further values for non-ionic stationary phases [114,115], liquid organic salts [68,116], cyclodextrins [117], and lanthanide chelates dissolved in a poly(dimethylsiloxane) [118] are summarized elsewhere. 

System constants derived from the solvation parameter model for packed column stationary phases at 121°C... 

Figure 2.7. Principal component score plot with the system constants from the solvation parameter model as variables for 52 non hydrogen-bond acid stationary phases at 121°C. Loading for PC 1 0.996 a + 0.059 s -0.054 / - 0.024 c - 0.014 r. Loading for PC 2 0.940 s + 0.328 / + 0.080 r + 0.027 c - 0.037 a.

Figure 2.8. Nearest neighbor complete link cluster dendrogram for the stationary phases in Table 2.6. The system constants from the solvation parameter model were used as variables.

The system of stationary phase constants introduced by Rohrschneider [282,283] and later modified by McReynolds [284] was the first widely adopted approach for the systematic organization of stationary phases based on their selectivity for specific solute interactions. Virtually all-popular stationary phases have been characterized by this method and compilations of phase constants are readily available [28,30]. Subsequent studies have demonstrated that the method is unsuitable for ranking stationary phases by their selectivity for specific interactions [29,102,285-287]. The solvation parameter model is suggested for this purpose (section 2.3.5). A brief summary of the model is presented here because of its historical significance and the fact that it provides a useful approach for the prediction of isothermal retention indices. 

Rohrschneider s approach is able to predict retention index values for solute s with known solute constants (a, through e) [283,288]. These are determined from AI values for the solute on at least five phases of known phase constants and solving the series of linear equations. The retention index of the solute on any phase of known phase constants (X through S ) can then be calculated from Eq. (2.8). The theoretical defects of the method for assigning intermolecular interactions do not apply to the prediction of retention index values. A mean error of about 6 index units was indicated in some calculations. The retention or retention index values for thousands of compounds can be calculated from literature compilations of solute descriptors and the system constants summarized in Tables 2.6 and 2.8 using the solvation parameter model [103]. The field of structure-driven prediction of retention in gas chromatography is not well developed at present and new approaches will likely emerge in the future. 

Influence of solvent type on the system constants of the solvation parameter model for a cyanopropylsiloxane-bonded silica sorbent in reversed-phase chromatography (r = 0 in all cases)... 

The solvation parameter model has been extended to include ternary solvent systems using a mixture-design approach [271,272]. The system maps are now smooth three-... 

The cavity model of solvation provides the basis for a number of additional models used to explain retention in reversed-phase chromatography. The main approaches are represented by solvophobic theory [282-286] and lattice theories based on statistical thermodynamics [287-291]. To a lesser extent classical thermodynamics combining partition and displacement models [292] and the phenomenological model of solvent effects [293] have also been used. Compared with the solvation parameter model all these models are mathematically complex, and often require the input of system variables that are either unknown or difficult to calculate, particularly for polar compounds. For this reason, and because of a failure to provide a simple conceptual picture of the retention process in familiar chromatographic terms, these models have largely remained the province of the physical chemist. 

System constants of the solvation parameter model for liquid-solid chromatography... 

For chromatographic applications, the most useful models of solvent properties are the solubility parameter concept, Snyder s solvent strength and selectivity parameters, solvatochromic parameters and the system constants of the solvation parameter model for gas to liquid transfer. The Hildebrand solubility parameter, 8h (total solubility parameter), is a rough measure of solvent strength, and is easily caleulated from the physical properties of the pure solvent. It is equivalent to the square root of the solvent vaporization energy divided by its molar volume. The original solubility parameter concept was developed from assumptions of regular solution behavior in which the principal intermolecular interactions were dominated by dispersion forces. 

System constants from the solvation parameter model for transfer from the gas phase to the solvent at 25°C... 

---