Influence of Charge Scaling on the Solvation Properties of Ionic Liquid Solutions

Effectively enhancing ion diffusion in superconcentrated ionic liquid electrolytes using co-solvent additives

The incorporation of high salt concentrations in ionic liquid (IL) electrolytes, forming superconcentrated ionic liquids, has been shown to improve Li-ion transference numbers and enhance cycling stability against lithium metal anodes. However, this benefit comes at the cost of significantly increased viscosity and reduced ionic conductivity due to the formation of large ion aggregates. To optimize conductivity further, a co-solvent can be introduced at an optimal concentration to enhance ion transport while preserving superior interfacial stability. The effectiveness of this approach depends on the solvent as it affects ion diffusion to varying degrees. This computational study examines how co-solvents can effectively enhance metal ion diffusion in superconcentrated ionic liquids by comparing two widely used organic solvents. We found that the key lies in their ability to effectively participate in Li solvation shells, disrupting the large Li-anion aggregates. Our results show that anion exchange in a Li(anion)x(solvent)y hybrid solvation shell occurs more rapidly than in a Li (anion)z solvation shell, facilitating Li diffusion through a structural diffusion mechanism. A co-solvent with a high donor number exhibits a stronger affinity for lithium ions, which is identified as a crucial factor in enhancing ion diffusion. This work provides valuable insights to guide the design of superconcentrated ionic liquid electrolytes for lithium-metal battery development.

Lessons Learned on Obtaining Reliable Dynamic Properties for Ionic Liquids

Ionic liquids are nowadays investigated with respect to their use as electrolytes for high-performance energy storage materials. In this study, we provide a tutorial on how to calculate dynamic properties such as self-diffusion coefficients, ionic conductivities, transference numbers, as well as ion pair and ion cage dynamics, that all play a role in judging the applicability of ionic liquids as electrolytes. For the case of the ionic liquid[ C 2 C 1 Im ] [ NTf 2 ] ${[{\rm{C}}_2 {\rm{C}}_1 {\rm{Im}}][{\rm{NTf}}_2 ]}$ , we investigate the performance of different force fields. Amongst them are non-polarizable models employing unity charges, a charge-scaled version of a non-polarizable model, a polarizable model and another non-polarizable model with refined Lennard-Jones parameters. We also study the influence of the system size on the dynamic properties. While all studied force field models capture qualitatively correct trends, only the polarizable force field and the non-polarizable force field with refined Lennard-Jones parameters provide quantitative agreement to reference data, making the latter model very attractive for the reason of lower computational costs.

Screening and Design of Aqueous Zinc Battery Electrolytes Based on the Multimodal Optimization of Molecular Simulation

Aqueous batteries, such as aqueous zinc-ion batteries (AZIB), have garnered significant attention because of their advantages in intrinsic safety, low cost, and eco-friendliness. However, aqueous electrolytes tend to freeze at low temperatures, which limits their potential industrial applications. Thus, one of the core challenges in aqueous electrolyte design is optimizing the formula to prevent freezing while maintaining good ion conductivity. However, the experimental trial-and-error approach is inefficient for this purpose, and existing simulation tools are either inaccurate or too expensive for high-throughput phase transition predictions. In this work, we employ a small amount of experimental data and differentiable simulation techniques to develop a multimodal optimization workflow. With minimal human intervention, this workflow significantly enhances the prediction power of classical force fields for electrical conductivity. Most importantly, the simulated electrical conductivity can serve as an effective predictor of electrolyte freezing at low temperatures. Generally, the workflow developed in this work introduces a new paradigm for electrolyte design. This paradigm leverages both easily measurable experimental data and fast simulation techniques to predict properties that are challenging to access by using either approach alone.

Predicting Ionic Conductivity of Imidazolium-Based Ionic Liquid Mixtures Using Quantum-Mechanically Derived Partial Charges in the Condensed Phase

