Simulating structure and dynamics in small droplets of 1-ethyl-3-methylimidazolium acetate

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.

Uncovering the Properties of Dicationic Ionic Liquid Nanodroplets through Ab Initio Molecular Dynamics Simulations

The behavior of nanodroplets of an imidazolium-based dicationic ionic liquid, i.e., [C1(mim)2][PF6]2, was investigated in this study using ab initio molecular dynamics simulations. The vibrational features as well as the structural, interfacial, and dynamical properties of different sized droplets were analyzed and compared to the bulk phase system. Structural properties of the droplets, such as π-π stacking, radial distribution functions, structure factors, combined distribution functions, and angular distribution functions were analyzed to understand the interactions and orientations of their ions. The vibrational features and hydrogen bonding strength of droplets were studied by calculating their infrared (IR) and power spectra, determining the contribution of different types of hydrogen bonding to each vibrational mode. The calculated spectra showed good overall agreement with the experimental results. The interfacial properties of the droplets and the orientation of their ions were analyzed using density profiles and an exposed surface. The results showed that, in all systems studied, cations and anions were equally likely to exist in both inner and outer layers, and the cations tended to be oriented toward the center of droplets with obtuse angles. Additionally, the droplet densities were extrapolated to predict the bulk phase density with less than 2% deviation. The dynamical properties of hydrogen bonds, mean square displacement, and van Hove correlations of cations and anions were also analyzed. The results indicated that there was no regular trend in the dynamic properties of droplets with an increasing system size.

Hydrogen Bonds in Lead Halide Perovskites: Insights from Ab Initio Molecular Dynamics

Hydrogen bonds (HBs) play an important role in the rotational dynamics of organic cations in hybrid organic/inorganic halide perovskites, thus affecting the structural and electronic properties of the perovskites. However, the properties and even the existence of HBs in these perovskites are not well established. In this study, we investigate HBs in perovskites MAPbBr3 (MA+ = CH3NH3+), FAPbI3 (FA+ = CH(NH2)2+), and their solid solution with composition (FAPbI3)7/8(MAPbBr3)1/8, using ab initio molecular dynamics and electronic structure calculations. We consider HBs donated by X-H fragments (X = N and C) of the organic cations and accepted by the halides (Y = Br and I) and characterize their properties based on pair distribution functions and on a combined distribution function of the hydrogen-acceptor distance with the donor-hydrogen-acceptor angle. By analyzing these functions, we establish geometrical criteria for HB existence based on the hydrogen-acceptor (H-Y) distance and donor-hydrogen-acceptor angle (X-H-Y). The distance condition is defined as d(H - Y) < 3 Å for N-H-donated HBs and d(H - Y) < 4 Å for C-H-donated HBs. The angular condition is 135° < (X - H - Y) < 180° for both types of HBs. A HB is considered to be formed when both angular and distance conditions are simultaneously satisfied. At the simulated temperature (350 K), the HBs dynamically break and form. We compute the time correlation functions of HB existence and HB lifetimes, which range between 0.1 and 0.3 ps at that temperature. The analysis of HB lifetimes indicates that N-H-Br bonds are relatively stronger than N-H-I bonds, while C-H-Y bonds are weaker, with a minimal influence from the halide and cation. To evaluate the impact of HBs on the vibrational spectra, we present the power spectrum in the region of N-H and C-H stretching modes, comparing them with the normal mode frequencies of isolated cations. We show that the peaks associated with N-H stretching modes in perovskites are redshifted and asymmetrically deformed, while the C-H peaks do not exhibit these effects.

A force field for bio-polymers in ionic liquids (BILFF) - part 2: cellulose in [EMIm][OAc]/water mixtures

We present the extension of our force field BILFF (Bio-Polymers in Ionic Liquids Force Field) to the bio-polymer cellulose. We already published BILFF parameters for mixtures of ionic liquid 1-ethyl-3-methylimidazolium acetate ([EMIm][OAc]) with water. Our all-atom force field focuses on a quantitative reproduction of the hydrogen bonds in the complex mixture of cellulose, [EMIm]⁺, [OAc]^- and water when compared to reference ab initio molecular dynamics (AIMD) simulations. To enhance the sampling, 50 individual AIMD simulations starting from different initial configurations were performed for cellulose in solvent instead of one long simulation, and the resulting averages were used for force field optimization. All cellulose force field parameters were iteratively adjusted starting from the literature force field of W. Damm et al. We were able to obtain a very good agreement with respect to both the microstructure of the reference AIMD simulations and experimental results such as the system density (even at higher temperatures) and the crystal structure. Our new force field allows performing very long simulations of large systems containing cellulose solvated in (aqueous) [EMIm][OAc] with almost ab initio accuracy.

Perspectives in the Computational Modeling of New Generation, Biocompatible Ionic Liquids

