Subscriber access provided by La Trobe University Library
Article
Removal of thiophene from mixtures with n-heptane by selective extraction using deep eutectic solvents Mohamed Kamel Hadj-Kali, Sarwono Mulyono, Hanee Farzana Hizaddin, Irfan Wazeer, Lahssen El-blidi, Emad Ali, Mohd Ali Hashim, and Inas Muen AlNashef Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.6b01654 • Publication Date (Web): 07 Jul 2016 Downloaded from http://pubs.acs.org on July 9, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Industrial & Engineering Chemistry Research is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Industrial & Engineering Chemistry Research
1
Removal of thiophene from mixtures with n-heptane by
2
selective extraction using deep eutectic solvents
3 4
Mohamed K. Hadj-Kali*a, Sarwono Mulyonoa, Hanee F. Hizaddinb, Irfan Wazeera, Lahssen El-
5
blidia, Emad Alia, Mohd. Ali Hashimb, Inas M. AlNashefc a
6 7
Chemical Engineering Department, College of Engineering King Saud University, P.O.Box 800, Riyadh, 11421, Saudi Arabia.
8 9
b
University of Malaya Center for Ionic Liquids (UMCiL), Department of Chemical Engineering, Faculty of Engineering, University of Malaya, 50603 Kuala Lumpur, Malaysia.
10 11
c
12
*Corresponding Author. Email address:
[email protected] Department of Chemical and Environmental Engineering, Masdar Institute of Science and Technology, Abu Dhabi, United Arab Emirates.
13 14
Abstract
15
This work investigates the use of deep eutectic solvents (DESs) to extract sulfur-based
16
compound from n-heptane as model diesel compound. Four DESs were prepared by combining
17
tetrabutylammonium bromide or methyltriphenylphosphonium bromide with ethylene glycol,
18
triethylene glycol or sulfolane. All DESs showed good ability to extract thiophene with the best
19
extraction efficiency, 35%, for the sulfolane-based DES. The extraction efficiency can be further
20
enhanced to reach 98% when five extraction cycles are performed. Moreover, the DESs were
21
easily regenerated using rotary evaporation. In addition,1H NMR analysis is used to elucidate the
22
extraction mechanism. Finally, the COSMO-RS model was used to predict the ternary tie lines
23
for the studied systems and the NRTL model allowed correlating the experimental data with an
24
average root mean square deviation lower than 2% for all DESs. These models can be utilized
25
for further simulation analysis of the extraction process.
26 27
Keywords: Desulfurization, Liquid-liquid Extraction, Deep Eutectic Solvents, COSMO-RS,
28
NRTL. 1 ACS Paragon Plus Environment
Industrial & Engineering Chemistry Research
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
29
Page 2 of 33
1. Introduction
30
The rapid growth of the world’s population accompanied by the increase in the standard of
31
living over the last several decades has caused environmental pollution by the emission of SOx
32
to continue increasing. Stringent legislation has been implemented all over the world to regulate
33
the sulfur content of transportation fuels. Hydrodesulfurization is the conventional industrial
34
process for removing sulfur containing compounds from liquid fuels.1 However, this process
35
operates under severe conditions with high pressure and temperature2 and is not efficient for
36
deep removal of refractory heterocyclic sulfur compounds such as thiophene, benzothiophene,
37
dibenzothiophene and its derivatives.3-6 Thus, developments of innovative complementary
38
technologies for fuel oil processing such as extraction, adsorption, oxidation and bio-
39
desulfurization have been sought during the last few decades. Extractive desulfurization is one of
40
the most promising desulfurization method because it operates under mild conditions and have
41
simple separation mechanism.7 Volatile organic compounds such as sulfolane, N-methyl-2-
42
pyrrolidone and dimethylformamide have been used as conventional solvents in the extraction
43
process in the last few decades.8 Hence, research efforts are focused on finding the most suitable
44
solvent for use in deep desulfurization process which has less impact on the environment.
45
Recently, the use of ionic liquids (ILs) to replace conventional solvents have been reported in
46
various research work. ILs are gaining increasing attention owing to their favorable chemical and
47
physical properties5 including high thermal stability, non-volatility and recyclability.9-13 A
48
comprehensive review on the application of ILs on desulfurization of fuel oils can be referred in
49
a review done by Francisco et al.14
50
In the past decade, ILs analogues called deep eutectic solvents (DESs) have emerged as
51
alternatives to conventional solvents and ILs.15,16 A DES is a mixture of two or more components
52
that has a melting point less than that of any of its individual constituent. DESs are usually 2 ACS Paragon Plus Environment
Page 3 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Industrial & Engineering Chemistry Research
53
composed of salts combined with a hydrogen-bond donor (HBD) or a complexing agent. DESs
54
have many advantages over ILs as they can be easily synthesized with high purity and they are
55
generally cheaper than ILs due to the simple synthesis process and their cost depend on the cost
56
of the constituent components. By carefully choosing the components of the DES, it is possible
57
to obtain non-toxic and biodegradable solvents. In addition, most DESs do not react with water.
58
DESs have been used extensively in various applications including material preparation,17
59
electrochemistry,18 purification of palm oil,19 enzyme catalysis20 and separation processes.21
60
DESs have good dissolution capabilities for different substances because they can form hydrogen
61
bond by receiving or giving electrons or protons. However, reports on the applications of DESs
62
in separation processes are not as widespread as those using ILs. Among the earliest report on
63
the use of DESs in separation processes is by Abbott et al. who used a mixture of quaternary
64
ammonium salts with glycerol to extract excess glycerol from soybean oil-based biodiesel.21
65
Later, Guo et al. used DESs to extract phenol from oil22 and achieved high phenol extraction
66
efficiencies by using quaternary ammonium salts based DESs which could reach up to 99.9%.
67
They also reported that the extraction efficiency was not sensitive to temperature. Shahbaz et al.
68
investigated eighteen DESs for the removal of residual catalyst content from palm oil-based
69
biodiesel23 and reported that the removal efficiency of the residual catalyst increased with the
70
increase in DES:biodiesel molar ratio.
71
More recently, sixteen DESs were experimentally screened by Li et al. to extract
72
benzothiophene from n-octane for use in desulfurization of liquid fuels.1 They stated that choline
73
chloride/propionate (1:2), tetrabutyl ammonium chloride/propionate (1:2) and tetrabutyl
74
ammonium chloride / polyethylene glycol (1:2) DESs reported the highest extraction efficiency
75
and that the extraction of benzothiophene from n-octane can be conducted up to five times
3 ACS Paragon Plus Environment
Industrial & Engineering Chemistry Research
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 33
76
without significantly compromising the extraction efficiency. However, after six times of
77
repeated use, the DESs became saturated and lost their extraction capacity.
78
Previously, our group has reported liquid-liquid equilibrium data on the use of imidazolium
79
and pyridinium-based ILs and ammonium- and phosphonium-based DESs to remove nitrogen
80
compounds from n-hexadecane as model diesel fuel.24,25 The DESs used in our previous work
81
are selected based on COSMO-RS screening and it was hypothesized that the DESs selected are
82
highly potential for use in removal of aromatic sulfur compounds as well.26 This work, thus, aims
83
at studying the effectiveness of using similar DESs to remove thiophene from liquid fuel. We
84
showed that low cost DESs can be used for the removal of sulfur containing compounds using
85
liquid-liquid extraction method. Four DESs have been tested where ethylene glycol (EG),
86
triethylene glycol (TrEG) and sulfolane (Sul) were used as complexing agents, while
87
tetrabutylammonium bromide (TBAB) and methyltriphenylphosphonium bromide (MeTPPBr)
88
were chosen as typical salts. The extraction effectiveness was assessed experimentally and also
89
numerically using two models: COSMO-RS and NRTL. Furthermore, the extraction mechanism
90
with the best DES was investigated using 1H NMR analysis.
91
2. Experimental Details
92
2.1.Chemicals and experimental protocol
93
All chemicals used in this work are listed in Table S1 in the Supporting Information and their
94
chemical structures are listed in Table 1. They were of high purity (>98 wt.%) and were used
95
without any further purification. The DESs were prepared according to the method described by
96
Abbott et al.16 The different combinations tested are listed in Table S2 in Supporting Information
97
along with their specific molar ratio. The mixtures were put in screw-capped bottles and then
98
stirred in an incubating-shaker equipped with temperature and speed control at a temperature of 4 ACS Paragon Plus Environment
Page 5 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Industrial & Engineering Chemistry Research
99
100°C (±0.1°C) with a rotational speed of 200 rpm until a clear liquid was formed. For the
100
mixture TBAB/sulfolane, the minimum molar ratio giving a liquid mixture at room temperature
101
was 1:7. However, mixing MeTPPBr with either sulfolane or TrEG led to a solid mixture at
102
room temperature even for higher molar ratios.
103
The thermal stability of three TBAB-based DESs investigated in this work has been checked
104
using TGA-50 thermogravimetric analyzer (Shimadzu, Kyoto, Japan) under air atmosphere with
105
flux of 100 mL/min and a heating rate of 10 °C/min. approximately 3 mg of each DES sample
106
were subjected to a temperature range of 25 to 400 °C. The analysis results are shown in fig. S1
107
(Supporting Information). TGA analysis confirms that all three DESs are thermally stable and
108
have similar behavior to their individual constituents. Moreover, it can be concluded from this
109
analysis that the hydrogen bond interactions are stronger between TBAB and EG compared with
110
TrEG which is also stronger than that between TBAB and sulfolane.
111
On the other hand, feed mixture containing 10 wt% of thiophene in n-heptane was prepared
112
by mixing weighed amounts of the chemicals using an analytical balance (±0.0001 g).The feed
113
was then mixed with the DESs in a mass ratio of 1:1. Each set of experiments was conducted at
114
25°C. The vials were placed in the incubator shaker. The shaking time was six hours followed by
115
a settling time of about 12 hours to guarantee that the equilibrium state was completely attained.
116
Samples were then taken from the top and bottom layers and analyzed using HPLC.
117
2.2. Analysis
118
Samples from the top and bottom layers were withdrawn using a syringe and then diluted using
119
acetonitrile. The samples were analyzed using a HPLC Agilent 1100 series (USA) with a zorbax
120
eclipse xdb-c8 column (Agilent, USA). The temperature of the column oven was set to 30°C.
121
The mobile phase was acetonitrile and distilled water with a volume ratio of 3:1. The flow rate of
122
the mobile phase was 1.4 ml min-1 with a pressure of 120 bars. Each sample was analyzed three 5 ACS Paragon Plus Environment
Industrial & Engineering Chemistry Research
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 33
123
times and the average is reported. The maximum uncertainty is estimated to be ±0.008 in mole
124
fraction, while the average uncertainty doesn't exceed ±0.002. Details about uncertainty
125
estimation and error propagation, as well as the calibration of thiophene are shown in the
126
supporting information.
127
Moreover, in order to highlight the extraction mechanism, both nuclear magnetic resonance
128
spectroscopy (NMR) and FT-IR were used. 1H NMR spectra were recorded using a JEOL
129
RESONANCE spectrome-ter ECX-500 II. Dimethyl sulfoxide (DMSO) was used as solvent. The
130
infrared spectra were recorded on a IRTracer-100 FT-IR spectrometer interfaced (Shimadzu,
131
Kyoto, Japan). The spectrum was obtained over 64 scans with a spectrum resolution of 4 cm−1 in
132
the region 4000 to 400 cm−1.
133
3. Results and discussion
134
3.1. Consistency tests
135
The Othmer-Tobias27 and Hand28 correlations were used to perform the consistency tests of
136
the experimental results as shown by Eq. 1 and 2, respectively.
137
138
Thio Thio
= +
hep
= +
DES
(Eq. 1)
(Eq. 2)
139
Here xhep, xDES and xThio refer to the concentrations of n-heptane, DES and thiophene,
140
respectively. Fitting parameters of the Othmer-Tobias correlation are a and b, while c and d are
141
the fitting parameters of the Hand correlation, whereas ' and '' refers to the extract and raffinate
142
phases, respectively. The value of R2 close to unity (the linearity of each plot) shows the
143
excellent degree of consistency of the data. The parameters of the Othmer-Tobias and Hand
144
correlations are listed in Table 2.
6 ACS Paragon Plus Environment
Page 7 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Industrial & Engineering Chemistry Research
145
3.2. Distribution ratio and selectivity
146
Table S3 in the Supporting Information shows the experimental tie lines with the values of
147
distribution coefficient D, and selectivity S, for the four ternary systems (thiophene (1) + n-
148
heptane (2) + DES (3)) where the compositions are expressed in mole fractions. Figure 1 depicts
149
the corresponding triangular phase diagrams for each system. The values of distribution
150
coefficient obtained in this work are all less than unity, which indicates that a relatively large
151
quantity of solvent is required in order to achieve high separation effectiveness. Nevertheless, the
152
values of selectivity reported here are encouraging. Figures. 2 and 3 show the plot of distribution
153
coefficient and selectivity as a function of thiophene concentration in the raffinate phase for
154
different DESs. TBAB/Sul (1:7) gave the highest values of distribution coefficient followed by
155
TBAB/TrEG (1:4), MeTPPBr/EG (1:4) and TBAB/EG (1:4). This can be attributed to the higher
156
molar ratio of the HBD in TBAB/Sul (1:7) DES compared to other DESs. For selectivity, at low
157
concentrations the TBAB/TrEG (1:4) gave the highest values followed by TBAB/Sul (1:7),
158
TBAB/EG (1:4) and MeTPPBr/EG (1:4). However, as the concentration of thiophene in the
159
raffinate increases, the TBAB/TrEG (1:4) gave lower value of S compared to the other DESs.
160
Besides, for the best DES (TBAB/Sul (1:7)), the 1H NMR of the raffinate phase has shown that
161
no sulfolane is detected in the top layer (see fig. S2 in the supporting material).
162
Moreover, we have compared the selectivity and distribution ratio results obtained in this
163
work with those available in the literature for more than 20 ionic liquids and pure sulfolane. This
164
comparison is provided in table S4 as supporting material. In overall, it is obvious that many
165
ionic liquids show better extraction performances than all DESs investigated in this work.
166
However, if we consider the relative high price and the complicated synthesis procedures of
167
ionic liquids compared to DESs, we can conclude that DESs could be considered as potential
168
alternatives especially for industrial scale applications. When compared to pure sulfolane, we can 7 ACS Paragon Plus Environment
Industrial & Engineering Chemistry Research
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 33
169
see that the salt has almost no effect on the extraction performance of thiophene (selectivity and
170
distribution ratio values are almost the same). However, the addition of TBAB has other
171
advantages compared to the pure sulfolane, since (i) no cross contamination is observed between
172
the two layers with TBAB/Sul DES and (ii) the operating temperature is lower than in the
173
classical process. So, low energy consumption is required when using the DES.
174
3.3. Thiophene extraction efficiency
175
The following equation is used to calculate the extraction efficiency.
176
Extraction efficiency '%) =
*1 −*2 *1
× 100
(Eq. 3)
177
Where C1 is the initial concentration of thiophene in fuel and C2 is the final concentration of
178
thiophene in the raffinate phase after extraction with the DES. The thiophene extraction
179
efficiency of each DES is given in Table 3. As expected the highest extraction efficiency is for
180
TBAB/Sul (1:7). It can be noted that both the salt and HBD affect the extraction efficiency.
181
However, in most cases the difference is not significant. The results obtained in this work are
182
somewhat comparable with those reported in the literature for the same type of oil mixtures. For
183
example, Jian-long et al.29 studied the extraction of thiophene from thiophene, n-heptane, xylol
184
mixture using five different ionic liquids. The reported sulfur removal percentage, with IL to oil
185
mass ratio of 1:1, ranges from 21.8% to 45.5% depending on the type of IL used. This means that
186
the DES investigated in this work has slightly lower extraction ability (ranging from 16% to
187
35%). However, we used binary oil mixture while ternary oil mixture was utilized in the
188
aforementioned work.29 Alonso et al. showed that using [C8mim][BF4] IL, thiophene recovery
189
can reach 79 wt% when three cycles are used.30,31 However, they used different oil mixtures
190
which contain thiophene in i-octane, in n–hexane, or in toluene.
8 ACS Paragon Plus Environment
Page 9 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Industrial & Engineering Chemistry Research
191
3.4. Extraction mechanism
192
1
H NMR and FT-IR analysis were used to elucidate the mechanism of extracting thiophene
193
using DESs. Fig. 4 shows the 1H NMR analysis of pure TBAB, TBAB with sulfolane in DES
194
form, and different molar ratios of DES/thiophene. First, the comparison of free TBAB and DES
195
showed that the chemical shifts of all peaks of TBAB decreased, this shielding effect is due to
196
interactions between oxygen atoms of sulfolane and the ammonium cation in TBAB. These
197
interactions lead to the decrease of ammonium cation electron withdrawing effect on all TBAB
198
hydrogen atoms, especially atoms adjacent to it such as 1 (moves from 3.360 ppm to 2.966 ppm).
199
The addition of thiophene to DES causes an opposite effect, that is to say, all peaks of TBAB
200
return to values very close to that in pure TBAB. This shows that the bromide anion in TBAB
201
interacts with thiophene aromatic ring. Moreover, it can be noted that oxygen atoms in sulfolane
202
created new interactions with the sulfur atom in thiophene. The result of these interactions is the
203
decreasing of the electron withdrawing power of sulfur atom over hydrogen atoms in thiophene
204
ring and increasing the electronic density. These effects are reflected in the values of chemical
205
shifts of hydrogen atoms of thiophene and sulfolane. Fig. 4 shows that the hydrogen atoms in
206
thiophene moved to higher field intensities with increasing of thiophene:DES ratio from 0.4:1to
207
1:1. For example, H atom in position 1 (7.437 ppm, 7.420, 7.255, 7.242) and for position 2
208
(7.232 ppm, 7.188, 7.034, 7.020). On the other hand, for sulfolane, these interactions resulted in
209
the increasing of electron withdrawing power of sulfur atom, and its hydrogen atoms moved in
210
lower field intensities (from 1.923 ppm, 2.700 in DES/thiophene: 1/0.4 to 2.050 ppm, 2.860 in
211
DES/thiophene: 1/1 respectively). Therefore, the results of 1H NMR analysis showed the
212
existence of interactions between DES and thiophene, and it can be concluded that these
213
interactions contributed to the extraction efficiency.
9 ACS Paragon Plus Environment
Industrial & Engineering Chemistry Research
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 33
214
The FT-IR analysis confirms the extraction mechanism. Three samples have been analyzed:
215
(i) pure thiophene, (ii) thiophene – sulfolane mixture and (iii) thiophene – (TBAB – Sulfolane)
216
DES mixture, as shown in fig. 5. As shown in this figure, FT-IR spectra confirms the existence
217
of interactions between thiophene and pure sulfolane (the green line) and also between thiophene
218
and DES system (the blue line). The main modifications due to the addition of sulfolane alone or
219
the DES mixture are observed on the signals relative to the symmetric (1408 cm-1) and the anti-
220
symmetric (1590 cm-1) carbon double bond (C=C) stretching vibration bands of the thiophene
221
rings. These bands were moved to higher wavenumbers (from 1408 to 1415 and from 1590 to
222
1604 (with sulfolane) then to 1625 cm-1 (with the DES)) and became less stretched and stretched.
223
These effects could be attributed to the interactions between thiophene ring and oxygen atoms of
224
sulfolane and these interactions were conserved by the addition of the salt in the DES mixture.
225
3.5. Number of extraction stages
226
In order to determine the number of cycles (stages) necessary to extract all thiophene from
227
thiophene/n-heptane mixture using TBAB:Sulfolane DES with molar ratio 1:7, we have
228
separated the top layer (n-heptane rich phase) from the bottom layer (DES rich phase) and we
229
have added the same amount of DES to the raffinate phase. This procedure was repeated many
230
times and the thiophene concentration for each step was determined using HPLC analysis.
231
Fig. 6 shows the HPLC chromatograms of thiophene in the raffinate phase after each
232
extraction step. As can be seen, the peak representing thiophene decreased gradually until
233
reaching a total removal of about 98% after six extraction cycles. Therefore, after only five steps,
234
it is possible to reduce the sulfur content of the model fuel (n-heptane) to less than 0.02 of its
235
initial value.
10 ACS Paragon Plus Environment
Page 11 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Industrial & Engineering Chemistry Research
236
3.6. Deep eutectic solvent regeneration
237
The regeneration of the extraction solvent and its reuse has economic and environmental
238
benefits. For this purpose, the back extraction and rotary evaporator process were used. In our
239
case, the regeneration of TBAB/Sul (1:7) DES was performed by using a rotary evaporator under
240
100 mbar vacuum and 40°C. After the first extraction, the two phases were separated using a
241
separating funnel and thiophene was removed from the bottom layer (DES-rich phase) using a
242
rotary evaporator. Then, the recovered DES was used for new extraction cycle. As can be seen
243
from Fig. 7, even after four cycles, the performance of recycled DES was as effective as that of a
244
fresh one. In this figure, the extraction performance reflects the quantity of thiophene trapped
245
and is expressed in (%) based on the first extraction cycle. Therefore, these advantages, i.e. easy
246
recycling and maintained extraction efficiency, make this DES an attractive solvent for
247
desulfurization.
248
On the other hand, based on the normal boiling point data of thiophene and different
249
hydrogen bond donors (listed in table S5), it is clear that the same technique (rotary evaporation)
250
can also be used for all DESs investigated in this work.
251
4. LLE Modeling
252
4.1. COSMO-RS predictions
253
The Conductor-like Screening Model for Real Solvents is useful to describe the
254
thermodynamic behavior of pure compounds and mixture of compounds. The only input required
255
in COSMO-RS calculation is the chemical structure of the species involved in order to predict
256
the properties and behavior of the compounds and mixture of compounds. In COSMO-RS,
257
molecules are placed in a virtual conductor environment and divided into segments. Interaction
258
energies are calculated based on the energetic cost of removing the virtual conductor 11 ACS Paragon Plus Environment
Industrial & Engineering Chemistry Research
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 33
259
environment which is described in terms of the screening charge density, σ. Klamt and co-
260
workers derived the mathematical expressions for the interaction energies, chemical potential,
261
and activity coefficient and these can be referred in their work32,33 and are also summarized in
262
our previous works. 26,34,35
263
In order to do predictions using COSMO-RS, a .cosmo file must be generated for all
264
compounds i.e. cations, anions, HBD, thiophene, and n-heptane. First, geometry optimization
265
was performed at Hartree-Fock level and 6-31G* basis set. This was followed by a single point
266
calculation using density functional theory combined with Becke-Perdew functional and triple
267
zeta valence potential (TZVP) basis set to generate the .cosmo files. All computational work
268
were performed using TURBOMOLE software package.36 Thereafter, ternary LLE tie lines
269
based on the experimental feed composition of each ternary system were generated using
270
COSMOthermX software package with the parameterization file BP_TZVP_C30_1401.ctd.37
271
Description of ILs in COSMO-RS approaches are i) the electro-neutral approach, ii) the ion
272
pair approach, and iii) the meta-file approach. The first approach i.e. the electroneutral approach
273
is widely accepted to be the closest representation of ILs in COSMO-RS because it describes the
274
cation and anion as separate species in a liquid mixture. In this work, the same approach is
275
considered to represent DESs in COSMO-RS. The mathematical adaptation can be referred to in
276
our previous works. 26,34,35
277
Figs. 1 (a) – (d) show the ternary tie lines calculated using COSMO-RS along with the
278
experimental tie lines. The goodness of fit between the COSMO-RS predicted tie lines and the
279
experimental tie lines was evaluated using the root mean square deviation (RMSD) as defined in
280
Equation 4, where x is the concentration of species in mole fraction, and the subscripts i, j, and k
12 ACS Paragon Plus Environment
Page 13 of 33
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Industrial & Engineering Chemistry Research
281
designate the component, phase, and tie line, respectively, and m is the number of tie lines, c the
282
number of components, and j refers to the phases. 678 :78 ; ? = 0123 '%) = 1004∑> CA ∑BA ∑@A 9
283
9
=>?