A considerable effort has been expended over the years to tune the properties of ionic liquids (ILs) by designing cations, anions, and pendant groups on the ions. A simple and effective approach to altering the properties of ILs is formulating IL-IL mixtures. However, the measurements and properties of such mixtures lag considerably behind those of pure ILs. From a molecular simulation point of view, binary IL mixtures have been investigated using charge distributions of pure ILs, which implicitly assumes that the ions of different polarizability do not influence the local electronic environment due to changing concentrations. To understand this effect, molecular dynamics (MD) simulations were conducted for a series of IL-IL mixtures containing the common cation 1-ethyl-3-methylimidazolium [C2mim] varying the composition of various combinations of anions (tetrafluoroborate [BF4] and dicyanamide [DCA], [BF4] and bis(trifluoromethanesulfonyl)imide [NTF2], [BF4] and trifluoromethanesulfonate [TFO], and [TFO] and [NTF2]). The effect of changing the electronic environment was evaluated by deriving partial charges using density functional theory (DFT) calculations in the condensed phase. It was observed that the overall charge on the cation and anion was a function of the cation-anion pairings for pure ILs. Moreover, the cation charge was found to vary linearly with anionic concentrations. Improved agreement of predicted density and ionic conductivity with experimental values was found for binary IL mixtures with this approach, in comparison to that when a fixed charge model is employed.

Evaluation of Density-Functional Tight-Binding Methods for Simulation of Protic Molecular Ion Pairs

In this work, we benchmark the accuracy of the density-functional tight-binding (DFTB) method, namely the long-range corrected second-order (LC-DFTB2) and third-order (DFTB3) models, for predicting energetics of imidazolium-based ionic liquid (IL) ion pairs. We compare the DFTB models against popular density functionals such as LC-ωPBE and B3LYP, using ab initio domain-based local pair-natural orbital coupled cluster (DLPNO-CC) energies as reference. Calculations were carried out in the gas phase, as well as in aqueous solution using implicit solvent methods. We find that the LC-DFTB2 model shows excellent performance in the gas phase and agrees well with reference energies in implicit solvent, often outperforming DFTB3 predictions for complexation energetics. Our study identifies a range of opportunities for use of the LC-DFTB method and quantifies its sensitivity to protonation states and the types of chemical interactions between ion pairs.

Molecular dynamics investigation of IEPOX chemical behavior at the interface and in the bulk phase of acidic aerosols

Isoprene epoxydiol (IEPOX) is an important reactive gas-phase intermediate produced by the photooxidation of isoprene under low NOx conditions, playing a key role in the formation of secondary organic aerosols (SOA). Previous studies have mostly focused on the liquid-phase reactions of IEPOX within aerosols; however, interfacial heterogeneous chemical reactions are equally important in SOA formation. This study systematically explores the reaction mechanisms of IEPOX at the acidic aerosol interface and in the bulk phase using classical molecular dynamics (MD) and ab initio molecular dynamics simulations (AIMD). The study found that the free energy of IEPOX at the aerosol interface significantly decreases, indicating that interfacial heterogeneous chemical reactions are indispensable for the formation of IEPOX-derived SOA. The research reveals the formation pathways of 2-methyltetrols (2-MTO) and 1,3,4-trihydroxy-3-methylbutan-2-yl sulfates (2-MTOOS), finding that the protonation of the epoxy O atom and the cleavage of the C-O bond are the rate-controlling steps, while the nucleophilic addition is a spontaneous process. Through multiple sets of simulations, it was observed that the formation frequency of 2-MTO at the acidic aerosol interface and in the bulk phase reached 53.8%, significantly higher than the 30.8% of 2-MTOOS, which is consistent with field observation data. Additionally, through metadynamics (MTD) simulations, it was suggested that IEPOX could undergoes acid-catalyzed ring-opening reactions at the interface, potentially followed by the transfer of H atoms from primary alcohols into the aerosol, leading to the possible formation of the intermediate product 3-methylbut-3-ene-1,2,4-triol (one of the proposed structures of C5-alkene triols). These findings provide new insights into the formation mechanism of IEPOX-derived SOA and offer a scientific basis for future studies on their physicochemical properties and atmospheric fate.

A Polarizable Forcefields for Glyoxal Acetals as Electrolyte Components for Lithium-Ion Batteries

