Virtual Site OPLS Force Field for Imidazolium-Based Ionic Liquids

Monte Carlo QM/MM simulation studies of the Cannizzaro reaction in ionic liquids for improved biofuel production

The conversion of biomass to 5-hydroxymethylfurfural (HMF) holds substantial promise as a renewable energy source. Notably, HMF can be transformed into 2,5-bis(hydroxymethyl)furan (BHMF), a crucial reactant in biofuel production, but requires harsh operating conditions, H2, and precious metal catalysts. A recently reported Cannizzaro reaction of HMF to BHMF, characterized by its efficiency, mild conditions, and eco-friendliness, instead employed ionic liquids (ILs) to achieve high yields. In this study, combined quantum mechanical and molecular mechanical (QM/MM) simulations in conjunction with Metropolis Monte Carlo statistical mechanics and free-energy perturbation theory utilized M06-2X/6-31+G(d), PDDG/PM3, and the OPLS-VSIL force field to uncover important solute-solvent interactions present in the HMF to BHMF reaction pathway. The Cannizzaro reaction was examined in water and in five ILs composed of the 1-butyl-3-methylimidazolium [BMIM] cation coupled to hexafluorophosphate, tetrafluoroborate, thiocyanate, chloride, and bromide. Energetic and structural analysis of the rate-determining hydride transfer between HMF and the hydride-donor anion HMFOH- attributed the enhanced reactivity to highly organized solvent interactions featuring (1) hydrogen bonding between the ring protons of [BMIM] and the negatively charged carbonyl oxygen atoms on the transition structure, (2) favorable electrostatic interactions between the IL anions and solute hydroxyl groups, and (3) beneficial π-π stacking interactions between [BMIM] and the two aromatic rings present in HMF and HMFOH-.

Competitive Effects of Anions on Protein Solvation by Aqueous Ionic Liquids

The present study utilizes molecular dynamics simulations to examine how different anions compete for protein solvation in aqueous solutions of ionic liquids (ILs). Ubiquitin is used as model protein and studied in IL mixtures sharing the same cation, 1-ethyl-3-methylimidazolium (EMIM), and two different anions in the same solution, from combinations of dicyanamide (DCA), chloride (Cl), nitrate (NO3), and tetrafluoroborate (BF4). Our findings reveal that specific interactions between anions and the protein are paramount in IL solvation, but that combinations of anions are not additive. For example, DCA exhibits a remarkable ability to form hydrogen bonds with the protein, resulting in a significantly stronger preferential binding to the protein than other anions. However, the combination of DCA with NO3, which also forms hydrogen bonds with the protein, results in a smaller preferential solvation of the protein than the combination of DCA with chloride ions, which are weaker binders. Thus, combining anions with varying affinities for the protein surface modulates the overall ion accumulation through nonadditive mechanisms, highlighting the importance of the understanding of competition for specific interaction sites, cooperative binding, bulk-solution affinity, and overall charge compensations, on the overall solvation capacity of the solution. Such knowledge may allow for the design of novel IL-based processes in biotechnology and material science, where fine-tuning protein solvation is crucial for optimizing performance and functionality.

Assessing OPLS-based force fields for investigating the characteristics of imidazolium-based dicationic ionic liquids: A comparative study with AIMD simulations and experimental findings

