1. Introduction
Aromatics produced by naphtha reforming and catalytic cracking are important feedstocks for many petrochemical applications [
1]. Nevertheless, olefin impurities are common in aromatic streams. Therefore, it is necessary to separate aromatics and olefins. For such separation, three techniques are commercially used depending on the solution concentration: (i) at low aromatic content (20~65 wt%), liquid–liquid extraction is usually used; (ii) at medium aromatic content (65 to 90 wt%), extractive distillation is usually used; and (iii) at high aromatic content (>90 wt%), azeotropic distillation is used. However, to date, there is no practical process for separation when the aromatic content in the feed mixture is less than 20 wt%, while liquid extraction is considered the most favorable process when the aromatic content is less than 20%.
The process of separating the components of a liquid stream by contacting that stream with another liquid stream, which may be insoluble or only partially soluble, is called liquid–liquid extraction. It is possible to separate the components because some of them have a preference to be more soluble in one of the liquid streams than in the others. According to Coquelet and Ramjugernath [
2], there are typically three distinct types of liquid–liquid equilibrium phase diagrams: (i) a binary component pair is partially miscible (type 1), (ii) two binary component pairs are somewhat miscible (type 2), and (iii) all three binary component pairs are somewhat miscible (type 3).
A critical step in the liquid–liquid extraction process is to find an efficient and cost-effective solvent. An ideal solvent should provide high extraction performance characterized by high solute selectivity, a high partition ratio, ease of regeneration, and a minimal feedstock to solvent ratio. From an environmental perspective, the solvent should also be environmentally friendly and non-toxic. From an economic point of view, the solvent should be available at low cost or should be able to be produced by a simple and cheap synthesis process. In addition, the physical and thermodynamic properties of the solvent, such as viscosity, thermal stability, density, and surface tension, are among the required information for industrial-scale applications. Traditional industrial-scale processes usually use organic solvents such as sulfolane (SUL), furfuryl alcohol, ethylene glycols, N-methylpyrrolidone and N-formylmorpholine.
However, the organic solvents have undesirable properties such as high toxicity, flammability, volatility and high cost of regeneration. Later, ionic liquid (IL) was introduced as a new advanced solvent and was extensively studied by researchers due to its invaluable advantage of negligible vapor pressure. To the best of our knowledge, the separation of toluene/1-hexene was described only by Meindersma et al. [
3] using 3-methyl-N-butyl-pyridinium dicyanamide ([3-mebupy]N(CN)
2) IL in a pilot plant with rotating disk contactor at T = 303.15 K. The selectivity obtained was also reported. The selectivity obtained ranged from 5 to 13, while the distribution ratio for toluene ranged from 0.258 to 0.350 [
3]. This result suggests that the IL has a higher affinity for toluene than for 1-hexene. This finding was also observed in aromatics/aliphatics separation, where the solvent tends to extract the aromatic compounds rather than the aliphatic ones [
4,
5,
6]. Even more interestingly, the tie lines reported were identical to those determined for toluene-heptane. Consequently, we can assume that most of the solvents used for the separation of toluene/heptane could also be used for the separation of toluene and 1-hexene.
However, despite the clear advantages of ILs, most are too expensive to be used at an industrial scale. They are also more difficult to synthesize than organic solvents and are not consistently environmentally friendly [
7,
8]. To overcome the limitations of ILs, deep eutectic solvents (DESs) have been explored as versatile substitutes for ILs and organic solvents. A DES usually consists of two or more components that combine through hydrogen bonds to form eutectic mixtures and are characterized by a melting point lower than that of the individual components. The components of DES are commonly referred to as hydrogen bond donor (HBD) and hydrogen bond acceptor (HBA). DESs represent an evolving class of green solvents. Most of them have negligible vapor pressure [
9], are biodegradable [
10], biocompatible, non-flammable, and non- or low-toxic [
11]. Several researchers have reported the use of different DESs for the separation of aromatic and aliphatic compounds [
12,
13] (
Table S1).
Kareem et al. [
14] used ETPI:SU and ETPI:EG DESs for the removal of toluene from heptane. However, it was found that EG or SU did not appear in the raffinate phase, indicating a critical problem in the removal of aromatics by liquid–liquid extraction. Wang et al. [
5] studied the effect of HBAs on the extraction of toluene. The DESs with bromide-based quaternary ammonium salts showed higher values in selectivity and the distribution ratio than those based on chloride. Hou et al. [
6] used tetrabutylphosphonium bromide (TBPBr) and tetrabutylammonium bromide (TBABr) as HBAs. The extraction rates of toluene for TBABr- and TBPBr-based DESs showed similar capacities, but TBPBr provided much higher selectivity than TBABr. In addition to HBAs, HBDs could also have a significant effect on extraction performance. Polyalcohols such as EG and TEG and carbonyl groups with carboxylic acid such as levulinic acid (LA) showed higher selectivities than sulfolane. Compared to EG, LA showed higher extraction efficiency, which was due to higher selectivity of toluene. This could be due to the fact that the carbon-oxygen double bond (C=O) present in the structure of LA could enhance its interaction with the aromatic ring of toluene by forming a π–π bond.
Numerous thermodynamic predictions and validations in critical steps such as denitrification [
15], desulfurization [
16,
17], and separation of aromatic and aliphatic mixtures [
18,
19,
20] have been performed using the Conductor-like Screening Model for Realistic Solvents (COSMO-RS). In this study, COSMO-RS screening and experimental liquid–liquid extraction were used to investigate and validate the extraction performance of DESs for the separation of toluene from 1-hexene.
2. COSMO-RS Screening
In the initial phase of this work, we performed a literature search and proposed more than 100 potential DESs for the separation of toluene from 1-hexene. These DESs were successfully synthesized, characterized, and used in many applications reported in the literature. They were mainly ammonium-, phosphonium-, and choline-based DESs with a variety of HBDs. However, when searching the database COSMO-RS provided by the supplier, it was found that some components of these DESs were not available. To address this issue, the molecular geometry of the unavailable compounds was calculated using the TurboMole program package, also provided by the same vendor. Geometry optimization of these compounds was successfully performed. Nevertheless, some compounds were omitted due to their polymer nature (PEG 200, PEG 400, PEG 600, PEG 1000, and PEG 4000) and the inability to import “
.cosmo” files into COSMOtherm (betaine, betaine hydrochloride, caffeic acid, choline acetyl chloride, p-toluenesulfonic acid, and trimethyl hydrochloride). These restrictions brought the DES candidates to a shortlist of 53 types (
Table 1) that were finally included in the COSMO-RS screening.
Screening Methodology
Geometry optimization of species not available in the database COSMO-RS was performed using the program package Turbomole (TmoleX). In this program, the chemical structure of the target molecule was first drawn. Then, geometry optimization was performed at the Hartree–Fock level and the 6-31G* basis set. To use COSMO-RS as a screening tool, you need to create the “.cosmo” files of the target components (salt cations, salt anions, and HBDs). The “.cosmo” file of a molecule contains information about the screening charge density (σ) of the segmented molecule in a virtual conductor environment. The “.cosmo” file was generated by a one-point calculation using density functional theory (DFT) with Becke–Perdew and the triple-ζ zeta valence potential (TZVP) basis set. Finally, the “.cosmo” files were exported to the COSMOthermX program with the parameterization BP _TZVP_C30_1301.ctd.
The activity coefficient at infinite dilution (
) of toluene and 1-hexene in each DES was predicted using the generated
“.cosmo” files in the COSMO-RS calculations. The selectivity of a DES to the solute compared to the carrier (
) can be expressed by the ratio of the activity coefficient for carrier and solute (Equation (1)).
In addition, the capacity at infinite dilution (C
∞) can be used to qualitatively determine the amount of a DES required for the extraction process. In this study, the capacity of a DES for solute (
indicates the maximum amount of solute that can be dissolved in the DES, and can be obtained by using the following equation:
The final parameter used to evaluate solvent properties in this extraction procedure is the performance index at infinite dilution (PI
∞). In this process, the two characteristics of capacity and selectivity are combined to estimate the overall performance of a DES. PI is simply expressed as the product of selectivity and capacity.
The performance of DESs for the separation of toluene and 1-hexene was compared using estimated C
∞ and S
∞ values. The BP_TZVPD_FINE_C30_1401.ctd parametrization file was used in the COSMOthermX program to perform the COSMO-RS calculations. DESs were represented using an electroneutral approach in the COSMO-RS approach as suggested by the developer of the COSMOtherm package for representing Ils. There are two other approaches that can be used in addition to the electroneutral approach to represent DESs in COSMO-RS, viz., metafile and ion-pair approaches. However, we chose the electroneutral approach because our previous studies have shown that it best describes the presence of Ils and DESs in the bulk mixture [
21,
22]. The COSMO-RS screening results with respect to C
∞ and S
∞ are shown in
Figure 1 and
Figure 2, respectively, and PI
∞ is represented in
Figure S1.
It is worth noting that the low C
∞ and S
∞ values were obtained in comparison to other systems in this work as well as in previous work with other aromatic-aliphatic separations [
18,
20,
23,
24]. This would support the extent of difficulty in separating 1-hexene and n-toluene. This would also necessitate experimental studies to validate the screening results, as the experimental results would reflect the actual extraction performance.
The COSMO-RS screening results in terms of capacity, selectivity, and the performance index at infinite dilution are reported in the (
supporting information Table S2).
Table 2 summarizes the top ten DESs from the COSMO-RS screening for each evaluation criteria (C
∞, S
∞, and PI
∞). It is noteworthy that TBABr:TEG (1:3) was not only highest for S∞, but was also among the top performers for C
∞ DES. It is expected that the highest value of C
∞ produced with BzTPPCl:TEG (1:8) is due to the high content of triethylene glycol (TEG), which increases the interaction with both 1-hexene and toluene. The potential DESs for experimental validation based on capacity and selectivity are TBABr:TEG (1:3), BzTPPCl:TEG (1:8), and TBABr:EG (1:4). Sigma surfaces of the HBAs, HBDs, as well as the sigma profile and sigma potential of three selected DESs for the toluene-1-hexene system are included in the (
supporting information Figures S2–S4). MTPPBr:TEG (1:4) was not selected because it had the lowest selectivity among the top ten DESs in the COSMO-RS screening. In addition, TBABr:LA (1:2) and TBABr:LA (1:3) were also rejected because levulinic acid appeared in the raffinate phase (
Figure S5).
Synthesis of these three potential candidates (BzTPPCl:TEG (1:8), TBABr:TEG (1:3), and TBABr:EG (1:4)) is then necessary to validate their actual performance. After they are successfully synthesized and characterized, the actual performance is determined in liquid–liquid extraction experiments. The liquid–liquid equilibria (LLE) data for each DES will be examined and supported with thermodynamic models and consistency tests.