A Refined All-Atom Potential for Imidazolium-Based Room Temperature Ionic Liquids: Acetate, Dicyanamide, and Thiocyanate Anions

Modeling High Concentration Bisalt-in-Sulfolane Electrolytes and the Observation of Ligand-Bridged Cation-Pair Complexes

Despite the abundance of sodium over lithium in Earth's crust and the copious amounts of expensive lithium salt required to make Li-ion high-concentration electrolytes (HCEs), studies of HCEs made from sodium salts remain sparse. A comparative molecular-level study of Li- and Na-ion HCEs and mixed cation or bisalt HCEs in an organic solvent is missing. To fill this gap, we studied model HCEs of pure and mixed Li and Na salts of bis(fluorosulfonyl)amide (FSI) in sulfolane using a confluence of classical molecular dynamics (MD), ab initio MD (AIMD) simulations, and quantum chemical cluster calculations. While Li-ion HCEs display transport properties superior to those of Na-ion HCEs, the latter's performance can be considerably improved by replacing even 25% of Na-ions with Li-ions. While the effects of doping are largely systemic, a larger sensitivity of the identity of solvation shells of Li-ions to the Li-content of the HCE is observed; in contrast, those of Na-ions are more oblivious to it. Fascinating ligand-bridged, short-distance cation pairs observed in the classical MD simulations are confirmed using density functional theory-based AIMD simulations. Quantum chemical calculations in the gas phase reveal the thermodynamic stability of such cation pairs complexed with multiple anions and solvent molecules.

Solvation Structures and Transport Mechanisms of Cations in Water-in-Bisalt Electrolytes for High-Concentration Lithium-Ion Batteries Revealed by Molecular Dynamics Simulations

The development of safe and cost-effective electrolytes for rechargeable batteries is currently underway. While water-based electrolytes hold promise, their restricted electrochemical stability window poses a challenge. Combining multiple ionic species emerges as a promising strategy to broaden this stability window and optimize Li-ion battery performance. This study focuses on dual-cation electrolytes, which blend lithium and potassium acetates to enhance the electrochemical characteristics of the solution at high concentrations. We investigated the solvation structure of each ion and its interactions on a molecular level. Our analysis reveals that ion clusters and aggregates are formed through shared acetate and water molecules at high salt concentrations. Furthermore, the residence time analyses of atom pairs indicate that cations diffuse in vehicular mode at low concentrations. In contrast, they switch to a structural mode at high concentrations due to diminishing water content. This study offers a comprehensive model for exploring diverse solvation structures of cations and gaining insights into their diffusion mechanisms within water-in-bisalt electrolytes for aqueous Li-ion batteries.

Insight to the Local Structure of Mixtures of Imidazolium-Based Ionic Liquids and Molecular Solvents from Molecular Dynamics Simulations and Voronoi Analysis

While the physicochemical properties as well as the NMR and vibration spectroscopic data of the mixtures of ionic liquids (ILs) with molecular solvents undergo a drastic change around the IL mole fraction of 0.2, the local structure of the mixtures pertaining to this behavior remains unclear. In this work, the local structure of 12 mixtures of 1-butyl-3-methylimidazolium cation (C₄mim⁺) combined with perfluorinated anions, such as tetrafluoroborate (BF₄^-), hexafluorophosphate (PF₆^-), trifluoromethylsulfonate (TFO^-), and bis(trifluoromethanesulfonyl)imide, (TFSI^-), and aprotic dipolar solvents, such as acetonitrile (AN), propylene carbonate (PC), and gamma butyrolactone (γ-BL) is studied by molecular dynamics simulations in the entire composition range, with an emphasis on the IL mole fractions around 0.2. Distributions of metric properties corresponding to the Voronoi polyhedra of the particles (volume assigned to the particles, local density, radius of spherical voids) are determined, using representative sites of the cations, anions, and the solvent molecules, to characterize the changes in the local structure of these mixtures. By analyzing the mole fraction dependence of the average value, fluctuation, and skewness parameter of these distributions, the present study reveals that, around the IL mole fraction of 0.2, the local structure of the mixture undergoes a transition between that determined by the interionic interactions and that determined by the interactions between the ions and solvent molecules. It should be noted that the strength of the interactions between the ions and the solvent molecules, modulated by the change in the composition of the mixture, plays an important role in the occurrence of this transition. The signature of the change in the local structure is traced back to the nonlinear change of the mean values, fluctuations, and skewness values of the metric Voronoi polyhedra distributions.

