Journal of Chemistry

Volume 2013 (2013), Article ID 838590, 7 pages

http://dx.doi.org/10.1155/2013/838590

Research Article

Investigation on Aromaticity Index and Double-Bond Equivalent of Aromatic Compounds and Ionic Liquids for Fuel Desulphurization

Syamsul B. Abdullah,¹ Z. Man,² and M. A. Bustam²

¹Faculty of Chemical and Natural Resources Engineering, Universiti Malaysia Pahang, Lebuhraya Tun Razak, Gambang, 26300 Pahang, Malaysia

²Chemical Engineering Department, PETRONAS Ionic Liquid Centre, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, Tronoh, 31750 Perak, Malaysia

Received 18 July 2013; Revised 5 November 2013; Accepted 15 November 2013

Academic Editor: Yoshihiro Kudo

Copyright © 2013 Syamsul B. Abdullah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Aromaticity index (AI) and double-bond equivalent (DBE) were studied to seek a simpler approach in identifying suitable ionic liquids for the desulphurization process. Manifestation of interaction mechanism by COSMO-RS and Raman spectroscopy is discussed. The findings demonstrate that AI and DBE calculations can be used to analyze and arrogate potential combination of cation and anion for fuel desulphurization purposes.

1. Introduction

Due to the limited efficiency of hydrodesulphurization towards removing aromatic organosulphur compounds, a number of recent researches had focused on extracting them using ionic liquids [1, 2]. Notable characteristics of ionic liquids (ILs) are their nonmeasurable vapour pressure, thermal stability with wide liquid range, and solvating properties for diverse materials. Because of these notable characteristics, ILs have drawn great attention in which their properties are believed to leave minimal impact on the environment as well as human health [3–9].

The prospect of ILs in becoming an industrial solvent for the desulphurization process may be illustrated from the solute-solvent properties’ point of view [7, 9–13]. These may include the polarity of the solute, the dipole moment, molecular size, and molar volume. As for the aromatic organosulfur functioning as a solute, its miscibility in ILs can mostly be described through its -electron density. The solute with a strong -electron density would be highly miscible in ILs and vice versa. This leads to our observation that -electron density could be a vital point in screening a potential IL for desulphurization [1, 11–14]. However, to determine the -electron density is not simple whereby it needs to be calculated using the quantum chemical calculations [13, 14]. Thus, the quantification of the -electron density for ILs using simple calculation needs immediate justification. Therefore, here, by suggesting a new complimentary method that estimates the aromaticity index and double-bond equivalent which are closely related to -electron density, the possibility of using it to preliminarily select an appropriate IL for fuel desulphurization purposes shall be explored.

2. Experimental and Theory

From Conductor-like Screening Model for Realistic Solvents (COSMO-RS) database, 407 combinations of organic ILs (11 of cations and 37 of anions) were investigated, while metal-based ILs were excluded. First, three important parameters that determine the extractive ability of the ILs which are predicted using COSMO-RS were basically calculated from the estimated activity coefficient () at infinite dilution. These three parameters are capacity (), selectivity (), and performance index (PI). Capacity is defined as the maximum amount of solute that can be dissolved in ILs, whereby selectivity is simplified from the ratio of the composition of solute in the extract (ILs phase) and its composition in raffinate (model fuel phase), while the resultant of multiplication between capacity and selectivity is defined as performance index. Numerically, these three parameters can be expressed by (1), (2), and (3), respectively [13]. These parameters have been estimated for a preliminary data comparison study against reported desulphurization performances. The three parameters can be expressed as follows:

Koch and Dittmar introduced an AI that can be calculated solely from the exact molecular mass of the organic compound. AI is a measure of the C–C double-bond “density” which additionally considers the contribution of -bonds by heteroatoms. Since all heteroatoms are taken into account as potential contributors, the actual aromaticity in a molecule can be higher and AI must be regarded as the most conservative approach, according to (4) [15]:

The number of each atom in the molecule, for example, carbon (C), oxygen (O), sulphur (S), hydrogen (H), nitrogen (N), and phosphorus (P), provided the value of AI. The atoms which are described in (4) are the most common elements in an organic compound [15].

