Mutual Solubility of Aromatic Hydrocarbons in ... - ACS Publications

Jan 10, 2017 - The experimental LLE data were modeled using the CPA-EoS originally developed by Kontogeorgis et al.(64) This model has demonstrated a ...
0 downloads 0 Views 1MB Size
Article pubs.acs.org/jced

Mutual Solubility of Aromatic Hydrocarbons in Pyrrolidinium and Ammonium-Based Ionic Liquids and Its Modeling Using the CubicPlus-Association (CPA) Equation of State Patricia Fernández Requejo,† Ismael Díaz,‡ Emilio J. González,‡ and Á ngeles Domínguez*,† †

Advanced Separation Processes Group, Department of Chemical Engineering, University of Vigo, Campus Lagoas-Marcosende, 36310 Vigo, Spain ‡ Department of Industrial Chemical Engineering and Environmental, Technical University of Madrid, E-28006 Madrid, Spain S Supporting Information *

ABSTRACT: In this work, the mutual solubility of 30 binary systems {aromatic hydrocarbon (1) + ionic liquid (2)} involving benzene, toluene, ethylbenzene, o-xylene, m-xylene, and p-xylene as aromatics and the ionic liquids, 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide, [BMpyr][NTf2], 1-butyl-1-methylpyrrolidinium trifluoromethanesulfonate, [BMpyr][TfO], 1-butyl1-methylpyrrolidinium dicyanamide, [BMpyr][DCA], butyltrimethylammonium bis(trifluoromethylsulfonyl)imide, [N4111][NTf2], and tributylmethylammonium bis(trifluoromethylsulfonyl)imide, [N4441][NTf2], were determined in the temperature range of T = (293.15 to 333.15) K at p = 0.1 MPa. The influence of the temperature and the structural characteristics of the aromatic hydrocarbons as well as the effect of the structure of the ionic liquids on the phase behavior was analyzed. A comparison between the experimental LLE of the binary systems presented in this work and those found in literature was also carried out. Finally, all the binary systems were satisfactorily modeled using the Cubic-Plus-Association (CPA) Equation of State (EoS).

1. INTRODUCTION The extraction of aromatic hydrocarbons from the petrochemical streams is one of the main targets of the petrochemical industry1 since the aromatics are used as basic raw materials for the production of many different chemical products.2 However, this separation is challenging because of the presence of azeotropes and components in the mixtures with close boiling points.3,4 Furthermore, the new environmental regulations demand a reduction of the level of sulfur and of the aromatic hydrocarbons in the petrochemical products (gasoline, diesel fuel, or engine oils) to minimize their environmental impact.5,6 Therefore, the design of processes more efficient from an energetic point of view and less dangerous for the environment as well as the search of new solvents to replace the volatile organic solvents used in many industrial processes are two of the targets tackled by the scientific community. In this background, the ionic liquids (ILs) have been proposed as a suitable alternative to the separation agents used in the petrochemical industry because of their advantageous properties such as a negligible vapor pressure which reduces odor problems and atmospheric pollution, high chemical and thermal stability, nonflammability, and the possibility of modifying their physicochemical properties by combining different anions and cations.7−9 Additionally, the ILs present a very low solubility in the hydrocarbons, both aliphatic and © 2017 American Chemical Society

aromatic hydrocarbons, which allows one to reduce the number of steps of the recovery solvent meaning a reduction of the energy consumption of the separation process. However, before applying the ILs as solvents it is essential to know the mutual solubility of the binary mixtures {aromatic hydrocarbon (1) + ionic liquid (2)} since this information is crucial to design any process on an industrial scale. Moreover, phase equilibrium results allow us to enhance our knowledge about interactions between the ILs and other organic compounds. A database of experimental binary liquid−liquid equilibria (LLE) data of different ILs with a broad range of hydrocarbons has been reported. By inspection of the literature available in the Scopus database, most of the references (25) are related to mixtures containing imidazolium-based ILs, while scarce publications about the phase behavior of binary mixtures with ILs based on other cations such as ammonium (3), phosphonium (3), pyridinium (7), pyrrolidinium (6), or quinolinium (7) have been found.10−54 Because of the large number of {aromatic hydrocarbon (1) + ionic liquid (2)} potential mixtures, new solubility studies are needed since the Received: July 20, 2016 Accepted: December 23, 2016 Published: January 10, 2017 633

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642

Journal of Chemical & Engineering Data

Article

Table 1. Specifications of the Pure Compounds Used in This Worka chemical name

CAS no

source

mass fraction purityg

benzene ethylbenzene toluene o-xylene m-xylene p-xylene [BMpyr][NTf2]b [BMpyr][TfO]c [BMpyr][DCA]d [N4111][NTf2]e [N4441][NTf2]f

71-43-2 100-41-4 108-88-3 95-47-6 108-38-3 106-42-3 223437-11-4 367522-96-1 370865-80-8 258273-75-5 405514-94-5

VWR Sigma-Aldrich VWR Sigma-Aldrich Sigma-Aldrich Fluka IoLiTec IoLiTec IoLiTec IoLiTec IoLiTec

≥0.999 ≥0.998 ≥0.998 ≥0.990 ≥0.990 ≥0.990 >0.99 >0.99 >0.98 >0.99 >0.99

purity method analysisg

purification method

GCi

molecular sieves and degassing

NMRj, ICk

vacuum desiccation

water contenth (ppm)

152 114 72 62 77

a Standard uncertainty: u(water content) = 5 ppm. b1-Butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide. c1-Butyl-1-methylpyrrolidinium trifluoromethanesulfonate. d1-Butyl-1-methylpyrrolidinium dicyanamide. eButyltrimethylammonium bis(trifluoromethylsulfonyl)imide. fTributylmethylammonium bis(trifluoromethylsulfonyl)imide. gGiven by the suppliers. hMeasured after the purification method. iGas chromatography. jNuclear magnetic resonance. kIon chromatography.

