Simulating the vibrational spectra of ionic liquid systems: 1-ethyl-3-methylimidazolium acetate and its mixtures

Infrared Spectroscopy of Li+ Solvation in Diglyme: Ab Initio Molecular Dynamics and Experiment

Infrared (IR) spectra of solutions of the lithium salt LiBF4 in diglyme, CH3O(CH2CH2O)2CH3, are studied via IR spectroscopy and ab initio molecular dynamics (AIMD) simulations. Experiments show that the major effects of LiBF4, compared to neat diglyme, are the appearance of a new broad band in the 250-500 cm-1 frequency region and a broadening and intensity enhancement of the diglyme band in the 900-1150 cm-1 region accompanied by a red-shift. Computational analysis indicates that hindered translational motions of Li+ in its solvation cage are mainly responsible for the new far-IR band, while the changes in the mid-IR are due to Li+-coordination-dependent B-F stretching vibrations of BF4- anions coupled with diglyme vibrations. Molecular motions in these and lower frequency regions are generally correlated, revealing the collective nature of the vibrational dynamics, which involve multiple ions/molecules. Herein, a detailed analysis of these features via AIMD simulations of the spectrum and its components, combined with analysis of the generalized normal modes of the solution components, is presented. Other minor spectral changes as well as diglyme conformational changes induced by the lithium salt are also discussed.

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.

DFT and ab initio molecular dynamics simulation study of the infrared spectrum of the protic ionic liquid 2-hydroxyethylammonium formate

Protic ionic liquids (PILs) typically show a complex band shape in their infrared (IR) spectra in the high-frequency range due to the hydrogen stretching vibrations of functional groups forming rather strong hydrogen bonds (H-bonds). In the low-frequency range, the intermolecular stretching mode of the H-bond leaves a mark in the far-IR spectrum of PILs. In this study, the IR spectrum of the PIL 2-hydroxyethylammonium formate, [HOCH2CH2NH3][HCOO], is investigated in order to identify the different modes that contribute to the high-frequency band shape, i.e. the cation ν(NH), ν(OH), and ν(CH) modes, and the anion ν(CH) mode, as well as the intermolecular mode of the strongest H-bond in the far-IR spectrum. The assignment is validated by quantum chemistry calculations of clusters at the density functional theory (DFT) level for four ionic pairs and by ab initio molecular dynamics (AIMD) simulations of ten ionic pairs. There is good agreement between the vibrational frequencies obtained from DFT and AIMD simulations for both the high- and low-frequency ranges. Based on the calculations, the strong H-bond interaction between the cation -NH3 group and [HCOO]- gives a broad band envelope associated with the ν(NH) mode in the high-frequency range of the IR spectrum on which there are narrower peaks corresponding to the ν(OH) and ν(CH) modes. In the far-IR (FIR) spectrum, the anions' rattling motion gives a broad feature with a maximum at 160 cm-1, while the H-bond's intermolecular NH⋯O stretching mode appears as a peak at 255 cm-1.

Understanding the role of imidazolium-based ionic liquids in the electrochemical CO2 reduction reaction

The development of efficient CO₂ capture and utilization technologies driven by renewable energy sources is mandatory to reduce the impact of climate change. Herein, seven imidazolium-based ionic liquids (ILs) with different anions and cations were tested as catholytes for the CO₂ electrocatalytic reduction to CO over Ag electrode. Relevant activity and stability, but different selectivities for CO₂ reduction or the side H₂ evolution were observed. Density functional theory results show that depending on the IL anions the CO₂ is captured or converted. Acetate anions (being strong Lewis bases) enhance CO₂ capture and H₂ evolution, while fluorinated anions (being weaker Lewis bases) favour the CO₂ electroreduction. Differently from the hydrolytically unstable 1-butyl-3-methylimidazolium tetrafluoroborate, 1-Butyl-3-Methylimidazolium Triflate was the most promising IL, showing the highest Faradaic efficiency to CO (>95%), and up to 8 h of stable operation at high current rates (-20 mA & -60 mA), which opens the way for a prospective process scale-up.

Intermolecular Interactions and Electrochemical Studies on Highly Concentrated Acetate-Based Water-in-Salt and Ionic Liquid Electrolytes

