Ternary Solid Polymer Electrolytes at the Electrochemical Interface: A Computational Study.

Polymer-based solid-like gel electrolytes have emerged as a promising alternative to improve battery performance. However, there is a scarcity of studies on the behavior of these media at the electrochemical interface. In this work, we report classical MD simulations of ternary polymer electrolytes composed of poly(ethylene oxide), a lithium salt [lithium bis(trifluoromethanesulfonyl)imide], and different ionic liquids [1-butyl-1-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide and 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide] confined between two charged and uncharged graphene-like surfaces. The molecular solvation of Li+ ions and their diffusion as well as the polymer conformational picture were characterized in terms of the radial distribution functions, coordination numbers, number density profiles, orientations, displacement variance, polymer radius of gyration, and polymer end-to-end distance. Our results show that the layering behavior of the ternary electrolyte in the interfacial region leads to a decrease of Li+ mobility in the direction perpendicular to the electrodes and high energy barriers that hinder lithium cations from coming into direct contact with the graphene-like surface. The nature of the ionic liquid and its concentration were found to influence the structural and dynamic properties at the electrode/electrolyte interface, the electrolyte with low amounts of the pyrrolidinium-based ionic liquid being that with the best performance since it favors the migration of Li+ cations toward the negative electrode when compared to the imidazolium-based one.

Size Matters: A Computational Study of Hydrogen Absorption in Ionic Liquids.

We combined both density functional theory and classical molecular dynamics simulations to investigate the molecular mechanisms governing hydrogen solvation in a total of 12 ionic liquids. Overall, the analysis of the structural properties under high temperature and pressure conditions revealed weak interactions between hydrogen and the ionic liquids, with a slight preference of this gas to be placed at the apolar domains. Interestingly, those ionic liquids comprising nitrate anions allow the accommodation of hydrogen molecules also in the polar areas. The study of the hydrogen velocity autocorrelation functions supports this observation. In addition, the structure of all of the tested ionic liquids was almost insensitive to the addition of hydrogen, so the available free volume and cavity formation are presumably the most important factors affecting solubility.

Ion-ion correlations across and between electrified graphene layers.

When an ionic liquid adsorbs onto a porous electrode, its ionic arrangement is deeply modified due to a screening of the Coulombic interactions by the metallic surface and by the confinement imposed upon it by the electrode's morphology. In particular, ions of the same charge can approach at close contact, leading to the formation of a superionic state. The impact of an electrified surface placed between two liquid phases is much less understood. Here we simulate a full supercapacitor made of the 1-butyl-3-methylimidazolium hexafluorophosphate and nanoporous graphene electrodes, with varying distances between the graphene sheets. The electrodes are held at constant potential by allowing the carbon charges to fluctuate. Under strong confinement conditions, we show that ions of the same charge tend to adsorb in front of each other across the graphene plane. These correlations are allowed by the formation of a highly localized image charge on the carbon atoms between the ions. They are suppressed in larger pores, when the liquid adopts a bilayer structure between the graphene sheets. These effects are qualitatively similar to the recent templating effects which have been reported during the growth of nanocrystals on a graphene substrate.

The effect of alkyl chain length on the structure and thermodynamics of protic-aprotic ionic liquid mixtures: a molecular dynamics study.

Mixtures of alkylammonium based protic ionic liquids and alkylmethylimidazolium based aprotic ionic liquids were studied by means of molecular dynamics simulations. Close to ideal mixing is observed in most studied magnitudes; however, the effect of increasing alkyl chain length in each of the cations is markedly different, with longer protic cations showing larger deviations, especially with regards to mixing enthalpy, which exhibits a strong compound forming tendency. The compound forming nature of these protic ionic liquids is shown to induce sharp changes in their local environment upon mixing.

Two-dimensional pattern formation in ionic liquids confined between graphene walls.