vacuum (p = 0.2 Pa) for at least 48 h prior to use. Once dried, their water content was measured by Karl Fisher tritation (Mettler Toledo C20 Coulometric KF Titrator) with an uncertainty in the measurement of 5 ppm. In order to prevent water absorption, the ILs were kept and manipulated in a glovebox under argon atmosphere. A summary of the specifications of all of the pure compounds used in this work is shown in Table 1. 2.2. Apparatus. All of the samples were prepared by weighing using a Mettler AX-205 Delta Range balance with an uncertainty in the measurement of 3 × 10−4 g. Densities of both pure compounds and mixtures were measured using an Anton Paar DSA-5000 M digital vibrating tube densimeter with an uncertainty in the measurement of 3 × 10−3 g·cm−3 for the pure ILs and 1 × 10−3 g·cm−3 for the remaining samples. These uncertainties were estimated according to the guidelines reported by Chirico et al.65 The calibration of the densimeter was checked with known densities of pure compounds, Mili-Q water, and some alcohols. To ensure the two-phase equilibrium, a PoliScience thermostatic bath with a digital temperature controller was used, and the temperature was controlled with a digital thermometer ASL model F200 with an uncertainty in the measurement of 0.01 K. 2.3. Experimental Procedure. 2.3.1. Solubility Study. In general terms, the solubility of the aromatic hydrocarbons in the ILs increases when the temperature increases; nevertheless, some binary systems show a decrease in the solubility of the aromatics in the ILs when the temperature increases.11−15,52 Therefore, in order to know the influence of the temperature on the solubility of the aromatics in the selected ILs, previously to the LLE determination, a preliminary study of their solubility was carried out. To this end, the solubilities of the six aromatics in the pyrrolidinium and ammonium-based ILs were determined at three temperatures: room temperature, T = 333.15 K and T = 283.15 K, using the “cloud point” method.66 For this, approximately 3 mL of ionic liquid was placed in a vial, and the aromatic was added dropwise until a slight turbidity in the samples was observed at room temperature; then, the compositions of the samples were determined by weighing. Next, the vial was placed into the thermostatic bath at T = 333.15 K. After a few minutes, if the turbidity of the samples disappeared, the corresponding aromatic hydrocarbon was added drop by drop until turbidity was observed again; then,

solubility of aromatic compounds in the ILs affects their use as extraction media. This article is a continuation of our previous works on solubility studies of aromatic hydrocarbons in different ILs.14,55 In this case, the LLE data of 30 binary systems {aromatic hydrocarbon (1) + ionic liquid (2)} were experimentally determined in the temperature range of (293.15 to 333.15) K at p = 0.1 MPa. The aromatic hydrocarbons used were benzene, toluene, ethylbenzene, and the three xylenes isomers, and the ILs selected were ILs containing 1-butyl-1-methylpyrrolidinium as cation combined with the anions bis(trifluoromethylsulfonyl)imide, [BMpyr][NTf2], trifluoromethanesulfonate, [BMpyr][TfO], and dicyanamide, [BMpyr][DCA], and two ammonium-based ILs, butyltrimethylammonium bis(trifluoromethylsulfonyl)imide, [N4111][NTf2], and tributylmethylammonium bis(trifluoromethylsulfonyl)imide, [N4441][NTf2]. These ILs were chosen since both the pyrrolidinium and the ammonium-based ILs present lower toxicity in comparison with that of the imidazolium and pyridinium-based ILs and because their application as extraction agent has been hardly studied in the literature up to date.56−63 Similar systems, some of them with the same ILs studied in this work, have been previously studied by other authors and, therefore, literature LLE data were also taken into account for comparison purposes.27,39,50,52−54 By comparison of the systems reported in this work, the influence of the structural characteristics of both the aromatic hydrocarbons and the ILs as well as the effect of the temperature on the phase behavior were analyzed. Finally, the experimental LLE data were modeled using the Cubic-Plus-Association (CPA) Equation of State (EoS).64

2. EXPERIMENTAL SECTION 2.1. Chemicals. Benzene, toluene, ethylbenzene, o-xylene, m-xylene, and p-xylene were procured from different suppliers with purities greater than 0.99 in mass fraction (provided by the suppliers) in all cases. They were degassed using an ultrasonic cleaner BRANSON 3510 and subsequently dried over molecular sieves of 3 × 10−10 m (Sigma-Aldrich). Then, they were kept in bottles under argon atmosphere without any additional treatment. All of the ILs were supplied by IoLiTec GmbH with high purity. The ILs were dried by heating at T = 343.15 K under 634

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642

Journal of Chemical & Engineering Data

Article

Figure 1. Experimental LLE data of the binary systems: (a) {benzene (1) + ionic liquid (2)}, (b) {toluene (1) + ionic liquid (2)}, (c) {ethylbenzene (1) + ionic liquid (2)}, (d) {o-xylene (1) + ionic liquid (2)}, (e) {m-xylene (1) + ionic liquid (2)}, and (f) {p-xylene (1) + ionic liquid (2)} as a function of the mass fraction of the corresponding aromatic hydrocarbon in the IL-rich phase at T = (293.15 to 333.15) K and at p = 0.1 MPa. Dashed lines represent calculated data using the CPA-EoS. Ionic liquids: (◇) [BMpyr][NTf2], (□) [BMpyr][TfO], (○) [BMpyr][DCA], (△) [N4111][NTf2], and (▲) [N4441][NTf2].

in the determination of density vs composition curves involved in the calculation of the phase compositions. 2.3.2. Liquid−Liquid Equilibrium Determination. In order to determine the experimental LLE data, an immiscible mixture of two compounds was placed inside 15 mL glass vials, the total volume of each sample being 12 mL (4 mL of ionic liquid + 8 mL of aromatic hydrocarbon). This mixture was prepared in a glovebox under argon atmosphere and sealed using a rubber cover to avoid losses by evaporation or absorbance of moisture. Then, the vials were vigorously stirred for at least 3 h in order to ensure intimate contact between both phases at T = 293.15 K. Thereafter, they were settled for at least 3 h in a thermostatic