Water-in-salt electrolytes constitute a new class of materials that have distinct properties relative to lower-concentration solutions. A recent approach to further increase the salt concentration and decrease the water content includes the addition of an ionic liquid to a highly concentrated aqueous solution. However, the physicochemical and electrochemical properties of aqueous lithium acetate-1-ethyl-3-methylimidazolium acetate solutions as well as the molecular interactions between electrolyte species have not been characterized. Here, we investigate these properties by evaluation of the ionic conductivity, viscosity, and thermal properties as well as the electrochemical behavior of various electrodes in these electrolytes. The intermolecular interactions are probed by nuclear magnetic resonance and infrared spectroscopies. We find that the addition of the ionic liquid increases the solubility limit of lithium acetate and that with an increase in both acetate salt and ionic liquid concentration in the electrolyte and decrease in water concentration, a strong acetate-water network is formed. The electrochemical stability window increases upon addition of the ionic liquid and reaches a value larger than 5 V for a set of negative Al and positive Ti electrodes in the highest acetate salt/ionic liquid concentration. Preliminary electrochemical charge storage performance measurements of a symmetric device based on two porous carbon electrodes cycled at a current density of 25 mA g^-1 delivered a specific capacitance of 20 F g^-1 with a Coulombic efficiency higher than 99% using a 1.8 V voltage window.

Hydrogen Bonding and Infrared Spectra of Ethyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide/Water Mixtures: A View from Molecular Dynamics Simulations

Simulations of ab initio molecular dynamics have been performed for mixtures of ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) ionic liquid and water. Statistics of donors and acceptors of hydrogen bonds has revealed that with increasing water content, hydrogen bonds between EMIM cations and TFSI anions are replaced by bonds to water molecules. In the mixture of liquids, the total number of bonds (from EMIM cations or water molecules) formed by TFSI acceptors increases. IR spectra obtained from ab initio molecular dynamics trajectories are in good agreement with literature data for ionic liquid/water systems. Analysis of oscillations of individual C-H and O-H bonds has shown correlations between vibrational frequencies and hydrogen bonds formed by an EMIM cation or water molecule and has indicated that the changes in the IR spectrum result from the decreased number of water-water hydrogen bonds in the mixture. The tests of DFTB methodology with tailored parameterizations have yielded reasonably good description of the IR spectrum of bulk water, whereas available parameterizations have failed in satisfactory reproduction of the IR spectrum of EMIM-TFSI/water mixtures in the region above 3000 cm^-1.

Cluster Analysis in Liquids: A Novel Tool in TRAVIS

We present a novel cluster analysis implemented in our open-source software TRAVIS and its application to realistic and complex chemical systems. The underlying algorithm is exclusively based on atom distances. Using a two-dimensional model system, we first introduce different cluster analysis functions and their application to single snapshots and trajectories including periodicity and temporal propagation. Using molecular dynamics simulations of pure water with varying system size, we show that our cluster analysis is size-independent. Furthermore, we observe a similar clustering behavior of pure water in classical and ab initio molecular dynamics simulations, showing that our cluster analysis is universal. In order to emphasize the application to more complex systems and mixtures, we additionally apply the cluster analysis to ab initio molecular dynamics simulations of the [C₂C₁Im][OAc] ionic liquid and its mixture with water. Using that, we show that our cluster analysis is able to analyze the clustering of the individual components in a mixture as well as the clustering of the ionic liquid with water.

Vibrational Circular Dichroism from DFT Molecular Dynamics: The AWV Method

The paper illustrates the Activity Weighted Velocities (AWV) methodology to compute Vibrational Circular Dichroism (VCD) anharmonic spectra from Density Functional Theory (DFT) molecular dynamics. AWV calculates the spectra by the Fourier Transform of the time correlation functions of velocities, weighted by specific observables: the Atomic Polar Tensors (APTs) and the Atomic Axial Tensors (AATs). Indeed, AWV shows to correctly reproduce the experimental spectra for systems in the gas and liquid phases, both in the case of weakly and strongly interacting systems. The comparison with the experimental spectra is striking especially in the fingerprint region, as demonstrated by the three benchmark systems discussed: (1S)-Fenchone in the gas phase, (S)-(−)-Propylene oxide in the liquid phase, and (R)-(−)-2-butanol in the liquid phase. The time evolution of APTs and AATs can be adequately described by a linear combination of the tensors of a small set of appropriate reference structures, strongly reducing the computational cost without compromising accuracy. Additionally, AWV allows the partition of the spectral signal in its molecular components without any expensive postprocessing and any localization of the charge density or the wave function.

