Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

Characterizing Microheterogeneity in Liquid Mixtures via Local Density Fluctuations

We present a novel approach to characterize and quantify microheterogeneity and microphase separation in computer simulations of complex liquid mixtures. Our post-processing method is based on local density fluctuations of the different constituents in sampling spheres of varying size. It can be easily applied to both molecular dynamics (MD) and Monte Carlo (MC) simulations, including periodic boundary conditions. Multidimensional correlation of the density distributions yields a clear picture of the domain formation due to the subtle balance of different interactions. We apply our approach to the example of force field molecular dynamics simulations of imidazolium-based ionic liquids with different side chain lengths at different temperatures, namely 1-ethyl-3-methylimidazolium chloride, 1-hexyl-3-methylimidazolium chloride, and 1-decyl-3-methylimidazolium chloride, which are known to form distinct liquid domains. We put the results into the context of existing microheterogeneity analyses and demonstrate the advantages and sensitivity of our novel method. Furthermore, we show how to estimate the configuration entropy from our analysis, and we investigate voids in the system. The analysis has been implemented into our program package TRAVIS and is thus available as free software.

Insight to the Local Structure of Mixtures of Imidazolium-Based Ionic Liquids and Molecular Solvents from Molecular Dynamics Simulations and Voronoi Analysis

While the physicochemical properties as well as the NMR and vibration spectroscopic data of the mixtures of ionic liquids (ILs) with molecular solvents undergo a drastic change around the IL mole fraction of 0.2, the local structure of the mixtures pertaining to this behavior remains unclear. In this work, the local structure of 12 mixtures of 1-butyl-3-methylimidazolium cation (C₄mim⁺) combined with perfluorinated anions, such as tetrafluoroborate (BF₄^-), hexafluorophosphate (PF₆^-), trifluoromethylsulfonate (TFO^-), and bis(trifluoromethanesulfonyl)imide, (TFSI^-), and aprotic dipolar solvents, such as acetonitrile (AN), propylene carbonate (PC), and gamma butyrolactone (γ-BL) is studied by molecular dynamics simulations in the entire composition range, with an emphasis on the IL mole fractions around 0.2. Distributions of metric properties corresponding to the Voronoi polyhedra of the particles (volume assigned to the particles, local density, radius of spherical voids) are determined, using representative sites of the cations, anions, and the solvent molecules, to characterize the changes in the local structure of these mixtures. By analyzing the mole fraction dependence of the average value, fluctuation, and skewness parameter of these distributions, the present study reveals that, around the IL mole fraction of 0.2, the local structure of the mixture undergoes a transition between that determined by the interionic interactions and that determined by the interactions between the ions and solvent molecules. It should be noted that the strength of the interactions between the ions and the solvent molecules, modulated by the change in the composition of the mixture, plays an important role in the occurrence of this transition. The signature of the change in the local structure is traced back to the nonlinear change of the mean values, fluctuations, and skewness values of the metric Voronoi polyhedra distributions.

Direct Correlation between Short-Range Vibrational Spectral Diffusion and Localized Ion-Cage Dynamics of Water-in-Salt Electrolytes

The molecular dynamics simulations of a "water-in-salt" electrolyte, lithium bis(trifluoromethyl sulfonyl) imide (LiNTf₂), with a varying concentration range of 3 to 20 m were performed to establish a direct connection between a dynamic property like the ion-cage lifetime with the short-range vibrational stretching frequency shift of the used probe, HOD. The properties reported here are compared to that obtained from experiments performed at the same concentrations. The time-series wavelet transform was adopted as a preferable mathematical tool for calculating the instantaneous fluctuating frequencies of the probe O-D stretch mode and the concentration-dependent vibrational stretch spectral signature based on the variable functions associated with a particular chemical bond derived from classical molecular dynamics trajectories. The decay time constants of frequency fluctuations and the lifetime of the ion cage (τ_IC) were estimated as a function of salt concentration. Herein, we emphasize the correlation between the slowest time constant (τ₃) of the decay of O-D stretch frequency fluctuations and the timescales associated with the lifetime of ion cages (τ_IC). The results exhibit that the existing relationships were also concentration-dependent. Therefore, this study highlights the connection between the ionic motions that regulate the overall system dynamics with the short-range vibrational frequency shift of the used probe, which was used similar to experiments. It also provides an understanding of the interionic interactions and the dynamical and spectral properties of the electrolytic mixtures. We establish a direct correlation between short-range frequency profile and localized ion-cage lifetime, which can fill the gap of understanding between viscosity, vibrational frequency, and ion-cage dynamics of electrolytes.

Molecular Simulation-Guided Spectroscopy of Imidazolium-Based Ionic Liquids and Effects of Methylation on Ion-Cage and -Pair Dynamics

Classical molecular dynamics simulations were performed to assess an atomistic interpretation of the ion-probe structural interactions in two typical ionic liquids (ILs), 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf₂] and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BDimIm][NTf₂] through computational ultrafast spectroscopy. The nitrile stretching vibrations of the thiocyanate anion, [SCN]^-, serve as the local mode of the ultrafast system dynamics within the imidazolium-based ionic liquid environment. The wavelet transform of classical trajectories determines the time-varying fluctuating frequencies and the stretch spectral signatures of SCN^- in the normalized distribution. However, computational modeling of the two-dimensional (2D) spectra from the wavelet-derived vibrational frequencies yields time evolution of the local molecular structure along with the varied time-dependent dynamics of the spectral diffusion process. We calculated the frequency-frequency correlation functions (FFCFs), time correlations associated with the ion-pair and -cage dynamics, and mean square displacements as a function of time, depicting diffusive dynamics. The calculated results based on the pair correlation functions and the distribution of atomic density suggest that the hydrogen and methylated carbon at the two-position of the imidazolium ring of [BMIm] and [BDimIm] cations, respectively, strongly interact with the probe through the N of the thiocyanate anion rather than the S atom. The center-of-mass center-of-mass (COM-COM) cation-probe radial distribution functions (RDFs) in conjunction with the site-specific structural analysis further reveal well-structured interactions of the thiocyanate ion and [BMIm]⁺ cation rather than the [BDimIm] cation. In contrast, the anion-probe COM-COM RDFs depict weak interactive associations within the vibrational probe [SCN]^- and [NTf₂]^- ions. Methylation at the two-position of the imidazolium ring predicts slower structural reorganization and breaking and reformation dynamics of the ion pairs and cages within the ionic liquid framework.

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.