We perform molecular dynamics simulations of ionic liquids confined between graphene walls under a large variety of conditions (pure ionic liquids, mixtures with water and alcohols, mixtures with lithium salts and defective graphene walls). Our results show that the formation of striped and hexagonal patterns in the Stern layer can be considered as a general feature of ionic liquids at electrochemical interfaces, the transition between patterns being controlled by the net balance of charge in the innermost layer of adsorbed molecules. This explains previously reported experimental and computational results and, for the first time, why these pattern changes are triggered by any perturbation of the charge density at the innermost layer of the electric double layer (voltage and composition changes, and vacancies at the electrode walls, among others), which may help tuning electrode-ionic liquid interfaces. Using Monte Carlo simulations we show that such structures can be reproduced by a simple two-dimensional lattice model with only nearest-neighbour interactions, governed by highly screened ionic interactions and short-range and excluded volume interactions. We also show that the results of our simulations are consistent with those inferred from the Landau-Brazovskii theory of pattern formation in self-assembling systems. The presence of these patterns at the ionic liquid graphene-electrode interfaces may have a strong impact on the process of ionic transfer from the bulk mixtures to the electrodes, on the differential capacitance of the electrode-electrolyte double layer or on the rates of redox reactions at the electrodes, among other physicochemical properties, and is therefore an effect of great technological interest.

Confinement Effects on an Electron Transfer Reaction in Nanoporous Carbon Electrodes.

Nanoconfinement generally leads to a drastic effect on the physical and chemical properties of ionic liquids. Here we investigate how the electrochemical reactivity in such media may be impacted inside of nanoporous carbon electrodes. To this end, we study a simple electron transfer reaction using molecular dynamics simulations. The electrodes are held at constant electric potential by allowing the atomic charges on the carbon atoms to fluctuate. We show that the Fe3+/Fe2+ couple dissolved in an ionic liquid exhibits a deviation with respect to Marcus theory. This behavior is rationalized by the stabilization of a solvation state of the Fe3+ cation in the disordered nanoporous electrode that is not observed in the bulk. The simulation results are fitted with a recently proposed two solvation state model, which allows us to estimate the effect of such a deviation on the kinetics of electron transfer inside of nanoporous electrodes.

Nanostructured solvation in mixtures of protic ionic liquids and long-chained alcohols.

The structural and dynamical properties of bulk mixtures of long-chained primary and secondary alcohols (propanol, butanol, and 2-pentanol) with protic ionic liquids (ethylammonium and butylammonium nitrate) were studied by means of molecular dynamics simulations and small angle X-ray scattering (SAXS). Changes in the structure with the alcohol concentration and with the alkyl chain length of the alcohol moieties were found, showing variations in the radial distribution function and in the number of hydrogen bonds in the bulk liquids. Moreover, the structural behaviour of the studied mixtures is further clarified with the spatial distribution functions. The global picture in the local scale is in good agreement with the nanostructured solvation paradigm [T. Méndez-Morales et al. Phys. Chem. B 118, 761 (2014)], according to which alcohols are accommodated into the hydrogen bonds' network of the ionic liquid instead of forming clusters in the bulk. Indeed, our study reveals that the alcohol molecules are placed with their polar heads at the interfaces between polar and nonpolar nanodomains in the ionic liquid, with their alkyl chains inside the nonpolar organic nanodomains. The influence of alcohol chain length in the single-particle dynamics of the mixtures is also reported calculating the velocity autocorrelation function and vibrational densities of states of the different species in the ionic liquid-alcohol mixtures, and a weak caging effect for the ethylammonium cations independent of the chain size of the alcohols was found. However, the SAXS data collected for the studied mixtures show an excess of the scattering intensities which indicates that there are also some structural heterogeneities at the nanoscale.

Molecular dynamics analysis of the effect of electronic polarization on the structure and single-particle dynamics of mixtures of ionic liquids and lithium salts.

