Phase Behaviour of 1-Ethyl-3-methylimidazolium Thiocyanate Ionic Liquid with Catalytic Deactivated Compounds and Water at Several Temperatures: Experiments and Theoretical Predictions

Density, surface tension and refractive index were determined for the binary mixture of catalytic deactivated compounds with 1-ethyl-3-methylimidazolium thiocyanate {[EMIM][SCN]} at temperature of (298.15 to 323.15) K. For all the compounds with ILs, the densities varied linearly in the entire mole fraction with increasing temperature. From the obtained data, the excess molar volume and deviation of surface tension and refractive index have been calculated. A strong interaction was found between similar (cation-thiophene or cation-pyrrole) compounds. The interaction of IL with dissimilar compounds such as indoline and quinoline and other multiple ring compounds was found to strongly depend on the composition of IL at any temperatures. For the mixtures, the surface tension decreases in the order of: thiophene > quinoline > pyridine > indoline > pyrrole > water. In general from the excess volume studies, the IL-sulphur/nitrogen mixture has stronger interaction as compared to IL-IL, thiophene-thiophene or pyrrole-pyrrole interaction. The deviation of surface tension was found to be inversely proportional to deviation of refractive index. The quantum chemical based COSMO-RS was used to predict the non-ideal liquid phase activity coefficient for all mixtures. It indicated an inverse relation between activity coefficient and excess molar volumes.


Introduction
Zero emission of nitrogen and sulphur in diesel oil production is one of the largest challenges for petroleum refineries in the world.However the hydrodenitrification (HDN) and hydrodesulphurization (HDS) processes represent an even more important challenge.Most of the refineries in the world are able to produce diesel oil with above 70 ppm of nitrogen [1] and 500 ppm of sulphur [2] by using high hydrogen consumption, large reactor volume, and severe operating conditions [3].Thus the refinery must remove the aromatic compounds of nitrogen/sulphur such as pyrrole, indole, indoline, carbazole, benzocarbazole, pyridine, quinoline, benzoquinoline, thiophene, benzothiophene, dibenzothiophene, and their derivatives [1,2,4].Moreover the catalysts having higher activity or selectivity do not interact favourably with such compounds, even at moderate temperatures [5].On the other hand the nature of the compounds coupled with its chemical structure posseses more sterical hindrance as compared to the catalysts.The HDN and HDS processes have greater ability in removing nonaromatic structure of nitrogen and sulphur molecules from diesel oil with high efficiency although not economically [6,7].
Over the past few years Ionic Liquid (ILs) are alternative solvent for organic solvents in various applications and are particularly used as an extracting solvent for the removal of sulphur from fuel oil [7][8][9][10][11].The physiochemical properties of ILs are easily adjustable by changing cation and anion combination.Further it has negligible vapour pressure, posses nonflammability and high thermal and chemical stability.ILs provide better solution for removal of aromatic skeleton structure with multiple rings of nitrogen and sulphur compounds in diesel oil via CH-π bond interaction, ππ stacking, and n-π bond interaction [12][13][14][15][16]. Additionally ILs have higher LUMO/HOMO energy and high polarity in mixtures with such compounds [17].There are several International Journal of Chemical Engineering works which report desulphurization using ionic liquid with a target of below 10 ppm of total sulphur content in diesel and other liquid fuels.However these works are carried out without proper screening of ionic liquid (combinations as large as 10 18 ILs!!).
In our previous work we have reported that smaller size of cation has higher selectivity but lower capacity for the removal of thiophene using COnductor-Like Screening MOdel for Real Solvent (COSMO-RS) predictions at T = 298.15K and atmospheric pressure [18].Further we have studied the removal of pyrrole, indole, indoline, carbazole, benzocarbazole, pyridine, quinoline, benzoquinoline, thiophene, benzothiophene, dibenzothiophene, and their derivatives from diesel oil using smaller size cation based on imidazolium, pyridinium, pyrrolidinium, piperidinium, morpholinuim, and triazolium with fluoro-and nonfluoro anions [13,14].According to these results we have chosen few ionic liquids for the removal of such compounds via Liquid-Liquid Extraction (LLE) at T = 298.15K.According to our work via COSMO-RS prediction [12][13][14][15][16] the IL: 1-ethyl-3-methylimidazolium thiocyanate ([EMIM][SCN]) was the most efficient IL for the removal of pyrrole, pyridine, indoline, quinoline, and thiophene molecules.Thus in this work we have studied the phase behaviour of [EMIM] [SCN] with such compounds at p = 1 bar.
In a recent work, Doma ńska et al. [19] studied the phase behaviour of the binary system of thiophene with 1-ethyl-3-methylimidazolium thiocyanate and other ionic liquids.In this work, they measured densities, surface tension, and refractive index of mixtures of [EMIM][SCN] with pyrrole, pyridine, indoline, quinoline, thiophene, and water.Finally the quantum chemical-based COSMO-RS (Conductor-like Screening MOdel for Real Solvents) [12][13][14][15] model has been used to predict the activity coefficient for known composition of mixtures.This model is independent of experimental data and uses the molecular structure to determine the Screening Charge Densities (SCDs) or the sigma profile of the pure molecule.This is the only indicator in computing the chemical potential of a component in the solution.The computational details of COSMO-RS will be described briefly in the later section.) ionic liquid was purchased from Fluka with a purity of >95%.In order to increase the purity and for reducing the water content and volatile compounds to negligible values, vacuum (0.1 Pa) for at least 48 h was applied to the IL sample prior to the measurements.A water content of 0.044% in [EMIM][SCN] was obtained by Karl Fischer titrator (KF 787, Metrohm).Pyrrole (C 4 H 5 N) and indoline (C 8 H 9 N) were supplied from Spectrochem, India with mass fraction of 99%.Pyridine (C 5 H 5 N), quinoline (C 9 H 7 N) and thiophene (C 4 H 4 S) were received from Sigma-Aldrich with mass fractions greater than 98% and were used without further purification.The pure component structures, experimental densities, refractive index, and surface tension are given in Tables 1 to 3. A comparison between reported and experimental measurements was carried out on all the properties.

