• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Solvatochromic Parameters and Preferential Solvation Behavior for Binary Mixtures of 1,3-Dialkylimidazolium Ionic Liquids with Water

    2016-09-23 06:06:24ShanDing,Li-gangWei,Kun-lanLi
    CHINESE JOURNAL OF CHEMICAL PHYSICS 2016年4期

    ?

    Solvatochromic Parameters and Preferential Solvation Behavior for Binary Mixtures of 1,3-Dialkylimidazolium Ionic Liquids with Water

    I.INTRODUCTION

    Ⅰonic liquids(ⅠLs),as potential environment-friendly solvents,have gained a great deal of interest with a wide range of applications because of their insignificant vapor pressure,high thermal and chemical stabilities,good solvating capabilities,and excellent recyclability [1,2].To expand the application field ofⅠLs,recent researchers have been focusing on the mixed systems of ⅠL and other solvent,such as alcohols[3]and dimethyl sulfoxide[4].These mixed systems have been credited with special physicochemical properties in comparison with pureⅠL.However,the study of physicochemical properties is more complex in the mixed solvent than that in pureⅠL[5].Ⅰntermolecular interactions occur betweenⅠL and other solvent in the mixture,as well as intramolecular interactions between the cation and anion ofⅠL[6,7].

    Ⅰnvestigating solvent-solvent interactions using solvatochromic indicators is simple and effective because the transition energy of the indicator depends on the solvating sphere composition and properties[8]. Solvatochromic indicator studies offer direct information on solvent properties,such as polarity,dipolarity/polarizability(π?),and hydrogen-bond donating/accepting capabilities.Reichardt proposed the solvatochromic zwitterionic pyridinium N-phenolate betaine dyes,which depended on the observed changes in the ultraviolet-visible(UV-Vis)spectral band shift,to indicate solvent polarity[8-10].Reichardt’s parameter ET(30)and its dimensionless normalized solvent polarity parameter ENTquantify the polarity and acidity,but are not dependent on polarizability.On this basic finding,Kamlet and Taft[11-14]exploited a multiparameter polarity system,solvatochromic parameters,consisting of hydrogen-bond donor(HBD)acidity(α) [12],hydrogen-bond acceptor(HBA)basicity(β)[13],and π?[14],to study solute-solvent and solvent-solvent interactions.Solvatochromicism has been proven to be successful in correlating intermolecular interactions in the binary mixtures[15].Beniwal et al.[16]explored the polarity parameters in both hydrophilic and hydrophobicⅠLs.The results revealed that synergistic effects were generated as a result of high β values of alcohols compared with those of theⅠLs.A strong correlation was observed in the β values and hydrophobicities of pureⅠLs,suggesting that β values can play a role in describing synergism.Similarly,Ali et al.investigated the role of the solvent microsphere in changing the behavior of the solvatochromic absorbance probe within binary solvent mixtures ofⅠLs with common molecular organic solvents[17].The polarity parameters,ENT,and π?of the mixtures,indicated the“hyperpolarity”behavior on addition of molecular solvents toⅠLs.

    The binary mixtures ofⅠL and water(ⅠL-water)have been widely used in separation[18],organic reaction[19],material synthesis[20],pretreatment and delignification of lignocellulosic biomass[21,22].Solvent properties ofⅠL-water mixtures have also been studied using solvatochromic indicator[23-26].Adam et al.explored the solvatochromic response of dipolarity micro-sensors based on Reichardt’s ET(30)and Kamlet-Taft solvent parameters and the application of the solvent exchange model[25].The results confirmed the formation of different intersolvent complexes in the binary mixtures of water and 1-butyl-3-methylimidazolium tetrafluoroborate([C4C1ⅠM]BF4)/1-butyl-3-methylimidazolium bromide([C4C1ⅠM]Br).These complexes provided hydrogen-bonds or electron pairs to the polar network,respectively.Nayak et al.investigated the solvatochromic characterization of the mixtures of a series ofⅠLs and water aimed to describe the changes in the solvent environment as a function of composition [26].They found that the critical solution behavior of poly(N-isopropyl acrylamide)inⅠL-water mixtures was correlated to the values of α and β.These specific properties ofⅠL-water mixtures make them more advantageous for a number of technical processes and then reduce the effect of theⅠLs on environment.To profit fromⅠL-water mixtures,their solvent properties should be fully understood.Sorting out,organizing,and analyzing the solvatochromic parameters of various ⅠL-water binary mixtures systematically are crucial.

    Ⅰn this study,the empirical solvatochromic parameters(ENT,α,β,and π?)of the binary mixtures of water and 1,3-dialkylimidazolium basedⅠLs with different cations and anions were measured over the entire range ofⅠL mole fractions(xIL)at room temperature(25?C),aiming to ascertain the discrepancies in mixed solvent properties.The dyes used in this study were Reichardt’s dye(RD),4-nitroaniline,and 4-nitroanisole,as shown in Fig.1.Preferential solvation model that measures the deviation in ET(30)from linearity with the addition of water inⅠLs was applied to fit the parameters as well.

    II.EXPERIMENTS

    A.Materials

    TheⅠLs used in this work were as follows:1-ethyl-3-methylimidazolium acetate([C2C1ⅠM]OAc),1-butyl-3-methylimidazolium acetate([C4C1ⅠM]OAc),1-hexyl-3-methylimidazolium acetate([C6C1ⅠM]OAc),1-octyl-3-methylimidazolium acetate([C8C1ⅠM]OAc),1-butyl-3-ethylimidazolium acetate([C4C2ⅠM]OAc),1-butyl-3-propylimidazolium acetate([C4C3ⅠM]OAc),1-butyl-3-methylimidazolium chloride([C4C1ⅠM]Cl),1-butyl-3-methylimidazoliumdibutylphosphate([C4C1ⅠM](C4H9O)2PO2),[C4C1ⅠM]Br,and[C4C1ⅠM]BF4. Molecular structures of the usedⅠLs are indicated in Fig.2.TheseⅠLs(purity>98.5%)were obtained from Lanzhou GreenchemⅠLs,LanzhouⅠnstitute of Chemical Physics,Chinese Academy of Sciences.TheseⅠLs were characterized by1H NMR spectroscopy and the spectra were consistent with reported data.The water content of theⅠLs was measured using a Karl Fischer moisture titrator(MKS-500),which didn’t exceed 47 ppm for the usedⅠLs.The residual content of chloride or bromide ions in theⅠLs was determined to be less than 50 ppm by using the Volhard titration method. The details about the impurities in theⅠLs are shown in Table S1 in supplementary materials.

    The RD dye was purchased from Sigma-Aldrich(>90%)and used without additional purification.The dyes of 4-nitroaniline(>99.5%)and 4-nitroanisole(>98%)were purchased from Aladdin Chemistry Co.,Shanghai,China.All reagents were kept out of light and stored at a desiccator.Deionized water with conductivity of 0.056μS/cm was used in this work.

    FⅠG.1Dyes used:(a)Reichardt’s dye(RD),(b)4-nitroaniline,and(c)4-nitroanisole.

    FⅠG.2 Molecular structures ofⅠLs used in this work.

    B.Methods

    UV-Vis spectroscopic measurements were performed with a Perkin-Elmer Lambda 35 spectrophotometer (USA).Polarity measurements were implemented in accordance with the recognized process[23-26].An appropriate concentration of acetone solution of indicators was added toⅠL and mixed adequately.The mixed reagents were dried under vacuum at 80?C for 5 h to remove the acetone.VariousⅠL-water binary mixtures were prepared by weight,with a systematic change in the xILfrom 0 to 1.0.All the indicators were dissolved in a concentration range of 0.5 mmol/L to 1.0 mmol/L. Temperature was maintained at 25±1?C.After 12 h standing to ensure that the mixtures were stable,all samples were tested for UV-Vis spectra with a meanvalue of three times in duplicate.The absorption spectra of the indicator were collected with scopes of λ=250 nm to 400 nm for the 4-nitroanisole,300 nm to 450 nm for the 4-nitroaniline,and 420 nm to 600 nm for the RD.

    C.Solvatochromic polarity analyses

    Solvatochromic polarity parameters were determined as follows.ET(30)polarity was defined as the electronic transition energy of the RD indicator in a particular solvent,which can be obtained from the spectroscopic shift according to Eq.(1).The value was always normalized using water and tetramethylsilane(TMS)as references,as defined by Eq.(2).

    Laurence et al.proposed to use 4-nitroanisole(OMe)as the indicator dye for determining π?[27],in accordance with Eq.(3).

    Hydrogen bond acceptor(HBA)basicity(β)was measured using the spectroscopic shift of the dye 4-nitroaniline(NH2)and can be calculated using Eq.(4).

    Hydrogen-bond donor(HBD)acidity(α)was derived from Reichardt’s scale and can be calculated using Eq.(5)in relation to ET(30)and π?.

    Ⅰdeal mixtures show a linear relationship among the solvatochromic parameter,the solvent composition,and the indicator that are equally solvated by the constituents of mixtures.The experiments performed for predicting the solvatochromic parameters of binary mixtures ofⅠLs and water indicated the measured values didn’t follow the expected ideal behavior.The measured values of solvatochromic parameters for binary mixtures derive from the predicted for linear mixing. Thus,the excess solvatochromic parameters(?SP)was defined as follows.

    Where SPmixture,SP0IL,SP0water,xILand xwaterare the the solvatochromic parameter of the mixture,pureⅠL,pure water,mole fraction ofⅠL and water in the binary mixture,respectively.The?SP can be employed to evaluate the deviations from the ideality.

    III.RESULTS AND DISCUSSION

    Given that the RD is almost insoluble in pure water[24,28],the solvatochromic polarities for pure water(ENT=1.0,β=0.14,α=1.17,π?=1.09)were adopted in this study according to well-accepted values[2,15]. The values of solvatochromic parameters(π?,β,α,ENT) for the binary mixtures ofⅠLs and water at 25?C over the entire range of xIL,are given in Fig.3-Fig.6 and Table S2 in the supplementary materials.The polarity values of the some pureⅠLs andⅠL-water mixtures basically conformed to the values reported in the literatures[23,25,29-33],as shown in Table S3 in the supplementary materials.

    A.The ETNparameter

    Reichardt’spolarityparameterETNreflectsthe chargetransferfromthephenolatetothepyridinium part of zwitterionic molecule and is affected bypolarizability,dipolemoment,andhydrogenbonding.AsshowninFig.3(a)and(c),for the pureⅠLs(xIL=1.0),the trends for ETNvalues were as follows:[C4C1ⅠM]BF4(0.67)>[C4C1ⅠM]Br (0.65)>[C2C1ⅠM]OAc(0.61)>[C4C2ⅠM]OAc(0.60)>

    [C4C1ⅠM]Cl(0.60)>[C4C1ⅠM]OAc(0.59)=[C8C1ⅠM]OAc (0.56)=[C6C1ⅠM]OAc(0.56)>[C4C3ⅠM]OAc(0.55)>

    [C4C1ⅠM](C4H9O)2PO2(0.54).

    The trends of ETNvalues for binary mixtures ofⅠLs and water were similar to that for the pureⅠLs.The presence and position of alkyl substitutes do not affect the variations to a very great extent in the ETNparameters of binary systems.Ⅰt should be noted that the ETNparameters didn’t show a synergistic behavior.Ⅰt was apparently opposite to the trend observed in the binary mixtures ofⅠLs and alcohol[16,17],suggesting that the ⅠL-water mixtures provided a different microenviroment for the RD indicator in comparison withⅠL-alcohol mixtures.

    Figure 3(a)and(c)show that the ENTvalues largely decreased at low xILregion(0-0.2)by adding small quantities of theⅠL.The RD dye has a highly dipolar ground state in comparison with its excited state and therefore exhibits a significant permanent dipole moment,a large polarizable π electron system,a substantial negative charge on the phenoxide O,and a positive charge on the pyridinium N.The strong negative charge on the phenoxide O of RD dye in the electronic ground state can acquire additional solvent stabilization by hydrogen-bonding interactions with theⅠL,which decreased ENTvalues.This suggested that the RD dye is much more preferentially solvated by theⅠL than water.Ⅰt was also observed in Fig.3(a)and(c)that the ENTvalues varied smoothly in the region of xILfrom 0.2 to 1.0,indicating that there wasn’t the nonspecific solvation effects due to the formation of the dye-ⅠL complex.As shown in Fig.3(b)and(d),the ENTparam-eters showed negative deviation.The reason for this nonlinear behavior(nonideal behavior)is“preferential solvation”of the indicator by one component of the binary mixtures.Preferential solvation arises whenever the bulk mole fraction solvent composition is different from the solvation microsphere[25].Ⅰn principle,this phenomenon includes the contribution of specific (e.g.hydrogen-bonding)and non-specific indicatorsolvent interactions.Ⅰt was worthwhile to note that the maximum?ENTwas obtained at xIL=0.1-0.3 for all the binary mixtures ofⅠLs and water,suggesting that there were the strongest interactions betweenⅠL and the RD dye in this xILregion.

    FⅠG.3 The ENTparameters(a,c)and?ENTvalues(b,d)for binary mixtures ofⅠLs and water at 25?C.The maximum ?ENTof[C4C1ⅠM](C4H9O)2PO2/H2O was≤±0.006.

    B.The α and β parameter

    As shown in Fig.4(a)and(c),for the pure ⅠLs(xIL=1.0),the order of α values was as follows.[C4C1ⅠM]BF4(0.66)>[C4C2ⅠM]OAc(0.48)>[C4C1ⅠM]Br(0.46)>[C2C1ⅠM]OAc(0.44)=[C8C1ⅠM]OAc (0.44)>[C4C1ⅠM]OAc(0.43)>[C4C1ⅠM](C4H9O)2PO2(0.42)>[C6C1ⅠM]OAc(0.40)>[C4C3ⅠM]OAc(0.39)>[C4C1ⅠM]Cl(0.38).

    Ⅰt was found that the parameter α was mainly affected by cation,but with an anion influence as well.The parameter β is a measure of the HBA basicity of a solvent.As shown in Fig.5(a)and (c),for the pureⅠLs,the order of β values was as follows.[C4C3ⅠM]OAc(1.14)>[C4C1ⅠM]OAc(1.13) =[C4C2ⅠM]OAc(1.13)>[C4C1ⅠM](C4H9O)2PO2(1.08) >[C6C1ⅠM]OAc(1.07)>[C2C1ⅠM]OAc(1.04)=[C8C1ⅠM] OAc(1.04)>[C4C1ⅠM]Cl(0.79)>[C4C1ⅠM]Br(0.70)>[C4C1ⅠM]BF4(0.51).

    Given uniform cation with 1-butyl-3-methyl imidazolium([C4C1ⅠM]+),β values obviously changed with different anions.Conversely,given uniform anion with acetate(OAc-),β values were relatively stable. The β values for[C4C3ⅠM]OAc,[C4C1ⅠM]OAc and [C4C2ⅠM]OAc were slightly higher than that for the otherⅠLs with acetate anions.The phenomenon mainly stemmed from the fact that the butyl substituent made the imidazolium rings highly symmetric.Generally,the parameter β of pureⅠL was largely affected by the nature of the anion,but there was a small cation effect. The polarity of solvent depends on both cations and anions ofⅠLs.

    Water is a good HBD because of its fairly high α value(α=1.17)and relatively low β scale(β=0.14).On the contrary,an opposite phenomenon of α and β scales for pureⅠLs was observed.ⅠLs,as HBA,made it easy to build hydrogen-bond with water.Variation in Fig.4 and Fig.5 of solvatochromic parameters for allⅠL-water binary systems provides valuable information that the solvent properties of binary mixtures can be affected by solvent composition(water content)distinctly.

    FⅠG.4 The parameters α(a,c)and?α values(b,d)for binary mixtures ofⅠLs and water at 25?C.The maximum?α of [C4C1ⅠM](C4H9O)2PO2/H2O was≤±0.021.

    FⅠG.5 The parameter β(a,c)and?β values(b,d)for the binary mixtures ofⅠLs and water at 25?C.The maximum?β of[C4C1ⅠM](C4H9O)2PO2/H2O was≤±0.022.

    The α and β scales changed dramatically asⅠL was initially added in water.The β values increased by almost 0.5 unit,while the α values decreased signifi-cantly.[C6C1ⅠM]OAc/H2O,[C8C1ⅠM]OAc/H2O,and [C4C1ⅠM](C4H9O)2PO2/H2O systems decreased initially,whereas others were the exact contrary.This phenomenon was reasonable because the addition ofⅠLs induced competitions in molecular interactions among ⅠLs and water.A commonly increasing trend of β values and a moderately decreasing trend of α values were observed with increasing xIL.ⅠnⅠL-rich regions(in the xILrange of 0.4 to 1.0),similar trends for the parameters α and β were followed;all values began to flatten and were basically steady,which can be due to the approximate composition of binary mixture by the time.This behavior was exactly in line with what was observed by Moita et al.[34],who considered that the competitions in molecular interactions of transition energies for pairs depended on composition range.

    The values α and β(ⅠLs-water)were between that of pureⅠLs and water.This result demonstrated no synergistic effect in these imidazolium-basedⅠL-water binary systems.The trends of parameter α for binary mixtures of water andⅠLs with acetate anion were tanglesome on the whole,as depicted in Fig.4 (a)and(c).[C8C1ⅠM]OAc/H2O system is on the high side,and[C4C3ⅠM]OAc/H2O system is the least. Sarkar et al.[24]deemed that the HBD acidity of [C4C1ⅠM]BF4/H2O system stemmed from the strong acidity for H on the two positions of imidazolium ring,which is in accordance with our results.The entire imidazolium ring protons are acidic,which provides a possibility to hydrogen-bonding with anions and dye indicator-solute as HBD.Ⅰn microstructural perspective,a highly symmetric cation leads to a high hydrophobicity degree of the cation cluster.As a result,this phenomenon is less likely to provide a proton to form hydrogen-bond.The α values for binary mixtures are determined by the nature of cations and anions together,but cation is the main factor.

    Considering that basicity is dominated by the characteristic of anions,the tendency of β values for the binary mixtures of water andⅠLs with different anions was specifically distinct as shown in Fig.5(a)and(c). For 1-butyl-3-methylimidazoliumⅠLs-water mixtures,[C4C1ⅠM](C4H9O)2PO2/H2O and[C4C1ⅠM]OAc/H2O systems show high HBA basicity because of the presence of oxyanions.The β values generally have a close relationship with the capability ofⅠL-water mixtures to dissolve lignin[35,36].However,the order of β values can’t exactly explain this dissolving behavior.

    As shown in Fig.4 and Fig.5,the β values displayed a positive deviation(?β>0),whereas the α values expressed a negative deviation(?α<0).Such a difference revealed preferential solvation of the indicator by the most polar constituent.Harifi-Mood et al.studied the solvatochromic behavior among water,methanol,ethanol and[C4C1ⅠM]BF4,and found a similar trend of solvent properties forⅠL-water mixtures[23].As depicted in Fig.4(b)and(d),the α values of[C4C3ⅠM]OAc/H2O showed large negative deviation in the case of xIL=0.2,which was due to high symmetric cation of[C4C3ⅠM]OAc.As shown in Fig.5 (b)and(d),the β values of[C2C1ⅠM]OAc/H2O and [C4C1ⅠM](C4H9O)2PO2/H2O presented large positive deviations in the case of around xIL=0.1,which indicated that the basic character of theseⅠLs resembled more different from that of water than otherⅠLs.The reason for this finding can be that theⅠL with small size of cation or large size of anion can donate their protons easier than other usedⅠLs.

    Ⅰt should be noted that the binary mixtures of [C4C1ⅠM]BF4/H2O owned the highest HBD acidity and the smallest deviations from the ideality,as shown in Fig.4.The solvatochromic parameters of [C4C1ⅠM]BF4/H2O measured in this work were similar to the data reported by Harifi-Mood et al.[23] However,F(xiàn)reire et al.[37]proposed that the hydrolysis of BF4-was observed under even at room temperature and the hydrolysis extent was increased with increasing time and temperature.Thus,in this work,the[C4C1ⅠM]BF4/H2O mixtures were measured for a short time for reducing the influence of the hydrolysis of BF4-.Even so,the influence of the hydrolysis of BF4-on the solvatochromic parameters should be further investigated.Here,the[C4C1ⅠM]BF4/H2O were only used in theoretically comparing the solvatochromic parameters of the binary mixtures,as did in other researches[16,25].

    C.Net basicity(β-α)

    Hauru et al.discussed the relationship between solvent parameters and cellulose regeneration fromⅠL solutions and found that the capability ofⅠL-water system to dissolve cellulose can be interpreted by net basicity (β-α)rather than β alone[38].Shi et al.reported a similar conclusion that theⅠLs with high net basicity could dissolve cellulose more efficiently[35].

    Net basicity(β-α)for all 1,3-dialkylimidazoliumbasedⅠL-water binary mixtures over the entire range of mole fractions is presented in TableⅠ.For theseⅠLs,while trace water was added,the value increased only from 0.1 unit to 0.4 unit.By contrast,with the addition of traceⅠLs,the value increased by almost 1.8 units at maximum.This phenomenon revealed that severe intermolecular interactions occurred in solution,which fitted with the former discussion.

    For pureⅠL,the cellulose dissolution located in a region is roughly defined by 0.35<(β-α)<0.9 and 0.80<β<1.20[21].Shi et al.reported that cellulose can still be dissolved in the mixture of[C2C1ⅠM]OAc/H2O(ⅠL content of 20wt%)[35].TableⅠindicates that most of(β-α)values of 1,3-dialkylimidazolium-based ⅠL-water systems are within the range of the solvent window,which gave a relatively unsophisticated account of the solvation behavior.Wang et al.reported that lignin could be dissolved in a series of 1,3-alkylimidazolium-basedⅠL-water mixtures except for [C4C1ⅠM]BF4/H2O[36].As shown in TableⅠ,β-α values of[C4C1ⅠM]BF4/H2O were negative,whereas those of most ofⅠL-water mixtures were positive.Hence,the lignin dissolution capacity of[C4C1ⅠM]BF4/H2O was weaker than that of otherⅠL-water mixtures.

    The preceding results indicated that the net basicity can be used to predict the dissolution behavior of cellulose and lignin inⅠL-water mixtures.Net basicity was considered as an important solvatochromic parameter for the binary mixtures ofⅠL and water.

    TABLEⅠNet basicity(β-α)for allⅠL-water binary mixtures.

    D.Dipolarity/polarisability(π?)

    Theparameterπ?isameasureofdipolarity/polarisability of the solvent,which is derived from the variation of energy corresponding to the dye OMe maximum absorption.Ⅰt is induced by the local electric field generated by the solvent.The π?values of pureⅠLs are depicted in Fig.6(a)and(c).The order of these π?was as follows:[C4C1ⅠM]Cl(1.19)=[C4C1ⅠM]Br(1.19) >[C2C1ⅠM]OAc(1.13)>[C4C1ⅠM]OAc(1.08)>[C6C1ⅠM] OAc(1.03)=[C4C2ⅠM]OAc(1.03)>[C4C3ⅠM]OAc(1.01) >[C8C1ⅠM]OAc(0.98)=[C4C1ⅠM]BF4(0.98)>[C4C1ⅠM] (C4H9O)2PO2(0.95).

    For 1,3-dialkylimidazolium acetate,the π?values decreased with increasing length of alkyl chain of R1 (C3-C8)or R3(C1-C3)(Fig.6(a)and(c)),suggesting that the π?can be affected by the position and size of substituent.The π?value for pure water(1.07)was higher than for some pureⅠLs,such as[C4C1ⅠM]BF4and[C4C1ⅠM](C4H9O)2PO2,whereas was lower than that for[C4C1ⅠM]Cl and[C4C1ⅠM]Br(Fig.6(b)).Theπ?values were related to the size of the anion ofⅠLs. With increasing in the size of the anion,the π?values decreased for the usedⅠLs.

    The π?values for theseⅠL-water mixtures are also shown in Fig.6(a)and(c).Ⅰt was not surprising that π? parameter was greatly influenced by increasing solutesolvent interactions when the xILwas increased.Ⅰt was noted that the dramatic change of π?values occurred in the xILrange from 0 to 0.2,suggesting that the micro environment of binary mixtures substantially changed due to the molecular interactions betweenⅠL and water.The π?values kept relative steady with increasing xILfrom 0.2 to 1.0.Ⅰt was worthwhile to be noted that binary mixtures of[C4C1ⅠM]OAc/H2O,[C4C2ⅠM]OAc/H2O and[C4C3ⅠM]OAc/H2O showed strong synergistic effects,in which π?values of the mixture is higher than that of either of the pure solvents.This synergistic effect was attributed to the formation of hydrogen-bond complexes,whose properties were different from either of the solvents.Very small levels of synergism were observed in the binary mixtures of[C2C1ⅠM]OAc/H2O,[C4C1ⅠM]Cl/H2O and [C4C1ⅠM]Br/H2O,which was due to strong intramolecular interactions between anions and cations of these ⅠLs.For[C6C1ⅠM]OAc/H2O,[C8C1ⅠM]OAc/H2O and [C4C1ⅠM](C4H9O)2PO2/H2O)systems,small levels of synergism were also found,which was probably due to the very large size of the cations and anions,suggesting that steric hindrance restricts interactions between the ⅠLs and the dye.No synergistic effects were found for the binary mixtures of[C4C1ⅠM]BF4/H2O.This finding was consistent with other reports[23,24].The synergistic behavior in case of some π?scales validated the existence of complex S12(displaying a mixed solvent formed by the interaction of pure solvents 1 and 2)with a higher polarity than pure components have.This phenomenon didn’t necessarily indicate that no S12 existed in other systems,but S12 was not significant.

    The ENT,α and β values didn’t present the synergy with varying xIL(Fig.3-Fig.5),whereas the values showed the synergistic effects for most of binary mixtures(Fig.6).Ⅰt suggested that the synergy weren’t entirely characteristic of the studied mixtures but determined by the systems employed to study the mixtures.

    FⅠG.6 The π?parameter(a,c)and?π?values(b,d)for the binary mixtures ofⅠLs and water at 25?C.The maximum error associated with π?of[C4C1ⅠM](C4H9O)2PO2/H2O was≤±0.026,and the details can be found in Table S2 in supplementary materials.

    The selection of the systems,for example the polarity indicator,was decisive in determining the nature of the molecular interactions and the extent of synergy occurred.Ⅰn this work,the π?parameters measured using the dye 4-nitroanisole could reflect the synergy.Ⅰt was difficult to speculate about the nature of interactions between the solvent and the dye in this stage.But such difference in synergy in these sameⅠL-water mixtures emphasized the requirements for more extensive studies in the future.

    As illustrated in Fig.6(b)and(d),the π?valuesmeasuredforallthebinarymixturesdidn’t follow the expected ideal behavior.The π?values produced a negative deviation for the binary mixturesof[C6C1ⅠM]OAc/H2O,[C8C1ⅠM]OAc/H2O,[C4C1ⅠM]BF4/H2O,[C4C1ⅠM](C4H9O)2PO2/H2O). For[C6C1ⅠM]OAc/H2Oand[C8C1ⅠM]OAc/H2O,asymmetrical cation and anion of theⅠL came into being clusterization,which was not conducive to deformability.Containing a tetrahedral structure anion group,the deformability of(C4H9O)2PO2-and BF4-is poor,which directly led to a negative deviation of π?values.Ⅰt may be attributed to the specific interactions between theⅠL and water.As shown in Fig.6(b)and (d),the π?values produced a positive deviation for the binary mixtures due to their high electronegativity anions and sphere microstructure.

    As depicted in Fig.3-Fig.6,the maximum deviation to ideality for the solvatochromic parameters can be obtained in the xILrange from 0.1 to 0.3.Ⅰts origin must be some interesting molecular interactions.Thermodynamic studies of imidazolium-basedⅠL-water systems have revealed that cation and counterion seem completely dissociated only below xILof 0.013 to 0.015.At higher xILthan this threshold,ions begin to interact with each other and organize themselves by mutual attraction.At a higher concentration,such as xIL>0.5 to 0.6,ⅠL molecules form clusters as the pure state,and water molecules interact withⅠL clusters without forming a high hydrogen-bond network among themselves [39].As indicated in Fig.3-Fig.6,the solvatochromic parameters dramatically changed in a range of low xIL(0-0.3)and kept the stable stage in a range of high xIL (0.3-1.0).The solvatochromic parameters seemed to be closely related to the microstructure of the binary mixtures.

    E.Preferential solvation model

    Considering that indicator dye can be preferentially solvated by either binary mixtures,water orⅠL,the selection must be made with the help of a preferential solvation model.Skwierczynski and Connors[40]proposed a phenomenological model of solvent effects,which were based on a simple two-step solvent exchange according to Eq.(7)and Eq.(8),to quantitatively describe solventsolvent(the general medium effect),solvent-solute(the solvation effect),and solute-solute(the intersolute ef-fect)interactions in aqueous-organic binary mixtures. Bosch et al.developed the model for a wide application successfully[41,42].Nevertheless,this model can’t explain synergetic behavior[43].

    As stated previously,allⅠL-water binary mixtures show a non-synergetic effect on the solvatochromic parameters of ENT,α,and β.The generalized model below is applied in this work toⅠL-water binary systems investigated.

    TABLEⅠⅠFitting parameters obtained from the preferential solvation models in the selected binary mixtures at 25?C.

    where S1 and S2 indicate two pure solvents,and S12 displays a mixed solvent formed by the interaction of solvents 1 and 2.Ⅰrepresents the corresponding RD indicator.Ⅰ(S1)2,Ⅰ(S2)2,Ⅰ(S12)2stand for the RD solvated by S1,S2,and S12 components,respectively.The equilibrium constants for two solvent-exchange processes,f2/1and f12/1,measure the tendency of the probe to be solvated by solvent S2 and“mixed solvent”S12 with reference to solvent S1.

    where xs1,xs2,and xs12are the mole fractions of solvents S1,S2,and S12 in the sphere of solvation of the indicator,respectively.x01and x02are the mole fractions of the two solvents in the bulk mixed solvent.

    According to this model,Reichardt’s polarity ETof binary mixtures as a function of mole fraction x2can be calculated.

    Ⅰn Eq.(11),a and c are the regression coefficients that minimize the square residuals of ETvalues.For the convenience of calculation,Skwierczynski and Connors[40] assumed ETpolarity of solvent S12 to be the average of ET1and ET2,as presented in Eq.(12).

    To understand the unique natures ofⅠL-water at the molecular-level aspect,two parameters,f2/1and f12/1,were measured by fitting the wavenumbers of maximum absorptions of the RD at different compositions of S1 (water)and S2(ⅠLs)to the proposed Eq.(11).The results are listed in TableⅠⅠ.For allⅠL-water binary mixtures,f2/1was generally greater than 1,which reflected that the indicator was more solvated byⅠL rather than water,and the solvation sphere was substantiated with ⅠL.Similarly,f12/1commonly exceeds 1,except for system[C8C1ⅠM]OAc/H2O.f2/1>f12/1exhibited thatⅠL built a closer solute-solvent interaction with the indicator in comparison with the mixture.Macrocosmically,the difference in the values of f2/1and f12/1demonstrated that the nature of anions and cations played an important role in the preferential solvation study. [C4C1ⅠM]BF4/H2O system,however,was not suitable for this preferential solvation model.

    The parameters f2/1and f12/1differed from each other for acetate-basedⅠL-water mixtures,which can affirm that the alkyl substitution on the imidazolium ring greatly affected the solute-solvent interactions of solvation sphere.With increasing alkyl chain of cation on one position,the value of f2/1first increased sharply and then decreased,whereas f12/1decreased continuously.Given that RD dye molecule is large,the phenomenon can be attributed to the steric-hindrance effect and symmetry of the large substitution on one position of the ring.Therefore,[C2C1ⅠM]OAc/H2O system held lower f2/1and higher f12/1than others do,which reflected that“mixed solvent”,combination of water and[C2C1ⅠM]OAc,was easier to form and more active to interact with the indicator than otherⅠL-watermixtures.On the contrary,the parameter f12/1for [C8C1ⅠM]OAc/H2O mixture was inferior to 1,which might be ascribed to the fact that large asymmetric cation putted a brake on the formation of the“mixed solvent”that directly led to a weak solute-solvent interaction.As for substituting proton on three positions,the values of f2/1and f12/1for[C4C3ⅠM]OAc/H2O system were higher than those for the two other systems. Hence,the symmetry of different alkyl group of cation tied the solute-solvent interactions up.From another aspect,the parameters f2/1and f12/1for the systems with butyl substitute on one position of cation were higher than the others.

    The preferential parameters f2/1and f12/1for 1-butyl-3-methylimidazolium-basedⅠL-watersystems were determined.The overall values followed f2/1>f12/1>1,which reflected that the dye is strongly preferentially solvated by pureⅠLs in comparison with the mixture with reference to water.Danten et al. [44]reported that an almost symmetric structure of anion···H-O-H···anion is formed by hydrogen-bonds inⅠL-water mixtures.When solute was added,similar to the dye indicator,a competition existed between anion and dye indicator,thereby indicating that the solvent-solvent interaction always strongly influenced solute-solvent interaction.For systems that shared the same cations,the observed trend for preferentialsolvationwas[C4C1ⅠM](C4H9O)2PO2/H2O>[C4C1ⅠM]Cl/H2O>[C4C1ⅠM]OAc/H2O>[C4C1ⅠM]Br/ H2O.This trend was closely related to the electronegativity and magnitude of anions.

    IV.CONCLUSION

    Theempiricalsolvatochromicparameters(ENT,α,β,andπ?)ofthebinarymixturesof1,3-dialkylimidazoliumⅠLs and water were investigated. This result indicated that solvent properties ofⅠL-water mixtures were influenced by both the xILandⅠL type. The behavior of the solvatochromic absorbance indicator obviously deviated from ideality in these mixtures. The maximum deviation to ideality was obtained in the xILrange from 0.1 to 0.3.Net basicity(β-α)for these ⅠL-water mixtures can predict the capacity of dissolving lignin.For most ofⅠL-water mixtures(except for [C4C1ⅠM]BF4/H2O),the π?parameters were found to be higher than that of pureⅠLs and water instead of other parameters(ENT,α and β).The synergy in π?parameters indicated interesting and strong solute-solvent and solvent-solvent interactions in these studied systems.With the preferential solvation model,the dye probe is preferentially solvated in the following order: pureⅠL>mixtures>water.This work will help to understand the physicochemical properties of theseⅠL-water mixtures and facilitate their practical applications.

    Supplementary materials:The impurities contents in the usedⅠLs,the measured solvatochromic parameters of the binary mixtures ofⅠLs and water,and comparison between reference data and the measured values in this work are given.

    V.ACKNOWLEDGMENTS

    This work is supported by the National Natural Science Foundation of China(No.21276034).

    [1]N.V.Plechkova and K.R.Seddon,Chem.Soc.Rev. 37,123(2008).

    [2]T.L.Greaves and C.J.Drummond,Chem.Rev.115,11379(2015).

    [3]O.Russina,A.Sferrazza,R.Caminiti,and A.Triolo,J.Phys.Chem.Lett.5,1738(2014).

    [4]F.Huo,Z.Liu,and W.Wang,J.Phys.Chem.B 117,11780(2013).

    [5]R.Sharma and R.K.Mahajan,RSC Adv.4,748 (2014).

    [6]S.Khatun and E.W.Castner Jr.,J.Phys.Chem.B 119,9225(2015).

    [7]Y.Zhao,J.Wang,H.Wang,Z.Li,X.Liu,and S.Zhang,J.Phys.Chem.B 119,6686(2015).

    [8]C.Reichardt,Chem.Rev.94,2319(1994).

    [9]C.Reichardt and T.Welton,Solvents and Solvent Effects in Organic Chemistry,4th Edn.,Weinheim: Wiley-VCH,132(2010).

    [10]C.Reichardt,Pure Appl.Chem.76,1903(2004).

    [11]R.W.Taft and M.J.Kamlet,J.Am.Chem.Soc.98,2886(1976).

    [12]M.J.Kamlet and R.W.Taft,J.Am.Chem.Soc.98,377(1976).

    [13]M.J.Kamlet,J.L.Addoud,and R.W.Taft,J.Am. Chem.Soc.99,6027(1977).

    [14]T.Yokoyama,R.W.Taft,and M.J.Kamlet,J.Am. Chem.Soc.98,3233(1976).

    [15]P.G.Jessop,D.A.Jessop,D.B.Fu,and L.Phan,Green Chem.14,1245(2012).

    [16]V.Beniwal and A.Kumar,ChemPhysChem 16,1026 (2015).

    [17]A.Ali,M.Ali,N.A.Malik,S.Uzair,and A.B.Khan,J.Chem.Eng.Data 59,1755(2014).

    [18]A.F.M.Cl′audio,M.C.Neves,K.Shimizu,J.N.C. Lopes,M.G.Freire,and J.A.Coutinho,Green Chem. 17,3948(2015).

    [19]D.W.Kim,D.J.Hong,J.W.Seo,H.S.Kim,H.K. Kim,C.E.Song,and D.Y.Chi,J.Org.Chem.69,3186(2004).

    [20]X.C.Kang,J.L.Zhang,W.T.Shang,T.B.Wu,P. Zhang,B.X.Han,Z.H.Wu,G.Mo,and X.Q.Xing,J.Am.Chem.Soc.136,3768(2014).

    [21]T.V.Doherty,M.Mora-Pale,S.E.Foley,R.J.Linhardt,and J.S.Sordick,Green Chem.12,1967(2010).

    [22]L.G.Wei,K.L.Li,Y.C.Ma,and X.Hou,Ⅰnd.Crop. Prod.37,227(2012).

    [23]A.R.Harifi-Mood,A.Habibi-Yangjeh,and M.R.Gholami,J.Phys.Chem.B 110,7073(2006).

    [24]A.Sarkar and S.Pandey,J.Chem.Eng.Data 51,2051 (2006).

    [25]C.G.Adam,M.V.Bravo,P.M.E.Mancini,and G. G.Fortunato,J.Phys.Org.Chem.27,841(2014).

    [26]P.K.Nayak,A.P.Hathorne,and H.Bermudez,Phys. Chem.Chem.Phys.15,1806(2013).

    [27]C.Laurence,P.Nicolet,M.T.Dalati,J.L.M.Abboud,and R.Notario,J.Phys.Chem.98,5807(1994).

    [28]S.N.Baker,G.A.Baker,and F.V.Bright,Green Chem.4,165(2002).

    [29]M.A.Ab Rani,A.Brant,L.Crowhurst,A.Dolan,M. Lui,N.H.Hassan,J.P.Hallett,P.A.Hunt,H.Niedermeyer,J.M.Perez-Arlandis,M.Schrems,T.Welton,and R.Wildinga,Phys.Chem.Chem.Phys.13,16831 (2011).

    [30]A.V.Ladesov,D.S.Kosyakov,K.G.Bogolitsyn,and N.S.Gorbova,Russ.J.Phys.Chem.A 89,1814(2015).

    [31]P.M.Mancini,G.G.Fortunato,C.G.Adam,and L. R.Vottero,J.Phys.Org.Chem.21,87(2008).

    [32]V.Strehmel,R.Lungwitz,H.Rexhausen,and S. Spange,New J.Chem.34,2125(2010).

    [33]Y.Cao,Y.Chen,X.Wang,and T.Mu,RSC Adv.4,5169(2014).

    [34]M.L.C.J.Moita,A.F.S.Santos,J.F.C.C.Silva,andⅠ.M.S.Lampreia,J.Chem.Eng.Data 57,2702 (2012).

    [35]J.Shi,K.Balamurugan,R.Parthasarathi,N.Sathitsuksanoh,S.Zhang,V.Stavila,V.Subramanian,B.A. Simmons,and S.Sinqh,Green Chem.16,3830(2014).

    [36]Y.T.Wang,L.G.Wei,K.L.Li,Y.C.Ma,N.N.Ma,S.Ding,L.L.Wang,D.Y.Zhao,B.Yan,W.Y.Wan,Q.Zhang,X.Wang,J.M.Wang,and H.Li,Bioresour. Technol.170,499(2014).

    [37]M.G.Freire,C.M.S.S.Neves,Ⅰ.M.Marrucho,J.A. P.Coutinho,and A.M.J.Fernandes,J.Phys.Chem. A 114,3744(2010).

    [38]L.K.J.Hauru,M.Hummel,A.W.T.King,Ⅰ. Kilpelainen,and H.Sixta,Biomacromolecules 13,2896 (2012).

    [39]L.Zhang,Z.Xu,Y.Wang,and H.Li,J.Phys.Chem. B 112,6411(2008).

    [40]R.D.Skwierczynski and K.A.Connors,J.Chem.Soc. Perkin Trans.2,467(1994).

    [41]U.Buhvestov,F(xiàn).Rived,C.Rafols,E.Bosch,and M. Roses,J.Phys.Org.Chem.11,185(1998).

    [42]M.Roses,U.Buhvestov,C.Rafols,F(xiàn).Rived,and E. Bosch,J.Chem.Soc.Perkin Trans.2,1341(1997).

    [43]N.D.Khupse and A.Kumar,J.Phys.Chem.B 115,711(2011).

    [44]Y.Danten,M.Ⅰ.Cabao,and M.Besnard,J.Phys. Chem.A 113,2873(2009).

    Shan Ding,Li-gang Wei?,Kun-lan Li,Ying-chong Ma
    School of Light Industry and Chemical Engineering,Dalian Polytechnic University,Dalian 116034,China
    (Dated:Received on January 6,2016;Accepted on March 15,2016)
    Binary mixtures of 1,3-dialkylimidazolium based ionic liquids(ⅠLs)and water were selected as solvent systems to investigate the solute-solvent and solvent-solvent interactions on the preferential solvation of solvatochromic indicators at 25?C.Empirical solvatochromic parameters,dipolarity/polarizability(π?),hydrogen-bond donor acidity(α),hydrogen-bond acceptor basicity(β),and Reichardt’s polarity parameters(ENT)were measured from the ultraviolet-visible spectral shifts of 4-nitroaniline,4-nitroanisole,and Reichardt’s dye.The solvent properties of theⅠL-water mixtures were found to be influenced byⅠL type andⅠL mole fraction(xIL).All these studied systems showed the non-ideal behavior.The maximum deviation to ideality for the solvatochromic parameters can be obtained in the xILrange from 0.1 to 0.3.For most of the binary mixtures,the π?values showed the synergistic effects instead of the ENT,α and β values.The observed synergy extent was dependent on the studied systems,such as the dye indicator andⅠL type.A preferential solvation model was utilized to gather information on the molecular interactions in the mixtures.The dye indicator was preferentially solvated on the following trend:ⅠL>ⅠL-water complex>water.
    Key words:Solvatochromic parameters,Preferential solvation,Ⅰonic liquid,Water

    ?

    Author to whom correspondence should be addressed.E-mail: wei-ligang@hotmail.com,Tel.:+86-411-86323726

    建设人人有责人人尽责人人享有的| 性高湖久久久久久久久免费观看| 免费观看人在逋| 男人操女人黄网站| 亚洲欧美日韩高清在线视频 | 伊人亚洲综合成人网| 日韩伦理黄色片| av网站在线播放免费| 亚洲欧美一区二区三区久久| 亚洲av在线观看美女高潮| 亚洲,欧美,日韩| 在线看a的网站| 在线观看免费午夜福利视频| 国精品久久久久久国模美| 亚洲成人免费电影在线观看 | 亚洲情色 制服丝袜| 不卡av一区二区三区| 国产在线一区二区三区精| 99国产精品免费福利视频| 电影成人av| 国产成人欧美在线观看 | 久久99一区二区三区| 老司机午夜十八禁免费视频| 美女主播在线视频| av在线播放精品| 日韩一区二区三区影片| 欧美老熟妇乱子伦牲交| 男女免费视频国产| 久久久国产一区二区| 一本综合久久免费| 一个人免费看片子| 亚洲国产精品一区二区三区在线| 又粗又硬又长又爽又黄的视频| 亚洲精品国产av成人精品| 韩国精品一区二区三区| 久久久国产精品麻豆| 成人黄色视频免费在线看| 搡老乐熟女国产| 九色亚洲精品在线播放| 热99久久久久精品小说推荐| 大香蕉久久成人网| 亚洲精品久久午夜乱码| 在线亚洲精品国产二区图片欧美| 久久久国产精品麻豆| 欧美日韩福利视频一区二区| 久久国产亚洲av麻豆专区| 亚洲精品第二区| 波多野结衣av一区二区av| 国产高清不卡午夜福利| 久久久久视频综合| 男人爽女人下面视频在线观看| 十八禁高潮呻吟视频| 亚洲中文av在线| 交换朋友夫妻互换小说| 亚洲天堂av无毛| 最新在线观看一区二区三区 | 中文字幕亚洲精品专区| 午夜福利免费观看在线| 女人精品久久久久毛片| 精品卡一卡二卡四卡免费| 18禁观看日本| 精品国产一区二区三区久久久樱花| 国产成人一区二区三区免费视频网站 | 久久久久视频综合| 黄色视频在线播放观看不卡| 亚洲视频免费观看视频| 在线观看免费日韩欧美大片| 亚洲色图综合在线观看| 亚洲精品一区蜜桃| 51午夜福利影视在线观看| 首页视频小说图片口味搜索 | 亚洲国产欧美一区二区综合| 久久久精品94久久精品| 丝袜在线中文字幕| 免费高清在线观看日韩| 成人三级做爰电影| 日韩 亚洲 欧美在线| 国产精品香港三级国产av潘金莲 | 久久精品久久久久久久性| 亚洲男人天堂网一区| 欧美精品一区二区免费开放| 啦啦啦视频在线资源免费观看| 午夜视频精品福利| 久久女婷五月综合色啪小说| 国产一区亚洲一区在线观看| 欧美日韩精品网址| 久久精品国产亚洲av涩爱| 久久女婷五月综合色啪小说| 国产成人一区二区三区免费视频网站 | av在线播放精品| 每晚都被弄得嗷嗷叫到高潮| 国产熟女欧美一区二区| 男女边吃奶边做爰视频| 国产无遮挡羞羞视频在线观看| 一本一本久久a久久精品综合妖精| 欧美日韩成人在线一区二区| 精品少妇久久久久久888优播| 婷婷成人精品国产| 成人亚洲精品一区在线观看| 中文精品一卡2卡3卡4更新| 制服人妻中文乱码| 丰满饥渴人妻一区二区三| 欧美人与善性xxx| 亚洲一卡2卡3卡4卡5卡精品中文| 一区二区三区乱码不卡18| 午夜福利视频在线观看免费| 亚洲一区二区三区欧美精品| 91精品国产国语对白视频| 女性被躁到高潮视频| 久久精品aⅴ一区二区三区四区| 亚洲一区二区三区欧美精品| 热re99久久国产66热| av线在线观看网站| 伊人亚洲综合成人网| av天堂久久9| 另类精品久久| 天堂中文最新版在线下载| 各种免费的搞黄视频| 欧美精品啪啪一区二区三区 | 欧美日本中文国产一区发布| 成人国产av品久久久| 熟女少妇亚洲综合色aaa.| 下体分泌物呈黄色| 日本黄色日本黄色录像| 久久久久视频综合| 精品国产超薄肉色丝袜足j| 日本黄色日本黄色录像| 精品欧美一区二区三区在线| 超碰97精品在线观看| 黄色一级大片看看| 大片电影免费在线观看免费| 亚洲熟女毛片儿| 99re6热这里在线精品视频| 操出白浆在线播放| 老司机亚洲免费影院| 久久久精品免费免费高清| 汤姆久久久久久久影院中文字幕| 精品高清国产在线一区| 人妻 亚洲 视频| 中文欧美无线码| 日日夜夜操网爽| 亚洲av在线观看美女高潮| 国产成人精品久久久久久| 多毛熟女@视频| 国产成人欧美| 男女免费视频国产| 啦啦啦在线免费观看视频4| 精品人妻1区二区| 亚洲国产精品一区二区三区在线| xxxhd国产人妻xxx| 天天影视国产精品| 观看av在线不卡| 好男人视频免费观看在线| 国产精品久久久av美女十八| 成人黄色视频免费在线看| 天天影视国产精品| 久久久亚洲精品成人影院| 天堂8中文在线网| 午夜视频精品福利| 天堂中文最新版在线下载| 男女国产视频网站| 乱人伦中国视频| 亚洲av日韩在线播放| 国产日韩一区二区三区精品不卡| 男女午夜视频在线观看| 免费高清在线观看视频在线观看| 丁香六月欧美| 丝袜在线中文字幕| 亚洲,欧美,日韩| 亚洲欧美一区二区三区黑人| 国产91精品成人一区二区三区 | 亚洲精品日韩在线中文字幕| 大片电影免费在线观看免费| 国产片特级美女逼逼视频| 日韩视频在线欧美| 欧美精品人与动牲交sv欧美| 欧美精品av麻豆av| 美女国产高潮福利片在线看| 国产一区有黄有色的免费视频| 欧美激情高清一区二区三区| 日本vs欧美在线观看视频| 好男人视频免费观看在线| 午夜免费成人在线视频| 久久人妻福利社区极品人妻图片 | 国产爽快片一区二区三区| 国产精品二区激情视频| 18在线观看网站| 国产高清videossex| 亚洲五月色婷婷综合| av在线老鸭窝| 亚洲七黄色美女视频| 99久久精品国产亚洲精品| 精品国产国语对白av| 国产爽快片一区二区三区| 欧美激情极品国产一区二区三区| 国产精品偷伦视频观看了| 搡老乐熟女国产| 亚洲专区国产一区二区| 久热这里只有精品99| 色94色欧美一区二区| 亚洲欧洲国产日韩| av天堂久久9| 午夜激情av网站| 波多野结衣av一区二区av| 女性生殖器流出的白浆| 欧美日韩成人在线一区二区| 精品一区在线观看国产| 色婷婷久久久亚洲欧美| 亚洲av在线观看美女高潮| 国产在线视频一区二区| 嫁个100分男人电影在线观看 | 国产日韩欧美在线精品| 免费观看人在逋| 99久久精品国产亚洲精品| 国产av国产精品国产| 水蜜桃什么品种好| 国产真人三级小视频在线观看| 国产成人精品在线电影| 久9热在线精品视频| 久久ye,这里只有精品| 日韩免费高清中文字幕av| 亚洲精品国产av蜜桃| 中文字幕人妻丝袜一区二区| 久久精品成人免费网站| 久久久久视频综合| 色播在线永久视频| 中文字幕亚洲精品专区| 黄色视频不卡| 国产成人影院久久av| 国产精品三级大全| 国产片特级美女逼逼视频| 好男人电影高清在线观看| 欧美日本中文国产一区发布| 日本av免费视频播放| 久久久久久久久免费视频了| 欧美精品亚洲一区二区| 免费少妇av软件| 99国产综合亚洲精品| 日韩大片免费观看网站| 亚洲国产欧美网| 女人被躁到高潮嗷嗷叫费观| av视频免费观看在线观看| 亚洲国产精品一区二区三区在线| av片东京热男人的天堂| 久久综合国产亚洲精品| 国产精品久久久久久精品电影小说| 涩涩av久久男人的天堂| 色网站视频免费| 一区二区av电影网| 成年女人毛片免费观看观看9 | 国产片内射在线| 亚洲一区中文字幕在线| 无遮挡黄片免费观看| 亚洲九九香蕉| 日韩伦理黄色片| 99国产精品99久久久久| 日本欧美国产在线视频| 精品一品国产午夜福利视频| 99精国产麻豆久久婷婷| 亚洲欧美精品自产自拍| 丰满人妻熟妇乱又伦精品不卡| 免费看av在线观看网站| 九草在线视频观看| 久久久精品国产亚洲av高清涩受| 欧美中文综合在线视频| 午夜激情av网站| avwww免费| 天堂中文最新版在线下载| 精品一区在线观看国产| 99热全是精品| 另类亚洲欧美激情| 美女大奶头黄色视频| 日韩av免费高清视频| 老司机影院成人| 久久精品熟女亚洲av麻豆精品| 99国产精品一区二区蜜桃av | 欧美黄色淫秽网站| 99国产精品99久久久久| 十八禁高潮呻吟视频| 国产伦人伦偷精品视频| 日韩中文字幕视频在线看片| 老司机亚洲免费影院| 亚洲,欧美精品.| 十八禁人妻一区二区| 啦啦啦中文免费视频观看日本| 老汉色∧v一级毛片| 欧美激情极品国产一区二区三区| 69精品国产乱码久久久| 国产99久久九九免费精品| 99久久人妻综合| 中文字幕人妻丝袜制服| 亚洲综合色网址| 中文乱码字字幕精品一区二区三区| 国产免费又黄又爽又色| 日韩人妻精品一区2区三区| 一级毛片 在线播放| 国产成人a∨麻豆精品| 亚洲av成人不卡在线观看播放网 | 免费高清在线观看日韩| 在线av久久热| 亚洲精品在线美女| 国产一区二区激情短视频 | 日本vs欧美在线观看视频| 国产成人免费无遮挡视频| 嫩草影视91久久| 看十八女毛片水多多多| 女人精品久久久久毛片| 又黄又粗又硬又大视频| 日韩精品免费视频一区二区三区| 亚洲av日韩在线播放| 久久久久视频综合| 亚洲精品日本国产第一区| 在线观看免费视频网站a站| 人成视频在线观看免费观看| 成人亚洲精品一区在线观看| 最新在线观看一区二区三区 | 午夜影院在线不卡| 成人免费观看视频高清| 国产欧美亚洲国产| 国产一区二区在线观看av| 亚洲欧美清纯卡通| 国产精品久久久av美女十八| www.av在线官网国产| 欧美日韩成人在线一区二区| 午夜免费鲁丝| 18禁观看日本| 人体艺术视频欧美日本| 日韩中文字幕视频在线看片| 久久久久国产精品人妻一区二区| 国产精品二区激情视频| 国产麻豆69| 色播在线永久视频| 老汉色∧v一级毛片| 国产精品.久久久| 久久精品亚洲av国产电影网| 高清不卡的av网站| 国产xxxxx性猛交| 国产精品偷伦视频观看了| 欧美日韩视频高清一区二区三区二| 国产日韩一区二区三区精品不卡| 999精品在线视频| 视频区欧美日本亚洲| 一本—道久久a久久精品蜜桃钙片| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲伊人色综图| 男的添女的下面高潮视频| 啦啦啦在线免费观看视频4| 免费少妇av软件| 精品熟女少妇八av免费久了| 免费日韩欧美在线观看| 日本黄色日本黄色录像| 国产精品人妻久久久影院| 欧美成人午夜精品| 91成人精品电影| 最近中文字幕2019免费版| 国产成人免费观看mmmm| 午夜福利影视在线免费观看| 久久精品熟女亚洲av麻豆精品| 免费人妻精品一区二区三区视频| 高潮久久久久久久久久久不卡| 91字幕亚洲| 女性生殖器流出的白浆| 久久久国产一区二区| 久久久精品区二区三区| 夫妻性生交免费视频一级片| 国产男女内射视频| 黄色 视频免费看| 精品久久久久久久毛片微露脸 | 欧美日韩视频精品一区| 欧美乱码精品一区二区三区| 免费看不卡的av| 视频区图区小说| 亚洲欧美中文字幕日韩二区| 亚洲av美国av| 性少妇av在线| 曰老女人黄片| 国产在线免费精品| 日韩大片免费观看网站| 一级黄片播放器| 亚洲国产日韩一区二区| 亚洲国产欧美在线一区| 久久久久久免费高清国产稀缺| 色94色欧美一区二区| 午夜久久久在线观看| 老司机靠b影院| cao死你这个sao货| 一边亲一边摸免费视频| 亚洲欧美中文字幕日韩二区| 午夜福利一区二区在线看| 国产成人av教育| 国产黄色视频一区二区在线观看| 欧美日本中文国产一区发布| 最新在线观看一区二区三区 | 欧美黑人欧美精品刺激| av片东京热男人的天堂| 国产精品免费视频内射| 波野结衣二区三区在线| 制服人妻中文乱码| 午夜免费成人在线视频| 又大又爽又粗| 丝袜美足系列| 国产不卡av网站在线观看| 少妇人妻 视频| 97精品久久久久久久久久精品| 老汉色av国产亚洲站长工具| 国产精品人妻久久久影院| 蜜桃国产av成人99| 色94色欧美一区二区| 国产高清国产精品国产三级| videos熟女内射| 免费高清在线观看视频在线观看| 国产免费一区二区三区四区乱码| avwww免费| av片东京热男人的天堂| 亚洲欧美激情在线| 岛国毛片在线播放| 大片电影免费在线观看免费| 久久人人97超碰香蕉20202| 大话2 男鬼变身卡| 久久天躁狠狠躁夜夜2o2o | 久久影院123| 欧美精品一区二区免费开放| 亚洲精品中文字幕在线视频| 久久久久久久久久久久大奶| 91成人精品电影| 无遮挡黄片免费观看| 久久精品人人爽人人爽视色| 亚洲av日韩在线播放| 欧美国产精品一级二级三级| 丝袜脚勾引网站| 大陆偷拍与自拍| 久久久国产一区二区| 伊人久久大香线蕉亚洲五| 精品福利永久在线观看| 一级毛片 在线播放| 99热全是精品| 天堂中文最新版在线下载| www.精华液| 欧美成人精品欧美一级黄| 国产欧美日韩精品亚洲av| 国产精品三级大全| 韩国高清视频一区二区三区| 午夜福利,免费看| 少妇的丰满在线观看| 亚洲精品国产av蜜桃| 婷婷色麻豆天堂久久| av在线播放精品| 国产精品国产av在线观看| 最新在线观看一区二区三区 | 啦啦啦啦在线视频资源| 性高湖久久久久久久久免费观看| 亚洲欧洲日产国产| 伊人久久大香线蕉亚洲五| 亚洲国产精品一区二区三区在线| 中国国产av一级| 超碰97精品在线观看| 国产高清不卡午夜福利| 成人手机av| 50天的宝宝边吃奶边哭怎么回事| 九草在线视频观看| 亚洲欧美成人综合另类久久久| 日日夜夜操网爽| 天天影视国产精品| 伊人亚洲综合成人网| 一级片'在线观看视频| 老司机在亚洲福利影院| 国产一级毛片在线| 老汉色av国产亚洲站长工具| a级毛片在线看网站| 亚洲综合色网址| 老司机影院成人| 久久国产精品男人的天堂亚洲| 黄频高清免费视频| 狂野欧美激情性xxxx| www日本在线高清视频| 国产免费又黄又爽又色| 中文字幕另类日韩欧美亚洲嫩草| 十分钟在线观看高清视频www| 亚洲av片天天在线观看| 亚洲欧美激情在线| 国产主播在线观看一区二区 | 女警被强在线播放| 免费在线观看日本一区| 久久久久久久大尺度免费视频| 中文字幕人妻丝袜一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久精品精品| 亚洲欧美激情在线| 国产av一区二区精品久久| 亚洲欧美精品自产自拍| av网站在线播放免费| 一二三四在线观看免费中文在| 亚洲成人国产一区在线观看 | 91精品国产国语对白视频| 午夜免费观看性视频| 又大又黄又爽视频免费| 曰老女人黄片| 久久久久精品人妻al黑| 国产高清videossex| 91字幕亚洲| 亚洲精品国产一区二区精华液| 成年av动漫网址| 国产精品 国内视频| 男人爽女人下面视频在线观看| 欧美黑人精品巨大| 国产一区亚洲一区在线观看| 亚洲成人免费av在线播放| 免费不卡黄色视频| www.自偷自拍.com| 欧美日韩福利视频一区二区| 欧美成人精品欧美一级黄| 麻豆国产av国片精品| 麻豆乱淫一区二区| 天天操日日干夜夜撸| 中文字幕精品免费在线观看视频| 国产精品av久久久久免费| 日日爽夜夜爽网站| 亚洲国产日韩一区二区| 国产99久久九九免费精品| 久久鲁丝午夜福利片| 亚洲情色 制服丝袜| kizo精华| 中国国产av一级| 国产又爽黄色视频| 可以免费在线观看a视频的电影网站| 免费在线观看影片大全网站 | 精品第一国产精品| 免费在线观看完整版高清| 黑人欧美特级aaaaaa片| 午夜免费观看性视频| 男女高潮啪啪啪动态图| 免费看不卡的av| 亚洲国产精品一区三区| 欧美黑人欧美精品刺激| 天堂8中文在线网| 亚洲精品美女久久久久99蜜臀 | 考比视频在线观看| av一本久久久久| 91麻豆精品激情在线观看国产 | 青草久久国产| av不卡在线播放| 精品久久久久久久毛片微露脸 | 亚洲精品国产一区二区精华液| 国产主播在线观看一区二区 | 成人午夜精彩视频在线观看| 免费av中文字幕在线| 国产精品久久久人人做人人爽| 国产亚洲午夜精品一区二区久久| 另类亚洲欧美激情| 国产亚洲欧美精品永久| 男女无遮挡免费网站观看| 午夜久久久在线观看| 久久毛片免费看一区二区三区| 国产精品 欧美亚洲| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 91精品三级在线观看| 国产av国产精品国产| 国产亚洲欧美在线一区二区| 美女福利国产在线| 国产一区二区激情短视频 | 精品高清国产在线一区| 亚洲成色77777| 一本综合久久免费| 大陆偷拍与自拍| 美国免费a级毛片| 国产福利在线免费观看视频| 97人妻天天添夜夜摸| 日韩av不卡免费在线播放| 欧美成人午夜精品| 制服诱惑二区| 亚洲 国产 在线| 日本vs欧美在线观看视频| 亚洲成人免费av在线播放| 国产男人的电影天堂91| 亚洲av成人不卡在线观看播放网 | 久久久久久亚洲精品国产蜜桃av| 国产成人系列免费观看| 国产激情久久老熟女| 国产91精品成人一区二区三区 | 免费高清在线观看日韩| 久久久久久久精品精品| 亚洲av日韩精品久久久久久密 | 免费黄频网站在线观看国产| av网站在线播放免费| 国产一卡二卡三卡精品| 日日夜夜操网爽| 国产精品av久久久久免费| 女性被躁到高潮视频| 18在线观看网站| 久久精品久久久久久久性| 女警被强在线播放| 男人舔女人的私密视频| 国产成人91sexporn| 美女中出高潮动态图| 亚洲欧洲日产国产| 男女床上黄色一级片免费看| 亚洲精品av麻豆狂野| 国产免费现黄频在线看| 高清欧美精品videossex| 99国产精品一区二区蜜桃av | 五月天丁香电影| 最新的欧美精品一区二区| 午夜视频精品福利| 久久久久久亚洲精品国产蜜桃av| av欧美777| 久久影院123| 视频区欧美日本亚洲| 热99国产精品久久久久久7| 男的添女的下面高潮视频| 精品国产超薄肉色丝袜足j| 男人操女人黄网站| 新久久久久国产一级毛片| 欧美另类一区| 精品少妇黑人巨大在线播放|