the sample was weighed, and its composition was calculated. If the turbidity of the samples does not disappear when the temperature rises, they were introduced in a thermostatic water bath at T = 283.15 K in order to determine if the solubility increases when the temperature decreases. Taking into account the results obtained in this preliminary study, it can be concluded that the solubility of the aromatic hydrocarbons in the ammonium-based ILs increases when the temperature is increased, while the solubility of the aromatics in the pyrrolidinium-based ILs decreases with the increase in the temperature. In both cases, this dependence of solubility with temperature is very smooth. These behaviors will be considered 635

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642

Journal of Chemical & Engineering Data

Article

bath to ensure complete phase separation. The temperature was controlled using a digital thermometer. Preliminary tests showed that the chosen time of stirring and of settling was enough to guarantee that equilibrium between phases was reached. Next, a sample of each phase was withdrawn using a syringe, the density of both phases was measured at the corresponding temperature, and the compositions were calculated using the polynomial expression previously obtained, as is described below. This procedure was repeated for each studied temperature from T = (293.15 K to 333.15) K, in 5 K steps. The determination of the IL-rich phase composition was carried out by means of measuring the densities and fitting these values versus the composition (ρ vs xi). Thus, prior to the experimental LLE determination, the density of the miscible binary mixtures with known compositions was determined. Bearing the above-mentioned results in mind, the density of the binary mixtures involving the ammonium-based ILs as solvents was measured at T = 343.15 K in order to ensure the complete miscibility of the binary mixtures, while the density of the binary mixtures containing the pyrrolidinium-based ILs was measured at T = 283.15 K and at p = 0.1 MPa for the same reason. The density values of binary mixtures as well as the polynomial expressions for ρ vs xi are shown in Table S1 and Table S2, respectively (Supporting Information). Besides, two binary mixtures were evaluated in order to obtain the error of the technique used for the determination of the solubility curves with the maximum error found being 0.005 in mole fraction. Finally, due to the fact that the tested ILs are practically immiscible in the aromatic hydrocarbons used in this work at the studied temperature range, the absence of ILs in the aromatic hydrocarbon-rich phase was assumed. This fact was checked with the values of the density of the aromatic hydrocarbon-rich phase since they were practically unchanged from those of the corresponding pure compound.

The effects of the structure of the studied aromatics as well as of ILs included in this work on the LLE data are analyzed below. 3.1.1. Influence of Aromatic Hydrocarbons. As can be inferred by comparing the different items in Figure 1a−e, the trend of the solubility of the aromatic hydrocarbons increases according to benzene > toluene > ethylbenzene in all of the studied ILs. This behavior could be due to the steric hindrance of the aromatics, given that the presence of one substituent as well as an increase in the alkyl chain length of the substituent on the benzene ring can increase steric hindrance and lead to a less effective package of the aromatic molecules in the IL-rich phase.52 Regarding the solubility of the xylene isomers in the ILs, the observed general trend is o-xylene > m-xylene ≥ pxylene. The solubility sequence of xylene isomers could be due to the impact on the aromatic induced dipole moment caused by the different positions of the second methyl group.14−17,52 This general trend can be observed in the LLE of all the presented binary systems with the exception of the binary systems involving [N4441][NTf2] in which the solubility of pxylene is much lower than those of the other xylene isomers in this IL. Finally, the general trend observed by many authors for the solubility of the aromatic hydrocarbons in the ILs is benzene > toluene > ethylbenzene > xylenes;14−17,52 nevertheless, the solubility sequence of the aromatics in the studied ILs follows in general the order benzene > toluene > o-xylene > ethylbenzene > m-xylene ≥ p-xylene (see Figure 1 or Table S2). For a better visualization of the conclusions drawn from the results, by way of example, the LLE data of the binary system {aromatic hydrocarbon (1) + [BMpyr][DCA] (2)} are displayed in Figure 2.

3. RESULTS AND DISCUSSION 3.1. Experimental Data. The experimental LLE data for the binary systems {aromatic hydrocarbon (1) + ionic liquid (2)} in the temperature range of T = (293.15 to 333.15) K, every 5 K, and at p = 0.1 MPa are displayed in Figure 1 and reported in Table S2 provided as Supporting Information. Because of the appreciable differences in the molar masses of the studied ILs, the graphical representation of the LLE data was performed in mass fraction, w, instead of molar fraction, x, in order to make a reliable analysis of the different influences. Moreover, as stated above, the studied ILs were assumed practically immiscible in all of the tested aromatic hydrocarbons at the studied temperature range; therefore, only the compositions of the corresponding aromatic hydrocarbon in the IL-rich phase are plotted in the figures presented in this work. The experimental LLE data presented in Figure 1 allow the analysis of the effect of the temperature on the solubility of the aromatics in the ILs. As has been mentioned in section 2.3.1, although the temperature effect on the phase behavior is small, a slight increase in the solubility of all the aromatics in the ammonium-based ILs with an increase in the temperature is observed at the studied temperature range, while a slight decrease in the solubility of the aromatics is found in the pyrrolidinium-based ILs. Despite the fact that this latter behavior is not typical for the binary systems, this trend has been previously reported for other mixtures.11−15,52

Figure 2. Experimental LLE data of the binary systems {aromatic hydrocarbon (1) + [BMpyr][DCA] (2)} as a function of the mass fraction of the corresponding aromatic hydrocarbon in the IL-rich phase: (●) benzene, (■) toluene, (○) ethylbenzene, (◇) o-xylene, (×) m-xylene, and (△) p-xylene at T = (293.15 to 333.15) K and at at p = 0.1 MPa.

3.1.2. Influence of Ionic Liquids. The effect of the structure of the anion and the cation of the ionic liquid on the phase equilibria can be analyzed from results shown in Figure 1. Comparing the solubilities of the aromatics in the pyrrolidinium-based ILs, two different phase behaviors are observed: when the aromatic hydrocarbon is benzene, the order followed by the solubilities is [BMpyr][DCA] > [BMpyr][NTf2] > 636

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642

