On the Structure of Ionic Liquids: Comparisons between Electronically Polarizable and Nonpolarizable Models I

Ultrafast and Reversible Superwettability Switching of 3D Graphene Foams via Solvent-Exclusive Plasma Treatments

For highly active electron transfer and ion diffusion, controlling the surface wettability of electrically and thermally conductive 3D graphene foams (3D GFs) is required. Here, we present ultrasimple and rapid superwettability switching of 3D GFs in a reversible and reproducible manner, mediated by solvent-exclusive microwave arcs. As the 3D GFs are prepared with vapors of nonpolar acetone or polar water exclusively, short microwave radiation (≤10 s) leads to plasma hotspot-mediated production of methyl and hydroxyl radicals, respectively. Upon immediate radical chemisorption, the 3D surfaces become either superhydrophobic (water contact angle = ∼170°) or superhydrophilic (∼0°), and interestingly, the wettability transition can be repeated many times due to the facile exchange between previously chemisorbed and newly introduced radicals via the formation of methanol-like intermediates. When 3D GFs of different surficial polarities are incorporated into electric double-layer capacitors with nonpolar ionic liquids or polar aqueous electrolytes, the polarity matching between graphene surfaces and electrolytes results in ≥548.0 times higher capacitance compared to its mismatching at ≥0.5 A g-1, demonstrating the significance of wettability-controlled 3D GFs.

Liquid Structure and Hydrogen Bonding in Aqueous Hydroxylammonium Nitrate

Hydroxylammonium nitrate (HAN) has emerged as a promising component in ionic liquid-based spacecraft propellants. However, the physicochemical and structural properties of aqueous HAN have been largely overlooked. The purpose of this study is to investigate the hydrogen bonding in aqueous HAN and understand its implications on these properties and the proton transfer mechanism as a function of concentration. Classical polarizable molecular dynamics simulations have been employed with the APPLE&P force field to analyze the geometry of individual hydrogen bonds and the overall hydrogen-bonding network in various concentrations of aqueous HAN. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) indicate the structural arrangement of the species and their hydrogen bonds. Projections of water density and the orientation of its electric dipole moment near the ions provide insight into the hydrogen-bonding network. The incorporation of water into the hydrogen-bonding network at high ion concentrations occurs via interstitial accommodation around the ions immediately outside the first solvation shell. While ion pairs are observed at all concentrations considered, the frequency of Ha···On hydrogen bonds increases substantially with the ion concentration. The findings contribute to a better fundamental understanding of HAN and the precursors of reactivity, crucial to the development of "green" spacecraft propellants.

GAFF-Based Polarizable Force Field Development and Validation for Ionic Liquids

Ionic liquids (ILs) have been used in many applications, including gas separations, electrochemistry, lubrication, and catalysis. Understanding how the different properties of ILs are related to their chemical structure and composition is crucial for these applications. Experimental investigations often provide limited insights and can be tedious in exploring a range of state points. Therefore, molecular simulations have emerged as a powerful tool that not only offers a microscopic perspective but also enables rapid screening and prediction of physical properties. The accuracy of these predictions, however, depends on the quality of the intermolecular potentials (force fields) used. The widely used classical fixed charge models, such as GAFF, OPLS, and CL&P, are popular due to their simplicity and computational efficiency. However, it has been shown that the use of integer charges with these classical models leads to sluggish dynamics. The use of scaled charge models can improve the dynamics, but these mean-field approaches are unable to account for polarization effects explicitly. Several different approaches have been proposed to include polarizability in IL force fields. In this work, we follow the protocol of the CL&Pol model to develop a Drude oscillator model based on the GAFF force field (Goloviznina, K., et al. J. Chem. Theory Comput. 2019, 15, 5858). We compare the performance of the model for eight imidazolium- and pyrrolidinium-based ILs against that of other models. We find that the new model provides reasonable estimations of density, self-diffusivity, and structural properties for these ILs and suggests a relatively simple way of extending the general GAFF model to more ILs.

Molecular Dynamics Simulation on the Charge Transport Properties in a Salt-in-Ionic Liquid Electrolyte

