Mechanisms of ion transport in lithium salt-doped polymeric ionic liquid electrolytes at higher salt concentrations

We used atomistic simulations to study the mechanisms of ion transport in salt-doped polymeric ionic liquid systems at higher salt concentrations. Consistent with the experimental observations, our simulations indicate that at higher salt concentrations, the anion mobilities become lower than that of the lithium cations. Further, the anion mobilities become relatively insensitive to the salt concentration, while the mobilities of lithium increase with increasing salt concentration. We rationalize the results for the anion mobilities by considering the fractions of anions which are exclusively coordinated with the polycations (Type1); co-coordinated with cations and lithium (Type2); and those exclusively coordinated with lithium (Type3). By considering the coordination characteristics of the different types of anions and their hopping motions, we demonstrate that the net anion mobilities results from a compensation effect of the salt concentration dependence of the mobilities of the different anions. With respect to the mobilities of the lithium ions, we demonstrate that the latter moves primarily by a structural diffusion mechanism involving refreshing of the solvation shell during hopping. Further, for the majority of the lithium ions, the solvation shell is shown to be comprised of co-coordinated Type2 anions, and that the number of polycations and the unique polymer chains involved in such coordination decreases with increasing salt concentration. Such changes are shown to weaken the solvation shell around the lithium, thereby facilitating faster ion motion. Together, our results suggest that systems in which the anion which exhibits a stronger coordination to the polycation in comparison to that of the lithium can facilitate higher transference numbers without a concomitant reduction in the mechanical strength.     

Ion and polymer dynamics in polymer electrolytes PPO-LiClO4. I. Insights from NMR line-shape analysis

We investigate ion and polymer dynamics in polymer electrolytes PPO-LiClO4 performing 2H and 7Li NMR line-shape analysis. Comparison of temperature dependent 7Li and 2H NMR spectra gives evidence for a coupling of ion and polymer dynamics. 2H NMR spectra for various salt concentrations reveal a strong slowdown of the polymer segmental motion when the salt content is increased. The 2H NMR line shape further indicates that the segmental motion is governed by dynamical heterogeneities. While the width of the distribution of correlation times G(log tau) is moderate for low and high salt content, an extremely broad distribution exists for an intermediate salt concentration of 15:1 PPO-LiClO4. For the latter composition, a weighted superposition of two spectral components, reflecting the fast and the slow polymer segments of the distribution, describes the 2H NMR line shape over a broad temperature range. Analysis of the temperature dependent relative intensity of both spectral components indicates the existence of a continuous rather than a discontinuous distribution G(log tau). Such continuous distribution is consistent with gradual fluctuations of the local salt concentration and, hence, of the local environments of the polymer segments, whereas it is at variance with the existence of large salt-depleted and salt-rich domains featuring fast and slow polymer dynamics, respectively. Finally, for all studied PPO-LiClO4 mixtures, the 2H NMR line shape strongly depends on the echo delay in the applied echo-pulse sequence, indicating that the structural relaxation of the polymer segments involves successive rotational jumps about small angles gamma < 20 degrees .     

Molecular dynamics simulation of polymer electrolytes based on poly(ethylene oxide) and ionic liquids. I. Structural properties

Molecular dynamics (MD) simulations have been performed for prototype models of polymer electrolytes in which the salt is an ionic liquid based on 1-alkyl-3-methylimidazolium cations and the polymer is poly(ethylene oxide), PEO. The MD simulations were performed by combining the previously proposed models for pure ionic liquids and polymer electrolytes containing simple inorganic ions. A systematic investigation of ionic liquid concentration, temperature, and the 1-alkyl- chain length, [1,3-dimethylimidazolium]PF6, and [1-butyl-3-methylimidazolium]PF6, effects on resulting equilibrium structure is provided. It is shown that the ionic liquid is dispersed in the polymeric matrix, but ionic pairs remain in the polymer electrolyte. Imidazolium cations are coordinated by both the anions and the oxygen atoms of PEO chains. Probability density maps of occurrences of nearest neighbors around imidazolium cations give a detailed physical picture of the environment experienced by cations. Conformational changes on PEO chains upon addition of the ionic liquid are identified. The equilibrium structure of simulated systems is also analyzed in reciprocal space by using the static structure factor, S(k). Calculated S(k) display a low wave-vector peak, indicating that spatial correlation in an extended-range order prevail in the ionic liquid polymer electrolytes. Long-range correlations are assigned to nonuniform distribution of ionic species within the simulation box.     

