the online version of this article ( doi:10.1007/s10953 - 014 - 0276-y ) contains supplementary material , which is available to authorized users .
in recent years , the deep desulfurization of diesel fuel has become the most studied process with different techniques ( extraction , liquid  liquid separation , oxidative desulfurization , adsorption ) .
the emission of sulfur from petrol and diesel oils , which is linked to acid rain , plays a crucial role in pollution problems of large conglomerates .
thus , the usa and european countries have issued regulations regarding sulfur content in fuels [ 1 , 2 ] . due to this situation ,
the european union approved a new directive stating that the content of total sulfur in european gasoline and diesel fuels from 2010 onwards must be at a maximum concentration level of 10  ppm .
ionic liquids ( ils ) have the ability to extract aromatic sulfur - containing compounds at ambient conditions .
additionally , ils are immiscible with the fuel , are non - volatile and can be regenerated and recycled by solvent washing .
oxidative desulfurization in future years probably will bring better results than simple liquid  liquid separation , however , first the best ils must be chosen . at present , the hydrodesulfurization ( hds ) processes is the established method used in some industrial technologies to remove organic sulfur from fuels .
however , to achieve low sulfur targets with current hds technology , higher temperatures , higher pressures , larger reactor volumes , and more active catalysts are needed .
the hds process does not purify fuels of polycyclic organic sulfides such as thiophene , benzothiophene , methyldibenzothiophenes , 4,6-dibenzothiophenethiols , thioethers , and disulfides .
extraction desulfurization , which has begun to be popular , especially with ils , has the potential for being an alternative and future complementary technology for deep desulfurization [ 410 ] . in order to solve this problem , extractive liquid
 liquid equilibrium ( lle ) desulfurization with ils has been proposed [ 7 , 1018 ] .
