Search

We use polarizable molecular dynamics simulations to study the thermal dependence of both structural and dynamic properties of two ionic liquids sharing the same cation (1-ethyl-3-methylimidazolium). The linear temperature trend in the structure is accompanied by an exponential Arrhenius-like behavior of the dynamics. Our parameter-free Voronoi tessellation analysis directly casts doubt on common concepts such as the alternating shells of cations and anions and the ionicity. The latter tries to explain the physico-chemical properties of the ionic liquids based on the association and dissociation of an ion pair. However, cations are in the majority of both ion cages, around cations and around anions. There is no preference of a cation for a single anion. Collectivity is a key factor in the dynamic properties of ionic liquids. Consequently, collective rotation relaxes faster than single-particle rotations, and the activation energies for collective translation and rotation are lower than those of the single molecules. [ABSTRACT FROM AUTHOR]

The Kamlet–Taft dipolarity/polarizability parameters π* for various ionic liquids were determined using 4-tert-butyl-2-((dicyanomethylene)-5-[4-N,N-diethylamino)-benzylidene]-Δ3-thiazoline and 5-(N,N-dimethylamino)-5′-nitro-2,2′-bithiophene as solvatochromic probes. In contrast to the established π*-probe N,N-diethylnitroaniline, the chromophores presented here show excellent agreement with polarity measurement using the chemical shift of 129Xe. They do not suffer from additional bathochromic UV/vis shifts caused by hydrogen-bonding resulting in too high π*-values for some ionic liquids. In combination with large sets of various ionic liquids, these new chromophores thereby allow for detailed analysis of the physical significance of π* and the comparison to quantum-mechanical methods. We find that π* correlates strongly with the ratio of molar refractivity to molar volume, and thus with the refractive index. [ABSTRACT FROM AUTHOR]

Different types of spectroscopy capture different aspects of dynamics and different ranges of intermolecular contributions. In this article, we investigate the dielectric relaxation spectroscopy (DRS) of collective nature and the time-dependent Stokes shift (TDSS) of disputed nature. Our computational study of unconfined and confined water clearly demonstrates that the TDSS reflects local, non-collective dynamics. Surprisingly, we found that the reaction field continuum model (RFCM) used to estimate TDSS curves solely from collective DRS spectra correctly transforms collective dynamics to local ones even in cases when the relaxation time trends are quite different. This correct transformation is possible due to structural information available in the DRS amplitude in a Kivelsen–Madden like context. [ABSTRACT FROM AUTHOR]

The inclusion of explicit polarization in molecular dynamics simulation has gained increasing interest during the last several years. An understudied area is the role of polarizability in computer simulations of solvation dynamics around chromophores, particularly for the large solutes used in experimental studies. In this work, we present fully polarizable ground and excited state force fields for the common fluorophores N-methyl-6-oxyquinolium betaine and Coumarin 153. While analyzing the solvation responses in water, methanol, and the highly viscous ionic liquid 1-ethyl-3-methylimidazolium trifluoromethanesulfonate we found that the inclusion of solute polarizability considerably increases the agreement of the obtained Stokes shift relaxation functions with experimental data. Solute polarizability slows down the inertial solvation response in the femtosecond time regime and enables the chromophore to adapt its dipole moment to the environment. Furthermore, the developed chromophore force field reproduces the solute dipole moments in both the electronic ground and excited state in environments ranging from gas phase to highly polar media correctly. Based on these studies it is anticipated that polarizable models of chromophores will lead to an improved understanding of the relationship of their environment to their spectroscopic properties. [ABSTRACT FROM AUTHOR]

The time-dependent Stokes shift (TDSS) has attracted increasing interest for measuring hydration dynamics around biomolecules during the last decades. Its ability to report on hydration dynamics around proteins, however, was questioned recently since the experimental signal stems from both water and protein motion with an unknown ratio of contribution. Using large-scale computer simulations, we examine the ability of the TDSS to capture local hydration dynamics at nine different sites around the protein ubiquitin. By computationally constraining protein motion, it is shown that the remaining water component is meaningful and in line with the picture of a heterogeneous yet overall mobile hydration layer. However, protein contributions are excessively large and cannot be removed in an experimental context, thus obscuring the water component. Consequently, we conclude that the experimental TDSS may not be suitable for the investigation of hydration dynamics around proteins. [ABSTRACT FROM AUTHOR]

The validity of linear response theory (LRT) in computer simulations of solvation dynamics, i.e. the time-dependent Stokes shift, has been debated widely during the last decades. Since the use of LRT is computationally less expensive than the calculation of the true nonequilibrium response, it is often invoked for large systems exhibiting a particularly slow solvation response, e.g. ionic liquids. In the case of ionic liquids, LRT does not only need to capture the correct overall dynamics of the system, but also the contributions and timescales of the respective cation and anion movement. We show by large scale computer simulations that the contribution of the permanent dipoles to the solvation response obeys LRT to some extent, whereas the induced contributions in polarizable simulations lead to a failure of LRT for the respective ion contributions. [ABSTRACT FROM AUTHOR]

