Pairing of 1-hexyl-3-methylimidazolium and tetrafluoroborate ions in n-pentanol

Polarizability Plays a Decisive Role in Modulating Association between Molecular Cations and Anions

Electrostatic interactions involving proteins depend on not only the ionic charges involved but also their chemical identities. Here we examine the origins of incompletely understood differences in the strength of association of different pairs of monovalent molecular ions that are relevant to protein-protein and protein-ligand interactions. Cationic analogues of the basic amino acid side chains are simulated, along with oxyanionic analogues of cation-exchange ligands and acidic amino acids. Experimentally observed association trends with respect to the cations, but not anions, are captured by a nonpolarizable model. An effective continuum correction to account for electronic polarizability can capture both trends better but at the expense of fidelity to the underlying free energy landscape for ion-pair association. A polarizable model proves decisive in capturing experimentally suggested trends with respect to both cations and anions; critically, the free energy landscape for ion-pair association is itself altered, thus altering configurational sampling.

Specific Ion Solvation and Pairing Effects in Glycerol Carbonate

Identifying the driving forces behind the solvation of inorganic salts by nonaqueous solvents is an important step in the development of green solvents. Here we focus on one promising solvent: glycerol carbonate (GC). Using ab initio molecular dynamics simulations, we build upon our previous work by detailing glycerol carbonate's interactions with a series of anions, a lithium ion, and the LiF ion pair. Through these investigations, we highlight the changes in solvation behavior as the anion size increases, the competition of binding shown by lithium for the oxygens of GC, and the behavior of the LiF ion pair in a GC solution. These results indicate the importance of the cation's identity in ion-pairing structure and dynamics and lend insight into the key factors behind the specific ion effects seen in GC.

Comparison of single-ion molecular dynamics in common solvents

Laying a basis for molecularly specific theory for the mobilities of ions in solutions of practical interest, we report a broad survey of velocity autocorrelation functions (VACFs) of Li⁺ and PF₆^- ions in water, ethylene carbonate, propylene carbonate, and acetonitrile solutions. We extract the memory function, γ(t), which characterizes the random forces governing the mobilities of ions. We provide comparisons controlling for the effects of electrolyte concentration and ion-pairing, van der Waals attractive interactions, and solvent molecular characteristics. For the heavier ion (PF₆^-), velocity relaxations are all similar: negative tail relaxations for the VACF and a clear second relaxation for γt, observed previously also for other molecular ions and with n-pentanol as the solvent. For the light Li⁺ ion, short time-scale oscillatory behavior masks simple, longer time-scale relaxation of γt. But the corresponding analysis of the solventberg Li⁺H₂O₄ does conform to the standard picture set by all the PF₆^- results.

Assessment of Simple Models for Molecular Simulation of Ethylene Carbonate and Propylene Carbonate as Solvents for Electrolyte Solutions

Progress in understanding liquid ethylene carbonate (EC) and propylene carbonate (PC) on the basis of molecular simulation, emphasizing simple models of interatomic forces, is reviewed. Results on the bulk liquids are examined from the perspective of anticipated applications to materials for electrical energy storage devices. Preliminary results on electrochemical double-layer capacitors based on carbon nanotube forests and on model solid-electrolyte interphase (SEI) layers of lithium ion batteries are considered as examples. The basic results discussed suggest that an empirically parameterized, non-polarizable force field can reproduce experimental structural, thermodynamic, and dielectric properties of EC and PC liquids with acceptable accuracy. More sophisticated force fields might include molecular polarizability and Buckingham-model description of inter-atomic overlap repulsions as extensions to Lennard-Jones models of van der Waals interactions. Simple approaches should be similarly successful also for applications to organic molecular ions in EC/PC solutions, but the important case of Li\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ deserves special attention because of the particularly strong interactions of that small ion with neighboring solvent molecules. To treat the Li\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ ions in liquid EC/PC solutions, we identify interaction models defined by empirically scaled partial charges for ion-solvent interactions. The empirical adjustments use more basic inputs, electronic structure calculations and ab initio molecular dynamics simulations, and also experimental results on Li\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ thermodynamics and transport in EC/PC solutions. Application of such models to the mechanism of Li\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ transport in glassy SEI models emphasizes the advantage of long time-scale molecular dynamics studies of these non-equilibrium materials.

Liquid/liquid interface layering of 1-butanol and [bmim]PF6 ionic liquid: a nonlinear vibrational spectroscopy and molecular dynamics simulation study

IV-SFG vibrational spectroscopy and MD simulation studies reveal a local polar/nonpolar layering structure at the interface of 1-butanol-d₉ and 1-butyl-3-methylimidazolium hexafluorophosphate ([bmim]PF₆).

State of hydrophobic and hydrophilic ionic liquids in aqueous solutions: are the ions fully dissociated?

Molecular dynamics simulations were performed for aqueous solutions of five ionic liquids (ILs): 1-ethyl-3-methylimidazolium ([C2mim]) bis(trifluoromethanesulfonyl) imide ([NTf2]), 1-n-butyl-3-methylimidazolium ([C4mim]) [NTf2], 1-n-hexyl-3-methylimidazolium ([C6mim]) [NTf2], [C2mim] ethylsulfate ([C2H5SO4]), and [C2mim] chloride (Cl) to determine whether the ions of these ILs are associated at relatively high dilutions and whether the association is governed by hydrophobicity/hydrophilicity of the ILs. The adaptive biasing force technique was applied to calculate the potential of mean force (PMF) for each IL ion pair. For all of the ILs, the PMF is characterized by two distinct contact minima in which the ions have different relative conformations. The hydrophobic ILs bearing the anion [NTf2](-) exist predominantly in the associative state; the strength of the association of these ILs increases with increase in the alkyl chain length. The most hydrophilic IL [C2mim] Cl was determined to be almost fully dissociated at the concentration examined in the study. [C2mim] [C2H5SO4] showed hydration behavior that was intermediate between that exhibited by the ILs in which the anion is substituted with either Cl(-) or [NTf2](-) paired with [C2mim](+). Association constants for these ILs were also computed. Radial distribution functions calculated by constraining the ions at the contact minima showed that hydration of the anion plays the dominant role in determining the microscopic behavior of these ILs in aqueous solutions.