Tailoring Intermolecular Interactions Towards High-Performance Thermoelectric Ionogels at Low Humidity

Development of ionic thermoelectric (iTE) materials is of immense interest for efficient heat-to-electricity conversion due to their giant ionic Seebeck coefficient (S_i ), but challenges remain in terms of relatively small S_i at low humidity, poor stretchability, and ambiguous interaction mechanism in ionogels. Herein, a novel ionogel is reported consisting of polyethylene oxide (PEO), polyethylene oxide-polypropylene oxide-polyethylene oxide (P123), and 1-ethyl-3-methylimidazolium acetate (Emim:OAC). By delicately designing the interactions between ions and polymers, the migration of anions is restricted due to their strong binding with the hydroxyl groups of polymers, while the transport of cations is facilitated through segmental motions due to the increased amorphous regions, thereby leading to enlarged diffusion difference between the cations and anions. Moreover, the plasticizing effect of P123 and Emim:OAC can increase the elongation at break. As a consequence, the ionogel exhibits excellent properties including high S_i (18 mV K^-1 at relative humidity of 60%), good ionic conductivity (1.1 mS cm^-1 ), superior stretchability (787%), and high stability (over 80% retention after 600 h). These findings show a promising strategy to obtain multifunctional iTE materials by engineering the intermolecular interactions and demonstrate the great potential of ionogels for harvesting low-grade heat in human-comfortable humidity environments.

Completely Computational Model Setup for Spectroscopic Techniques: The Ab Initio Molecular Dynamics Indirect Hard Modeling Approach

The spectroscopic quantification of mixture compositions usually requires pure compounds and mixtures of known compositions for calibration. Since they are not always available, methods to fill such gaps have evolved, which are, however, not generally applicable. Therefore, calibration can be extremely challenging, especially when multiple unstable species, for example, intermediates, exist in a system. This study presents a new calibration approach that uses ab initio molecular dynamics (AIMD)-simulated spectra to set up and calibrate models for the physics-based spectral analysis method indirect hard modeling (IHM). To demonstrate our approach called AIMD-IHM, we analyze Raman spectra of ternary hydrogen-bonding mixtures of acetone, methanol, and ethanol. The derived AIMD-IHM pure-component models and calibration coefficients are in good agreement with conventionally generated experimental results. The method yields compositions with prediction errors of less than 5% without any experimental calibration input. Our approach can be extended, in principle, to infrared and NMR spectroscopy and allows for the analysis of systems that were hitherto inaccessible to quantitative spectroscopic analysis.

Fe4(OAc)10[EMIM]2: Novel Iron-Based Acetate EMIM Ionic Compound

We synthesized and characterized a novel iron(II) aceto EMIM coordination compound, which has a simplified empirical formula Fe4(OAc)10[EMIM]2, in two different hydration forms: as anhydrous monoclinic compound and triclinic dihydrate Fe4(OAc)10[EMIM]2·2H2O. The dihydrate compound is isostructural with recently reported Mn4(OAc)10[EMIM]2·2H2O, while the anhydrate is a superstructure of the Mn counterpart, suggesting the existence of solid solutions. Both new Fe compounds contain chains of Fe2+ octahedrally coordinated exclusively by acetate groups. The EMIM moieties do not interact directly with the Fe2+ and contribute to the structural framework of the compound through van der Waals forces and C-H···O hydrogen bonds with the acetate anions. The compounds have a melting temperature of ∼94 °C; therefore, they can be considered metal-containing ionic liquids. Differential thermal analysis indicates three endothermic transitions associated with melting, structural rearrangement in the molten state at about 157 °C, and finally, thermal decomposition of the Fe4(OAc)10[EMIM]2. Thermogravimetric analyses indicate an ∼72 wt % mass loss during the decomposition at 280-325 °C. The Fe4(OAc)10[EMIM]2 compounds have higher thermal stability than their Mn counterparts and [EMIM][OAc] but lower compared to iron(II) acetate. Temperature-programmed desorption coupled with mass spectrometry shows that the decomposition pathway of the Fe4(OAc)10[EMIM]2 involves four distinct regimes with peak temperatures at 88, 200, 267, and 345 °C. The main species observed in the decomposition of the compound are CH3, H2O, N2, CO, OC-CH3, OH-CO, H3C-CO-CH3, and H3C-O-CO-CH3. Variable-temperature infrared vibrational spectroscopy indicates that the phase transition at 160-180 °C is associated with a reorientation of the acetate ions, which may lead to a lower interaction with the [EMIM]+ before the decomposition of the Fe4(OAc)10[EMIM]2 upon further heating. The Fe4(OAc)10[EMIM]2 compounds are porous, plausibly capable of accommodating other types of molecules.