A clear picture of charge transport properties in salt-in-ionic liquid electrolyte (SILE) is indispensable for the applications in lithium-ion batteries. In this study, we applied molecular dynamics (MD) simulations on a typical SILE system, composed of lithium bis(fluorosulfonyl)imide (LiFSI) with a molar fraction of 0.3 doped in 1-ethyl-3-methylimidazolium bis(fluorosulfonyl)imide (EMIMFSI). Based on the MD simulations, we calculated conductivity spectra from 108 Hz to 1014 Hz, charge current correlation functions, and charge mean square displacements, based on the center-of-mass (COM) velocities of the ions. The conductivity spectra show a bimodal feature between 1012 Hz and 1013 Hz, attributed to the interionic vibrations of the EMIM+-FSI- and Li+-FSI- contact ion pairs, respectively. Structural relaxation is observed between 109 Hz and 1012 Hz, and a flat plateau below 109 Hz, attributed to the direct current (DC) conductivity. For this SILE composed of three constituent ions, i.e., Li+, EMIM+, and FSI-, the above transport properties are further partitioned to the contributions of the individual constituent ions, including self, distinct contribution of the same constituent ions, and also the cross correlation between them. Detailed analyses on the individual contributions reveal strongly correlated motions in this complex ionic system.

Drought timing and severity affect radial growth of Picea crassifolia at different elevations in the western Qilian Mountains

In the context of continued global climate change, the intensity and frequency of droughts have increased to varying degrees in many places. Due to the complexity of drought events, the mechanisms by which trees respond to drought are not well understood. In this study, we analyzed the growth trends of Qinghai spruce (Picea crassifolia) at different elevations in the western part of Qilian Mountains and the dynamic response to climate change. We also compared the differences in radial growth of trees at different elevations in response to drought events in the growing and non-growing seasons based on resistance (Rt), recovery (Rc), and resilience (Rs). The results showed that (1) trees at all three elevations were limited by drought stress and the lower the elevation the more sensitive the trees were to drought. (2) The response of middle- and low-elevation trees to the standardized precipitation evaporation index in June of that year was stable. (3) Growing season drought limits radial growth of trees more than non-growing season drought, and Rt is smaller and Rc is larger at low elevations. With increasing drought severity, trees at all three elevations exhibited a trend of decreasing Rt and Rs and increasing Rc. (4) There were significant differences in the growth trends of trees at the three elevations. Therefore, we should continuously pay attention to the dynamics of the forest ecosystem in the western part of Qilian Mountains and take improved measures to cope with the adverse effects of drought on Qinghai spruce.

A Unified Approach to Derive Atomic Partial Charges and Polarizabilities of Ionic Liquids

We propose a unified approach to fit simultaneously a set of atomic partial charges and polarizabilities of the polarizable model against the ab initio electrostatic potential (ESP) and polarizability. The polarizable model is represented with interactive atomic dipoles with distance-dependent attenuation. For the polarizable model employed in this study, the internal electric field on the polarization sites is fully turned on, and thus allows self-induced dipoles, which persist even for an isolated molecule/ion. By such treatment, the contribution of ESP stems not only from the partial charges but also from the self-induced dipoles, and the atomic partial charges and polarizabilities can be fitted simultaneously against ESP in a unified manner. The fitting with 1-ethyl-3-methylimidazolium (EMIM⁺) and nitrate (NO₃^-), a prototypical organic cation and inorganic anion, respectively, that can form ionic liquid, demonstrates that allowance of the self-induced dipoles gives much better fitness. Moreover, test on the total dipole of an EMIM⁺/NO₃^- ion pair shows that the agreement with the ab initio dipole is also much improved for the polarizable model, which highlights the importance of the polarization effects of ionic liquids.

Long-Range Organization Study of Piperidinium-Based Ionic Liquids by Polarizable Molecular Dynamics Simulations