Non-monotonic composition dependence of the breakdown of Stokes-Einstein relation for water in aqueous solutions of ethanol and 1-propanol: explanation using translational jump-diffusion approach

A series of experimental and simulation studies examined the validity of the Stokes-Einstein relationship (SER) of water in binary water/alcohol mixtures of different mixture compositions. These studies revealed a strong non-monotonic composition dependence of the SER with maxima at the specific alcohol mole fraction where the non-idealities of the thermodynamic and transport properties are observed. The translational jump-diffusion (TJD) approach elucidated the breakdown of the SER in pure supercooled water as caused by the jump translation of molecules. The breakdown of SER in the supercooled water/methanol binary mixture was successfully explained using the same TJD approach. To further generalize the picture, here we focus on the non-monotonic composition dependence of SER breakdown of water in two water/alcohol mixtures (water/ethanol and water/propanol) for a broad temperature range. In agreement with previous studies, maximum breakdown of SER is observed for the mixture with alcohol mole fraction x = 0.2. Diffusion of the water molecules at the maximum SER breakdown point is largely contributed by jump-diffusion. The residual-diffusion, obtained by subtracting the jump-diffusion from the total diffusion, approximately follows the SER for different compositions and temperatures. We also performed hydrogen (H-)bond dynamics and observed that the contribution of jump-diffusion is proportional to the total free energy of activation of breaking all H-bonds that exist around a molecule. This study, therefore, suggests that the more a molecule is trapped by H-bonding, the more likely it is to diffuse through the jump-diffusion mechanism, eventually leading to an increasing degree of SER breakdown.

Structures and properties of ionic crystals and condensed phase ionic liquids predicted with the generalized energy-based fragmentation method

The generalized energy-based fragmentation (GEBF) approach is extended to facilitate ab initio investigations of structures, lattice energies, vibrational spectra and ¹ H NMR chemical shifts of ionic crystals and condensed-phase ionic liquids (ILs) with the periodic boundary conditions (PBC). For selected periodic systems, our results demonstrate that the so-called PBC-GEBF approach can provide satisfactory descriptions on ground-state energies, structures, and vibrational spectra of ionic crystals and IL crystals. The PBC-GEBF approach is then applied to three realistic condensed phase systems. For three ionic crystals (LiCl, NaCl, and KCl), we apply the PBC-GEBF approach with MP2 theory as well as some popular DFT methods to investigate their crystal structures and lattice energies. Our calculations indicate that the crystal structures obtained with PBC-GEBF-MP2/6-311 + G** are very close to the corresponding X-ray structures, while PBC-GEBF-ωB97X-D/6-311 + G** provides satisfactory prediction for crystal structures and lattice energies. For two polymorphs of [n-C₄ mim][Cl] crystals, we find that the PBC-GEBF approach at the M06-2X/6-311 + G** level can give a satisfactory descriptions on structures and Raman spectra of these two crystals. Furthermore, for [C₂ mim][BF₄ ] ILs, we demonstrate that their ¹ H NMR chemical shifts can be estimated from averaging over 5 typical snapshots (extracted from MD simulations) with the PBC-GEBF approach at the B97-2/pcSseg-2 level. The calculated results account for the observed experimental data quite well. Therefore, we expect that the PBC-GEBF approach, combined with various quantum chemistry methods, will become an effective tool in predicting structures and properties of ionic crystals and condensed-phase ILs.

Xenon Diffusion in Ionic Liquids with Blurred Nanodomain Separation