In this work we have derived the parameters of an AMOEBA-like polarizable forcefield for electrolytes based on tetramethoxy and tetraethoxy-glyoxal acetals, and propylene carbonate. The resulting forcefield has been validated using both ab-initio data and the experimental properties of the fluids. Using molecular dynamics simulations, we have investigated the structural features and the solvation properties of both the neat liquids and of the corresponding 1 M LiTFSI electrolytes at the molecular level. We present a detailed analysis of the Li ion solvation shells, of their structure and highlight the different behavior of the solvents in terms of their molecular structure and coordinating features.

Temperature-Dependent Dielectric Relaxation Measurements of (Betaine + Urea + Water) Deep Eutectic Solvent in Hz-GHz Frequency Window: Microscopic Insights into Constituent Contributions and Relaxation Mechanisms

A combined experimental and simulation study of dielectric relaxation (DR) of a deep eutectic solvent (DES) composed of betaine, urea, and water with the composition [Betaine:Urea:Water = 11.7:12:1 (weight ratio) and 9:18:5 (molar ratio)] was performed to explore and understand the interaction and dynamics of this system. Temperature-dependent (303 ≤ T/K ≤ 343) measurements were performed over 9 decades of frequency, combining three different measurement setups. Measured DR, comprising four distinct steps with relaxation times spreading over a few picoseconds to several nanoseconds, was found to agree well with simulations. The simulated total DR spectra, upon dissection into three self (intraspecies) and three cross (interspecies) interaction contributions, revealed that the betaine-betaine self-term dominated (∼65%) the relaxation, while the urea-urea and water-water interactions contributed only ∼7% and ∼1%, respectively. The cross-terms (betaine-urea, betaine-water, and urea-water) together accounted for <30% of the total DR. The slowest DR component with a time constant of ∼1-10 ns derived dominant contribution from betaine-betaine interactions, where betaine-water and urea-water interactions also contributed. The subnanosecond (0.1-0.6 ns) time scale originated from all interactions except betaine-water interaction. An extensive interaction of water with betaine and urea severely reduced the average number of water-water H-bonds (∼0.7) and heavily decreased the static dielectric constant of water in this DES (εs ∼ 2). Furthermore, simulated first rank collective single particle reorientational relaxations (C1(t)) and the structural H-bond fluctuation dynamics (CHB (t)) exhibited multiexponential kinetics with time scales that corresponded well with those found both in the simulated and measured DR.

Anion-Dependent Strength Scale of Interactions in Ionic Liquids from X-ray Photoelectron Spectroscopy, Ab Initio Molecular Dynamics, and Density Functional Theory

Using a combination of experiments and calculations, we have gained new insights into the nature of anion-cation interactions in ionic liquids (ILs). An X-ray photoelectron spectroscopy (XPS)-derived anion-dependent electrostatic interaction strength scale, determined using XPS core-level binding energies for IL cations, is presented here for 39 different anions, with at least 18 new anions included. Linear correlations of experimental XPS core-level binding energies for IL cations with (a) calculated core binding energies (ab initio molecular dynamics (AIMD) simulations were used to generate high-quality model IL structures followed by single-point density functional theory (DFT) to obtain calculated core binding energies), (b) experimental XPS core-level binding energies for IL anions, and (c) other anion-dependent interaction strength scales led to three main conclusions. First, the effect of different anions on the cation can be related to ground-state interactions. Second, the variations of anion-dependent interactions with the identity of the anion are best rationalized in terms of electrostatic interactions and not occupied valence state/unoccupied valence state interactions or polarizability-driven interactions. Therefore, the XPS-derived anion-dependent interaction strength scale can be explained using a simple electrostatic model based on electrostatic site potentials. Third, anion-probe interactions, irrespective of the identity of the probe, are primarily electrostatic, meaning that our electrostatic interaction strength scale captures some inherent, intrinsic property of anions independent of the probe used to measure the interaction strength scale.

Adsorption Isotherm and Mechanism of Ca2+ Binding to Polyelectrolyte