Ionic Liquids and Water: Hydrophobicity vs. Hydrophilicity

Many chemical processes rely extensively on organic solvents posing safety and environmental concerns. For a successful transfer of some of those chemical processes and reactions to aqueous media, agents acting as solubilizers, or phase-modifiers, are of central importance. In the present work, the structure of aqueous solutions of several ionic liquid systems capable of forming multiple solubilizing environments were modeled by molecular dynamics simulations. The effect of small aliphatic chains on solutions of hydrophobic 1-alkyl-3-methylimidazolium bis(trifluoromethyl)sulfonylimide ionic liquids (with alkyl = propyl [C₃C₁im][NTf₂], butyl [C₄C₁im][NTf₂] and isobutyl [iC₄C₁im][NTf₂]) are covered first. Next, we focus on the interactions of sulphonate- and carboxylate-based anions with different hydrogenated and perfluorinated alkyl side chains in solutions of [C₂C₁im][C_nF_2n+1SO₃], [C₂C₁im][C_nH_2n+1SO₃], [C₂C₁im][CF₃CO₂] and [C₂C₁im][CH₃CO₂] (n = 1, 4, 8). The last system considered is an ionic liquid completely miscible with water that combines the cation N-methyl-N,N,N-tris(2-hydroxyethyl)ammonium [N_{1 2OH 2OH 2OH}]⁺, with high hydrogen-bonding capability, and the hydrophobic anion [NTf₂]^-. The interplay between short- and long-range interactions, clustering of alkyl and perfluoroalkyl tails, and hydrogen bonding enables a wealth of possibilities in tailoring an ionic liquid solution according to the needs.

Spectroscopic Signatures of Hydrogen-Bonding Motifs in Protonic Ionic Liquid Systems: Insights from Diethylammonium Nitrate in the Solid State

Diethylammonium nitrate, [N0 0 2 2][NO3], and its perdeuterated analogue, [N D D 2 2] [NO3], were structurally characterized and studied by infrared, Raman, and inelastic neutron scattering (INS) spectroscopy. Using these experimental data along with state-of-the-art computational materials modeling, we report unambiguous spectroscopic signatures of hydrogen-bonding interactions between the two counterions. An exhaustive assignment of the spectral features observed with each technique has been provided, and a number of distinct modes related to NH···O dynamics have been identified. We put a particular emphasis on a detailed interpretation of the high-resolution, broadband INS experiments. In particular, the INS data highlight the importance of conformational degrees of freedom within the alkyl chains, a ubiquitous feature of ionic liquid (IL) systems. These findings also enable an in-depth physicochemical understanding of protonic IL systems, a first and necessary step to the tailoring of hydrogen-bonding networks in this important class of materials.

Practical Online Monitoring of Ionic Liquid Fiber Welding Solvent

Ionic liquids (ILs) are becoming important solvents in commerce, but monitoring their purity and performance in industrial applications presents new challenges. Fiber welding technology utilizes ILs to mold and shape natural fibers (cotton, hemp, flax, silk, and wool) into morphologies that are typically attained only using synthetic, petroleum-based non-biodegradable plastics. The result is an atom-efficient process that up-converts fibrous substrates to value-added products and materials. A key aspect of bringing this and other IL-enabled technologies to market relies on efficient monitoring and recycling of IL-based solvents. Implementing online IL quality monitoring enhances the unit economics of these processes. Here, we characterize and report conductivity measurements, refractometry, and ATR-FTIR spectroscopy techniques for online IL monitoring during an industrial fiber welding process. The online analysis enables more efficient recycling of the IL solvent, increasing the process efficiency and product quality.

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.