Journal of Chemical & Engineering Data

Article

[BMpyr][TfO], while when the aromatic hydrocarbon is other than benzene, the sequence followed is [BMpyr][NTf2] > [BMpyr][DCA] > [BMpyr][TfO]. This unusual behavior observed for the binary system involving benzene may be related to the lower molecular mass and molar volume of the [DCA]− anion since the number of ions of ionic liquid in the vicinity of benzene, which is the smallest aromatic hydrocarbons, is higher than that when the ionic liquid contains a [NTf2]− or [TfO]− anion, with higher molar volumes. However, the comparison of the values of the solubility of the aromatic hydrocarbons in the ammonium-based ILs shows that the solubility sequence is [N4441][NTf2] > [N4111][NTf2] for all the studied aromatics. A possible explanation for this behavior is that the interaction between the IL cation and the aromatic is enhanced by the decrease of the IL cation polarity due to the increase of the alkyl chain length. Furthermore, the entropic effects that contribute to increase the solubility (asymmetry and free volume of the IL) are also enhanced by a longer alkyl chain.52 3.2. Comparison with the Literature. The comparison between the experimental LLE data of the binary systems in the temperature range of T = (293.15 to 333.15) K and at p = 0.1 MPa presented in this work and those data found in the literature27,39,50,52−54,56,62,63 are displayed in Figure 3. To our knowledge, only LLE data for the binary systems {benzene, or toluene or ethylbenzene (1) + [BMpyr][NTf2] (2)},56,62 {benzene, or toluene (1) + [BMpyr][TfO] (2)},27 and {benzene, or toluene, or ethylbenzene (1) + [BMpyr][DCA] (2)}52,63 were previously published. 3.2.1. Comparison of Binary Systems with the Same IL. As can be observed in Figure 3, the experimental LLE data of the systems {benzene or toluene (1) + [BMpyr][TfO] (2)} are in agreement with those found in the literature.27 The same conclusion is possible to draw from the comparison of the experimental and reported LLE data of the systems {benzene or toluene or ethylbenzene (1) + [BMpyr][NTf2] (2)} and {benzene (1) + [BMpyr][DCA] (2)}.56,62,63 Otherwise, some discrepancies are observed in the comparison of the experimental solubility data of the binary systems {benzene or toluene or ethylbenzene (1) + [BMpyr][DCA] (2)} with those reported in literature.52 These differences can be due to the different experimental methodology employed to determine the phase compositions. Paduszyński et al.52 determined the compositions of equilibrium using gravimetry. 3.2.2. Comparison of Binary Systems with the Same Alkyl Chain Length of the IL Cation. The comparison of the experimental LLE data for the systems {benzene or toluene (1) + [BMpyr][TfO] (2)} with those reported for the systems {benzene or toluene (1) + [BMim][TfO] (2)},27 {benzene or toluene (1) + [BM3py][TfO] (2)}27 are presented in Figure 3a,b. The solubilities of the aromatic hydrocarbons follow the order [BMpyr][TfO] ≈ [BMim][TfO] < [BM3py][TfO], both for benzene and for toluene. In the case of the comparison of the systems {benzene or toluene or ethylbenzene (1) + IL (2)} the ILs being [BMpyr][NTf2] and [BM4py][NTf2],39 it can be concluded that the solubility of the three aromatics is lower in [BMpyr][NTf2] than in its pyridinium-based homologue. 3.2.3. Comparison of Binary Systems with Different IL Anions. The influence of the anionic nature of the IL on the solubility of the aromatic hydrocarbons can be analyzed comparing the systems containing benzene, toluene, and ethylbenzene as aromatics and the ILs [BMpyr][NTf2], [BMpyr][TfO], [BMpyr][DCA], [BMpyr][TCM],50 [BMpyr]-

Figure 3. Experimental LLE data of the binary systems: (a) {benzene (1) + ionic liquid (2)}, (b) {toluene (1) + ionic liquid (2)}, and (c) {ethylbenzene (1) + ionic liquid (2)} as a function of the mass fraction of the corresponding aromatic hydrocarbon in the IL-rich phase at the 637

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642

Journal of Chemical & Engineering Data

Article

mixtures {aromatic hydrocarbon (1) + ionic liquid (2)} were obtained from the experimental LLE data. Parameter estimation was carried out using the maximum-likelihood objective function and the Britt-Luecke minimization algorithm.71 The CPA-EoS combines a typical cubic term (SRK) with an association term for the calculation of compressibility factor (Z):

Figure 3. continued studied temperatures and at p = 0.1 MPa. Ionic liquids: (□) [BMpyr][TfO] (this work), (○) [BMpyr][DCA] (this work), (◇) [BMpyr][NTf2] (this work), (■, gray) [BMpyr][TfO],27 (●, gray) [BMpyr][DCA],52 (●) [BMpyr][DCA],63 (▲, gray) [BMpyr][NTf2],62 (◆) [BMpyr][NTf2],56 (■) [BMim][TfO],27 (◪) [BM3py][TfO],27 (◆, gray) [BM4py][NTf2],39 (▽) [BMpyr][TCM],50 (△) [BMpyr][FPA],53 and (▼, gray) [BMpyr][TCB].54

Z = Z cubic + Z assoc

