No Arabic abstract
To understand the non-exponential relaxation associated with solvation dynamics in the ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate, we study power spectra of the fluctuating Franck-Condon energy gap of a diatomic probe solute via molecular dynamics simulations. Results show 1/f dependence in a wide frequency range over 2 to 3 decades, indicating distributed relaxation times. We analyze the memory function and solvation time in the framework of the generalized Langevin equation using a simple model description for the power spectrum. It is found that the crossover frequency toward the white noise plateau is directly related to the time scale for the memory function and thus the solvation time. Specifically, the low crossover frequency observed in the ionic liquid leads to a slowly-decaying tail in its memory function and long solvation time. By contrast, acetonitrile characterized by a high crossover frequency and (near) absence of 1/f behavior in its power spectra shows fast relaxation of the memory function and single-exponential decay of solvation dynamics in the long-time regime.
Computer simulations of (i) a [C12mim][Tf2N] film of nanometric thickness squeezed at kbar pressure by a piecewise parabolic confining potential reveal a mesoscopic in-plane density and composition modulation reminiscent of mesophases seen in 3D samples of the same room-temperature ionic liquid (RTIL). Near 2D confinement, enforced by a high normal load, relatively long aliphatic chains are strictly required for the mesophase formation, as confirmed by computations for two related systems made of (ii) the same [C12mim][Tf2N] adsorbed at a neutral solid surface and (iii) a shorter-chain RTIL ([C4mim][Tf2N]) trapped in the potential well of part i. No in-plane modulation is seen for ii and iii. In case ii, the optimal arrangement of charge and neutral tails is achieved by layering parallel to the surface, while, in case iii, weaker dispersion and packing interactions are unable to bring aliphatic tails together into mesoscopic islands, against overwhelming entropy and Coulomb forces. The onset of in-plane mesophases could greatly affect the properties of long-chain RTILs used as lubricants.
Room Temperature Ionic Liquids (RTILs) have attracted much of the attention of the scientific community in the past decade due the 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 (MD) calculations are reported. 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 and choice of the interaction potential, is analyzed.
Ionic liquids constrained at interfaces or restricted in subnanometric pores are increasingly employed in modern technologies, including energy applications. Understanding the details of their behavior in these conditions is therefore critical. By using molecular dynamics simulation, we clarify theoretically and numerically the effect of confinement at the nanoscale on the static and dynamic properties of an ionic liquid. In particular, we focus on the interplay among the size of the ions, the slit pore width, and the length scale associated to the long-range organization of polar and apolar domains present in the bulk material. By modulating both the temperature and the extent of the confinement, we demonstrate the existence of a complex reentrant phase behavior, including isotropic liquid and liquid-crystal-like phases with different symmetries. We show how these changes impact the relative organization of the ions, with substantial modifications of the Coulombic ordering, and their dynamical state. In this respect, we reveal a remarkable decoupling of the dynamics of the counterions, pointing to very different roles played by these in charge transport under confinement. We finally discuss our findings in connection with very recent experimental and theoretical work.
In V-T theory the atomic motion is harmonic vibrations in a liquid-specific potential energy valley, plus transits, which move the system rapidly among the multitude of such valleys. In its first application to the self intermediate scattering function (SISF), V-T theory produced an accurate account of molecular dynamics (MD) data at all wave numbers q and time t. Recently, analysis of the mean square displacement (MSD) resolved a crossover behavior that was not observed in the SISF study. Our purpose here is to apply the more accurate MSD calibration to the SISF, and assess the results. We derive and discuss the theoretical equations for vibrational and transit contributions to the SISF. The time evolution is divided into three successive intervals: the vibrational interval when the vibrational contribution alone accurately accounts for the MD data; the crossover when the vibrational contribution saturates and the transit contribution becomes resolved; and the diffusive interval when the transit contribution alone accurately accounts for the MD data. The resulting theoretical error is extremely small at all q and t. Comparison of V-T and mode-coupling theories for the MSD and SISF reveals that, while their formulations differ substantially, their underlying atomic motions are in logical correspondence.
We investigate the structure of the [bmim][Tf2N]/silica interface by simulating the indentation of a thin (4 nm) [bmim][Tf2N] film by a hard nanometric tip. The ionic liquid/silica interface is represented in atomistic detail, while the tip is modelled by a spherical mesoscopic particle interacting via an effective short-range potential. Plots of the normal force (Fz) on the tip as a function of its distance from the silica surface highlight the effect of weak layering in the ionic liquid structure, as well as the progressive loss of fluidity in approaching the silica surface. The simulation results for Fz are in near-quantitative agreement with new AFM data measured on the same [bmim][Tf2N]/silica interface at comparable thermodynamic conditions.