On the role of water in the hydrogen bond network in DESs: an ab initio molecular dynamics and quantum mechanical study on the urea–betaine system

We herein report an ab initio molecular dynamics study on a natural DES composed of urea and betaine in a 3 : 2 ratio, as a test case for evaluating the water effect.

Synthesis, Characterization, and Crystal Structures of Two New Manganese Aceto EMIM Ionic Compounds with Chains of Mn2+ Ions Coordinated Exclusively by Acetate

We synthesized and determined crystal structures of two manganese(II) aceto EMIM coordination compounds with simplified empirical formulas Mn₄(OAc)₁₀[EMIM]₂ and Mn₄(OAc)₁₀[EMIM]₂·2H₂O. Both compounds feature extended chains of Mn²⁺ octahedrally coordinated exclusively by acetate anions, which has been observed for the first time. The EMIM moieties and water molecules participate in hydrogen bonding with acetate anions but do not directly interact with the metal cation. Both compounds have melting temperatures around 120 °C and can be considered as (non-room-temperature) ionic liquids. The structural arrangement represented by the two title compounds is robust in terms of accommodating other types of cations and allows for tuning of physical properties of the ionic liquid by means of cation substitution. Thermal analysis results obtained using TGA-DSC and VT IR suggest melting phase transitions around 120 °C, followed by structural rearrangement in the molten state taking place around 140-160 °C. Compounds I and II have a higher thermal stability range compared to [EMIM][OAc] ionic liquid, with an onset decomposition temperature above 260 °C.

Solvent Effects and Aggregation Phenomena Studied by Vibrational Optical Activity and Molecular Dynamics: The Case of Pantolactone

Raman and Raman optical activity (ROA), IR, and vibrational circular dichroism (VCD) spectra of (R)- and (S)-pantolactone have been recorded in three solvents. ROA has been employed on water and DMSO solutions, VCD on DMSO and CCl₄ solutions. In the last solvent, monomer-dimer equilibrium is present. Due to the low conformational flexibility of the isolated molecule and to the possibility of aggregation, this compound has been used here to test different protocols for computation of the spectroscopic responses taking into account solvent effects. Molecular dynamics (MD) simulations have been carried out together with statistical clustering methods based on collective variables to extract the structures needed to calculate the spectra. Quantum mechanical DFT calculations based on PCM are compared with approaches based on different representations of the solvent shell (MM or QM level). Appropriate treatment of the solvent permits obtaining of good band-shapes, with the added advantage that the MD analysis allows one to take into account flexibility of dimeric structures justifying the broadness of observed bands and the absence of intense VCD couplets in the carbonyl and OH stretching regions.

The infrared spectra of protic ionic liquids: performances of different computational models to predict hydrogen bonds and conformer evolution

The temperature dependence of the far- and mid-infrared spectrum of two prototypical protic ionic liquids (PILs) sharing a common trialkylammonium cation, but having different anions, is investigated. The exploitation of both the FIR and MIR ranges provides complementary information about the microscopic configurations and the intermolecular interactions, which determine the structure and the properties of ILs. The analysis of the data collected for all the measured frequencies in a wide temperature range reveals several phase transitions and allows the evaluation of the conformer distribution in the different physical states. The difference in the average energy between the H-bonded configurations and the dispersion-governed ones was also determined for the two PILs. Moreover, a computational model for ionic couples based on the ωB97X-D functional and a polar solvent is here successfully exploited for the description of the hydrogen bonding between anion and cation. For the attribution of vibrational lines of the conformers of the cation, the picture based on single ion calculations at the B3LYP level is more valuable and provides better agreement with the experiments.

Computing Bulk Phase Resonance Raman Spectra from ab Initio Molecular Dynamics and Real-Time TDDFT

We present our novel approach for computing resonance Raman (RR) spectra of periodic bulk phase systems from ab initio molecular dynamics, including solvent influence and some anharmonic effects. Based on real-time time-dependent density functional theory, we obtain the RR spectra for all laser wavelengths in one pass. We compute the RR spectrum of uracil in aqueous solution, which is in good agreement with experiment. This is the first simulation of a bulk phase RR spectrum.

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

