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

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

We develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

Horn-like Pore Entrance Boosts Charging Dynamics and Charge Storage of Nanoporous Supercapacitors

Optimizing the synergy between nanoporous carbons and ionic liquids can significantly enhance the energy density of supercapacitors. The highest energy density has been obtained as the size of porous carbon matches the size of ionic liquids, while it may result in slower charging dynamics and thus reduce the power density. Enhancing energy storage without retarding charging dynamics remains challenging. Herein, we designed porous electrodes by introducing an optimized horn-like entrance to the nanopore, which can concurrently improve supercapacitors' charging dynamics and energy storage. Our results revealed the mechanism of improved charging lies in the gradual desolvation process and optimized ion motion paths: the former expedites the adsorption of the counterion by reducing the transitional energy barrier for ions entering the pores, and the latter accelerates the co-ion desorption and eliminates ion overfilling. Meanwhile, the enhancement of energy density could be attributed to the multi-ion coordinated migration.

A molecular dynamics study of a fully zwitterionic copolymer/ionic liquid-based electrolyte: Li+ transport mechanisms and ionic interactions

The development of polymer electrolytes (PEs) is crucial for advancing safe, high-energy density batteries, such as lithium-metal and other beyond lithium-ion chemistries. However, reaching the optimum balance between mechanical stiffness and ionic conductivity is not a straightforward task. Zwitterionic (ZI) gel electrolytes comprising lithium salt and ionic liquid (IL) solutions within a fully ZI polymer network can, in this context, provide useful properties. Although such materials have shown compatibility with lithium metal in batteries, several fundamental structure-dynamic relationships regarding ionic transport and the Li⁺ coordination environment remain unclear. To better resolve such issues, molecular dynamics simulations were carried out for two IL-based electrolyte systems, N-butyl-N-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide ([BMP][TFSI]) with 1 M LiTFSI salt and a ZI gel electrolyte containing the IL and a ZI copolymer: poly(2-methacryloyloxyethyl phosphorylcholine-co-sulfobetaine vinylimidazole), poly(MPC-co-SBVI). The addition of ZI polymer decreases the [TFSI]^- -[Li]⁺ interactions and increases the IL ion diffusivities, and consequently, the overall ZI gel ionic conductivity. The structural analyses showed a large preference for lithium-ion interactions with the polymer phosphonate groups, while the [TFSI]^- anions interact directly with the sulfonate group and the [BMP]⁺ cations only display secondary interactions with the polymer. In contrast to previous experimental data on the same system, the simulated transference numbers showed smaller [Li]⁺ contributions to the overall ionic conductivities, mainly due to negatively charged lithium aggregates and the strong lithium-ion interactions in the systems.

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.

Spatial-Decomposition Analysis of Electrical Conductivity in Mixtures of Ionic Liquid and Sodium Salt for Sodium-Ion Battery Electrolytes

Ionic liquid (IL)-based electrolytes are a promising material for the development of sodium-ion batteries, and their performance can be quantified by electrical conductivity. In this highly concentrated ionic system, the correlated motions of ion pairs are influential on the ionic transport properties. Herein, all-atom analyses are conducted through molecular dynamics simulations to bridge the macroscopically observable electrical conductivity with the molecular pictures of correlated motion of ion pairs. The analysis is applied to three mixtures of IL with sodium salt that are relevant to the electrolyte for a sodium-ion battery: [1-ethyl-3-methylimidazolium, Na][bis(fluorosulfonyl)amide] ([C₂C₁im, Na][FSA]), [N-methyl-N-propylpyrrolidinium, Na][FSA] ([C₃C₁pyrr, Na][FSA]), and [K, Na][FSA]. The computational results on electrical conductivities are in agreement with the experimental reports, and their dependency on temperature and sodium-ion composition is reproduced well. The overall contributions from cross-correlated motions are found to be negative in all the IL mixtures; thus, the total conductivities are less than their Nernst-Einstein estimates. The spatial view of cross-correlated motions is further obtained by decomposing the time correlation functions of velocities according to the distances between ion pairs. It is observed that ion pairs are moving in the same direction for ∼0.3 ps when they were initially within the first coordination shell, followed by motions toward opposite directions. The cross-correlation terms are also dissected into local and nonlocal components, and it is commonly seen for all the ion pairs that the local component is negative for cation-anion pairs and is positive for cation-cation and anion-anion pairs. The motions of ion pairs are accompanied by a "backflow" that manifests in the form of the nonlocal component whose sign is opposite to the corresponding local component. In fact, the contributions of the correlated motions of ions to the electrical conductivity are not localized to contact pairs and extend spatially beyond the first coordination shell of the cation-anion pairs.