In this study, we extended the optimized potentials for liquid simulation-ionic-liquid virtual site (OPLS-VSIL) force field (FF) to imidazolium-based dicationic ionic liquids (DILs) and evaluated the ability of different OPLS-based FFs (i.e., OPLS-2009IL, 0.8*OPLS-2009IL, and OPLS-VSIL) in predicting different properties of the studied DIL by comparing their results with ab initio molecular dynamics (AIMD) simulation and experimental results. To achieve this purpose, MD simulations with three different OPLS-based FFs as well as AIMD simulation were performed for [C3(mim)2][NTF2]2 DIL and its structural, dynamical, vibrational, and volumetric properties were analyzed. Structural properties of the studied DIL, i.e., radial distribution functions (RDFs), structure factor, and hydrogen-bond network, showed that compared to 0.8*OPLS-2009IL FF, there is a much better agreement between the results of both OPLS-2009IL and OPLS-VSIL FFs with the AIMD simulation. On the other hand, the results of dynamical properties, such as mean square displacements, van Hove correlation functions as well as hydrogen bond, ion pair, and ion cage dynamics, depicted that in both 0.8*OPLS-2009IL and OPLS-VSIL FFs, the dynamics of the system is almost similar, and compared to OPLS-2009IL FF, they have better agreements with experimental results where they exist. So, it can be seen that although reducing the total charge of studied DIL by 20% leads to an increase in the dynamics of the system, the type distribution of partial charges on each atom does not significantly affect the system's dynamics. The calculated infrared (IR) and power spectra showed that the vibrational features of studied DIL in three OPLS-based FFs are mostly the same and reducing total charge and different type distribution of partial charges have no significant effect on the studied system. Furthermore, in volumetric properties, OPLS-VSIL FF shows somehow better agreement with experimental results. Overall, the evaluation of different structural, dynamical, vibrational, and volumetric properties of [C3(mim)2][NTF2]2 DIL shows that the OPLS-VSIL FF may be the best choice among the different studied OPLS FFs.

Molecular dynamics study of fluorosulfonyl ionic liquids as electrolyte for electrical double layer capacitors

The development of high-performance supercapacitors is an important goal in the field of energy storage. Ionic liquids (ILs) are promising electrolyte materials for efficient energy storage in supercapacitors, because of the high stability, low volatility, and wider electrochemical stability window than traditional electrolytes. However, ILs-based supercapacitors usually show a relatively lower power density owing to the inherent viscosity-induced low electrical conductivity. Fluorosulfonyl ILs have aroused much attention in energy storage devices due to its low toxicity and excellent stability. Here, we propose that structural modification is an effective way to improve the energy storage performance of fluorosulfonyl ILs through the classical molecular dynamics (MD) method. Four fluorosulfonyl ILs with different sizes and symmetries were considered. Series of properties including conductivity, interface structure, and double-layer capacitance curves were systematically investigated. The results show that smaller size and more asymmetric structure can enhance self-diffusion coefficient and conductivity, and improve the electrochemical performance. Appropriate modification of the electrodes can further enhance the capacitive performance. Our work provides an opportunity to further understand and develop the fluorosulfonyl ILs electrolyte in supercapacitors.

Uncovering the Critical Factors that Enable Extractive Desulfurization of Fuels in Ionic Liquids and Deep Eutectic Solvents from Simulations

Environmental regulatory agencies have implemented stringent restrictions on the permissible levels of sulfur compounds in fuel to reduce harmful emissions and improve air quality. Problematically, traditional desulfurization methods have shown low effectiveness in the removal of refractory sulfur compounds, e.g., thiophene (TS), dibenzothiophene (DBT), and 4-methyldibenzothiophene (MDBT). In this work, molecular dynamics (MD) simulations and free energy perturbation (FEP) have been applied to investigate the use of ionic liquids (ILs) and deep eutectic solvents (DESs) as efficient TS/DBT/MDBT extractants. For the IL simulations, the selected cation was 1-butyl-3-methylimidazolium [BMIM] and the anions included chloride [Cl], thiocyanate [SCN], tetrafluoroborate [BF₄], hexafluorophosphate [PF₆], and bis(trifluoromethylsulfonyl)amide [NTf₂]. The DESs were composed of choline chloride with ethylene glycol (CCEtg) or with glycerol (CCGly). Calculation of excess chemical potentials predicted the ILs to be more promising extractants with energies lower by 1-3 kcal/mol compared to DESs. Increasing IL anion size was positively correlated to enhanced solvation of S-compounds, which was influenced by energetically dominant solute-anion interactions and favorable solute-[BMIM] π-π stacking. For the DESs, the solvent components offered a range of synergistic, yet comparatively weaker, electrostatic interactions that included hydrogen bonding and cation-π interactions. An in-depth analysis of the structure of IL and DES systems is presented, along with a discussion of the critical factors behind experimental trends of S-compound extraction efficiency.

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.

The effect of various partial atomic charges on the bulk and liquid/vacuum interface properties of iodobenzene derivatives at their melting points