Ion and polymer dynamics in polymer electrolytes PPO-LiClO4. II. 2H and 7Li NMR stimulated-echo experiments

We use 2H NMR stimulated-echo spectroscopy to measure two-time correlation functions characterizing the polymer segmental motion in polymer electrolytes PPO-LiClO4 near the glass transition temperature Tg. To investigate effects of the salt on the polymer dynamics, we compare results for different ether oxygen to lithium ratios, namely, 6:1, 15:1, 30:1, and infinity. For all compositions, we find nonexponential correlation functions, which can be described by a Kohlrausch function. The mean correlation times show quantitatively that an increase of the salt concentration results in a strong slowing down of the segmental motion. Consistently, for the high 6:1 salt concentration, a high apparent activation energy Ea=4.1 eV characterizes the temperature dependence of the mean correlation times at Tg相似文献   

The study of structure formation in a polymer-containing ionic liquid in terms of the integral equation theory

The structural properties of a polymer-containing ionic liquid under the conditions of good solubility of a flexible polymer are studied theoretically. Two systems are discussed: In one, polymer solubility is due to the presence of specific interaction between polymer chains and solvent cations; in the other, polymer solubility is due to the presence of specific interactions between the polymer and solvent anions. The dependences of the structural characteristics of a solution on the polymer concentration and the energy of attraction between polymer chains and solvent ions are calculated. In a semidilute polymer solution, long-range correlations of polymer chains with a power dependence of the characteristic scale of ordering on the polymer density appear. The conditions under which, along with the intermediate ordering typical of a pure ionic liquid, the long-range ordering of the solvent cations and anions occur after addition of a polymer to the ionic liquid are studied.     

Dependence of ion hydration on the sign of the ion's charge

The solvation of simple ions in water is studied using molecular dynamics simulations with a polarizable force field. Previous simulations using this potential demonstrated that anions are more favorably solvated in water than cations. The present work is an attempt to explain this result by examining the effects of ions on the surrounding water structure, with particular focus on the first solvation shell and its interactions with the surrounding water. We conclude that while the first solvation shell surrounding cations is frustrated by competition between ion-water and water-water interactions, solvation of anions is compatible with good water-water interactions.     

Simulation Study of the Conformational Properties of Diblock Polyelectrolytes in Salt Solutions

Coarse-grained molecular dynamics simulations are performed to understand the behavior of diblock polyelectrolytes in solutions of divalent salt by studying the conformations of chains over a wide range of salt concentrations. The polymer molecules are modeled as bead spring chains with different charged fractions and the counterions and salt ions are incorporated explicitly. Upon addition of a divalent salt, the salt cations replace the monovalent counterions, and the condensation of divalent salt cations onto the polyelectrolyte increases, and the chains favor to collapse. The condensation of ions changes with the salt concentration and depends on the charged fraction. Also, the degree of collapse at a given salt concentration changes with the increasing valency of the counterion due to the bridging effect. As a quantitative measure of the distribution of counterions around the polyelectrolyte chain, we study the radial distribution function between monomers on different polyelectrolytes and the counterions inside the counterion worm surrounding a polymer chain at different concentrations of the divalent salt. Our simulation results show a strong dependence of salt concentration on the conformational properties of diblock copolymers and indicate that it can tune the self-assembly behaviors of such charged polyelectrolyte block copolymers.     

Aqueous solutions of ionic liquids: study of the solution/vapor interface using molecular dynamics simulations

We performed a detailed molecular dynamics study of the interfacial structure of aqueous solutions of 1-butyl-3-methylimidazolium tetrafluoroborate in order to explain the anomalous dependence of the surface tension on concentration. At low concentrations the surface tension decreases with concentration. At higher concentrations the surface becomes saturated; a plateau is observed in simulations with a non-polarizable force field while a possible increase is detected in simulations with a polarizable force field. The structure is characterized by a surplus of cations at the surface (with hydrophobic butyl chains pointing toward vacuum) which at low concentrations is only partly compensated by the anions because of asymmetric solvation. A more hydrophobic 1-butyl-3-methylimidazolium hexafluorophosphate is also simulated for comparison.     

Highly concentrated polycarbonate‐based solid polymer electrolytes having extraordinary electrochemical stability

