Article pubs.acs.org/JPCB
Structure of 1‑Alkyl-1-methylpyrrolidinium Bis(trifluoromethylsulfonyl)amide Ionic Liquids with Linear, Branched, and Cyclic Alkyl Groups Hemant K. Kashyap,† Cherry S. Santos,‡ N. Sanjeeva Murthy,§ Jeevapani J. Hettige,† Kijana Kerr,∥ Sharon Ramati,⊥ JinHee Gwon,⊥ Masao Gohdo,∥ Sharon I. Lall-Ramnarine,⊥ James F. Wishart,∥ Claudio J. Margulis,*,† and Edward W. Castner, Jr.*,‡ †
Department of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States § New Jersey Center for Biomaterials, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States ∥ Chemistry Department, Brookhaven National Laboratory, Upton, New York 11973-5000, United States ⊥ Department of Chemistry, Queensborough Community College-CUNY, Bayside, New York 11364, United States ‡
S Supporting Information *
ABSTRACT: X-ray scattering and molecular dynamics simulations have been carried out to investigate structural differences and similarities in the condensed phase between pyrrolidinium-based ionic liquids paired with the bis(trifluoromethylsulfonyl)amide (NTf2−) anion where the cationic tail is linear, branched, or cyclic. This is important in light of the charge and polarity type alternations that have recently been shown to be present in the case of liquids with cations of moderately long linear tails. For this study, we have chosen to use the 1-alkyl-1-methylpyrrolidinium, Pyrr1,n+ with n = 5 or 7, as systems with linear tails, 1-(2ethylhexyl)-1-methylpyrrolidinium, Pyrr1,EtHx+, as a system with a branched tail, and 1-(cyclohexylmethyl)-1-methylpyrrolidinium, Pyrr1,ChxMe+, as a system with a cyclic tail. We put these results into context by comparing these data with recently published results for the Pyrr1,n+/NTf2− ionic liquids with n = 4, 6, 8, and 10.1,2 General methods for interpreting the structure function S(q) in terms of q-dependent natural partitionings are described. This allows for an in-depth analysis of the scattering data based on molecular dynamics (MD) trajectories that highlight the effect of modifying the cationic tail.
■
INTRODUCTION Modern room-temperature ionic liquids (RTILs) are novel materials that are being pursued for a variety of applications including energy storage in solar cells, supercapacitor devices, lithium ion batteries, as well as catalysis and separation processes.3−8 In the past decade, a significant number of Xray1,9−34 and neutron35−37 scattering experiments as well as MD simulations23,28,31,32,38−51 have been carried out to unravel the local and mesoscale structure of RTILs. As our group has indicated in several recent publications,2,23,28,50 the X-ray structure function S(q) of RTILs typically has either two or three peaks in the relevant intermolecular region at q values Pyrr1,8+ > Pyrr1,10+. This ordering will become more clear as we discuss the different partitionings of S(q) in subsequent sections. Since we plan to interpret the experimental S(q) by using computer simulations, we first show in Figure 3 the excellent agreement between computed and experimental results. The high level of agreement between simulations and experiments permits us to use the computational structure function and different C
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX
The Journal of Physical Chemistry B
Article
partitioning of S(q) sums the contributions from the cationic and anionic charged parts together in the polar component. This results in destructive interference in the charge alternation region of S(q). In essence, the “natural partition” to study polarity alternation is that for which polar and apolar groups are separately identified as different contributors to the overall S(q). Conversely, the natural partitioning to study charge alternation is one for which cations and anions are not grouped together but instead constitute mutually exclusive S(q) subcomponents. Figure 5 shows that for the case of Pyrr1,10+ a large peak near q = 0.35 Å−1 is observed both in the polar−polar and apolar− apolar subcomponents of S(q). An antipeak at the same value of scattering vector is observed for the polar−apolar subcomponent. This is the signature of polarity alternation. The signature of alternation becomes less clear when the cationic alkyl groups are shorter, branched, or cyclic, since both peaks and antipeaks become broad and poorly defined for this partitioning of S(q). Ionic Partitioning of S(q). Figure 6 depicts a Pyrr1,7+/ NTf2− simulation box color coded according to the ionic
Figure 4. (top) A snapshot of the Pyrr1,7+/NTf2− system colored based on the polarity partitioning of S(q). (bottom) A single ion pair with the same color coding. Green corresponds to polar and white to apolar.
should be positive going, whereas the cross term (polar− apolar) + (apolar−polar) contribution to S(q) should be negative going.2,54 Peaks and antipeaks in the structure function S(q) are a signature of alternation because they correspond to features of the same spatial periodicity but with an origin offset.2 Figure 5
Figure 6. (top) A snapshot of the Pyrr1,7+/NTf2− system colored on the basis of the ionic partitioning of S(q). (bottom) A single ion pair with the same color coding. Blue corresponds to cation and green to anion.
partitioning scheme. The most obvious use of the ionic partitioning of S(q) is to identify the region in reciprocal space in which charge alternation occurs. We expect all ionic liquids to show some type of charge alternation in the ionic partitioning of S(q) independent of their nonpolar tails even when this is absent because of cancellations in the overall S(q). Figure 7 shows exactly this for q ≈ 0.85 Å−1. The cation−cation and anion−anion contributions in Figure 7 lead to peaks for all systems studied, whereas the (cation−anion) + (anion−cation) contributions lead to antipeaks. This is the ubiquitous signature of charge alternation in ionic liquids.
Figure 5. Polar−apolar partitioning of the structure function S(q). Each of the partial structure functions is spaced by +2 units vertically.
shows a plot of the polarity partitioning of S(q) highlighting this effect.54 When analyzing Figure 5, the reader should only focus on the prepeak region (i.e., q < 0.5 Å−1) for which the polar/apolar partitioning of S(q) holds special significance. On the contrary, the polar/apolar partitioning is not useful to interpret features such as charge alternation that occur at around q = 0.85 Å−1 in S(q). This is because the polar/apolar D
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX
The Journal of Physical Chemistry B
Article
In the same way, Pyrr1,ChxMe+ appears similar to Pyrr1,5+ and not to cations with longer linear tails. Sub-Ionic Partitioning of S(q). In the previous subsection, we indicated that an ionic partitioning of S(q) was advantageous in order to better understand the charge alternation region of S(q). However, with the exception of the anionic component, the ionic partitioning is of little value to study polarity alternation in the prepeak region of S(q). The appreciation of the importance of anionic contribution to the prepeak region and the apparent lack of importance of cationic contribution to the same region has been reported before by our group.50 The underlying reason for this correct but perhaps misleading observation is that the cationic contribution to the prepeak in the ionic partitioning of S(q) is convoluted. This is because both polar and apolar components of S(q) project onto the correlations of cations with cations and those of cations with anions, leading to net cancellations in S(q). In order to deconvolute this information, one must resort to a subionic partitioning of S(q). Figure 9 depicts a Pyrr1,7+/NTf2− simulation box color coded according to the subionic partitioning scheme. As one would
Figure 7. Cation−cation (left), cation−anion (middle), and anion− anion (right) partial structure functions for all the ionic liquids studied here. Curves are offset vertically by 2 units for visual clarity.
Note that, while we cannot extract much useful information from this partitioning in the case of cation−cation or (cation− anion) + (anion−cation) interactions in the prepeak region, much can be learned instead by looking at the anion−anion subcomponent of S(q). This is because both polar and apolar components of S(q) are included in the cationic and (cation− anion) + (anion−cation) partitions but only purely polar components of S(q) are projected onto the anionic partition. Therefore, in this particular partitioning of S(q), information related to polarity alternation that occurs at the prepeak region is not convoluted only for the case of the purely anionic subcomponent of S(q). In the ionic partitioning of S(q), if one is interested in the prepeak, one should only analyze Figure 7 (right). At the value of q for the prepeak, the left and center graphs show data that is convoluted. Figure 7 (right) shows as expected a positive going peak due to polarity alternation. This peak is better defined and more prominent when tails are longer. The purely anion−anion component of the prepeak in Figure 8a and b reveals the reciprocal space value at which polarity alternation occurs. These figures show that at this q value Pyrr1,EtHx+ is much more comparable to Pyrr1,6+ (similar in length) and not to Pyrr1,8+ (same number of tail carbon atoms). Figure 9. (top) A snapshot of the Pyrr1,7+/NTf2− system colored on the basis of the subionic partitioning of S(q). (bottom) A single ion pair with the same color coding. Blue corresponds to cation head, green to anion, and white to cation tail. Notice that head is defined as including the ring, the methyl group, and up to the second methylene group of the longer alkyl tail.
predict, for the case of the longer cations, Figure 10 unambiguously demonstrates that cation−cation polar subcomponents (i.e., head−head) show positive going peaks at the prepeak region. At the same time, polar−apolar antipeaks can be seen in the same region when considering head−tail correlations. Furthermore, the apolar−apolar correlation of tails also gives rise to positive going peaks. These findings are again the signature of polarity alternation and are in complete agreement with matching observations in Figure 5. With increasing length of the cationic alkyl groups, cancellations in the prepeak region due to convolution of
Figure 8. The anion−anion subcomponents of total S(q) for (a) Pyrr1,ChxMe+ (black), Pyrr1,5+ (red), and Pyrr1,7+ (green) and (b) Pyrr1,EtHx+ (black), Pyrr1,6+ (red), and Pyrr1,8+ (green) RTILs. E
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX
The Journal of Physical Chemistry B
Article
Figure 10. Cation−cation partial structure function split into head− head (left), head−tail (middle), and tail−tail (right) subcomponents. Curves are vertically offset by 2 units for visual clarity.
Figure 11. Cation−anion contributions to the structure function S(q) partitioned into cation head−anion and cation tail−anion subcomponents. Curves are offset vertically by 2 units for visual clarity.
polar and apolar components in the ionic partitioning of S(q) are partially overcome. However, it is much better to examine the polarity partitioning or the subionic partitioning of S(q) to properly understand longer-range polarity alternations in the prepeak region. Because of its striking simplicity, perhaps the most elegant way73 to concurrently check for the existence of charge and polarity alternations in ionic liquids is to plot the cationic head−anion subionic partitioning of S(q). Since both the cationic head and the anion are polar, their correlation should always give rise to a positive going peak in the prepeak region if polarity alternation is indeed present in the system. Furthermore, because cation heads and anions carry charges of opposite sign, their contribution to S(q) should always appear as an antipeak in the charge alternation region. This function, as opposed to others described in previous sections, does not suffer from convolution problems either at the prepeak region or at the charge alternation region in reciprocal space. We expect these general principles to hold true for most if not all ionic liquids. A caveat is in place here; our line of reasoning in this and previous articles is based on the well established model of cations with nonpolar tails and anions that are completely polar. The situation will become more complex for ILs comprising cations and anions that both have longer alkyl tail substituents such as 1-octyl-1-methylpyrrolidinium octylsulfate or for systems in which tails are polar.55 Figure 11 (left) demonstrates that the cationic head−anion subionic partitioning of S(q) carries all the desirable information needed to establish polarity alternation and charge alternation regions. Systems with longer cationic alkyl groups show both charge and polarity alternation, whereas those with shorter tails only show charge alternation. Figure 11 (right) only teaches us about polarity alternation but not about charge alternation. Systems with longer cationic alkyl groups show antipeaks in the cation-tail−anion subionic partitioning of S(q). Figures 12 and 13 provide a clear comparison of the subionic contributions of cyclic and branched systems against those with linear tails. In all cases, the left panel shows that for branched and cyclic tailed cations polar−polar components scatter in a way that is similar to systems with linear tails of approximately the same spatial length but not of the same number of carbons.
Figure 12. Cation−cation partial structure function split into head− head, head−tail, and tail−tail terms for (a) Pyrr1,ChxMe+ (black), Pyrr1,5+ (red), and Pyrr1,7+ (green) and (b) Pyrr1,EtHx+ (black), Pyrr1,6+ (red), and Pyrr1,8+ (green) RTILs.
This analogy is not as clear for the other subionic components of S(q). Bulk Densities. Figure 14 shows that, in the case of the systems in which cations have linear alkyl tail groups, experimental and simulated densities (see Table S1, Supporting Information) follow very similar trends with density decreasing F
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX
The Journal of Physical Chemistry B
Article
with increasing number of tail carbons. For these liquids, computational densities are about 5% higher than the corresponding experimental densities. Density deviations between simulation and experiment in the case of the branched and cyclic cationic systems are less than 3.5%. In our experience from previous studies and from the data presented in this article, excellent agreement between experimental and computational S(q) can be obtained even if small density discrepancies are observed. S(q) does not appear to be sensitive to these small density differences. For fixed n = 7, the density of the system with the cyclic chain is significantly larger by about 5 or 6% with respect to that of the linear chain cation system. Instead, the density difference between the system with linear cationic tail of n = 8 and the branched system is much less significant. Notably, an earlier study75 found that the branched sec-butyl isomer of Pyrr1,4+/NTf2− was denser than the linear isomer (1.46 vs 1.39 g/cm3), in keeping with the trend observed here.
■
CONCLUSIONS In this article, we compare the structure of linear, branched, and cyclic alkyl tail pyrrolidinium bis(trifluoromethylsulfonyl)amide RTILs. From a methodological perspective, we found that the ideal way to study these and other ionic liquids is by using qdependent partitionings of the structure function. These qdependent partitionings enable us to clearly identify the physical origin of the different intermolecular features of S(q). As expected, ILs having cations with long linear alkyl groups show both the signs of polarity alternation and charge alternation. ILs comprising cations with shorter alkyl tails show charge alternation but no polarity alternation. The behavior of our branched and cyclic systems at the prepeak region is from the perspective of the charged groups similar to that observed for systems where the alkyl tail in the cations is of similar length but not necessarily of similar number of carbon atoms. The behavior of the polar−apolar components and apolar−apolar components in these two systems is somewhat more complex and not necessarily similar to that of systems with the same tail length.
Figure 13. Cation−anion partial structure function split into cation head−anion and cation tail−anion terms for (a) Pyrr1,ChxMe+ (black), Pyrr1,5+ (red), and Pyrr1,7+ (green) and (b) Pyrr1,EtHx+ (black), Pyrr1,6+ (red), and Pyrr1,8+ (green) RTILs.
■
ASSOCIATED CONTENT
S Supporting Information *
Details on the synthesis of the compounds and a table of measured and simulated densities. This material is available free of charge via the Internet at http://pubs.acs.org.
■
AUTHOR INFORMATION
Corresponding Author
*E-mail:
[email protected] (C.J.M.); ed.castner@ rutgers.edu (E.W.C.). Author Contributions
H.K.K. and C.J.M. completed the molecular simulations and theoretical analysis of the structure factor. J.J.H. worked with H.K.K. in checking the correctness of all force field parameters in various simulation input files. C.S.S., N.S.M., and E.W.C. carried out the X-ray experiments at APS and analyzed the Xray data. K.K., S.R., J.G., M.G., S.L.-R., and J.F.W. prepared and characterized three ionic liquids. H.K.K. drafted most of the manuscript in collaboration with C.J.M. and E.W.C.
Figure 14. Experimental and simulated bulk densities of the RTILs as a function of the number of alkyl chain carbon atoms, n. The densities of Pyrr1,n+/NTf2− RTILs are shown with black circles and those of Pyrr1,ChxMe+/NTf2− and Pyrr1,EtHx+/NTf2− with red circles. All simulated densities are at 295 K. Reported experimental densities74−76 were measured between 293 and 303 K, as indicated in Table S1 of the Supporting Information.
Notes
The authors declare no competing financial interest. G
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX
The Journal of Physical Chemistry B
■
Article
(13) Atkin, R.; Warr, G. G. The Smallest Amphiphiles: Nanostructure in Protic Room-Temperature Ionic Liquids with Short Alkyl Groups. J. Phys. Chem. B 2008, 112, 4164−4166. (14) Fukuda, S.; Takeuchi, M.; Fujii, K.; Kanzaki, R.; Takamuku, T.; Chiba, K.; Yamamoto, H.; Umebayashi, Y.; Ishiguro, S.-i. Liquid Structure of N-Butyl-N-methylpyrrolidinium Bis(trifluoromethanesulfonyl) Amide Ionic Liquid Studied by Large Angle X-ray Scattering and Molecular Dynamics Simulations. J. Mol. Liq. 2008, 143, 2−7. (15) Fujii, K.; Seki, S.; Fukuda, S.; Takamuku, T.; Kohara, S.; Kameda, Y.; Umebayashi, Y.; Ishiguro, S.-i. Liquid Structure and Conformation of a Low-Viscosity Ionic Liquid, N-Methyl-N-propylpyrrolidinium Bis(fluorosulfonyl) Imide Studied by High-Energy X-ray Scattering. J. Mol. Liq. 2008, 143, 64−69. (16) Xiao, D.; Hines, L. G., Jr.; Li, S.; Bartsch, R. A.; Quitevis, E. L.; Russina, O.; Triolo, A. Effect of Cation Symmetry and Alkyl Chain Length on the Structure and Intermolecular Dynamics of 1,3Dialkylimidazolium Bis(trifluoromethanesulfonyl)amide Ionic Liquids. J. Phys. Chem. B 2009, 113, 6426−6433. (17) Triolo, A.; Russina, O.; Fazio, B.; Appetecchi, G. B.; Carewska, M.; Passerini, S. Nanoscale Organization in Piperidinium-Based Room Temperature Ionic Liquids. J. Chem. Phys. 2009, 130, 164521. (18) Russina, O.; Triolo, A.; Gontrani, L.; Caminiti, R.; Xiao, D.; Hines, L. G., Jr.; Bartsch, R. A.; Quitevis, E. L.; Pleckhova, N.; Seddon, K. R. Morphology and Intermolecular Dynamics of 1-Alkyl-3methylimidazolium Bis(trifluoromethane)sulfonylamide Ionic Liquids: Structural and Dynamic Evidence of Nanoscale Segregation. J. Phys.: Condens. Matter 2009, 21, 424121. (19) Gontrani, L.; Russina, O.; Celso, F. L.; Caminiti, R.; Annat, G.; Triolo, A. Liquid Structure of Trihexyltetradecylphosphonium Chloride at Ambient Temperature: an X-ray Scattering and Simulation Study. J. Phys. Chem. B 2009, 113, 9235−9240. (20) Pott, T.; Méléard, P. New Insight into the Nanostructure of Ionic Liquids: a Small Angle X-ray Scattering (SAXS) Study on Liquid Tri-alkyl-methyl-ammonium Bis(trifluoromethanesulfonyl)amides and Their Mixtures. Phys. Chem. Chem. Phys. 2009, 11, 5469−5475. (21) Fujii, K.; Mitsugi, T.; Takamuku, T.; Yanaguchi, T.; Umebayashi, Y.; Ishiguro, S. Effect of Methylation at the C2 Position of Imidazolium on the Structure of Ionic Liquids Revealed by Large Angle X-ray Scattering Experiments and MD Simulations. Chem. Lett. 2009, 38, 340−341. (22) Zheng, W.; Mohammed, A.; Hines, L. G.; Xiao, D.; Martinez, O. J.; Bartsch, R. A.; Simon, S. L.; Russina, O.; Triolo, A.; Quitevis, E. L. Effect of Cation Symmetry on the Morphology and Physicochemical Properties of Imidazolium Ionic Liquids. J. Phys. Chem. B 2011, 115, 6572−6584. (23) Santos, C. S.; Annapureddy, H. V. R.; Murthy, N. S.; Kashyap, H. K.; Castner, E. W., Jr.; Margulis, C. J. Temperature-Dependent Structure of Methyltributylammonium Bis(trifluoromethylsulfonyl)amide: X-ray Scattering and Simulations. J. Chem. Phys. 2011, 134, 064501. (24) Aoun, B.; Goldbach, A.; Gonzalez, M. A.; Kohara, S.; Price, D. L.; Saboungi, M.-L. Nanoscale Heterogeneity in Alkyl-Methylimidazolium Bromide Ionic Liquids. J. Chem. Phys. 2011, 134, 104509. (25) Fujii, K.; Kanzaki, R.; Takamuku, T.; Kameda, Y.; Kohara, S.; Kanakubo, M.; Shibayama, M.; ichi Ishiguro, S.; Umebayashi, Y. Experimental Evidences for Molecular Origin of Low-Q Peak in Neutron/X-ray Scattering of 1-Alkyl-3-methylimidazolium Bis(trifluoromethanesulfonyl)amide Ionic Liquids. J. Chem. Phys. 2011, 135, 244502. (26) Hayes, R.; Imberti, S.; Warr, G. G.; Atkin, R. Pronounced Sponge-Like Nanostructure in Propylammonium Nitrate. Phys. Chem. Chem. Phys. 2011, 13, 13544−13551. (27) Umebayashi, Y.; Hamano, H.; Seki, S.; Minofar, B.; Fujii, K.; Hayamizu, K.; Tsuzuki, S.; Kameda, Y.; Kohara, S.; Watanabe, M. Liquid Structure of and Li+ Ion Solvation in Bis(trifluoromethanesulfonyl)amide Based Ionic Liquids Composed of 1-Ethyl-3-methylimidazolium and N-Methyl-N-propylpyrrolidinium Cations. J. Phys. Chem. B 2011, 115, 12179−12191.
ACKNOWLEDGMENTS The authors thank Prof. Mark Maroncelli for the sample of Pyrr1,5+/NTf2−, Dr. Tomasz Szreder for measuring the density of Pyrr1,4+/NTf2−, Dr. Alison Funston for measuring the density of Pyrr1,5+/NTf2−, and Jasmine Hatcher for assistance preparing scattering samples. We thank Dr. Chris Benmore and Dr. Yang Ren for help in the X-ray measurements and data analysis at APS beamline 11-ID-C. This work was supported by the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences under contracts DE-SC0001780 (E.W.C.), DE-SC0008644 (C.J.M.), and DE-AC02-98CH10886 (J.F.W.). Preparation of Pyrr1,ChxMe+ and Pyrr1,7+ salts by K.K. was supported at BNL by the DOE Office of Nuclear Energy. Researchers (S.L.-R., S.R., J.G.) from Queensborough Community College were supported by PSC-CUNY Research Award Grants and the Louis Stokes Alliance for Minority Participation. Use of the Advanced Photon Source at Argonne National Laboratory was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC0206CH11357.
■
REFERENCES
(1) Santos, C. S.; Murthy, N. S.; Baker, G. A.; Castner, E. W., Jr. Communication: X-ray Scattering from Ionic Liquids with Pyrrolidinium Cations. J. Chem. Phys. 2011, 134, 121101. (2) Kashyap, H. K.; Hettige, J. J.; Annapureddy, H. V. R.; Margulis, C. J. SAXS Anti-Peaks Reveal the Length-Scales of Dual PositiveNegative and Polar-Apolar Ordering in Room-Temperature Ionic Liquids. Chem. Commun. 2012, 48, 5103−5105. (3) Wishart, J. F. Energy Applications of Ionic Liquids. Energy Environ. Sci. 2009, 2, 956−961. (4) Castner, E. W., Jr.; Wishart, J. F. Spotlight on Ionic Liquids. J. Chem. Phys. 2010, 132, 120901. (5) MacFarlane, D. R.; Huang, J.; Forsyth, M. Lithium-Doped Plastic Crystal Electrolytes Exhibiting Fast Ion Conduction for Secondary Batteries. Nature 1999, 402, 792−794. (6) Wang, P.; Wenger, B.; Humphry-Baker, R.; Moser, J.-E.; Teuscher, J.; Kantlehner, W.; Mezger, J.; Stoyanov, E. V.; Zakeeruddin, S. M.; Grätzel, M. Charge Separation and Efficient Light Energy Conversion in Sensitized Mesoscopic Solar Cells Based on Binary Ionic Liquids. J. Am. Chem. Soc. 2005, 127, 6850−6856. (7) Kuang, D.; Wang, P.; Ito, S.; Zakeeruddin, S. M.; Grätzel, M. Stable Mesoscopic Dye-Sensitized Solar Cells Based on Tetracyanoborate Ionic Liquid Electrolyte. J. Am. Chem. Soc. 2006, 128, 7732− 7733. (8) Han, X.; Armstrong, D. W. Ionic Liquids in Separations. Acc. Chem. Res. 2007, 40, 1079−1086. (9) Bradley, A. E.; Hardacre, C.; Holbrey, J. D.; Johnston, S.; McMath, S. E. J.; Nieuwenhuyzen, M. Small-Angle X-ray Scattering Studies of Liquid Crystalline 1-Alkyl-3-methylimidazolium Salts. Chem. Mater. 2002, 14, 629−635. (10) Triolo, A.; Mandanici, A.; Russina, O.; Rodriguez-Mora, V.; Cutroni, M.; Hardacre, C.; Nieuwenhuyzen, M.; Bleif, H.-J.; Keller, L.; Ramos, M. A. Thermodynamics, Structure, and Dynamics in Room Temperature Ionic Liquids: The Case of 1-Butyl-3-methyl Imidazolium Hexafluorophosphate ([bmim][PF6]). J. Phys. Chem. B 2006, 110, 21357−21364. (11) Triolo, A.; Russina, O.; Bleif, H.-J.; Di Cola, E. Nanoscale Segregation in Room Temperature Ionic Liquids? J. Phys. Chem. B 2007, 111, 4641−4644. (12) Triolo, A.; Russina, O.; Fazio, B.; Triolo, R.; DiCola, E. Morphology of 1-Alkyl-3-methylimidazolium Hexafluorophosphate Room Temperature Ionic Liquids. Chem. Phys. Lett. 2008, 457, 362−365. H
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX
The Journal of Physical Chemistry B
Article
(28) Kashyap, H. K.; Santos, C. S.; Annapureddy, H. V. R.; Murthy, N. S.; Margulis, C. J.; Castner, E. W., Jr. Temperature-Dependent Structure of Ionic Liquids: X-ray Scattering and Simulations. Faraday Discuss. 2012, 154, 133−143. (29) Russina, O.; Triolo, A. New Experimental Evidence Supporting the Mesoscopic Segregation Model in Room Temperature Ionic Liquids. Faraday Discuss. 2012, 154, 97−109. (30) Russina, O.; Triolo, A.; Gontrani, L.; Caminiti, R. Mesoscopic Structural Heterogeneities in Room-Temperature Ionic Liquids. J. Phys. Chem. Lett. 2012, 3, 27−33. (31) Macchiagodena, M.; Ramondo, F.; Triolo, A.; Gontrani, L.; Caminiti, R. Liquid Structure of 1-Ethyl-3-methylimidazolium Alkyl Sulfates by X-ray Scattering and Molecular Dynamics. J. Phys. Chem. B 2012, 116, 13448−13458. (32) Li, S.; Banuelos, J. L.; Guo, J.; Anovitz, L.; Rother, G.; Shaw, R. W.; Hillesheim, P. C.; Dai, S.; Baker, G. A.; Cummings, P. T. Alkyl Chain Length and Temperature Effects on Structural Properties of Pyrrolidinium-Based Ionic Liquids: A Combined Atomistic Simulation and Small-Angle X-ray Scattering Study. J. Phys. Chem. Lett. 2012, 3, 125−130. (33) Triolo, A.; Russina, O.; Caminiti, R.; Shirota, H.; Lee, H. Y.; Santos, C. S.; Murthy, N. S.; Castner, E. W., Jr. Comparing Intermediate Range Order for Alkyl- vs. Ether-Substituted Cations in Ionic Liquids. Chem. Commun. 2012, 48, 4959−4961. (34) Kofu, M.; Nagao, M.; Ueki, T.; Kitazawa, Y.; Nakamura, Y.; Sawamura, S.; Watanabe, M.; Yamamuro, O. Heterogeneous Slow Dynamics of Imidazolium-Based Ionic Liquids Studied by Neutron Spin Echo. J. Phys. Chem. B 2013, 117, 2773−2781. (35) Triolo, A.; Russina, O.; Arrighi, V.; Juranyi, F.; Janssen, S.; Gordon, C. Quasielastic Neutron Scattering Characterization of the Relaxation Processes in a Room Temperature Ionic Liquid. J. Chem. Phys. 2003, 119, 8549−8557. (36) Hardacre, C.; Holbrey, J. D.; Nieuwenhuyzen, M.; Youngs, T. G. A. Structure and Solvation in Ionic Liquids. Acc. Chem. Res. 2007, 40, 1146−1155. (37) Hardacre, C.; Holbrey, J. D.; Mullan, C. L.; Youngs, T. G. A.; Bowron, D. T. Small Angle Neutron Scattering from 1-Alkyl-3methylimidazolium Hexafluorophosphate Ionic Liquids ([Cnmim][PF6], n = 4, 6, and 8). J. Chem. Phys. 2010, 133, 74510−74517. (38) Urahata, S. M.; Ribeiro, M. C. Structure of Ionic Liquids of 1Alkyl-3-methylimidazolium Cations: A Systematic Computer Simulation Study. J. Chem. Phys. 2004, 120, 1855 −1863. (39) Del Popolo, M. G.; Voth, G. A. On the Structure and Dynamics of Ionic Liquids. J. Phys. Chem. B 2004, 108, 1744−1752. (40) Wang, Y.; Voth, G. A. Unique Spatial Heterogeneity in Ionic Liquids. J. Am. Chem. Soc. 2005, 127, 12192−12193. (41) Wang, Y.; Voth, G. A. Tail Aggregation and Domain Diffusion in Ionic Liquids. J. Phys. Chem. B 2006, 110, 18601−18608. (42) Borodin, O.; Smith, G. D. Structure and Dynamics of N-MethylN-propylpyrrolidinium Bis(trifluoromethanesulfonyl)imide Ionic Liquid from Molecular Dynamics Simulations. J. Phys. Chem. B 2006, 110, 11481−11490. (43) Bhargava, B. L.; Devane, R.; Klein, M. L.; Balasubramanian, S. Nanoscale Organization in Room Temperature Ionic Liquids: a Coarse Grained Molecular Dynamics Simulation Study. Soft Matter 2007, 3, 1395−1400. (44) Bhargava, B. L.; Balasubramanian, S.; Klein, M. L. Modelling Room Temperature Ionic Liquids. Chem. Commun. 2008, 29, 3339− 3351. (45) Bhargava, B.; Klein, M.; Balasubramanian, S. Structural Correlations and Charge Ordering in a Room-Temperature Ionic Liquid. ChemPhysChem 2008, 9, 67−70. (46) Lopes, J. N. C.; Pádua, A. A. H. Nanostructural Organization in Ionic Liquids. J. Phys. Chem. B 2006, 110, 3330−3335. (47) Canongia Lopes, J. N.; Shimizu, K.; Pádua, A. A. H.; Umebayashi, Y.; Fukuda, S.; Fujii, K.; Ishiguro, S.-i. A Tale of Two Ions: The Conformational Landscapes of Bis(trifluoromethanesulfonyl)amide and N,N-Dialkylpyrrolidinium. J. Phys. Chem. B 2008, 112, 1465−1472.
(48) Morrow, T. I.; Maginn, E. J. Molecular Dynamics Study of the Ionic Liquid 1-n-Butyl-3-methylimidazolium Hexafluorophosphate. J. Phys. Chem. B 2002, 106, 12807−12813. (49) Maginn, E. J. Molecular Simulation of Ionic Liquids: Current Status and Future Opportunities. J. Phys.: Condens. Matter 2009, 21, 373101. (50) Annapureddy, H. V. R.; Kashyap, H. K.; De Biase, P. M.; Margulis, C. J. What Is the Origin of the Prepeak in the X-ray Scattering of Imidazolium-Based Room-Temperature Ionic Liquids? J. Phys. Chem. B 2010, 114, 16838−16846. (51) Siqueira, L. J. A.; Ribeiro, M. C. C. Charge Ordering and Intermediate Range Order in Ammonium Ionic Liquids. J. Chem. Phys. 2011, 135, 204506. (52) Castner, E. W., Jr.; Margulis, C. J.; Maroncelli, M.; Wishart, J. F. Ionic Liquids: Structure and Photochemical Reactions. Annu. Rev. Phys. Chem. 2011, 62, 85−105. (53) Discussion, G. Faraday Discuss. 2012, 154, 189−220. (54) Kashyap, H. K.; Margulis, C. J. Theoretical Deconstruction of the X-ray Structure Function Exposes Polarity Alternations in Room Temperature Ionic Liquids. ECS Trans. 2013, 50, 301−307. (55) Kashyap, H. K.; Santos, C. S.; Daly, R. P.; Hettige, J. J.; Murthy, N. S.; Shirota, H.; Castner, E. W.; Margulis, C. J. How Does the Ionic Liquid Organizational Landscape Change When Nonpolar Cationic Alkyl Groups Are Replaced by Polar Isoelectronic Diethers? J. Phys. Chem. B 2013, 117, 1130−1135. (56) Hammersley, A.; Svensson, S.; Hanfland, M.; Fitch, A.; Hausermann, D. Two-Dimensional Detector Software: from Real Detector to Idealised Image or Two-Theta Scan. High Pressure Res. 1996, 14, 235−248. (57) Hammersley, A. P. FIT2D V10.3 Reference Manual V4.0; European Synchrotron Radiation Facility: 1998. (58) Qui, X.; Thompson, J. W.; Billinge, S. J. L. PDFgetX2: a GUIDriven Program to Obtain the Pair Distribution Function from X-ray Powder Diffraction Data. J. Appl. Crystallogr. 2004, 37, 678. (59) Lopes, J. N. C.; Pádua, A. A. H. Molecular Force Field for Ionic Liquids Composed of Triflate or Bistriflylimide Anions. J. Phys. Chem. B 2004, 108, 16893−16898. (60) Rizzo, R. C.; Jorgensen, W. L. OPLS All-Atom Model for Amines: Resolution of the Amine Hydration Problem. J. Am. Chem. Soc. 1999, 121, 4827−4836. (61) Price, M. L. P.; Ostrovsky, D.; Jorgensen, W. L. Gas-Phase and Liquid-State Properties of Esters, Nitriles, and Nitro Compounds with the OPLS-AA Force Field. J. Comput. Chem. 2001, 22, 1340−1352. (62) Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435−447. (63) van der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A. E.; Berendsen, H. J. C. GROMACS: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701−1718. (64) Nosé, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511−519. (65) Nosé, S. A Molecular Dynamics Method for Simulations in the Canonical Ensemble. Mol. Phys. 1984, 52, 255 −268. (66) Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52, 7182−7190. (67) Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An Nlog(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089−10092. (68) Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577−8593. (69) Prince, E., Ed. International Tables for Crystallography; International Union of Crystallography: 2006; Vol. C. (70) Lorch, E. Neutron Diffraction by Germania, Silica and Radiation-Damaged Silica Glasses. J. Phys. C: Solid State Phys. 1969, 2, 229−237. (71) Du, J.; Benmore, C. J.; Corrales, R.; Hart, R. T.; Weber, J. K. R. A Molecular Dynamics Simulation Interpretation of Neutron and XI
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX
The Journal of Physical Chemistry B
Article
ray Diffraction Measurements on Single Phase Y(2)O(3)-Al(2)O(3) Glasses. J. Phys.: Condens. Matter 2009, 21, 205102. (72) Lado, F. Numerical Fourier Transforms in One, Two, and Three Dimensions for Liquid State Calculations. J. Comput. Phys. 1971, 8, 417−432. (73) Hettige, J. J.; Kashyap, H. K.; Annapureddy, H. V. R.; Margulis, C. J. Anions, the Reporters of Structure in Ionic Liquids. J. Phys. Chem. Lett. 2013, 4, 105−110. (74) MacFarlane, D. R.; Meakin, P.; Sun, J.; Amini, N.; Forsyth, M. Pyrrolidinium Imides: A New Family of Molten Salts and Conductive Plastic Crystal Phases. J. Phys. Chem. B 1999, 103, 4164−4170. (75) Jin, H.; OʼHare, B.; Dong, J.; Arzhantsev, S.; Baker, G. A.; Wishart, J. F.; Benesi, A.; Maroncelli, M. Physical Properties of Ionic Liquids Consisting of the 1-Butyl-3-methylimidazolium Cation with Various Anions and the Bis(trifluoromethylsulfonyl)imide Anion with Various Cations. J. Phys. Chem. B 2008, 112, 81−92. (76) Appetecchi, G. B.; Montanino, M.; Zane, D.; Carewska, M.; Alessandrini, F.; Passerini, S. Effect of the Alkyl Group on the Synthesis and the Electrochemical Properties of N-Alkyl-N-methylpyrrolidinium Bis(trifluoromethanesulfonyl)imide Ionic Liquids. Electrochim. Acta 2009, 54, 1325−1332. (77) Cang, H.; Li, J.; Fayer, M. D. Orientational Dynamics of the Ionic Organic Liquid 1-Ethyl-3-methylimidazolium Nitrate. J. Chem. Phys. 2003, 119, 13017−13023. (78) Li, J.; Wang, I.; Fruchey, K.; Fayer, M. D. Dynamics in Supercooled Ionic Organic Liquids and Mode Coupling Theory Analysis. J. Phys. Chem. A 2006, 110, 10384−10391. (79) Fruchey, K.; Fayer, M. D. Dynamics in Organic Ionic Liquids in Distinct Regions Using Charged and Uncharged Orientational Relaxation Probes. J. Phys. Chem. B 2010, 114, 2840−2845. (80) Nicolau, B. G.; Sturlaugson, A.; Fruchey, K.; Ribeiro, M. C. C.; Fayer, M. D. Room Temperature Ionic Liquid-Lithium Salt Mixtures: Optical Kerr Effect Dynamical Measurements. J. Phys. Chem. B 2010, 114, 8350−8356. (81) Fruchey, K.; Lawler, C. M.; Fayer, M. D. Investigation of Nanostructure in Room Temperature Ionic Liquids Using Electronic Excitation Transfer. J. Phys. Chem. B 2012, 116, 3054−3064. (82) Sturlaugson, A. L.; Fruchey, K. S.; Fayer, M. D. Orientational Dynamics of Room Temperature Ionic Liquid/Water Mixtures: Water-Induced Structure. J. Phys. Chem. B 2012, 116, 1777−1787. (83) Wong, D. B.; Giammanco, C. H.; Fenn, E. E.; Fayer, M. D. Dynamics of Isolated Water Molecules in a Sea of Ions in a Room Temperature Ionic Liquid. J. Phys. Chem. B 2013, 117, 623−635.
J
dx.doi.org/10.1021/jp403518j | J. Phys. Chem. B XXXX, XXX, XXX−XXX