Aqueous interfaces with hydrophobic room-temperature ionic liquids: a molecular dynamics study

We report a molecular dynamics study of the interface between water and (macroscopically) water-immiscible room-temperature ionic liquids "ILs", composed of PF6(-) anions and butyl- versus octyl-substituted methylimidazolium+ cations (noted BMI+ and OMI+). Because the parameters used to simulate the pure ILs were found to exaggerate the water/IL mixing, they have been modified by scaling down the atomic charges, leading to better agreement with the experiment. The comparison of [OMI][PF6] versus [BMI][PF6] ILs demonstrates the importance of the N-alkyl substituent on the extent of solvent mixing and on the nature of the interface. With the most hydrophobic [OMI][PF6] liquid, the "bulk" IL phase is dryer than with the [BMI][PF6] liquid. At the interface, the OMI+ cations retain direct contacts with the bulk IL, whereas the more hydrophilic PF6(-) anions gradually dilute in the local water micro-environment and are thus isolated from the "bulk" IL. The interfacial OMI+ cations are ordered with their imidazolium moiety pointing toward the aqueous side and their octyl chains toward the IL side of the interface. With the [BMI][PF6] liquid, the system gradually evolves from an IL-rich to a water-rich medium, leading to an ill-defined interfacial domain with high intersolvent mixing. As a result, the BMI+ cations are isotropically oriented "at the interface". Because the imidazolium cations are more hydrophobic than the PF6(-) anions, the charge distribution at the interface is heterogeneous, leading to a positive electrostatic potential at the interface with the two studied ILs. Mixing-demixing simulations on [BMI][PF6]/water mixtures are also reported, comparing Ewald versus reaction field treatments of electrostatics. Phase separation is very slow (at least 30 ns), in marked contrast with mixtures involving classical organic liquids, which separate in less than 0.5 ns at the microscopic level. The results allow us to better understand the specificity of the aqueous interfaces with hydrophobic ionic liquids, compared with classical organic solvents, which has important implications as far as the mechanism of liquid-liquid ion extraction is concerned.     

Structure of the liquid-vacuum interface of room-temperature ionic liquids: a molecular dynamics study

Molecular dynamics simulations for the liquid-vacuum interface of the ionic liquid 1-ethyl-3-methylimidazolium nitrate (EMIM+/NO3-) were performed for both electronically polarizable and nonpolarizable potential energy surfaces. The interfacial structural properties, such as the oscillation in the number density profile, the orientational ordering, and the local clustering in the interfacial region, were calculated. The simulations with both the polarizable and nonpolarizable model demonstrate the existence of an inhomogeneous interfacial structure normal to the surface layer. It was found for both models that the ethyl tail group on EMIM+ is likely to protrude outward from the surface. In the outmost surface layer, the cation is likely to lie on the surface with the imidazolium ring parallel to the interface, while there is a second region with enhanced density from that in the bulk where the cation preferably slants with the imidazolium ring tending to be perpendicular to the surface. The results also reveal that the electronic polarization effect is important for the ionic liquid interface. It is found that the cation is likely to be segregated at the ionic liquid surface for the polarizable model, while for the nonpolarizable model, the anion is found to be more likely to exhibit such behavior. The surface tension of the polarizable model (58.5 +/- 0.5 mN/m) is much smaller than that of the nonpolarizable model (82.7 +/- 0.6 mN/m), in better agreement with extrapolated experimental measurements on similar ionic liquid systems.     

A comparison of ether- and alkyl-derivatized imidazolium-based room-temperature ionic liquids: a molecular dynamics simulation study

Molecular dynamics simulations of ether-derivatized imidazolium-based room-temperature ionic liquids (EDI-RTILs), [C(5)O(2)mim][TFSI] and [C(5)O(2)mim][BF(4)], have been performed and compared with simulations of alkyl-derivatized analogues (ADI-RTILs). Simulations yield RTIL densities, self-diffusion coefficients and viscosity in excellent agreement with experimental data. Simulations reveal that structure in the EDI-RTILs, quantified by the extent of nanoscale segregation of tails as well as cation-ion and cation-cation correlations, is reduced compared to that observed in the ADI-RTILs. Significant correlation between ether tail oxygen atoms and imidazolium ring hydrogen atoms was observed in the EDI-RTILs. This correlation is primarily intramolecular in origin but has a significant intermolecular component. Competition of ether oxygen atoms with oxygen atoms of TFSI(-) or fluorine atoms of BF(4)(-) for coordination of the ring hydrogen atoms was found to reduce the extent of cation-anion correlation in the EDI-RTILs compared to the ADI-RTILs. The reduction in intermolecular correlation, particularly tail-tail segregation, as well as weakening of cation-anion specific interactions due to the ether tail, may account for the faster dynamics observed in the EDI-RTILs compared to ADI-RTILs.     