In this Perspective, I review the current state of computational simulations on ionic liquids with an emphasis on the recent biocompatible variants. These materials are used here as an example of relatively complex systems that highlights the limits of some of the approaches commonly used to study their structure and dynamics. The source of these limits consists of the coexistence of nontrivial electrostatic, many-body quantum effects, strong hydrogen bonds, and chemical processes affecting the mutual protonation state of the constituent molecular ions. I also provide examples on how it is possible to overcome these problems using suitable simulation paradigms and recently improved techniques that, I expect, will be gradually introduced in the state-of-the-art of computational simulations of ionic liquids.

Modelling biocompatible ionic liquids based on organic acids and amino acids: challenges for computational models and future perspectives

In this short review I shall highlight the basic principle and the difficulties that arise in attempting the computational modeling of seemingly simple systems which hide an unexpected complexity. Biocompatible ionic liquids which are based on the coupling of organic or amino acid anions with metabolic cations such as cholinium are the target of this review. These substances have been the subject of intense research activities in the last few years and have attracted the attention of computational chemists. I shall show that the computational description of these substances is far from trivial and requires the use of sophisticated techniques in order to account for a surprisingly rich chemistry that is due to several phenomena such as polarization, charge transfer, proton transfer equilibria and tautomerization reactions.

Optimized Atomic Partial Charges and Radii Defined by Radical Voronoi Tessellation of Bulk Phase Simulations

We present a novel method for the computation of well-defined optimized atomic partial charges and radii from the total electron density. Our method is based on a two-step radical Voronoi tessellation of the (possibly periodic) system and subsequent integration of the total electron density within each Voronoi cell. First, the total electron density is partitioned into the contributions of each molecule, and subsequently the electron density within each molecule is assigned to the individual atoms using a second set of atomic radii for the radical Voronoi tessellation. The radii are optimized on-the-fly to minimize the fluctuation (variance) of molecular and atomic charges. Therefore, our method is completely free of empirical parameters. As a by-product, two sets of optimized atomic radii are produced in each run, which take into account many specific properties of the system investigated. The application of an on-the-fly interpolation scheme reduces discretization noise in the Voronoi integration. The approach is particularly well suited for the calculation of partial charges in periodic bulk phase systems. We apply the method to five exemplary liquid phase simulations and show how the optimized charges can help to understand the interactions in the systems. Well-known effects such as reduced ion charges below unity in ionic liquid systems are correctly predicted without any tuning, empiricism, or rescaling. We show that the basis set dependence of our method is very small. Only the total electron density is evaluated, and thus, the approach can be combined with any electronic structure method that provides volumetric total electron densities-it is not limited to Hartree-Fock or density functional theory (DFT). We have implemented the method into our open-source software tool TRAVIS.

A force field for bio-polymers in ionic liquids (BILFF) - part 1: [EMIm][OAc]/water mixtures

We present BILFF, a novel force field for bio-polymers in ionic liquids. In the first part of our study, we introduce optimized force field parameters for mixtures of the ionic liquid (IL) 1-ethyl-3-methylimidazolium acetate ([EMIm][OAc]) with water. This imidazolium-based IL is of particular practical importance as it can dissolve significant amounts of cellulose even at room temperature. An understanding of this dissolution process via molecular dynamics simulations requires a quantitative description of the microscopic structure and the strong hydrogen bonds with a method able of simulating at least several dozen nanoseconds, which is the main aim of our novel force field. To reach this goal, we optimize the force field parameters to reproduce radial, spatial, and combined distribution functions, hydrogen bond lifetimes, diffusion coefficients, and several other quantities from reference ab initio molecular dynamics (AIMD) simulations. Non-trivial effects such as dispersion interactions between the side chains and π-π stacking of the cations are reproduced very well. We further validate the force field by comparison to experimental data such as thermal expansion coefficients, bulk modulus, and density at different temperatures, which yields good agreement and correct trends. No other force field with optimized parameters for mixtures of [EMIm][OAc] and water has been presented in the literature yet. Optimized force field parameters for cellulose and other ILs will be published in upcoming articles.

Ionic Liquid Clusters Generated from Electrospray Thrusters: Cold Ion Spectroscopic Signatures of Size-Dependent Acid-Base Interactions

We determine the intramolecular distortions at play in the 2-hydroxyethylhydrazinium nitrate (HEHN) ionic liquid (IL) propellant, which presents the interesting case that the HEH⁺ cation has multiple sites (i.e., hydroxy, primary amine, and secondary ammonium groups) available for H-bonding with the nitrate anion. These interactions are quantified by analyzing the vibrational band patterns displayed by cold cationic clusters, (HEH⁺)_n(NO₃^-)_n-1, n = 2-6, which are obtained using IR photodissociation of the cryogenically cooled, mass-selected ions. The strong interaction involving partial proton transfer of the acidic N-H proton in HEH⁺ cation to the nitrate anion is strongly enhanced in the ternary n = 2 cluster but is suppressed with increasing cluster size. The cluster spectra recover the bands displayed by the bulk liquid by n = 5, thus establishing the minimum domain required to capture this aspect of macroscopic behavior.

Exploring Free Energy Profiles of Enantioselective Organocatalytic Aldol Reactions under Full Solvent Influence

