Exploring the capabilities and limitations of the Van Hove function to understand directional correlations in ion movements within Li-ion battery electrolytes

Understanding microscopic directional correlations in ion movements within lithium-ion battery (LIB) electrolytes is important because these correlations directly affect the ionic conductivity. Onsager transport coefficients are widely used to understand these correlations. On the other hand, the Van Hove function (VHF) is also capable of determining correlated motions. However, identifying various types of ion correlated motions in LIB electrolytes using VHF is not well explored. Here, we have conducted molecular dynamics simulations of a representative experimental LIB electrolyte system-lithium hexafluorophosphate (LiPF6)-at different concentrations in a (9:1 wt. %) mixture of ethyl methyl carbonate and fluoroethylene carbonate in order to explore the capabilities and limitations of using VHF to understand different types of ion correlations. We conclude that analysis of VHF can qualitatively describe both the positive correlation between cation-anion at different salt concentrations and the negative correlation between cation-cation and anion-anion present in high salt concentration, but it cannot foretell which correlation is dominating at any given electrolyte concentration. This type of quantitative information can be obtained only via Onsager's approach. This could be seen as a limitation of relying solely on VHF to fully understand ion correlation in electrolyte media.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Tracing the Effect of Replacing [Gly]- with [Ala]- and Hydroxylation of [emim]+ on the Fine-Tuning of the Transport Properties of the Corresponding Amino Acid-Based Ionic Liquids Using MD Simulation

Natural amino acid-based ionic liquids (AAILs) composed of deprotonated amino acids, [AA]^-, as anions and hydroxylated imidazolium cations provide an eco-friendly nontoxic IL family with the growing number of chemical and biochemical revolutionary applications. In this paper, the transport properties of four AAILs composed of 1-(2-hydroxyethyl)-3-methylimidazolium ([HOemim]⁺) and 1-ethyl-3-methylimidazolium ([emim]⁺) cations with alaninate and glycinate anions were studied by molecular dynamics (MD) simulations. A nonpolarizable all-atom force field with the scaled charge (±0.8e) on each of the ions was applied and compared with the unit charge model in some cases. The tunable effects of the presence of the hydroxyl group in the side chain of the imidazolium cation, the type of amino acid anion, and the varied temperature on the dynamical behavior of AAILs were investigated in detail. The experimentally compatible trends of the simulated ionic self-diffusion coefficients, ionic conductivity, and ionicity were found to be inverse to the viscosity and ionic association of these ILs as [emim][Gly] > [emim][Ala] > [HOemim][Gly] > [HOemim][Ala]. The main reason behind these trends is the higher ability of the hydroxylated cation for the hydrogen-bond formation with [AA]^-. The mean square displacement (MSD), self-diffusion, and transference number of imidazolium cations are larger than those of [AA]^- anions, except in the case of [HOemim][Gly]. It was found that the activation energy for diffusion of [AA]^- is lower than that of [HOemim]⁺ but higher than that of [emim]⁺ in [HOemim][AA] and [emim][AA] ILs, respectively. The computed velocity autocorrelation function (VACF) showed that [Gly]^-, as the lightest ion, has the shortest mean collision time and velocity randomization time among the ions, especially in the [HOemim][Gly] IL. Replacing [emim]⁺ with [HOemim]⁺, similar to the effect of decreasing temperature, causes significant decreasing of the ionic self-diffusion and increasing of the well depth of the first minimum of the ionic VACFs. Current findings show that introducing suitable functional groups in the side chain of imidazolium cations can be a viable approach for efficient engineering design and fine-tuning of the transport properties of these AAILs.

The effect of explicit counterion binding on the transference number of polyelectrolyte solutions