Polyelectrolytes, such as poly(acrylic acid) (PAA), can effectively mitigate CaCO3 scale formation. Despite their success as antiscalants, the underlying mechanism of binding of Ca2+ to polyelectrolyte chains remains unresolved. Through all-atom molecular dynamics simulations, we constructed an adsorption isotherm of Ca2+ binding to sodium polyacrylate (NaPAA) and investigated the associated binding mechanism. We find that the number of calcium ions adsorbed [Ca2+]ads to the polymer saturates at moderately high concentrations of free calcium ions [Ca2+]aq in the solution. This saturation value is intricately connected with the binding modes accessible to Ca2+ ions when they bind to the polyelectrolyte chain. We identify two dominant binding modes: the first involves binding to at most two carboxylate oxygens on a polyacrylate chain, and the second, termed the high binding mode, involves binding to four or more carboxylate oxygens. As the concentration of free calcium ions [Ca2+]aq increases from low to moderate levels, the polyelectrolyte chain undergoes a conformational transition from an extended coil to a hairpin-like structure, enhancing the accessibility to the high binding mode. At moderate concentrations of [Ca2+]aq, the high binding mode accounts for at least one-third of all binding events. The chain's conformational change and its consequent access to the high binding mode are found to increase the overall Ca2+ ion binding capacity of the polyelectrolyte chain.

Nonaqueous Ion Pairing Exemplifies the Case for Including Electronic Polarization in Molecular Dynamics Simulations

The inclusion of electronic polarization is of crucial importance in molecular simulations of systems containing charged moieties. When neglected, as often done in force field simulations, charge-charge interactions in solution may become severely overestimated, leading to unrealistically strong bindings of ions to biomolecules. The electronic continuum correction introduces electronic polarization in a mean-field way via scaling of charges by the reciprocal of the square root of the high-frequency dielectric constant of the solvent environment. Here, we use ab initio molecular dynamics simulations to quantify the effect of electronic polarization on pairs of like-charged ions in a model nonaqueous environment where electronic polarization is the only dielectric response. Our findings confirm the conceptual validity of this approach, underlining its applicability to complex aqueous biomolecular systems. Simultaneously, the results presented here justify the potential employment of weaker charge scaling factors in force field development.

Molecular Modeling of Water-in-Salt Electrolytes: A Comprehensive Analysis of Polarization Effects and Force Field Parameters in Molecular Dynamics Simulations

Accurate modeling of highly concentrated aqueous solutions, such as water-in-salt (WiS) electrolytes in battery applications, requires proper consideration of polarization contributions to atomic interactions. Within the force field molecular dynamics (MD) simulations, the atomic polarization can be accounted for at various levels. Nonpolarizable force fields implicitly account for polarization effects by incorporating them into their van der Waals interaction parameters. They can additionally mimic electron polarization within a mean-field approximation through ionic charge scaling. Alternatively, explicit polarization description methods, such as the Drude oscillator model, can be selectively applied to either a subset of polarizable atoms or all polarizable atoms to enhance simulation accuracy. The trade-off between simulation accuracy and computational efficiency highlights the importance of determining an optimal level of accounting for atomic polarization. In this study, we analyze different approaches to include polarization effects in MD simulations of WiS electrolytes, with an example of a Na-OTF solution. These approaches range from a nonpolarizable to a fully polarizable force field. After careful examination of computational costs, simulation stability, and feasibility of controlling the electrolyte properties, we identify an efficient combination of force fields: the Drude polarizable force field for salt ions and non-polarizable models for water. This cost-effective combination is sufficiently flexible to reproduce a broad range of electrolyte properties, while ensuring simulation stability over a relatively wide range of force field parameters. Furthermore, we conduct a thorough evaluation of the influence of various force field parameters on both the simulation results and technical requirements, with the aim of establishing a general framework for force field optimization and facilitating parametrization of similar systems.

Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

Toward Bottom-Up Understanding of Transport in Concentrated Battery Electrolytes

Bottom-up understanding of transport describes how molecular changes alter species concentrations and electrolyte voltage drops in operating batteries. Such an understanding is essential to predictively design electrolytes for desired transport behavior. We herein advocate building a structure-property-performance relationship as a systematic approach to accurate bottom-up understanding. To ensure generalization across salt concentrations as well as different electrolyte types and cell configurations, the property-performance relation must be described using Newman's concentrated solution theory. It uses Stefan-Maxwell diffusivity, ij , to describe the role of molecular motions at the continuum scale. The key challenge is to connect ij to the structure. We discuss existing methods for making such a connection, their peculiarities, and future directions to advance our understanding of electrolyte transport.