Experimental Section
Densities of the pure components and binary mixture were measured at atmospheric pressure with Anton Paar DSA-4500MA digital vibrating U-tube densimeter.The densitometer has a well-defined thermoelectric temperature control system.The uncertainty in the density measurement is ±0.0001 g•cm −3 .The apparatus was calibrated by measuring the density of Millipore quality water and ambient air.The surface tension of the pure components and binary mixtures was measured with a tensiometer by plate type method (Hanging drop tensiometer method, kruss k9, Germany) with a precision of 0.01 mN/m.The apparatus was calibrated by measuring the surface tension of Millipore quality water at ambient temperature.Refractive index of the pure components and binary mixtures was determined at ambient temperature using an automatic refractometer AD-13 model (ABBEMAT-WR Dr. Kernchen), with an uncertainty of ±0.00004.
Samples were prepared by transferring known mass of the pure liquids into stoppered bottles via syringe.The stoppered bottles are closed with screw caps to seal and prevent evaporation.All weighing was carried out in a balance (Mettler Toledo AT 261) with an accuracy of ±10 −4 gm.Previous experiments [18] showed that equilibrium was attained after 6 hours of stirring at 100 rpm at temperature at 298.15 K, using circulating water bath along with automatic controller.Equilibrium was attained by keeping the mixture settled for 12 hours during which good contact was obtained between two pure components.Samples from the homogeneous mixture were withdrawn using syringes for the thermodynamic phase behavior studies.All the samples were prepared immediately prior to performing density, surface tension, and refractive index measurements so as to prevent variation in composition due to water/air retention via the hygroscopic IL.

Computational Details
The sigma profile generation and methodologies are given in detail in our previously published work on ionic liquid systems [12][13][14][15].The SCDs obtained in COSMO calculations are given for each segment.Thereafter the expression for chemical potential of segment in the mixture phase is derived in which interaction energies between segments are calculated from COSMO-RS.The chemical potential of each molecule is then obtained by summing the contributions of the segments.For the prediction, a complete dissociation of ionic liquids into cations and anions was assumed as in our earlier work [13][14][15].The activity coefficient of any component in the mixture is then given by Here ln Γ S (σ) and ln Γ i (σ) refer to the activity coefficient of the segment in the mixture and compound i, respectively.p i (σ) refers to the sigma profile of the component i, and n i implies the contribution of number of segments for component i. Thereafter the Staverman-Guggenheim term for the combinatorial contribution [ln γ SG i/S ] is calculated using the following equation: where Here r i and q i are normalized using the volume (66.69A 2 ) and surface area (79.53A 2 ) of a functional group, that is, r i = COSMO Volume/66.69,q i = COSMO Area/79.53.The COSMO volume and COSMO area are taken from the output of the COSMO file.COSMO-RS thereafter estimates the activity coefficients of the components in the liquid-liquid system according to (2).  and benchmarked at different temperatures (Tables 2 and 3).

