- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Advances in Physical Chemistry
Volume 2012 (2012), Article ID 829523, 7 pages
Conformational Analysis of Quaternary Ammonium-Type Ionic Liquid Cation, N,N-Diethyl-N-methyl-N-(2-methoxyethyl) Ammonium Cation
1Department of Applied Chemistry, National Defense Academy, 1-10-20, Hashirimizu, Yokosuka, Kanagawa 239-8686, Japan
2Department of Materials Science and Engineering, National Defense Academy, 1-10-20, Hashirimizu, Yokosuka, Kanagawa 239-8686, Japan
Received 17 April 2012; Accepted 13 June 2012
Academic Editor: Ranko Richert
Copyright © 2012 Takahiro Takekiyo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Conformational preference of N,N-diethyl-N-methyl-N-(2-methoxyethyl) ammonium cation ([DEME]+), which is one of the quaternary ammonium-based ionic liquid cation, in the gas phase has been investigated using a density functional theory (DFT) calculation. Eight candidates for the stable conformers of [DEME]+ exist in the gas phase, and can it energetically classify into two groups. One is a five conformers group, which has the intramolecular attractive interaction form (the folded form). The other is a three conformers group, which is the noninteraction form (the extended form). The transformation from the folded form to the extended form induces large changes in the dipole moment and partial charges of N and O atoms. Here we show that the difference in the dipole moment and partial charges of N and O atoms associated with the conformational change of [DEME]+ are closely related to the molecular orientation of [DEME]-based ionic liquids in the liquid state.
The conformational analysis of room temperature ionic liquids (RTILs) using a vibrational spectroscopy combined with the quantum chemical calculation has been studied because the investigation of the conformational equilibria in RTILs is a good indicator of the molecular orientation for the phase transition and solution structure of RTILs [1–3].
Recently, the phase behaviors and solution structures of N,N-diethyl-N-methyl-N-(2-methoxyethyl) ammonium tetrafluoroborate ([DEME][BF4]) and iodide ([DEME][I]), which are typical quaternary ammonium-based RTILs, in the pure solution and also in the aqueous solutions have been investigated using various experimental techniques such as Raman, deferential thermal analysis (DTA), and X-ray diffraction methods [4–12]. Imai et al.  and Abe et al.  reported the phase diagram of [DEME][BF4]-water mixed solutions as functions of temperature and water concentration. They found that [DEME][BF4]-water mixed solutions show complicated phase behavior. Moreover, Yoshimura et al.  found the existence of the nearly free hydrogen bonds of water in the [DEME][BF4]-water mixed solutions at 298 and 77 K using a Raman spectroscopy. In this relation, Imai et al. [9, 12] reported that [DEME] ( = BF4 and I)-water mixed solutions have double glass transition temperatures () due to the liquid-liquid immiscibility using a DTA method. It is, however, still unclear the relationship between the molecular orientation such as conformation and the complicated phase behavior in the aqueous [DEME](= BF4 and I) solutions.
Although the conformations of 1-alkyl-3-methylimidazolium-based RTILs have been well studied [2, 3, 13–18], there has been no report on the conformation of [DEME]+. As one of the reasons, many free internal rotations in the [DEME]+ complicates the conformation of [DEME]+ in the liquid states rather than those of 1-ethyl-3-methylimidazolium cation ([emim]+) and 1-butyl-3-methylimidazolium cation ([bmim]+). Clear identification of the stable conformers of [DEME]+ in the liquid state is difficult, but, as a first step, it is important to reveal the stable conformers of [DEME]+ in the gas phase for determining the complicated conformation in the liquid state.
In this study, we have attempted to determine the conformational preference of [DEME]+ in the gas phase using a density functional theory (DFT) calculation. Our calculation results show that [DEME]+ in the isolated system classifies into the intramolecular attractive interaction form (the folded form) and the noninteraction form (the extended form). These forms have large differences in the dipole moments and partial charges.
Density functional theory (DFT) calculations were carried out using GAUSSIAN03 program . For the present calculation of [DEME]+, we used the Becke’s three-parameter (B3) exchange function . The B3 exchange function was combined with Lee-Yang-Peer correlation function (B3LYP) . Calculations of geometries, total energies, and atomic charges of [DEME]+ were performed using the 6–311++G(d,p) basis set. The normal coordinate calculations using the B3LYP method were performed by uniform wavenumber scaling factor of 0.9613. The basis set superposition errors (BSSEs) were confirmed to be negligible by the counterpoise method [22, 23]. The initial structure of [DEME]+ was started from the all conformer shown in Figure 1. The Breneman and Wiberg CHELPG scheme [24, 25] was used to determine the atomic charges from the electrostatic potentials.
3. Results and Discussion
We consider the conformations on (i) the CN+CCOC chain and (ii) the ethyl-chain (CCN+C) to determine the conformational preference of [DEME]+. First, we focused on the conformational preference of the CN+CCO−C chain. The changes in potential energy surfaces (PES) for the CN+CC, N+CCO−, and CCO−C angles were performed by the B3LYP/6-311++G(d,p) basis set as shown in Figure 2. As seen in Figure 2(a), the CN+CC angle has double energy minima at ±30°. The energy differences () between at ±180° and ±30° of the CN+CC angle showed a large value (over +20 kJ mol−1). The CN+CC angle energetically prefers to the angle at 180° (the form). A remarkable result is the change in the PES for the N+CCO− angle. The N+CCO− angle showed the smallest PES change among three angles, and has double energy minima at ±85°. between at ±180° and ±85° (the form) of the N+CCO− angle are about −1.3 kJ/mol, and are below a RT order (2.48 kJ/mol). The N+CCO− angle takes the form (180°) and the form (±85°), and the form is energetically more stable than the form. On the other hand, the CCO−C angle having the nonenergy minimum only takes the form because the PES shows the lowest at 180°. On the basis of these results, the CN+CC and N+CCO− angles may be important key to determine the conformation of the CN+CCO−C chain. Therefore, the CN+CCO−C chain takes the () and () conformers.
Here, we must discuss why the form of the N+CCO− angle is stable. One possible explanation is the electron polarizabilities of N and O atoms. According to the DFT and Raman studies by Yoshida et al. [26, 27], the OCCO angle of 1,2-dimethoxyethane (DME: CH3OCH2CH2OCH3) prefers the form to the form owing to the electron polarizability of O atoms. Thus, we presume that the higher stability of the form for the N+CCO− angle of [DEME]+ might arise from the intramolecular attractive interaction due to the electron polarizabilities of N and O atoms because the structure of N+COO− angle of [DEME]+ resembles that of the OCCO angle of DME.
To assess this speculation, we estimated the partial atomic charges ( and ) of the N and O atoms, the distance () between the N and O atoms, and dipole moment (μ) for the [DEME]+ conformers using the B3LYP/6-311++G(d,p) basis set. When the N+CCO− angle takes the form, the and values are +0.19 and −0.40, respectively, and the and μ values are 3.72 Å and 4.12 D, respectively. On the other hand, when the N+CCO− angle takes the form, the and values are about +0.12 and −0.30, respectively, and the and μ values 3.06 Å, and 2.67 D. Thus, with transforming from the to forms, the and values become smaller and the value of becomes shorter. The μ value becomes larger ( D). Accordingly, we conclude that the intramolecular attractive interaction due from the electron polarizabilities of N and O atoms is the main cause for the higher stability of the N+CCO− angle of [DEME]+.
Next, we must consider an affection from the ethyl-chain for the conformational stability of the CN+CCO−C chain since the steric effect owing to the conformational preference of ethyl-chain group would change the conformational stability of the CN+CCO−C chain. In an effort to ensure the affection from the ethyl-chain for the conformational stability of the CN+CCO−C chain, we calculated the PES of CCN+C dihedral angles of the ethyl-chain. Seven conformers of ethyl-chain are found to be energetically stable; the , , , , , , and conformers. Here, the large letters of “”, “”, and “” represent the , (+: the positive angle for the form) and (−: the negative angle for the form) forms, respectively. From the results of ( = , , , , , and ) as shown in Table 1, the and conformers have higher energetic stability rather than other conformers. Previous quantum chemistry calculations indicated that the form of tetraethylammonium cation (Et4N+) having the similar structure of [DEME]+ is less stable than the form due to the steric effect in the gas phase [28, 29]. The energy calculation of the ethyl-chain conformers of [DEME]+ qualitatively agrees with the quantum chemistry calculation results of Et4N+ [28, 29].
We checked the conformational preference of the CN+CCO−C chain again, when the ethyl-chain takes the conformer (Figure 3). Compared with Figures 2 and 3, the stable conformers of the N+CCO− and CCO−C angles are independent of the ethyl-chain conformers. On the other hand, the change in the PES of CN+CC angle is largely dependent of the ethyl-chain conformers since the position of the CN+CC angle is very close with the ethyl-chain. Thus, transformation from the conformer to conformer of the ethyl-chain causes the change in double energy minima of the CN+CC angle. Especially, the energy minimum of the CN+CC angle at +30° moved to the angle at +60° (the form) by a steric effect between the ethyl-chain and the CN+CCO−C chain, and between at 180° and +60° for the CN+CC angle is +8.6 kJ/mol. The and forms of CN+CC angle are energetically stable when the ethyl-chain takes the conformer. Therefore, when the CN+CC angle takes the conformer, the CN+CCO−C chain takes the six conformers; the ,,,,, and conformers, respectively.
In view of the results from the conformational combinations of the CN+CCO−C chain and ethyl-chain, [DEME]+ can take eight conformers in the gas phase as shown in Figure 4. The energy difference and structural parameters for eight conformers of [DEME]+ were done as shown in Table 2. The rank ordering of the conformational stability of [DEME]+ is > > > > ≫ > >. (Here the large and small letters (-) represent the conformers for the ethyl-chain () and the CN+CCO−C chain (), resp.) A notable point is the balance of among the conformers. Eight conformers can energetically classify into two types. One is a group taking the form for the N+CCO− angle (, , and ). The other is a group taking the form for the N+CCO− angle (, , , , and ). Probably, the higher stability of five conformers results from the intramolecular attractive interaction between the N and O atoms. Interestingly, this classification is independent on the ethyl-chain conformer. The former group structurally may be depicted as an extended form and the later group as a folded form. Therefore, [DEME]+ in the isolated system has the intramolecular attractive interaction form (the folded form) and the noninteraction form (the extended form).
Another notable point is the difference in the average dipole moments of the folded and extended forms of [DEME]+. Here, we mention the dipole moment (μ), which is expressed by the product of the total charges and distance between the positive and negative charges, of [DEME]+. [DEME]+ is a monovalent cation but has a dipole moment (μ) due to the electron polarizabilities of N and O atoms in the N+CCO−C chain. As mentioned above, [DEME]+ takes many conformers for the N+CCO−C chain. Thus, the electron polarizabilities of N and O atoms and distance between the N and O atoms of each conformation are different as shown in Table 2. Consequently, each conformer has a different dipole moment ( Debye (D)) despite that [DEME]+ is a monovalent cation. According to previous calculation results, two conformers of other cations (e.g., [emim]+ and [bmim]+) have a different dipole moment [15, 18]. The values of μ for [emim]+ are 1.13 D (planar) and 1.64 D (nonplanar), and those for [bmim]+ are 5.76 D () and 4.96 D (), respectively. Interestingly, the μ values of [DEME]+ showed the same order with those of imidazolium cations.
The difference in the average dipole moments (Δμ) between the folded and extended forms is 1.45 D. The values of Δμ between the planar and nonplanar conformers for [emim]+  and between the and conformers of [bmim]+  are about 0.51 and 0.80 D, respectively. Δμ among the conformers of [DEME]+ obtained by the present calculation is larger than those of other imidazolium cations such as [emim]+ and [bmim]+. This is clearly due to the electron polarizabilities of N and O atoms in the CN+CCO−C chain. Our calculation results indicated that the electron polarizabities of the N and O atoms in the CN+CCO−C chain are the most important key for determining the intermolecular interaction between the [DEME]+ and other molecular liquids or RTILs in the liquid state. In other words, the dipole moments and partial charges of N and O atoms associated with the conformational change of [DEME]+ have close connection with various interesting phenomena such as complicated phase transitions, in [DEME] ( = BF4 and I)-water mixed solutions [4–12].
Finally, we mention the relation of the present calculation of [DEME]+ monomer in the gas phase and in the liquid state such as [DEME]( = BF4 and I). The average between the folded and extended forms of [DEME]+ is about 10 kJ/mol, and this energy corresponds to the energy of the hydrogen-bonded formation of molecules in the liquid state . Therefore, the energy balance between the conformational equilibrium of [DEME]+ in the liquid state may be different with that in the gas phase. It is known that interactions in RTILs are not weak, but strong and mostly electrostatic . Actually some researchers have mentioned that around 50% to 15% of ions in RTILs exist as ion pairs [32, 33]. Thus, the conformation of [DEME]+ in the liquid state is more complicated and the calculations might include larger size clusters such as ion pair, dimmer.
In summary, conformational preference of [DEME]+ in the N,N-diethyl-N-methyl-N-(2-methoxyethyl) ammonium-based RTIL in the gas phase has been investigated by a DFT calculation. Present calculations showed that the [DEME]+ consists of eight candidates for the conformers in the gas phase. We also found that [DEME]+ has an equilibrium between the intramolecular attractive interaction form (the folded from) and the noninteraction form (the extended form) in the gas phase. We found that the CN+CCO−C chain is the most important key for determining the molecular orientation of [DEME]+ in the liquid state. The present calculations determined the structural parameters and conformational identification of the complicated cation at the B3LYP level. Although a set of gas phase calculation on a single ion is simple, but this work may prove important to the community of force field developers. Further study such as larger RTIL cluster size calculation including the ion par and dimmer may facilitate to understand more details of the complicated conformations of [DEME]+ in the liquid state.
- H. O. Hamaguchi and R. Ozawa, “Structure of ionic liquids and ionic liquid compounds: are ionic liquids genuine liquids in the conventional sense?” Advances in Chemical Physics, vol. 131, pp. 85–104, 2005.
- T. Endo, T. Kato, K. I. Tozaki, and K. Nishikawa, “Phase behaviors of room temperature ionic liquid linked with cation conformational changes: 1-Butyl-3-methylimidazolium Hexafluorophosphate,” Journal of Physical Chemistry B, vol. 114, no. 1, pp. 407–411, 2010.
- R. W. Berg, “Raman spectroscopy and ab-initio model calculations on ionic liquids,” Monatshefte fur Chemie, vol. 138, no. 11, pp. 1045–1075, 2007.
- Y. Imai, H. Abe, T. Goto, Y. Yoshimura, Y. Michishita, and H. Matsumoto, “Structure and thermal property of N,N-diethyl-N-methyl-N-2-methoxyethyl ammonium tetrafluoroborate-H2O mixtures,” Chemical Physics, vol. 352, no. 1–3, pp. 224–230, 2008.
- Y. Imai, H. Abe, T. Goto, Y. Yoshimura, S. Kushiyama, and H. Matsumoto, “Orientational ordering of crystal domains in ionic liquid based mixtures,” Journal of Physical Chemistry B, vol. 112, no. 32, pp. 9841–9846, 2008.
- Y. Imai, H. Abe, and Y. Yoshimura, “X-ray diffraction study of ionic liquid based mixtures,” Journal of Physical Chemistry B, vol. 113, no. 7, pp. 2013–2018, 2009.
- Y. Yoshimura, T. Goto, H. Abe, and Y. Imai, “Existence of nearly-free hydrogen bonds in an ionic liquid, N,N -diethyl-N-methyl-N-(2-methoxyethyl) ammonium tetrafluoroborate-water at 77 K,” Journal of Physical Chemistry B, vol. 113, no. 23, pp. 8091–8095, 2009.
- H. Abe, Y. Yoshimura, Y. Imai, T. Goto, and H. Matsumoto, “Phase behavior of room temperature ionic liquid—H2O mixtures: N, N-diethyl-N-methyl-N-2-methoxyethyl ammonium tetrafluoroborate,” Journal of Molecular Liquids, vol. 150, no. 1–3, pp. 16–21, 2009.
- Y. Imai, H. Abe, T. Miyashita, T. Goto, H. Matsumoto, and Y. Yoshimura, “Two glass transitions in N, N-diethyl-N-methyl-N-(2-methoxyethyl) ammonium tetrafluoroborate-H2O mixed solutions,” Chemical Physics Letters, vol. 486, no. 1-3, pp. 37–39, 2010.
- H. Abe, Y. Imai, T. Takekiyo, and Y. Yoshimura, “Deuterated water effect in a room temperature ionic liquid: N,N-diethyl-N-methyl-2-methoxyethyl ammonium tetrafluoroborate,” Journal of Physical Chemistry B, vol. 114, no. 8, pp. 2834–2839, 2010.
- H. Abe, Y. Imai, T. Goto et al., “Water network in room-temperature ionic liquid: N,N-diethyl-N-methyl-N-2-methoxyethyl ammonium tetrafluoroborate,” Metallurgical and Materials Transactions A, vol. 41, no. 5, pp. 1137–1143, 2010.
- Y. Imai, H. Abe, H. Matsumoto, O. Shimada, T. Hanasaki, and Y. Yoshimura, “Glass transition behaviour of the quaternary ammonium type ionic liquid, mixtures,” Journal of Chemical Thermodynamics, vol. 43, no. 3, pp. 319–322, 2011.
- R. Ozawa, S. Hayashi, S. Saha, A. Kobayashi, and H. O. Hamaguchi, “Rotational isomerism and structure of the 1-butyl-3-methylimidazolium cation in the ionic liquid state,” Chemistry Letters, vol. 32, no. 10, pp. 948–949, 2003.
- K. Fujii, T. Fujimori, T. Takamuku, R. Kanzaki, Y. Umebayashi, and S. I. Ishiguro, “Conformational equilibrium of bis(trifluoromethanesulfonyl) imide anion of a room-temperature ionic liquid: raman spectroscopic study and DFT calculations,” Journal of Physical Chemistry B, vol. 110, no. 16, pp. 8179–8183, 2006.
- Y. Umebayashi, T. Fujimori, T. Sukizaki et al., “Evidence of conformational equilibrium of 1-ethyl-3-methylimidazolium in its ionic liquid salts: raman spectroscopic study and quantum chemical calculations,” Journal of Physical Chemistry A, vol. 109, no. 40, pp. 8976–8982, 2005.
- J. C. Lassègues, J. Grondin, R. Holomb, and P. Johansson, “Raman and ab initio study of the conformational isomerism in the 1-ethyl-3-methyl-imidazolium bis(trifluoromethanesulfonyl)imide ionic liquid,” Journal of Raman Spectroscopy, vol. 38, no. 5, pp. 551–558, 2007.
- Y. Umebayashi, T. Mitsugi, S. Fukuda et al., “Lithium ion solvation in room-temperature ionic liquids involving bis(trifluoromethanesulfonyl) imide anion studied by Raman spectroscopy and DFT calculations,” Journal of Physical Chemistry B, vol. 111, no. 45, pp. 13028–13032, 2007.
- R. Holomb, A. Martinelli, I. Albinsson, J. C. Lassègues, P. Johansson, and P. Jacobsson, “Ionic liquid structure: the conformational isomerism in 1-butyl-3-methyl-imidazolium tetrafluoroborate ([bmim][BF4]),” Journal of Raman Spectroscopy, vol. 39, no. 7, pp. 793–805, 2008.
- M. J. Frish, G. W. Trucks, H. B. Schlegel, et al., GAUSSIAN 03, Gaussian, Pittsburgh, Pa, USA, 2003.
- A. D. Becke, “Density-functional exchange-energy approximation with correct asymptotic behavior,” Physical Review A, vol. 38, no. 6, pp. 3098–3100, 1988.
- C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B, vol. 37, no. 2, pp. 785–789, 1988.
- S. F. Boys and F. Bernardi, “The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors,” Molecular Physics, vol. 19, no. 4, pp. 553–566, 1970.
- S. Simon, M. Duran, and J. J. Dannenberg, “How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?” Journal of Chemical Physics, vol. 105, no. 24, pp. 11024–11031, 1996.
- C. M. Breneman and K. B. Wiberg, “Determining atom-centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis,” Journal of Computational Chemistry, vol. 11, no. 3, pp. 361–373, 1990.
- K. B. Wiberg, T. A. Keith, M. J. Frisch, and M. Murcko, “Solvent effects on 1,2-dihaloethane gauche/trans ratios,” Journal of Physical Chemistry, vol. 99, no. 22, pp. 9072–9079, 1995.
- H. Yoshida and H. Matsuura, “Density functional study of the conformations and vibrations of 1,2-dimethoxyethane,” Journal of Physical Chemistry A, vol. 102, no. 16, pp. 2691–2699, 1998.
- N. Goutev, K. Ohno, and H. Matsuura, “Raman spectroscopic study on the conformation of 1,2-dimethoxyethane in the liquid phase and in aqueous solutions,” Journal of Physical Chemistry A, vol. 104, no. 40, pp. 9226–9232, 2000.
- V. B. Luzhkov, F. Österberg, P. Acharya, J. Chattopadhyaya, and J. Åqvist, “Computational and NMR study of quaternary ammonium ion conformations in solution,” Physical Chemistry Chemical Physics, vol. 4, no. 19, pp. 4640–4647, 2002.
- T. Takekiyo and Y. Yoshimura, “Drastic change in the conformational equilibria of tetraethylammonium bromide in the glassy aqueous solution,” Chemical Physics Letters, vol. 420, no. 1–3, pp. 8–11, 2006.
- T. Yasuda, M. Okuyama, N. Tanimoto, S. Sekiguchi, and S. I. Ikawa, “Intramolecular hydrogen bonding and molecular conformation. Part 4.-Corrections of thermodynamic properties obtained by IR measurements at different temperatures and pressures,” Journal of the Chemical Society, Faraday Transactions, vol. 91, no. 19, pp. 3379–3383, 1995.
- T. Welton, “Room-temperature ionic liquids. Solvents for synthesis and catalysis,” Chemical Reviews, vol. 99, no. 8, pp. 2071–2083, 1999.
- R. Katoh, M. Hara, and S. Tsuzuki, “Ion pair formation in [bmim]I ionic liquids,” Journal of Physical Chemistry B, vol. 112, no. 48, pp. 15426–15430, 2008.
- S. G. Raju and S. Balasubramanian, “Aqueous solution of [bmim][PF6]: ion and solvent effects on structure and dynamics,” Journal of Physical Chemistry B, vol. 113, no. 14, pp. 4799–4806, 2009.