Polyelectrolyte solutions have been proposed as a method to improve the efficiency of lithium-ion batteries by increasing the cation transference number because the polymer self-diffusion coefficient is much lower than that of the counterion. However, this is not necessarily true for the polymer mobility. In some cases, negative transference numbers have been reported, which implies that the lithium ions are transporting to the same electrode as the anion, behavior that is often attributed to a binding of counterions to the polyion. We use a simple model where we bind some counterions to the polymer via harmonic springs to investigate this phenomenon. We find that both the number of bound counterions and the strength of their binding alter the transference number, and, in some cases, the transference number is negative. We also investigate how the transference number depends on the Manning parameter, the ratio of the Bjerrum length to charge separation along the chain. By altering the Manning parameter, the transference number can almost be doubled, which suggests that charge spacing could be a way to increase the transference number of polyelectrolyte solutions.

Structure and Capacitance of Electrical Double Layers in Tricationic Ionic Liquids with Organic Solvents

Tricationic ionic liquid (TIL) electrolytes have been successfully employed in supercapacitors with graphene electrodes, but the low power density of the TILs-based supercapacitors caused by strong cations-anions associations requires enhancement by adding organic solvents to the liquid. In this paper, the role of the solvents acetonitrile (ACN) and ethylene carbonate (EC) on the ion diffusion, the conductivity of the TIL [C₆(mim)₃](Tf₂N)₃, and the structures and the capacitances of the electrical double layers (EDLs) in TIL/ACN and TIL/EC electrolytes were probed by molecular dynamics (MD) simulations. The results indicate that adding organic solvents to the liquid significantly reduces interactions between ions, thereby greatly improving the ion diffusion coefficients and the conductivity of the TIL, and the maximum conductivity is found at the 0.55 M TIL/ACN electrolyte concentration. Moreover, the reduced packing of counterions and the strong expulsion of coions near charged electrodes are observed in the organic electrolytes, especially in the TIL/EC electrolyte. Further analyses on EDLs affirm that the asymmetric camel-shaped differential capacitance-voltage (C-V) curve in the pure TIL electrolyte is weakly changed by the solvent ACN or EC. Besides, the EDL capacitance in the TIL-based hybrid electrolytes is improved slightly by the organic solvents. Comparing the integral capacitances in TIL/ACN and TIL/EC with different solvent contents, it is found that reducing the solvent polarity may be more beneficial to promote the EDL capacitance. Comprehensively, in this work, the 0.55 M TIL/ACN electrolyte is the optimal choice for the high-performance supercapacitor. Hence, solvating TIL electrolytes in supercapacitors by suitable solvents can effectively enhance the power density without compromising energy density.

Molecular Dynamics and Emerging Network Graphs of Interactions in Dinitrile-Based Li-Ion Battery Electrolytes

Advancements in battery research have shown interesting formulations of battery electrolytes that have helped improve the efficiency of Li-ion batteries over the decades. However, the quest for a safer and affordable battery electrolyte still proceeds with more unique formulations reported in the literature regularly. The dinitriles, especially adiponitrile and glutaronitrile, have caught the attention of the research community as part of this quest. In this work, we performed molecular dynamics simulations of dinitrile electrolytes with lithium bistrifluorosulfonimide (LiTFSI) as the electrolyte salt at varying concentrations and temperatures. On analysis of our simulations, we find that the densities of the mixtures follow the same trend as that of experimental values. The solvation properties were explored using the radial distribution functions. The connectivity of the Li⁺ with the dinitrile molecules and anions is established for all of the electrolyte concentrations using network graphs. We observe that the electrolytes form highly networked structures as the concentration increases without being affected by the rise in temperature. The networking of ionic interactions was quantified by calculating the average degree of each graph. Ionic conductivity calculations were computed using three methods: Nernst-Einstein relation, correlated method, and current autocorrelation function. We report the importance of accounting for the correlated motion of ions while estimating the ionic conductivity. The correlated conductivity and current autocorrelation function calculations provide a satisfactory estimation of the ionic conductivity compared to the experimental values.

Impact of Anion Asymmetry on Local Structure and Supercooling Behavior of Water-in-Salt Electrolytes