Microscopic origins of conductivity in molten salts unraveled by computer simulations

Molten salts are crucial materials in energy applications, such as batteries, thermal energy storage systems or concentrated solar power plants. Still, the determination and interpretation of basic physico-chemical properties like ionic conductivity, mobilities and transference numbers cause debate. Here, we explore a method for determination of ionic electrical mobilities based on non-equilibrium computer simulations. Partial conductivities are then determined as a function of system composition and temperature from simulations of molten LiF_αCl_βI_γ (with α + β + γ = 1). High conductivity does not necessarily coincide with high Li⁺ mobility for molten LiF_αCl_βI_γ systems at a given temperature. In salt mixtures, the lighter anions on average drift along with Li⁺ towards the negative electrode when applying an electric field and only the heavier anions move towards the positive electrode. In conclusion, the microscopic origin of conductivity in molten salts is unraveled here based on accurate ionic electrical mobilities and an analysis of the local structure and kinetics of the materials.

Insights into the Stabilization of Fluoride Ions in Ionic Liquids: Pointers to Better Fluorinating Agents

The fluorination efficiency of a fluorinating agent depends on the free availability of the fluoride ions, which in turn depends on its interaction with its solvation shell. A stable fluoride-based poor solvate ionic liquid (SIL) comprising 1-ethyl-3-methylimidazolium (EMIM) cation and ethylene glycol (EG) was recently reported and demonstrated as a fluorinating agent. Herein, we performed ab initio calculations and ab initio molecular dynamics simulations to gain a microscopic understanding of the intermolecular interactions in this SIL in gas, liquid, and crystalline phases. Ethylene glycol (EG), being capable of forming hydrogen bond(s) with the fluoride ion, prevents the latter from reacting with the EMIM cation. Fluoride forms hydrogen bonds with both the cation and the EG molecule, but it was found to have more affinity toward EG, forming a stronger hydrogen bond with its hydroxyl proton than with the acidic proton of the cation. An optimal concentration of EG in the SIL balances its contribution to stabilizing the fluoride ion and yet making fluoride available for fluorination.

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.

Aqueous Nanoclusters Govern Ion Partitioning in Dense Polymer Membranes

The uptake and sorption of charged molecules by responsive polymer membranes and hydrogels in aqueous solutions is of key importance for the development of soft functional materials. Here, we investigate the partitioning of simple monatomic (Na⁺, K⁺, Cs⁺, Cl^-, I^-) and one molecular ion (4-nitrophenolate; NP^-) within a dense, electroneutral poly(N-isopropylacrylamide) membrane using explicit-water molecular dynamics simulations. Inside the predominantly hydrophobic environment, water distributes in a network of polydisperse water nanoclusters. The average cluster size determines the mean electrostatic self-energy of the simple ions, which preferably reside deeply inside them; we therefore find substantially larger partition ratios K ≃10^-1 than expected from a simple Born picture using a uniform dielectric constant. Despite their irregular shapes, we observe that the water clusters possess a universal negative electrostatic potential with respect to their surroundings, as is known for aqueous liquid-vapor interfaces. This potential, which we find concealed in cases of symmetric monatomic salts, can dramatically impact the transfer free energies of larger charged molecules because of their weak hydration and increased affinity to interfaces. Consequently, and in stark contrast to the simple ions, the molecular ion NP^- can have a partition ratio much larger than unity, K ≃10-30 (depending on the cation type) or even 10³ in excess of monovalent salt, which explains recent observations of enhanced reaction kinetics of NP^- reduction catalyzed within dense polymer networks. These results also suggest that ionizing a molecule can even enhance the partitioning in a collapsed, rather hydrophobic gel, which strongly challenges the traditional simplistic reasoning.

Atomistic Insights into the Layered Microstructure and Time-Dependent Stability of [BMIM][PF6] Confined within the Meso-Slit of Carbon

Clarifying the microstructures and time-dependent stability of ionic liquids (ILs) within the confinement of the meso-slit of carbon is the first step to understand the intrinsic synergy effect between ILs and a promising mesoporous carbon electrode. In this work, we adopted molecular dynamics to systematically investigate the behavior of [BMIM][PF₆] in the 2.8 nm-wide meso-slit of carbon. The confined ILs formed a pronounced layered spatial distribution and can be divided into three distinct regions, namely, com-, sub-, and cen-layer, according to valley coordinates in the number density profiles. In the com-layer region, the imidazolium rings of ILs possess two dominant orientations, namely, "parallel" and "tilted standing". The rotation ability of all the ions is highly restrained. In the sub-layer and cen-layer regions, a part of the [BMIM] imidazolium ring has a preferred "tilted standing" orientation. The [BMIM] cations are still in a rotational restrain state and show a preferred rotation motion along the x-y plane. The hydrogen bond between [BMIM] cations and [PF₆] anions play a crucial role in determining the confined multilayered spatial distribution and distinctive orientation properties of ILs.

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.