Nonpolarizable Force Fields through the Self-Consistent Modeling Scheme with MD and DFT Methods: From Ionic Liquids to Self-Assembled Ionic Liquid Crystals

A key to achieve the accuracy of molecular dynamics (MD) simulation is the set of force fields used to express the atomistic interactions. In particular, the electrostatic interaction remains the main issue for the precise simulation of various ionic soft materials from ionic liquids to their supramolecular compounds. In this study, we test the nonpolarizable force fields of ionic liquids (ILs) and self-assembled ionic liquid crystals (ILCs) for which the intermolecular charge transfer and intramolecular polarization are significant. The self-consistent modeling scheme is adopted to refine the atomic charges of ionic species in a condensed state through the use of density functional theory (DFT) under the periodic boundary condition. The atomic charges of the generalized amber force field (GAFF) are effectively updated to express the electrostatic properties of ionic molecules obtained by the DFT calculation in condensed phase, which improves the prediction accuracy of ionic conductivity with the obtained force field (GAFF-DFT). The derived DFT charges then suggest that the substitution of a hydrophobic liquid-crystalline moiety into IL-based cations enhances the charge localization of ionic groups in the amphiphilic molecules, leading to the amplification of the electrostatic interactions among the hydrophilic/ionic groups in the presence of hydrophobic moieties. In addition, we focus on an ion-conductive pathway hidden in the self-assembled nanostructure. The MD results indicate that the ionic groups of cation and anion interact strongly for keeping the bicontinuous nanosegregation of ionic nanochannel. The partial fractions of hydrophilic/ionic and hydrophobic nanodomains are then quantified with the volume difference from referenced IL systems, while the calculated ionic conductivity decreases in the self-assembled ILCs more than the occupied volume of ionic nanodomains. These analyses suggest that the mobility of ions in the self-assembled ILCs remains quite restricted even with small tetrafluoroborate anions because of strong attractive interaction among ionic moieties.

Recent Developments in Polarizable Molecular Dynamics Simulations of Electrolyte Solutions

Polarizable molecular dynamics simulations are a fast progressing field in the scientific research of ionic liquids. The fundamentals of polarizable simulations, as well as their application to ionic liquids, were summarized in a review [Bedrov, D.; Piquemal, J.-P.; Borodin, O.; MacKerell, Jr., A. D.; Roux, B.; Schröder, C. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields. Chem. Rev. 2019, 119, 7940–7995] in 2019. Since then, new methods to treat intermolecular interaction of induced dipoles in these highly charged systems were developed. This concerns the damping of these interactions and additional charge transfer as well as the prediction of ionic materials with ultrahigh refractive indices. In addition to the progress of the polarizable force fields, also thermostats and barostats for polarizable simulations evolved recently.

Ion Correlation in Choline Chloride-Urea Deep Eutectic Solvent (Reline) from Polarizable Molecular Dynamics Simulations

In recent years, deep eutectic solvents (DESs) emerged as highly tunable and environmentally friendly alternatives to common ionic liquids and organic solvents. In this work, a polarizable model based on the CHARMM Drude polarizable force field is developed for a 1:2 ratio mixture of choline chloride/urea (reline) DES. To successfully reproduce the structure of the liquid as compared to first-principles molecular dynamics simulations, a damping factor was introduced to correct the observed over-binding between the chloride and the hydrogen bonding site of choline. Investigated radial distributions reveal the formation of hydrogen bonds between all the constituents of reline and similar interactions for chloride and urea's oxygen atoms, which could contribute to the melting point depression of the mixture. Predicted dynamic properties from our polarizable force field are in good agreement with experiments, showing significant improvements over nonpolarizable models. Similar to some ionic liquids, an oscillatory behavior in the velocity autocorrelation function of the anion is visible, which can be interpreted as a rattling motion of the lighter anion surrounded by the heavier cations. The obtained results for ionic conductivity of reline show some degree of correlated ion motion in this DES. However, a joint diffusion of ion pairs cannot be observed during the simulations.

Electronic Polarization Is Essential for the Stabilization and Dynamics of Buried Ion Pairs in Staphylococcal Nuclease Mutants

