Abstract
This study separates tetrahydrofuran (THF)-water mixtures containing varying THF percentages, using sorbitan esters (Spans) with different tail characteristics. We probe the separation mechanisms using attenuated total reflectance-Fourier transform infrared spectroscopy and small angle X ray scattering (SAXS). THF and water are miscible and interact through hydrogen bonds. Water splits the COC absorbance band of THF into a peak at ≈1,070 cm−1 (crystalline THF) and a dominant peak at ≈1,050 cm−1 (glassy THF), indicating disorder. Depending on the Span, separation occurs for mixtures containing up to 70% THF (v/v, relative to water). Spans with unsaturated tails separate the lowest THF percentages. Tail length and number of Span tails enhances ordering of THF, and the crystalline THF peak at ≈1,070 cm−1 dominates. Spans interact with THF through hydrogen bonds, as reflected in the splitting of the COC band of THF. Furthermore, C-H…O hydrogen bonds cause a blueshift in the νas(CH2) band of Spans with increasing THF. This effect is most significant in Span 40 and 60, indicating that they interact with THF more strongly than Span 20, Span 80 and Span 85. In contrast, they interact with water less strongly than Span 20, Span 80 and Span 85, as indicated by their flocculation at low THF percentages. Therefore, we propose that separation between THF and water occurs primarily through two mechanisms: 1) Span 20, Span 80 and Span 85 compete against THF for interactions with water through their hydrophilic head, and 2) Span 40 and Span 60 preferentially interact with THF through their tails. Nonetheless, water also interacts with the heads of Span 40 and Span 60, as indicated by SAXS. SAXS shows that in THF Spans self-assemble into micelles, which aggregate into either surface fractals or mass fractals. There are two persistence lengths because of the limited order in THF. Water orders self-assembled structures, likely by favoring the formation of micelles which host water in their interior. Therefore, we identify a single persistence length (≈25 Å), representative of the distance between the micelle centers.
1 Introduction
THF, acetonitrile, dimethylformamide or methanol are water miscible pollutants used in industrial processes (; ; ; ; ; ; ). They are found in industrial wastewater as well as in groundwater, in the case of spills such as those caused by the notable Ohio trail derailment accident (2023). Our study focuses on THF, which is used as a stabilizer for chlorinated solvents, as well as to produce pharmaceuticals and pesticides (; ).
Beyond their intrinsic toxicity, spills of water miscible solvents are particularly problematic because they can rapidly migrate in impacted aquifers, harming downstream receptors (). Furthermore, their separation from water presents greater challenges compared to water-immiscible pollutants, such as hydrocarbons (). These can be removed from water through mechanical separation methods (e.g., coalescers ()) or chemical methods (such as oxidation (; ; )). Also, when they are present in groundwater below levels at which they are toxic to the bacteria, they can biodegraded (; ; ). Since biodegradation tends to be a fairly slow process, there are limitations to its use in the case of rapidly migrating pollutants. These consideration also apply to phytoremediation using plants (; ). Moreover, phytoremediation is only effective for shallower contamination, which can be reached by the roots of the plants. In the case of miscible organic solvents, pervaporation is often used (). Beyond its high energy cost, its implementation for the treatment of impacted aquifers is not feasible. This is because complex treatment units could be integrated in a chemical plant, but they cannot be built as part of a remediation plant on polluted sites. Therefore, it is necessary to identify cost-effective, facile methods to separate miscible solvents from either industrial wastewater or groundwater. These methods were scant until our most recent research.
Previous studies used CO2 switchable species (Jessop and Cunningham) or increased the hydrophobicity of the organic solvent () to increase their incompatibility with water, thereby leading to separation. Other studies used species that preferentially interact with water, at the expense of its interactions with the solvent. Examples of these species are sugars and their alcohols (; ; ) as well as salts, such as sodium chloride and choline salts) (; ).
In our recent research, we have used amphiphiles to separate miscible solvents from water. Examples include sodium lauroyl lactylate (), glycerol monooleate (), fatty acids (; ; ) and a sorbitan ester (Span 80, (),). In these studies, we showed that separation involves preferential interactions either between water and the amphiphile [e.g., in the case of Span 80, ()] or preferential interactions between the solvent and the amphiphile [e.g., in the case of stearic acid, ()]. To investigate interactions between water and other components in the mixture (solvent and amphiphiles), we probed the OH stretch band of water using Attenuated Total Reflectance–Fourier Transform Infrared (ATR-FTIR) spectroscopy. We deconvolved the OH stretch peak into four peaks, representing water species with different coordination to other molecules in the mixtures. Specifically, we considered double (DD) and single (SD) hydrogen (H) bond donors, and double (DA) and single (SA) H bond acceptors. Based on this analysis, we concluded that separation is promoted when amphiphiles interact with water the same way the solvent does, i.e., when amphiphiles compete with the solvent for interactions with the same water species. We also analyzed H bonding of solvent with hydrophobic amphiphiles, such as stearic acid (). We showed that the COC band of THF contains information about H bonding, by comparing changes of the OH stretch band of water and the COC band of THF, in THF-water mixtures (). By analyzing this band, we demonstrated that stearic acid separates water from THF because it competes for interactions with the same THF species, i.e., THF molecules that H bond and coordinate the same way with other molecules in solution ().
Furthermore, we highlighted that the tail plays an important role in separation. For example, separation between THF and water does not occur upon adding acetic acid, although the structure of acetic acid is the same as the head of stearic acid (). Given this similarity, both acetic acid and stearic acid H bond with THF the same way (). We proposed that the tail of stearic acid promotes H bonding between THF and the hydrophilic tail at the expense of THF-water interactions. In contrast, in the absence of the tail, acetic acid can H bond either water or THF, thereby enabling the existence of a ternary solution. Indeed, solubility plays a key role. Specifically, separation is favored when the amphiphile added has intermediate solubility in the mixture, and is soluble in either water or the solvent, but not both (). In ternary mixtures where all components are miscible two by two, there are instances where separation may still occur, but predictions of the mixing behaviour are notably more complex (). Solubility of the amphiphile in each component and in the mixture is clearly affected by the tail. While our previous studies started highlighting the role of the tail, we had not yet conducted a systematic analysis of the role of the tail on solvent separation. The current study aims to fill this research gap.
In particular, the objective of this study is to expand on our previous research, comparing the separation between THF and water using sorbitan esters with a different tail structure (e.g., tail length and hence hydrophobicity, and double vs. single bonds in the tail). Specifically, we investigate interactions (H bonding) between THF and different sorbitan esters, comparing them with interactions between THF and water. We also analyze the role of the tail on THF ordering and its separation from water.
This objective supports the overall goal of our research, which is to develop a robust approach to predict the separation of water from miscible solvents, and to identify the most effective amphiphiles for water purification. The toolbox we are building relies on calculated values for H bonds between solvent-water couples, amphiphile-water couples, as well as amphiphile-solvent couples. We are calculating these values in our research. Values estimated on representative amphiphile and solvents will be used to provide initial estimated for other similar compounds, predicting the mixing behaviour of ternary mixtures. Our toolbox also relies on building a systematic and quantitative understanding of the effect of the tail characteristics and amphiphile solubility. The current study will be followed by others, where broader groups of amphiphiles will be used to consolidate the portion of our predictive toolbox focused on the tails. Our approach will be complementary to models such as COSMO, which we recently used to predict activity coefficients, because they are correlated to solvent separation ().
2 Materials and methods
2.1 Materials
Sorbitan esters (Span 20, Span 40, Span 60, Span 65, Span 80 and Span 85) were purchased from Sigma Aldrich (Canada). Their structure and hydrophilic lipophilic balance (HLB) number are given in Table 1 and Scheme 1. THF (reagent grade, Caledon laboratory chemicals) was purchased from Fisher Scientific (Canada). All experiments were conducted using deionised (DI) water.
TABLE 1
| Surfactant name | HLB | Tail length |
|---|---|---|
| Span 20 (sorbitan monolaurate) | 8.6 | 11 |
| Span 40 (sorbitan monopalmitate) | 6.7 | 15 |
| Span 60 (sorbitan monostearate) | 4.7 | 17 |
| Span 80 (sorbitan monooleate) | 4.3 | 17 |
| Span 65 (sorbitan tristearate) | 2.1 | 17 on each tail, 3 tails |
| Span 85 (sorbitan trioleate) | 1.8 | 17 on each tail, 3 tails |
HLB number and tail length of the Spans used in this study.
SCHEME 1
2.2 Bottle tests
Samples were prepared using 40 g/L of sorbitan esters and mixtures containing varying percentages of DI water and THF (0%–100% THF, v/v relative to water). Samples were agitated by hand for 10 s after heating to 65°C for 2 min, using a hotplate. Samples were then cooled and observed at 20°C under quiescent conditions for 24 h, to qualitatively determine their phase behavior. These same samples were also imaged using optical microscopy (Section 2.3) and analyzed using ATR-FTIR (Section 2.4).
2.3 Optical microscopy
A VHX-5000 digital microscope (Keyence) was used to image THF-water emulsions stabilized with different sorbitan esters.
2.4 Attenuated total reflectance–Fourier transform infrared spectroscopy (ATR-FTIR)
Absorbance spectra were collected using an ATR-FTIR spectrometer (Thermoscientific Nicolet Summit FTIR spectrometer with an Everest ATR), with an accompanying IR solution software (Omnic 9, Thermo Fisher Scientific). Each spectrum was the average of 10 scans, with a resolution of 2 cm−1. The wavenumbers ranged from 400 cm−1 to 4,000 cm−1. Measurements were conducted at least in duplicate.
Spectra were analyzed using Quasar 1.5.0 (Orange-Spectroscopy software), using a rubber band baseline correction. Additionally, the spectra where normalized using min-max normalization. To this end, each spectrum was divided by its Absmax–Absmin range.
2.5 Synchrotron based small angle X-ray scattering (SAXS)
SAXS experiments were carried out at the Canadian Light Source Synchrotron (CLS) on the Brockhouse Diffraction Sector Undulator Beamline (BXDS-IVU) (), to identify the self-assembly of sorbitan esters into either pure THF or THF-water mixtures containing 2% water (relative to THF). In this study, samples were introduced into quartz tubing (Polyamide 0.0575″ ID x 0.0615″ OD), which was sealed to avoid solvent evaporation during data collection. A clay compound was used for sealing. SAXS data was collected with a Rayonix MX300HE detector (8,192 × 8,192 pixels) with 2 × 2 binning (4,096 × 4,096 pixel) for an effective pixel size of 73.242 µm with background stability mode ON. SAXS patterns were collected with a photon energy of 12.18 keV and sample-to-detector distance of ∼233 cm. SAXS patterns were collected in transmission geometry with a 120 s dwell time. Patterns were processed with GSASII (Argonne National Laboratory (C), 2010). This product includes software developed by the UChicago Argonne, LLC) (; ). SAXS data was calibrated with silver behenate (AgBeh) and instrument parameters such as sample-to-detector distances, detector tilt, beam centre, were refined as described elsewhere (). SAXS patterns were integrated from q = 0.012 to q = 0.637 Å−1. Measurements were done in duplicate for each sample.
Absolute intensity with background removed as a function of the scattering vector, q, were fitted to a model recently developed by our group. It has the general form,
The scattering intensity is interpreted as arising from scale-independent fractal scattering in the low q-range, superimposed on top Guinier’s form factor, namely a zero-mean Gaussian. Any periodicities arising from other correlations in the system are captured by non-zero mean Gaussians. In our case, we only required one zero mean Gaussian and one non-zero mean Gaussian.
The fractal dimension is captured in the modeling process by writing the scattering intensity, such as
This is a simple power-law decay on a double logarithmic scale within a given q-range limited by reciprocal fractal outer and inner cutoffs, qoc and qic, respectively, i.e., when qoc < q < qic. It is also known in the literature as the fractal region. Here, q is the module of the scattering vector, Ds, Dm and Dp are the surface, mass and respectively the pore fractal dimensions, and d is the Euclidean dimension of the space in which the fractal is embedded. In the case of surface fractals d−1 < Ds < d, Dm = Dp = d and the scattering intensity reduces to I(q)∼qDs−2d. For mass fractals 0 < Dm < d, Dm = Ds and Dp = d, and I(q)∼q-Dm, while for pore fractals Dp = Ds and Dm = d, and thus I(q)∼q−Dp. Therefore, one can obtain the fractal dimension from the slope (s) of the experimental scattering curve. Moreover, one can differentiate between mass and surface fractals, i.e., if the measured slope is s < d then the sample is a mass fractal, while if d < s < d + 1, the sample is a surface fractal. For a surface fractal, Ds is equal to 2 for a perfectly smooth surface, and approaches 3 for a highly folded/convoluted surface.
The term qp corresponds to the center of the Gaussian peak, while σ corresponds to the standard deviation of the Gaussian, which characterizes the width of the distribution.
The model was fitted to the data by non-linear regression using Prism 10.1 (GraphPad Software, San Diego, CA). We used a standard nonlinear regression routines (; ), and data were weighted by 1/Y2. This was done to ensure that solutions were stable and in order to determine confidence intervals.
The center position of the Gaussian peaks in reciprocal space was converted to its corresponding value in direct space, dp, by dp = 2π/qp. To determine the domain size, correlation length, coherence length, sometimes referred to the persistence length of the scattering event related to the Gaussian peaks we used an approximation based on the standard Scherrer model:
The term refers to the full width half maximum of a peak. The conversion of a standard deviation to a full-width half maximum is , while the shape factor is usually assumed to be K = for spherical domains of cubic symmetry. The correlation length can then be calculated from the standard deviation of the Gaussian using:
2.6 FTIR spectromicroscopy
FTIR spectromicroscopy (in ATR mode) was used to analyze emulsions prepared using different sorbitan esters and varying percentages of THF and water.
Measurements were conducted at the Canadian Light Source (CLS) in Saskatoon (Canada) using an Agilent Cary 670 FTIR Interferometer with Cary 620 Microscope, employing a 2D Focal Plane Array MCT-A detector. The hyperspectra were collected in transmission mode by placing 0.5 µL in between two calcium fluoride windows separated by a 10 μm Teflon spacer. This technique was used to map the samples, characterizing compositional differences between the interior and the exterior of emulsified droplets.
3 Results and discussion
THF and water are freely miscible in one another in binary solutions. Our previous study revealed that Span 80, a sorbitan ester, separates THF from water (). Here, we examine THF-water separation with different sorbitan esters (cf. Table 1), to assess the effect of their structure, and their tail in particular, on the mechanisms of solvent-water separation.
Note that in our study, separation was qualitatively assessed through bottle tests, optical microscopy and FTIR spectromicroscopy (Table 2). Our previous research conducted using Span 80 quantitatively showed that this sorbitan ester increases the water purity, but does not reduce THF to values compatible with legislative limits. These are for instance 50–1,300 μg/L in different states in the United States (). Span 80 is also one of the sorbitan esters analyzed in our current study. A modeling study reports that swing distillation could theoretically achieve up to 99.9% removal of THF from water (). Swing distillation is, however, energy intensive. Another study conducted using electrooxidation reports that THF removal was higher than 90% with suitable electrodes (). However, electrokinetic methods also use energy and require the installation of electrodes in impacted aquifers, which can be challenging and costly. The goal of our study is to reduce the concentrations of THF in water using low energy approaches, to then allow for subsequent treatments to be effective. As an example, a study reported the successful removal of THF from water using anaerobic bacteria and membranes (). Bacterial remdiation is energy effective, but bacteria can remediate pollutants only below concentrations at which they are toxic to them. Therefore, our approach can represent the first step in a treatment train for the removal of THF from water. Identifying the mechansism of separation and the surfactant characteristics that lead to effective separation will allow us to optimize this first step.
TABLE 2
| Tail | THF % (v/v, relative to water) | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Span | HLB | Number of carbons | Number of tails | Bonds | 10 | 20 | 30 | 40 (L) | 50 (L) | 60 (L) | 70 (L) | 80 (L) | 90 (L) |
| Span 20 | 8.6 | 11 | 1 | Single | E | E | E | 2 | 2 | 2 | 1 | 1 | 1 |
| Span 40 | 6.7 | 15 | 1 | Single | 1L + F | 1L + F | 1L + F | 2 | 2 | 2 | 2 | 1 | 1 |
| Span 60 | 4.7 | 17 | 1 | Single | 1L + F | 1L + F | 1L + F | 2 | 2 | 2 | 2 | 1 | 1 |
| Span 80 | 4.3 | 17 | 1 | Double | E | E | E | 2 | 2 | 2 | 2 | 1 | 1 |
| Span 65 | 2.1 | 17 | 3 | Single | 1L + F | E | E | 2 | 2 | 2 | 2 | 2 | 1 |
| Span 85 | 1.8 | 17 | 3 | Double | E | E | E | 2 | 2 | 2 | 2 | 2 | 1 |
Summary of the different phases observed in THF-water mixtures, with different THF to water ratios, and with different sorbitan esters. This table reflects observations with the naked eye and with an optical microscope, over a 24 h period. The colors of the cells distinguish different phase behaviours, and are used as visual aid. 1L = one liquid phase, based on observations with the naked eye and under an optical microscope; 1L + F = one liquid phase, based on observations with the naked eye and under an optical microscope, mixed with flocculated surfactant; E= emulsion, as seen under an optical microscope; 2L = separation into bulk phases, as seen with the naked eye in bottle tests.
Bottle tests and optical microscopy show that separation occurs at THF to water ratios dependent on the type of sorbitan ester (cf. Table 2; Figure 1; Supplementary Figures S1.1–S1.3). Specifically, we observe that both the hydrophobic lipophilic balance (HLB) number and saturation of the tail affect the THF to water ratio at which separation occurs.
FIGURE 1
While all sorbitan esters are hydrophobic (soluble in THF, insoluble in water), Span 20 has the highest HLB number (HLB = 8.6). It yields emulsions with 10%–30% THF (v/v, relative to water), free layers with 40%–60% THF and one liquid phase with 70%–80% THF. Span 40 (HLB = 6.7), Span 60 (HLB = 4.7) and Span 80 (HLB = 4.3) have intermediate HLB numbers. They separate THF from water into free layers with THF percentages ranging from 40% to 70%, while yielding one liquid phase only at 80% THF. With Span 40 and Span 60, separation does not occur with 30% THF or less, where the observed turbidity is due to flocs (Supplementary Figure S1.2). With Span 80, emulsions are observed at 30% THF or below (Supplementary Figure S1.2). Although Span 80 and Span 60 have similar HLB and the same tail length, Span 80 differs from Span 60 in that it has a double bond in its tail (i.e., its tail is unsaturated). Finally, Span 65 (HLB = 2.1) and Span 85 (HLB = 1.8) have the lowest HLB number, and they yield free phases with THF percentages from 40% to 80% THF. With 20%–30% THF, they yield emulsions. Finally, with 10% THF and Span 65, the observed turbidity is due to flocs (Supplementary Figure S1.2), while emulsions form with Span 85. Note that, similar to Span 80 and Span 60, Span 65 and 85 have similar HLB and the same tail length for each of their three tails. However, Span 85 has double bonds in its tails, while Span 65 does not. In both cases, double bonds in the tail (i.e., unsaturated tails) promote separation at the lowest THF percentages.
In our previous study, we explained water-THF separation by Span 80 based on the ability of Span 80 to compete against THF for hydrogen (H) bonding with water (). Here, we explore the alternative hypothesis that sorbitan esters interact with THF, competing against water. In our previous studies, we investigated H bonding between water and THF by probing the OH stretch band of water (; ). Water species coordinate differently with other molecules in solution. In particular, they can donate one (single donors, SD) or two (double donors, DD) hydrogens, and accept one (single acceptor, SA) or two (double acceptor, DA) hydrogens (; ). Therefore, water species include the following: SD-DA (coordination = 3, at ∼2,850–3,200 cm−1), SD-SA (coordination = 2, at ∼3,200–3,300 cm−1), DD-SA (coordination = 3, at ∼3,400–3,550 cm−1) and DD-DA (coordination = 4, at ∼3,300–3,400 cm-1). In some instances, water species are intermediate between these four main species, as indicated by their peak position. The relative abundance of each water species can be estimated from the area of each of the peaks representing them. In our previous research, we found that in pure water SD-SA and DD-DA are the most abundant species, and have similar abundance (; ). THF alters the distribution of water species, because it H bonds water and mainly accepts hydrogens (). Therefore, DD-SA dominate at high THF percentages (above 70%) (). Instead, DD-DA is dominant up to 40% THF. At intermediate THF percentages (50%–70% THF), the dominant water species is intermediate between DD-DA and SD-SA. The COC peak of THF at ≈1,065 cm−1 can also reveal H bonding between THF and water. Upon mixing with water, the COC peak of THF shifts and splits, as seen in Figure 2. Our recent study compared the OH stretch band of water and the COC band of THF, in THF-water mixtures (). This comparison showed that THF species that interact with DD-SA water species correspond to the 1,070 cm−1 peak in the COC band. The peak at 1,050 cm−1 dominates with >30% THF (v/v, relative to water), and it is related to THF species that mainly accept a single H. These water species correspond to SD-SA below 30% THF, and to species intermediate between SD-SA and DD-DA above 30% THF. Finally, the peak at 1,040 cm−1 corresponds to THF species that interact with SD-DA and DD-DA water species. In light of our previous study, we now use the COC band to study H bonding between THF and water, and between THF and sorbitan esters.
FIGURE 2
In particular, we analyze the COC band of THF to probe the hypothesis that sorbitan esters H bond THF. We also compare how H bonding of THF with water and with sorbitan esters affects the COC band of THF. With pure THF, we observe a dominant COC stretch peak at ≈1,070 cm−1 (
Sorbitan esters also interact with THF through H bonds, as indicated by splitting of the COC band in THF-sorbitan ester mixtures. We propose that H bonding occurs between the oxygen on the ether ring of THF and the OH groups of the sorbitan esters. Qualitatively, we observe that this interaction results in a shift of the OH stretch band of sorbitan esters (Supplementary Figure S1.4). Figure 2 shows the COC band of THF, which can be deconvolved into three peaks at ≈1,050 cm−1 (glassy THF), at ≈1,070 cm−1 (crystalline THF) and at ≈1,030 cm−1. Figure 3 shows the percent area under each peak, with the different sorbitan esters. These areas represent the relative abundance of the different THF types, e.g., crystalline vs. glassy. Dissimilar to observations with water, in mixtures of THF and sorbitan esters the peak at ≈1,050 cm−1 (glassy THF) is less prominent than the peak at the peak at 1,070 cm−1 (crystalline or ordered THF). This result suggests that THF is more ordered upon mixing with sorbitan esters compared to water, thereby leading to a more crystalline-like structure. THF ordering depends on the length and number of the sorbitan esters tails. THF is least ordered with Span 20, with the highest HLB number and the greatest 1,050–1,070 cm−1 ratio. Ordering increases when increasing the length of a single saturated tail from eleven carbons (for Span 20) to seventeen carbons (for Span 60), as seen in Figure 2A. Also, THF is less ordered with Span 80 (one unsaturated tail with 17 carbons) than with Span 85 (three unsaturated tails with 17 carbons), as seen in Figure 2B. While Span 65 was not analyzed in detail and is not shown in Figures 2, 3, its effect on the COC band is qualitatively similar to Span 85, as it also has a low HLB number (Supplementary Figure S1.5). Analyzing in detail the effect of Span 65 on the COC band of THF will be the objective of future research.
FIGURE 3

Percent areas under the COC band. These areas were obtained from the deconvolution of the dry reference subtracted COC band located at 1,030 (black), 1,055 (red), and 1,070 (green) cm−1 for (A) THF + Span 20 (squares) and THF + Span 60 (circles) and (B) THF + Span 40 (triangles), THF + Span 80 (diamonds) and THF + Span 85 (stars).
Moreover, in pure THF, we observe a peak at 905 cm−1, which corresponds to δ(C-C-O) (
FIGURE 4

δ(C-C-O) band of THF upon mixing with water and with the different sorbitan esters. The colors match the concentrations indicated in Figure 2.
The δ(C-C-O) band of THF is also affected by mixing THF with sorbitan esters (cf. Figure 4). In this case, the peak at 890 cm−1 (glassy THF) appears as a shoulder while the dominant peak is at 910 cm−1 (crystalline THF), in contrast with THF-water mixtures. This effect is most marked with the most hydrophobic sorbitan esters, Span 65 and Span 85, as qualitatively seen in Figure 4. We also quantitatively analyze the overall position of the δ(C-C-O) peak, rather than the two peaks convolved under it (Figure 5). Figure 5 shows that the shift of the peak to higher wavenumbers is greatest with the most hydrophobic sorbitan ester (e.g., Span 85), in line with the greater ordering of THF, whereas the shift is smallest with the most hydrophilic sorbitan ester (Span 20). Note that while we have not analyzed Span 65 in detail, its effect on the δ(C-C-O) is qualitatively similar to Span 85 (Figure 5), with a similar HLB number. In summary, the analysis of the δ(C-C-O) confirms the hypothesis that THF interacts with sorbitan esters through H bonds. It also indicates that THF is more ordered upon mixing with sorbitan esters than with water, in agreement with the analysis of the COC stretch band of THF.
FIGURE 5

Center position of δ(C-C-O) peak upon mixing different mole percentages of THF with different sorbitan esters.
So far, we have discussed H bonding between sorbitan esters and THF. We also highlighted that THF ordering upon mixing with sorbitan esters is correlated to their HLB numbers. Now, we compare the interactions between the different sorbitan ester tails and THF. Recall that sorbitan esters with saturated tails flocculate at low THF percentages, whereas this does not occur when tails are unsaturated. This result suggests that the characteristics of the tails influence THF-sorbitan ester interactions. To probe this aspect, we examine the CH band in the spectrum of sorbitan esters (Figure 6). In our previous study, we analyzed the CH bands of THF, to study THF interactions with water (
FIGURE 6

Comparison between the νas(CH2) band of different sorbitan esters, with varying percentages of THF. Black: neat THF; grey: 66 mole% sorbitan ester; red: 50 mole% sorbitan ester; green: 27 mole % sorbitan ester; blue: 17 mole% sorbitan ester; pink: neat sorbitan esters.
In particular, we use SAXS to compare the self-assembly of sorbitan esters in pure THF and with low (2%) water percentages relative to THF. While higher water percentages lead to emulsification of larger droplets or bulk separation, at low water percentages we can observe the onset of THF-water separation. In pure THF, we propose that sorbitan esters can interact with the solvent through either their heads or tails, leading to their self-assembly into micelles with their heads pointing outwards or inwards. This hypothesis is confirmed by SAXS data (Figures 7, 8). The SAXS patterns of all sorbitan esters in pure THF display broad peaks, which are characteristic of micelles with their heads pointing outwards or inwards (
TABLE 3
| THF | ||||
|---|---|---|---|---|
| D | Fractal type | Persistence length 1 | Persistence length 2 | |
| Span 20 | 2.13 | Surface | 25.1 | 11.9 |
| Span 40 | 2.78 | Mass | 25.5 | 13.7 |
| Span 60 | N/A | |||
| Span 65 | 2.41 | Mass | 27.2 | 12.1 |
| Span 80 | 2.38 | Surface | 26.6 | 15.9 |
| Span 85 | 2.38 | Surface | 27.4 | 21.6 |
| THF +2% water | ||||
| D | Fractal type | Persistence length 1 | Persistence length 2 | |
| Span 20 | 2.75 | Surface | 27.4 | 19.3 |
| Span 40 | 2.18 | Surface | 27.2 | 24.3 |
| Span 40, replicate | 2.21 | Surface | 27.3 | 23.9 |
| Span 60 | 2.46 | Surface | 26.9 | 23.3 |
| Span 65 | 2.27 | Surface | 23.9 | 23.2 |
| Span 80 | 2.72 | Surface | 23.2 | 24.8 |
| Span 85 | 2.15 | Mass | 24.5 | 31.6 |
| Span 85, replicate | 2.17 | Mass | 24.4 | 32.1 |
D and fractal type of sorbitan ester self assembled structures, for different sorbitan esters.
FIGURE 7

XRD patterns in the SAXS region for different sorbitan esters in either THF (A) or THF-DI water mixtures (B) containing 98% THF (v/v, relative to water). In (A), Span40 THF and Span40 THF B represent replicate experiments. Similarly, in (B), Span40 THF DI and Span40 THF DI Bj represent replicate experiments. These replicates are shown to allow the readers to compare replicate experiments. The comparison of the XRD of Span 85 in THF vs. THF and DI water is shown separately in Figure 8.
FIGURE 8

XRD patterns in the SAXS region for Span 85 in either pure THF or THF and DI water. In the figure legend, Span85THFDI and Span85THFDIJ represent replicate experiments, shown to highlight the reproducibility of the data. In (A) the peak at 42.7 Å has been removed for the purpose of the fit, while in (B) the peak is shown.
In summary, IR data suggest that the longer, saturated tail of some sorbitan esters (e.g., Span 40 and Span 60) interact with THF more strongly than the shorter, saturated tail of Span 20, or the unsaturated tails of Span 80 and Span 85. SAXS data complement these findings. They show that interactions between all sorbitan esters and water play a role in THF-water separation, beyond differences in the way they interact with THF. Indeed, in all cases 2% water enhances ordering, since it favors the formation of micelles with water in their interior.
4 Conclusion
Bottle tests, optical microscopy and FTIR spectromicroscopy show that sorbitan esters (Span 20, Span 40, Span 60, Span 80, Span 65 and Span 85) separate THF and water, which would otherwise be miscible due to H bonding. Separation occurs with up to 70% THF (v/v, relative to water), at THF to water ratios dependent on the sorbitan ester structure. Span 80 and Span 85 have double bonds in their tails. Double bonds in the tail (i.e., unsaturated tails) promote separation at the lowest THF percentages. In contrast, Span 40 and Span 60 (with saturated tails) flocculate at the lowest THF percentages (10%–30% THF) and do not effectively separate these THF-water mixtures. This suggests weaker interactions with water, and hence lower solubility at high water percentages. Saturated, but shorter tails (e.g., Span 20) still allow separation at low THF percentages.
We use ATR-FTIR to study separation mechanisms between THF and water by sorbitan esters. ATR-FTIR reveals that mixing with water disorders THF, resulting in an increase of the percentage of glassy to crystalline THF.
Sorbitan esters impede THF-water interactions. Similar to water, they H bond with THF. However, dissimilar to water, they order THF. Ordering increases with tail length (e.g., from Span 20 to Span 60) and with increasing number of tails (e.g., from Span 80 to Span 85). This is reflected in the COC and δ(C-C-O) band of THF.
Furthermore, we observe a blueshift in the νas(CH2) band of sorbitan esters upon mixing with increasing THF percentages, because of the formation of C-H…O H bonds. The blueshift is most marked for Span 40 and 60. This result indicates that Span 40 and Span 60 interact with THF more strongly than Span 20, Span 80 and Span 85. In contrast, Span 40 and Span 60 interact with water less strongly than Span 20, Span 80 and Span 85, as indicated by their flocculation at the lowest THF percentages (10%–30% THF, v/v relative to water).
Based on these results, we propose that separation between THF and water occurs due to two primary mechanisms, depending on the sorbitan ester characteristics. Specifically, separation occurs primarily (but not exclusively) because the hydrophilic head of Span 20, Span 80 and Span 85 compete for interactions with water against THF. In contrast, separation occurs primarily (but not exclusively) due to preferential interactions between the tail of Span 40 and Span 60 and THF. Note that we are not contending that these are the only mechanisms involved. Indeed, THF interacts with the tails of all sorbitan esters, and water interacts with their heads.
The role of water-sorbitan ester interactions on THF-water separation is highlighted by SAXS data. In THF, sorbitan esters self-assemble into micelles with their heads pointing outwards or inwards, because both the heads and the tails of sorbitan esters interact with THF. Both types of micelles further assemble to form either closed-packed clusters (surface fractals, in the case of Span 20, Span 80, and Span 85), or open-packed clusters (mass fractals, in the case of Span 40, Span 60 and Span 65). In THF, we identify two different persistence lengths, because of the limited order of sorbitan ester self-assembled structures. With 2% water, we observe a transition from mass to surface fractals in all cases except Span 85, possibly due to its incomplete hydration. Under these conditions, sorbitan esters self assemble into micelles with their heads pointing inwards, which host water in their interior. We identify a single persistence length (≈25 Å), which corresponds to the distance between the micelle centers. This shows enhanced ordering.
Our future research will focus on further elucidating how separation mechanisms are affected by the tail structure, for different amphiphiles, including ones the head of which does not interact with THF, as well as others with varying heads that H bond THF. This research is part of our overall goal of developing a predictive toolbox to optimize solvent separation from water.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
BB: Data curation, Formal Analysis, Investigation, Validation, Visualization, Writing–original draft, Writing–review and editing. AM: Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. TL: Data curation, Formal Analysis, Software, Visualization, Writing–review and editing. EP: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. Part of the research described in this paper was performed at the Canadian Light Source, a national research facility of the University of Saskatchewan, which is supported by the Canada Foundation for Innovation (CFI), the Natural Sciences and Engineering Research Council (NSERC), the National Research Council (NRC), the Canadian Institutes of Health Research (CIHR), the Government of Saskatchewan, and the University of Saskatchewan. The authors acknowledge the support of the Natural Sciences and Engineering Research Council of Canada (provided through an NSERC Discovery grant, awarded to EP, PIN 537871).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frsfm.2023.1329058/full#supplementary-material
References
1
AliagaC.TsungC. K.AlayogluS.KomvopoulosK.YangP.SomorjaiG. A. (2011). Sum frequency generation vibrational spectroscopy and kinetic study of 2-methylfuran and 2, 5-dimethylfuran hydrogenation over 7 nm platinum cubic nanoparticles. J. Phys. Chem. C115, 8104–8109. 10.1021/jp111343j
2
BarmanS. R.BanerjeeP.MukhopadhayayA.DasP. (2017). Biodegradation of acenapthene and naphthalene by Pseudomonas mendocina: process optimization, and toxicity evaluation. J. Environ. Chem. Eng.5, 4803–4812. 10.1016/j.jece.2017.09.012
3
BartokovaB.LaredoT.MarangoniA. G.PensiniE. (2023a). Mechanism of tetrahydrofuran separation from water by stearic acid. J. Mol. Liq.391, 123262. 10.1016/j.molliq.2023.123262
4
BartokovaB.MarangoniA. G.LaredoT.PensiniE. (2023b). Role of hydrogen bonding on solvent separation using amphiphilic sorbitan ester. Colloids and Surfaces C. 1, 100004.
5
BartokovaB.MarangoniA. G.LaredoT.StobbsJ.MeszarosP.PensiniE. (2023c). Effect of hydrogen bonding on the mixing behaviour of ternary aqueous mixtures. J. Mol. Liq.283, 122124. 10.1016/j.molliq.2023.122124
6
BelhajA.DesnouesN.ElmerichC. (2002). Alkane biodegradation in Pseudomonas aeruginosa strains isolated from a polluted zone: identification of alkB and alkB-related genes. Res. Microbiol.153, 339–344. 10.1016/s0923-2508(02)01333-5
7
BeshaA. T.BekeleD. N.NaiduR.ChadalavadaS. (2018). Recent advances in surfactant-enhanced In-Situ Chemical Oxidation for the remediation of non-aqueous phase liquid contaminated soils and aquifers. Environ. Technol. Innovation9, 303–322. 10.1016/j.eti.2017.08.004
8
BoonchanS.BritzM. L.StanleyG. A. (1998). Surfactant‐enhanced biodegradation of high molecular weight polycyclic aromatic hydrocarbons by Stenotrophomonas maltophilia. Biotechnol. Bioeng.59, 482–494. 10.1002/(sici)1097-0290(19980820)59:4<482::aid-bit11>3.0.co;2-c
9
de Brito CardosoG.Nascimento SouzaI.MourãoT.FreireM. G.SoaresC. M. F.Silva LimaÁ. (2014). Novel aqueous two-phase systems composed of acetonitrile and polyols: phase diagrams and extractive performance. Sep. Purif. Technol.124, 54–60. 10.1016/j.seppur.2014.01.004
10
DhamoleP. B.MahajanP.FengH. (2010). Phase separation conditions for sugaring-out in Acetonitrile−Water systems. J. Chem. Eng. Data55 (9), 3803–3806. 10.1021/je1003115
11
EarndenL.MarangoniA. G.LaredoT.StobbsJ.MarshallT.PensiniE. (2022a). Decontamination of water co-polluted by copper, toluene and tetrahydrofuran using lauric acid. Sci. Rep.12, 15832–15920. 10.1038/s41598-022-20241-4
12
EarndenL.MarangoniA. G.LaredoT.StobbsJ.PensiniE. (2022b). Mechanisms of separation between tetrahydrofuran and water using hydroxystearic acid. Phys. Fluids34, 097119. 10.1063/5.0108008
13
EarndenL.MarangoniA. G.LaredoT.StobbsJ.PensiniE. (2022c). Self-Assembled glycerol monooleate demixes miscible liquids through selective hydrogen bonding to water. J. Mol. Liq.367, 120551. 10.1016/j.molliq.2022.120551
14
FleyfelF.DevlinJ. P. (1991). Carbon dioxide clathrate hydrate epitaxial growth: spectroscopic evidence for formation of the simple type-II carbon dioxide hydrate. J. Phys. Chem.95, 3811–3815. 10.1021/j100162a068
15
HallJ.SooleK.BenthamR. (2011). Hydrocarbon phytoremediation in the family Fabacea—a review. Int. J. phytoremediation13, 317–332. 10.1080/15226514.2010.495143
16
HarthcockM. A. (1989). Probing the complex hydrogen bonding structure of urethane block copolymers and various acid containing copolymers using infra-red spectroscopy. Polymer30, 1234–1242. 10.1016/0032-3861(89)90041-4
17
HuD.LiX.ChenZ.CuiY.GuF.JiaF.et al (2018). Performance and extracellular polymers substance analysis of a pilot scale anaerobic membrane bioreactor for treating tetrahydrofuran pharmaceutical wastewater at different HRTs. J. Hazard. Mater.342, 383–391. 10.1016/j.jhazmat.2017.08.028
18
IsaacsonC.MohrT. K.FieldJ. A. (2006). Quantitative determination of 1, 4-dioxane and tetrahydrofuran in groundwater by solid phase extraction GC/MS/MS. Environ. Sci. Technol.40, 7305–7311. 10.1021/es0615270
19
JessopP.CunninghamM. (20206). CO2-switchable materials. London, United Kingdom: Royal Society of Chemistry.
20
KolehmainenE.TurunenI. (2007). Micro-scale liquid–liquid separation in a plate-type coalescer. Chem. Eng. Process. Process Intensif.46, 834–839. 10.1016/j.cep.2007.05.027
21
LeeJ.ChoJ.KimD. M.ParkS. (2011). Separation of tetrahydrofuran and water using pressure swing distillation: modeling and optimization. Korean J. Chem. Eng.28, 591–596. 10.1007/s11814-010-0467-1
22
LeontowichA. F.GomezA.Diaz MorenoB.MuirD.SpasyukD.KingG.et al (2021). The lower energy diffraction and scattering side-bounce beamline for materials science at the Canadian Light Source. J. Synchrotron Radiat.28, 961–969. 10.1107/s1600577521002496
23
LevenbergK. (1944). A method for the solution of certain non-linear problems in least squares. Q. Appl. Math.2, 164–168. 10.1090/qam/10666
24
LiH.StanwixP.AmanZ.JohnsM.MayE.WangL. (2016). Raman spectroscopic studies of clathrate hydrate formation in the presence of hydrophobized particles. J. Phys. Chem. A120, 417–424. 10.1021/acs.jpca.5b11247
25
LiuL.XiaY.WangL.ShiR.YanS.ZhaoX.et al (2022). Cyanate ester resin with high heat‐resistance and degradable diacetal structure: synthesis, polymerization, and properties. Macromol. Mater. Eng.307, 2200423. 10.1002/mame.202200423
26
LongerasO.GautierA.Ballerat-BusserollesK.AndansonJ. M. (2020). Deep eutectic solvent with thermo-switchable hydrophobicity. ACS Sustain. Chem. Eng.8, 12516–12520. 10.1021/acssuschemeng.0c03478
27
MalakahmadA.HoL. L. H. (2017). UV/H2O2 oxidation process optimization by response surface methodology for removal of polycyclic aromatic hydrocarbons (PAHs) from water. Desalination Water Treat.65, 408–417. 10.5004/dwt.2017.20278
28
MammoneJ. F.SharmaS. K.NicolM. (1980). Raman spectra of methanol and ethanol at pressures up to 100 kbar. J. Phys. Chem.84, 3130–3134. 10.1021/j100460a032
29
MarquardtD. W. (1963). An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Industrial Appl. Math.11, 431–441. 10.1137/0111030
30
MarshallT.EarndenL.MarangoniA. G.LaredoT.PensiniE. (2022). Cubic mesophases of self-assembled amphiphiles separate miscible solvents. Colloids Surfaces A650, 129548. 10.1016/j.colsurfa.2022.129548
31
MizunoK.ImafujiS.FujiwaraT.OhtaT.TamiyaY. (2003). Hydration of the CH groups in 1,4-dioxane probed by NMR and IR: contribution of blue-shifting CH···OH 2 hydrogen bonds. J. Phys. Chem. B107, 3972–3978. 10.1021/jp021712+
32
PatelV.MarangoniA. G.Mirzaee GhazaniS.LaredoT.StobbsJ.PensiniP. (2023). Effect of bacterial surfactants on the phase behavior of miscible pollutants in water. Colloids Surf. C1, 100013. 10.1016/j.colsuc.2023.100013
33
PrasadP. S. R.PrasadK. S.ThakurN. K. (2007). Laser Raman spectroscopy of THF clathrate hydrate in the temperature range 90–300 K. Spectrochim. Acta Part A68, 1096–1100. 10.1016/j.saa.2007.06.049
34
SalesC. M.GrosternA.ParalesJ. V.ParalesR. E.Alvarez-CohenL. (2013). Oxidation of the cyclic ethers 1, 4-dioxane and tetrahydrofuran by a monooxygenase in two Pseudonocardia species. Appl. Environ. Microbiol.79, 7702–7708. 10.1128/aem.02418-13
35
SionkowskaA.PłaneckaA. (2013). Preparation and characterization of silk fibroin/chitosan composite sponges for tissue engineering. J. Mol. Liq.178, 5–14. 10.1016/j.molliq.2012.10.042
36
SouzaR. L.LimaR. A.CoutinhoJ. A.SoaresC. M.LimaÁ. S. (2015). Aqueous two-phase systems based on cholinium salts and tetrahydrofuran and their use for lipase purification. Sep. Purif. Technol.155, 118–126. 10.1016/j.seppur.2015.05.021
37
SunY.ZhouQ.XuY.WangL.LiangX. (2011). Phytoremediation for co-contaminated soils of benzo [a] pyrene (B [a] P) and heavy metals using ornamental plant Tagetes patula. J. Hazard. Mater.186, 2075–2082. 10.1016/j.jhazmat.2010.12.116
38
TabataM.KumamotoM.NishimotoJ. (1994). Chemical properties of water-miscible solvents separated by salting-out and their application to solvent extraction. Anal. Sci.10, 383–388. 10.2116/analsci.10.383
39
TobyB. H.Von DreeleR. B. (2013). GSAS-II: the genesis of a modern open-source all purpose crystallography software package. J. Appl. Crystallogr.46, 544–549. 10.1107/s0021889813003531
40
UrtiagaA.GómezP.ArrutiA.OrtizI. (2014). Electrochemical removal of tetrahydrofuran from industrial wastewaters: anode selection and process scale‐up. J. Chem. Technol. Biotechnol.89, 1243–1250. 10.1002/jctb.4384
41
Von DreeleR. B. (2014). Small-angle scattering data analysis in GSAS-II. J. Appl. Crystallogr.47, 1784–1789. 10.1107/s1600576714018366
42
WangB.EzejiasT.FengH.BlaschekH. (2008). Sugaring-out: a novel phase separation and extraction system. Chem. Eng. Sci.63, 2595–2600. 10.1016/j.ces.2008.02.004
43
WangW. H.HoagG. E.CollinsJ. B.NaiduR. (2013). Evaluation of surfactant-enhanced in situ chemical oxidation (S-isco) in contaminated soil. Soil Pollut.224, 1713. 10.1007/s11270-013-1713-z
44
YinY.YangY.de Lourdes MendozaM.ZhaiS.FengW.WangY.et al (2017). Progressive freezing and suspension crystallization methods for tetrahydrofuran recovery from Grignard reagent wastewater. J. Clean. Prod.144, 180–186. 10.1016/j.jclepro.2017.01.012
45
YoshikawaM.MasakiK.IshikawaM. (2002). Pervaporation separation of aqueous organic mixtures through agarose membranes. J. Membr. Sci.205, 293–300. 10.1016/s0376-7388(02)00131-x
Summary
Keywords
mixing behaviour, separation, sorbitan esters, hydrogen bonding, solvent structure
Citation
Bartokova B, Marangoni AG, Laredo T and Pensini E (2023) Effect of sorbitan ester structure on the separation between tetrahydrofuran and water. Front. Soft Matter 3:1329058. doi: 10.3389/frsfm.2023.1329058
Received
27 October 2023
Accepted
29 November 2023
Published
18 December 2023
Volume
3 - 2023
Edited by
Mohd Sajid Ali, King Saud University, Saudi Arabia
Reviewed by
Animesh Pan, University of Rhode Island, United States
Ajaya Bhattarai, Tribhuvan University, Nepal
Updates

Check for updates
Copyright
© 2023 Bartokova, Marangoni, Laredo and Pensini.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Erica Pensini, epensini@uoguelph.ca
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.