Charge Environment and Hydrogen Bond Dynamics in Binary Ionic Liquid Mixtures: A Computational Study

The extent of charge transfer between the cation and the anion in a room-temperature ionic liquid depends on the basicity of the anion. Ion charges determined in the condensed state via density functional theory calculations capture this effect rather well, and charges derived in such a manner have been employed in force field-based molecular dynamics simulations to quantitatively reproduce several physical properties of the liquids. However, the issue of transferability of cation charges in mixtures of ionic liquids, say with one type of cation and two different anion types needs to be addressed. Herein, we demonstrate that the cation charge in such a mixture varies linearly with anion composition, a result that ties in rather well with X-ray photoelectron spectroscopic experiments. The variation in cation charge with bulk anion composition is shown to be a result of changes in its coordination environment. Cations surrounded by a higher proportion of more basic anions possess lower charges than those surrounded by less basic anions. Time scales for the exchange of anion types for the occupation of hydrogen bonding sites around the cation have been determined and are seen to be constituted by three processes-breakage of existing hydrogen bond, diffusion to the hydrogen bonding site and displacement of the incumbent anion from its site in the cation coordination shell. These time scales explain the differences observed between infrared and NMR spectroscopic experiments in ionic liquid mixtures rather well.

Preferential solvation and ion association properties in aqueous dimethyl sulfoxide solutions

We study the solvation and the association properties of ion pairs in aqueous dimethyl sulfoxide (DMSO) solution by atomistic molecular dynamics (MD) simulations. The ion pair is composed of two lithium and a single sulfonated diphenyl sulfone ion whose properties are studied under the influence of different DMSO concentrations. For increasing mole fractions of DMSO, we observe a non-ideal behavior of the solution as indicated by the derivatives of the chemical activity. Our findings are complemented by dielectric spectra, which also verify a complex DMSO-water mixing behavior. In agreement with these results, further simulation outcomes reveal an aqueous homoselective solvation of the ion species which fosters the occurrence of pronounced ion association constants at higher DMSO mole fractions. The consequences of this finding are demonstrated by lower ionic conductivities for increasing concentrations of DMSO.

Local environment structure and dynamics of CO2 in the 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide and related ionic liquids

Despite the innumerous papers regarding the study of the ionic liquids as a potential candidate for CO₂ capture, many details concerning the structure and dynamics of CO₂ in the system are still to be revealed, i.e., the correlation between the local environment structure and the dynamic properties of the substance. This present work relied on the performance of molecular dynamics both for the neat [C₂mim][Tf₂N] and [C₂mim][Tf₂N]/CO₂ mixtures in an attempt to elucidate the local environment of CO₂ and their effects on the dynamic properties of [C₂mim][Tf₂N]. A slight change in the orientation of the cation and anion could be observed, which was correlated to the cation and anion moving away from each other in order to receive the carbon dioxide. The gas molecules pushed both the cation and the anion away to create sufficient void to its accommodation. The diffusion coefficient of [C₂mim]⁺ is higher than [Tf₂N]^- regardless the increase of the CO₂ concentration. The addition of CO₂ in the ionic liquid has shown an increase of 4-5 times for the diffusivity of ions, which was related to the decrease of cation-anion interaction strength. The transport properties' results showed that the addition of CO₂ in the ionic liquid generates the fluidization of the system, decreasing the viscosity as a consequence of the local environment structure changing. Likewise, the effect of the type of anion and cation on the system properties was studied considering [Ac]^- and [BMpyr]⁺ ions, showing large effects by the change of anion to [Ac]^- which rise from the strong [C₂mim]⁺-[Ac]^- interaction, which conditions the solvation of ions by CO₂ molecules.

Significance of weak interactions in imidazolium picrate ionic liquids: spectroscopic and theoretical studies for molecular level understanding