The dynamics of xenon gas, loaded in a series of 1-alkyl-3-methylimidazolium based ionic liquids, probes the formation of increasingly blurred polar/apolar nanodomains as a function of the anion type and the cation chain length. Exploiting ¹²⁹ Xe NMR spectroscopy techniques, like Pulse Gradient Spin Echo (PGSE) and inversion recovery (IR), the diffusion motion and relaxation times are determined for 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [C_n C₁ im][TFSI]. A correlation between the ILs nano-structure and both xenon diffusivity and relaxation times, as well as chemical shifts, is outlined. Interestingly, comparison with previous results of the same properties in the homologous imidazolium chlorides and hexafluorophospate shows an opposite trend with the alkyl chain length. Classical molecular dynamics (MD) simulations are used to calculate the xenon and cation and anion diffusion coefficients in the same systems, including imidazolium cations with longer chains (n=4, 6, 8 … 20). An almost quantitative agreement with the experiments validates the MD simulations and, at the same time, provides the necessary structural and dynamic microscopic insights on the nano-segregation and diffusion of xenon in bistriflimide, chloride and hexafluorphosphate salts allowing to observe and rationalize the shaping effect of the cation in the nanostructure.

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.

Molecular simulation of osmometry in aqueous solutions of the BMIMCl ionic liquid: a potential route to force field parameterization of liquid mixtures

Despite widespread development and use of ionic liquids (ILs) in both academic and industrial research, computational force fields (FFs) for most of those are not available for a precise description of inter-species interactions in aqueous environments. In the scope of this study, by means of molecular simulations, the osmotic coefficient of an aqueous solution of an IL is calculated and used as a basis to reparameterize popular IL-FFs existing in the literature. We first calculate the osmotic coefficients (at 298.15 K and 1 atm pressure) of aqueous solutions of 1-butyl-3-methylimidazolium chloride (BMIMCl), a generic IL, popularly used in biomass processing and the subsequent conversion to value-added intermediates. The performance of two popular atomic, nonpolarizable FFs developed for BMIMCl, one by Lopes, Pádua, and coworkers (FF-LP) and the other by Sambasivarao, Acevedo, and coworkers (FF-SA), when mixed with the SPC/E water model, is tested with respect to their ability to reproduce the experimental osmotic coefficient data. Interestingly, the osmotic coefficient is found to be increasing with a gradual increase in IL molality within the concentration range of our investigation, which is contrary to the experimental trend reported in the literature for the same IL-water mixture. Henceforth, necessary corrections to the nonbonded ion-ion and ion-water interactions are made to match the experimental osmotic coefficient. To further assess the reliability of the new FF, we extensively explore the thermodynamic (density, isothermal compressibility, and thermal expansion coefficient), dynamic (diffusivity and viscosity), and association/dissociation properties (rationalized with the help of radial distribution functions) with both the original and reparameterized FF for a wider range of concentrations up to a molality of 18.50 mol kg-1. The calculated quantities are compared against experimental data wherever available. The modified FF parameters exhibit significant improvements in terms of its ability to match experimental solution properties, such as density, viscosity, association/dissociation, etc. We report that excessive dissociation of BMIMCl in water is responsible for the shortcomings observed in the original FFs and improved prediction of physicochemical properties could be achieved using the modified FFs.

Voronoi Polyhedra as a Tool for the Characterization of Inhomogeneous Distribution in 1-Butyl-3-methylimidazolium Cation-Based Ionic Liquids