To investigate the structure and dynamics of small ionic liquid droplets in gas phase, we performed a DFT-based ab initio molecular dynamics study of several 1-ethyl-3-methylimidazolium acetate clusters in vacuum as well as a bulk phase simulation. We introduce an unbiased criterion for average droplet diameter and density. By extrapolation of the droplet densities, we predict the experimental bulk phase density with a deviation of only a few percent. The hydrogen bond geometry between cations and anions is very similar in droplets and bulk, but the hydrogen bond dynamics is significantly slower in the droplets, becoming slower with increasing system size, with hydrogen bond lifetimes up to 2000 ps. From a normal mode analysis of the trajectories, we identify the modes of the ring proton C-H stretching, which are strongly affected by hydrogen bonding. From analyzing these, we find that the hydrogen bond becomes weaker with increasing system size. The cations possess an increased concentration inside the clusters, whereas the anions show an excess concentration on the outside. Almost all anions point towards the droplet center with their carboxylic groups. Ring stacking is found to be a very important structural motif in the droplets (as in the bulk), but side chain interactions are only of minor importance. By using Voronoi tessellation, we define the exposed droplet surface and find that it consists mainly of hydrogen atoms from the cation's and anion's methyl and ethyl groups. Polar atoms are rarely found on the surface, such that the droplets appear completely hydrophobic on the outside.

Finding the best density functional approximation to describe interaction energies and structures of ionic liquids in molecular dynamics studies

Ionic liquids raise interesting but complicated questions for theoretical investigations due to the fact that a number of different inter-molecular interactions, e.g., hydrogen bonding, long-range Coulomb interactions, and dispersion interactions, need to be described properly. Here, we present a detailed study on the ionic liquids ethylammonium nitrate and 1-ethyl-3-methylimidazolium acetate, in which we compare different dispersion corrected density functional approximations to accurate local coupled cluster data in static calculations on ionic liquid clusters. The efficient new composite method B97-3c is tested and has been implemented in CP2K for future studies. Furthermore, tight-binding based approaches which may be used in large scale simulations are assessed. Subsequently, ab initio as well as classical molecular dynamics simulations are conducted and structural analyses are presented in order to shed light on the different short- and long-range structural patterns depending on the method and the system size considered in the simulation. Our results indicate the presence of strong hydrogen bonds in ionic liquids as well as the aggregation of alkyl side chains due to dispersion interactions.

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.

Polar solvent fluctuations drive proton transfer in hydrogen bonded complexes of carboxylic acid with pyridines: NMR, IR and ab initio MD study

We study a series of intermolecular hydrogen-bonded 1 : 1 complexes formed by chloroacetic acid with 19 substituted pyridines and one aliphatic amine dissolved in CD₂Cl₂ at low temperature by ¹H and ¹³C NMR and FTIR spectroscopy. The hydrogen bond geometries in these complexes vary from molecular (O-HN) to zwitterionic (O^-H-N⁺) ones, while NMR spectra show the formation of short strong hydrogen bonds in intermediate cases. Analysis of C[double bond, length as m-dash]O stretching and asymmetric CO₂^- stretching bands in FTIR spectra reveal the presence of proton tautomerism. On the basis of these data, we construct the overall proton transfer pathway. In addition to that, we also study by use of ab initio molecular dynamics the complex formed by chloroacetic acid with 2-methylpyridine, surrounded by 71 CD₂Cl₂ molecules, revealing a dual-maximum distribution of hydrogen bond geometries in solution. The analysis of the calculated trajectory shows that the proton jumps between molecular and zwitterionic forms are indeed driven by dipole-dipole solvent-solute interactions, but the primary cause of the jumps is the formation/breaking of weak CHO bonds from solvent molecules to oxygen atoms of the carboxylate group.

CO2 capture in ionic liquid 1-alkyl-3-methylimidazolium acetate: a concerted mechanism without carbene