Transport properties of room-temperature ionic liquids from classical molecular dynamics

Room-temperature ionic liquids (RTILs) have attracted much attention in the scientific community in the past decade due their novel and highly customizable properties. Nonetheless, their high viscosities pose serious limitations to the use of RTILs in practical applications. To elucidate some of the physical aspects behind transport properties of RTILs, extensive classical molecular dynamics calculations are reported. Here, in particular, bulk viscosities and ionic conductivities of butyl-methyl-imidazole based RTILs are presented over a wide range of temperatures. The dependence of the properties of the liquids on simulation parameters, e.g., system-size effects or the choice of the interaction potential, is analyzed in detail.     

A comparative study of solvation dynamics in room-temperature ionic liquids

The solvation dynamics of ionic liquids have been the subject of many experimental and theoretical studies but remain poorly understood. We analyze these dynamics by modeling the time-resolved fluorescence response of coumarin 153 in two room-temperature ionic liquids: 1-butyl-1-methylpyrrolidinium bromide and 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide. Our results demonstrate that phenomena such as electrostatic screening operate significantly differently in the two liquids, and the relative importance of translational and rovibrational components of the ionic response depends significantly on the character of the ions involved. However, collective motion dominates the response of both ionic liquids, and the qualitative features of this collective behavior are strikingly similar in both cases.     

A molecular dynamics computer simulation study of room-temperature ionic liquids. II. Equilibrium and nonequilibrium solvation dynamics

The molecular dynamics (MD) simulation study of solvation structure and free energetics in 1-ethyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium hexafluorophosphate using a probe solute in the preceding article [Y. Shim, M. Y. Choi and H. J. Kim, J. Chem. Phys. 122, 044510 (2005)] is extended to investigate dynamic properties of these liquids. Solvent fluctuation dynamics near equilibrium are studied via MD and associated time-dependent friction is analyzed via the generalized Langevin equation. Nonequilibrium solvent relaxation following an instantaneous change in the solute charge distribution and accompanying solvent structure reorganization are also investigated. Both equilibrium and nonequilibrium solvation dynamics are characterized by at least two vastly different time scales--a subpicosecond inertial regime followed by a slow diffusive regime. Solvent regions contributing to the subpicosecond nonequilibrium relaxation are found to vary significantly with initial solvation configurations, especially near the solute. If the solvent density near the solute is sufficiently high at the outset of the relaxation, subpicosecond dynamics are mainly governed by the motions of a few ions close to the solute. By contrast, in the case of a low local density, solvent ions located not only close to but also relatively far from the solute participate in the subpicosecond relaxation. Despite this difference, linear response holds reasonably well in both ionic liquids.     

Uranyl and strontium salt solvation in room-temperature ionic liquids. A molecular dynamics investigation

Using molecular dynamics simulations, we compare the solvation of uranyl and strontium nitrates and uranyl chlorides in two room-temperature ionic liquids (ILs): [BMI][PF(6)] based on 1-butyl-3-methylimidazolium(+),PF(6)(-) and [EMI][TCA] based on 1-ethyl-3-methylimidazolium(+),AlCl(4)(-). Both dissociated M(2+),2NO(3)(-) and associated M(NO(3))(2) states of the salts are considered for the two cations, as well as the UO(2)Cl(2) and UO(2)Cl(4)(2)(-) uranyl complexes. In a [BMI][PF(6)] solution, the "naked" UO(2)(2+) and Sr(2+) ions are surrounded by 5.8 and 10.1 F atoms, respectively. The first-shell PF(6)(-) anions rotate markedly during the dynamics and are coordinated, on the average, monodentate to UO(2)(2+) and bidentate to Sr(2+). In an [EMI][TCA] solution, UO(2)(2+) and Sr(2+) coordinate 5.0 and 7.4 Cl atoms of AlCl(4)(-), respectively, which display more restricted motions. Four Cl atoms sit on a least motion pathway of transfer to uranyl, to form the UO(2)Cl(4)(2)(-) complex. The free NO(3)(-) anions and the UO(2)Cl(4)(2)(-) complex are surrounded by imidazolium(+) cations ( approximately 4 and 6-9, respectively). The first shell of the M(NO(3))(2) and UO(2)Cl(2) neutral complexes is mostly completed by the anionic components of the IL, with different contributions depending on the solvent, the M(2+) cation, and its counterions. Insights into energy components of solvation are given for the different systems.     