The inhomogeneity distribution in four imidazolium-based ionic liquids (ILs) containing the 1-butyl-3-methylimidazolium (C₄mim) cation, coupled with tetrafluoroborate (BF₄), hexafluorophosphate (PF₆), bis(trifluoromethanesulfonyl)amide (TFSA), and trifluoromethanesulfonate (TfO) anions, was characterized using Voronoi polyhedra. For this purpose, molecular dynamic simulations have been performed on the isothermal-isobaric (NpT) ensemble. We checked the ability of the potential models to reproduce the experimental density, heat of vaporization, and transport properties (diffusion and viscosity) of these ionic liquids. The inhomogeneity distribution of ions around the ring, methyl, and butyl chain terminal hydrogen atoms of the C₄mim cation was investigated by means of Voronoi polyhedra analysis. For this purpose, the position of the C₄mim cation was described successively by the ring, methyl, and butyl chain terminal hydrogen atoms, while that of the anions was described by their F or O atom. We calculated the Voronoi polyhedra distributions of the volume, the density, and the asphericity parameters as well as that of the radius of the spherical intermolecular voids. We carried out the analysis in two steps. In the first step, both ions were taken into account. The calculated distributions gave information on the neighboring ions around a reference one. In the second step, to distinguish between like and oppositely charged ions and then to get information on the inhomogeneity distribution of the like ions, we repeated the same calculations on the same sample configurations and removed one of the ions and considered only the other one. Detailed analysis of these distributions has revealed that the ring hydrogen atoms are mainly solvated by the anions, while the methyl and butyl terminal H atoms are surrounded by like atoms. The extent of this inhomogeneity was assessed by calculating the cluster size distribution that shows that the dimers are the most abundant ones.

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.

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.

New scaling relations to compute atom-in-material polarizabilities and dispersion coefficients: part 1. Theory and accuracy

Polarizabilities and London dispersion forces are important to many chemical processes. Force fields for classical atomistic simulations can be constructed using atom-in-material polarizabilities and C n (n = 6, 8, 9, 10…) dispersion coefficients. This article addresses the key question of how to efficiently assign these parameters to constituent atoms in a material so that properties of the whole material are better reproduced. We develop a new set of scaling laws and computational algorithms (called MCLF) to do this in an accurate and computationally efficient manner across diverse material types. We introduce a conduction limit upper bound and m-scaling to describe the different behaviors of surface and buried atoms. We validate MCLF by comparing results to high-level benchmarks for isolated neutral and charged atoms, diverse diatomic molecules, various polyatomic molecules (e.g., polyacenes, fullerenes, and small organic and inorganic molecules), and dense solids (including metallic, covalent, and ionic). We also present results for the HIV reverse transcriptase enzyme complexed with an inhibitor molecule. MCLF provides the non-directionally screened polarizabilities required to construct force fields, the directionally-screened static polarizability tensor components and eigenvalues, and environmentally screened C6 coefficients. Overall, MCLF has improved accuracy compared to the TS-SCS method. For TS-SCS, we compared charge partitioning methods and show DDEC6 partitioning yields more accurate results than Hirshfeld partitioning. MCLF also gives approximations for C8, C9, and C10 dispersion coefficients and quantum Drude oscillator parameters. This method should find widespread applications to parameterize classical force fields and density functional theory (DFT) + dispersion methods.

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.

Effective charges of ionic liquid determined self-consistently through combination of molecular dynamics simulation and density-functional theory

A self-consistent scheme combining the molecular dynamics (MD) simulation and density functional theory (DFT) was recently proposed to incorporate the effects of the charge transfer and polarization of ions into non-poralizable force fields of ionic liquids for improved description of energetics and dynamics. The purpose of the present work is to analyze the detailed setups of the MD/DFT scheme by focusing on how the basis set, exchange-correlation (XC) functional, charge-fitting method or force field for the intramolecular and Lennard-Jones interactions affects the MD/DFT results of 1,3-dimethylimidazolium bis(trifluoromethylsulfonyl) imide ( [C1mim][NTf2]) and 1-ethyl-3-methylimidazolium glycinate ( [C2mim][Gly]). It was found that the double-zeta valence polarized or larger size of basis set is required for the convergence of the effective charge of the ion. The choice of the XC functional was further not influential as far as the generalized gradient approximation is used. The charge-fitting method and force field govern the accuracy of the MD/DFT scheme, on the other hand. We examined the charge-fitting methods of Blöchl, the iterative Hirshfeld (Hirshfeld-I), and REPEAT in combination with Lopes et al.'s force field and general AMBER force field. There is no single combination of charge fitting and force field that provides good agreements with the experiments, while the MD/DFT scheme reduces the effective charges of the ions and leads to better description of energetics and dynamics compared to the original force field with unit charges. © 2017 Wiley Periodicals, Inc.

Local Structure in Terms of Nearest-Neighbor Approach in 1-Butyl-3-methylimidazolium-Based Ionic Liquids: MD Simulations

