Research Article  Open Access
Takahiro Takekiyo, Yusuke Imai, Hiroshi Abe, Yukihiro Yoshimura, "Conformational Analysis of Quaternary AmmoniumType Ionic Liquid Cation, N,NDiethylNmethylN(2methoxyethyl) Ammonium Cation", Advances in Physical Chemistry, vol. 2012, Article ID 829523, 7 pages, 2012. https://doi.org/10.1155/2012/829523
Conformational Analysis of Quaternary AmmoniumType Ionic Liquid Cation, N,NDiethylNmethylN(2methoxyethyl) Ammonium Cation
Abstract
Conformational preference of N,NdiethylNmethylN(2methoxyethyl) ammonium cation ([DEME]^{+}), which is one of the quaternary ammoniumbased 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 [1–3].
Recently, the phase behaviors and solution structures of N,NdiethylNmethylN(2methoxyethyl) ammonium tetrafluoroborate ([DEME][BF_{4}]) and iodide ([DEME][I]), which are typical quaternary ammoniumbased 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 Xray diffraction methods [4–12]. Imai et al. [4] and Abe et al. [8] reported the phase diagram of [DEME][BF_{4}]water mixed solutions as functions of temperature and water concentration. They found that [DEME][BF_{4}]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][BF_{4}]water mixed solutions at 298 and 77 K using a Raman spectroscopy. In this relation, Imai et al. [9, 12] reported that [DEME][] ( = BF_{4} and I)water mixed solutions have double glass transition temperatures () due to the liquidliquid 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][](= BF_{4} and I) solutions.
Although the conformations of 1alkyl3methylimidazoliumbased 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 1ethyl3methylimidazolium cation ([emim]^{+}) and 1butyl3methylimidazolium 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 threeparameter (B3) exchange function [20]. The B3 exchange function was combined with LeeYangPeer 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.
(a)
(b)
3. Results and Discussion
We consider the conformations on (i) the CN^{+}CCOC chain and (ii) the ethylchain (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/6311++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.
(a)
(b)
(c)
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,2dimethoxyethane (DME: CH_{3}OCH_{2}CH_{2}OCH_{3}) 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/6311++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 ethylchain for the conformational stability of the CN^{+}CCO^{−}C chain since the steric effect owing to the conformational preference of ethylchain group would change the conformational stability of the CN^{+}CCO^{−}C chain. In an effort to ensure the affection from the ethylchain for the conformational stability of the CN^{+}CCO^{−}C chain, we calculated the PES of CCN^{+}C dihedral angles of the ethylchain. Seven conformers of ethylchain 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 (Et_{4}N^{+}) 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 ethylchain conformers of [DEME]^{+} qualitatively agrees with the quantum chemistry calculation results of Et_{4}N^{+} [28, 29].

We checked the conformational preference of the CN^{+}CCO^{−}C chain again, when the ethylchain 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 ethylchain conformers. On the other hand, the change in the PES of CN^{+}CC angle is largely dependent of the ethylchain conformers since the position of the CN^{+}CC angle is very close with the ethylchain. Thus, transformation from the conformer to conformer of the ethylchain 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 ethylchain 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 ethylchain 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.
(a)
(b)
(c)
In view of the results from the conformational combinations of the CN^{+}CCO^{−}C chain and ethylchain, [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 ethylchain () 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 ethylchain 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]^{+} [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^{+}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][] ( = BF_{4} 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][]( = BF_{4} 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 hydrogenbonded 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,NdiethylNmethylN(2methoxyethyl) ammoniumbased 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.
References
 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. View at: Google Scholar
 T. Endo, T. Kato, K. I. Tozaki, and K. Nishikawa, “Phase behaviors of room temperature ionic liquid linked with cation conformational changes: 1Butyl3methylimidazolium Hexafluorophosphate,” Journal of Physical Chemistry B, vol. 114, no. 1, pp. 407–411, 2010. View at: Publisher Site  Google Scholar
 R. W. Berg, “Raman spectroscopy and abinitio model calculations on ionic liquids,” Monatshefte fur Chemie, vol. 138, no. 11, pp. 1045–1075, 2007. View at: Publisher Site  Google Scholar
 Y. Imai, H. Abe, T. Goto, Y. Yoshimura, Y. Michishita, and H. Matsumoto, “Structure and thermal property of N,NdiethylNmethylN2methoxyethyl ammonium tetrafluoroborateH_{2}O mixtures,” Chemical Physics, vol. 352, no. 1–3, pp. 224–230, 2008. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 Y. Imai, H. Abe, and Y. Yoshimura, “Xray diffraction study of ionic liquid based mixtures,” Journal of Physical Chemistry B, vol. 113, no. 7, pp. 2013–2018, 2009. View at: Publisher Site  Google Scholar
 Y. Yoshimura, T. Goto, H. Abe, and Y. Imai, “Existence of nearlyfree hydrogen bonds in an ionic liquid, N,N diethylNmethylN(2methoxyethyl) ammonium tetrafluoroboratewater at 77 K,” Journal of Physical Chemistry B, vol. 113, no. 23, pp. 8091–8095, 2009. View at: Publisher Site  Google Scholar
 H. Abe, Y. Yoshimura, Y. Imai, T. Goto, and H. Matsumoto, “Phase behavior of room temperature ionic liquid—H_{2}O mixtures: N, NdiethylNmethylN2methoxyethyl ammonium tetrafluoroborate,” Journal of Molecular Liquids, vol. 150, no. 1–3, pp. 16–21, 2009. View at: Publisher Site  Google Scholar
 Y. Imai, H. Abe, T. Miyashita, T. Goto, H. Matsumoto, and Y. Yoshimura, “Two glass transitions in N, NdiethylNmethylN(2methoxyethyl) ammonium tetrafluoroborateH_{2}O mixed solutions,” Chemical Physics Letters, vol. 486, no. 13, pp. 37–39, 2010. View at: Publisher Site  Google Scholar
 H. Abe, Y. Imai, T. Takekiyo, and Y. Yoshimura, “Deuterated water effect in a room temperature ionic liquid: N,NdiethylNmethyl2methoxyethyl ammonium tetrafluoroborate,” Journal of Physical Chemistry B, vol. 114, no. 8, pp. 2834–2839, 2010. View at: Publisher Site  Google Scholar
 H. Abe, Y. Imai, T. Goto et al., “Water network in roomtemperature ionic liquid: N,NdiethylNmethylN2methoxyethyl ammonium tetrafluoroborate,” Metallurgical and Materials Transactions A, vol. 41, no. 5, pp. 1137–1143, 2010. View at: Publisher Site  Google Scholar
 Y. Imai, H. Abe, H. Matsumoto, O. Shimada, T. Hanasaki, and Y. Yoshimura, “Glass transition behaviour of the quaternary ammonium type ionic liquid, $\{[\text{DEME}][\text{I}]+{\text{H}}_{2}\text{O}}$ mixtures,” Journal of Chemical Thermodynamics, vol. 43, no. 3, pp. 319–322, 2011. View at: Publisher Site  Google Scholar
 R. Ozawa, S. Hayashi, S. Saha, A. Kobayashi, and H. O. Hamaguchi, “Rotational isomerism and structure of the 1butyl3methylimidazolium cation in the ionic liquid state,” Chemistry Letters, vol. 32, no. 10, pp. 948–949, 2003. View at: Publisher Site  Google Scholar
 K. Fujii, T. Fujimori, T. Takamuku, R. Kanzaki, Y. Umebayashi, and S. I. Ishiguro, “Conformational equilibrium of bis(trifluoromethanesulfonyl) imide anion of a roomtemperature ionic liquid: raman spectroscopic study and DFT calculations,” Journal of Physical Chemistry B, vol. 110, no. 16, pp. 8179–8183, 2006. View at: Publisher Site  Google Scholar
 Y. Umebayashi, T. Fujimori, T. Sukizaki et al., “Evidence of conformational equilibrium of 1ethyl3methylimidazolium 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. View at: Publisher Site  Google Scholar
 J. C. Lassègues, J. Grondin, R. Holomb, and P. Johansson, “Raman and ab initio study of the conformational isomerism in the 1ethyl3methylimidazolium bis(trifluoromethanesulfonyl)imide ionic liquid,” Journal of Raman Spectroscopy, vol. 38, no. 5, pp. 551–558, 2007. View at: Publisher Site  Google Scholar
 Y. Umebayashi, T. Mitsugi, S. Fukuda et al., “Lithium ion solvation in roomtemperature 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. View at: Publisher Site  Google Scholar
 R. Holomb, A. Martinelli, I. Albinsson, J. C. Lassègues, P. Johansson, and P. Jacobsson, “Ionic liquid structure: the conformational isomerism in 1butyl3methylimidazolium tetrafluoroborate ([bmim][BF4]),” Journal of Raman Spectroscopy, vol. 39, no. 7, pp. 793–805, 2008. View at: Publisher Site  Google Scholar
 M. J. Frish, G. W. Trucks, H. B. Schlegel et al., GAUSSIAN 03, Gaussian, Pittsburgh, Pa, USA, 2003.
 A. D. Becke, “Densityfunctional exchangeenergy approximation with correct asymptotic behavior,” Physical Review A, vol. 38, no. 6, pp. 3098–3100, 1988. View at: Publisher Site  Google Scholar
 C. Lee, W. Yang, and R. G. Parr, “Development of the ColleSalvetti correlationenergy formula into a functional of the electron density,” Physical Review B, vol. 37, no. 2, pp. 785–789, 1988. View at: Publisher Site  Google Scholar
 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. View at: Google Scholar
 S. Simon, M. Duran, and J. J. Dannenberg, “How does basis set superposition error change the potential surfaces for hydrogenbonded dimers?” Journal of Chemical Physics, vol. 105, no. 24, pp. 11024–11031, 1996. View at: Google Scholar
 C. M. Breneman and K. B. Wiberg, “Determining atomcentered 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. View at: Google Scholar
 K. B. Wiberg, T. A. Keith, M. J. Frisch, and M. Murcko, “Solvent effects on 1,2dihaloethane gauche/trans ratios,” Journal of Physical Chemistry, vol. 99, no. 22, pp. 9072–9079, 1995. View at: Google Scholar
 H. Yoshida and H. Matsuura, “Density functional study of the conformations and vibrations of 1,2dimethoxyethane,” Journal of Physical Chemistry A, vol. 102, no. 16, pp. 2691–2699, 1998. View at: Google Scholar
 N. Goutev, K. Ohno, and H. Matsuura, “Raman spectroscopic study on the conformation of 1,2dimethoxyethane in the liquid phase and in aqueous solutions,” Journal of Physical Chemistry A, vol. 104, no. 40, pp. 9226–9232, 2000. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 T. Welton, “Roomtemperature ionic liquids. Solvents for synthesis and catalysis,” Chemical Reviews, vol. 99, no. 8, pp. 2071–2083, 1999. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 S. G. Raju and S. Balasubramanian, “Aqueous solution of [bmim][PF_{6}]: ion and solvent effects on structure and dynamics,” Journal of Physical Chemistry B, vol. 113, no. 14, pp. 4799–4806, 2009. View at: Publisher Site  Google Scholar
Copyright
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.