Solvation of uranium hexachloro complexes in room-temperature ionic liquids. A molecular dynamics investigation in two liquids

We report a molecular dynamics study of the solvation of UCl(6)(-), UCl(6)(2-), and UCl(6)(3-) complexes in the [BMI][Tf(2)N] and [MeBu(3)N][Tf(2)N] ionic liquid cations based on the same anion (bis(trifluoromethylsulfonyl)imide (Tf(2)N-)) and the butyl-3-methyl-imidazolium+ (BMI+) or methyl-tri-n-butyl-ammonium (MeBu(3)N+) cation, respectively. The comparison of two electrostatic models of the complexes (ionic model with -1 charged halides versus quantum mechanically derived charges) yields similar solvation features of a given solute. In the two liquids, the first solvation shell of the complexes is positively charged and evolves from purely cationic in the case of UCl(6)(3-) to a mixture of cations and anions in the case of UCl(6)(-). UCl(6)(3-) is exclusively "coordinated" to BMI+ or MeBu(3)N+ solvent cations that mainly interact via their CH aromatic protons or their N-Me group, respectively. Around the less charged UCl(6)(-) complex, the cations interact via the less polar moieties (butyl chains of BMI+ or MeBu(3)N+) and the anions display nonspecific interactions. In no case does the uranium atom further coordinate solvent ions. According to an energy components analysis, UCl(6)(3-) interacts more attractively with the [BMI][Tf(2)N] liquid than with [MeBu(3)N][Tf(2)N], while UCl(6)(-) does not show any preference, suggesting a significant solvation effect of the redox properties of uranium, also supported by free energy perturbation simulations. The effect of ionic liquid (IL) humidity is investigated by simulating the three complexes in 1:8 water/IL mixtures. In contrast to the case of "naked" ions (e.g., lanthanide(3+), UO2(2+), alkali, or halides), water has little influence on the solvation of the UCl(6)(n-) complexes in the two simulated ILs, as indicated by structural and energy analysis. This is in full agreement with the experimental observations (Nikitenko, S. I.; et al. Inorg. Chem. 2005, 44, 9497).     

A molecular dynamics computer simulation study of room-temperature ionic liquids. I. Equilibrium solvation structure and free energetics

Solvation in 1-ethyl-3-methylmidazolium chloride and in 1-ethyl-3-methylimidazolium hexafluorophosphate near equilibrium is investigated via molecular dynamics computer simulations with diatomic and benzenelike molecules employed as probe solutes. It is found that electrostriction plays an important role in both solvation structure and free energetics. The angular and radial distributions of cations and anions become more structured and their densities near the solute become enhanced as the solute charge separation grows. Due to the enhancement in structural rigidity induced by electrostriction, the force constant associated with solvent configuration fluctuations relevant to charge shift and transfer processes is also found to increase. The effective polarity and reorganization free energies of these ionic liquids are analyzed and compared with those of highly polar acetonitrile. Their screening behavior of electric charges is also investigated.     

Europium(III) and its halides in anhydrous room-temperature imidazolium-based ionic liquids: a combined TRES, EXAFS, and molecular dynamics study

Combining spectroscopic techniques (TRES and EXAFS) and molecular dynamics simulations, we have investigated the state of trivalent europium dissolved in room-temperature ionic liquids (RTILs), as a function of the RTIL anion and in the presence of added chloride anions. The studied RTILs are based on the 1-butyl-3-methyl-imidazolium (Bumim+) cation and differ by their anionic counterparts: BF4-, PF6-, Tf- (triflate, CF3SO3-), and Tf2N- [(CF3SO2)2N-]. The results show the strong influence of the RTIL nature on the first solvation shell of europium and on its complexation with chloride. Depending on the RTIL, europium(III), which was introduced in solution as a triflate salt, is found to be solvated either by RTIL anions only or as neutral undissociated EuTf3 moieties completed by solvent anions. Kinetic effects, related to the viscosity of the RTIL and the nature of the europium salt, also markedly influence the coordination of added Cl- or F- anions to the metal.     