the 1-alkylpiperidinium - based , or pyrrolidinium - based ils with different anions , or 1-alkylcyanopyridinium - based ils , have been recently studied in our laboratory in ternary lle { il  +  thiophene , or benzothiophene  +  heptane ) with high selectivities .
attractive extraction parameters were presented as well for 1-ethyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}imide , [ emim][ntf2 ] ( , and references cited therein ) , 1-ethyl-3-methylimidazolium acetate , [ emim][oac ] , 1-ethyl-3-methylimidazolium thiocyanate , [ emim][scn ] , and 1,3-dimethylimidazolium methylphosphonate [ dmim][mp ] .
this work is a continuation of our systematic studies on the physicochemical properties and the extraction abilities of piperidinium - based ils ( and references cited therein ) .
proposed by us are new interaction parameters for the group contribution method modified unifac for the piperidnium - based ils , predicted attractive infinite dilution selectivity , and capacity of piperidinium - based ils ( alkane chain , n  =  36 ) in the thiophene / heptane separation problem at t  =  328.15  k. to our best knowledge , the phosphonium - based il ( tributyl - methylphosphonium methylsulfate , [ p1,4,4,4][ch3so4 ] ) was measured in ternary lle for the separation of thiophene from cyclohexane at t  =  298.15  k with very low selectivities in a range of 1.5 to 5.4 .
better results were obtained with deep eutectic solvents , des , containing phosphonium - based ils with ethylene glycol .
des is composed of methyltripentylphosphonium bromide , [ p1,5,5,5][br ] , and ethylene glycol ( as a hydrogen bond donor ) showed selectivities of about s  =  60100 for ternary lle at t  =  318  k for benzene / hexane separation .
poorer results were estimated with des composed of tetrabutylphosphonium bromide , [ p4,4,4,4][br ] , and ethylene glycol for the separation of toluene / heptane .
usually , the results of separation processes for aliphatic / aromatic hydrocarbons provide good information for the separation of aliphatic / aromatic sulfur compounds .
we can expect similar or even better results for the chosen il . on the other side ,
there are very good results obtained with methylphosphonate   and diethylphosphate   anions of ils in ternary lle ( il  +  thiophene  +  heptane ) mixtures . in this work
we report experimental ternary lle data for one additional piperidinium - based il , { 1-pentyl-1-methylpiperidinium bis{(trifluoromethyl)sylfonyl}imide [ c5mpip][ntf2 ] , for comparison with measured earlier 1-propyl , or 1-butyl- , or 1-hexyl-1-methylpiperidinium bis{(trifluoromethyl)sylfonyl}imide .
moreover , tributylethylphosphonium diethylphosphate , [ p2,4,4,4][dep ] was chosen to check the influence of the anion .
the [ dep ] anion in [ emim][dep ] shows interesting results for thiophene extraction from hexane at t  =  298.15  k .
the solvents heptane , thiophene , and benzothiophene used in this work are model compounds for fuel and sulfur organic hydrocarbons , respectively .
the ternary systems { il ( 1 )  +  thiophene , or benzothiophene ( 2 )  +  heptane ( 3 ) } were investigated at t  =  308.15  k and p  =  101.33  kpa . the experimental tie - lines for four ternary mixtures
the solute distribution ratio and the extractive selectivity were determined from the experimental data , and are compared to the literature data .
the ils studied , [ c1c5pip][ntf2 ] and [ p2,4,4,4][dep ] , were purchased from iolitec .
the names , abbreviations , structures , measured densities and mass fraction of ils are listed in table  1 .
the names , cas numbers , sources , mass fraction purities , purification method , water content , and measured and literature densities of all chemicals used are shown in table  1s in the supplementary material information .
the samples of ils were dried for 24  h at 300  k under reduced pressure to remove volatile impurities and trace amounts of water .
thiophene and benzothiophene were stored over freshly activated molecular sieves of type 4   ( union carbide ) .
the densities for all substances were measured at t  =  298.15 and 101.33  kpa . the method and uncertainties
have been described previously .table  1list of investigated ionic liquids : structure , name , abbreviation of name , molar mass ( m ) and densitycationanionname , abbreviationm/(gmol)exp .
density /gcm ( 298.15  k ; 101.33  kpa ) 



 1-pentyl-1-methylpiperidinium bis{(trifluoromethyl ) sulfonyl}imide [ c1c5pip][ntf2]450.461.35016 



 tributylethylphosphonium diethylphosphate , [ p2,4,4,4][dep]384.471.0089 list of investigated ionic liquids : structure , name , abbreviation of name , molar mass ( m ) and density the water content was analyzed by the karl - fischer titration ( method titroline kf ) .
the sample of il , or solvent , was dissolved in methanol and titrated in steps of 0.0025  cm .
the error in the water content is 10    10 in mass fraction for the 3  cm of injected il .
the water content in solvents used was less than 350    10 in mass fraction . to obtain the experimental lle tie - lines , mixtures with compositions inside the immiscible region of the systems
the vessel was tightly closed to avoid losses by evaporation or pickup of moisture from the atmosphere .
the jackets were connected to a thermostatic water bath ( lauda alpha ) to maintain a constant temperature of t  =  308.15  k ( 0.05 ) .
the mixtures were stirred for 6  h to reach thermodynamic equilibrium and after a minimum of 12  h were analyzed .
after the phase separation , samples of about ( 0.10.3 )    10 cm were taken from both phases using glass syringes with coupled stainless steel needles .
a sample of the phase was placed in an ampoule with a capacity of 2    10 cm .
next , acetone ( 1.0  cm ) was added to the samples to avoid phase splitting and to maintain a homogeneous mixture .
because of the low vapor pressure , the ils used in this work can not be analyzed by gc .
thus , only thiophene or benzothiophene and heptane were analyzed ; the mass fraction of the third component , the il , was determined by subtracting the mole fractions of the two other components from unity .
the compositions were analyzed by gas chromatography ( perkinelmer clarus 580 gc equipped with auto sampler and fid and tcd detectors ) .
the capillary column of the chromatograph was protected with a pre - column to avoid the non - volatile ionic liquid reaching the column in the case of a leak from the glass wool in the liner .
the totalchrom workstation software was used to obtain the chromatographic areas for the thiophene , or benzothiophene , heptane and the internal standard propan-1-ol .
details of the operational conditions of the apparatus are reported in table  2s in the supplementary material .
the estimated uncertainty in the determination of mole fraction compositions is    0.003 for compositions of the hydrocarbon - rich phase and 0.005 for compositions of the il - rich phase .
the equilibrium compositions of the experimental tie - line ends in ternary systems of four mixtures { il ( 1 )  +  thiophene or benzothiophene ( 2 )  +  heptane ( 3 ) } , at t  =  308.15  k and p  =  101.33  kpa are reported in table  2 .
experimental solubilities for [ c1c5pip][ntf2 ] and [ p2,4,4,4][dep ] in heptane at t  =  308.15  k are totally different from each other . in the binary { il ( 1 )  +  heptane ( 3 )
} system complete liquid miscibility ( solubility of heptane in the il ) is up to mole fraction of heptane \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{3}^{\text{il } } $ $ \end{document}x3il  =  0.089 and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{3}^{\text{il } } $ $ \end{document}x3il  =  0.456 for [ c1c5pip][ntf2 ] and [ p2,4,4,4][dep ] , respectively .
the solubility of heptane is much larger in [ p2,4,4,4][dep ] than that in [ c1c5pip][ntf2 ] .
the piperidinium - based il shows much lower solubility of heptane in the il . in comparison with piperidinium - based il measured by us earlier
, heptane shows higher solubility in [ c1c5pip][ntf2 ] than in [ c1c3pip][ntf2 ] ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{3}^{\text{il } } $ $ \end{document}x3il  =  0.051 , at t  =  308.15  k ) .
this effect is due to an increase in the van der waals interactions between the hydrocarbon chain of the cation and heptane.table  2compositions of experimental tie lines , solute distribution ratios ,  , and selectivity , s , for ternary systems { [ c1c5pip][ntf2 ] or [ p2,4,4,4][dep ] ( 1 )  +  thiophene or benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k , p  =  101.33  kpahydrocarbon - rich phaseil - rich phase 
  

 s 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{1}^{\text{i } } $ $ \end{document}x1i 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{2}^{\text{i } } $ $ \end{document}x2i 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{1}^{\text{ii } } $ $ \end{document}x1ii 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{2}^{\text{ii } } $ $ \end{document}x2ii 
 [ c1c5pip][ntf2 ]  +  thiophene  +  heptane  0.0000.0000.9110.000  0.0000.0570.7700.1422.4926.7  0.0000.1110.6630.2502.2523.0  0.0000.1910.5640.3521.8417.7  0.0000.2300.5230.3921.7015.4  0.0000.3530.4240.4901.3910.4  0.0000.4470.3760.5421.218.2  0.0000.5300.3370.5821.106.4  0.0000.6960.2810.6500.934.1 
0.0000.8000.2410.7020.883.1  0.0000.8910.2130.7500.842.5  0.0001.0000.1860.8140.81[c1c5pip][ntf2 ]  +  benzothiophene  +  heptane  0.0000.0000.9110.000  0.0000.0250.8030.1094.3648.3  0.0000.0520.7050.2104.0445.0  0.0000.0920.5890.3263.5437.9  0.0000.1610.4440.4722.9329.3  0.0000.2430.3830.5312.1919.2  0.0000.3460.2990.6111.7712.8  0.0000.4600.2360.6741.478.8  0.0000.5530.2080.6961.265.9 
0.0000.7660.1440.7580.992.4  0.0000.8530.0960.8190.961.7  0.0000.9150.0600.8730.951.2[p2,4,4,4][dep ]  +  thiophene  +  heptane  0.0020.0000.5440.000  0.0020.0160.5150.0432.696.0  0.0040.0430.4430.1042.425.1  0.0050.0790.3760.1702.154.3  0.0080.1120.3230.2141.913.6  0.0090.1480.2740.2571.743.1  0.0080.1720.2440.2801.632.8  0.0100.1990.2010.3011.512.4  0.0130.2120.1870.3091.462.2  0.0130.2390.1600.3261.362.0[p2,4,4,4][dep ]  +  benzothiophene  +  heptane  0.0020.0000.5480.000  0.0060.0150.4950.0624.139.1  0.0070.0440.4180.1633.708.4  0.0090.0680.3690.2243.297.5  0.0060.0920.3270.2773.016.9  0.0110.1200.2870.3232.696.0  0.0060.1670.2180.3802.284.7  0.0140.2410.1710.4221.753.2  0.0140.2710.1470.4351.612.7 
0.0150.3110.1150.4531.462.3standard uncertainties are : u(x )  <  0.003 , u(t )  =  0.05  k compositions of experimental tie lines , solute distribution ratios ,  , and selectivity , s , for ternary systems { [ c1c5pip][ntf2 ] or [ p2,4,4,4][dep ] ( 1 )  +  thiophene or benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k , p  =  101.33  kpa standard uncertainties are : u(x )  <  0.003 , u(t )  =  0.05  k the solubility of thiophene at t  =  308.15  k is equal to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{2}^{\text{il } } $ $ \end{document}x2il  =  0.814 for [ c1c5pip][ntf2 ] ( x2il  =  0.797 for [ c1c3pip][ntf2 ] at t  =  298.15  k , the influence of temperature is minimal ; the largest solubility of thiophene in the piperidinium - based il was observed for c1c6pip][ntf2 ] ) .
complete miscibility with thiophene was observed for [ p2,4,4,4][dep ] . in the binary system with benzothiophene ,
the solubility of benzothiophene in [ c1c5pip][ntf2 ] at t  =  308.15  k is equal to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{2}^{\text{il } } $ $ \end{document}x2il  =  0.873 ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{2}^{\text{il } } $ $ \end{document}x2il  =  0.945 for [ c1c3pip][ntf2 ] at t  =  308.15  k .
immiscibility is observed in the { thiophene , or benzothiophene ( 2 )  +  heptane ( 3 ) } binary mixture , as was reported previously .
the determined experimental tie - lines for the ternary lle systems are plotted in figs .  1 , 2 , 3 , 4 for thiophene and benzothiophene , respectively
. figures  1 , 2 , 3 , 4 show that the two - phase region is much larger for [ c1c5pip][ntf2 ] than that for [ p2,4,4,4][dep].fig . 
1plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ c1c5pip][ntf2 ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  kfig . 
2plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ c1c5pip][ntf2 ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  kfig . 
3plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ p2,4,4,4][dep ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  kfig . 
4plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ p2,4,4,4][dep ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ c1c5pip][ntf2 ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ c1c5pip][ntf2 ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ p2,4,4,4][dep ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k plot of the experimental ( filled circle , gray solid lines ) results versus values calculated with the nrtl equation ( square , black dotted lines ) for the composition tie lines of the ternary system { [ p2,4,4,4][dep ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k the results obtained in this work show that the more suitable il for the separation of thiophene , or benzothiophene from heptane , is [ c1c5pip][ntf2 ] because of its much larger selectivity ( s ) and the comparable solute distribution ratio (  ) .
these parameters are defined as follows:1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \beta = \frac{{x _ { 2}^{\text{ii } } } } { { x _ { 2}^{\text{i } } } } $ $ \end{document}=x2iix2i2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ s = \frac{{x _ { 2}^{\text{ii } } \cdot x _ { 3}^{\text{i } } } } { { x _ { 2}^{\text{i } } \cdot x _ { 3}^{\text{ii } } } } $ $ \end{document}s = x2iix3ix2ix3iiwhere x is the mole fraction ; superscripts i and ii refer to the heptane - rich phase and the il - rich phase , respectively .
the values of  and s are listed in table  2 for thiophene and benzothiophene .
figures  5 and 6 present measured values of  and s for ils for thiophene and benzothiophene.fig . 
5plot of the selectivity ( s ) as a function of the mole fraction of solute in the hydrocarbon - rich phase for the ternary systems : closed circle { [ c1c5pip][ntf2 ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } , filled square { [ c1c5pip][ntf2 ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } , open circle { [ p2,4,4,4][dep ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } , and open square { [ p2,4,4,4][dep ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } , at t  =  308.15  kfig . 
6plot of the solute distribution ratio (  ) as a function of the mole fraction of solute in the hydrocarbon - rich phase for the ternary systems : closed circle { [ c1c5pip][ntf2 ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } filled square { [ c1c5pip][ntf2 ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } open circle { [ p2,4,4,4][dep ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } , and open square { [ p2,4,4,4][dep ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } , at t  =  308.15  k plot of the selectivity ( s ) as a function of the mole fraction of solute in the hydrocarbon - rich phase for the ternary systems : closed circle { [ c1c5pip][ntf2 ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } , filled square { [ c1c5pip][ntf2 ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } , open circle { [ p2,4,4,4][dep ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } , and open square { [ p2,4,4,4][dep ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } , at t  =  308.15  k plot of the solute distribution ratio (  ) as a function of the mole fraction of solute in the hydrocarbon - rich phase for the ternary systems : closed circle { [ c1c5pip][ntf2 ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } filled square { [ c1c5pip][ntf2 ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } open circle { [ p2,4,4,4][dep ] ( 1 )  +  thiophene ( 2 )  +  heptane ( 3 ) } , and open square { [ p2,4,4,4][dep ] ( 1 )  +  benzothiophene ( 2 )  +  heptane ( 3 ) } , at t  =  308.15  k the values presented in table  2 show that the distribution ratio coefficient are in the range of 0.812.41 , 0.954.36 , 1.362.69 and 1.464.13 for [ c1c5pip][ntf2]/thiophene , [ c1c5pip][ntf2]/benzothiophene , [ p2,4,4,4][dep]/thiophene and [ p2,4,4,4][dep]/benzothiophene , respectively .
the selectivities of the separation in the system thiophene or benzothiophene / heptane is quite high for [ c1c5pip][ntf2 ] and very low for [ p2,4,4,4][dep ] .
the values listed in table  2 for the best tie - lines are : 26.7 , 48.3 , 6.0 and 9.1 for [ c1c5pip][ntf2]/thiophene , [ c1c5pip][ntf2]/benzothiophene , [ p2,4,4,4][dep]/thiophene and [ p2,4,4,4][dep]/benzothiophene , respectively . in this work
the effect of the alkane chain length on the cation of the piperidinium - based il was examined for comparison with previously measured data for [ c1c3pip][ntf2]/thiophene , [ c1c4pip][ntf2]/thiophene , and [ c1c6pip][ntf2]/thiophene at t  =  298.15  k , and of [ c1c3pip][ntf2]/benzothiophene at t  =  308.15  k   ( the influence of temperature is not large ) .
the characteristic extraction parameters obtained in this work are compared to the few previously described in the open literature in table  3 .
unfortunately , the selectivity for [ c1c5pip][ntf2 ] obtained in this work is slightly worse than that for [ c1c3pip][ntf2 ] measured by us earlier .
the values of selectivity presented for 1-alkylcyanopyridinium - based ils at t  =  308.15  k measured in our earlier work   are also larger than those for piperidinium - based ils   ( see table  3).table  3comparison of solute distribution ratio (  ) and selectivity ( s ) for sulfur compounds extractionillle system 

t / k 
  
 max 

 s 
 max 
 ref.[c1c5pip][ntf2]il  +  thiophene  +  heptane3082.4926.7this work[c1c5pip][ntf2]il  +  benzothiophene  +  heptane3084.3648.3this work[c1c3pip][ntf2]il  +  thiophene  +  heptane2982.5060.3[c1c3pip][ntf2]il  +  benzothiophene  +  heptane3085.3693.0[coc2mpip][ntf2 ] 
 il  +  thiophene  +  heptane2982.6462.9[coc2mpip][fap ] 
 il  +  thiophene  +  heptane2984.0056.8[bcnpy][ntf2 ] 
 il  +  thiophene  +  heptane3081.9362.2[bcnpy][ntf2 ] 
 il  +  benzothiophene  +  heptane3083.50117.1[p2,4,4,4][dep]il  +  thiophene  +  heptane3082.696.0this work[p2,4,4,4][dep]il  +  benzothiophene  +  heptane3084.139.1this work[emim][dep ] 
 il  +  thiophene  +  hexane2982.8848.8[p1,4,4,4][ch3so4]il  +  thiophene  +  cyclohexane2981.385.38[dmim][mp ] 
 il  +  thiophene  +  heptane2980.421,756 
 1-(2-methoxyethyl)-1-methylpiperidinium bis{(trifluoromethyl)sulfonyl}imide 
 1-(2-methoxyethyl)-1-methylpiperidinium trifluorotris(perfluoroethyl)phosphate 
 1-butyl-4-cyanopyridinium bis{(trifluoromethyl)sulfonyl}imide 
 1-ethyl-3-methylimidazolium diethylphosphate 
 1,3-dimethylimidazolium methylphosphonate comparison of solute distribution ratio (  ) and selectivity ( s ) for sulfur compounds extraction 
 1-(2-methoxyethyl)-1-methylpiperidinium bis{(trifluoromethyl)sulfonyl}imide 
 1-(2-methoxyethyl)-1-methylpiperidinium trifluorotris(perfluoroethyl)phosphate 
 1-butyl-4-cyanopyridinium bis{(trifluoromethyl)sulfonyl}imide 
 1-ethyl-3-methylimidazolium diethylphosphate 
 1,3-dimethylimidazolium methylphosphonate the selectivities for [ c1c3pip][ntf2 ] are comparable to those for 4-(2-methoxyethyl)-4-methylpiperidinium trifluorotris(perfluoroethyl)phosphate [ coc2mpip][fap ] ils at t  =  298.15  k , or to 4-(2-methoxyethyl)-4-methylpiperidinium bis{(trifluoromethyl)sulfonyl}imide [ coc2mpip][ntf2 ] .
the extraction results for [ p2,4,4,4][dep ] are very low and similar to [ p1,4,4,4][ch3so4 ] .
it can be definitely concluded that phosphonium - based cations are not suitable for these separation processes .
however , for the diethylphosphate anion [ dep ] and imidazolium - based cation [ emim ] , the results are comparable to those obtained in this work with [ c1c5pip][ntf2 ] but with a lower  value   ( see table  3 ) .
it can be also seen from figs .  5 and 6 that  and s decrease as the solute mole fraction ( thiophene , or benzothiophene ) in the heptane phase increases , for all systems , when going through the tie - line end compositions .
the ternary lle data measured in this study were correlated ( the tie - line correlation ) using the well known non - random liquid equation , nrtl .
the equations and algorithms used for the calculation of the compositions in both phases follow the method described by walas .
the objective function f(p ) was used to minimize the difference between the experimental and calculated compositions:3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ f(p ) = \sum\limits_{i = 1}^{n } { \left [ { x_{2}^{{{\text{i,}}exp } } - x _ { 2}^{\text{i , calc } } \left ( { pt } \right ) } \right]^{2 } + \left [ { x_{3}^{{{\text{i,}}exp } } - x_{3}^{\text{i , calc } } \left ( { pt } \right ) } \right]}^{2 } + \left [ { x_{2}^{{{\text{ii,}}exp } } - x _ { 2}^{\text{ii , calc } } \left ( { pt } \right ) } \right]^{2 } + \left [ { x_{3}^{{{\text{ii,}}exp } } - x _ { 3}^{\text{ii , calc } } \left ( { pt } \right ) } \right]^{2 } $ $ \end{document}f(p)=i=1nx2i , exp - x2i , calcpt2+x3i , exp - x3i , calcpt2+x2ii , exp - x2ii , calcpt2+x3ii , exp - x3ii , calcpt2where p is the set of parameters vector , n is the number of experimental points , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{2}^{{{\text{i}},exp } } $ $ \end{document}x2i , exp , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{3}^{{{\text{i}},exp } } $ $ \end{document}x3i , exp and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x _ { 2i}^{{{\text{i}},{\text{calc } } } } \left ( { pt } \right ) $ $ \end{document}x2ii , calcpt , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{3}^{\text{i , calc } } \left ( { pt } \right ) $ $ \end{document}x3i , calcpt are the experimental and calculated mole fractions of one phase , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{2}^{{{\text{ii}},exp } } $ $ \end{document}x2ii , exp , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x_{3}^{{{\text{ii,}}exp } } $ $ \end{document}x3ii , exp\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x _ { 2}^{\text{ii , calc } } \left ( { pt } \right ) $ $ \end{document}x2ii , calcpt , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ x _ { 3}^{\text{ii , calc } } \left ( { pt } \right ) $ $ \end{document}x3ii , calcpt are the experimental and calculated mole fractions of the second phase .
the binary parameters of each constituent were regressed by minimizing the sum of the squares of the differences between the experimental and calculated mole fractions of each component of both liquid phases for each ternary system .
the value of the non - randomness parameter , ij , was optimized in order to obtain the best model fit .
the correlated parameters are given in table  4 along with the root mean square deviations ( rmsd ) .
the rmsd values , which are a measure of the precision of the correlation , were calculated according the equation:4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \text{rmsd } } = \left ( { \sum\limits_{i } { \sum\limits_{l } { \sum\limits_{m } { \left [ { x_{ilm}^{exp } - x_{ilm}^{\text{calc } } } \right]^{2 } /6k } } } } \right)^{1/2 } $ $ \end{document}rmsd=ilmxilmexp - xilmcalc2/6k1/2where x is the mole fraction and the subscripts i , l , and m designate the component , phase , and tie - line , respectively .
the experimental and calculated lle data agreed relatively well.table  4binary interaction parameters , parameter  
 ij and root mean square deviation (  
 x ) for the nrtl equation for the ternary systems { [ c1c5pip][ntf2 ] or [ p2,4,4,4][dep ] ( 1 )  +  thiophene or benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k , p  =  101.33  kpa 
 ij 
 g 
 12/(jmol)g 
 21/(jmol ) 
  
 ij 
 rmsd  
 x 
 [ c1c5pip][ntf2 ]  +  thiophene  +  heptane  128261.1123225.940.20.007  132536.4315361.68  23460.97475.45[c1c5pip][ntf2 ]  +  benzothiophene  +  heptane  125698.7312608.560.20.007  135142.9534672.57  232498.842232.66[p2,4,4,4][dep ]  +  thiophene  +  heptane  125566.067459.650.20.003  132785.5519670.14  233351.226008.81[p2,4,4,4][dep ]  +  benzothiophene  +  heptane  129265.5911427.530.20.004  132823.0616531.68  231048.012808.72 binary interaction parameters , parameter  
 ij and root mean square deviation (  
 x ) for the nrtl equation for the ternary systems { [ c1c5pip][ntf2 ] or [ p2,4,4,4][dep ] ( 1 )  +  thiophene or benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k , p  =  101.33  kpa
liquid phase equilibrium data were measured in this study for the extraction of thiophene or benzothiophene from heptane using two ils . four ternary systems { il  +  thiophene or benzothiophene  +  heptane } were analytically determined using gc for the composition analysis at temperature t  =  308.15  k at ambient pressure .
it has been demonstrated that the 1-pentyl-1-methylpiperidinium bis{(trifluoromethyl)sulfonyl}imide il is much more effective than the phosphonium - based il for extraction of thiophene or benzothiophene from alkanes .
sulfur compounds can be extracted easily by piperidinium - based ils , leading to low sulfur content in fuels .
our earlier experimental results revealed that the solubility of sulfur compounds in the il increases as the alkyl chain length increases .
the capacity of extraction , described in terms of the selectivity and the solute distribution ratio coefficients , was calculated for all ternary systems and compared to the published data used in similar extraction problems .
based on the values obtained , [ c1c5pip][ntf2 ] was found to be useful for the extraction of sulfur compounds from alkanes ; however , it is not as good as [ c1c3pip][ntf2 ] measured previously . the selectivity and the solute distribution ratio decrease as the mole fraction of thiophene or benzothiophene in the heptane - rich phase increases . the best selectivity ( s )
is observed for very low mole fractions of s - compounds in the hydrocarbon - rich phase \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( x_{2}^{\text{hc } } \ , = \,0.0 5 ) $ $ \end{document}(x2hc=0.05 ) ( see fig . 
5 ) , which may be compared with the results of the hds method for the removal of the s - compounds . the experimental data in this work
the non - randomness parameter was also determined through the reduction of the experimental data .
the model exhibited an excellent fit to the data with the average rmsd values between 0.003 and 0.007 . 

 in this work , the desulfurization ability of alkyl - piperidinium - based and phosphonium - based ionic liquids ( ils ) for ( thiophene or benzothiophene  +  heptane ) mixtures are studied . with this aim , ternary liquid 
  liquid phase equilibrium data ( lle ) have been obtained for mixtures of { il ( 1 )  +  thiophene , or benzothiophene ( 2 )  +  heptane ( 3 ) } at t  =  308.15  k and p  =  101.33  kpa . for this study 1-pentyl-1-methylpiperidinium bis{(trifluoromethyl)sulfonyl}imide , [ c1c5pip][ntf2 ] , and tributylethylphosphonium diethylphosphate , [ p2,4,4,4][dep ] , 
 were used . 
 the suitability of these ils as solvents for extractive desulfurization has been evaluated in terms of the solute distribution ratio and selectivity . 
 immiscibility was observed in the binary liquid systems of ( thiophene , or benzothiophene  +  heptane ) with both ils . 
 one of the studied ils , [ c1c5pip][ntf2 ] , shows high distribution ratios and high selectivities for extraction of sulfur compounds . 
 the data obtained have been correlated with the non - random two liquid nrtl model . 
 the experimental tie - lines and the phase compositions in mole fractions in the ternary systems were calculated with an average root mean square deviation of 0.0052.electronic supplementary materialthe online version of this article ( doi:10.1007/s10953 - 014 - 0276-y ) contains supplementary material , which is available to authorized users . 
Electronic supplementary material
Introduction
Experimental
Results and Discussion
Data Correlation
Conclusions
Electronic Supplementary Material