Solid polymer electrolytes (SPEs) are compounds of great interest as safe and flexible alternative ionics materials, particularly suitable for energy storage devices. We study an unusual dependence on the salt concentration of the ionic conductivity in an SPE system based on poly(ethylene carbonate) (PEC). Dielectric relaxation spectroscopy reveals that the ionic conductivity of PEC/lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) electrolyte continues to increase with increasing salt concentration because the segmental motion of the polymer chains is enhanced by the plasticizing effect of the imide anion. Fourier transfer‐infrared (FTIR) spectroscopy suggests that this unusual phenomenon arises because of a relatively loose coordination structure having moderately aggregated ions, in contrast to polyether‐based systems. Comparative FTIR study against PEC/lithium perchlorate (LiClO₄) electrolytes suggests that weak ionic interaction between Li and TFSI ions is also important. Highly concentrated electrolytes with both reasonable conductivity and high lithium transference number (t₊) can be obtained in the PEC/LiTFSI system as a result of the unusual salt concentration dependence of the conductivity and the ionic solvation structure. The resulting concentrated PEC/LiTFSI electrolytes have extraordinary oxidation stability and prevent any Al corrosion reaction in a cyclic voltammetry. These are inherent effects of the highly concentrated salt. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 2442–2447     

Effects of alkali metal halide salts on the hydrogen bond network of liquid water

Measurements of the oxygen K-edge X-ray absorption spectrum (XAS) of aqueous sodium halide solutions demonstrate that ions significantly perturb the electronic structure of adjacent water molecules. The addition of halide salts to water engenders an increase in the preedge intensity and a decrease in the postedge intensity of the XAS, analogous to those observed when increasing the temperature of pure water. The main-edge feature exhibits unique behavior and becomes more intense when salt is added. Density functional theory calculations of the XAS indicate that the observed red shift of the water transitions as a function of salt concentration arises from a strong, direct perturbation of the unoccupied molecular orbitals on water by anions, and does not require significant distortion of the hydrogen bond network beyond the first solvation shell. This contrasts the temperature-dependent spectral variations, which result primarily from intensity changes of specific transitions due to geometric rearrangement of the hydrogen bond network.     

Specific Ion and Concentration Effects in Acetate Solutions with Na+, K+ and Cs+

How salt ions affect solutes and the water beyond the solvation shell is not well understood. Molecular dynamics simulations of alkali-acetate solutions were analysed here in order to examine if, and how, different cations and solute concentrations affect the water structure and the interactions between water and acetates. The results revealed that water structure is perturbed to more than 1 nm away from the acetates and that this effect is more pronounced in physiological than in molar electrolyte concentrations. Analysis of simulations of a soluble protein revealed that the water orientation is perturbed to at least 1.5 nm from the protein structure. Furthermore, modifications to the orientation of water around carboxylate side chains were shown to depend on the local environment on the protein surface, and could extend to well over 1 nm, which may have an effect on protein dynamics during MD simulations in small water boxes.     

How ions affect the structure of water