The effects of interionic hydrogen bonding and π-π stacking interactions on the physical properties of a new series of picrate anion based ionic liquids (ILs) have been investigated experimentally and theoretically. The existence of aromatic (C2-HO) and aliphatic (C7-HO-N22 and C6-HO-N20) hydrogen bonding and π-π stacking interactions in these ILs has been observed using various spectroscopic techniques. The aromatic and aliphatic C-HO hydrogen bonding interactions are found to have a crucial role in binding the imidazolium cation and picrate anion together. However, the π-π stacking interactions between two successive layers are found to play a decisive role in tight packing in ILs leading to differences in physical properties. The drastic difference in the melting points of the methyl and propyl derivatives (mmimPic and pmimPic respectively) have been found to be primarily due to the difference in the strength and varieties of π-π stacking interactions. While in mmimPic, several different types of π-π stacking interactions between the aromatic rings (such as picrate-picrate, picrate-imidazole and imidazolium-imidazolium cation rings) are observed, only one type of π-π stacking interaction (picrate-picrate rings) is found to exist in the pmimPic IL. NMR spectroscopic studies reveal that the interaction of these ILs with solvent molecules is different and depends on the dielectric constant of the solvent. While an ion solvation model explains the solvation in high dielectric solvents, an ion-pair solvation model is found to be more appropriate for low dielectric constant solvents. The enhanced stability of these investigated picrate ILs compared with that of inorganic picrate salts under high doses of γ radiation clearly indicates the importance of weak interionic interactions in ILs, and also opens up the possibility of the application of picrate ILs as prospective diluents in nuclear separation for advanced fuel cycling process.

Computational approaches to understanding reaction outcomes of organic processes in ionic liquids

The utility of using a combined experimental and computational approach for understanding ionic liquid media, and their effect on reaction outcome, is highlighted through a number of case studies.

Solvent effects of 1-ethyl-3-methylimidazolium acetate: solvation and dynamic behavior of polar and apolar solutes

We study the solvation properties of the ionic liquid 1-ethyl-3-methylimidazolium acetate ([eMIM]⁺[ACE]⁻) and the resulting dynamic behavior for differently charged model solutes at room temperature via atomistic molecular dynamics (MD) simulations of 500 ns length.

Quantitative prediction of physical properties of imidazolium based room temperature ionic liquids through determination of condensed phase site charges: a refined force field

Quantitative prediction of physical properties of room temperature ionic liquids through nonpolarizable force field based molecular dynamics simulations is a challenging task. The challenge lies in the fact that mean ion charges in the condensed phase can be less than unity due to polarization and charge transfer effects whose magnitude cannot be fully captured through quantum chemical calculations conducted in the gas phase. The present work employed the density-derived electrostatic and chemical (DDEC/c3) charge partitioning method to calculate site charges of ions using electronic charge densities obtained from periodic density functional theory (DFT) calculations of their crystalline phases. The total ion charges obtained thus range between -0.6e for chloride and -0.8e for the PF6 ion. The mean value of the ion charges obtained from DFT calculations of an ionic liquid closely matches that obtained from the corresponding crystal thus confirming the suitability of using crystal site charges in simulations of liquids. These partial charges were deployed within the well-established force field developed by Lopes et al., and consequently, parameters of its nonbonded and torsional interactions were refined to ensure that they reproduced quantum potential energy scans for ion pairs in the gas phase. The refined force field was employed in simulations of seven ionic liquids with six different anions. Nearly quantitative agreement with experimental measurements was obtained for the density, surface tension, enthalpy of vaporization, and ion diffusion coefficients.

Probing molecular interaction in ionic liquids by low frequency spectroscopy: Coulomb energy, hydrogen bonding and dispersion forces

Low vibrational spectroscopy provides detailed information on the strength and type of interaction and their influence on the properties of ionic liquids.

Ionic liquids studied across different scales: a computational perspective

For theoreticians, ionic liquids represent a major challenge. This is due to the fact that intermolecular interactions are particularly strong because of ionic liquids' ionicity. This, in turn, causes a subtle interplay between different scales which is encoded in the measured macro- and mesoscopic properties and also in the molecular electrostatic characteristics. Therefore, force fields have to describe the microscopic processes correctly in order to reproduce macroscopic properties accurately over a large range of state variables. Herein, imidazolium-based ionic liquids were studied at different scales, going from the detailed quantum electronic scale to the classical atomistic scale. It is indicated how the information gained at each level could be used for the other scales. In particular, the issue of deriving suitable partial charges for use in classical force fields is addressed. The Blöchl method was employed to generate partial charges reproducing the multipole distribution accurately for bulk systems. This led naturally to absolute ionic charges of less than /l e/, i.e., charge scaling. So, the monopole structure of the herein introduced force field mimics the quantum chemical behaviour observed in the liquid phase. This led to a substantial improvement in the description of dynamical properties of immediate experimental interest, such as electric conductivity. For further insight, the electric dipole moment of the ions was taken as physical indicator of their electronic structure. The electric dipole moment was found to fluctuate strongly and to depend on polarisation. Hence, our scale-combined study offers a gateway to rational design of models, based on the relevant underlying physics rather than on mere numerical parameterisation, and thereby to (possibly) more direct physical interpretation of experimental results.