Description of the local microscopic structure in ionic liquids (ILs) is a prerequisite to obtain a comprehensive understanding of the influence of the nature of ions on the properties of ILs. The local structure is mainly determined by the spatial arrangement of the nearest neighboring ions. Therefore, the main interaction patterns in ILs, such as cation-anion H-bond-like motifs, cation-cation alkyl tail aggregation, and ring stacking, were considered within the framework of the nearest-neighbor approach with respect to each particular interaction site. We employed classical molecular dynamics (MD) simulations to study in detail the spatial, radial, and orientational relative distribution of ions in a set of imidazolium-based ILs, in which the 1-butyl-3-methylimidazolium (C4mim(+)) cation is coupled with the acetate (OAc(-)), chloride (Cl(-)), tetrafluoroborate (BF4(-)), hexafluorophosphate (PF6(-)), trifluoromethanesulfonate (TfO(-)), or bis(trifluoromethanesulfonyl)amide (TFSA(-)) anion. It was established that several structural properties are strongly anion-specific, while some can be treated as universally applicable to ILs, regardless of the nature of the anion. Namely, strongly basic anions, such as OAc(-) and Cl(-), prefer to be located in the imidazolium ring plane next to the C-H(2/4-5) sites. By contrast, the other four bulky and weakly coordinating anions tend to occupy positions above/below the plane. Similarly, the H-bond-like interactions involving the H(2) site are found to be particularly enhanced in comparison with the ones at H(4-5) in the case of asymmetric and/or more basic anions (C4mimOAc, C4mimCl, C4mimTfO, and C4mimTFSA), in accordance with recent spectroscopic and theoretical findings. Other IL-specific details related to the multiple H-bond-like binding and cation stacking issues are also discussed in this paper. The secondary H-bonding of anions with the alkyl hydrogen atoms of cations as well as the cation-cation alkyl chain aggregation turned out to be poorly sensitive to the nature of the anion.

Understanding SO2 Capture by Ionic Liquids

Ionic liquids have generated interest for efficient SO2 absorption due to their low vapor pressure and versatility. In this work, a systematic investigation of the structure, thermodynamics, and dynamics of SO2 absorption by ionic liquids has been carried out through quantum chemical calculations and molecular dynamics (MD) simulations. MP2 level calculations of several ion pairs complexed with SO2 reveal its preferential interaction with the anion. Results of condensed phase MD simulations of SO2-IL mixtures manifested the essential role of both cations and anions in the solvation of SO2, where the solute is surrounded by the "cage" formed by the cations (primarily its alkyl tail) through dispersion interactions. These structural effects of gas absorption are substantiated by calculated Gibbs free energy of solvation; the dissolution is demonstrated to be enthalpy driven. The entropic loss of SO2 absorption in ionic liquids with a larger anion such as [NTf2](-) has been quantified and has been attributed to the conformational restriction of the anion imposed by its interaction with SO2. SO2 loading IL decreases its shear viscosity and enhances the electrical conductivity. This systematic study provides a molecular level understanding which can aid the design of task-specific ILs as electrolytes for efficient SO2 absorption.

Capacitive Energy Storage: Current and Future Challenges

Capacitive energy storage devices are receiving increasing experimental and theoretical attention due to their enormous potential for energy applications. Current research in this field is focused on the improvement of both the energy and the power density of supercapacitors by optimizing the nanostructure of porous electrodes and the chemical structure/composition of the electrolytes. However, the understanding of the underlying correlations and the mechanisms of electric double layer formation near charged surfaces and inside nanoporous electrodes is complicated by the complex interplay of several molecular scale phenomena. This Perspective presents several aspects regarding the experimental and theoretical research in the field, discusses the current atomistic and molecular scale understanding of the mechanisms of energy and charge storage, and provides a brief outlook to the future developments and applications of these devices.

Divergent trend in density versus viscosity of ionic liquid/water mixtures: a molecular view from guanidinium ionic liquids

Ionic liquids (ILs) have shown great potential in the dissolution and stability of biomolecules when a low-to-moderate quantity of water is added.