Results and Discussion
Besides the effect of temperature and experimental method, the presence of trace amount of impurities such as water or ions provides a remarkable effect on thermodynamic properties [20].Table 2 shows the comparison of measured densities for pure [EMIM][SCN] with temperature.Table 3 shows the comparison between experimental and literature data of the pure aromatic nitrogen/sulphur and water at 298.15 K.This agrees will with data reported for pyridine, pyrrole, and quinoline [21].This also matches with the experimental data obtained earlier for thiophene [22] and water [23].It was observed that the density ρ of pure ionic liquid is higher than that of water and also than that of aromatic nitrogen/sulphur compounds.Table 3 shows a comparison for the surface tension σ and refractive index nD of pure components with literature data at 298.15 K.The surface tension and refractive index were compared successfully against the known value of water [21] and thiophene [22], respectively.The difference between experimental and literature data for pure [EMIM][SCN] can be explained by the thermodynamic properties which strongly depend on water and other impurities.Experimental densities over the entire composition range for  with increasing temperature.The value of α = 0.117 K −1 (ionic liquid) was calculated from

Thermal Expansion "α". The coefficient of thermal expansion of the [EMIM][SCN] α is shown in
where V is the molar volume of the pure fluid, ρ is the density of the pure fluid, and subscript P indicates constant pressure.The volume expansivity of [EMIM][SCN] was calculated using the measured density observed at different temperatures.The volume expansivity of [EMIM][SCN] was found to be independent of temperature (Figure 1).These observations are more consistent with the studies of Rodríguez and Brennecke [24], Wandschneider et al. [25], and Pereiro et al. [26,27].