In this work, various atomic charge schemes including natural bond orbital (NBO), electrostatic potential based (CHELPG), and σ-hole model charges were applied in the OPLS-AA force field to investigate their effects on the thermophysical and structural properties of iodobenzene and its derivatives. Molecular dynamics simulations presented in this work show that the studied structural and thermophysical properties are in good agreement with experiments when the CHELPG charge was coupled with the OPLS-AA force field. Also, the arrangement of iodobenzene derivatives in the liquid phase was investigated via combined radial/angular distribution functions (CDFs) analyses and halogen-bonding theory. The most probable orientation of iodobenzene derivatives at the liquid/vacuum interface was assigned by atom density profile and bivariate orientational distribution maps. For all studied iodobenzene derivatives, benzene rings are oriented such that the iodine atoms tend toward the vacuum phase.

Ionic Liquids as Organocatalysts for Nucleophilic Fluorination: Concepts and Perspectives

Besides their extremely useful properties as solvent, ionic liquids (ILs) are now considered to be highly instructive tools for enhancing the rates of chemical reactions. The ionic nature of the IL anion and cation seems to be the origin of this fascinating function of ILs as organocatalyst/promoter through their strong Coulombic forces on other ionic species in the reaction and also through the formation of hydrogen bonds with various functional groups in substrates. It is now possible to tailor-make ILs for specific purposes as solvent/promoters in a variety of situations by carefully monitoring these interactions. Despite the enormous potentiality, it seems that the application of ILs as organocatalysts/promoters for chemical reactions have not been fully achieved so far. Herein, we review recent developments of ILs for promoting the nucleophilic reactions, focusing on fluorination. Various aspects of the processes, such as organocatalytic capability, reaction mechanisms and salt effects, are discussed.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

Hydrogen-bonding and "π-π" interaction promoted solution-processable mixed matrix membranes for aromatic amines detection

Detection of hazardous compounds can alleviate risk to human health. However, it remains a challenge to develop easy-to-use testing tools for carcinogenic aromatic amines. Herein, we presented a conjugated molecule-based aniline detector, mixed matrix membranes (MMMs), through the solution-processable strategy. The pentacene-based dispersed phase is achieved using the state-of-the-art ionic liquids (ILs) as the continuous phase, based on which MMMs are easily manufactured by a solution process. Moreover, molecular dynamics (MD) simulations and quantum mechanical calculations suggested that hydrogen bonding and π-π interaction between ILs cations and pentacene could promote the dissolution. These prepared MMMs can offer easy-operation and on-site detection of carcinogenic primary aromatic amines with eye-readable fluorescence signal. This work provides a paradigm for the design of a portable testing device for various hazardous compounds.

Molecular-Level Understanding of Surface Roughness Boosting Segregation Behavior at the ZIF-8/Ionic Liquid Interfaces

Here, we perform a series of classical molecular dynamics simulations for two different [HEMIM][DCA] and [BMIM][BF₄] ionic liquids (ILs) on the ZIF-8 surface to explore the interfacial properties of metal-organic framework (MOFs)/IL composite materials at the molecular level. Our simulation results reveal that the interfacial structures of anions and cations on the ZIF-8 surface are dominated by the surface roughness due to the steric hindrance, which is extremely different from the driving mechanism based on solid-ion interactions of ILs on flat solid surfaces. At the ZIF-8/IL interfaces, the open sodalite (SOD) cages of the ZIF-8 surface can block most of the large-size cations outside and significantly boost the segregation behavior of anions and cations. In comparison with the [BMIM][BF₄] IL, the [HEMIM][DCA] IL has much more anions entering into the open SOD cages owing to the combination of stronger ZIF-8-[DCA]^- interactions and more ordered arrangement of [DCA]^- anions on the ZIF-8 surface. Furthermore, more and stronger ZIF-8-[BF₄]^- hydrogen bonds (HBs) are found to exist on the cage edges than the ZIF-8-[DCA]^- HBs, further preventing [BF₄]^- anions from entering into SOD cages. By more detailed analyses, we find that the hydrophobic interaction has an important influence on the interfacial structures of the side chains of [HEMIM]⁺ and [BMIM]⁺ cations, while the π-π stacking interaction plays a key role in determining the interfacial structures of the imidazolium rings of both cations. Our simulation results in this work provide a molecular-level understanding of the underlying driving mechanism on segregation behavior at the ZIF-8/IL interfaces.

