Solute Rotation in Ionic Liquids: Size, Shape, and Electrostatic Effects

Do Electrostatics Control the Diffusive Dynamics of Solitary Water? NMR and MD Studies of Water Translation and Rotation in Dipolar and Ionic Solvents

NMR-based measurements of the diffusion coefficients and rotation times of solitary water and benzene at 300 K are reported in a diverse collection of 13 conventional organic solvents and 10 imidazolium ionic liquids. Proton chemical shifts of water are found to be correlated to water OH-stretching frequencies, confirming the importance of electrostatic interactions in these shifts. However, the influence of magnetic interactions in aromatic solvents renders chemical shifts a less reliable indicator of electrostatics. Diffusion coefficients (DB) and rotational correlation times (τB) of benzene in the solvents examined are accurately described as functions of viscosity (η) by DB ∝ η-0.81 and τB ∝ η0.64. Literature values of DB and τB in alkane and normal alcohols, which were not included among the solvents studied here, are systematically faster than predicted by these correlations, indicating that factors beyond solvent viscosity play a role in determining the friction on benzene. In contrast to benzene, water diffusion and rotation are poorly described in terms of viscosity alone, even in the dipolar and ionic solvents measured here. The present data and the substantial literature data already available on dilute water diffusion show a systematic dependence of DW on solvent polarity among isoviscous solvents. The aspect of solvent polarity most relevant to water dynamics is the ability of a solvent to accept hydrogen bonds from water, as conveniently quantified by the frequency of water's OH stretching band, ΔνOH. The friction on translation, ζtr = kBT/DW, and rotation, ζrot = kBTτW, are both well correlated by functions of the form ζ(η, ΔνOH) = a1ηa2 exp (a3ΔνOH), where the ai are adjustable parameters. Molecular dynamics simulations reveal a strong coupling between electrostatic and nonelectrostatic water-solvent interactions, which makes it impossible to dissect the friction on water into additive dielectric and hydrodynamic components. Simulations also provide a tentative explanation for the unusual form of the correlating function ζ(η, ΔνOH), at least in the case of ζrot.

Solvent Role of Ionic Liquids in Fundamental Chemical Reaction Dynamics Analyzed by Time-Resolved Spectroscopy

Ionic liquids (ILs), which are used as solvents for chemical reactions, are different from conventional organic solvents owing to their designability. Physicochemical parameters of the ILs, such as polarity and viscosity, that affect chemical equilibria and reaction kinetics can be tuned by changing the combination of anions and cations or by varying the lengths of the alkyl chains present in the cations. We were interested in knowing how these physicochemical parameters affect fundamental chemical reactions in ILs. Therefore, in this personal account, we investigate our recent work on two different photochemical reactions in ILs, namely excited-state intramolecular proton transfer of hydroxyflavone and photodissociation of aminodisulfide, using time-resolved spectroscopic techniques. Interestingly, the roles of the ILs in these chemical reactions are quite different. The effect of the cationic species of the ILs (i. e., the head groups and number of alkyl carbons) on the solvation environment upon photoexcitation and reaction rate are discussed.

Simulations of electric field gradient fluctuations and dynamics around sodium ions in ionic liquids

The T₁ relaxation time measured in nuclear magnetic resonance experiments contains information about electric field gradient (EFG) fluctuations around a nucleus, but computer simulations are typically required to interpret the underlying dynamics. This study uses classical molecular dynamics (MD) simulations and quantum chemical calculations, to investigate EFG fluctuations around a Na⁺ ion dissolved in the ionic liquid 1-ethyl 3-methylimidazolium tetrafluoroborate, [Im₂₁][BF₄], to provide a framework for future interpretation of NMR experiments. Our calculations demonstrate that the Sternheimer approximation holds for Na⁺ in [Im₂₁][BF₄], and the anti-shielding coefficient is comparable to its value in water. EFG correlation functions, C_EFG(t), calculated using quantum mechanical methods or from force field charges are roughly equivalent after 200 fs, supporting the use of classical MD for estimating T₁ times of monatomic ions in this ionic liquid. The EFG dynamics are strongly bi-modal, with 75%-90% of the de-correlation attributable to inertial solvent motion and the remainder to a highly distributed diffusional processes. Integral relaxation times, ⟨τ_EFG⟩, were found to deviate from hydrodynamic predictions and were non-linearly coupled to solvent viscosity. Further investigation showed that Na⁺ is solvated by four tetrahedrally arranged [BF₄]^- anions and directly coordinated by ∼6 fluorine atoms. Exchange of [BF₄]^- anions is rare on the 25-50 ns timescale and suggests that motion of solvent-shell [BF₄]^- is the primary mechanism for the EFG fluctuations. Different couplings of [BF₄]^- translational and rotational diffusion to viscosity are shown to be the source of the non-hydrodynamic scaling of ⟨τ_EFG⟩.