DBE has been introduced through organic matters characterization studies and its calculation is a well-established tool in mass spectrometry. In short, DBE can be defined as a summation of the unsaturations plus rings in a molecule which basically can be calculated via a molecular formula, using (5), where and are the number of atom and valence of each atom, respectively [15]:

In order to identify the interaction mechanism, a vial containing a sample of BT/[BMIM][OSO₄] mixture was placed on the top of the laser spot. The total amount of the sample mixture was fixed accordingly for interaction purposes between benzothiophene (BT) and 1-butyl-3-methylimidazolium octylsulfate ([BMIM][OSO₄]). A series of BT/[BMIM][OSO₄] binary system with mole fraction in the range of , and 0.72 was studied.

The investigated ILs were synthesized in-house and some were purchased from Merck with stated halides content of less than 1%. Allowable water content in ILs for extraction has to be maintained below 500 ppm. In the present extractive desulphurization process, a mixture of -dodecane (-C₁₂) and BT (1000 µg/mL) was used as model fuel. The extraction was carried out in 25 mL screw-cap vials. The mixture was vigorously stirred at room temperature for 10 min and another 10 min for settling to attain equilibrium, which has been previously identified [2]. The ILs have higher densities from the model fuel. The upper phase which is the model fuel could be withdrawn easily and analyzed for its sulphur content using the HPLC (DAD detector, reverser-phase ZORBAX SB-C18 column), and, based on the values, a material balance was calculated. At equilibrium, the sulphur partition coefficient determines the extent of the mobility of sulphur into ILs phase. The sulphur partition coefficient is defined as the ratio of equilibrium concentration of sulphur in ILs phase () to that in model fuel phase (), as shown in (6) below. The details of the procedures and analysis are provided elsewhere [2]

3. Result and Discussion

In the electroneutral mixture method of the COSMO-RS, in order to estimate the activity coefficient (), a mixture of ILs and solute was treated as ternary system (1 : 1 IL, cation + anion + solute) that can be converted into a binary system. The activity coefficient value is estimated for calculating the capacity and selectivity values. Later, the prediction extractability or rather known as performance index of the ILs at infinite dilution was determined using calculated capacity and selectivity values. Meanwhile, the experimental extractability which was collected from the literature is plotted against performance index for interaction mechanism identification, as shown in Figure 1. In general, it can be seen that the experimental extractability of sulphur compounds which are extracted from -C₁₂ by [BMIM][OSO₄] decreased in the order of DBT > BT > 4-MDBT > 4,6-DMDBT > TS > tetrahydrothiophene > dodecanethiol while, on the other hand, the prediction extractability by COSMO-RS depicted a different order with TS > tetrahydrothiophene > BT > DBT > 4-MDBT > 4,6-DMDBT > dodecanethiol. In the experimental approach, the researchers have already explained that the order of sulphur compounds extracted mostly depended upon the -electron density of the compounds but, for the prediction approach via COSMO-RS, the explanations came from the sigma profile plot, shown in Figure 2. On the left-hand side of the histogram is the hydrogen bond donor region (;  e/Å²) and on the right-hand side is the hydrogen acceptor region (;  e/Å²). From this figure, it can be seen that the cation part of [BMIM][OSO₄] lies in the donor region while the anion part lies in the acceptor region. A small fraction of the sigma profile of all predicted sulphur compounds stretched out into the donor and acceptor region, thus indicating that these sulphur compounds favour a very weak hydrogen bonding. The overlapping of the sigma profile between the sulphur compounds and [BMIM][OSO₄] represents that miscibility between these compounds is possible. However, thiophene (TS) showed more overlaps compared to the other sulphur compounds which proves that more thiophene will be extracted from -C₁₂ as per estimated. With these inconsistent observations between experimental and prediction approach by COSMO-RS, we were determined to investigate the interaction mechanism between BT molecules and [BMIM][OSO₄] using Raman spectroscopy.