On the diffusion of ferrocenemethanol in room-temperature ionic liquids: an electrochemical study

The electrochemical behaviour of ferrocenemethanol (FcMeOH) has been studied in a range of room-temperature ionic liquids (RTILs) using cyclic voltammetry, chronoamperomery and scanning electrochemical microscopy (SECM). The diffusion coefficient of FcMeOH, measured using chronoamperometry, decreased with increasing RTIL viscosity. Analysis of the mass transport properties of the RTILs revealed that the Stokes-Einstein equation did not apply to our data. The "correlation length" was estimated from diffusion coefficient data and corresponded well to the average size of holes (voids) in the liquid, suggesting that a model in which the diffusing species jumps between holes in the liquid is appropriate in these liquids. Cyclic voltammetry at ultramicroelectrodes demonstrated that the ability to record steady-state voltammograms during ferrocenemethanol oxidation depended on the voltammetric scan rate, the electrode dimensions and the RTIL viscosity. Similarly, the ability to record steady-state SECM feedback approach curves depended on the RTIL viscosity, the SECM tip radius and the tip approach speed. Using 1.3 μm Pt SECM tips, steady-state SECM feedback approach curves were obtained in RTILs, provided that the tip approach speed was low enough to maintain steady-state diffusion at the SECM tip. In the case where tip-induced convection contributed significantly to the SECM tip current, this effect could be accounted for theoretically using mass transport equations that include diffusive and convective terms. Finally, the rate of heterogeneous electron transfer across the electrode/RTIL interface during ferrocenemethanol oxidation was estimated using SECM, and k(0) was at least 0.1 cm s(-1) in one of the least viscous RTILs studied.     

How ion properties determine the stability of a lipase enzyme in ionic liquids: a molecular dynamics study

The influence of eight different ionic liquid (IL) solvents on the stability of the lipase Candida antarctica lipase B (CAL-B) is investigated with molecular dynamics (MD) simulations. Considered ILs contain cations that are based either on imidazolium or guanidinium as well as nitrate, tetrafluoroborate or hexafluorophosphate anions. Partial unfolding of CAL-B is observed during high-temperature MD simulations and related changes of CAL-B regarding its radius of gyration, surface area, secondary structure, amount of solvent close to the backbone and interaction strength with the ILs are evaluated. CAL-B stability is influenced primarily by anions in the order NO(3)(-)? BF(4)(-) < PF(6)(-) of increasing stability, which agrees with experiments. Cations influence protein stability less than anions but still substantially. Long decyl side chains, polar methoxy groups and guanidinium-based cations destabilize CAL-B more than short methyl groups, other non-polar groups and imidazolium-based cations, respectively. Two distinct causes for CAL-B destabilization are identified: a destabilization of the protein surface is facilitated mostly by strong Coulomb interactions of CAL-B with anions that exhibit a localized charge and strong polarization as well as with polar cation groups. Surface instability is characterized by an unraveling of α-helices and an increase of surface area, radius of gyration and protein-IL total interaction strength of CAL-B, all of which describe a destabilization of the folded protein state. On the other hand, a destabilization of the protein core is facilitated when direct core-IL interactions are feasible. This is the case when long alkyl chains are involved or when particularly hydrophobic ILs induce major conformational changes that enable ILs direct access to the protein core. This core instability is characterized by a disintegration of β-sheets, diffusion of ions into CAL-B and increasing protein-IL van der Waals interactions. This process describes a stabilization of the unfolded protein state. Both of these processes reduce the folding free energy and thus destabilize CAL-B. The results of this work clarify the impact of ions on CAL-B stabilization. An extrapolation of the observed trends enables proposing novel ILs in which protein stability could be enhanced further.     

Sonochemistry and sonoluminescence of room-temperature ionic liquids

The sonication of ionic organic liquids leads to decomposition of the liquids. Multibubble sonoluminescence spectra and headgas analysis reveal a variety of decomposition products from the sonolysis of N,N'-dialkylimidazolium ionic liquids. The decomposition is a result of acoustic cavitation, which generates localized hot spots from the implosive collapse of bubbles in the ionic liquids. Despite the negligible vapor pressure of the ionic liquids, reaction still occurs in a heated shell of the bubbles or from microdroplets thrown into the collapsing bubbles.     