Buried charged residues play important roles in the modulation of protein stabilities and conformational dynamics and make crucial contributions to protein functions. Considering the generally nonpolar nature of protein interior, a key question concerns the contribution of electronic polarization to the stabilization and properties of buried charges. We answer this question by conducting free energy simulations using the latest polarizable CHARMM force field based on Drude oscillators for a series of Staphylococcal nuclease mutants that involve a buried Glu-Lys pair in different titration states and orientations. While a nonpolarizable model suggests that the ionized form of the buried Glu-Lys pair is more than 40 kcal/mol less stable than the charge-neutral form, the two titration states are comparable in stability when electronic polarization is included explicitly, a result better reconcilable with available experimental data. Analysis of free energy components suggests that additional stabilization of the ionized Glu-Lys pair has contributions from both the enhanced salt-bridge strength and stronger interaction between the ion-pair and surrounding protein residues and penetrated water. Despite the stronger direct interaction between Glu and Lys, the ion-pair exhibits considerably larger and faster structural fluctuations when polarization is included, due to compensation of interactions in the cavity. Collectively, observations from this work provide compelling evidence that electronic polarization is essential to the stability, hydration, dynamics, and therefore function of buried charges in proteins. Therefore, our study advocates for the explicit consideration of electronic polarization for mechanistic and engineering studies that implicate buried charged residues, such as enzymes and ion transporters.

Charge transfer and polarisability in ionic liquids: a case study

The practical use of ionic liquids (ILs) is benefiting from a growing understanding of the underpinning structural and dynamic properties, facilitated through classical molecular dynamics (MD) simulations. The predictive and explanatory power of a classical MD simulation is inextricably linked to the underlying force field. A key aspect of the forcefield for ILs is the ability to recover charge based interactions. Our focus in this paper is on the description and recovery of charge transfer and polarisability effects, demonstrated through MD simulations of the widely used 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C₄C₁im][NTf₂] IL. We study the charge distributions generated by a range of ab initio methods, and present an interpolation method for determining atom-wise scaled partial charges. Two novel methods for determining the mean field (total) charge transfer from anion to cation are presented. The impact of using different charge models and different partial charge scaling (unscaled, uniformly scaled, atom-wise scaled) are compared to fully polarisable simulations. We study a range of Drude particle explicitly polarisable potentials and shed light on the performance of current approaches to counter known problems. It is demonstrated that small changes in the charge description and MD methodology can have a significant impact; biasing some properties, while leaving others unaffected within the structural and dynamic domains.

Interfacial Properties of Hydrophobic Deep Eutectic Solvents with Water

Hydrophobic deep eutectic solvents (DESs) have recently gained much attention as water-immiscible solvents for a wide range of applications. However, very few studies exist in which the hydrophobicity of these DESs is quantified. In this work, the interfacial properties of hydrophobic DESs with water were computed at various temperatures using molecular dynamics simulations. The considered DESs were tetrabutylammonium chloride–decanoic acid (TBAC–dec) with a molar ratio of 1:2, thymol–decanoic acid (Thy–dec) with a molar ratio of 1:2, and dl-menthol–decanoic acid (Men–dec) with a molar ratio of 2:1. The following properties were investigated in detail: interfacial tensions, water-in-DES solubilities (and salt-in-water solubilities for TBAC–dec/water), density profiles, and the number densities of hydrogen bonds. Different ionic charge scaling factors were used for TBAC–dec. Thy–dec and Men–dec showed a high level of hydrophobicity with negligible computed water-in-DES solubilities. For charge scaling factors of 0.7 and 1 for the thymol and decanoic acid components of Thy–dec, the computed interfacial tensions of the DESs are in the following order: TBAC–dec (ca. 4 mN m^–1) < Thy–dec (20 mN m^–1) < Men–dec (26 mN m^–1). The two sets of charge scaling factors for Thy–dec did not lead to different density profiles but resulted in considerable differences in the DES/water interfacial tensions due to different numbers of decanoic acid–water hydrogen bonds at the interfaces. Large peaks were observed for the density profiles of (the hydroxyl oxygen of) decanoic acid at the interfaces of all DES/water mixtures, indicating a preferential alignment of the oxygen atoms of decanoic acid toward the aqueous phase.