Figure 1: Comparison of extractability ( value) of the sulfur species from -C₁₂ collected from the literature [8, 10] and COSMO-RS prediction.

Figure 2: Sigma profile for [BMIM][OSO₄] and respective sulfur species.

Raman spectroscopy is a vibrational molecular spectroscopy which derives from an inelastic light scattering process. In Raman spectroscopy, sample preparation is simpler as compared to other Fourier transforms spectroscopy techniques. Basically, with Raman spectroscopy, a laser photon is scattered by a sample molecule and losses or gains of energy during the process. The amount of energy lost is seen as a change in energy (wavelength) of the irradiating photon. This energy loss is characteristic for a particular bond in the molecule. Figure 3 shows the comparison in Raman spectra in the 1600–1400 cm⁻¹ region where shifting occurred for two identified peaks at 1560 and 1420 cm⁻¹, respectively. By observing these peaks, the shifted peaks have been assigned to stretching and the C–H vibrations of the aromatic ring system were identified. It is probable that the interaction largely involved the -bond between BT and C–H of imidazolium cation of [BMIM][OSO₄]. Irrevocably, it proved that -electron density and aromaticity played a crucial role in extracting aromatic organosulphur which ruled out the conceptual observation explicated by the sigma profile of COSMO-RS. While looking at this mechanism which may fit any kind of aromatic compounds, we were driven to quantify potential ILs for fuel desulphurization by -electron and aromaticity aspect using AI and DBE calculation approach, instead of relying on -electron density calculations which are hard to be estimated and are complicated calculations. But first, preassessment is needed in order to justify this screening method.

Figure 3: Raman spectra of over the range of 1600–1400 cm⁻¹ with the mole fraction increasing from 0.22 to 0.72.

The AI values of five aromatic organosulphurs (TS, BT, DBT, 4-MDBT, and 4,6-DMDBT) and three aromatic hydrocarbons (benzene, toluene, and xylene) were calculated using (4) and illustrated in Figure 4. With the aim of looking for some pattern within each aromatic species, the experimental extractability has been added for comparison study. The AI values of eight studied compounds were found to be relatively proportional with the experimental extractability whereby they followed the order of DBT > BT > TS > 4-MDBT > benzene > 4,6-DMDBT > toluene > xylene. For example, based on the miscibility study of aromatic hydrocarbons alone, as collected from the literature, the most miscible compound in imidazolium-based ILs is benzene followed by toluene and xylene. Apart from considering the molecular size or alkyl chain length, there is an alternative explanation in describing the miscibility of the aromatic hydrocarbons in imidazolium-based ILs: miscibility decreases as the value of calculated AI decreases. This fact indicates that calculating AI value could become a prediction tool for miscibility study of aromatic compounds in imidazolium-based ILs. Due to this, it is particularly exciting that the relationship between AI and the experimental extractability of aromatic compounds may open the possibility for guesstimating and perhaps screening potential ILs for the desulphurization process.

Figure 4: Comparison of extractability ( value) organosulphur and aromatic hydrocarbon compounds collected from the literature [8, 10–12] and AI value.

Table 1 shows the calculated AI and DBE values for common cation and anion. The AI values were calculated based on the root structure of cation alone by ignoring the alkyl chain attached prior to determine the pattern toward experimental extractability for extractive desulphurization process. Five ILs have been synthesized accordingly (refer to Table 2) and their desulphurization performance was determined in terms of sulfur partition coefficient ( value). Basically, it is apparent from the pattern of AI values that the order was as follows: imidazole ≈ pyrazole > benzimidazole > pyridine > pyrrolidone. Meanwhile, the results from extractive desulphurization depicted the order differently, as shown in Figure 5. The inconsistent order between AI and value was found to be due to the physical phase of ILs whereby benzimidazole-based ILs are in the solid form at room temperature, resulting in poor desulphurization performance ( value =0.124). As part of our study on screening the ILs for desulphurization, physical phase is an important aspect that should be considered first before evaluating AI value for cation selection.

Table 1: Calculated aromaticity index (AI) and double-bond equivalent (DBE) value for common cation and anion.