We report a molecular dynamics study on the effect of electronic polarization on the structure and single-particle dynamics of mixtures of the aprotic ionic liquid 1-ethyl-3-methylimidazolium bis-(trifluoromethylsulfonyl)-imide ([EMIM][TFSI]) doped with a lithium salt with the same anion at 298 K and 1 bar. In particular, we analyze the effect of electron density fluctuations on radial distribution functions, velocity autocorrelation functions, cage correlation functions, mean-squared displacements, and vibrational densities of states, comparing the predictions of the quantum-chemistry-based Atomistic Polarizable Potential for Liquids, Electrolytes, & Polymers (APPLE&P) with those of its nonpolarizable version and those of the standard non-polarizable Optimized Potentials for Liquid Simulations-All Atom (OPLS-AA). We found that the structure of the mixture is scarcely modified by the fluctuations in electron charge of their constituents, but their transport properties are indeed quite drastically changed, with larger mobilities being predicted for the different species in the bulk mixtures with the polarizable force field. Specifically, the mean-squared displacements are larger for the polarizable potentials at identical time intervals and the intermediate subdiffusive plateaus are greatly reduced, so the transition to the diffusive regime takes place much earlier than in the non-polarizable media. Moreover, the correlations of the added cations inside their cages are weakened out earlier and their vibrational densities of states are slightly red-shifted, reflecting the weakening effect of the electronic polarization on the Coulomb coupling in these dense ionic media. The comparison of OPLS-AA with non-polarizable APPLE&P indicates that adding polarization to OPLS-AA is not sufficient to achieve results close to experiments.

Molecular dynamics simulation of the structure and interfacial free energy barriers of mixtures of ionic liquids and divalent salts near a graphene wall.

A molecular dynamics study of mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIm][BF4]) with magnesium tetrafluoroborate (Mg[BF4]2) confined between two parallel graphene walls is reported. The structure of the system is analyzed by means of ionic density profiles, lateral structure of the first layer close to the graphene surface and angular orientations of imidazolium cations. Free energy profiles for divalent magnesium cations are calculated using two different methods in order to evaluate the height of the potential barriers near the walls, and the results are compared with those of mixtures of the same ionic liquid and a lithium salt (Li[BF4]). Preferential adsorption of magnesium cations is analyzed using a simple model and compared to that of lithium cations, and vibrational densities of states are calculated for the cations close to the walls analyzing the influence of the graphene surface charge. Our results indicate that magnesium cations next to the graphene wall have a roughly similar environment to that in the bulk. Moreover, they face higher potential barriers and are less adsorbed on the charged graphene walls than lithium cations. In other words, magnesium cations have a more stable solvation shell than lithium ones.

Molecular dynamics simulations of mixtures of protic and aprotic ionic liquids.

Molecular dynamics simulations of mixtures of the protic ionic liquid ethylammonium nitrate (EAN) and the aprotic 1-ethyl-3-methylimidazolium tetrafluoroborate ([EMIM][BF4]) are reported and the results are compared with experimental density and electrical conductivity measurements. Essentially ideal mixing of the ionic liquids is seen to take place by means of experimental and simulated excess molar volumes, whose very low values suggest a gradual transition between the structures of the two end constituents of the mixture. A weak dominance of the structure of the protic ionic liquid is nevertheless registered, due to a slight preferential formation of the network of hydrogen bonds, as reflected in the coordination number and the number of hydrogen bonds in the mixture. A novel conductivity curve showing pronounced deviations from the simple ideal mixing rule is reported, with three different regions defined by a local maximum - reflecting enhanced translational dynamics relative to ideal mixture behaviour - and a global minimum at intermediate concentrations. The physical origin of this behaviour is discussed along with the structure and single-particle dynamics of the mixture, and it is seen that these regions are defined by the onset of the formation of the EAN hydrogen bonded network (xEAN = 0.2) and the virtual disappearance of the structure of the aprotic ionic liquid at xEAN = 0.7. It is concluded that the delicate interplay between both networks has a deep effect on the placement and mobility of [EMIM](+) cations in the mixture all throughout the different stages of the structural transition, which seems to be the driving force behind the reported transport properties of the mixture at intermediate to high EAN concentrations.

Molecular origin of high free energy barriers for alkali metal ion transfer through ionic liquid-graphene electrode interfaces.