The effect of hydration on the stability of ionic liquid crystals: MD simulations of [C14C1im]Cl and [C14C1im]Cl·H2O

The thermal range of the stability of Ionic Liquid Crystal (ILC) phases of imidazolium ILCs, and the type of the mesophase itself are affected by several molecular structural features, the two prominent ones being the alkyl chain length and the counter-anion. Hydration is also very important: monohydrate samples of 1-alkyl-3-methylimidazolium halides have a higher clearing point and a wider thermal range of the stability of the ionic smectic phase, compared with the analogous anhydrous sample. To understand the reasons, at a microscopic level, for such increased stability due to hydration, we run classical Molecular Dynamics (MD) simulations of a typical ionic liquid crystal, 1-tetradecyl-3-methylimidazolium chloride, and of its monohydrate form. We tested a full-charge non-polarizable force field and a scaled-charge version having the total charge of the ions scaled by a factor of 0.80. Comparison of the structural and dynamic properties with available experimental data reveals that the scaling of the charge by a factor of 0.80 results in a good agreement between simulated and experimental data and it sheds light on the microscopic mechanism responsible for the increased stability of the monohydrated phase. A hydrogen-bond network between water and the chloride anion is established in the ionic layer which increases the stability of the ionic layer; this in turn increases the nano-segregation between the ionic and hydrophobic layers which eventually produce an increased order of the alkylic layer as well.

Deep Eutectic Solvents: Molecular Simulations with a First-Principles Polarizable Force Field

The unique properties of deep eutectic solvents make them useful in a variety of applications. In this work we develop a first-principles force field for reline, which is composed of choline chloride and urea in the molar ratio 1:2. We start with the symmetry adapted perturbation theory (SAPT) protocol and then make adjustments to better reproduce the structure and dynamics of the liquid when compared to first-principles molecular dynamics (FPMD) simulations. The resulting force field is in good agreement with experiments in addition to being consistent with the FPMD simulations. The simulations show that primitive molecular clusters are preferentially formed with choline-chloride ionic pairs bound with a hydrogen bond in the hydroxyl group and that urea molecules coordinate the chloride mainly via the trans-H chelating hydrogen bonds. Incorporating polarizability qualitatively influences the radial distributions and lifetimes of hydrogen bonds and affects long-range structural order and dynamics. The polarizable force field predicts a diffusion constant about an order of magnitude larger than the nonpolarizable force field and is therefore less computationally intensive. We hope this study paves the way for studying complex hydrogen-bonding liquids from a first-principles approach.

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.

Why Lithium Ions Stick to Some Anions and Not Others

Designing battery electrolytes for lithium-ion batteries has been a topic of extensive research for decades. The ideal electrolyte must have a large conductivity as well as high Li⁺ transference number. The conductivity is very sensitive to the nature of the anions and dynamical correlations between ions. For example, lithium bis(trifluoromethane)sulfonimide (LiTFSI) has a large conductivity, but the chemically similar lithium trifluoromethanesulfonate (LiOTf) shows poor conductivity. In this work, we study the binding of Li⁺ to these anions in an ethylene carbonate (EC) solvent using enhanced sampling metadynamics. The evaluated free energies display a large dissociation barrier for LiOTf compared to LiTFSI, suggesting long-lived ion-pair formation in the former but not the latter. We probe these observations via unbiased molecular dynamics simulations and metadynamics simulations of TFSI with a hypothetical OTF-like partial charge model indicating an electrostatic origin for those differences. Our results highlight the deleterious impact of sulfonate groups in lithium-ion battery electrolytes and provide a new basis for the assessment of electrolyte designs.

Quantum Simulations of Hydrogen Bonding Effects in Glycerol Carbonate Electrolyte Solutions