Cation (basic structure)
Anion

Chemical formula
AI
Chemical formula
DBE

Benzimidazole, C₇H₆N₂
1
Acetate, C₂H₃O₂
1.5

Benzyltriphenylphosphine,
0.6458
Benzoate, C₇H₅O₂
5.5

C₂₅H₂₁P

Bis(tpM)phos., C₁₆H₃₄O₂P
0.5

Guanidine, CH₅N₃
—
Butylsulfate,C₄H₉O₄S
0.5

Imidazole, C₃H₄N₂
2
Decanoate, C₁₀H₁₉O₂
1.5

Oxazole, C₃H₃NO
1.5
DBP, C₈H₁₈O₄P
0.5

Piperidine, C₄H₁₀N₂
0
Dicyanamide, C₂N₃
4.5

Pyrazole, C₃H₄N₂
2
DEP, C₄H₁₀O₄P
0.5

Pyridine, C₅H₅N
0.875
DHP, H₂O₄P
0.5

Pyrrolidone, C₅H₁₁N
0.125
DMP, C₂H₆O₄P
0.5

Quinoline, C₉H₇N
1.083
EESO₄, C₄H₉O₅S
0.5

Thiazole, C₃H₃NS
1.5
ESO₄, C₂H₅O₄S
0.5

Hydrogensulfate,HSO₄
0.5

Methanesulfonate, CH₃O₃S
0.5

MESO₄, C₃H₇O₅S
0.5

MSO₄, CH₃O₄S
0.5

Nitrate, NO₃
1.5

OSO₄, C₈H₁₇O₄S
0.5

Salicylate, C₇H₅O₃
5.5

Thiocyanate, CNS
2.5

TOS, C₇H₇O₃S
4.5

TCM, C₄N₃
6.5

(ME)ESO₄, C₅H₁₁O₆S
0.5

Imidazolide or pyrazolide, C₃H₃N₂
3.5

Bis(PFE)phos., C₄F₁₀O₂P
0.5

Bis(PFES)ami., C₄F₁₀NO₄S₂
0.5

Bis(TFM)imi., C₂F₆N
0.5

Bis(TFMS)meth., C₃HF₆O₄S₂
0.5

Perchlorate, ClO₄
0.5

HFB, C₄F₇O₂
0.5

IOPF, F₄NO₂P₂
0.5

pFBS, C₄F₉O₃S
0.5

NTf₂, C₂F₆NO₄S₂
0.5

TFA, C₂F₃O₂
1.5

OTf, CF₃O₃S
0.5

Tris(tFMS)meth., C₄F₉O₆S₃
0.5

Table 2: Selected combination of cation and anion for ionic liquids with their respective calculated AI and DBE value.

Cation
Anion
value

1, 3-Dimethylimidazolium, AI = 2
MSO₄, DBE = 0.5
0.409

1, 2-Dimethylpyrazolium, AI = 2
MSO₄, DBE = 0.5
0.366

1, 3-Dimethylbenzimidazolium, AI = 1
MSO₄, DBE = 0.5
0.124

1, 4-Dimethylpyridinium, AI = 0.875
MSO₄, DBE = 0.5
0.334

1, 1-Dimethylpyrrolidinium, AI = 0.125
MSO₄, DBE = 0.5
0.306

1-Butyl-3-methylimidazolium, AI = 2
DHP, DBE = 0.5
0.191

1-Butyl-3-methylimidazolium, AI = 2
DBP, DBE = 0.5
1.703

1-Butyl-3-methylimidazolium, AI = 2
HSO₄, DBE = 0.5
0.205

1-Butyl-3-methylimidazolium, AI = 2
MSO₄, DBE = 0.5
1.222

1-Butyl-3-methylimidazolium, AI = 2
OSO₄, DBE = 0.5
2.030

1-Butyl-3-methylimidazolium, AI = 2
OTf, DBE = 0.5
0.942

1-Butyl-3-methylimidazolium, AI = 2
Acetate, DBE = 1.5
1.564