In this work we study mechanisms of solvent-mediated ion interactions with charged surfaces in ionic liquids by molecular dynamics simulations, in an attempt to reveal the main trends that determine ion-electrode interactions in ionic liquids. We compare the interfacial behaviour of Li(+) and K(+) at a charged graphene sheet in a room temperature ionic liquid, 1-butyl-3-methylimidazolium tetrafluoroborate, and its mixtures with lithium and potassium tetrafluoroborate salts. Our results show that there are dense interfacial solvation structures in these electrolytes that lead to the formation of high free energy barriers for these alkali metal cations between the bulk and direct contact with the negatively charged surface. We show that the stronger solvation of Li(+) in the ionic liquid leads to the formation of significantly higher interfacial free energy barriers for Li(+) than for K(+). The high free energy barriers observed in our simulations can explain the generally high interfacial resistance in electrochemical storage devices that use ionic liquid-based electrolytes. Overcoming these barriers is the rate-limiting step in the interfacial transport of alkali metal ions and, hence, appears to be a major drawback for a generalised application of ionic liquids in electrochemistry. Some plausible strategies for future theoretical and experimental work for tuning them are suggested.

Nanostructure of mixtures of protic ionic liquids and lithium salts: effect of alkyl chain length.

The bulk structure of mixtures of two protic ionic liquids, propylammonium nitrate and butylammonium nitrate, with a salt with a common anion, is analyzed at room temperature by means of small angle X-ray scattering and classical molecular dynamics simulations. The study of several structural properties, such as density, radial distribution functions, spatial distribution functions, hydrogen bonds, coordination numbers and velocity autocorrelation functions, demonstrates that increasing the alkyl chain length of the alkylammonium cation results in more segregated, better defined polar and apolar domains, the latter having a larger size. This increase, ascribed to the erosion of the H-bond network in the ionic liquid polar regions as salt is added, is confirmed by means of small angle X-ray scattering measurements, which show a clear linear increase of the characteristic spatial sizes of the studied protic ionic liquids with salt concentration, similar to that previously reported for ethylammonium nitrate (J. Phys. Chem. B, 2014, 118, 761-770). In addition, larger ionic liquid cations lead to a lower degree of hydrogen bonding and to more sparsely packed three-dimensional structures, which are more easily perturbed by the addition of lithium salts.

Mixtures of protic ionic liquids and molecular cosolvents: a molecular dynamics simulation.

In this work, the effect of molecular cosolvents (water, ethanol, and methanol) on the structure of mixtures of these compounds with a protic ionic liquid (ethylammonium nitrate) is analyzed by means of classical molecular dynamics simulations. Included are as-yet-unreported measurements of the densities of these mixtures, used to test our parameterized potential. The evolution of the structure of the mixtures throughout the concentration range is reported by means of the calculation of coordination numbers and the fraction of hydrogen bonds in the system, together with radial and spatial distribution functions for the various molecular species and molecular ions in the mixture. The overall picture indicates a homogeneous mixing process of added cosolvent molecules, which progressively accommodate themselves in the network of hydrogen bonds of the protic ionic liquid, contrarily to what has been reported for their aprotic counterparts. Moreover, no water clustering similar to that in aprotic mixtures is detected in protic aqueous mixtures, but a somehow abrupt replacing of [NO3](-) anions in the first hydration shell of the polar heads of the ionic liquid cations is registered around 60% water molar concentration. The spatial distribution functions of water and alcohols differ in the coordination type, since water coordinates with [NO3](-) in a bidentate fashion in the equatorial plane of the anion, while alcohols do it in a monodentate fashion, competing for the oxygen atoms of the anion. Finally, the collision times of the different cosolvent molecules are also reported by calculating their velocity autocorrelation functions, and a caging effect is observed for water molecules but not in alcohol mixtures.

Molecular dynamics simulations of the structure of the graphene-ionic liquid/alkali salt mixtures interface.