Salts with asymmetric (fluorosulfonyl)(trifluoromethanesulfonyl)imide (FTFSI) anions have recently been shown to suppress crystallization of water-in-salt electrolytes, enabling low-temperature operation of high-voltage aqueous rechargeable batteries. To clarify the underlying mechanism for the kinetic suppression of crystallization, we investigate the local solution structures and dynamic behaviors of water-in-salt electrolytes based on the asymmetric FTFSI anion and its symmetric anion analogues by Raman spectroscopy and molecular dynamics simulations. We find that monodentate coordination of FTFSI to cations leads to high rotational mobility of the uncoordinated SO₂CF₃ group. We conclude that the peculiar, coordination-dependent, local dynamics in the asymmetric FTFSI anion, manifested by enhanced intramolecular bond rotation, enables the strong supercooling behavior.

Fine probing the effect of replacing [PF6]- with [PF3(C2F5)3]- on the local structure and nanoscale organization of [bmim]+-based ionic liquids using MD simulation

Comparative all-atom molecular dynamics simulations are used to study the microscopic local structure and interionic interactions of two ionic liquids (ILs) composed of the 1-butyl-3-methylimidazolium cation, [bmim]+, coupled with the hexafluorophosphate, [PF6]-, or tris(pentafluoroethyl)trifluorophosphate, [FAP]-, anions. Respective distribution functions clearly reveal that the structural correlations between the cation and anion decrease when (i) replacing [PF6]- with [FAP]-, (ii) scaling the partial atomic charges, and (iii) considering the anion's structural flexibility versus rigidity. Replacement of [PF6]- with [FAP]- expands the nonpolar domains totally and causes the decreasing of the three-dimensional polar networks as well as the diminishing of the nano-aggregation of cation side chains. Current simulations show that with increasing the anion size and its charge delocalization, the probability of the in-plane cation-anion conformation, its related hydrogen bond acceptor ability, and the cation-cation π-π interaction decreases in accordance with the fluidity enhancements of the corresponding imidazolium-based IL. Hence, structural findings can explain and justify rationally the origins of the observed trends in the simulated dynamical properties of these ILs in our previous report. A complete understanding of the microscopic structure of ILs is necessary to control the outstanding properties of ILs as designer solvents that will support experimentalists for the best engineering design and a breakthrough efficiency of IL-related processes.

Insights into High Conductivity of the Two-Dimensional Iodine-Oxidized sp2-c-COF

A recent experiment [ Jin , E. ; Science 2017 , 357 , 673 - 676 ] shows that the conductivity of a two-dimensional (2D) sp²-carbon-hybridized π-conjugated covalent organic framework (sp²-c-COF) can be enhanced by as much as 12 orders of magnitude after iodine oxidation processing. To understand the physical mechanism underlying such a huge increase in the conductivity, we perform multiscale computations and find that the high conductivity of the iodine-oxidized 2D COF can be attributed to both hole transfer and ion transfer within the 2D COF. The computed dominant charge distribution corresponding to the valence band maximum (VBM) suggests that the delocalized π electrons occur mostly at the active reaction sites. The computed low ionization energy at the active reaction sites further supports that the 2D COF tends to lose electrons during iodine oxidation and to yield cationic COF and anionic triiodide I₃^-. Complementary classical molecular dynamics simulation shows a relatively high anion conductivity of 13.63 × 10^-2 S m^-1, consistent with the high conductivity measured from the experiment (7.1 × 10^-2 S m^-1). Meanwhile, we find that the cations in 2D COF can also induce a shift of the Fermi level to cross the valence band, thereby enhancing the hole mobility to 86.75 cm² V^-1 s^-1. For proposing a potential application of the highly conductive iodine-oxidized 2D sp²-c-COF, we design a prototypical model of the 2D spirally wound lithium-ion battery and find that it exhibits enhanced stability than a typical electrolyte material.

Capturing the effect of [PF3(C2F5)3]−vs. [PF6]−, flexible anion vs. rigid, and scaled charge vs. unit on the transport properties of [bmim]+-based ionic liquids: a comparative MD study

Effect of [PF₃(C₂F₅)₃]⁻vs. [PF₆]⁻, flexible anion vs. rigid, and scaled charge vs. unit on the transport properties of ILs.