1-Butyl-3-methylimidazolium, AI = 2
NO₃, DBE = 1.5
1.500

1-Butyl-3-methylimidazolium, AI = 2
CNS, DBE = 2.5
1.857

1-Butyl-3-methylimidazolium, AI = 2
Imidazolide, DBE = 3.5
2.125

1-Butyl-3-methylimidazolium, AI = 2
Pyrazolide, DBE = 3.5
2.030

1-Butyl-3-methylimidazolium, AI = 2
TOS, DBE = 4.5
0.471

1-Butyl-3-methylimidazolium, AI = 2
Dicyanamide, DBE = 4.5
2.448

1-Butyl-3-methylimidazolium, AI = 2
Benzoate, DBE = 5.5
2.846

1-Butyl-3-methylimidazolium, AI = 2
Salicylate, DBE = 5.5
2.571

1-Butyl-3-methylimidazolium, AI = 2
TCM, DBE = 6.5
3.762

Figure 5: Comparison of AI value and desulphurization partition coefficient ( value) of cations solely with MSO₄.

In the second part of predicting ILs for extractive desulphurization, it should be noted that all of the ILs used are composed of relatively similar cations. Table 2 reports the value of 16 ILs at room temperature of experimental conditions along with their respective calculated DBE values, which are also illustrated in Figure 6. DBE value was investigated instead of AI due to the limited number of anions that are currently available. The DBE value agrees reasonably well with the experimental findings if 1-butyl-3-methylimidazolium tosylate (TOS) is excluded. The solid phase of this IL attributes lavishly in decreasing the desulphurization ability on model fuel. Within the same root of two types of anions, namely, phosphate-based ILs (referred to as dihydrogenphosphate (DHP) and dibutylphosphate (DBP)) and sulfate-based ILs (referred to as hydrogensulfate (HSO₄), methylsulfate (MSO₄), and octylsulfate (OSO₄)), anions that have longer alkyl chains gave higher performances. Weakening of anions due to alkyl chain length can be correlated to the strength of cation-anion interactions whereby the shorter the alkyl chain, the stronger the interaction between cation and anion. Mochizuki and Sugawara [1] have observed a similar trend in their work where the extraction of sulphur compound increases linearly with the increase in alkyl chain length of imidazolium-based sulfate ILs. Among the ILs, tricyanomethane (TCM) stood out with a value of 3.762. This is closely followed by benzoate and salicylate anion with values of 2.846 and 2.571, respectively. This proves that, by determining AI and DBE values, the appropriate ILs for fuel desulphurization can be classified.

Figure 6: Comparison of DBE value and desulphurization partition coefficient ( value) of anions solely with BMIM.

4. Conclusion

Systematic investigations on AI and DBE of ILs have been presented. The predicted values obtained using COSMO-RS did not converge well with values obtained from the literature due to the interaction theory consideration. Using Raman spectroscopy, interaction mechanism has been revealed where the extractability of BT by [BMIM] [OSO₄] is through aromaticity and -electron density. Further study then confirmed that the extractability of aromatic compounds is proportionate to -electron which was closely related to AI and DBE. In short, directly simple AI and DBE calculations can be used to analyze and arrogate potential combination of cation and anion ILs for fuel desulphurization purposes.

Abbreviations

AI: | Aromaticity index
BT: | Benzothiophene
Bz: | Benzene
COSMO-RS: | Conductor-like Screening Model for Realistic Solvents
DBT: | Dibenzothiophene
DBE: | Double bond equivalent
IL: | Ionic liquid
mmBz: | 1,3-Dimethylbenzimidazolium
mmim: | 1,3-Dimethylimidazolium
mmPy: | 1,4-Dimethylpyridinium
mmPyrr: | 1,1-Dimethylpyrrolidinium
mmPyz: | 1,2-Dimethylpyrazolium
Tl: | Toluene
TS: | Thiophene
Xy: | Xylene
: | 1-Butyl-3-methylimidazolium octylsulfate
4-MDBT: | 4-Methyldibenzothiophene
4,6-DMDBT: | 4,6-Dimethyldibenzothiophene.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