Thermophoresis in liquids: a molecular dynamics simulation study

Thermophoresis in liquids is studied by molecular dynamics simulation (MD). A theory is developed that divides the problem in the way consistent with the characteristic scales. MD is then conducted to obtain the solution of each problem, which is to be all combined for macroscopic predictions. It is shown that when the temperature gradient is applied to the nonconducting liquid bath that contains neutral particles, there occurs a pressure gradient tangential to the particle surface at the particle-liquid interface. This may induce the flow in the interfacial region and eventually the particle to move. This applies to the material system that interacts through van der Waals forces and may be a general source of the thermophoresis phenomenon in liquids. The particle velocity is linearly proportional to the temperature gradient. And, in a large part of the given temperature range, the particle motion is in the direction toward the cold end and decreases with respect to the temperature. It is also shown that the particle velocity decreases or even reverses its sign in the lowest limit of the temperature range or with a particle of relatively weak molecular interactions with the liquid. The characteristics of the phenomenon are analyzed in molecular details.     

Picosecond time-resolved fluorescence study on solute-solvent interaction of 2-aminoquinoline in room-temperature ionic liquids: aromaticity of imidazolium-based ionic liquids

Time-resolved fluorescence spectra and fluorescence anisotropy decay of 2-aminoquinoline (2AQ) have been measured in eight room-temperature ionic liquids, including five imidazolium-based aromatic ionic liquids and three nonaromatic ionic liquids. The same experiments have also been carried out in several ordinary molecular liquids for comparison. The observed time-resolved fluorescence spectra indicate the formation of pi-pi aromatic complexes of 2AQ in some of the aromatic ionic liquids but not in the nonaromatic ionic liquids. The fluorescence anisotropy decay data show unusually slow rotational diffusion of 2AQ in the aromatic ionic liquids, suggesting the formation of solute-solvent complexes. The probe 2AQ molecule is likely to be incorporated in the possible local structure of ionic liquids, and hence the anisotropy decays only through the rotation of the whole local structure, making the apparent rotational diffusion of 2AQ slow. The rotational diffusion time decreases rapidly by adding a small amount of acetonitrile to the solution. This observation is interpreted in terms of the local structure formation in the aromatic ionic liquids and its destruction by acetonitrile. No unusual behavior upon addition of acetonitrile has been found for the nonaromatic ionic liquids. It is argued that the aromaticity of the imidazolium cation plays a key role in the local structure formation in imidazolium-based ionic liquids.     

Room temperature ionic liquids containing low water concentrations-a molecular dynamics study

We have performed classical molecular dynamics to study the properties of a water-miscible and a water-immiscible room-temperature ionic liquid when mixed with small quantities of water. The two ionic liquids consist of the same 1-ethyl-3-methylimidazolium ([EMIM]) cation combined with either the boron tetrafluoride ([BF(4)]) or bis(trifluoromethylsulfonyl)imide ([NTf(2)]) anion. It is found that, in both ionic liquids, water clusters of varying sizes are typically hydrogen bonded to two anions with the cation playing a minor role. We also highlight the difficulties of obtaining dynamic quantities such as self-diffusion coefficients from simulations of such viscous systems.     

Ab initio molecular dynamics simulation of ionic liquids

Ab initio Car-Parinnello molecular dynamics is used to simulate the structure and the dynamics of 1-butyl-3-methylimidazolium iodide ([bmim]I) ionic liquid at 300 K. Site-site pair correlation functions reveal that the anion has a strong interaction with any three C-H's of the imidazolium ring. The ring bends over and wraps around the anion such that the two nitrogen atoms take a distance to the anion. Electron donating butyl group contributes the electronic polarization in addition to geometrical (out-of-plane) polarization of the ring due to the liquid environment. This facilitates bending of the ring along the axis passing through nitrogen atoms. The average bending angle depends largely on the alkyl chain length and slightly on the anion type. Redistribution of electron density over the ring caused by the electron donating alkyl group provides additional independent evidence to the instability of lattice structure, hence the low melting point of the ionic liquid. Simulated viscosity and diffusion coefficients of [bmim]I are in quite agreement with the experiments.     

Limiting diffusion coefficients of ionic liquids in water and methanol: a combined experimental and molecular dynamics study

