Christian Krekeler

Effect of Anions on Static Orientational Correlations, Hydrogen Bonds, and Dynamics in Ionic Liquids: A Simulational Study

Three different ionic liquids are investigated via atomistic molecular dynamics simulations using the force field of Lopes and PAdua (J. Phys. Chem. B 2006, 110, 19586). In particular, the 1-ethyl-3-methylimidazolium cation EMIM+ is studied in the presence of three different anions, namely, chloride Cl-, tetrafluoroborate BF(4)(-), and bis((trifluoromethyl)sulfonyly)imide TF2N-. In the focus of the present study are the static distributions of anions and cations around a cation as a function of anion size. It is found that the preferred positions of the anions change from being close to the imidazolium hydrogens to being above and below the imidazolium rings. Lifetimes of hydrogen bonds are calculated and found to be of the same order of magnitude as those of pure liquid water and of some small primary alcohols. Three kinds of short-range cation-cation orderings are studied, among which the offset stacking dominates in all of the investigated ionic liquids. The offset stacking becomes weaker from [EMIM][Cl] to [EMIM][BF4] to [EMIM][TF2N]. Further investigation of the dynamical behavior reveals that cations in [EMIM][TF2N] have a slower tumbling motion compared with those in [EMIM][Cl] and [EMIM][BF4] and that pure diffusive behavior can be observed after 1.5 ns for all three systems at temperatures 90 K above the corresponding melting temperatures.

Ionic Charge Reduction and Atomic Partial Charges from First-Principles Calculations of 1,3-Dimethylimidazolium Chloride

We present a detailed calculation of partial charges for the 1,3-dimethylimidazolium chloride ionic liquid. We first analyze MP2 electronic structure calculations and DFT results on isolated ion pairs with various methods of assigning partial charges to the atomic centers. In a second run we analyze the trajectory of a 25 ps long Car-Parrinello MD run of 30 ion pairs under bulk conditions using a charge fitting procedure due to Blöchl. Both, the single ion pair and the bulk system, provide us with a similar total ionic charge considerably less than unity. Especially the liquid state DFT results give convincing evidence for a reduced ionic charge on the ions. The similarity of both results suggest that the delocalization of the Cl charge is due only to local interactions. The relevance of our results is 2-fold; on the one hand they shed light on the basic property of the liquid and its reduced ionic character, and on the other hand, the ab initio derived partial charges provide a fundamental theoretical basis for the recent attempts to use the total ionic charge as an adjustable parameter. Furthermore, all our partial charges are subject to large fluctuations, hinting to the importance of polarization effects.

A comparative study of two classical force fields on statics and dynamics of [EMIM][BF4] investigated via molecular dynamics simulations

The influences of two different commonly employed force fields on statical and dynamical properties of ionic liquids are investigated for [EMIM][BF(4)]. The force fields compared in this work are the one of Canongia Lopes and Padua [J. Phys. Chem. B 110, 19586 (2006)] and that of Liu et al. [J. Phys. Chem. B 108, 12978 (2004)]. Differences in the strengths of hydrogen bonds are found, which are also reflected in the static ion distributions around the cation. Moreover, due to the stronger hydrogen bonding in the force field of Liu et al., the diffusive motions of cations and anions and the rotational behavior of the cations are slower compared with those obtained with the force field of Canongia Lopes and Padua. Both force fields underestimate the zero-field electrical conductivity, while the experimental dielectric constant can be reproduced within the expected statistical error boundaries.

Density functional study of ion hydration for the alkali metal ions (Li+,Na+,K+) and the halide ions (F-,Br-,Cl-)

