Removal of Thiophene from Mixtures with n-Heptane by Selective

Jul 7, 2016 - This work investigates the use of deep eutectic solvents (DESs) to extract sulfur-based compounds from n-heptane as model diesel compoun...
0 downloads 12 Views 977KB Size
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

=>?