Recent advances in NMR spectroscopy of ionic liquids

This review presents recent developments in the application of NMR spectroscopic techniques in the study of ionic liquids. NMR has been the primary tool not only for the structural characterization of ionic liquids, but also for the study of dynamics. The presence of a host of NMR active nuclei in ionic liquids permits widespread use of multinuclear NMR experiments. Chemical shifts and multinuclear coupling constants are used routinely for the structure elucidation of ionic liquids and of products formed by their covalent interactions with other materials. Also, the availability of a multitude of NMR techniques has facilitated the study of dynamical processes in them. These include the use of NOESY to study inter-ionic interactions, pulsed-field gradient techniques for probing transport properties, and relaxation measurements to elucidate rotational dynamics. This review will focus on the application of each of these techniques to investigate ionic liquids.

Molecular Dynamics Simulations of Bimolecular Electron Transfer: the Distance-Dependent Electronic Coupling

Understanding the distance dependence of the parameters underpinning Marcus theory is imperative when interpreting the results of experiments on electron transfer (ET). Unfortunately, most of these parameters are difficult or impossible to access directly with experiments, necessitating the use of computer simulations to model them. In this work, we use molecular dynamics simulations in conjunction with constrained density functional theory calculations to study the distance dependence of the electronic coupling matrix element, |H_RP|, for bimolecular ET. Contrary to what is typically assumed for such intermolecular reactions, we find that the magnitude of |H_RP| does not decay exponentially with the center-of-mass separation of the reactants, r_COM. The addition of other simple measures of donor/acceptor (D/A) orientation did not improve the correlation of |H_RP| with r_COM. Using the minimum distance separation, r_min, of the reactants as the structural descriptor allowed the system to be partitioned into high-coupling/close-contact and low-coupling/non-contact regimes, but large fluctuations of |H_RP| were still found for the close-contact reactant pairs. Despite the persistent large fluctuations of |H_RP|, its mean value was found to decay piecewise exponentially with increasing r_min, which was attributed to significant changes in the average D/A pair structure. The results herein advise one to use caution when interpreting the experimental results derived from spherical reactant models of bimolecular ET.

Heterogeneous Structures of Ionic Liquids as Probed by CO Rotation with Nuclear Magnetic Resonance Relaxation Analysis and Molecular Dynamics Simulations

The rotational dynamics of carbon monoxide (CO) in ionic liquids (ILs) was investigated by nuclear magnetic resonance (NMR) relaxation measurements and molecular dynamics (MD) simulations. NMR spin-lattice relaxation time measurements were performed for ¹⁷O-enriched CO in 10 ILs (four imidazolium-cation-based, four phosphonium-cation-based, and two ammonium-cation-based ILs, all paired with the bis(trifluorosulfonylmethane)imide anion). In combination with previously reported data for five ILs and viscosity data, our results indicated that the obtained rotational relaxation times (τ_2R) were much smaller than those predicted using the Stokes-Einstein-Debye (SED) theory. For the same viscosity/temperature values, the τ_2R^-1 value increased linearly with increasing carbon number of the alkyl group in the cation. The deviation from the SED equation was due to the insensitivity of τ_2R to the carbon number, even though a higher carbon number generally leads to higher viscosity values for ILs. To investigate the unique rotational properties of CO in the ILs, MD simulations were performed on five representative ILs (two imidazolium, two phosphonium, and one ammonium) containing CO solutes. From rotational correlation function analyses, the CO rotation mainly occurred in a free rotation-like manner within 1 ps, which explained the relative insensitivity of CO rotation to viscosity. In the subsequent time scale (>1 ps), the minor component of the CO rotation was discriminated among different ILs. It was strongly suggested that, because CO preferably locates in the outer part of the alkyl groups in the cation, the slow CO rotation is correlated with the outer alkyl dynamics, which are decoupled from the whole cation rotation.

Dynamics in the BMIM PF6/acetonitrile mixtures observed by femtosecond optical Kerr effect and molecular dynamics simulations