The nanoscale organization of some classes of ionic liquids is responsible for their singular properties. In this paper, we use polarizable molecular dynamics simulations and small-angle X-ray scattering to probe the structure of two piperidinium- and (trifluoromethylsulfonyl)imide-based ionic liquids ([EBPip⁺][NTf₂^-] and [EOPip⁺][NTf₂^-]) that differ in the alkyl chain length of their cation. The X-ray scattering intensities calculated numerically, from the radial distribution functions, are in excellent agreement with the experimental data. The analysis of the different contributions of the X-ray scattering data allowed us to highlight the correlations responsible for the low q peak observed for the long-chain alkyl cations. New angular analyses showed that anions were more likely to align with alkyl chains as their size increased, inducing angular correlation between anions at larger distances. They also showed that the long alkyl chains of the cations aligned more with each other than the short ones. These more aligned alkyl chains induce a smaller volume of the apolar microdomains compared to the well-studied imidazolium-based ionic liquids, leading to the smaller correlation distance for piperidinium-based ionic liquids.

Hydrophobicity-Dependent Heterogeneous Nanoaggregates and Fluorescence Dynamics in Room-Temperature Ionic Liquids

The hydrophobicity of room-temperature ionic liquids (RTILs) has been shown to have a very significant effect on the optical and structural properties of and in RTILs. The average excited state lifetime of neat RTILs has been shown to be increasing with increasing hydrophobicity of the RTILs. By employing pico-nanosecond-based fluorescence anisotropy decay, the volume of the nanoaggregates in neat RTILs have been calculated. The volume of these nanoaggregates have been shown to be decreasing with increase in hydrophobicity of the RTILs. Thus, hydrophobicity has been shown to have an important role, i.e., hydrophobicity can be used as a handle to tune the properties of RTILs as designer solvents. Moreover, the excited-state lifetime of red-emitting fluorophores, i.e., whose fluorescence emission is not perturbed by the inherent emission of RTILs, has been shown to increase with the increasing hydrophobicity of the RTILs. Highly hydrophobic RTILs have been shown to exhibit positive deviation and highly hydrophilic RTIL has been shown to exhibit negative deviation from the linear correlation between average solvation time (τ_s) versus viscosity/temperature (η/T).

Current Status of AMOEBA-IL: A Multipolar/Polarizable Force Field for Ionic Liquids

Computational simulations of ionic liquid solutions have become a useful tool to investigate various physical, chemical and catalytic properties of systems involving these solvents. Classical molecular dynamics and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations of IL systems have provided significant insights at the atomic level. Here, we present a review of the development and application of the multipolar and polarizable force field AMOEBA for ionic liquid systems, termed AMOEBA–IL. The parametrization approach for AMOEBA–IL relies on the reproduction of total quantum mechanical (QM) intermolecular interaction energies and QM energy decomposition analysis. This approach has been used to develop parameters for imidazolium– and pyrrolidinium–based ILs coupled with various inorganic anions. AMOEBA–IL has been used to investigate and predict the properties of a variety of systems including neat ILs and IL mixtures, water exchange reactions on lanthanide ions in IL mixtures, IL–based liquid–liquid extraction, and effects of ILs on an aniline protection reaction.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Ionic Liquid-Microwave-Based Extraction of Biflavonoids from Selaginella sinensis

Selaginella sinensis (Desv.) Spring has been used for many years as traditional Chinese medicine (TCM) for many years. Recently, ionic liquids (ILs) have attracted great attentions in extraction and separation technology of TCM as a new green solvent. In this paper, microwave assisted extraction-IL (MAE-IL) that extracted amentoflavone (AME) and hinokiflavone (HIN) from Selaginella sinensis was reported for the first time. The contents of two biflavonoids were simultaneously determined by a high performance liquid chromatography (HPLC) method. After different ionic liquids were compared, it was found [C₆mim]BF₄ had a high selectivity and efficiency. Moreover, the important extraction conditions, including solid-liquid ratio, IL concentration, extraction time, microwave power and radiation temperature, were also investigated and optimized by response surface methodology (RSM) using AME and HIN yields as index. The results showed that the extraction yields of AME and HIN from S. sinensis were 1.96 mg/g and 0.79 mg/g, respectively, under the optimal process parameters (0.55 mmol/L, 300 W, 40 min, 1:11 g/mL and 48 °C). Compared with the conventional extraction methods, MAE-IL could not only achieve higher yield in shorter time, but also could reduce the consumption of solvent. This effective, rapid and green MAE-IL method was suitable for the extraction of AME and HIN.

