Solvents solvation parameter model

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]. 

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. 

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. 

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... 

In order to quantify intermolecular solute-IL interactions, (Abraham et al., 2003 Acree Abraham, 2006) 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, (Sprunger et al, 2007 Sprunger et al., 2008 Sprunger et al., 2009a Sprunger et al., 2009b) modified the Abraham solvation parameter model ... 

It should be born in mind, however, that the activation parameters calculated refer to the sum of several reactions, whose enthalpy and/or entropy changes may have different signs from those of the decrystalUzation proper. Specifically, the contribution to the activation parameters of the interactions that occur in the solvent system should be taken into account. Consider the energetics of association of the solvated ions with the AGU. We may employ the extra-thermodynamic quantities of transfer of single ions from aprotic to protic solvents as a model for the reaction under consideration. This use is appropriate because recent measurements (using solvatochromic indicators) have indicated that the polarity at the surface of cellulose is akin to that of aliphatic alcohols [99]. Single-ion enthalpies of transfer indicate that Li+ is more efficiently solvated by DMAc than by alcohols, hence by cellulose. That is, the equilibrium shown in Eq. 7 is endothermic ... 

Here Vij denotes the distance between atoms i and j and g(i) the type of the amino acid i. The Leonard-Jones parameters Vij,Rij for potential depths and equilibrium distance) depend on the type of the atom pair and were adjusted to satisfy constraints derived from as a set of 138 proteins of the PDB database [18, 17, 19]. The non-trivial electrostatic interactions in proteins are represented via group-specific dielectric constants ig(i),g(j) depending on the amino-acid to which atom i belongs). The partial charges qi and the dielectric constants were derived in a potential-of-mean-force approach [20]. Interactions with the solvent were first fit in a minimal solvent accessible surface model [21] parameterized by free energies per unit area (7j to reproduce the enthalpies of solvation of the Gly-X-Gly family of peptides [22]. Ai corresponds to the area of atom i that is in contact with a ficticious solvent. Hydrogen bonds are described via dipole-dipole interactions included in the electrostatic terms... 

For molecular systems in the vacuum, exact analytical derivatives of the total energy with respect to the nuclear coordinates are available [22] and lead to very efficient local optimization methods [23], The situation is more involved for solvated systems modelled within the implicit solvent framework. The total energy indeed contains reaction field contributions of the form ER(p,p ), which are not calculated analytically, but are replaced by numerical approximations Efp(p,p ), as described in Section 1.2.5. We assume from now on that both the interface Y and the charge distributions p and p depend on n real parameters (A, , A ). In the geometry optimization problem, the A, are the cartesian coordinates of the nuclei. There are several nonequivalent ways to construct approximations of the derivatives of the reaction field energy with respect to the parameters (A1 , A ) ... 

Simple electrostatic models can be used to interpret the activity coefficients of polar molecules in terms of just three parameters a radius, the dipole moment of the solute, and the dielectric constant of the solvent. The continuum model of the solvent can be used to deduce a value for the free energy of solvation of a spherical molecule of radius r containing a point dipole at its center. The value obtained by Kirkwood from electrostatic theory is... 

Both AP and Ga have a tightly bound hydrate shell in aqueous solution and both are prone to hydrolysis. In terms of the Hertz electrostatic model for quadrupolar relaxation of ionic nuclei in electrolyte solution (see Section III.C) one therefore expects effective quenching of the electric field gradient caused by the surrounding water dipoles, due to a nearly perfect coordination symmetry. Any contribution to the e.f.g. should therefore arise from outer-sphere solvent dipoles. In terms of the fully orientated solvation (FOS) model this would correspond to a distribution width parameter approaching zero (/. -> 0) with the first term in equation (4) vanishing. This is indeed what Hertz (24) found for both AF" and Ga ", and the experimental infinite dilution relaxation rates ( AP" 7-5 s Ga 350 s ) are remarkably well matched by the computed ones... 

The various parameters have been fit to reproduce experimentap34-337 aqueous solvation data. Much like the earlier quantum models, the primary dependence of the ENP terms is on the solvent dielectric constant, which is taken from experiment. Cavity definition, regardless of shape, is parametric in every model, although many researchers avoid the term nevertheless, van der Waals radii, isodensity surface values, and so on are parametric choices. The more important point is that the cavity parameters are not expected to show much sensitivity to solvent in any model. 

The value of X for a typical polar solvent is approximately two. This equation was introduced earlier in the development of the MSA for ion-solvent interactions (section 3.5). It was seen that the MSA gives an improved description of ion solvation parameters with respect to the Born model. However, it fails to distinguish between the solvation of cations and anions of the same size. In other words, it fails to distinguish between the short-range chemical interactions which stabilize ions of differing charge. 

---