We model ion solvation in water. We use the MB model of water, a simple two-dimensional statistical mechanical model in which waters are represented as Lennard-Jones disks having Gaussian hydrogen-bonding arms. We introduce a charge dipole into MB waters. We perform (NPT) Monte Carlo simulations to explore how water molecules are organized around ions and around nonpolar solutes in salt solutions. The model gives good qualitative agreement with experiments, including Jones-Dole viscosity B coefficients, Samoilov and Hirata ion hydration activation energies, ion solvation thermodynamics, and Setschenow coefficients for Hofmeister series ions, which describe the salt concentration dependence of the solubilities of hydrophobic solutes. The two main ideas captured here are (1) that charge densities govern the interactions of ions with water, and (2) that a balance of forces determines water structure: electrostatics (water's dipole interacting with ions) and hydrogen bonding (water interacting with neighboring waters). Small ions (kosmotropes) have high charge densities so they cause strong electrostatic ordering of nearby waters, breaking hydrogen bonds. In contrast, large ions (chaotropes) have low charge densities, and surrounding water molecules are largely hydrogen bonded.     

Effects of ion solvation on phase equilibrium and interfacial tension of liquid mixtures

We study the bulk thermodynamics and interfacial properties of electrolyte solution mixtures by accounting for electrostatic interaction, ion solvation, and inhomogeneity in the dielectric medium in the mean-field framework. Difference in the solvation energy between the cations and anions is shown to give rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The ion solvation affects the phase equilibrium of the solvent mixture, depending on the dielectric constants of the solvents, reflecting the competition between the solvation energy and translation entropy of the ions. Miscibility is decreased if both solvents have low dielectric constants and is enhanced if both solvents have high dielectric constant. At the mean-field level, the ion distribution near the interface is determined by two competing effects: accumulation in the electrostatic double layer and depletion in a diffuse interface. The interfacial tension shows a nonmonotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations and decreases approximately as the square root of the salt concentration for dilute solutions, reaching a minimum near 1 mM. We also find that, for a fixed cation type, the interfacial tension decreases as the size of anion increases. These results offer qualitative explanations within one unified framework for the long-known concentration and ion size effects on the interfacial tension of electrolyte solutions.     

Laser Raman and FTIR studies on Li+ interaction in PVAc-LiClO4 polymer electrolytes

The polymer electrolytes composed of poly(vinyl acetate) (PVAc) with various stoichiometric ratios of lithium perchlorate (LiClO(4)) salt have been prepared by solution casting method. The techniques Fourier transform infra-red (FTIR) and Laser Raman spectroscopy have been used to monitor polymer-salt complex formation, ion-ion and ion-polymer interactions as a function of salt concentration. Significant changes in both Laser Raman and FTIR spectra are observed which reveals an interaction between ester oxygens with lithium cation coordination. These results strongly suggest the interaction of lithium cation and network polymer chains. When the salt content is increased, the intensity of the internal Raman modes of the ClO(4)(-) increases. The ClO(4)(-) stretching mode observed at 934 cm(-1) in Laser Raman shows some additional shoulder peaks with increase in salt concentration. This reveals the presence of free anions, ion contact pairs and higher order ionic clusters. From the FTIR and Laser Raman results the transport mechanism of ions in PVAc:LiClO(4) polymer electrolytes has been discussed.     

Nematogens containing oxyethylene units at a lateral or terminal position and their mixtures with salts

Short terminal polyoxyethylene chains may be introduced into mesogens containing lateral substituents such as two hexyloxy chains or a crown-ether fragment. These compounds have a large nematic range near room temperature and can dissolve large amounts of the salt LiBF₄ without destroying the nematic arrangment. In the nematic phase, the ions are ordered and the quadrupolar splitting associated with these ions can allow this ordering to be monitored. Increasing the salt concentration does not seem to change the ordering of the individual ions. The ¹³C thermal evolution of the field-induced chemical shift in the oxyethylene (OE) units with or without salt does not show any real difference, indicating that the interaction between the ions and the mesogen is small. This means that there is very little change in conformation in the OE unit when adding salt to the mesogen.     

Nematogens containing oxyethylene units at a lateral or terminal position and their mixtures with salts

Short terminal polyoxyethylene chains may be introduced into mesogens containing lateral substituents such as two hexyloxy chains or a crown-ether fragment. These compounds have a large nematic range near room temperature and can dissolve large amounts of the salt LiBF₄ without destroying the nematic arrangment. In the nematic phase, the ions are ordered and the quadrupolar splitting associated with these ions can allow this ordering to be monitored. Increasing the salt concentration does not seem to change the ordering of the individual ions. The ¹³C thermal evolution of the field-induced chemical shift in the oxyethylene (OE) units with or without salt does not show any real difference, indicating that the interaction between the ions and the mesogen is small. This means that there is very little change in conformation in the OE unit when adding salt to the mesogen.     

On the molecular mechanism of ion specific Hofmeister series

Hofmeister series ranks the ability of salt ions in influencing a variety of properties and processes in aqueous solutions.In this review,we reexamine how these ions and some other small molecules affect water structure and thermodynamic properties,such as surface tension and protein backbone solvation.We illustrate the difficulties in interpreting the thermodynamic information based on structural and dynamic arguments.As an alternative,we show that the solvation properties of ions and proteins/small molecules can be used to explain the salt effects on the thermodynamic properties of the solutions.Our analysis shows that the often neglected cation-anion cooperativity plays a very important role in these effects.We also argue that the change of hydrogen donor/acceptor equilibrium by added cosolutes/cosolvents can be used to explain their effects on protein secondary structure denaturation/protection:those increase hydrogen donor concentrations such as urea and salts with strongly solvated cations/weakly hydrated anions tend to dissolve protein backbone acting as secondary structure denaturants,whereas those lack of hydrogen donors but rich in acceptors have the opposite effect.     

Combined depletion and electrostatic forces in polymer-induced membrane adhesion: a theoretical model

We develop a semi-quantitative analytical theory to describe adhesion between two identical planar charged surfaces embedded in a polymer-containing electrolyte solution. Polymer chains are uncharged and differ from the solvent by their lower dielectric permittivity. The solution mimics physiological fluids: It contains 0.1 M of monovalent ions and a small number of divalent cations that form tight bonds with the headgroups of charged lipids. The components have heterogeneous spatial distributions. The model was derived self-consistently by combining: (a) a Poisson-Boltzmann like equation for the charge densities, (b) a continuum mean-field theory for the polymer profile, (c) a solvation energy forcing the ions toward the polymer-poor regions, and (d) surface interactions of polymers and electrolytes. We validated the theory via extensive coarse-grained Molecular Dynamics (MD) simulations. The results confirm our analytical model and reveal interesting details not detected by the theory. At high surface charges, polymer chains are mainly excluded from the gap region, while the concentration of ions increases. The model shows a strong coupling between osmotic forces, surface potential and salting-out effects of the slightly polar polymer chains. It highlights some of the key differences in the behaviour of monomeric and polymeric mixed solvents and their responses to Coulomb interactions. Our main findings are: (a) the onset of long-ranged ion-induced polymer depletion force that increases with surface charge density and (b) a polymer-modified repulsive Coulomb force that increases with surface charge density. Overall, the system exhibits homeostatic behaviour, resulting in robustness against variations in the amount of charges. Applications and extensions of the model are briefly discussed.     

On the different roles of anions and cations in the solvation of enzymes in ionic liquids

The solvation of the enzyme Candida antarctica lipase B (CAL-B) was studied in eight different ionic liquids (ILs). The influence of enzyme-ion interactions on the solvation of CAL-B and the structure of the enzyme-IL interface are analyzed. CAL-B and ILs are described with molecular dynamics (MD) simulations in combination with an atomistic empirical force field. The considered cations are based on imidazolium or guanidinium that are paired with nitrate, tetrafluoroborate or hexafluorophosphate anions. The interactions of CAL-B with ILs are dominated by Coulomb interactions with anions, while the second largest contribution stems from van der Waals interactions with cations. The enzyme-ion interaction strength is determined by the ion size and the magnitude of the ion surface charge. The solvation of CAL-B in ILs is unfavorable compared to water because of large formation energies for the CAL-B solute cages in ILs. The internal energy in the IL and of CAL-B increases linearly with the enzyme-ion interaction strength. The average electrostatic potential on the surface of CAL-B is larger in ILs than in water, due to a weaker screening of charged enzyme residues. Ion densities increased moderately in the vicinity of charged residues and decreased close to non-polar residues. An aggregation of long alkyl chains close to non-polar regions and the active site entrance of CAL-B are observed in one IL that involved long non-polar decyl groups. In ILs that contain 1-butyl-3-methylimidazolium cations, the diffusion of one or two cations into the active site of CAL-B occurs during MD simulations. This suggests a possible obstruction of the active site in these ILs. Overall, the results indicate that small ions lead to a stronger electrostatic screening within the solvent and stronger interactions with the enzyme. Also a large ion surface charge, when more hydrophilic ions are used, increases enzyme-IL interactions. An increase of these interactions destabilizes the enzyme and impedes enzyme solvation due to an increase in solute cage formation energies.     

"Swiss-cheese" polyelectrolyte gels as media with extremely inhomogeneous distribution of charged species

"Swiss-cheese" polyelectrolyte gels (i.e., gels containing a regular set of closed spherical pores) are considered as a suitable system for modeling of a medium with extremely inhomogeneous distribution of charged species. It is shown that the inhomogeneous distribution of ions in Swiss-cheese polyelectrolyte gels can be reached simply by immersion of the gels in an aqueous solution of charged species (e.g., low-molecular 1-1 salt or multivalent ions and macroions charged likely to the gel chains). If a polymer gel is kept in such a solution for a long time, the concentration of ions within relatively big voids becomes equal to that in external solution. On the other hand, due to the Donnan effect the ion's concentration in polymer matrix is always lower than that in external solution. As a result the multivalent ions distribute between water voids and polymer matrix. The extent of this distribution is characterized by partition coefficient kD (determined as ratio kD = n(s)(void)/n(s)(mat) of the concentrations n(s)(void) and n(s)(mat) of ions in water voids and in polymer matrix, correspondingly). It is shown that the partition coefficient kD can be larger than 10 for low-molecular salt, reaches 10(3) for bivalent ions, and is higher than 10(6) for tetravalent ions. In the case of polymer macroions the partition coefficient kD tends to infinity. Our calculations show that the lower limit of characteristic scales of heterogeneity (determined by water voids size starting from which the condition of total electroneutrality is fulfilled and effect of partition is the most pronounced) can be equal to tens of nanometers.