[FPA],53 and [BMpyr][TCB].54 By inspection of Figure 3, it can be inferred that the solubility of benzene follows the order [BMpyr][FPA] < [BMpyr][TfO] < [BMpyr][NTf 2 ] < [BMpyr][DCA] < [BMpyr][TCM] < [BMpyr][TCB], while the order followed for the solubilities of toluene and ethylbenzene is [BMpyr][FPA] < [BMpyr][TfO] < [BMpyr][DCA] < [BMpyr][NTf2] < [BMpyr][TCM] < [BMpyr][TCB]. Furthermore, it is important to comment that the effect of the temperature on the LLE data is rather small for all the studied binary systems and the systems found in the literature. Nevertheless, the dependence of the solubility of the aromatics with the temperature in the binary systems involving the IL [BMpyr][FPA] is significant, given that an increase in the temperature leads to an increase in the solubility of benzene, toluene, and ethylbenzene in this IL. As can be observed in Figure 3, the solubility of benzene and toluene is lower in [BMpyr][FPA] than in the remaining ILs at low temperatures, whereas higher solubilities are found when the temperature is increased. When the aromatic hydrocarbon is ethylbenzene, the solubility values in [BMpyr][FPA] are lower than those obtained for the rest ILs in the comparison. 3.3. Thermodynamic Modeling Using the CPA-EoS. The experimental LLE data were modeled using the CPA-EoS originally developed by Kontogeorgis et al.64 This model has demonstrated a good ability to correlate phase equilibria data for mixtures involving ILs.67−69 Data fitting was carried out using the data regression tool implemented in a commercial process flow sheet simulator (Aspen Plus). First, pure component parameters were obtained by fitting density and vapor pressure, and then binary interaction parameters of the

(1)

Z cubic =

vm a − vm − b RT (vm + b)

Z assoc =

∑ xi ∑ ρi ∑ ⎢⎢⎜⎜ i

Ai

(3)

∑ ∑ xixjaij i

b=

⎤ 1 1 ⎞ δX A i ⎥ − ⎟⎟ 2 ⎠ δρi ⎥⎦ ⎣⎝ X A i ⎡⎛

i

a=

(2)

(4)

j

∑ xibi

(5)

i

aij =

aiaj (1 − kij)

ai = a0i(1 + c1i(1 −

(6)

Tri ))2

(7)

where ρ is the molar density, R is the universal gas constant, a and b are the SRK parameters, vm is the mixture molar volume, xi is the mole fraction of component i, XAi is the fraction of sites of type A at molecule i that does not form any bonds with other active sites, and Tri is the reduced temperature of the component i. This fraction is calculated as X Ai =

1 (1 + ρ ∑i xi ∑ B X BjΔA iBj )

(8)

j

with ⎤ ⎡ ⎛ ε A iBj ⎞ ΔA iBj = bijβ A iBjg (ρ)⎢exp⎜ ⎟ − 1⎥ ⎦ ⎣ ⎝ RT ⎠

(9)

Table 2. CPA Pure Component Parameters Aromatic Hydrocarbon parameter

benzene

toluene

mm Pcm (bar) Tcm (K) εAiBj (K) βAiBj ρ deviation Pv deviation

0.7703 55.02 572.8 0.0 0 toluene > o-xylene > ethylbenzene > m-xylene ≥ p-xylene. Moreover, the analysis of the influence of the change of the anion of the ionic liquid shows that the solubility decrease in the order [BMpyr][NTf2] > [BMpyr][DCA] > [BMpyr][TfO] for all of the studied aromatics with the exception of the solubility of benzene that follows the order [BMpyr][DCA] > [BMpyr][NTf2] > [BMpyr][TfO]. Regarding the effect of the

REFERENCES

(1) Gaile, A. A.; Zalishchevskii, G. D.; Gafur, N. N.; Semenov, L. V.; Varshavski, O. M.; Fedyanin, N. P.; Koldobskaya, L. L. Removal of aromatic hydrocarbons from reforming naphtha. Combined extraction − Extractive-azeotropic distillation process. Chem. Technol. Fuels Oils 2004, 40, 215−221. (2) Weissermel, K.; Arpe, H. J. Industrial Organic Chemistry, 4th ed.; Wiley-WCH: Weinheim, Germany, 2003. (3) Meindersma, G. W.; Podt, A. J. G.; de Haan, A. B. Selection of ionic liquids for the extraction of aromatic hydrocarbons from aromatic/aliphatic mixtures. Fuel Process. Technol. 2005, 87, 59−70. (4) Domańska, U.; Pobudkowska, A.; Królikowski. Separation of aromatic hydrocarbons from alkanes using ammonium ionic liquid C2NTf2 at T = 298.15 K. Fluid Phase Equilib. 2007, 259, 173−179. (5) Dixson-Decléve, S. Global fuel quality trends. Proceedings of the 6th International Colloquium Fuels; Bartz, W. J., Ed.; Esslingen, Germany, 2007; pp 3−5. (6) Kędra-Królik, K.; Fabrice, M.; Jaubert, J.-N. Extraction of thiophene or pyridine from n-heptane using ionic liquids gasoline and diesel desulfurization. Ind. Eng. Chem. Res. 2011, 50, 2296−2306. (7) Earle, M. J.; Esperança, J. M. S. S.; Gilea, M. A.; Canongia Lopes, J. N.; Rebelo, L. P. N.; Magee, J. W.; Seddon, K. R.; Widegren, J. A. The distillation and volatility of ionic liquids. Nature 2006, 439, 831− 834. 640

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642

Journal of Chemical & Engineering Data

Article