Mutual diffusion coefficients D(12) of the ionic liquids 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([C(2)MIM][NTf(2)]) and [C(4)MIM][NTf(2)] in highly diluted solutions of water and methanol have been measured at different temperatures between 288 K and 313 K using the Taylor dispersion technique. Tracer diffusion coefficients of the two cations [C(2)MIM](+) and [C(4)MIM](+) as well as the anion [NTf(2)](-) in these solutions have been obtained by molecular dynamics (MD) simulations. For our simulations we used well established force fields for the solvents water and methanol and a recently developed force field for imidazolium-based ionic liquid [C(n)MIM][NTf(2)]. Mutual diffusion coefficients D(12) have been calculated from the tracer diffusion coefficients using the Nernst-Hartley equation strictly valid only at low ionic concentration. The agreement between the diffusion coefficients reported in the literature, the experimental data obtained in this work and the MD results is excellent.     

Behavior of electrogenerated bases in room-temperature ionic liquids

The reductive electrochemistry of substituted benzophenones in the aprotic room-temperature ionic liquid (RTIL) 1-butyl-1-methylpyrrolidinium bistriflimide occurs via two consecutive one-electron processes leading to the radical anion and dianion, respectively. The radical anion exhibited electrochemical reversibility at all time-scales whereas the dianion exhibited reversibility at potential sweep rates of >or=10 V s(-1), collectively indicating the absence of strong ion-paring with the RTIL cation. In contrast, reduction in 1-butyl-3-methylimidazolium bistriflimide is complicated by proton-transfer from the [Bmim] cation. At low potential sweep rates, reduction involves a single two-electron process characteristic of either an electrochemical, chemical, electrochemical (ECE) or disproportion-type (DISP1) mechanism. The rate of radical anion protonation in [Bmim] is governed by basicity and conforms to the Hammett free-energy relation. At higher potential sweep rates in [Bmim][NTf2], reduction occurs via two consecutive one-electron processes, giving rise to the partially reversible generation of the radical anion and the irreversible generation of the dianion, respectively. Also, the redox potentials for the reversible parent/radical anion couples were found to be a linear function of Hammett substituent constants in both RTIL media and exhibited effectively equivalent solvent-dependent reaction constants, which are similar to those for reduction in polar molecular solvents such as acetonitrile or alcohols.     

Sum frequency generation study of the room-temperature ionic liquids/quartz interface

The purpose of this investigation is to study the ionic liquid/quartz interface with sum frequency generation vibrational spectroscopy (SFG). SFG spectroscopy was chosen for this study because of its unique ability to yield vibrational spectra of molecules at an interface. Different polarization combinations are used, which probe different susceptibilities, giving SFG the ability to determine molecular orientation at the interface. The ionic liquids used were 1-butyl-3-methylimidazolium tetrafluoroborate, [BMIM][BF(4)], and 1-butyl-3-methylimidazolium hexafluorophosphate, [BMIM][PF(6)]. To determine the influence of the molecular structure and charge on orientation at the interface, neutral, 1-methylimidazole, and 1-butylimidazole were also studied. Raman spectra and depolarization ratios were obtained for neat samples of 1-methylimidazole, 1-butylimidazole, and 1-butyl-3-methylimidazolium tetrafluoroborate recorded from 2700 to 3300 cm(-1). SFG spectra of the 1-methylimidazole/quartz interface showed both methyl and aromatic C-H vibrations. Orientation calculations determined that the ring of the molecule is tilted 45-68 degrees from normal, with the methyl group oriented 32-35 degrees from normal. The SFG spectra of 1-butylimidazole contain several resonances from the alkyl chain with only one weak resonance from the aromatic ring. Orientation calculations suggest that the ring is lying in the plane of the surface with the methyl group pointing 43-47 degrees from normal. The orientation of the [BMIM][PF(6)] ionic liquid was sensitive to trace amounts of water and had to be evacuated to <3 x 10(-5) Torr for the water to be removed. SFG spectra of both ionic liquids were similar, displaying resonances from the alkyl chain as well as the aromatic ring. Orientation analysis suggests the aromatic ring was tilted 45-90 degrees from normal for [BMIM][BF(4)], while the ring for [BMIM][PF(6)] was tilted 38-58 degrees from normal. This suggests the orientation of the molecule is influenced by the size of the anion.