Ion Transport in the EMITFSI/PVDF System at Different Temperatures: A Molecular Dynamics Simulation

We used all-atom molecular dynamics simulations to study the ion transport in the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide/poly(vinylidene fluoride) (EMITFSI/PVDF) system with 40.05 wt % EMITFSI at different temperatures. The glass-transition temperature (T g = 204 K) of this system shows a good agreement with the experimental value (200 K). With the increase of temperature, the peaks of the pair correlation function show an increasing trend. Interestingly, the coordination numbers of ion pairs and the degree of independent ion motion are mainly affected by the binding energy between ion pairs as the temperature increases. In addition, the ion transport properties with increasing temperature can be studied by the ion-pair relaxation times, ion-pair lifetimes, and diffusion coefficients. The simulation results illustrate that the ion transport is intensified. Especially, the cations can always diffuse faster than the anions. The power law shows that mobilities of anions and cations are seen to exhibit a "superionic" behavior. With the increase of temperature, transference numbers of anions decrease first and then increase and transference numbers of cations show the opposite changes; ionic conductivity increases gradually; and viscosity decreases gradually, indicating that the diffusion resistance of ions decreases. In general, after adding PVDF into the EMITFSI system, the glass-transition temperature and viscosity increase, the ionic conductivity and degree of independent ion motion decrease, and diffusion coefficients of cations decrease faster than those of the anions.

Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective

The theoretical treatment of ionic liquids must focus now on more realistic models while at the same time keeping an accurate methodology when following recent ionic liquids research trends or allowing predictability to come to the foreground. In this Perspective, we summarize in three cases of advanced ionic liquid research what methodological progress has been made and point out difficulties that need to be overcome. As particular examples to discuss we choose reactions, chirality, and radicals in ionic liquids. All these topics have in common that an explicit or accurate treatment of the electronic structure and/or intermolecular interactions is required (accurate methodology), while at the same time system size and complexity as well as simulation time (realistic model) play an important role and must be covered as well.

Two Faces of the Two-Phase Thermodynamic Model

The quantum harmonic model and the two-phase thermodynamic method (2PT) are widely used to obtain quantum-corrected properties such as isobaric heat capacities or molar entropies. 2PT heat capacities were calculated inconsistently in the literature. For water, the classical heat capacity was also considered, but for organic liquids, it was omitted. We reanalyzed the performance of different quantum corrections on the heat capacities of common organic solvents against experimental data. We have pointed out serious flaws in previous 2PT studies. The vibrational density of states was calculated incorrectly causing a 39% relative error in diffusion coefficients and 45% error in the 2PT heat capacities. The wrong conversion of isobaric and isochoric heat capacities also caused about 40% error but in the other direction. We have introduced the concept of anharmonic correction (AC), which is simply the deviation of the classical heat capacity from that of the harmonic oscillator model. This anharmonic contribution is around +30 to 40 J/(mol K) for water depending on the water model and −8 to −10 J/(mol K) for hydrocarbons and halocarbons. AC is unrealistically large, +40 J/(K mol) for alcohols and amines, indicating some deficiency of the OPLS force field. The accuracy of the computations was also assessed with the determination of the self-diffusion coefficients.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N_112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.

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.

Partial Charges Optimized by Genetic Algorithms for Deep Eutectic Solvent Simulations