Influence of Electronic Polarization on the Structure of Ionic Liquids

The liquid structure and electrical screening ability of ionic liquids are fundamentally intertwined. The molecular nature of the charge carriers means that screening distances of external fields depend sensitively on the ion packing and structure of the ionic liquid. In this work, we quantitatively illustrate how the liquid structure itself is directly modulated by electrostatic screening conditions. In particular, electronic polarization fundamentally relaxes long-range ion structuring in asymmetric ionic liquids such as [BMIM⁺][BF₄^-], with the influence propagating to short-range ion-ion correlation. A consequence of the exact Stillinger-Lovett second moment condition is that, at fixed density, any pairwise-additive, nonpolarizable force field will necessarily predict artificially enhanced long-range ion structuring. This is because the screening condition is set by the infinite-frequency dielectric response. There is no ad-hoc fix: One has to use polarizable force fields to correctly reproduce the optical dielectric constant. Our illustration of this fundamental effect significantly clarifies interpretation of previous work comparing property prediction using polarizable and nonpolarizable force fields.

Charge-Induced Long-Range Order in a Room-Temperature Ionic Liquid

We report direct evidence for charge-induced long-range (ca. 100 μm) order in the room-temperature ionic liquid (IL) 1-butyl-3-methylimidazolium tetrafluoroborate (BMIM(+)BF4(-)), supported on a silica surface. We have measured the rotational diffusion dynamics of anionic, cationic, and neutral chromophores as a function of distance from a silica surface. The results reflect the excess charge density gradient induced in the IL by the (negative) charge present on the silica surface. Identical measurements in ethylene glycol reveal spatially invariant reorientation dynamics for all chromophores. Capping the silica support with Me2SiCl2 results in spatially invariant reorientation dynamics in the IL. We understand these data in the context of the IL exhibiting a spatially damped piezoelectric response mediated by IL fluidity and disorder.

Ionic Liquid Crystals: Versatile Materials

This Review covers the recent developments (2005-2015) in the design, synthesis, characterization, and application of thermotropic ionic liquid crystals. It was designed to give a comprehensive overview of the "state-of-the-art" in the field. The discussion is focused on low molar mass and dendrimeric thermotropic ionic mesogens, as well as selected metal-containing compounds (metallomesogens), but some references to polymeric and/or lyotropic ionic liquid crystals and particularly to ionic liquids will also be provided. Although zwitterionic and mesoionic mesogens are also treated to some extent, emphasis will be directed toward liquid-crystalline materials consisting of organic cations and organic/inorganic anions that are not covalently bound but interact via electrostatic and other noncovalent interactions.

What type of nanoscopic environment does a cationic fluorophore experience in room temperature ionic liquids?

In the presence of a cationic fluorophore (rhodamine 6G) whose absorption has a significant spectral overlap with the emission of a room temperature ionic liquid (RTIL), the emission of the latter gets quenched, and the quenching has been shown to be dynamic in nature. It has been shown that resonance energy transfer (RET) indeed happens between the RTIL (donor) and rhodamine 6G (cationic acceptor), and RET is the reason for the quenching of the RTIL emission. The spectral and temporal aspects of the RET (between neat RTILs as the donors and rhodamine 6G as the acceptor) were closely studied by steady-state and picosecond time-resolved fluorescence spectroscopy. The influence of the alkyl chain length of the cation, size of the anion, excitation wavelength and concentration of the acceptor on the RET dynamics were also investigated. The energy transfer time (obtained from the rise time of the acceptor) was noted to vary from 2.5 ns to 4.1 ns. By employing the Förster formulation, the donor-acceptor distance was obtained, and its magnitude was found to vary between 31.8 and 37.1 Å. The magnitude of the donor-acceptor distance was shown to be independent of the alkyl chain length of the cation but dependent on the size of the anion of the RTIL. Moreover, the donor-acceptor distance was observed to be independent of the excitation wavelength or concentration of the acceptor. It was shown that the Förster formulation can possibly account for the mechanism and hence can explain the experimental observables in the RET phenomenon. Following the detailed experiments and rigorous analysis, a model has been put forward, which can successfully explain the nanoscopic environment that a cationic fluorophore experiences in an RTIL. Moreover, the nanoscopic environment experienced by the cationic probe has been noted to be different from that experienced by a neutral fluorophore.