Nanoscopic Study on Aliphatic Choline-Based Naphthenic Acid Ionic Liquids: Structural and Dynamical Properties

The structural and dynamical properties of two cholinium-based naphthenic acid ionic liquids (CBNAILs), namely, cholinium cyclopentane carboxylate ([CH][CPC]) and cholinium cyclohexane carboxylate ([CH][CHC]) have been studied using the classical molecular dynamics simulations and quantum mechanical methods. The results have also been compared with corresponding aromatic ones. Interaction energies, charge transfers, and hydrogen bondings between the studied ion pairs were investigated by density functional theory calculations and also the theory of atoms in molecules. Density, mean-square displacement, self-diffusivity, viscosity, electrical conductivity, transference number, ionicity, and fragility have been computed for the studied CBNAILs in the temperature range of 298.15-450 K and at 0.1 MPa. The simulated values were in good agreement with experimental data where they exist. The structural features of these CBNAILs were characterized by calculating the partial site-site radial distribution functions and spatial distribution functions. The results show a density cap of hydrogen atoms of hydroxyl groups of cations extended along the COO^- groups of anions. In these ILs, [CPC]^- and [CHC]^- are connected by hydrogen bonding to [CH]⁺ mainly through the carboxylate group. Increasing the number of carbon atoms in the ring of [CHC]^- with respect to [CPC]^- makes remarkable changes in self-diffusions, electrical conductivities, and viscosities. Altogether, this work gives a better insight into the dynamics and structuring of this class of biodegradable ILs at a molecular level.

Protic ammonium carboxylate ionic liquids: insight into structure, dynamics and thermophysical properties by alkyl group functionalization

This study is aimed at characterising the structure, dynamics and thermophysical properties of five alkylammonium carboxylate ionic liquids (ILs) from classical molecular dynamics simulations.

N,N-Dimethylpyridin-4-amine (DMAP) based ionic liquids: evaluation of physical properties via molecular dynamics simulations and application as a catalyst for Fisher indole and 1H-tetrazole synthesis

N,N-dimethylpyridin-4-amine (DMAP) based ionic liquids (ILs) as new and efficient catalysts for the facile synthesis of indoles (via Fischer indole synthesis), and 1H-tetrazoles (via click chemistry).

Capacitive performance of amino acid ionic liquid electrolyte-based supercapacitors by molecular dynamics simulation

Capacitance–electric potential curves of amino acid ionic liquid electrolyte-based supercapacitors.

Molecular dynamics simulation of amino acid ionic liquids near a graphene electrode: effects of alkyl side-chain length

The nanoscopic structure of amino acid ionic liquids (AAILs) as biodegradable electrolytes near a neutral graphene surface was studied by molecular dynamics (MD) simulation.

Viscosity scaling of the self-diffusion and velocity cross-correlation coefficients of two functionalised ionic liquids and of their non-functionalized analogues

Functionalised and unfunctionalised ammonium based ionic liquid analogues show a common Stokes–Einstein–Sutherland scaling with viscosity.

Thermodynamic studies of ionic hydration and interactions for amino acid ionic liquids in aqueous solutions at 298.15 K