We present large scale polarizable simulations of mixtures of the ionic liquids 1-ethyl-3-methylimidazolium trifluoromethanesulfonate and 1-ethyl-3-methylimidazolium dicyanamide with water, where the dielectric spectra, the ion hydration and the conductivity were evaluated. The dielectric decrement, the depression of the dielectric constant of water upon addition of ions, is found to follow a universal functional of Langevin type. Only three physical properties need to be known to describe the complete range of possible concentrations, namely the dielectric constant of pure water, of pure ionic liquid and the linear slope of the dielectric decrement at low ionic liquid concentrations. Both the generalized dielectric constant, as well as the water contribution to the dielectric permittivity follow the functional dependence. We furthermore find that a scaling of van der Waals parameters upon addition of polarizable forces to the force field is necessary to correctly describe the frequency dependent dielectric conductivity and its contribution to the dielectric spectrum, as well as the static electric conductivity, which is also treated in the framework of a pseudolattice theory. [ABSTRACT FROM AUTHOR]

The distribution of a molecule's polarizability to individual atomic sites is inevitable to develop accurate polarizable force fields. We present the direct quantum mechanical calculation of atomic polarizabilities of 27 common ionic liquids. The method is superior to previously published distribution routines based on large databases of the molecular polarizability, and enables the correct description of any ionic liquid and its peculiarities within the quantum mechanical framework. [ABSTRACT FROM AUTHOR]

Ground and excited state dipoles and polarizabilities of the chromophores N-methyl-6-oxyquinolinium betaine (MQ) and coumarin 153 (C153) in solution have been evaluated using time-dependent density functional theory (TD-DFT). A method for determining the atomic polarizabilities has been developed; the molecular dipole has been decomposed into atomic charge transfer and polarizability terms, and variation in the presence of an electric field has been used to evaluate atomic polarizabilities. On excitation, MQ undergoes very site-specific changes in polarizability while C153 shows significantly less variation. We also conclude that MQ cannot be adequately described by standard atomic polarizabilities based on atomic number and hybridization state. Changes in the molecular polarizability of MQ (on excitation) are not representative of the local site-specific changes in atomic polarizability, thus the overall molecular polarizability ratio UGRAPHIC DISPLAY="INLINE" ID="UGT1" SRC="UGT1"/ does not provide a good approximation for local atom-specific polarizability changes on excitation. Accurate excited state force fields are needed for computer simulation of solvation dynamics. The chromophores considered in this study are often used as molecular probes. The methods and data reported here can be used for the construction of polarizable ground and excited state force fields. Atomic and molecular polarizabilities (ground and excited states) have been evaluated over a range of functionals and basis sets. Different mechanisms for including solvation effects have been examined; using a polarizable continuum model, explicit solvation and via sampling of clusters extracted from a MD simulation. A range of different solvents have also been considered. [ABSTRACT FROM AUTHOR]

This study presents the large scale computer simulations of two common fluorophores, N-methyl-6-oxyquinolinium betaine and coumarin 153, in five polar or ionic solvents. The validity of linear response approximations to calculate the time-dependent Stokes shift is evaluated in each system. In most studied systems linear response theory fails. In ionic liquids the magnitude of the overall response is largely overestimated, and linear response theory is not able to capture the individual contributions of cations and anions. In polar liquids, the timescales of solvation dynamics are often not correctly reproduced. These observations are complemented by a detailed analysis of Gaussian statistics including higher order correlation functions, variance of the energy gap distribution and its time evolution. The analysis of higher order correlation functions was found to be not suitable to predict a failure of linear response theory. Further analysis of radial distribution functions and hydrogen bonds in the ground and excited state, as well as the time evolution of the number of hydrogen bonds after solute excitation reveal an influence of solvent structure in some of the studied systems. [ABSTRACT FROM AUTHOR]

The time-dependent fluorescence of a chromophore can be calculated from either nonequilibrium simulations, or, as long as linear response theory holds true, from equilibrium solvent fluctuations in the ground or excited state if the perturbation inflicted by the chromophore is small. The assumption of Gaussian statistics, in contrast, links the nonequilibrium dynamics to solvent fluctuations solely in the excited state, as long as the energy gap distribution is Gaussian throughout the process. The validity of linear response theories on the ground and excited state surface as well as Gaussian statistics is thoroughly tested in this study by calculating the time-dependent Stokes shift of different benzene-like solutes. The effect of the size of change in partial charges of the solute, the multipolar order of charge distribution, the direction of change, as well as the influence of different solvents on the validity of linear response theory is examined by simulating 54 different systems. Calculation of the Gaussian character of the energy distribution in equilibrium, as well as the time-evolution of the peak width in the nonequilibrium simulation sheds light on the validity of Gaussian statistics in a nonstationary regime. We observed that a large intermediate broadening of the width of the energy distribution correlates with a failure of correlation functions to describe the nonequilibrium event. These results are accompanied by analysis of higher order correlation functions, as well as the structure of the solvents water, acetonitrile and methanol around the solute, to yield a comprehensive view, as well as general guidelines, on when and why equilibrium solvent fluctuations can correctly depict solvation dynamics. [ABSTRACT FROM AUTHOR]