Deep eutectic solvents (DESs) are a class of solvents often composed of ammonium-based chloride salts and a neutral hydrogen bond donor (HBD) at specific ratios. These cost-effective and environmentally friendly solvents have seen significant growth in multiple fields, including organic synthesis, and in materials and extractions because of their desirable properties. In the present work, a new software called genetic algorithm machine learning (GAML) was developed that utilizes a genetic algorithm (GA) approach to facilitate the development of optimized potentials for liquid simulation (OPLS)-based force field (FF) parameters for eight unique DESs based on three ammonium-based salts and five HBDs at multiple salt:HBD ratios. As an initial test of GAML, partial charges were created for 86 conventional solvents based on neutral organic molecules that yielded excellent overall mean absolute deviations (MADs) of 0.021 g/cm³, 0.63 kcal/mol, and 0.20 kcal/mol compared to experimental densities, heats of vaporization (ΔH_vap), and free energies of hydration (ΔG_hyd), respectively. FFs for DESs constructed from ethylammonium, N,N-diethylethanolammonium, and N-ethyl-N,N-dimethylethanolammonium chloride salts were then parameterized using GAML with exceptional agreement achieved at multiple temperatures for experimental densities, surface tensions, and viscosities with MADs of 0.024 g/cm³, 4.2 mN/m, and 5.3 cP, respectively.

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.

Excess chemical potential of thiophene in [C4MIM] [BF4, Cl, Br, CH3COO] ionic liquids, determined by molecular simulations

The excess chemical potential of thiophene in imidazolium-based ionic liquids [C₄mim][BF₄], [C₄mim][Cl], [C₄mim][Br], and [C₄mim][CH₃COO] were determined by means of molecular dynamics in conjunction with free energy perturbation techniques employing non-polarizable force fields at 300 K and 343.15 K. In addition, energetic and structural analysis were performed such as: interaction energies, averaged noncovalent interactions, radial, and combined distribution functions. The results from this work revealed that the ionic liquids (ILs) presenting the most favorable excess chemical potentials ([C₄mim][BF₄], [C₄mim][CH₃COO]) are associated with the strongest energetic interaction between the thiophene molecule and the ionic liquid anion, and with the weakest energetic interaction between the thiophene molecule and the ionic liquid cation.

Excess chemical potential of thiophene in imidazolium-based ionic liquids [C₄mim][BF₄], [C₄mim][Cl], [C₄mim][Br], and [C₄mim][CH₃COO] determined by molecular simulations.

Dynamics in the BMIM PF6/acetonitrile mixtures observed by femtosecond optical Kerr effect and molecular dynamics simulations

We have performed the measurements of the optical Kerr effect signal time evolution up to 4 ns for a mixture of 1-alkyl-3-methyl-imidazolium hexafluorophosphate (BMIM PF₆) ionic liquid and acetonitrile in the whole mole fractions range. The long delay line in our experimental setup allowed us to capture the complete reorientational dynamics of the ionic liquid. We have analysed the optical Kerr effect signal in the time and frequency domains with help of molecular dynamics simulations. In our approximation of the slow picosecond dynamics with a multi-exponential decay, we distinguish three relaxation times. The highest two are assigned to the reorientation of the cation and acetonitrile molecules that are in the vicinity of the imidazolium ring. The third one is recognized as originating from cation rotations and reorientation of acetonitrile molecules in the bulk or in the vicinity of the aliphatic chains of the cation. With help of the simulation we interpret the intermolecular band in the reduced spectral density, obtained from Kerr signal, as follows: its low-frequency side results from oscillations of one of the components in the cage formed by its neighbors, while the high-frequency side is attributed to the librations of the cation and acetonitrile molecule as well as the intermolecular oscillations of system components involved in specific interactions. We use this assignment and concentration dependence of the spectra obtained from velocity and angular velocity correlations to explain the mole fraction dependence of Kerr reduced spectral density.

Structural investigations of molecular solutes within nanostructured ionic liquids

Ionic liquids (ILs) containing sufficiently long alkyl chains form amphiphilic nanostructures with well-defined polar and non-polar domains. Here we have explored the robustness of these amphiphilic nanostructures to added solutes and gained insight into how the nature of the solute and IL ions affect the partitioning of these solutes within the nanostructured domains of ILs. To achieve this, small angle X-ray scattering (SAXS) investigations were performed and discussed for mixtures of 9 different molecular compounds with 6 different ILs containing imidazolium cations. The amphiphilic nanostructure of ILs persisted to high solute concentrations, over 50 mol% of added solute for most 1-butyl-3-methylimidazolium ILs and above 80 mol% for most 1-decyl-3-methylimidazolium ILs. Solute partitioning within these domains was found to be controlled by the inherent polarity and size of the solute, as well as specific interactions between the solute and IL ions, with SAXS results corroborated with IR spectroscopy and molecular dynamics simulations. Molecular dynamics simulations also revealed the ability to induce π⁺-π⁺ stacking between imidazolium cations through the use of these added molecular compounds. Collectively, these results provide scope for the selection of IL ions to rationally influence and control the partitioning behaviour of given solutes within the amphiphilic nanostructure of ILs.