Amino acid ionic liquids are a special class of ionic liquids due to their unique acid-base behavior, biological significance, and applications in different fields such as templates in synthetic chemistry, stabilizers for biological macromolecules, etc. The physicochemical properties of these ionic liquids can easily be altered by making the different combinations of amino acids as anion along with possible cation modification which makes amino acid ionic liquids more suitable to understand the different kinds of molecular and ionic interactions with sufficient depth so that they can provide fruitful information for a molecular level understanding of more complicated biological processes. In this context, volumetric and osmotic coefficient measurements for aqueous solutions containing 1-ethyl-3-methylimidazolium ([Emim]) based amino acid ionic liquids of glycine, alanine, valine, leucine, and isoleucine are reported at 298.15 K. From experimental osmotic coefficient data, mean molal activity coefficients of ionic liquids were estimated and analyzed using the Debye-Hückel and Pitzer models. The hydration numbers of ionic liquids in aqueous solutions were obtained using activity data. Pitzer ion interaction parameters are estimated and compared with other electrolytes reported in the literature. The nonelectrolyte contribution to the aqueous solutions containing ionic liquids was studied by calculating the osmotic second virial coefficient through an application of the McMillan-Mayer theory of solution. It has been found that the second osmotic virial coefficient which includes volume effects correlates linearly with the Pitzer ion interaction parameter estimated independently from osmotic data as well as the hydrophobicity of ionic liquids. The enthalpy-entropy compensation effect, explained using the Starikov-Nordén model of enthalpy-entropy compensation, and partial molar entropy analysis for aqueous [Emim][Gly] solutions are made by using experimental Gibb's free energy data and literature enthalpy data. This study highlights that the hydrophobic interaction persists even in the limit of infinite dilution where the hydration effects are usually dominant, implying importance of hydrophobic hydration. Analysis of the results further shows that the hydration of amino acid ionic liquids occurs through the cooperative H-bond formation with the kosmotropic effect in contrast to the usual inorganic salts or hydrophobic salts like tetraalkylammonium halides.

Amino acid ionic liquid modified mesoporous carbon: a tailor-made nanostructure biosensing platform

A novel nanocomposite based on ordered graphitized mesoporous carbon (GMC) and amino acid ionic liquids (AAIL) is obtained through controlled surface modification of GMC with hydrophilic AAILs (1-ethyl-3-methylimidazolium alanine, EMIM[Ala]), which is used as a platform for a tyrosinase biosensor to detect phenol. The GMC-AAIL nanocomposite possesses a better biocompatibility and improved aqueous-phase dispersion than hydrophobic GMC alone, owing to the introduction of hydrophilic and biocompatible AAILs. Comparative studies revealed that the catalytic activity of tyrosinase for phenol in phosphate buffer solution (PBS) containing EMIM[Ala] was about ten times higher than that in pure PBS. By entrapping tyrosinase molecules into the mesopores of GMC, making use of the synergy effect of GMC and AAIL (the "interspace confinement effect", the anti-fouling ability, and the biocompatible microenvironment), the GMC-AAIL-based biosensors display superior analytical performance to GMC-based ones in terms of signal-to-noise ratio, stability, repeatability, and working life. After 21-day storage, the electrode retained more than 90% of its initial response, indicating that surface modification of GMC with hydrophilic and biocompatible AAILs could significantly prolong the life of tyrosinase in vitro. The GMC10-EMIM[Ala]-based biosensor demonstrates a linear response for phenol concentrations from 0.1 to 10 µmol L(-1) with a low detection limit of 20 nmol L(-1) and sensitivity of 1385 mA cm(-2) M(-1). The GMC-AAIL nanocomposite proves to be a promising platform for enzyme-based biosensors and biocatalysis.

Acetonitrile boosts conductivity of imidazolium ionic liquids

We apply a new methodology in the force field generation (Phys. Chem. Chem. Phys.2011, 13, 7910) to study binary mixtures of five imidazolium-based room-temperature ionic liquids (RTILs) with acetonitrile (ACN). Each RTIL is composed of tetrafluoroborate (BF(4)) anion and dialkylimidazolium (MMIM) cations. The first alkyl group of MIM is methyl, and the other group is ethyl (EMIM), butyl (BMIM), hexyl (HMIM), octyl (OMIM), and decyl (DMIM). Upon addition of ACN, the ionic conductivity of RTILs increases by more than 50 times. It significantly exceeds an impact of most known solvents. Unexpectedly, long-tailed imidazolium cations demonstrate the sharpest conductivity boost. This finding motivates us to revisit an application of RTIL/ACN binary systems as advanced electrolyte solutions. The conductivity correlates with a composition of ion aggregates simplifying its predictability. Addition of ACN exponentially increases diffusion and decreases viscosity of the RTIL/ACN mixtures. Large amounts of ACN stabilize ion pairs, although they ruin greater ion aggregates.