Abstract

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.

1. Introduction

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 [13].

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 [412]. Imai et al. [4] and Abe et al. [8] 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. [7] 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, 1318], 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.

2. Experimental

Density functional theory (DFT) calculations were carried out using GAUSSIAN03 program [19]. For the present calculation of [DEME]+, we used the Becke’s three-parameter (B3) exchange function [20]. The B3 exchange function was combined with Lee-Yang-Peer correlation function (B3LYP) [21]. 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+CCOC chain. The changes in potential energy surfaces (PES) for the CN+CC, N+CCO, and CCOC 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 CCOC 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+CCOC chain. Therefore, the CN+CCOC 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+CCOC chain since the steric effect owing to the conformational preference of ethyl-chain group would change the conformational stability of the CN+CCOC chain. In an effort to ensure the affection from the ethyl-chain for the conformational stability of the CN+CCOC 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+CCOC 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 CCOC 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+CCOC 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+CCOC chain takes the six conformers; the , , , , , and conformers, respectively.

In view of the results from the conformational combinations of the CN+CCOC 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+CCOC 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+CCOC chain. As mentioned above, [DEME]+ takes many conformers for the N+CCOC 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]+ [15] and between the and conformers of [bmim]+ [18] 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+CCOC chain. Our calculation results indicated that the electron polarizabities of the N and O atoms in the CN+CCOC 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 [412].

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 [30]. 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 [31]. 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+CCOC 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.