Accurate Diels-Alder Energies and Endo Selectivity in Ionic Liquids Using the OPLS-VSIL Force Field

Our recently developed optimized potentials for liquid simulations-virtual site ionic liquid (OPLS-VSIL) force field has been shown to provide accurate bulk phase properties and local ion-ion interactions for a wide variety of imidazolium-based ionic liquids. The force field features a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to hydrogen bonding interactions. In this study, the Diels-Alder reaction between cyclopentadiene and methyl acrylate was computationally investigated in the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, [BMIM][PF₆], as a basis for the validation of the OPLS-VSIL to properly reproduce a reaction medium environment. Mixed ab initio quantum mechanics and molecular mechanics (QM/MM) calculations coupled to free energy perturbation and Monte Carlo sampling (FEP/MC) that utilized M06-2X/6-31G(d) and OPLS-VSIL gave activation free energy barriers of 14.9 and 16.0 kcal/mol for the endo-cis and exo-cis Diels-Alder reaction pathways, respectively (exptl. ΔH^‡ of 14.6 kcal/mol). The endo selectivity trend was correctly predicted with a calculated 73% endo preference. The rate and selectivity enhancements present in the endo conformation were found to arise from preferential hydrogen bonding with the exposed C4 ring hydrogen on the BMIM cation. Weaker electronic stabilization of the exo transition state was predicted. For comparison, our earlier ±0.8 charge-scaled OPLS-2009IL force field also yielded a ΔG^‡ of 14.9 kcal/mol for the favorable endo reaction pathway but did not adequately capture the highly organized solvent interactions present between the cation and Diels-Alder transition state.

Thermodynamic, structural and dynamic properties of ionic liquids [C4mim][CF3COO], [C4mim][Br] in the condensed phase, using molecular simulations

In this work a series of thermodynamic, structural, and dynamical properties for the 1-butyl-3-methylimidazolium trifluoroacetate ([C4mim][CF3COO]) and 1-butyl-3-methylimidazolium bromide, ([C4mim][Br]) ionic liquids (ILs) were calculated using Non-polarizable Force Fields (FF), parameterized using a methodology developed previously within the research group, for condensed phase applications. Properties such as the Vapor-Liquid Equilibrium (VLE) curve, critical points (ρ c, T c), Radial, Spatial and Combined Distribution Functions and self-diffusion coefficients were calculated using Equilibrium Molecular Dynamics simulations (EMD); other properties such as shear viscosities and thermal conductivities were calculated using Non-Equilibrium Molecular Dynamics simulations (NEMD). The results obtained in this work indicated that the calculated critical points are comparable with those available in the literature. The calculated structural information for these two ILs indicated that the anions interact mainly with hydrogen atoms from both the imidazolium ring and the methyl chain; the bromide anion displays twice the hydrogen coordination number than the oxygen atoms from the trifluoroacetate anion. Furthermore, Non-Covalent interactions (NCI index), determined by DFT calculations, revealed that some hydrogen bonds in the [C4mim][Br] IL displayed similar strength to those in the [C4mim][CF3COO] IL, in spite of the shorter O--H distances found in the latter IL. The majority of the calculated transport properties presented reasonable agreement with the experimental available data. Nonetheless, the self-diffusion coefficients determined in this work are under-estimated with respect to experimental values; however, by escalating the electrostatic atomic charges for the anion and cation to ±0.8e, only for this property, a remarkable improvement was obtained. Experimental evidence was recovered for most of the calculated properties and to the best of our knowledge, some new predictions were done mainly in thermodynamic states where data are not available. To validate the FF, developed previously within the research group, dynamic properties were also evaluated for a series of ILs such as [C4mim][PF6], [C4mim][BF4], [C4mim][OMs], and [C4mim][NTf2] ILs.