We performed molecular dynamics simulations of mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate with lithium tetrafluoroborate and potassium tetrafluoroborate between two charged and uncharged graphene walls, in order to analyze the structure of the well-known formation of layers that takes place on liquids under confinement. For this purpose, we studied the molecular density profiles, free energy profiles for bringing lithium and potassium cations from the bulk mixture to the graphene wall and the orientational distributions of imidazolium rings within the first adsorbed layer as a function of salt concentration and electrode potential. The charge densities in the electrodes were chosen to be zero and ±1 e nm(-2), and the salt molar percentages were %salt = 0, 10 and 25. We found that the layered structure extends up to 1-2 nm, where the bulk behaviour is recovered. In addition, whereas for the neutral surface the layers are composed of both ionic species, increasing the electrode potential, the structure changes to alternating cationic and anionic layers leading to an overcompensation of the charge of the previous layer. We also calculated the distribution of angles of imidazolium rings near neutral and charged graphene walls, finding a limited influence of the added salt. In addition, the average tilt of the imidazolium ring within the first layer goes from 36° with respect to a normal vector to the uncharged graphene wall to 62° in the presence of charged walls. The free energy profiles revealed that lithium and potassium ions are adsorbed on the negative surface only for the highest amount of salt, since the free energy barriers for approaching this electrode are considerably higher than kBT.

Solvation of lithium salts in protic ionic liquids: a molecular dynamics study.

The structure of solutions of lithium nitrate in a protic ionic liquid with a common anion, ethylammonium nitrate, at room temperature is investigated by means of molecular dynamics simulations. Several structural properties, such as density, radial distribution functions, hydrogen bonds, spatial distribution functions, and coordination numbers, are analyzed in order to get a picture of the solvation of lithium cations in this hydrogen-bonded, amphiphilically nanostructured environment. The results reveal that the ionic liquid mainly retains its structure upon salt addition, the interaction between the ammonium group of the cation and the nitrate anion being only slightly perturbed by the addition of the salt. Lithium cations are solvated by embedding them in the polar nanodomains of the solution formed by the anions, where they coordinate with the latter in a solid-like fashion reminiscent of a pseudolattice structure. Furthermore, it is shown that the average coordination number of [Li](+) with the anions is 4, nitrate coordinating [Li](+) in both monodentate and bidentate ways, and that in the second coordination layer both ethylammonium cations and other lithiums are also found. Additionally, the rattling motion of lithium ions inside the cages formed by their neighboring anions, indicative of the so-called caging effect, is confirmed by the analysis of the [Li](+) velocity autocorrelation functions. The overall picture indicates that the solvation of [Li](+) cations in this amphiphilically nanostructured environment takes place by means of a sort of inhomogeneous nanostructural solvation, which we could refer to as nanostructured solvation, and which could be a universal solvation mechanism in ionic liquids.

MD simulations of the formation of stable clusters in mixtures of alkaline salts and imidazolium-based ionic liquids.

Structural and dynamical properties of room-temperature ionic liquids containing the cation 1-butyl-3-methylimidazolium ([BMIM](+)) and three different anions (hexafluorophosphate, [PF6](-), tetrafluoroborate, [BF4](-), and bis(trifluoromethylsulfonyl)imide, [NTf2](-)) doped with several molar fractions of lithium salts with a common anion at 298.15 K and 1 atm were investigated by means of molecular dynamics simulations. The effect of the size of the salt cation was also analyzed by comparing these results with those for mixtures of [BMIM][PF6] with NaPF6. Lithium/sodium solvation and ionic mobilities were analyzed via the study of radial distribution functions, coordination numbers, cage autocorrelation functions, mean-square displacements (including the analysis of both ballistic and diffusive regimes), self-diffusion coefficients of all the ionic species, velocity and current autocorrelation functions, and ionic conductivity in all the ionic liquid/salt systems. We found that lithium and sodium cations are strongly coordinated in two different positions with the anion present in the mixture. Moreover, [Li](+) and [Na](+) cations were found to form bonded-like, long-lived aggregates with the anions in their first solvation shell, which act as very stable kinetic entities within which a marked rattling motion of salt ions takes place. With very long MD simulation runs, this phenomenon is proved to be on the basis of the decrease of self-diffusion coefficients and ionic conductivities previously reported in experimental and computational results.

Investigation of the local structure of mixtures of an ionic liquid with polar molecular species through molecular dynamics: cluster formation and angular distributions.