We performed first principles density functional calculations to study the effect of monovalent ions M+ (M = Li,Na,K) and A- (A = F,Cl,Br) in water with the aim of characterizing the local molecular properties of hydration. For this reason, several ion-water clusters, up to five or six water molecules were considered; such structures were optimized, and the Wannier analysis was then applied to determine the average molecular dipole moment of water. We found that with an increasing number of water molecules, the molecular polarization is determined by the water-water interaction rather than the water-ion interaction, as one would intuitively expect. These results are consistent with those obtained in previous density functional calculations and with other results obtained by employing classical polarizable water models. The main message of this work is that as one increases the number of molecules the average dipole moment of all water molecules and the ones in the first shell tends to the same value as the average of a similar sized cluster of pure water. This supports the use of nonpolarizable classical models of water in classical atomistic simulations.

Study of 1,3-dimethylimidazolium chloride with electronic structure methods and force field approaches.

The 1,3-dimethyl imidazolium chloride [MMIM]Cl is an example of ionic liquid and frequently studied in literature. In this article [MMIM]Cl is studied using an ab initio method [second order Moller-Plesset perturbation theory (MP2), density functional theory (DFT)] and classical force field approach with the aim of looking at some properties on different scales. Selected properties are studied with the different methods and compared to each other. The comparison between the results obtained with MP2 and the DFT approach allows us to comment on the validity of this latter and thus on its employment in larger systems. On the other hand, the comparison between the DFT results and those of the classical approach allows us to test the reproducibility of electrostatic properties by this latter approach. As the results show the used DFT setup is rather satisfactory, while the classical force fields are describing the electrostatic properties in an insufficient way. A revision (improvement) of the classical force fields is at this stage necessary in order to capture the electrostatic properties in a proper way.

Probing spatial locality in ionic liquids with the grand canonical adaptive resolution molecular dynamics technique

We employ the Grand Canonical Adaptive Resolution Simulation (GC-AdResS) molecular dynamics technique to test the spatial locality of the 1-ethyl 3-methyl imidazolium chloride liquid. In GC-AdResS, atomistic details are kept only in an open sub-region of the system while the environment is treated at coarse-grained level; thus, if spatial quantities calculated in such a sub-region agree with the equivalent quantities calculated in a full atomistic simulation, then the atomistic degrees of freedom outside the sub-region play a negligible role. The size of the sub-region fixes the degree of spatial locality of a certain quantity. We show that even for sub-regions whose radius corresponds to the size of a few molecules, spatial properties are reasonably reproduced thus suggesting a higher degree of spatial locality, a hypothesis put forward also by other researchers and that seems to play an important role for the characterization of fundamental properties of a large class of ionic liquids.

Towards open boundary molecular dynamics simulation of ionic liquids

We extend the use of the adaptive resolution (AdResS) method in its grand canonical-like version (GC-AdResS) to the molecular dynamics simulation of 1,3-dimethylimidazolium chloride. We show that the partitioning of the total system in a subsystem of interest with atomistic details and a reservoir of coarse-grained particles leads to satisfactory results. The challenging aspect of this study, compared to previous AdResS simulations, is the presence of charged particles and the necessity of addressing the question about the minimal physical input needed to model the coarse-grained particles in the reservoir. We propose two different approaches and show that in both cases they are sufficient to capture the decisive physical characteristics that allow a valid system-reservoir coupling. The technically satisfactory results pave the way for the multiscale analysis of ionic liquids and truly open boundary molecular simulations.

Adaptive resolution molecular dynamics technique: Down to the essential

We investigate the role of the thermodynamic (TD) force as an essential and sufficient technical ingredient for an efficient and accurate adaptive resolution algorithm. Such a force applied in the coupling region of an adaptive resolution molecular dynamics setup assures thermodynamic equilibrium between atomistically resolved and coarse-grained regions, allowing the proper exchange of molecules. We numerically prove that indeed for systems as relevant as liquid water and 1,3-dimethylimidazolium chloride ionic liquid, the combined action of the TD force and thermostat allows for computationally efficient and numerically accurate simulations, beyond the current capabilities of adaptive resolution setups, which employ switching functions in the coupling region.