Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Molecular ionic liquids are typically characterized by strong electrostatic interactions resulting in a charge ordering and retardation of their translational and rotational behaviour. Unfortunately, this effect is often overestimated in classical molecular dynamics simulations. This can be circumvented in a twofold way: the easiest way is to reduce the partial charges of the ions to sub-integer values of ±0.7-0.9 e. The more realistic model is to include polarizable forces, e.g. Drude-oscillators, but it comes along with an increasing computational effort. On the other hand, charge-scaled models are claimed to take an average polarizability into account. But do both models have the same impact on structure and dynamics of molecular ionic liquids? In the present study several molecular dynamics simulations of 1-ethyl-3-methylimidazolium trifluoromethanesulfonate are performed with different levels of polarization as well as with varying charge scaling factors of 0.74 to 0.90. The analysis of the structural and dynamical results are performed in different levels: from the atomic point of view over the molecular level to collective properties determined by the complete sample.

In this study we present the results of the molecular dynamics simulation of the ionic liquids: 1-butyl-3-methyl-imidazolium tetrafluoroborate and trifluoromethylacetate as well as 1-ethyl-3-methyl-imidazolium dicyanamide. Ionic liquids are characterized by both a molecular dipole moment and a net charge. Thus, in contrast to a solution of simple ions in a (non-) polar solvent, rotational and translational effects influence the very same molecule. This study works out the theoretical framework necessary to compute the conductivity spectrum and its low frequency limit of ionic liquids. Merging these computed conductivity spectra with previous simulation results on the dielectric spectra of ionic liquids yields the spectrum of the generalized dielectric constant, which may be compared to experiments. This spectrum was calculated for the three ionic liquids over six orders of magnitude in frequency ranging from 10 MHz to 50 THz. The role of rotation and translation and their coupling term on the generalized dielectric constant is discussed in detail with a special emphasis on the zero-frequency limit. Thereby, the frequency dependence of the cross correlation between the collective rotational dipole moment and the current is discussed.

The Drude oscillator model is applied to the molecular ionic liquid 1-ethyl-3-methyl-imidazolium triflate. The range of manageable Drude charges is tested. The strength of the polarizability is systematically varied from 0% to 100%. The influence on the structure, single particle dynamics, and collective dielectric properties is investigated. The generalized dielectric constant can be decomposed into a dielectric permittivity, a dielectric conductivity, and an optical dielectric constant ɛ(∞). The major part of the static generalized dielectric constant comes from the collective rotation of the ions, i.e., the dielectric permittivity. The translational contribution from the dielectric conductivity is about 58% of the dielectric permittivity. For the evaluation of the optical dielectric contribution, the computational dielectric theory was adapted to the case of heterogeneous polarizabilities. In case of 100% polarizability, it reaches a value of approximately 2.

The aim of this study is the analysis of the rotational motion in ionic liquids, in particular, 1-butyl-3-methyl-imidazolium tetrafluoroborate. By comparing single-particle and collective motion it is found that the Madden-Kivelson relation is fairly fulfilled in long-term simulation studies (>100 ns), i.e., the collective reorientation can be predicted by the corresponding single-particle property and the static dipolar correlation factor, GK. Furthermore, simulated reorientation is in accordance with hydrodynamic theories yielding hydrodynamic radii comparable to van der Waals radii. Since viscosity is the central quantity entering hydrodynamic formulas, we calculated and measured the viscosity of our system in order to have two independent cycles of hydrodynamic evaluation, a computational and an experimental one. While the static dielectric constant agrees with dielectric reflectance experiment, the hydrodynamic radii derived from the experiments are much lower as a consequence of enhanced rotational motion. Even more, a considerable dynamic broadening is observed in the experiments.

Intramolecular vibrational energy flow in excited bridged azulene-anthracene compounds is investigated by time-resolved pump-probe laser spectroscopy. The bridges consist of molecular chains and are of the type (CH(2))(m) with m up to 6 as well as (CH(2)OCH(2))(n) (n=1,2) and CH(2)SCH(2). After light absorption into the azulene S(1) band and subsequent fast internal conversion, excited molecules are formed where the vibrational energy is localized at the azulene side. The vibrational energy transfer through the molecular bridge to the anthracene side and, finally, to the surrounding medium is followed by probing the red edge of the azulene S(3) absorption band at 300 nm and/or the anthracene S(1) absorption band at 400 nm. In order to separate the time scales for intramolecular and intermolecular energy transfer, most of the experiments were performed in supercritical xenon where vibrational energy transfer to the bath is comparably slow. The intramolecular equilibration proceeds in two steps. About 15%-20% of the excitation energy leaves the azulene side within a short period of 300 fs. This component accompanies the intramolecular vibrational energy redistribution (IVR) within the azulene chromophore and it is caused by dephasing of normal modes contributing to the initial local excitation of the azulene side and extending over large parts of the molecule. Later, IVR in the whole molecule takes place transferring vibrational energy from the azulene through the bridge to the anthracene side and thereby leading to microcanonical equilibrium. The corresponding time constants tau(IVR) for short bridges increase with the chain length. For longer bridges consisting of more than three elements, however, tau(IVR) is constant at around 4-5 ps. Comparison with molecular dynamics simulations suggests that the coupling of these chains to the two chromophores limits the rate of intramolecular vibrational energy transfer. Inside the bridges the energy transport is essentially ballistic and, therefore, tau(IVR) is independent on the length.

In this work, the collective structure of aqueous solutions of ionic liquids was studied by means of molecular dynamics simulations. Various concentrations of 1-butyl-3-methyl-imidazolium tetrafluoroborate and TIP3P water were simulated at the very same size of the simulation box. For the analysis, the ternary system cation/anion/water was subdivided into binary networks. The local structure of each of these six networks is investigated by atom-atom radial distribution functions as well as by the so-called g coefficients, which reveal the mutual orientation of the network constituting partners. Furthermore, the collective structure of the whole samples was characterized by the contribution of each species to the static dielectric constant epsilon(omega=0) and to the Kirkwood G(K) factor. The combination of the analysis tools mentioned above provides knowledge about the cross-linking of the ionic species with the dipolar water. Thereby, the interplay between charge-charge and hydrogen bond networks is analyzed in detail.

In the last two decades, while searching for interesting applications of ionic liquids as potent solvents, their solvation properties and their general impact on biomolecules, and in particular on proteins, gained interest. It turned out that ionic liquids are excellent solvents for protein refolding and crystallization. Biomolecules showed increased solubilities and stabilities, both operational and thermal, in ionic liquids, which also seem to prevent self-aggregation during solubilization. Biomolecules can be immobilized, e.g. in highly viscous ionic liquids, for particular biochemical processes and can be designed to some extent by the proper choice of the ionic liquid cations and anions, which can be characterized by the Hofmeister series.