In this work, we used molecular dynamics simulations to analyze in detail the spatial distributions of the different constituents in mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate with three polar molecular species: water and two alcohols of different chain lengths (methanol and ethanol). In particular, we report results regarding the influence of the chosen species and its concentration on the formation of ionic and molecular clusters over the whole miscibility range, as well as on the angular distribution of polar molecules around the anion and the cation in these systems. Both analyses showed that addition of a molecular species breaks down the polar network of the pure ionic liquid in clusters whose mean size decreases progressively as more molecules are added. At very high concentrations of the molecular species, the ions are found to be isolated in mixtures with water and methanol, but they tend to form pairs in ethanol. In mixtures with water we identified large clusters that form a water network at very high water concentrations, while at low water concentrations polar molecules tend to form smaller aggregates. In contrast, in mixtures with alkanols there is no evidence of the formation of large alcohol clusters at any concentration. Spatial order in alcohol was also studied by means of the Kirkwood G factor, reaching the conclusion that the angular correlations which appear in pure alcohols due to dipole interactions are destroyed by the ionic liquid, even when present only in tiny amounts.

Dynamical properties of alcohol + 1-hexyl-3-methylimidazolium ionic liquid mixtures: a computer simulation study.

In this work, extensive molecular dynamics simulations of the dynamics of mixtures of ionic liquids (ILs) composed of the cation 1-hexyl-3-methylimidazolium and several anions of different hydrophobicity degrees (Cl(-), BF(4)(-), PF(6)(-)) with alcohols of different chain lengths (methanol and ethanol) are reported. We evaluated the influence of the nature of the anion, the length of the molecular chain of the alcohol, and the alcohol concentration on some dynamical properties of the mixtures, such as self-diffusion coefficients of all the species, mean square displacements (with an analysis of both ballistic and diffusive regimes), and velocity autocorrelation functions of alcohol molecules. The diffusivity of the mixtures was found to be highly dependent on the nature of the anion since the interaction between chloride and alcohols is greater than that with fluorinated anions and leads to slower dynamics. Additionally, our results show that self-diffusion coefficients increase with alcohol concentration. On the other hand, a subdiffusive regime over thousands of picoseconds was detected at intermediate times through analysis of the center-of-mass mean square displacements of alcohol molecules, a region that becomes narrower as alcohol concentration increases. Finally, the study of the role of the anion and of solvent concentration on velocity autocorrelation functions reflects an increase in mean collision times as the amount of alcohol increases until the value of pure alcohols is reached. These collision times are smaller in mixtures with halogenated ILs.

Molecular dynamics simulation of the structure and dynamics of water-1-alkyl-3-methylimidazolium ionic liquid mixtures.

We have performed extensive molecular dynamic simulations to analyze the influence of cation and anion natures, and of water concentration, on the structure and dynamics of water-1-alkyl-3-methylimidazolium ionic liquid mixtures. The dependence on water concentration of the radial distribution functions, coordination numbers, and hydrogen bonding degree between the different species has been systematically analyzed for different lengths of the cation alkyl chain (alkyl = ethyl, butyl, hexyl, and octyl) and several counterions. These include two halogens of different sizes and positions in Hoffmeister series, Cl(-) and Br(-), and the highly hydrophobic inorganic anion PF(6)(-) throughout its whole solubility regime. The formation of water clusters in the mixture has been verified, and the influences of both anion hydrophobicity and cation chain length on the structure and size of these clusters have been analyzed. The water cluster size is shown to be relatively independent of the cation chain length, but strongly dependent on the hydrophobicity of the anion, which also determines critically the network formation of water and therefore the miscibility of the ionic liquid. The greater influence of the anion relative to the cation one is seen to be reflected in all the analyzed physical properties. Finally, single-particle dynamics in IL-water mixtures is considered, obtaining the self-diffusion coefficients and the velocity autocorrelation functions of water molecules in the mixture, and analyzing the effect of cation, anion, and water concentration on the duration of the ballistic regime and on the time of transition to the diffusive regime. Complex non-Markovian behavior was detected at intermediate times within an interval progressively shorter as water concentration increases.