Figures
Abstract
In this work, we have explored the thermophysical properties of tetraalkylammonium hydroxide ionic liquids (ILs) such as tetrapropylammonium hydroxide (TPAH) and tetrabutylammonium hydroxide (TBAH) with isomers of butanol (1-butanol, 2-butanol and 2-methyl-2-propanol) within the temperature range 293.15–313.15 K, with interval of 5 K and over the varied concentration range of ILs. The molecular interactions between ILs and butanol isomers are essential for understanding the function of ILs in related measures and excess functions are sensitive probe for the molecular interactions. Therefore, we calculated the excess molar volume (VE) and the deviation in isentropic compressibility (Δκs) using the experimental values such as densities (ρ) and ultrasonic sound velocities (u) that are measured over the whole compositions range at five different temperatures (293.15, 298.15, 303.15, 308.15 and 313.15 K) and atmospheric pressure. These excess functions were adequately correlated by using the Redlich–Kister polynomial equation. It was observed that for all studied systems, the VE and Δκs values are negative for the whole composition range at 293.15 K. And, the excess function follows the sequence: 2-butanol>1-butanol>2-methyl-2-propanol, which reveals that (primary or secondary or tertiary) position of hydroxyl group influence the magnitude of interactions with ILs. The negative values of excess functions are contributions from the ion-dipole interaction, hydrogen bonding and packing efficiency between the ILs and butanol isomers. Hence, the position of hydroxyl group plays an important role in the interactions with ILs. The hydrogen bonding features between ILs and alcohols were analysed using molecular modelling program by using HyperChem 7.
Citation: Attri P, Baik KY, Venkatesu P, Kim IT, Choi EH (2014) Influence of Hydroxyl Group Position and Temperature on Thermophysical Properties of Tetraalkylammonium Hydroxide Ionic Liquids with Alcohols. PLoS ONE 9(1): e86530. https://doi.org/10.1371/journal.pone.0086530
Editor: Danilo Roccatano, Jacobs University Bremen, Germany
Received: August 30, 2013; Accepted: December 11, 2013; Published: January 29, 2014
Copyright: © 2014 Attri et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: We gratefully acknowledge SRC program of National Research Foundation of Korea (NRF) Grant funded by the Korean Government (MEST) (No. 20100029418) and in part by Kwangwoon University 2013. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Till date numerous research groups have focused their work on the study of fascinating physical properties of ionic liquids (ILs), due to their wide variety of applications in industries and applied chemistry [1]–[16]. For the applications in chemical and industrial processes, the knowledge of the thermophysical properties of IL is essential, as they represent the basis for the chemical and biological process [17]–[28]. Many of these studies have led to new possible applications for ILs [18]–[21]. Knowledge of structure and properties of ILs is essential for the understanding of their molecular interactions in the binary mixtures [4]–[6], [8]–[16]. Nevertheless, in order to transfer the ILs from laboratory to industry, designing future processes and equipment involving these ionic compounds, an accurate knowledge about their physical properties, either for pure ILs or mixed with other solvents, is crucial. Therefore, a deep knowledge of thermophysical properties of ILs and their liquid mixtures are essentially required for scientific community. Apparently, the physicochemical properties of ILs are quite sensitive toward the structure and nature of cations and anions [7]–[11]. The variations in thermophysical properties of ILs, such as density (ρ) and speed of sound (u) are observed to be very sensitive to the change in ion, mainly due to the microscopic level interactions between solvent molecules [12]–[16].
Binary mixtures of ILs with other solvents can also improve the thermodynamic and transport properties of working fluids as well as the efficiency of the chemical equipments such as batteries, photoelectrical cells, and other electrochemical apparatus. The use of the binary mixtures of ILs with polar compounds such as alcohols allows the change and control of the properties of the mixtures to suit a given situation [29].Thermodynamic properties of mixtures containing ILs and alcohols are important for both the design of many technological processes and an understanding of the solute–solvent interactions in the mixtures. These properties are required in the development of models for process design, energy efficiency, and in the evaluation of possible environmental impacts [30]. Regarding the study of physical properties for binary mixtures of alcohol+ILs, a large number of works have been published in recent years [31]–[41], showing the interest of the scientific community for this field. While, still there is no experimental or theoretical results are available for the thermophysical properties between the tetraalkylammonium hydroxide and butanol isomers. Additionally, there is no study to show the interactions between the hydroxide anion of the IL and hydroxyl group of alcohols.
In this research to study these interaction, we explore and compare the measurements of two thermophysical properties such as ρ and u of binary mixtures involving 1-butanol, 2-butanol and 2-methyl-2-propanol with tetrapropylammonium hydroxide [(C3H7)4N][OH] (TPAH) and tetrabutylammonium hydroxide [(C4H9)4N][OH] (TBAH) ILs over a complete mole fraction range at various temperatures from 293.15 to 313.15 K, with interval of 5 K. Further, the excess molar volume (VE), and deviation in isentropic compressibilities (Δκs) were calculated using experimental data. The resulting VE and Δκs values were found to be strongly dependent on the place of hydroxyl group attached in the chain and also on the interactions between the hydroxide anion of ILs and hydroxyl group of the alcohols. These deviations in physical parameters have been explained in terms of intermolecular interactions between alcohols and ILs. Additionally, the temperature also plays an important role in interaction studies. Moreover, the hydrogen bonding features between ILs and alcohols were carried out to get a deep insight into intermolecular interactions for the studied compounds. These studies were performed according to the semi-empirical calculations by using HyperChem 7.
Materials and Methods
Materials
1-butanol, 2-butanol and 2-methyl-2-propanol were obtained from Merck >99% of purity and stored over freshly activated 3 Å molecular sieves and were purified by the standard method described by Riddick et al [42]. A comparison is made for the pure alcohols in Table 1 between the experimental ρ and u values determined in the present study and those reported in the literature [31], [38], [41]–[50]. ILs were synthesized in laboratory and analysed using 1H-NMR, the preparation is given below.
Synthesis of ILs
Synthesis of Tetrapropylammonium Hydroxide (TPAH).
The synthesis of this IL was carried out in a 250 mL round bottomed flask, which was immersed in a water-bath, fitted with a reflux condenser. Solid potassium hydroxide (40 mmol) was added to a solution of tetrapropylammonium bromide [(C3H7)4N][Br] (40 mmol) in dry methylene chloride (20 mL), and the mixture was stirred vigorously at room temperature for 10 h. The precipitated KBr was filtered off, and the filtrate was then evaporated to leave the crude [(C3H7)4N][OH] as a viscous liquid that was washed with ether (2×20 mL) and dried at 343.15 K for 5 h to obtain the pure IL. The sample was analyzed by Karl Fisher titration and revealed very low levels of water (below 70 ppm). The yield of TPAH was 82%. 1H NMR (DMSOd6): δ (ppm) 0.8 (t, 12H), 1.46 (m, 8H), 2.92 (t, 8H), 4.56 (s, OH). HRMS calculated for C12H29NO (M+ - OH) 203.36, found 203.25.
Synthesis of Tetrabutylammonium Hydroxide (TBAH).
A procedure similar to that above for [(C3H7)4N][OH] was followed with the exception of the use of [(C4H9)4N][Br] ([cation]) instead of [(C3H7)4N][Br]. The yield of TBAH was 82%. 1H NMR (DMSOd6): δ (ppm) 0.94 (t, 12H), 1.37 (m, 8H), 1.96 (m, 8H), 3.43 (t, 8H), 4.78 (s, OH). HRMS calculated for C16H37NO (M+ - OH) 259.47 found out to be 259.34.
Experimental Procedure
Density (ρ) and speed of sound (u) measurements.
The density (ρ) and speed of sound (u) measurements were performed with an Anton-Paar DSA 5000 with an accuracy of temperature of ±0.01 K. The uncertainties in the density and speed of sound measurements were ±0.00005 g cm−3 and 0.01 m s−1 respectively. Prior to measurements, the instrument was calibrated with deionized water and dry air as standards at 293.15 K.
The binary mixtures of butanol isomers and IL were prepared by mass using a high-precision analytical balance with an uncertainty of ±1×10−4 g. All of the samples were prepared immediately before the measurements to avoid variations in composition due to evaporation of the solution. Clear and air bubble free solutions were used to perform the ρ and u experiments at different temperatures. The detailed measurement procedures used were described in detail in our previous research papers [12]–[15].
Hydrogen Bonding through Simulation Program.
The structures of ILs and alcohols were optimized based on molecular mechanics and semi-empirical calculations using the HyperChem 7 molecular visualization and simulation program [51]–[54]. Initial molecular geometry of butanol isomers and ILs were optimized with the PM3 semi-empirical calculations and single point calculations were carried out to determine the total energies. Now the optimized molecules, alcohols and IL were chosen and then placed on top of each other symmetrically (parallel) with a starting interplanar distance of 2.3 Å and the angle made by covalent bonds to the donor and acceptor atoms less than 1200 was fulfilled. Further, the geometries were optimized using geometry optimizations based on molecular mechanics (using the MM+force field) and PM3 semi-empirical calculations, the Polak-Ribiere routine with rms gradient of 0.01 as the termination condition was used. PM3 uses a set of parameters derived from a variety of experimental versus calculated molecular properties, as compared to other semiempirical methods, including the AM1 procedure [53]. Typically, nonbonded interactions are less repulsive in the PM3 procedure [54]. Hydrogen bonds were displayed using HyperChem “show hydrogen bonds” and “recompute hydrogen bond” options.
Results and Discussion
In order to have the better understanding of the molecular interactions between tetraalkylammonium hydroxide ILs with polar solvents such as 1-butanol, 2-butanol and 2-methyl-2-propanol, we have measured ρ and u properties over the whole composition range at various temperatures such as 293.15, 298.15, 303.15, 308.15 and 313.15 K under atmospheric pressure. The experimental ρ and u values of ILs with alcohols are presented as a function of IL concentration in Table S1 in File S1. Further, Figures 1–4 show the measured ρ and u values for the binary mixtures of different butanol isomers with both ILs (TPAH and TBAH) at all the studied temperatures. Figures 1 and 2 reveal that the variation of ρ values of TPAH or TBAH with 1-butanol, 2-butanol and 2-methyl-2-propanol, shows similar trends. It has been found that the ρ of the mixtures increased with the increasing concentrations of the ILs in alcohols. The effect of the ILs on the ρ in the alcohols has been examined at various temperatures. It has been observed that the ρ values decreased as temperature increased in the all systems. The results in Figure 1a clearly reveal that the ρ values of the TPAH+1-butanol mixture increase sharply up to x1≈0.8000 and later become almost constant at all the temperatures. While for TPAH+2-butanol, ρ values increase up to very rich IL concentration x1≈0.9900, as shown in Figure 1b. On the other hand, ρ values for TPAH+2-methyl-2-propanol increase sharply up to x1≈0.6400 and, no prominent changes have been observed afterwards. The increase in ρ values for TPAH+alcohols mixtures is possibly due to increase in the ion pair interactions between TPAH and alcohols. This shows that the density values of the TPAH+alcohol mixtures are not affected much due to change in the position of hydroxyl group in different isomers.
The solid line represents the smoothness of these data.
The solid line represents the smoothness of these data.
The solid line represents the smoothness of these data.
The solid line represents the smoothness of these data.
Whereas, in Figure 2a the ρ values for the TBAH+1-butanol mixtures increases sharply up to x1≈0.6200, later the increase is marginally very less at high concentration of IL region. Whereas, TBAH+2-butanol mixtures have similar trends as shown earlier by TPAH+2-butanol, ρ values increase up to very rich IL concentration x1≈0.9900, as depicted by Figure 2b. Further, TBAH+2-methyl-2-propanol mixture shows the increase in ρ values, while ρ doesn't increase sharply at mole fraction 0.5000–0.9900, which may be due to decrease in ion-pair interactions between TBAH and 2-methyl-2-propanol, as shown in Figure 2c. From Table S1 in File S1, we observed that the densities of investigated systems increase with increasing the length of alkyl chain in IL. It was found that ρ values to be higher in the TBAH+butanol isomers as compared to TPAH+butanol isomers at equimolar mixture. Whereas, according to early documented research articles the density decreases with increase in alkyl chain in a cation or anion [55], [56]. These discrepancies vary from IL to IL and solvent to solvent and also depend on the nature as well as structural arrangement of IL and solvent. Moreover, from close look on the Table S1 in File S1, we observed that with increase in temperature, ρ values of TBAH+butanol isomers decreases more as compared to TPAH+butanol isomers. This might be due to the assumption that the ion-pair interaction decreases more for high alkyl chain+butanol isomers as compared to lower alkyl chain+butanol isomers with the increase in temperature.
Ultrasonic sound velocities (u) prove to be an informative source regarding the properties of different solvents and their mixture. The values of u were found to decrease with an increase in temperature while u values increased with increasing in mole fraction of IL. As noted from Figures 3 and 4, there is a sharp increase of u in all ILs, except in the mixture of TPAH with 1-butanol at 293.15 K, in the mole fraction range from 0.8000 to 0.9900 of IL. Over this range, the u values decrease slightly for the mixtures of TPAH with 1-butanol at 293.15 K. Whereas, no change is observed in rest of the IL+butanol isomers at all investigated temperatures. This u value is significantly increased in IL-solvent interactions when the mole fraction of IL was increased. If we compare the TBAH+1-butanol to TPAH+1-butanol, it has been observed that the u values slightly decrease when the alkyl substituents size of cation increases. Whereas, the same trend was observed on comparing the u values of TBAH+2-methyl-2-propanol to TPAH+2-methyl-2-propanol. It has been found that the u values slightly decrease when the size of cation increases. While the u value for TPAH+2-butanol are lower than TBAH+2-butanol, which again reveal that u values slightly increase as the size of cation increases. Hence, our results lead to conclusion that interactions of ILs with alcohols, depends upon the position of the hydroxyl group.
Thermophysical properties of mixed solvents of ILs with butanol isomers can be tunable. The extent of deviation of liquid mixtures from ideal behavior is best expressed by excess functions. Excess molar volumes (VE) as well as ultrasonic studies are known to provide useful insights into solution structural effects and intermolecular interactions between component molecules. The extent of deviation of liquid mixtures from ideal behavior is best expressed by excess functions. Volumetric properties of binary mixtures of ILs with polar compounds are contributed to the clarification of the various intermolecular interactions existing between the different species found in solution. The excess volumes are determined from the density of pure compounds (ρ1 and ρ2) and mixture (ρm) using a standard equation [15]. The ultrasonic studies have been adequately employed in understanding the nature of molecular interaction in solvent mixed systems. In the chemical industry, knowledge of the ultrasonic and its related properties of solutions are essential in the design involving chemical separation, heat transfer, mass transfer, and fluid flow. Isentropic compressibilities (κs) of the binary mixtures were calculated using the relation from ρ and u. The composition dependence of the VE and Δκs properties represents the deviation from ideal behavior of the mixtures and provides an indication of the interactions between IL and alcohols. These properties were mathematically fitted by variable degree functions using the Redlich-Kister expression:(1)
(2)Where Y refers to VE or Δκs. ai are adjustable parameters and can be obtained by least-squares analysis. Values of the fitted parameters are listed in Table S2 in File S1, along with the standard deviations of the fit. The values of VE and Δκs for the binary mixtures at various temperatures as function of ILs concentrations are included in Table S1 in File S1. Figures 5 to 10 display the experimental data for the binary mixtures, and the fitted curves, along with the excess properties of VE and Δκs for the butanol isomers with ILs as function of IL concentrations at different temperatures.
Solid lines correlated by the Redlich-Kister equation.
Solid lines correlated by the Redlich-Kister equation.
Solid lines correlated by the Redlich-Kister equation.
Solid lines correlated by the Redlich-Kister equation.
Solid lines correlated by the Redlich-Kister equation.
Solid lines correlated by the Redlich-Kister equation.
From Figure 5 one can note that the values of VE are negative for all TPAH+butanol isomers systems at all measured temperatures over whole composition range. We have observed that the excess molar volumes present a minimum at x1≈0.5634 for TPAH+1-butanol at all investigated temperatures, whereas we obtained that the VE values present a minimum at x1≈0.4143 for the TPAH+2-butanol system. Further, minimum VE values lie at x1≈0.4244 for TPAH+2-methyl-2-propanol system at all investigated temperatures. The minimum VE values could be due to hydrogen bonds between alcohols and TPAH IL. The decrease in the magnitude of the negative VE values with an increase in the IL composition can be attributed to the decrease of hydrogen bonding. In other words, due to increase in the concentration of the IL results in decrease of packing efficiency. Further, with increase in temperature the magnitude of the negative VE values decreases in all the TPAH+butanol isomers systems. This can again be due to decrease in magnitude of the hydrogen bonding with increase in the temperature.
It is interesting to note that the VE values in TPAH+2-methyl-2-propanol mixture show more negative values of VE at the alcohol-rich composition than the TPAH+1-butanol and TPAH+2-butanol mixtures at 293.15 K (Figure 7a), implying that in the TPAH+2-methyl-2-propanol, there are ion-dipole interactions and packing effects with 2-methyl-2-propanol which are stronger than those in the 2-butanol and 1-butanol solution at x1≈0.4200. A comparison between the negative deviation of VE of TPAH+2-methyl-2-propanol, TPAH+1-butanol and TPAH+2-butanol suggests that there is a difference of the hydroxyl position in the alkyl chain, leading to variation in the interactions between the alcohols and TPAH. Molecular interaction between TPAH and alcohols follows the following order at 293.15 K, 2-methyl-2-propanol>1-butanol>2-butanol.
Further, the negative VE values are observed for TBAH+butanol isomers at all measured temperatures over whole composition range, except 2-butanol at higher temperatures 293.15 to 313.15 K. The VE values for 1-butanol with TBAH are as represented in Figure 6a. And, we found that negative VE values are observed over the entire mole fraction range at all investigated temperatures. These negative VE values reveal that a more efficient packing or attractive interaction occurred between the TBAH and 1-butanol. 1-Butanol forms a hydrogen bond with the alkyl chain cation, while the interactions decrease at higher temperatures. The interactions between the 1-butanol molecules and the alkyl chain of TBAH are due to ion-dipole or hydrogen bonding interactions. This will reduce the interactions between the tetrabutylammonium cation and hydroxide anion in the IL, which contributes to the negative VE values. Furthermore, the observed positive VE values for TBAH+2-butanol at higher temperatures show that there exist no specific interactions between unlike molecules, as displayed in Figure 6b. The magnitude and sign of VE values are a reflection of the type of interaction staking place in the mixture, which are the result of different effects containing the loss of the dipole interaction from each other and the breakdown of the IL ion pair (positive VE). The interaction between the ion pair of ILs increases as compared to IL+2-butanol interactions, which leads to positive contribution at higher temperature over whole composition range. While, Figure 6c shows the negative VE values for TBAH+2-methyl-2-propanol at all measured temperatures over whole composition range. This might be due to the large difference between the molar volumes of the 2-methyl-2-propanol and TPAH implying that it is possibly due to the fact that the relatively small organic molecules fit into the interstices upon mixing. Therefore, the filling effect of organic molecular liquids in the interstices of ILs, and the ion-dipole interactions between organic molecular liquid and alkyl cation of ammonium ILs, all contribute to the negative values of VE.
Clearly, the observed negative VE values increases further with increasing the temperature in the entire mole fraction range for all IL systems. It is interesting to note that the VE values in the ILs+2-methyl-2-propanol mixture shows more negative values of VE than the IL+2-butanol and IL+1-butanol mixtures at 298.15 K over the alcohol rich concentration range (Figure 7), implying that in the 2-methyl-2-propanol there is strong ion-dipole interactions and packing effects with ILs as compared to 2-butanol and 1-butanol. The magnitude and sign of VE values are a reflection of the type of interactions taking place in the mixture, which reveals that VE values are more negative for TBAH (VE = −8.149 cm3.mol−1 at x1 = 0.2831 for TBAH+2-methyl-2-propanoland VE = −5.787 cm3.mol−1 at x1 = 0.5428 for TBAH+2-butanol) than TPAH (VE = −7.547 cm3.mol−1 at x1 = 0.4143 for TPAH+2-methyl-2-propanol and VE = −5.702 cm3.mol−1 at x1 = 0.5634 for TPAH+2-butanol) in all systems except in ILs+1-butanol systems. For TBAH+1-butanols (VE = −2.448 cm3.mol−1 at x1 = 0.5993) and VE = −3.175 cm3.mol−1 at x1 = 0.4139 for TPAH+1-butanol, hence the VE is more negative for the TPAH+1-butanol than TBAH+1-butanol due to steric hindrance created by the long chain cation of TBAH, that reduces the interaction magnitude between TBAH and 1-butanol. Interestingly, the hydrogen bonding between ILs and butanol isomers has predicted using semiempirical calculations with the help of Hyperchem 7, and those interactions are explicitly elucidated in Figures S1, S2, S3, S4, S5, S6.
Using semiempirical calculations for the hydrogen bonding between ILs and butanol isomers, displayed in Figures S1, S2, S3, S4, S5, S6, we calculated heat of formation of the complexes and compared the values with those of the ILs and butanol isomers (as displayed in Table 2). In all the cases, ΔHf of the complex resulting from hydrogen bonding was higher than the sum of ΔHf's of butanol isomers and ILs. It is reasonable to assume that these differences (ΔΔHf), calculated according to equation 3, represent the energies of the hydrogen bond. The energies of the hydrogen bonding can also obtained by using the total binding energies of butanol isomers [54] and ILs, presented in Table 2, instead of ΔHf's for these calculations. The results in Table 2 indicate that the energy required for the formation of a weak hydrogen bond is less than required for the formation of a stronger hydrogen bond:(3)Where ΔHf (1) is the heat of formation of the butanol isomers, ΔHf (2) the heat of formation of the ILs and ΔHf (3) the heat of formation of the complex (butanol isomers and ILs).
A glance at the Figure S1, illustrated the hydrogen bonding between nitrogen group of TPAH IL with the “-OH” group of 1-buatnol. The binding energy of the TPAH is found to be −3784.96 kcal/mol and that of 1-butanol is found to be −1330.63 kcal/mol, but after the hydrogen bonding occur between the TPAH and 1-butanol, the binding energy of complete system comes out to be −5117.47 kcal/mol (Table 2). Hence, the estimated hydrogen bond energy of the above system is ≈1.87 kcal/mol, which could probably due to the interaction of TPAH with 1-butanol that leads to decrease in energy (less than sum of individual energy of TPAH and 1-butanol), and increase in the strength of hydrogen bonding. Similarly, the Figures S2 and S3, depict the possibility of hydrogen bonding between the nitrogen group of TPAH IL with hydroxyl group of 2-butanol and 2-methyl-2-propanol and the estimated hydrogen bond energies are ≈4.14 and ≈0.64 kcal/mol respectively (Table 2). This shows that the strength of hydrogen bonding is more for TPAH+2-butanol as compared to other butanol isomers. Further, the Figures S4 to S6, clearly again show the hydrogen bonding between the nitrogen group of TBAH with hydroxyl group of 1-butanol, 2-butanol and 2-methyl-2-propanol and now the hydrogen bond energies are ≈0.39, ≈2.29 and ≈0.08 respectively. Hence, we may conclude that the hydrogen bond in case of TBAH+2-butanol is stronger as compared to other butanol isomers.
Our interpretation of hydrogen bonding between of IL and butanol isomers (based VE data) is quite corroborated with our theoretical calculation of hydrogen bonding of IL+butanol isomers. It is noteworthy that the hydroxyl groups of alcohols are interacting with the nitrogen group of ILs (TPAH and TBAH) (Figures S1, S2, S3, S4, S5, S6). According to literature, the negative VE values are a result of contributions from both the accommodation of organic molecules in the interstice of the IL networks and the ion–dipole interactions between the organic molecules and cation of the ionic liquid [38], [57]. Our experimental results reveal that the negative VE values for entire composition and theoretical calculation suggested that hydroxyl groups of alcohols are interacting with the cation of ILs, as illustrated in Figures S1, S2, S3, S4, S5, S6. Hence, our results are very well correlated with literature results. Therefore, the absolute value of VE is an indicative to the difference in the packing efficiency and the interaction intensity. As can be seen from Figure 7, the VE values for the studied systems follow the sequence: 2-methyl-2-propanol>1-butanol>2-butanol. If only ion–dipole interactions are taken into consideration, the order 1-butanol>2-butanol>2-methyl-2-propanol is understandable. The decreased dielectric constant from 1-butanol (17.8), 2-butanol (16.6) and 2-methyl-2-propanol (10.9) leads to the weaker ion–dipole interaction and in turn resulting in the smaller VE values. Whereas, if we consider the energies of the hydrogen bond (Table 2), the order for ILs+butanol isomers follows: 2-butanol>1-butanol>2-methyl-2-propanol. Our experimental results suggest the order 2-methyl-2-propanol>1-butanol>2-butanol, which reveal that the interactions are not only due to individual contribution of ion-dipole interaction or H-bonding, but it is the combined effect of both the factors. Whereas, another plausible reason is that the butanol isomers makes it easy to accommodate in the interstice of the IL network, and the higher packing efficiency also leads to the larger VE values. While, with increase in temperature there is decrease in VE values in all the systems because at higher temperature the packing efficiency decreases of ILs. On the other hand, the ion-dipole interactions also decrease with the increase in temperature that leads to decrease in VE values.
Further, for better understanding of the interactions between the tetraalkylammonium hydroxide ILs, we have calculated the Δκs. As seen in Figure 8, Δκs values of tetraalkylammonium hydroxide ILs+butanol isomers are negative over the full composition range at 293.15 K as a function of ILs concentration.The behavior of Δκs, implies that these mixtures are less compressible than the ideal mixture. This is due to closer approach of unlike molecules and a stronger interaction between components of mixtures that leads to a decrease in the compressibility. From Figure 8a, it can be seen that the minimum Δκs values are observed at mole fraction of IL ≈0.4141 for the TPAH+1-butanol system. The negative Δκs values of TPAH+1-butanol are attributed to the strong attractive interactions due to the solvation of the ions in these solvents, over the complete composition range and at all studied temperatures. Similarly, the curves in Figure 8 (c and d), show that the Δκs values for the 2-butanol or 2-methyl-2-propanol+TPAH systems are negative over the complete composition range and at all studied temperatures, except at the higher temperatures (308.15 and 313.15 K) for ≈0.8100 to 0.9999 composition range. The minimum is approached at mole fraction of IL ≈0.4171, ≈0.3272 and ≈0.3138 for the TPAH+1-butanol, TPAH+2-butanol and TPAH+2-methyl-2-propanol systems at all temperatures, respectively. The negative Δκs values attributed to the strong attractive interactions between the molecules of the components. The negative values of Δκs of the TPAH with butanol isomers imply that solvent molecules around solute are less compressible than the solvent molecules in the bulk solutions. Whereas on further addition of IL, there is decrease in the compressibility graph at all studied temperature ranges. This might be due to the decreased attraction between IL and butanol isomers in IL rich concentration region. Additionally, for the 2-butanol and 2-methyl-2-propanol, there are positive Δκs values at higher temperatures, this is might be again due to decrease in attraction of TPAH and alcohol molecules in the IL-rich concentration region, since the interaction between the ILs increases and whereas decreases in case of IL and alcohols.
Figure 9 depicts the negative Δκs values of all TBAH+butanol isomers over the full composition range at 293.15 K. The curves in Figure 9 show that the Δκs values for the 1-butanol or 2-butanol or 2-methyl-2-propanol systems are negative over the complete composition range at low temperature. The minimum is approached at mole fractions of IL ≈0.4787, ≈0.4604 and ≈0.1985 for the TBAH+1-butanol, TBAH+2-butanol and TBAH+2-methyl-2-propanol systems, respectively. Our results show that for all the system, TBAH+butanol isomers shows the positive Δκs values in the IL-rich region at the higher temperatures. These results are very similar with the TPAH+alcohols at higher temperature, which might be due to the decrease in the attraction of TBAH and alcohol molecules in the IL-rich concentration region due to the increased interaction between the ILs and the decreased interaction between IL and alcohols. Obviously, the Δκs values in the ILs+2-methyl-2-propanol mixture shows more negative values of Δκs than the ILs+2-butanol and ILs+1-butanol mixtures at 293.15 K over the entire concentration range (Figure 10), implying that in the 2-methyl-2-propanol there is strong ion-dipole interactions and packing effects with ILs as compared to 2-butanol and 1-butanol. The magnitude and sign of Δκs values are a reflection of the type of interactions taking place in the mixture, which reveals that Δκs values are more negative for TBAH+2-methyl-2-propanol (Δκs = −211.351TPa−1 at x1 = 0.1985), than TBAH+1-butanol (Δκs = −168.519TPa−1 at x1 = 0.4787), and least is for TBAH+2-butanol (Δκs = −144.998TPa−1 at x1 = 0.4604). For TPAH+2-methyl-2-propanol (Δκs = −264.189TPa−1 at x1 = 0.3115), TPAH+1-butanol (Δκs = −222.322TPa−1 at x1 = 0.3218), and TPAH+2-butanol (Δκs = −153.949TPa−1 at x1 = 0.0.3244). TPAH+butanol isomers have more negative Δκs values then TBAH+butanol isomers this is might be due to the steric hindrance created by the long chain cation of TBAH, which reduces the interaction magnitude between TBAH and alcohols.
However, after close look about the physical properties of alcohols with ILs, we observed that hydroxyl position of alcohols are playing important role in addition to the cation chain length of the ILs. ILs (TPAH or TBAH) interact strongly with the 2-methyl-2-propanol as compared to the 2-butanol and 1-butanol, this is might be due to more+I-effect of 2-methyl-2-propanol, which increases its tendency to interact with ILs more strongly as compared to 2-butanol and 1-butanol. Also, 1-butanol interacts more strongly as compared to 2-butanol with ILs, which might be due to the steric hindrance the interaction of 2-butanol decreases as compared to 1-butanol. However, there is decrease in VE and Δκs values in all the system due to increase in temperature; this might be due to strong self-association between the alcohol molecules that prevents the alcohol-IL strong interactions. Our experimental results of VE values are very well supported with literature [38], [46], [47]. Wen-Lu Weng [47], showed the interactions of anisole with 2-butanol and 2-methyl-2-propanol, author observed that VE values of the 2-methyl-2-propanol is more negative than 2-butanol. Additionally, the VE values increases (less negative) with increase in temperature. Further, Qian et al. [38], showed that 1-methylimidazolium acetate IL interacts with methanol, ethanol, 1-propanol and 1-butanol, it was observed that in all the systems VE values increase with increase in the temperature. Moreover, during the interaction of formamide with 1-butanol negative VE values have been observed, whereas interaction of formamide with 2-butanol results in positive VE values [46]. These all results by various authors support our above results explanation that interactions between ILs+alcohols depend upon the position of hydroxyl group. Therefore, the physicochemical properties of ILs are quite sensitive toward the structure and nature of interacting molecules.
Conclusion
We have performed and compared thermophysical properties of binary mixtures of tetraalkylammonium hydroxide based ILs with butanol isomers over the whole composition range at various temperatures (293.15 to 313.15 K, in steps of 5 K). To obtain a more detailed picture of the molecular interactions, we measured temperature dependence properties of ρ and u for ILs with butanol isomers over the whole composition range at various temperatures. The ρ and u values increase with the increasing the cation alkyl chain length of ILs. Our results reveal that the position of hydroxyl group in alcohols leads to alteration of the thermophysical properties of ILs. To measure the non-ideality of the mixtures, we determined VE and Δκs at each temperature as a function of IL concentration. The predicted properties were correlated by the Redlich-Kister type equation. Our studies demonstrate that there is decrease in VE and Δκs values in all the systems due to increase in temperature; this might be due to strong self-association between the alcohol molecules that prevents the alcohol-IL strong interactions. Additionally, according to the theoretical calculations obtained by HyperChem 7, the energy of hydrogen bond is more for low alkyl chain ILs (TPAH) as compared to higher alkyl chain ILs (TBAH) with alcohols. Molecular interactions such as ion-dipole and hydrogen bonding between the butanol isomers and alkyl chain of ILs are suggested to be mainly responsible for variation in the thermophysical parameters. Our findings provide better molecular interactions for the mixing of the solvents and better analysis of the solvation process.
Supporting Information
Figure S1.
Schematic depiction of the hydrogen bonding interaction between TPAH and 1-butanol molecules, which is predicted by a semiempirical calculation with the help of HyperChem 7.
https://doi.org/10.1371/journal.pone.0086530.s001
(TIF)
Figure S2.
Schematic depiction of the hydrogen bonding interaction between TPAH and 2-butanol molecules, which is predicted by a semiempirical calculation with the help of HyperChem 7.
https://doi.org/10.1371/journal.pone.0086530.s002
(TIF)
Figure S3.
Schematic depiction of the hydrogen bonding interaction between TPAH and 2-methyl-2-propanol molecules, which is predicted by a semiempirical calculation with the help of HyperChem 7.
https://doi.org/10.1371/journal.pone.0086530.s003
(TIF)
Figure S4.
Schematic depiction of the hydrogen bonding interaction between TBAH and 1-butanol molecules, which is predicted by a semiempirical calculation with the help of HyperChem 7.
https://doi.org/10.1371/journal.pone.0086530.s004
(TIF)
Figure S5.
Schematic depiction of the hydrogen bonding interaction between TBAH and 2-butanol molecules, which is predicted by a semiempirical calculation with the help of HyperChem 7.
https://doi.org/10.1371/journal.pone.0086530.s005
(TIF)
Figure S6.
Schematic depiction of the hydrogen bonding interaction between TBAH and 2-methyl-2-propanolmolecules, which is predicted by a semiempirical calculation with the help of HyperChem 7.
https://doi.org/10.1371/journal.pone.0086530.s006
(TIF)
File S1.
Supporting Tables. Table S1. Mole fraction (x1) of IL, density (ρ), ultrasonic sound velocity (u), Excess molar volumes (VE), isentropic compressibility (κs), and deviation in isentropic compressibility (Δκs) for the systems of tetraalkylammonium hydroxide IL with butanol isomers at T = 293.15, 298.15, 303.15, 308.15 and 313.15 K and at atmospheric pressure. Table S2. Estimated Parameters of eq 1 and Standard Deviation σ, for the Systems of ILs with butanol isomers as Function of Temperature.
https://doi.org/10.1371/journal.pone.0086530.s007
(DOCX)
Author Contributions
Conceived and designed the experiments: PA. Performed the experiments: PA. Analyzed the data: PA PV. Contributed reagents/materials/analysis tools: KYB ITK EHC. Wrote the paper: PA.
References
- 1. Janikowski J, Razali MR, Forsyth CM, Nairn KM, Batten SR, et al. (2013) Physical Properties and Structural Characterization of Ionic Liquids and Solid Electrolytes Utilizing the Carbamoyl cyano(nitroso) methanide Anion. Chem Plus Chem 78: 486–497.
- 2. Capelo SB, Méndez-Morales T, Carrete J, Lago EL, Vila J, et al. (2012) Effect of Temperature and Cationic Chain Length on the Physical Properties of Ammonium Nitrate-Based Protic Ionic Liquids. J Phys Chem B 116: 11302–11312.
- 3. Liu ZP, Wu XP, Wang WC (2006) A Novel United-Atom Force Field for Imidazolium-Based Ionic Liquids. Phys Chem Chem Phys 8: 1096–1104.
- 4. Chen T, Chidambaram M, Liu Z, Berend SB, Bell AT (2010) Viscosities of the Mixtures of 1-Ethyl-3-Methylimidazolium Chloride with Water, Acetonitrile and Glucose: A Molecular Dynamics Simulation and Experimental Study. J Phys Chem B 114: 5790–5794.
- 5. Couadou E, Jacquemin J, Galiano H, Hardacre C, Anouti MA (2013) Comparative Study on the Thermophysical Properties for Two Bis[(trifluoromethyl) sulfonyl] imide-Based Ionic Liquids Containing the Trimethyl-Sulfonium or the Trimethyl-Ammonium Cation in Molecular Solvents. J Phys Chem B 117: 1389–1402.
- 6. Fox ET, Paillard E, Borodin O, Henderson WA (2013) Physicochemical Properties of Binary Ionic Liquid–Aprotic Solvent Electrolyte Mixtures. J Phys Chem C 117: 78–84.
- 7. Zhang S, Sun N, He X, Lu X, Zhang X (2006) Physical Properties of Ionic Liquids: Database and Evaluation. J Phys Chem Ref Data 35: 1475–1517.
- 8. González B, Calvar N, González E, Domínguez A (2008) Density and Viscosity Experimental Data of the Ternary Mixtures 1-Propanol or 2-Propanol+Water+1-Ethyl-3-Methylimidazolium Ethylsulfate. Correlation and Prediction of Physical Properties of the Ternary Systems. J Chem Eng Data 53: 881–887.
- 9. Pereiro AB, Legido JL, Rodríguez A (2007) Physical Properties Of Ionic Liquids Based on 1-Alkyl-3-Methylimidazolium Cation and Hexafluorophosphate as Anion and Temperature Dependence. J Chem Thermodynamics 39: 1168–1175.
- 10. Annat G, Forsyth M, MacFarlane DR (2012) Ionic Liquid Mixtures Variations in Physical Properties and Their Origins in Molecular Structure. J Phys Chem B 116: 8251–8258.
- 11. Deetlefs M, Seddon KR, Shara M (2006) Predicting Physical Properties of Ionic Liquids. Phys Chem Chem Phys 8: 642–649.
- 12. Kavitha T, Attri P, Venkatesu P, Ramadevi RS, Hofman T (2012) Influence of Alkyl Chain Length and Temperature on Thermophysical Properties of Ammonium-Based Ionic Liquids with Molecular Solvent. J Phys Chem B 116: 4561–4574.
- 13. Attri P, Venkatesu P, Hofman T (2011) Temperature Dependence Measurements and Structural Characterization of Trimethyl Ammonium Ionic Liquids with a Highly Polar Solvent. J Phys Chem B 115: 10086–10097.
- 14. Attri P, Venkatesu P, Kumar A (2010) Temperature Effect on the Molecular Interactions Between Ammonium Ionic Liquids and N,N-Dimethylformamide. J Phys Chem B 114: 13415–13425.
- 15. Attri P, Reddy PM, Venkatesu P, Kumar A, Hofman T (2010) Measurements and Molecular Interactions for N,N-Dimethylformamide with Ionic Liquid Mixed Solvents. J Phys Chem B 114: 6126–6133.
- 16. Hou M, Xu Y, Han Y, Chen B, Zhang W, et al. (2013) Thermodynamic Properties of Aqueous Solutions of Two Ammonium-Based Protic Ionic Liquids at 298.15 K. J Mol Liqs 178: 149–155.
- 17. Attri P, Venkatesu P, Kumar A (2011) Activity and Stability of α-Chymotrypsin in Biocompatible Ionic Liquids: Enzyme Refolding by Triethyl Ammonium Acetate. Phys Chem Chem Phys 13: 2788–2796.
- 18. Seddon KR (1997) Ionic Liquids for Clean Technology. J Chem Technol Biotechnol 68: 351–356.
- 19. Greaves TL, Drummond C (2008) Protic Ionic Liquids: Properties and Applications. J Chem Rev 108: 206–237.
- 20. Rogers RD, Seddon KR (2003) Ionic Liquids-Solvents of the Future?. Science 302: 792–793.
- 21. Davis JH (2004) Task-Specific Ionic Liquids. Chem Lett 33: 1072–1077.
- 22. Attri P, Venkatesu P, Kumar A (2011) Activity and Stability of α-Chymotrypsin in Biocompatible Ionic Liquids: Enzyme Refolding by Triethyl Ammonium Acetate. Phys Chem Chem Phys 13: 2788–2796.
- 23. Chiappe C, Pomelli CS, Rajamani S (2011) Influence of Structural Variations in Cationic and Anionic Moieties on the Polarity of Ionic Liquids. J Phys Chem B 115: 9653–9661.
- 24. Attri P, Venkatesu P, Kumar A (2012) Water and a Protic Ionic Liquid Acted as Refolding Additives for Chemically Denatured Enzymes. Org Biomol Chem 10: 7475–7478.
- 25. Noda A, Bin Hasan Susan A, Kudo K, Mitsushima S, Hayamizu K, et al. (2003) Brønsted Acid–Base Ionic Liquids as Proton-Conducting Nonaqueous Electrolytes. J Phys Chem B 107: 4024–4033.
- 26. Attri P, Lee SH, Hwang SW, IL Kim J, Lee SW, et al. (2013) Elucidating Interactions and Conductivity of Newly Synthesised Low Bandgap Polymer with Protic and Aprotic Ionic Liquids. PloS One 2013, 88: e68970
- 27. Attri P, Choi EH (2013) Influence of Reactive Oxygen Species on the Enzyme Stability and Activity in the Presence of Ionic Liquids. PloS One 2013, 8: e75096
- 28. Reddy PM, Venkatesu P (2011) Ionic Liquid Modifies the Lower Critical Solution Temperature (LCST) of Poly(N-isopropylacrylamide) in Aqueous Solution. J Phys Chem B 115: 4752–4757.
- 29. Abareshi M, Goharshadi EK, Zebarjad SM (2009) Thermodynamic Properties of the Mixtures of Some Ionic Liquids with Alcohols using a Simple Equation of State. J Mole Liq 149: 66–73.
- 30. Arce A, Rodil E, Soto A (2006) Physical and Excess Properties for Binary Mixtures of 1-Methyl-3-Octylimidazolium Tetrafluoroborate, [Omim][BF4], Ionic Liquid with Different Alcohols. J Solution Chem 35: 63–78.
- 31. González EJ, González B, Macedo EA (2013) Thermophysical Properties of the Pure Ionic Liquid 1-Butyl-1-methylpyrrolidinium Dicyanamide and Its Binary Mixtures with Alcohols. J Chem Eng Data 58: 1440–1448.
- 32. González EJ, González B, Calvar N, Domínguez A (2007) Physical Properties of Binary Mixtures of the Ionic Liquid 1-Ethyl-3-Methylimidazolium Ethylsulfate with Several Alcohols at T = (298.15, 313.15 K, and 328.15) K and Atmospheric Pressure. J Chem Eng Data 52: 1641–1648.
- 33. Lehmann J, Rausch MH, Leipertz A, Fröba AP (2010) Density and Excess Molar Volumes for Binary Mixtures of Ionic Liquid 1-Ethyl-3-Methylimidazolium Ethylsulfate with Solvents. J Chem Eng Data 55: 4068–4074.
- 34. Vercher E, Orchillés AV, Miguel PJ, Martínez-Andreu A (2007) Volumetric and Ultrasonic Studies of 1-Ethyl-3-Methylimidazolium Trifluoromethane sulfonate Ionic Liquid with Methanol, Ethanol, 1- Propanol, and Water at Several Temperatures. J Chem Eng Data 52: 1468–1482.
- 35. Rilo E, Ferreira AGM, Fonseca IMA, Cabeza O (2010) Densities and Derived Thermodynamic Properties of Ternary Mixtures 1-Butyl-3-Methyl-Imidazolium Tetrafluoroborate+Ethanol+Water at Seven Pressures and Two Temperatures. Fluid Phase Equilib 296: 53–59.
- 36. Anouti M, Jacquemin J, Lemordant D (2010) Volumetric Properties, Viscosities, and Isobaric Heat Capacities of Imidazolium Octanoate Protic Ionic Liquid in Molecular Solvents. J Chem Eng Data 55: 5719–5728.
- 37. Alvarez VH, Mattedi S, Martin-Pastor M, Aznar M, Iglesias M (2011) Thermophysical Properties of Binary Mixtures of {Ionic Liquid 2-Hydroxy Ethylammonium Acetate+(Water, Methanol, or Ethanol)}. J Chem Thermodyn 43: 997–1010.
- 38. Qian W, Xu Y, Zhu H, Yu C (2012) Properties of Pure 1-Methylimidazolium Acetate Ionic Liquid and its Binary Mixtures with Alcohols. J Chem Thermodyn 49: 87–94.
- 39. Gómez E, Calvar N, Macedo EA, Domínguez A (2012) Effect of the Temperature on the physical Properties of Pure 1-Propyl-3- MethylimidazoliumBis(Trifluoromethylsulfonyl)Imide and Characterization of its Binary Mixtures with Alcohols. J Chem Thermodyn 45: 9–15.
- 40. Jiang H, Wang J, Zhao F, Qi G, Hu Y (2012) Volumetric and Surface Properties of Pure Ionic Liquid n-Octylpyridinium Nitrate and its Binary Mixtures with Alcohols. J Chem Thermodyn 47: 203–208.
- 41. Domanska U, Zawadzki M, Lewandrowska A (2012) Effect of Temperature and Composition on the Density, Viscosity, Surface Tension, and Thermodynamic Properties of Binary Mixtures of N-OctylisoquinoliniumBis(Trifluoromethylsulfonyl)Imide with Alcohols. J Chem Thermodyn 48: 101–111.
- 42.
Riddick JA, Bunger WB, Sakano TK (1986) Organic Solvents 4th ed., Wiley-Inter Science, New York.
- 43. Mokhtarani B, Sharifi A, Mortaheb HR, Mirzaei M, Mafi M, et al. (2009) Density and Viscosity of 1-Butyl-3-Methylimidazolium Nitrate with Ethanol, 1-Propanol, or 1-Butanol at Several Temperatures. J Chem Thermodyn 41: 1432–1438.
- 44. De Cominges BE, Piñeiro MM, Mosteiro L, Mascato1 E, Mato MM, et al. (2002) Temperature Dependence of Thermophysical Properties of Octane+1-Butanol System. J Therm Anal Calorim 70: 217–227.
- 45. Gnanakumari P, Venkatesu P, Mohan KR, Rao MVP, Prasad DHL (2007) Excess Volumes and Excess Enthalpies of N-Methyl-2-Pyrrolidone with Branched Alcohols. Fluid Phase Equilibria 252: 137–142.
- 46. Nain AK (2007) Densities and Volumetric Properties of Binary Mixtures of Formamide with 1-Butanol, 2-Butanol, 1,3-Butanediol and 1,4-Butanediol at Temperatures between 293.15 and 318.15 K. J Solution Chem 36: 497–516.
- 47. Weng WL (1999) Densities and Viscosities for Binary Mixtures of Anisole with 2-Butanol, 2-Methyl-1-propanol, and 2-Methyl-2-propanol. J Chem Eng Data 44: 788–791.
- 48. Bravo-Sanchez MG, Iglesias-Silva GA, Estrada-Baltazar A (2010) Densities and Viscosities of Binary Mixtures of n-Butanol with 2-Butanol, Isobutanol, and tert-Butanol from (303.15 to 343.15) K. J Chem Eng Data 55: 2310–2315.
- 49. Łachwa J, Morgado P, Esperancüa José MSS, Guedes HJR, Lopes JNC, et al. (2006) Fluid-Phase Behavior of {1-Hexyl-3-methylimidazolium Bis(trifluoromethylsulfonyl) Imide, [C6mim][NTf2], +C2-C8 n-Alcohol} Mixtures: Liquid-Liquid Equilibrium and Excess Volumes. J Chem Eng Data 51: 2215–2221.
- 50. Domınguez M, Gascon I, Valen A, Royo FM, Urieta JS (2000) Densities of (2-Butanol+n-Hexane+1-Butylamine) at T = 298∶15 And T = 313∶15 K: Excess and Partial Excess Molar Volumes and Application of the Eras Model. J Chem Thermodyn 32: 1551–1568.
- 51. Huq F, Yu JQ (2002) Molecular Modeling Analysis: “Why is 2-Hydroxypyridine Soluble in Water but not 3-Hydroxypyridine?”. J Mol Model 8: 81–86.
- 52. Jorgensen WL, Chandrasekhas J, Madura JD, Impey RW, Klein ML (1983) Comparison of Simple Potential Functions for Simulating Liquid Water. J Chem Phys 79: 926–935.
- 53. Dewar MJS, Storch DM (1985) Development and Use of Quantum Molecular Models. 75. Comparative Tests of Theoretical Procedures for Studying Chemical Reactions. J Am Chem Soc 107: 3898–3902.
- 54. Pop E, Brewster ME (1997) Dimerization of Dexanabinol by Hydrogen Bonding Accounts for Its Hydrophobic Character. Intern J Quantum Chem 65: 1057–1064.
- 55. Wilkes JS (2004) Properties of Ionic Liquid Solvents for Catalysis. J Mol Catal A: Chem 214: 11–17.
- 56. Marsh KN, Boxall JA, Lichtenthaler R (2004) Room Temperature Ionic Liquids and Their Mixtures—a Review. Fluid Phase Equilib 219: 93–98.
- 57. Zhu A, Wang J, Liu R (2011) A Volumetric and Viscosity Study for the Binary Mixtures of 1-Hexyl-3-Methylimidazolium Tetrafluoroborate with some Molecular Solvents. J Chem Thermodyn 43: 796–799.