(8) Baranyai, K. J.; Deacon, G. B.; MacFarlane, D. R.; Pringle, J. M.; Scott, J. L. Thermal degradation of ionic liquids at elevated temperatures. Aust. J. Chem. 2004, 57, 145−147. (9) Smiglak, M.; Reichert, W. M.; Holbrey, J. D.; Wilkes, J. S.; Sun, L. Y.; Thrasher, J. S.; Kirichenko, K.; Singh, S.; Katritzky, A. R.; Rogers, R. D. Combustible ionic liquids by design: is laboratory safety another ionic liquid myth? Chem. Commun. 2006, 24, 2554−2556. (10) Domańska, U.; Marciniak, A. Solubility of 1-alkyl-3-methylimidazolium hexafluorophosphate in hydrocarbons. J. Chem. Eng. Data 2003, 48, 451−456. (11) Łachwa, J.; Szydłowski, J.; Makowska, A.; Seddon, K. R.; Esperança, J. M. S. S.; Guedes, H. J. R.; Rebelo, L. P. N. Changing from an unusual high-temperature demixing to a UCST-type in mixtures of 1-alkyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}amide and arenes. Green Chem. 2006, 8, 262−267. (12) Domańska, U.; Laskowska, M.; Pobudkowska, A. Phase equilibria study of the binary Systems (1-butyl-3-methylimidazolium thiocyanate ionic liquid + organic solvent or water). J. Phys. Chem. B 2009, 113, 6397−6404. (13) Domańska, U.; Królikowski, M.; Arasimowicz, M. Phase equilibria of (1-hexyl-3-methylimidazolium thiocyanate + water. alcohol. or hydrocarbon) binary systems. J. Chem. Eng. Data 2010, 55, 773−777. (14) González, E. J.; Requejo, P. F.; Domínguez, Á .; Macedo, E. A. Phase equilibria of binary mixtures (ionic liquid + aromatic hydrocarbon): Effect of the structure of the components on the solubility. Fluid Phase Equilib. 2013, 360, 416−422. (15) Shang, H.-T.; Wu, J.-S.; Zhou, Q.; Wang, L.-S. Solubilities of benzene. toluene. and xylene isomers in 1-butyl-3-methylimidazolium tetrafluoroborate. J. Chem. Eng. Data 2006, 51, 1286−1288. (16) Shiflett, M. B.; Niehaus, A. M. S. Liquid−liquid equilibria in binary mixtures containing substituted benzenes with ionic liquid 1ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide. J. Chem. Eng. Data 2010, 55, 346−353. (17) García, J.; Torrecilla, J. S.; Fernández, A.; Oliet, M.; Rodríguez, F. (Liquid + liquid) equilibria in the binary systems (aliphatic. or aromatic hydrocarbons + 1-ethyl-3-methylimidazolium ethylsulfate. or 1-butyl-3-methylimidazolium methylsulfate ionic liquids). J. Chem. Thermodyn. 2010, 42, 144−150. (18) Domańska, U.; Pobudkowska, A.; Eckert, F. Liquid−liquid equilibria in the binary systems (1.3-dimethylimidazolium. or 1-butyl3-methylimidazolium methylsulfate + hydrocarbons). Green Chem. 2006, 8, 268−276. (19) González, E. J.; Bottini, S. B.; Pereda, S.; Macedo, E. A. Liquid− liquid equilibria of binary systems {benzene + [x-Mim][NTf2] ionic liquid}: Experimental data and thermodynamic modeling using a group contribution equation of state. Fluid Phase Equilib. 2014, 362, 163−169. (20) Domańska, U.; Laskowska, M.; Marcianiak, A. Phase equilibria of (1-ethyl-3-methylimidazolium ethylsulfate + hydrocarbon. + ketone. and + ether) binary systems. J. Chem. Eng. Data 2008, 53, 498−502. (21) Yokozeki, A.; Shiflett, M. B. Binary and ternary phase diagrams of benzene. hexafluorobenzene. and ionic liquid [emim][Tf2N] using equations of state. Ind. Eng. Chem. Res. 2008, 47, 8389−8395. (22) Bendová, M.; Wagner, Z. Thermodynamic description of liquid−liquid equilibria in systems 1-ethyl-3-methylimidazolium ethylsulfate + C7-hydrocarbons by polymer-solution models. Fluid Phase Equilib. 2009, 284, 80−85. (23) Domańska, U.; Królikowska, M. Phase equilibria study of the binary systems (1-butyl-3-methylimidazolium tosylate ionic liquid + water. or organic solvent). J. Chem. Thermodyn. 2010, 42, 355−362. (24) Marciniak, A. Influence of anion structure on the liquid−liquid equilibria of 1-ethyl-3-methyl-imidazolium cation based ionic liquidhydrocarbon binary systems. J. Chem. Eng. Data 2011, 56, 368−374. (25) Calado, M. S.; Petrovski, Z.; Manic, M. S.; Najdanovic-Visak, V.; Macedo, E. A.; Visak, Z. P. Liquid−liquid equilibria of imidazolium ionic liquids having bistriflamide or triflate anions with aromatic nonpolar and polar compounds. Fluid Phase Equilib. 2013, 337, 67−72.