Ionic liquids (ILs) provide a promising medium for CO₂ capture. Recently, the family of ILs comprising imidazolium-based cations and acetate anions, such as 1-ethyl-3-methylimidazolium acetate (EMI⁺OAc^-), has been found to react with CO₂ and form carboxylate compounds. N-Heterocyclic carbene (NHC) is widely assumed to be responsible by directly reacting with CO₂ though NHC has not been detected in these ILs. Herein, a computational analysis of CO₂ capture in EMI⁺OAc^- is presented. Quantum chemistry calculations predict that NHC is unstable in a polar environment, suggesting that NHC is not formed in EMI⁺OAc^-. Ab initio molecular dynamics simulations indicate that an EMI⁺ ion "activated" by the approach of a CO₂ molecule can donate its acidic proton to a neighboring OAc^- anion and form a carboxylate compound with the CO₂ molecule. Analysis of this termolecular process indicates that the EMI⁺-to-OAc^- proton transfer and the formation of 1-ethyl-3-methylimidazolium-2-carboxylate occur essentially concurrently. Based on these findings, a novel concerted mechanism that does not involve NHC is proposed for CO₂ capture.

Spin-state dependence of the structural and vibrational properties of solvated iron(ii) polypyridyl complexes from AIMD simulations: aqueous [Fe(bpy)3]Cl2, a case study

LS and HS IR spectra of aqueous [Fe(bpy)₃]²⁺ and corresponding HS–LS difference IR spectrum as obtained from state-of-the-art ab initio molecular dynamics simulations applied to the determination of the structural and vibrational properties of the solvated complex.

A Molecular Level Understanding of Template Effects in Ionic Liquids

The structure-directing or template effect has been invoked several times for ionic liquids to explain the different outcome in material synthesis, namely, different scaffolds or geometrical arrangements with varying ionic liquids. It is obvious to assume that such an effect can originate from the most likely complex microstructure, being present within the ionic liquid itself. In that regard, ionic liquids have already been shown to undergo a nanosegregation into polar and nonpolar phases, which is commonly known and denoted as microheterogeneity. In order to provide detailed insight on the molecular level and to understand the effects rising from this structuring, we performed molecular dynamics simulations on selected very simple model systems composed of 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, considering ethyl, butyl, hexyl, and octyl side chains attached to the cations, mixed with either n-dodecanol or n-butanol. By analyzing snapshots of the simulation boxes and calculating spatial distribution functions, we can visualize that with increasing side chains, the systems show considerable nanosegregation into polar and nonpolar domains. Combined angular and distance distribution functions show that in case of the nanosegregating systems the side chains of the cations are preferentially arranged in a parallel fashion, which indicates a micelle-like structure for the ionic liquids. The alcohol molecules participate in and are, therefore, influenced by this microheterogeneity. It can be shown that in the case of the short IL alkyl side chains, the self-aggregation of the nonpolar units of the alcohols is much stronger, while for the long chain cations, the nonpolar entities of the alcohols are most often connected to the nonpolar units of the ionic liquids. Using our domain analysis tool, we can quantify these observations by tracking the number, size, and shape of the polar and nonpolar entities present in the different investigated systems. The aforementioned combined angular-distance distribution functions reveal a structure-directing effect of the ionic liquids on the alcohol molecules within our simple model systems. The ionic liquids act as template and order the alcohol molecules according to their own structure, resulting in a parallel alignment of the alkyl side chains of the alcohols and ionic liquid cations, with both polar groups being at the same side. These observations show that the microheterogeneous structure of ionic liquids can indeed be applied to order substrates with respect to each other or, for example, to catalysts in a predetermined fashion, opening new possibilities for explaining or enhancing selectivities of chemical reactions in ionic liquids.

Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations

We present our novel methodology for computing Raman optical activity (ROA) spectra of liquid systems from ab initio molecular dynamics (AIMD) simulations. The method is built upon the recent developments to obtain magnetic dipole moments from AIMD and to integrate molecular properties by using radical Voronoi tessellation. These techniques are used to calculate optical activity tensors for large and complex periodic bulk phase systems. Only AIMD simulations are required as input, and no time-consuming perturbation theory is involved. The approach relies only on the total electron density in each time step and can readily be combined with a wide range of electronic structure methods. To the best of our knowledge, these are the first computed ROA spectra for a periodic bulk phase system. As an example, the experimental ROA spectrum of liquid (R)-propylene oxide is reproduced very well.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Vibrational Spectroscopy of Ionic Liquids