On the origin of ionicity in ionic liquids. Ion pairing versus charge transfer

In this paper we show by using static DFT calculations and classical molecular dynamics simulations that the charge transfer between ionic liquid ions plays a major role in the observed discrepancies between the overall mobility of the ions and the observed conductivities of the corresponding ionic liquids, while it also directly suppresses the association of oppositely charged ions, thus the ion pairing. Accordingly, in electrochemical applications of these materials it is important to consider this reduction of the total charges on the ions, which can greatly affect the performance of the given process or device in which the ionic liquid is used. By slightly shifting from the salt-like to a molecular liquid-like system via the decreased charges, the charge transfer also fluidizes the ionic liquid. We believe that this vital information on the molecular level structure of ionic liquids offers a better understanding of these materials, and allows us to improve the a priori design of ionic liquids for any given purpose.

Effects of repulsive interaction on the electric double layer of an imidazolium-based ionic liquid by molecular dynamics simulation

Effects of repulsion between BMIM⁺/PF₆⁻ and the graphite electrode on electric double layer was studied by molecular dynamics simulation.

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.

Comparing induced point-dipoles and Drude oscillators

This work compares the two major methods to introduce polarisability in Molecular Dynamics simulations: induced point-dipoles and Drude oscillators.

Distributed polarizability models for imidazolium-based ionic liquids

Quantum chemical calculations are used to derive distributed polarizability models sufficiently accurate and compact to be used in classical molecular dynamics simulations of imidazolium-based room temperature ionic liquids. Two distributed polarizability models are fitted to reproduce the induction energy of three imidazolium cations (1,3-dimethyl-, 1-ethyl-3-methyl-, and 1-butyl-3-methylimidazolium) and four anions (tetrafluoroborate, hexafluorophosphate, nitrate, and thiocyanate) polarized by a point charge located successively on a grid of surrounding points. The first model includes charge-flow polarizabilities between first-neighbor atoms and isotropic dipolar polarizability on all atoms (except H), while the second model includes anisotropic dipolar polarizabilities on all atoms (except H). For the imidazolium cations, particular attention is given to the transferability of the distributed polarizability sets. The molecular polarizability and its anisotropy rebuilt by the distributed models are found to be in good agreement with the exact ab initio values for the three cations and 23 additional conformers of 1-ethyl-3-methyl-, 1-butyl-3-methyl-, 1-pentyl-3-methyl-, and 1-hexyl-3-methylimidazolium cations.

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.

Stochastic structure determination for conformationally flexible heterogenous molecular clusters: application to ionic liquids

We present a novel method that enables accurate and efficient computational determination of conformationally flexible clusters, "Kick(3)" This method uses stochastically generated structures in combination with fast quantum mechanical methods. We demonstrate the power of this method by elucidating the structure of ionic liquid (IL) ([xMIM(+)][NO3(-)])n clusters (x = E, B, D, n = 1-10,15). Dispersion-corrected, third-order self-consistent-charge density-functional tight-binding (DFTB3) is shown to be a computationally efficient, yet reliable approximation to density functional theory for predicting and understanding IL structure and stability. The presented approach, therefore, enables the accurate and efficient screening of ILs with high potential toward practical applications, without recourse to more expensive quantum chemical methods.

Electrostatics interactions in classical simulations

Electrostatic interactions are crucial for both the accuracy and performance of atomistic biomolecular simulations. In this chapter we review well-established methods and current developments aiming at efficiency and accuracy. Specifically, we review the classical Ewald summations, particle-particle particle-method particle-method Ewald algorithms, multigrid, fast multipole, and local methods. We also highlight some recent developments targeting more accurate, yet classical, representation of the molecular charge distribution.