We have performed the measurements of the optical Kerr effect signal time evolution up to 4 ns for a mixture of 1-alkyl-3-methyl-imidazolium hexafluorophosphate (BMIM PF₆) ionic liquid and acetonitrile in the whole mole fractions range. The long delay line in our experimental setup allowed us to capture the complete reorientational dynamics of the ionic liquid. We have analysed the optical Kerr effect signal in the time and frequency domains with help of molecular dynamics simulations. In our approximation of the slow picosecond dynamics with a multi-exponential decay, we distinguish three relaxation times. The highest two are assigned to the reorientation of the cation and acetonitrile molecules that are in the vicinity of the imidazolium ring. The third one is recognized as originating from cation rotations and reorientation of acetonitrile molecules in the bulk or in the vicinity of the aliphatic chains of the cation. With help of the simulation we interpret the intermolecular band in the reduced spectral density, obtained from Kerr signal, as follows: its low-frequency side results from oscillations of one of the components in the cage formed by its neighbors, while the high-frequency side is attributed to the librations of the cation and acetonitrile molecule as well as the intermolecular oscillations of system components involved in specific interactions. We use this assignment and concentration dependence of the spectra obtained from velocity and angular velocity correlations to explain the mole fraction dependence of Kerr reduced spectral density.

The effect of structural heterogeneity upon the microviscosity of ionic liquids

The behaviour of two molecular rotors, one charged - 3,3'-diethylthiacarbocyanine iodide (Cy3) and one neutral - 8-[4-decyloxyphenyl]-4,4-difluoro-4-bora-3a,4a-diaza-s-indacene (BODIPY-C10), have been studied in various ionic liquids. The fluorescent decay lifetime has been used to elucidate the structure of the immediate region around the rotor. The neutral BODIPY-C10 was found to prefer the non-polar alkyl chain environment, leading to two trends in the lifetime of the dye: one when it was fully partitioned into the non-polar domain, and one when it also sampled polar moieties. The positively charged Cy3 dye showed a complex relationship between the bulk viscosity of the ionic liquid and lifetime of the molecular rotor. This was attributed to a combination of polarity related spectral changes, changes in anion cages around the dye, and temperature dependent fluorescent lifetimes alongside the dependence of the rotor upon the viscosity.

Micelle Heterogeneity from the 2D Kinetics of Solute Rotation

The chemical and physical properties of microstructured materials vary with position. The photophysics of solute molecules can measure these local properties, but they often show multiple rates (rate dispersion), which complicates the interpretation. In the case of micelles, rate dispersion in a solute's anisotropy decay has been assigned to either local anisotropy or heterogeneity in the local viscosity. To resolve this conflict, the rotation of PM597 molecules in SDS micelles has been measured by polarized MUPPETS (multiple population-period transient spectroscopy). This 2D technique shows that heterogeneity is strong and that local anisotropy is minimal. The results suggest that on a subnanosecond time scale, the solute sees only one strong fluctuation of the micelle structure. The anisotropic, average structure emerges on longer time scales.

Dynamics of Dicyanamide in Ionic Liquids is Dominated by Local Interactions

The dynamics of probe molecules is commonly used to investigate the structural dynamics of room-temperature ionic liquids; however, the extent to which this dynamics reflects the dynamics of the ionic liquids or is probe specific has remained debated. Here, we explore to what extent the vibrational and rotational dynamics of the dicyanamide anion, a common ionic liquid anion, correlates with the structural relaxation of ionic liquids. We use polarization-resolved, ultrafast infrared spectroscopy to probe the temperature- and probe-concentration-dependent dynamics of samples with small amounts of 1-ethyl-3-methylimidazolium ([emim]⁺) dicyanamide ([DCA]⁻) dissolved in four [emim]⁺-based ionic liquids with tetrafluoroborate ([BF₄]⁻), bis(trifluoromethylsulfonyl)imide ([NTf₂]⁻), ethylsufate ([EtSO₄]⁻), and triflate ([OTf]⁻) as anions. The transient spectra after broad-band excitation at 2000–2300 cm^–1, resonant with the symmetric and antisymmetric C≡N stretching vibrations, initially contain oscillatory signatures due to the vibrational coherence between both modes. Vibrational population relaxation occurs on two distinct time scales, ∼6–7 and ∼15–20 ps. The vibrational dynamics is rather insensitive to the details of the ionic liquid anion and temperature, except for the slow vibrational relaxation component. The decay of the excitation anisotropy, a measure of the rotational dynamics of [DCA]⁻, markedly depends on temperature, and the obtained decay time exhibits an activation energy of ∼15–21 kJ/mol. Remarkably, neither the rotation time nor the activation energy can be simply explained by the variation of the macroscopic viscosity. Hence, our results suggest that the dynamics of dicyanamide is only in part representative of the ionic liquid structural dynamics. Rather, the dynamics of the probe anion seems to be determined by the specific interaction of [DCA]⁻ with the ionic liquid’s ions for the class of [emim]⁺-based ionic liquids studied here.