Vibrational spectroscopy has continued use as a powerful tool to characterize ionic liquids since the literature on room temperature molten salts experienced the rapid increase in number of publications in the 1990's. In the past years, infrared (IR) and Raman spectroscopies have provided insights on ionic interactions and the resulting liquid structure in ionic liquids. A large body of information is now available concerning vibrational spectra of ionic liquids made of many different combinations of anions and cations, but reviews on this literature are scarce. This review is an attempt at filling this gap. Some basic care needed while recording IR or Raman spectra of ionic liquids is explained. We have reviewed the conceptual basis of theoretical frameworks which have been used to interpret vibrational spectra of ionic liquids, helping the reader to distinguish the scope of application of different methods of calculation. Vibrational frequencies observed in IR and Raman spectra of ionic liquids based on different anions and cations are discussed and eventual disagreements between different sources are critically reviewed. The aim is that the reader can use this information while assigning vibrational spectra of an ionic liquid containing another particular combination of anions and cations. Different applications of IR and Raman spectroscopies are given for both pure ionic liquids and solutions. Further issues addressed in this review are the intermolecular vibrations that are more directly probed by the low-frequency range of IR and Raman spectra and the applications of vibrational spectroscopy in studying phase transitions of ionic liquids.

Symmetry and dynamics of FHF− anion in vacuum, in CD2Cl2 and in CCl4. Ab initio MD study of fluctuating solvent–solute hydrogen and halogen bonds

Asymmetric solvation of FHF⁻ by halogen- and hydrogen-bonding solvents breaks the symmetry of the anion.

Resolving the Origins of Crystalline Anharmonicity Using Terahertz Time-Domain Spectroscopy and ab Initio Simulations

Anharmonicity has been shown to be an important piece of the fundamental framework that dictates numerous observable phenomena. In particular, anharmonicity is the driving force of vibrational relaxation processes, mechanisms that are integral to the proper function of numerous chemical processes. However, elucidating its origins has proven difficult due to experimental and theoretical challenges, specifically related to separating the anharmonic contributions from other unrelated effects. While no one technique is particularly suited for providing a complete picture of anharmonicity, by combining multiple complementary methods such a characterization can be made. In this study the role of individual atomic interactions on the anharmonic properties of crystalline purine, the building block of many DNA and RNA nucleobases, is studied by experimental terahertz time-domain spectroscopy and first-principles density functional theory (DFT) and ab initio molecular dynamics simulations (AIMD). In particular, the detailed vibrational information provided by the DFT calculations is used to interpret the atomic origins of anharmonic-related effects as determined by the AIMD calculations, which are in good agreement with the experimental data. The results highlight that anharmonicity is especially pronounced in the intermolecular interactions, particularly along the amine hydrogen bond coordinate, and yields valuable insight into what is similarly observed complex biosystems and crystalline solids.

Classical Magnetic Dipole Moments for the Simulation of Vibrational Circular Dichroism by ab Initio Molecular Dynamics

We present a new approach for calculating vibrational circular dichroism spectra by ab initio molecular dynamics. In the context of molecular dynamics, these spectra are given by the Fourier transform of the cross-correlation function of magnetic dipole moment and electric dipole moment. We obtain the magnetic dipole moment from the electric current density according to the classical definition. The electric current density is computed by solving a partial differential equation derived from the continuity equation and the condition that eddy currents should be absent. In combination with a radical Voronoi tessellation, this yields an individual magnetic dipole moment for each molecule in a bulk phase simulation. Using the chiral alcohol 2-butanol as an example, we show that experimental spectra are reproduced very well. Our approach requires knowing only the electron density in each simulation step, and it is not restricted to any particular electronic structure method.

Ion pairing in ionic liquids

In the present article we briefly review the extensive discussion in literature about the presence or absence of ion pair-like aggregates in ionic liquids. While some experimental studies point towards the presence of neutral subunits in ionic liquids, many other experiments cannot confirm or even contradict their existence. Ion pairs can be detected directly in the gas phase, but no direct method is available to observe such association behavior in the liquid, and the corresponding indirect experimental proofs are based on such assumptions as unity charges at the ions. However, we have shown by calculating ionic liquid clusters of different sizes that assuming unity charges for ILs is erroneous, because a substantial charge transfer is taking place between the ionic liquid ions that reduce their total charge. Considering these effects might establish a bridge between the contradicting experimental results on this matter. Beside these results, according to molecular dynamics simulations the lifetimes of ion-ion contacts and their joint motions are far too short to verify the existence of neutral units in these materials.