(26) Larriba, M.; Navarro, P.; García, J.; Rodríguez, F. Liquid−liquid extraction of toluene from n-heptane by {[emim][TCM] + [emim][DCA]} binary ionic liquid mixtures. Fluid Phase Equilib. 2014, 364, 48−54. (27) Marciniak, A.; Karczemna, E. Influence of cation structure on binary liquid−liquid equilibria for systems containing ionic liquids based on trifluoromethanesulfonate anion with hydrocarbons. J. Phys. Chem. B 2010, 114, 5470−5474. (28) Domańska, U.; Pobudkowska, A.; Wisniewska, A. Solubility and excess molar properties of 1.3-dimethylimidazolium methylsulfate. or 1-butyl-3-methylimidazolium methylsulfate. or 1-butyl-3-methylimidazolium octylsulfate ionic liquids with n-alkanes and alcohols: Analysis in terms of the PFP and FBT models. J. Solution Chem. 2006, 35, 311− 334. (29) Domańska, U.; Marciniak, A. Liquid phase behaviour of 1-butyl3-methylimidazolium 2-(2-methoxyethoxy)-ethylsulfate with organic solvents and water. Green Chem. 2007, 9, 262−266. (30) Pereiro, A. B.; Rodríguez, A. Binary mixtures containing OMIM PF6: Density. speed of sound. refractive index and LLE with hexane. heptane and 2-propanol at several temperatures. Phys. Chem. Liq. 2008, 46, 162−174. (31) Domańska, U.; Marciniak, A. Liquid phase behavior of 1hexyloxymethyl- 3-methyl-imidazolium-based ionic liquids with hydrocarbons: The influence of anion. J. Chem. Thermodyn. 2005, 37, 577− 585. (32) Domań s ka, U.; Marciniak, A. Phase behavior of 1hexyloxymethyl-3-methyl-imidazolium and 1.3-dihexyloxymethyl-imidazolium based ionic liquids with alcohols. water. ketones and hydrocarbons: The effect of cation and anion on solubility. Fluid Phase Equilib. 2007, 260, 9−18. (33) Domańska, U.; Mazurowska, L. Solubility of 1.3-dialkylimidazolium chloride or hexafluorophosphate or methylsulfonate in organic solvents: Effect of the anions on solubility. Fluid Phase Equilib. 2004, 221, 73−82. (34) Domańska, U.; Bakala, I.; Pernak, J. Phase equilibria of an ammonium ionic liquid with organic solvents and water. J. Chem. Eng. Data 2007, 52, 309−314. (35) Domańska, U.; Lugowska, K.; Pernak, J. Phase equilibria of didecyldimethylammonium nitrate ionic liquid with water and organic solvents. J. Chem. Thermodyn. 2007, 39, 729−736. (36) Domańska, U.; Paduszyński, K. (Solid + liquid) and (liquid + liquid) phase equilibria measurements and correlation of the binary systems {tri-iso-butyl(methyl)phosphonium tosylate + alcohol. or + hydrocarbon}. Fluid Phase Equilib. 2009, 278, 90−96. (37) Blesic, M.; Lopes, J. N. C.; Gomes, M. F. C.; Rebelo, L. S. P. N. Solubility of alkanes. alkanols. and their fluorinated counterparts in tetraalkylphosphonium ionic liquids. Phys. Chem. Chem. Phys. 2010, 12, 9685−9692. (38) Domańska, U.; Królikowski, M.; Pobudkowska, A.; Letcher, T. M. Phase equilibria study of the binary systems (N-butyl-4methylpyridinium tosylate ionic liquid + organic solvent. or water). J. Chem. Eng. Data 2009, 54, 1435−1441. (39) Domańska, U.; Królikowski, M.; Ramjugernath, D.; Letcher, T. M.; Tumba, K. Phase equilibria and modeling of pyridinium-based ionic liquid solutions. J. Phys. Chem. B 2010, 114, 15011−15017. (40) Dreiseitlová, J.; Reháh, K.; Vreekamp, R. Mutual solubility of pyridinium-based tetrafluoroborates and toluene. J. Chem. Eng. Data 2010, 55, 3051−3054. (41) Domań s ka, U.; Roguszewska, M.; Kró l ikowski, M.; Ramjugernath, D.; Naidoo, P. Phase equilibria study of binary systems comprising an (ionic liquid + hydrocarbon). J. Chem. Thermodyn. 2015, 83, 90−96. (42) Domańska, U.; Królikowski, M.; Paduszynski, K. Phase equilibria study of the binary systems (N-butyl-3-methylpyridinium tosylate ionic liquid + an alcohol). J. Chem. Thermodyn. 2009, 41, 932−938. (43) Królikowski, M.; Karpińska, M.; Zawadzki, M. Phase equilibria study of the binary systems (n-hexylisoquinolinium thiocyanate ionic 641

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642

Journal of Chemical & Engineering Data

Article

liquid + organic solvent or water). J. Phys. Chem. B 2012, 116, 4292− 4299. (44) Królikowska, M.; Karpińska, M. Phase equilibria study of the (N-octylisoquinolinium thiocyanate ionic liquid + aliphatic and aromatic hydrocarbon. or thiophene) binary systems. J. Chem. Thermodyn. 2013, 63, 128−134. (45) Zawadzki, M.; Domańska, U. Thermodynamic properties of the N-octylquinolinium bis{(trifluoromethyl)sulfonyl}imide. J. Chem. Thermodyn. 2012, 48, 276−283. (46) Domańska, U.; Zawadzki, M. Thermodynamic properties of the N-butylisoquinolinium bis(trifluoromethylsulfonyl)imide. J. Chem. Thermodyn. 2011, 43, 989−995. (47) Domańska, U.; Zawadzki, M.; Zwolińska, M. Thermophysical properties and phase equilibria study of the binary systems {Nhexylquinolinium bis(trifluoromethylsulfonyl)imide + aromatic hydrocarbons. or an alcohol}. J. Chem. Thermodyn. 2011, 43, 775−781. (48) Domańska, U.; Zawadzki, M.; Królikowski, M.; Lewandrowska, A. Phase equilibria study of binary and ternary mixtures of {Noctylisoquinolinium bis{(trifluoromethyl)sulfonyl}imide + hydrocarbon. or an alcohol. or water}. Chem. Eng. J. 2012, 181−182, 63−71. (49) Domańska, U.; Zawadzki, M.; Tshibangu, M. M.; Ramjugernath, D.; Letcher, T. M. Phase equilibria study of {N-butylquinolinium bis{(trifluoromethyl)sulfonyl}imide + aromatic hydrocarbons. or an alcohol} binary systems. J. Chem. Thermodyn. 2010, 42, 1180−1186. (50) Domańska, U.; Lukoshko, E. V.; Szydłowski, J. Phase behavior of tricyanomethanide-based ionic liquids with alcohols and hydrocarbons. Fluid Phase Equilib. 2015, 387, 18−23. (51) Ferreira, A. R.; Freire, M. G.; Ribeiro, J. C.; Lopes, F. M.; Crespo, J. G.; Coutinho, J. A. P. An overview of the liquid-liquid equilibria of (ionic liquid + hydrocarbon) binary systems and their modeling by the Conductor-like Screening Model for real solvents. Ind. Eng. Chem. Res. 2011, 50, 5279−5294. (52) Paduszyński, K.; Lukoshko, E. V.; Królikowski, M.; Domańska, U.; Szydłowski, J. Thermodynamic study of binary mixtures of 1-butyl1-methylpyrrolidinium dicyanamide ionic liquid with molecular solvents: new experimental data and model with PC-SAFT equation of state. J. Phys. Chem. B 2015, 119, 543−551. (53) Domańska, U.; Lukoshko, E. V.; Królikowski, M. Phase behaviour of ionic liquid 1-butyl-1-methylpyrrolidinium tris(pentafluoroethyl)trifluorophosphate with alcohols. water and aromatic hydrocarbons. Fluid Phase Equilib. 2013, 345, 18−22. (54) Paduszyński, K.; Lukoshko, E. V.; Królikowski, M.; Domańska, U. Measurements and equation-of-state modelling of thermodynamic properties of binary mixtures of 1-butyl-1-methylpyrrolidinium tetracyanoborate ionic liquid with molecular compounds. J. Chem. Thermodyn. 2015, 90, 317−326. (55) González, E. J.; Requejo, P. F.; Maia, F. M.; Domínguez, Á .; Macedo, E. A. Solubility. density and excess molar volume of binary mixtures of aromatic compounds and common ionic liquids at T = 283.15 K and atmospheric pressure. Phys. Chem. Liq. 2015, 53, 419− 428. (56) Pereiro, A. B.; Veiga, H. I. M.; Esperança, J. M. S. S.; Rodríguez, A. Effect of temperature on the physical properties of two ionic liquids. J. Chem. Thermodyn. 2009, 41, 1419−1423. (57) Marciniak, A. Influence of cation and anion structure of the ionic liquid on extraction processes based on activity coefficients at infinite dilution. A review. Fluid Phase Equilib. 2010, 294, 213−233. (58) Marciniak, A.; Królikowski, M. Ternary liquid−liquid equilibria of bis(trifluoromethylsulfonyl)-amide based ionic liquids + thiophene + n-heptane. The influence of cation structure. Fluid Phase Equilib. 2012, 321, 59−63. (59) Petkovic, M.; Seddon, K. R.; Rebelo, L. P. N.; Silva Pereira, C. Ionic liquids: a pathway to environmental acceptability. Chem. Soc. Rev. 2011, 40, 1383−1403. (60) Frade, R. F. M.; Rosatella, A. A.; Marques, C. S.; Branco, L. C.; Kulkarni, P. S.; Mateus, N. M. M.; Afonso, C. A. M.; Duarte, C. M. M. Toxicological evaluation on human colon carcinoma cell line (CaCo2) of ionic liquids based on imidazolium. guanidinium. ammonium.