Differences in the behavior of dicationic and monocationic ionic liquids as revealed by time resolved-fluorescence, NMR and fluorescence correlation spectroscopy

With an aim to understand the behavior in terms of the intermolecular interactions, structure and dynamics of dicationic and monocationic ionic liquids (ILs), two imidazolium-based dicationic ionic liquids (DILs), 1,8-bis-(3-methylimidazolium-1-yl)octane bis-(trifluoromethylsulfonyl)amide ([C₈(mim)₂][NTf₂]₂), 1,9-bis-(3-methylimidazolium-1-yl)nonane bis-(trifluoromethylsulfonyl)amide ([C₉(mim)₂][NTf₂]₂), and one monocationic ionic liquid (MIL), 1-butyl-3-methyl imidazolium bis(trifluoromethylsulfonyl)amide ([C₄(mim)][NTf₂]), have been investigated through combined fluorescence, electron paramagnetic resonance (EPR), NMR and fluorescence correlation spectroscopy (FCS). The DILs were synthesized by following a standard synthetic protocol and subsequently characterized by different analytical techniques. Steady state absorption, emission and EPR spectroscopic data reveal that DILs are less polar compared to MIL. The polarities of the DILs and MIL were found to be close to those of acetonitrile and short chain alcohols, respectively. The excitation wavelength dependent emission data reveals that DILs are more micro-heterogeneous in nature than MIL. The rotational diffusion of two organic solutes, perylene and 8-methoxypyrene-1,3,6-sulfonate (MPTS), were examined in the DILs and MIL. The rotational diffusion data for perylene and MPTS were analyzed in light of the Stokes-Einstein-Debye (SED) hydrodynamic theory. The rotation of perylene in the DILs was observed to be relatively faster to that in the MIL, and it goes beyond the limit predicted by the SED theory. In order to explain the rotational motion of perylene in DILs, the data was analyzed further by invoking quasi-hydrodynamic theory. The observed rotational behavior of perylene has been explained by considering the fact that perylene is located in the nonpolar region of ILs, and larger solvent molecules (DILs) induce a lower friction to the rotating solute. Interestingly, unlike perylene, rotations of MPTS in both of the ILs were observed to be much hindered indicating a relatively stronger MPTS-IL interaction than perylene-IL interaction. More interestingly, rotation of MPTS was observed to be faster in the DILs than that in the MIL despite the fact that DILs are more viscous than MILs. Relatively faster rotation of MPTS in DILs has been explained by resorting to NMR and FCS studies. The outcomes of the NMR and FCS studies revealed that DILs in the experimental condition exist in their folded form and because of this structural restriction of DILs it becomes difficult for the bulky MPTS to make stronger hydrogen bonding interactions with DILs, which eventually makes the rotation of MPTS in DILs faster. Essentially, the outcomes of all of these studies have demonstrated that the behavior of DILs is quite different to that of the usual MILs.

Photoinduced Bimolecular Electron Transfer in Ionic Liquids: Cationic Electron Donors

Recently, we have reported a systematic study of photoinduced electron-transfer reactions in ionic liquid solvents using neutral and anionic electron donors and a series of cyano-substituted anthracene acceptors [ Wu , B. ; Maroncelli , M. ; Castner , E. W. Jr Photoinduced Bimolecular Electron Transfer in Ionic Liquids . J. Am. Chem. Soc. 139 , 2017 , 14568 ]. Herein, we report complementary results for a cationic class of 1-alkyl-4-dimethylaminopyridinium electron donors. Reductive quenching of cyano-substituted anthracene fluorophores by these cationic quenchers is studied in solutions of acetonitrile and the ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide. Varying the length of the alkyl chain permits tuning of the quencher diffusivities in solution. The observed quenching kinetics are interpreted using a diffusion-reaction analysis. Together with results from the prior study, these results show that the intrinsic electron-transfer rate constant does not depend on the quencher charge in this family of reactions.