Scott K.
Spear
*a,
Scott T.
Griffin
b,
Kimberly S.
Granger
b,
Jonathan G.
Huddleston
c and
Robin D.
Rogers
*b
aAlabama Institute for Manufacturing Excellence, The University of Alabama, Tuscaloosa, AL 35487, USA. E-mail: SSpear@bama.ua.edu; Fax: +1 205 3483510; Tel: +1 205 3480464
bCenter for Green Manufacturing and Department of Chemistry, The University of Alabama, Tuscaloosa, AL 35487, USA. E-mail: RDRogers@bama.ua.edu.; Fax: +1 205 3480823; Tel: +1 205 3484323
cMillipore Bioprocessing Ltd., Medomsley Road, Consett, County Durham, UK DH8 6SZ
First published on 24th May 2007
The physico-chemical properties of soybean oil methyl ester (SBME), better known as biodiesel, of importance to its use as a solvent in liquid–liquid separations have been examined. Partition coefficients of several organic species between SBME–water have been determined and compared to log P (1-octanol–water). The free energy of transfer of a methylene group has been obtained and the solvent properties of the SBME–water system determined from distribution data of a small solute set using Abraham's generalized solvation equation. Solute distribution behavior is similar to that found for conventional organic solvent–water systems, but is most similar to other vegetable oils such as olive oil. When ionizable solutes are partitioned in the SBME–water system at differing pH, the neutral species show the highest distribution. Partitioning is dependent on the solute's ability to form hydrogen bonds between water and its charged state. Metal ions (e.g., Fe3+, Co2+, and Ni2+) exhibit moderate partitioning to the SBME phase from water only in the presence of extractants. Actinides (UO22+, Am3+) exhibit significant partitioning to the SBME from aqueous solutions with the use of octyl(phenyl)-N,N-diisobutylcarbamoylmethyl phosphine oxide (CMPO). Soybean oil methyl ester may be a suitable “green” alternative for the replacement of volatile organic solvents in liquid–liquid extractions in selected applications.
World production of vegetable oils is recorded at 88.7 million metric tons, with over 10 million metric tons of vegetable oils produced in the United States.7 U. S. oilseed crops include soybean (by far the predominant seed oil crop), corn, cotton, sunflower, flax, and rapeseed (canola). The major use for vegetable oils is, obviously, in food products, however, some vegetable oils and their ester derivatives are finding increasing industrial usage. Two important uses of transesterified soybean oil are in soybean oil inks and in diesel fuel.8 Also, vegetable oils have recently found use as a cheap, non-toxic solvent for the synthesis of nanoparticles9 and quantum dots.10 Regulatory pressure to reduce VOC emissions is leading to an increasing use of such solvents in the coatings industry.11
Soybean oil methyl esters (SBMEs, e.g., SoyGold® 1000) are composed of a mixture of fatty acids (Table 1). Palmitic and stearic acid are saturated fatty acids and comprise approximately 15% of the final mixture. The remaining three major fatty acids, oleic acid, linoleic acid, and linolenic acid have varying degrees of unsaturation and comprise approximately 83% of the final mixture. By far the major fatty acid present in soybean oil is linoleic acid. Linoleic acid is an 18-carbon long-chain fatty acid with double bonds at the 9- and 12-carbon atoms.
SBME is produced by the transesterification of soybean oil and methanol. It is non-toxic, non-hazardous, and biodegradable. It is neither a hazardous air pollutant (HAP), nor an ozone-depleting chemical (ODC).
Transesterification chemistry (used in the production of SoyGold® 1000) is a more modern innovation, and involves the combination of organically derived oils or fats with a small chain alcohol such as methanol or ethanol in the presence of a catalyst (most often the catalyst is a base) to form fatty acid esters. This process is akin to saponification chemistry, the alkaline hydrolysis of fatty acid esters, in which triglycerides are reacted with sodium or potassium hydroxide to produce glycerol and fatty acid salts (soap). SBME has a high flash point (218 °C) and low vapor pressure (<2 mm Hg).13
Materials such as SBME may have an important role to play in replacing the previous generation of solvents with materials of a more benign nature derived from renewable resources. However, there is a need to understand the physico-chemical properties of these materials before they can find widespread application. The solvent property data we report here is intended to inform that understanding, and begin to place SBME in context with other common solvents.
The 14C-labeled organic solutes were purchased from Sigma (St. Louis, MO, USA). Upon receipt, the tracers were diluted to an activity of ca. 0.06–0.08 µCi µL−1 for use as the ‘spike’ in the partitioning experiments. The hydrophilic tracers were diluted in deionized water and the hydrophobic tracers were diluted in their unlabelled form.
The cobalt tracer (as 60CoCl2 in 0.5 M HCl) was obtained from New England Nuclear (Boston, MA, USA). The iron tracer (as 59FeCl3 in 0.1 M HCl) and the nickel tracer (as 63NiCl2 in 0.1 M HCl) were both obtained from Amersham Pharmacia Biotech (Piscataway, NJ, USA). The actinide tracers (233UO2Cl2 and 241AmCl3) were obtained from Isotope Products Laboratories (Valencia, CA, USA). All metal ion radiotracers, except for 233U, were used as received for partitioning studies, or diluted with deionized water to an activity of ca. 0.03 µCi µL−1.
Stock solutions of the uranyl nitrate radiotracer solution were cleaned of impurities (228Th and its daughters from decay of 232U impurities) using TEVA resin (EIChrom Industries, Darien, IL, USA).14 20 µL of 233UO22+ stock solution was diluted in 10 mL of 5 M HCl and the solution was added to a column comprised of 1 g TEVA resin preconditioned with 5 M HCl. After loading, 20 mL of 0.5 M HCl was added to the column to elute 233U, which was collected, evaporated to dryness and dissolved in 0.1 M HNO3.
The distribution ratios reported here are the average of at least two measurements and are typically accurate to ±5%. Where necessary, the aqueous phase pH was measured using a Fisher Accumet pH meter and adjusted to acidic conditions using sulfuric acid or to basic conditions using sodium hydroxide. Multiple linear regression analyses were performed using Stat Box v2.5 (Grimmer Logiciels, 34 Rue de Dunkerque, Paris, France).
The relationship between the octanol–water partition coefficient and the distribution in the SBME–water system can be summarized in the form of eqn (1):
Log D (SBME) = b + m(log P) | (1) |
Fig. 2 shows the distribution of several n-alcohols in relation to their alkyl chain length in the SBME–water system, as well as literature data for their distribution in an olive oil–water system16 and the conventional 1-octanol–water system. These relationships can be described by eqn (2):
ln K = C + Enc | (2) |
ΔGCH2 = −RTE | (3) |
![]() | ||
Fig. 2 Distribution of short chain alcohols in the SBME–water (▲) system compared to distribution in the olive oil–water (▼) and 1-octanol–water (●) systems. |
Table 2 presents ΔGCH2 values for various solvent systems to allow comparison with the present SBME system and with the olive oil–water system. Also included in Table 2 are ΔGCH2 values for selected ionic liquid (IL) salt–salt20 and polyethylene glycol (PEG)–salt aqueous biphasic systems (ABS)4 that have been studied in our laboratory. From Table 2, one notes that based on the ΔGCH2 for the various solvents listed, SBME is very near to that observed for xylene and would suggest that SBME is a likely alternative to xylene as a solvent. Indeed, Hendrickson et al. have discovered that SBME is a suitable replacement for xylene as a developing solvent for photopolymerizable printing plates.21
Solvent–water system | −ΔGCH2/kcal mol−1 | Ref. | Solvent–water system | −ΔGCH2/kcal mol−1 | Ref. |
---|---|---|---|---|---|
a Present work. b Calculated in present work with data from Ref. 16. | |||||
Hexane | 1.10 | 18 | Diisopropyl ether | 0.68 | 18 |
Chloroform | 0.85 | 18 | SBME | 0.67 | —a |
Benzene | 0.84 | 18 | Xylene | 0.64 | 18 |
Octanol | 0.79 | 18 | Olive oil | 0.61 | —b |
Octane | 0.77 | 18 | n-Butanol | 0.54 | 18 |
Dodecane | 0.77 | 18 | Methyl ethyl ketone (MEK) | 0.43 | 18 |
Methyl isobutyl ketone (MIBK) | 0.72 | 18 | PEG–salt (general) | 0.142–0.696 | 4 |
IL salt–salt : 1-butyl-3-methylimidazolium chloride/K3PO4 or K2HPO4 or K2CO3 | 0.304–0.728 | 20 | PEG-2000 (40% w/w)/(NH4)2SO4 (1.68 M) | 0.214 | 4 |
Table 2 and Fig. 2 illustrate the similarity between the solvent properties of the currently investigated SBME–water system and literature data for the olive oil–water system. The free energy of transfer between the phases of the SBME system is much less than that typical of solvents dominated by van der Waals forces and much more similar to solvents typified by molecular interactions or by a significant water content at equilibrium. However, the equilibrium water content of the SBME phase is relatively low (0.35% wt/wt), and since it might be expected that the predominant intermolecular forces in an SBME phase would be dispersion forces, it is possible that molecular packing features contribute to the relative weakness of cavity forces in the olive oil and SBME phases implied by the relatively low free energy of CH2 transfer value obtained in these systems.
Some Property (SP) = cavity terms + polarity terms + hydrogen bonding terms + constant | (4) |
Linear free energy relationships (LFER) based upon the generalized solvation equation have been widely used to model many processes, such as partitioning in aqueous–organic systems, solubility, and transport across biological membranes.22–24 Abraham's generalized solvation equation, eqn (5):
log SP = c + rR2 + sπ2H + a∑α2H + b∑β2H + vVx | (5) |
In essence, the log of SP, in this case the partition coefficient of a series of solutes in a given system, can be related to the solute property descriptors of each solute through eqn (5). The solute distribution ratios for a very limited set of solutes determined in the SBME–water system and their corresponding Gibbs's energy related solute property descriptors are shown in Table 3.24,26,31
Solute | Abraham parameters | log D | |||||
---|---|---|---|---|---|---|---|
π 2 H | ∑α2H | ∑β2H | Vx | R 2 | Ref. | ||
1,2,4-Trichlorobenzene | 0.081 | 0 | 0 | 1.0836 | 0.98 | 23 | 2.99 |
1,4-Dichlorobenzene | 0.75 | 0 | 0.02 | 0.9612 | 0.825 | 23 | 2.75 |
Chlorobenzene | 0.65 | 0 | 0.07 | 0.8388 | 0.718 | 23 | 2.41 |
Toluene | 0.52 | 0 | 0.14 | 0.8573 | 0.601 | 23 | 2.33 |
p-Toluic acid | 0.90 | 0.60 | 0.38 | 1.073 | 0.73 | 25 | 1.28 |
Salicyclic acid | 0.84 | 0.71 | 0.38 | 0.9904 | 0.89 | 26 | 1.39 |
Benzene | 0.52 | 0 | 0.14 | 0.7164 | 0.61 | 23 | 2.23 |
Acetophenone | 1.01 | 0 | 0.49 | 1.0139 | 0.818 | 23 | 1.54 |
Methyl iodide | 0.43 | 0 | 0.13 | 0.5077 | 0.676 | 23 | 1.27 |
Pentanol | 0.42 | 0.37 | 0.48 | 0.8718 | 0.219 | 23 | 0.52 |
Butanol | 0.42 | 0.37 | 0.48 | 0.7309 | 0.224 | 23 | −0.096 |
Phthalic acid | 1.60 | 0.82 | 0.75 | 1.147 | 0.85 | 26 | −1.24 |
Propylamine | 0.35 | 0.16 | 0.61 | 0.631 | 0.225 | 23 | 0.41 |
Propanol | 0.42 | 0.37 | 0.48 | 0.59 | 0.236 | 23 | −0.39 |
Isopropanol | 0.36 | 0.33 | 0.56 | 0.59 | 0.212 | 23 | −0.87 |
Acetic acid | 0.65 | 0.61 | 0.45 | 0.4648 | 0.265 | 23 | −1.26 |
Ethanol | 0.42 | 0.37 | 0.48 | 0.4491 | 0.246 | 23 | −1.07 |
Acetonitrile | 0.90 | 0.04 | 0.33 | 0.4042 | 0.237 | 23 | −0.81 |
Methanol | 0.44 | 0.43 | 0.47 | 0.3082 | 0.278 | 23 | −1.44 |
The descriptors are the McGowan volume (Vx), the excess molar refraction (R2), the solute dipolarity–polarizability (π2H), and the solutes' overall hydrogen-bond acceptor basicity and hydrogen-bond donor acidity, ∑β2H and ∑α2H, respectively. The sign and magnitude of the regression coefficients obtained from eqn (5) (r, s, a, b, and v), by multiple linear regression of the solute descriptors on the log of the partition coefficient, may be considered proportional to the solvent properties of the phases corresponding to the appropriate solute descriptor. Thus, r corresponds to the relative strengths of the solute–solvent interactions determined by the excess molar refractivity of the solute, and a corresponds to the relative solvent hydrogen-bond basicity of the phases, and so on. Finally, c is a constant of proportionality.
The use of as few as twenty solutes in a multiple linear regression is unlikely to give very satisfactory results, and thus the results we report can only be considered a preliminary characterization of the solvent properties of the SBME–water partitioning system. Nevertheless, Table 4 shows the relative contributions of the regression coefficients to the partitioning of the solutes in the SBME–water system. Only solute acidity, solute basicity, and solute volume are significant at the 95% level. The most significant factor is solute basicity, for which the regression coefficient is large and negative. The SBME phase is significantly less acidic than the aqueous phase and basic (in the sense of Lewis) solutes have a strong preference for the aqueous phase.
Variable | Coefficient | Correlation/Y | Student's t | Probability |
---|---|---|---|---|
a Adjusted R2 = 0.90; multiple R2 = 0.96; F statistic = 55.34 | ||||
a | −1.67 | −0.61 | 2.87 | 0.0117 |
b | −3.61 | −0.80 | 4.80 | 0.0002 |
v | 3.20 | 0.59 | 6.98 | 4.41× 10−6 |
Constant | −0.012 |
The next most significant parameter is volume, which is comparatively small and (typical of solvent–water partitioning) positive. The v parameter may be considered to be a measure of the relative hydrophobicity of the system and is similar to ΔGCH2 in that it reflects the difference in the free energy required for cavity formation between the water and SBME phases. Typically, this parameter is large in solvent–water systems (Table 4) reflecting the fact that for most solvents, there are few intermolecular forces hindering cavity formation other than dispersion forces, whereas the free energy of cavity formation in aqueous phases is dominated by its extensively hydrogen bonded nature. There is a significant penalty to cavity formation in the aqueous phase due to the structuredness of water. However, the magnitude of this parameter is not as great as for most pure van der Waals solvents, which usually indicates either there is significant structure to the SBME phase, or that (and of course it amounts to the same thing) there is significant water present in the organic phase at equilibrium; for example, in the 1-octanol–water system. However, in the case of SBME, we have measured an equilibrium water content of only 0.35% wt/wt. It may be thought in the case of olive oil and SBME that molecular shape features play a role in determining the magnitude of this parameter.
Finally, the a coefficient is relatively small and negative, suggesting that the oil phase is less basic than the aqueous phase, it is also the least significant factor of those which are significant at the 95% confidence level. Again, it should be remembered that the solute set is small and the above interpretation tentative.
The full solvation equation for the 20 solute SBME–water partition is given in eqn (6):
Log D(SBME) = −0.11 − 0.78π2H − 1.58∑α2H − 2.90∑β2H + 3.19Vx − 0.59R2 F = 35.5 | (6) |
The coefficients are slightly different from those shown in Table 4 because in Table 4 only three terms are used (a, b, and v), and in eqn (6) all 5 terms have been used (a, b, v, s, and r). Both equations give the same result for Log D within statistical limits; however, the effect of the extra terms is to change the values of the terms in the first equation. The effect ought to be small because the extra terms are of small statistical significance, nevertheless, they do have an effect.
These results may be directly compared to a similar LSER derived for the olive oil–water system27 shown in eqn (7).
Log D(OO) = −0.011 − 0.8π2H − 1.47∑α2H − 4.92∑β2H + 4.17Vx + 0.58R2 F = 5841 | (7) |
The similarity in the roles of solute acidity, basicity, and volume in the SBME–water systems and the olive oil–water system is readily apparent. Surprisingly, even terms which were least significant in the regression of eqn (6), follow the olive oil–water regression quite closely.
![]() | ||
Fig. 3 pH dependent partitioning of ionizable solutes: 1 acetic acid (pKa = 4.75) ; 2 phthalic acid (pKa1 = 2.89, pKa2 = 5.51); 3 propylamine (pKa = 10.71); 4 salicylic acid (pKa1 = 2.97, pKa2 = 13.40); 5 p-toluic acid (pKa = 4.36). |
![]() | ||
Fig. 4 The distribution of metal ions (Co2+, Fe3+, Ni2+) in the absence and presence of PAN and TAN (each at 0.1 mM) vs. pH. |
For the transition metal ions (Co2+, Ni2+, Fe3+) examined, no useful extraction is achieved in the absence of an extractant. In the presence of the extractants PAN and TAN, the distribution coefficient of all the metal ions increases, except at the lowest pH. This increase is pH dependent, being highest under the most basic conditions. This is commonly found in the extraction of metal ion species by solvent extraction and is due to the increase in the complexation constant with pH and forms a useful technical means of metal ion purification. Extractants and organic phases may be loaded at high pH and unloaded at low pH in a similar forward and back extraction process as mentioned above for the extraction of organic acids and bases. Additionally, some selectivity in metal ion separation may be secured by this means. In the SBME–water system, there appears to be little difference between the performance of PAN and TAN which are most effective in the extraction of Co2+, but appear not practically useful for the extraction of Ni2+ and Fe3+. Undoubtedly, suitable extractants could be found for the extraction of these and other metal ion species using the SBME–water system.
In the TRUEX process for the extraction of uranium and transuranic ions, 0.2 M CMPO is used as the extractant, 1.2–1.4 M tri-n-butylphosphate (TBP) as a phase modifier, and the diluent is paraffinic hydrocarbons.34 Schultz and Horwitz,34 using the TRUEX formulation, reported the highest distribution values for UO22+ and Am3+ to be 103 and 10, respectively, at 6 M HNO3. Mathur et al., using dodecane as the diluent, 0.2 M CMPO and 1.2 M TBP, reported distribution ratios of 102 for UO22+ and 10 for Am3+ at 1–6 M HNO3.35 We have examined the partitioning of UO22+ (Fig. 5) and Am3+ (Fig. 6) in a similar system but with SBME as a replacement for the paraffinic hydrocarbon diluent. Without CMPO, the distribution ratios fall well below 1, indicating the preference for the aqueous phase. For the uranyl cation, addition of 0.2 M CMPO significantly improves the distribution ratios up to 102 at all nitric acid and NaNO3 concentrations studied.
![]() | ||
Fig. 5 The distribution of UO22+ in the SBME–water system without extractant (●), with 0.2 M CMPO (■), and with 0.2 M CMPO–1.2 M TBP (▲) vs. HNO3 (filled symbols, solid lines) or NaNO3 (open symbols, dashed lines). Distribution ratio from SBME–water = 0.0044. |
![]() | ||
Fig. 6 The distribution of Am3+ in the SBME–water system without extractant (●), with 0.2 M CMPO (■), and with 0.2 M CMPO–1.2 M TBP (▲) vs. HNO3 (filled symbols, solid lines) or NaNO3 (open symbols, dashed lines). Distribution ratio from SBME–water = 0.0082. |
As with the uranyl cation, the distribution ratios for Am3+ are well below 1 and although there is an increase in preference for the SBME phase with the addition of CMPO, the distribution ratios are still below 1 and favor the aqueous phase up to 0.1 M NO3−, where they then begin to increase above 1 and favor the SBME phase approaching 102 at 5 M NO3−. This behavior would allow the separation of UO22+ from Am3+ in this system.
The inclusion of 1.2 M TBP for either UO22+ or Am3+ only made small improvements in the distribution ratios and does not appear to be needed as a phase modifier in the SBME–water systems as observed for traditional dodecane–water systems;34,35 the distribution values remain comparable without its use. The elimination of TBP could simplify and reduce the costs of these systems.
This journal is © The Royal Society of Chemistry 2007 |