phosphonium. pyridinium and pyrrolidinium cations. Green Chem. 2009, 11, 1660−1665. (61) Á lvarez-Guerra, M.; Irabien, A. Design of ionic liquids: an ecotoxicity (Vibrio f ischeri) discrimination approach. Green Chem. 2011, 13, 1507−1516. (62) Requejo, P. F.; Calvar, N.; Dominguez, A.; Gomez, E. Comparative study of the LLE of the quaternary and ternary systems involving benzene, n-octane, n-decane and the ionic liquid [BMpyr] [NTf2]. J. Chem. Thermodyn. 2016, 98, 56−61. (63) Requejo, P. F.; Calvar, N.; Dominguez, A.; Gomez, E. Determination and correlation of (liquid + liquid) equilibria of ternary and quaternary systems with octane, decane, benzene and [BMpyr][DCA] at T = 298.15 K and atmospheric pressure. J. Chem. Thermodyn. 2016, 94, 197−203. (64) Kontogeorgis, G.M.; Voutsas, E.C.; Yakoumis, I.V.; Tassios, D. P. An Equation of State for Associating Fluids. Ind. Eng. Chem. Res. 1996, 35, 4310−4318. (65) Chirico, R. D.; Frenkel, M.; Magee, J. W.; Diky, V.; Muzny, C. D.; Kazakov, A. F.; Kroenlein, K.; Abdulagatov, I.; Hardin, G. R.; Acree, W. E., Jr.; et al. Improvement of quality in publication of experimental thermophysical property data: challenges. Assessment tools. global implementation. and online support. J. Chem. Eng. Data 2013, 58, 2699−2716. (66) Letcher, T. M.; Siswana, P. M. Liquid−liquid equilibria for mixtures of an alkanol + water + a methyl substituted benzene at 25°C. Fluid Phase Equilib. 1992, 74, 203−217. (67) Haghtalab, A.; Kheiri, A. High pressure measurement and CPA equation of state for solubility of carbon dioxide and hydrogen sulfide in 1-butyl-3-methylimidazolium acetate. J. Chem. Thermodyn. 2015, 89, 41−50. (68) Manic, M. S.; Queimada, A. J.; Macedo, E. A.; Najdanovic-Visak, V. High-pressure solubilities of carbon dioxide in ionic liquids based on bis(trifluoromethylsulfonyl)imide and chloride. J. Supercrit. Fluids 2012, 65, 1−10. (69) Sousa, J. M. M. V.; Granjo, J. F. O.; Queimada, A. J.; Ferreira, A. G. M.; Oliveira, N. M. C.; Fonseca, I. M. A. Solubilities of hydrofluorocarbons in ionic liquids: Experimental and modelling study. J. Chem. Thermodyn. 2014, 73, 36−43. (70) Maia, F. M.; Tsivintzelis, I.; Rodriguez, O.; Macedo, E. A.; Kontogeorgis, G. M. Equation of state modelling of systems with ionic liquids: Literature review and application with the Cubic Plus Association (CPA) model. Fluid Phase Equilib. 2012, 332, 128−143. (71) AspenTech. Aspen Plus v8.8 Software Documentation; Aspen Technology: Germany, 2015. (72) Kontogeorgis, G. M.; Folas, G. K. Thermodynamic Models for Industrial Applications: From Classical and Advanced Mixing Rules to Association Theories; Wiley: Chichester, UK, 2009. (73) Valderrama, J. O.; Rojas, R. E. Critical Properties of Ionic Liquids. Revisited. Ind. Eng. Chem. Res. 2009, 48, 6890−6900.

642

DOI: 10.1021/acs.jced.6b00655 J. Chem. Eng. Data 2017, 62, 633−642