The need for environmentally friendly nonaqueous solvents in electrochemistry and other fields has motivated recent research into the molecular-level solvation structure, thermodynamics, and dynamics of candidate organic liquids. In this paper, we present the results of quantum density functional theory simulations of glycerol carbonate (GC), a molecule that has been proposed as a solvent for green industrial chemistry, nonaqueous alternatives for biocatalytic reactions, and liquid media in energy storage devices. We investigate the structure and dynamics of both the pure GC liquid and electrolyte solutions containing KF and KCl ion pairs. These simulations reveal the importance of hydrogen bonding that controls the structural and dynamic behavior of the pure liquid and ion association in the electrolyte solutions. The results illustrate the difficulties associated with classical modeling of complex organic solvents. The simulations lead to a better understanding of the underlying mechanisms behind the previously observed peculiar ion-specific behavior in GC electrolyte solutions.

Realistic Ion Dynamics through Charge Renormalization in Nonaqueous Electrolytes

While many practically important electrolytes contain lithium ions, interactions of these ions are particularly difficult to probe experimentally because of their small X-ray and neutron scattering cross sections and large neutron absorption cross sections. Molecular dynamics (MD) is a powerful tool for understanding the properties of nonaqueous electrolyte solutions from the atomic level, but the accuracy of this computational method crucially depends on the physics built into the classical force field. Here, we demonstrate that several force fields for lithium bistriflimide (LiTFSI) in acetonitrile yield a solution structure that is consistent with the neutron scattering experiments, yet these models produce dramatically different ion dynamics in solution. Such glaring discrepancies indicate that inadequate representation of long-range interactions leads to excessive ionic association and ion-pair clustering. We show that reasonable agreement with the experimental observations can be achieved by renormalization of the ion charges using a "titration" method suggested herewith. This simple modification produces realistic concentration dependencies for ionic diffusion and conductivity in <2 M solutions, without loss in quality for simulation of the structure.

Molecular dynamics simulation of imidazolium CnMIM-BF4 ionic liquids using a coarse grained force-field

Ionic liquids feature thermophysical properties that are of interest in solvents, energy storage materials and tunable lubrication applications. Here we use new Coarse Grained (CG) models to investigate the structure, dynamics and interfacial properties of the [C_2-8MIM][BF₄] family of ionic liquids (ILs). The simulated equation of state and diffusion coefficients are in good agreement with experimental data and with all-atom force-fields. We quantify the nano-structure and liquid-vapour interfacial properties of the ILs as a function of the size of the imidazolium cation. The computational efficiency of the CG models enables the simulation of very long time scales (100's of nanoseconds), which are needed to resolve the dynamic and interfacial properties of ILs containing cations with long aliphatic chains. For [C_>4MIM] [BF₄] the break in symmetry associated to the liquid-vapour interface induces nanostructuring of polar and non-polar domains in the direction perpendicular to the interface plane, with the inhomogeneous regions penetrating deep inside the bulk liquid, typically 5 nm for C₈MIM cations.

Self-Consistent Scheme Combining MD and Order-N DFT Methods: An Improved Set of Nonpolarizable Force Fields for Ionic Liquids

The nonpolarizable force field of ionic liquids is tuned by using the self-consistent scheme of molecular dynamics (MD) simulation and first-principles calculation based on the order-N density functional theory (DFT). The atomic charges are determined by using the whole MD cell for DFT calculation and accounts effectively for the many-body effects of charge transfer and intramolecular polarization. The charges represent effective interactions in the condensed phase within the framework of the nonpolarizable force field and can be an alternative for an explicitly many-body model incorporating, for example, polarizability. Here we demonstrate the performance of nonpolarizable force field determined with the MD-DFT self-consistent scheme in imidazolium-, pyrrolidinium-, and ammonium-based ionic liquids. The variation ranges of molecular charges are much larger with the compositions of the ionic liquid than with the thermodynamic conditions, and the charge-ordering structures become systematically weaker with the effective charges. For energetic properties, while the calculated heat of vaporization depends on the atomic and molecular charges, the corresponding heat capacity is not strongly affected by the DFT-based variation. For transport properties, the self-diffusion coefficient, electrical conductivity, and viscosity vary much more in the self-consistent scheme. The effective DFT charge is observed to enhance the fluidity of ionic liquids and improve the accuracy of electrical conductivity and viscosity. This is due to the weakened interactions among the ions, and the too slow motions observed with a full-charge model are well corrected through the iteration of MD and DFT. We therefore conclude that the set of nonpolarizable force fields obtained with the MD-DFT self-consistent scheme leads to better description of transport properties of ionic liquids.