Y. Mochizuki and K. Sugawara, “Removal of organic sulfur from hydrocarbon resources using ionic liquids,” Energy and Fuels, vol. 22, no. 5, pp. 3303–3307, 2008.

B. A. Syamsul, Z. Man, L. Ismail, A. Maulud, M. I. A. Mutalib, and M. A. Bustam, “An optimization study on the extraction of benzothiophene from model fuel using the ionic liquids 1-butyl-3-methylimidazolium octylsulfate,” Research Journal of Chemistry and Environment, vol. 15, no. 2, pp. 510–518, 2011.

A. Marciniak, “Influence of anion structure on the liquid-liquid equilibria of 1-ethyl-3-methyl-imidazolium cation based ionic liquid-hydrocarbon binary systems,” Journal of Chemical and Engineering Data, vol. 56, no. 3, pp. 368–374, 2011.

M. B. Shiflett and A. M. S. Niehaus, “Liquid-liquid equilibria in binary mixtures containing substituted benzenes with ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide,” Journal of Chemical and Engineering Data, vol. 55, no. 1, pp. 346–353, 2010.

U. Domańska and M. Królikowski, “Determination of activity coefficients at infinite dilution of 35 solutes in the ionic liquid, 1-butyl-3-methylimidazolium tosylate, using gas-liquid chromatography,” Journal of Chemical and Engineering Data, vol. 55, no. 11, pp. 4817–4822, 2010.

M. B. Shiflett and A. Yokozeki, “Liquid-liquid equilibria in binary mixtures containing fluorinated benzenes and ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide,” Journal of Chemical and Engineering Data, vol. 53, no. 11, pp. 2683–2691, 2008.

Y. Huo, S. Xia, and P. Ma, “Solubility of alcohols and aromatic compounds in imidazolium-based ionic liquids,” Journal of Chemical and Engineering Data, vol. 53, no. 11, pp. 2535–2539, 2008.

C. C. Cassol, A. P. Umpierre, G. Ebeling, B. Ferrera, S. S. X. Chiaro, and J. Dupont, “On the extraction of aromatic compounds from hydrocarbons by imidazolium ionic liquids,” International Journal of Molecular Sciences, vol. 8, no. 7, pp. 593–605, 2007.

A. Arce, M. J. Earle, H. Rodríguez, and K. R. Seddon, “Separation of benzene and hexane by solvent extraction with 1-alkyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}amide ionic liquids: effect of the alkyl-substituent length,” Journal of Physical Chemistry B, vol. 111, no. 18, pp. 4732–4736, 2007.

J. Eßer, P. Wasserscheid, and A. Jess, “Deep desulfurization of oil refinery streams by extraction with ionic liquids,” Green Chemistry, vol. 6, no. 7, pp. 316–322, 2004.

Y. Nie, C. Li, A. Sun, H. Meng, and Z. Wang, “Extractive desulfurization of gasoline using imidazolium-based phosphoric ionic liquids,” Energy and Fuels, vol. 20, no. 5, pp. 2083–2087, 2006.

X. Jiang, Y. Nie, C. Li, and Z. Wang, “Imidazolium-based alkylphosphate ionic liquids—a potential solvent for extractive desulfurization of fuel,” Fuel, vol. 87, no. 1, pp. 79–84, 2008.

A. A. P. Kumar and T. Banerjee, “Thiophene separation with ionic liquids for desulphurization: a quantum chemical approach,” Fluid Phase Equilibria, vol. 278, no. 1-2, pp. 1–8, 2009.

B.-M. Su, S. Zhang, and Z. C. Zhang, “Structural elucidation of thiophene interaction with ionic liquids by multinuclear NMR spectroscopy,” Journal of Physical Chemistry B, vol. 108, no. 50, pp. 19510–19517, 2004.

B. P. Koch and T. Dittmar, “From mass to structure: an aromaticity index for high-resolution mass data of natural organic matter,” Rapid Communications in Mass Spectrometry, vol. 20, no. 5, pp. 926–932, 2006.