Effect of Composition on Density and Excess Molar
Volume.[18] and anion providing better interaction with aromatic compounds.A small size of IL molecule has a significant influence on the changing densities of the mixture with an increase in temperature over the studied mole fraction (0.2 to 0.8) of IL.Besides the hydrogen bonding strongly depends on the temperature [30] as compared to CH-π bond interaction [12][13][14][31][32][33] and π-π bond interactions [13,14,[27][28][29].These similar trends were observed for other binary systems [24,28,29].
The measured densities of binary mixtures over the studied mole fraction (0.  x x x x x x
Geng et al. [20] had earlier reported that the ratio of van der Waals and electrostatic force plays a significant role in the extraction process.It was found that mixtures of IL with similar or dissimilar structure of molecules for mole fraction 0.3 and above plays a dominating role at temperatures 293.15 + thiophene shows a similar trend at above 0.3 mole fraction of IL in the mixture.

Effect of Composition on Surface Tension and Refractive
Index.The measured surface tension data for all the studied systems are listed in Table 5 and illustrated in Figure 3.At

tension σ and deviation of surface tension for the binary systems: (a) [EMIM][SCN] + Pyridine, (b) [EMIM][SCN] + Pyrrole, (c) [EMIM][SCN] + Indoline, (d) [EMIM][SCN] + Quinoline, (e) [EMIM][SCN] + Thiophene, and (f) [EMIM]
[SCN] + Water at 298.15 K.   3 that the surface tension increases in the order of thiophene > quinoline > pyridine > indoline > pyrrole > water.Thus an increasing mole fraction of [EMIM][SCN] causes a decrease in the strength of sterical hindrance due to dissimilar molecules or addition of benzene ring in pyrrole or pyridine structure The van der Waals force between the ethyl substitution on the [EMIM] and the aromatic structure of nitrogen/sulphur is not much significant in liquid phase.Further when the similar or dissimilar aromatic structures of molecules are mixed together, then the interaction between [EMIM] cation and [SCN] anion decreases while the CH-π bond interaction, π-π stacking, and dispersion forces get increases particularly for the systems having similar or dissimilar aromatic structure of molecules [37].
Figure 4 shows the composition dependence of refractive index for the above-mentioned binary systems (seen in Table 6).

Effect of Composition on Excess Molar
Volume.The excess molar volume can be defined by [38,39] where V is the molar volume of the pure fluid, ρ is the density of the pure fluid, and subscript P indicates constant pressure.
V E m is the excess molar volume of the binary mixture, ρ mix is density of the mixture, and density of the pure components is denoted as ρ 1 and ρ 2 , respectively.x 1 and x 2 are the mole fraction of components 1 and 2. M 1 , M 2 are the molecular weight of components 1 and 2, respectively.
The calculated excess molar volume values are listed in Tables 4(a) to 4(f).[EMIM][SCN](1) + pyridine(2) (Figure 5) binary mixture exhibits a linear decrease with an increase in the mole fraction of [EMIM][SCN](1) at all studied temperatures.At temperature T = 298.15K a negative excess molar volume for all the systems is obtained within 0.2 to 0.8 mole fraction.This is due to fact that the specific interactions include the partial charge transfer, formation of hydrogen bond, CH-π bond, π-π stacking, and n-π bond interaction.The weak H-bonds are easily broken up when the temperature increases as compared to CH-π bond, π-π stacking, and n-π bond.This indicates that the IL-sulphur/nitrogen mixture has stronger interaction as compared to IL-IL, thiophene-thiophene, or pyrrole-pyrrole interaction.A similar phenomenon is applicable for all other studied systems.This can also be explained by the stronger association occurring between similar or dissimilar aromatic structure of molecules in liquid phase at any composition or temperatures without any temperature and composition influence.
In case of [EMIM][SCN](1) + water(2) mixture, a sudden decrease is seen at initial mole fraction, thereafter decreasing linearly as the composition increases.This is due to the fact that at individual concentration (higher dilution rate), the oxygen sites are stabilized by electron or partial charge transfer.Thus [EMIM][SCN] has strong attraction with water at low concentrations.Furthermore the temperature effect on the excess molar volume is negligible since the rate of kinetic energy is not enough to change the phase behavior upon mixing with aromatic structure of the nitrogen/sulphur molecules.It means that the structural effect becomes more important than the dispersion or dipolar interactions.Hence the composition has significant influence on excess molar volume of aromatic structure with multiple ring compounds like indoline and quinoline.

Effect of Composition on Deviation of Surface Tension
and Refractive Index.The deviation of surface tension and refractive index can be defined by [35][36][37][38][39] where Δσ and ΔnD are the deviation of surface tension and refractive index.The surface tension and refractive index of the binary mixture are denoted as σ mix and nD mix , respectively.x 1 and x 2 are the mole fraction of the components 1 and 2. The deviation of surface tension and refractive index are presented in Figures 6 and 7.The deviation of surface tension for [EMIM][SCN](1) + pyrrole(2) indicates an increasing trend with increasing mole fraction of [EMIM][SCN] while the refractive index deviation shows a decreasing trend at T = 298.15K.It can be observed that the strong electrostatic interaction takes place between the two similar aromatic structures of the molecules (IL-nitrogen/sulphur compound), even though molecules posses additional benzene ring such as indoline and quinoline.The deviation of surface tension was found to be inversely proportional to deviation of refractive index.This is due to the fact that the deviation of surface tension of the binary mixture strongly depends on the cohesive interaction at the surface level while the refractive index deviation depends on the dispersion interaction upon mixing.However both depend on the size and shape of the molecules in a liquid phase.These tendencies agree well with the reported data of Nain [40] and Freire et al. [41].In [EMIM][SCN] + Water mixture, the deviation of surface tension increases linearly with composition (Figure 7) while the deviation of refractive index decreases slightly (Figure 8).The molecular interaction is very strong between similar and dissimilar aromatic structures [31,34,42,43] of sulphur/nitrogen compounds and [EMIM] cation while the mobility of ions is weakened upon mixing.However [EMIM][SCN] interacts with the compounds in liquid phase via n-π bonding (where n: N,O from both the side).It can be observed that the minimum amount of [EMIM][SCN] ionic liquid is enough for separation of similar structure of pyrrole and thiophene while for pyridine, indoline, and quinoline more amount is required as compared to pyrrole and thiophene.Negative deviation was observed for the mole fraction (0.2 to 0.8) of IL at temperature from 298.15 to 323.15 K for both deviation in surface tension and refractive index studies.The observed trends are as follows.

Relation of Activity Coefficient with Excess Molar Volume.
In order to substantiate the findings we have predicted the activity coefficients of solutes (nitrogen and sulphur compounds) in mixture (Tables 4(a) to 4(f)).It is interesting to know that that pyrrole (Table 4(b)) and water (Table 4(c)) have activity coefficients less than one among all the studied systems.It implies that water is soluble in IL, thus the IL is hydrophilic in nature.Among all the compounds it seems that pyrrole is the most soluble.Figure 8 shows a comparison of predicted activity coefficient and excess volume for ILpyridine mixture.The predicted activity coefficient increases with increasing mole fraction of ionic liquid in the mixture, except for [EMIM][SCN] + water systems.However the excess molar volumes are seen to decrease with mole fraction.Thus it is fair to say that the activity coefficient is inversely proportional to excess volume.
It is observed that the properties of mixture are influenced by the aromatic ring structure of thiophene, pyridine, pyrrole, indoline, and quinoline (Tables 4(a) to 4(f)).The activity coefficient also gives us an indication of the strength of solubility.For example, [EMIM][SCN](1) + Pyridine(2) system (Figure 8) shows low activity coefficient (i.e., high solubility with mole fraction of IL) as compared to [EMIM][SCN] + thiophene (Table 4(e)) which implies that the solubility of pyridine in [EMIM][SCN] is higher than thiophene.In case of excess molar volume both the systems show an increasing excess molar volume with increasing mole fraction of IL.

Conclusion
Excess molar volume and deviation of surface tension and refractive index have been calculated from measured densities, surface tension, and refractive index for the systems of [EMIM][SCN] ionic liquid with pyrrole, pyridine, indoline, quinoline, thaiophene, and water in a liquid phase at temperature of 298.15 to 323.15 K and atmospheric pressure.The surface tension increases with increasing mole fraction in the order of thiophene > quinoline > pyridine > indoline > pyrrole > water, while the refractive index increases in the order of [EMIM][SCN]: pyrrole < Pyridine < thiophene < indoline <quinoline.Furthermore the deviation of surface tension and refractive index show a similar negative trend which follows: for surface tension: Water > quinoline > thiophene > pyrrole > pyridine >indoline, for refractive index: Indoline > thiophene > pyrrole > pyridine > quinoline > water.Finally it was concluded that the similar aromatic structures of nitrogen and/or sulphur compounds are removed via CH-π bond interaction and π-π stacking.The aromatic heterocycle with multiple ring of nitrogen and sulphur compounds removal was found to strongly depend on the composition of IL as compared to temperature effect.It is also observed that the activity coefficient of IL increases while the excess molar volumes decreases with increasing mole fraction of IL.

− 1 )Figure 3 :
Figure 3: Composition dependence of surface tension of the binary mixture plotted against mole fraction of IL at 298.15 K.

Figure 4 :
Figure 4: Composition dependence of surface tension of the binary mixture plotted against mole fraction of IL at 298.15 K.

Figure 7 :
Figure 7: Refractive index deviation against mole fraction of IL at 298.15 K.

Table 1 :
Name and structure of the studied ionic liquid.

Table 2 :
Comparison of measured density ρ as a function of temperature with literature data for pure compounds.

Table 3 :
Surface tension σ and refractive index nD for pure compounds of the binary system studied in this work at 298.15K and atmospheric pressure.

Table 2
shows the experimental densities of pure [EMIM][SCN], aromatic nitrogen, aromatic sulphur, and water at temperatures of T = (298.15to 323.15) K.It can be seen that the density of all the pure components including ionic liquid decreases linearly with increasing temperature.Tables 4(a) to 4(f) show the experimental densities for six binary systems studied as a function of temperature over the studied mole fraction (0.2 to 0.8) of [EMIM][SCN].