We present a computational study on the enantioselectivity of organocatalytic proline-catalyzed aldol reactions between aldehydes in dimethylformamide (DMF). To explore the free energy surface of the reaction, we apply two-dimensional metadynamics on top of ab initio molecular dynamics (AIMD) simulations with explicit solvent description on the DFT level of theory. We avoid unwanted side reactions by utilizing our newly developed hybrid AIMD (HyAIMD) simulation scheme, which adds a simple force field to the AIMD simulation to prevent unwanted bond breaking and formation. Our condensed phase simulation results are able to nicely reproduce the experimental findings, including the main stereoisomer that is formed, and give a correct qualitative prediction of the change in syn:anti product ratio with different substituents. Furthermore, we give a microscopic explanation for the selectivity. We show that both the explicit description of the solvent and the inclusion of entropic effects are vital to a good outcome-metadynamics simulations in vacuum and static nudged elastic band (NEB) calculations yield significantly worse predictions when compared to the experiment. The approach described here can be applied to a plethora of other enantioselective or organocatalytic reactions, enabling us to tune the catalyst or determine the solvent with the highest stereoselectivity.

Understanding the Complex Surface Interplay for Extraction: A Molecular Dynamics Study

By means of classical molecular dynamics simulation the interfacial properties of methanol and n‐dodecane, which are two potential candidate solvents for use in non‐aqueous liquid–liquid extraction, were assessed. The question of how the interface changes depending on the concentration of extractant (tri‐n‐butyl phosphate) and salt (LiCl) is addressed. Two different models to represent systems were used to evaluate how LiCl and tri‐n‐butyl phosphate affect mutual miscibility, and how the last‐named behaves depending on the chemical environment. Tri‐n‐butyl phosphate increases the mutual solubility of the solvents, whereas LiCl counteracts it. The extractant was found to be mostly adsorbed on the interface between the solvents, and therefore the structural features of the adsorption were investigated. Adsorption of tri‐n‐butyl phosphate changes depending on its concentration and the presence of LiCl. It exhibits a preferential orientation in which the butyl chains point at the n‐dodecane phase and the phosphate group points at the methanol phase. For high concentrations of tri‐n‐butyl phosphate, its molecular orientation is preserved by diffusion of the excess molecules into both the methanol and n‐dodecane phases. However, LiCl hinders the diffusion into the methanol phase, and thus increases the concentration of tri‐n‐butyl phosphate at the interface and forces a rearrangement with subsequent loss of orientation.

Dissolving Cellulose in 1,2,3-Triazolium- and Imidazolium-Based Ionic Liquids with Aromatic Anions

We present 1,2,3-triazolium- and imidazolium-based ionic liquids (ILs) with aromatic anions as a new class of cellulose solvents. The two anions in our study, benzoate and salicylate, possess a lower basicity when compared to acetate and therefore should lead to a lower amount of N-heterocyclic carbenes (NHCs) in the ILs. We characterize their physicochemical properties and find that all of them are liquids at room temperature. By applying force field molecular dynamics (MD) simulations, we investigate the structure and dynamics of the liquids and find strong and long-lived hydrogen bonds, as well as significant π-π stacking between the aromatic anion and cation. Our ILs dissolve up to 8.5 wt.-% cellulose. Via NMR spectroscopy of the solution, we rule out chain degradation or derivatization, even after several weeks at elevated temperature. Based on our MD simulations, we estimate the enthalpy of solvation and derive a simple model for semi-quantitative prediction of cellulose solubility in ILs. With the help of Sankey diagrams, we illustrate the hydrogen bond network topology of the solutions, which is characterized by competing hydrogen bond donors and acceptors. The hydrogen bonds between cellulose and the anions possess average lifetimes in the nanosecond range, which is longer than found in common pure ILs.

An Efficient Lossless Compression Algorithm for Trajectories of Atom Positions and Volumetric Data

We present our newly developed and highly efficient lossless compression algorithm for trajectories of atom positions and volumetric data. The algorithm is designed as a two-step approach. In the first step, efficient polynomial extrapolation schemes reduce the information entropy of the data by exploiting both spatial and temporal continuity. The second step processes the data by a series of transformations (Burrows-Wheeler, move-to-front, run length encoding) and finally compresses the stream with multitable canonical Huffman coding. Our approach reaches a compression ratio of around 15:1 for typical position trajectories in the XYZ format. For volumetric data trajectories in Gaussian Cube format (such as electron density), even a compression ratio of around 35:1 is yielded, which is by far the smallest size of all formats compared here. At the same time, compression and decompression are still reasonably fast for everyday use. The precision of the data can be selected by the user. For storage of the compressed data, we introduce the BQB file format, which is very robust, flexible, and efficient. In contrast to most archiving formats, it allows fast random access to individual trajectory frames. Our method is implemented in C++ and provided as free software under the GNU LGPL license. It has been included in the TRAVIS program package but is also available as stand-alone tool and as a library ("libbqb") for use in other projects.