Ion-pair Configurations Research Articles

The Many-Body Expansion for Aqueous Systems Revisited: IV. Stabilization of Halide-Anion Pairs in Small Water Clusters.

We report the structures, energetics, many-body effects, and vibrational spectra of water clusters stabilizing pairs of halide-anions, X-(H2O)kY- (k = 2-6, X/Y = F, Cl, Br, I) as well as their stability in the gas phase relative to fragmentation. We find that these metastable cluster structures mimicking the solvent-separated ion pair (SSIP) configurations in aqueous solutions are less stable relative to fragmentation into smaller ionic aqueous clusters containing a single halide anion. The many-body expansion (MBE) at these geometries was found to converge at the 4-body term, which is, however, significant, amounting to >20% of the total binding energy in several instances. The binding motif of these ion pair aqueous clusters starts as networks in which the water molecules form a "bridge" between the two halide-anions. As the cluster grows, these structures become destabilized by a more repulsive 3-body term for distances R(O-O) < 3.45 Å with respect to networks in which water molecules move outside the bridge, solvating the other side of the anion.

Structures and bonding characteristics of KCl(H2O)n clusters with n = 1-10 based on density functional theory.

Aqueous inorganic salt solutions play a prominent role in both physiological and chemical experiments, and significant attention has been directed toward understanding the mechanisms underlying salt dissolution. In our effort to elucidate the hydration process of potassium chloride, we employed a comprehensive genetic algorithm to explore the structures of KCl(H2O)n (n = 1-10). A series of stable structures were identified by high-level abinitio optimization and single-point energy calculations with a zero-point energy correction. An analysis of the probability distribution of KCl(H2O)1-10 revealed that clusters with high probability at low temperatures exhibit reduced probabilities at higher temperatures, while others become more prevalent. When n = 1-9, the contact ion pair configurations or partially dissociated structures dominate in the system, and the probability distribution plot shows that the proportion of the solvent-separated ion pair (SSIP) structures of KCl(H2O)n is very small, while the SSIP configuration in KCl(H2O)10 becomes a stable structure with increasing temperature. The results from natural bond orbital analysis reveal a stronger interaction between chloride ions and water molecules. These findings provide valuable insights for a more comprehensive understanding of the intricacies of potassium chloride dissolution in water.

Collective Proton Transfers in Cyclic Water–Ammonia Tetramers: A Path Integral Machine-Learning Study

We present results from machine-learning-based path integral molecular dynamics simulations that describe isomerization paths articulated via collective proton transfers along mixed, cyclic tetramers combining water and ammonia at cryogenic conditions. The net result of such isomerizations is a reverse of the chirality of the global hydrogen-bonding architecture along the different cyclic moieties. In monocomponent tetramers, the classical free energy profiles associated with these isomerizations present the usual symmetric double-well characteristics whereas the reactive paths exhibit full concertedness among the different intermolecular transfer processes. Contrastingly, in mixed water/ammonia tetramers, the incorporation of a second component introduces imbalances in the strengths of the different hydrogen bonds leading to a partial loss of concertedness, most notably at the vicinity of the transition state. As such, the highest and lowest degrees of progression are registered along OH···N and O···HN coordinations, respectively. These characteristics lead to polarized transition state scenarios akin to solvent-separated ion-pair configurations. The explicit incorporation of nuclear quantum effects promotes drastic depletions in the activation free energies and modifications in the overall shape of the profiles which include central plateau-like stages, indicating the prevalence of deep tunneling regimes. On the other hand, the quantum treatment of the nuclei partially restores the degree of concertedness among the evolutions of the individual transfers.

Adsorption and Oxidation of NO2 on Anatase TiO2: Concerted Nitrate Interaction and Photon-Stimulated Reaction

The catalytic and photon-induced oxidation of NO2 on anatase TiO2 has been studied and compared with the surface nitrate species obtained after adsorption of HNO3. Using a combination of in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), density functional theory (DFT), and temperature-programmed desorption (TPD), it is shown that identical products are obtained in all reaction systems but that their formation rates differ significantly. The surface reaction products are identified as combinations of surface–NO3– species, where NO2 bonds to the lattice oxygen, and freely adsorbed NO3– ions. These products can be obtained either by dissociative adsorption of HNO3 or by catalytic/photocatalytic oxidation of NO2, which is facilitated by UV light. A concerted reaction mechanism is unraveled that involves reorientation of bidentate nitrate that pushes out a neighboring protonated lattice oxygen to form a surface–NO3– species and a terminal OH group. The thermal stability of these surface species has been studied by means of TPD and simultaneous in situ DRIFTS measurements that reveal a main desorption peak (m/z = 46) at around 430 °C, which is attributed to concerted nitrate desorption through pentoxide (N2O5) formation. A weaker and broader TPD peak is found at about 185 °C and is attributed to desorption of nitrate species bonded in a compressed configuration. The experimental results can be explained by the changing stability of the identified nitrate products, which depends strongly on the surface chemical environment and the surface coverage. The DFT results show that the stabilization of intermediate NO2 adsorbates and the final nitrate reaction products occurs through a bifunctional charge exchange mechanism that is mediated by the TiO2 crystal. In particular, a stable surface–NO3– and NO3– ion pair configuration has been identified. This mechanism explains both the thermal and photoinduced oxidation of NO2 and their thermal stability and different formation rates, yielding high photoinduced oxidation reaction rates. Our results provide insights into the structure and chemical stability of nitrate surface products on TiO2 particles and their formation mechanism, which is important for understanding their catalyzed transformation into the harmful compounds HONO and N2O during continued UV light illumination.

Solvation and Ion-Pairing Effects of Choline Acetate Electrolyte in Protic and Aprotic Solvents Studied by NMR Titrations.

Choline‐based electrolytes have been proposed as environmentally friendly and low‐cost alternatives for secondary zinc air batteries. Choline acetate [Ch]+[OAc]− in protic (D2O) and aprotic (DMSO‐d6) solvents has been studied by means of concentration‐dependent 1H NMR, viscosity, and density measurements. The viscosities have been calculated on the basis of the Jones‐Dole equation and showed that the dominant contribution originates from short‐range ion‐solvent interactions. Site‐specific association affinities were assigned from NMR chemical shift titrations. In DMSO‐d6, the hydroxyl group of choline was found to have the smallest dissociation constant followed by the methyl group of acetate. The corresponding Gibbs energies at low concentration were found to be in agreement with a solvent‐separated ion pair (2SIP) configuration, whereas at concentrations above 300 mM, a solvent‐shared ion pair (SIP) configuration was assigned. For [Ch]+[OAc]− in D2O, association effects were found to be weaker, attributed to the high dielectric constant of the solvent. On time scales on the order of 100 ms, NMR linewidth perturbations indicated a change in the local rotational dynamics of the ions, attributed to short‐range cation‐solvent interactions and not to solvent viscosity. At 184 mM, ∼ 40 % of the cations in DMSO‐d6 and ∼ 10 % in D2O were found to exhibit short‐range interactions, as indicated by the linewidth perturbations. It was found that at about 300 mM, the ions in DMSO‐d6 exhibit a transition from free to collective translational dynamics on time scales on the order of 400 ms. In DMSO‐d6, both ions were found to be almost equally solvated, whereas in D2O solvation of acetate was stronger, as indicated by the obtained effective hydrodynamic radii. For [Ch]+[OAc]− in DMSO‐d6, the results suggest a solvent‐shared ion association with weak H‐bonding interactions for concentrations between 0.3–1 M. Overall, the extent of ion association in solvents such as DMSO is not expected to significantly limit charge transport and hinder the performance of choline‐based electrolytes.

Insights into Spontaneous Solid Electrolyte Interphase Formation at Magnesium Metal Anode Surface from Ab Initio Molecular Dynamics Simulations.

Spontaneous chemical reactivity at multivalent (Mg, Ca, Zn, Al) electrode surfaces is critical to solid electrolyte interphase (SEI) formation, and hence, directly affects the longevity of batteries. Here, we report an investigation of the reactivity of 0.5 M Mg(TFSI)2 in 1,2-dimethoxyethane (DME) solvent at a Mg(0001) surface using ab initio molecular dynamics (AIMD) simulations and detailed Bader charge analysis. Based on the simulations, the initial degradation reactions of the electrolyte strongly depend on the structure of the Mg(TFSI)2 species near the anode surface. At the surface, the dissociation of Mg(TFSI)2 species occurs via cleavage of the N-S bond for the solvent separated ion pair (SSIP) and via cleavage of the C-S bond for the contact ion pair (CIP) configuration. In the case of the CIP, both TFSI anions undergo spontaneous bond dissociation reactions to form atomic O, C, S, F, and N species adsorbed on the surface of the Mg anode. These products indicate that the initial SEI layer formed on the surface of the pristine Mg anode consists of a complex mixture of multiple components such as oxides, carbides, sulfides, fluorides, and nitrides. We believe that the atomic-level insights gained from these simulations will lay the groundwork for the rational design of tailored and functional interphases that are critical for the success of multivalent battery technology.

Benzene Excimer and Excited Multimers: Electronic Character, Interaction Nature, and Aromaticity.

In this Letter we analyze the forces involved in the formation of the benzene excimer and its electron structure, and (anti)aromatic character. We extend our study to excited states in molecular aggregates, the triplet excimer and the benzene-tricyanobenzene exciplex. Electronic wave functions are decomposed in terms of localized excitations and ion-pair configurations through diabatization, and we show that excimer (anti)aromaticity can be described as the linear combination of ground, excited, and ionic molecular states. Our analysis concludes that the benzene excimer must be characterized as antiaromatic, with weaker antiaromaticity than the molecular excited singlet. Moreover, we define a model electronic Hamiltonian for the excimer state and we use it as a building block for the extrapolation of electronic Hamiltonians in molecular aggregates. Benzene multimers present a nonuniform (anti)aromatic character, with the center of the column being antiaromatic and the edges behaving as aromatic. The implications of this work go beyond the study of the excimer, providing a general framework for the calculation and characterization of excited states in aggregates.

Size-Dependent Onset of Nitric Acid Dissociation in Cs+·(HNO3)(H2O)n=0-11 Clusters at 20 K.

We report the water-mediated charge separation of nitric acid upon incorporation into size-selected Cs+·(HNO3)(H2O)n=0-11 clusters at 20 K. Dramatic spectral changes are observed in the n = 7-9 range that are traced to the formation of many isomeric structures associated with intermediate transfer of the acidic proton to the water network. This transfer is complete by n = 10, which exhibits much simpler vibrational band patterns consistent with those expected for a tricoordinated hydronium ion (the Eigen motif) along with the NO stretching bands predicted for a hydrated NO3- anion that is directly complexed to the Cs+ cation. Theoretical analysis of the n = 10 spectrum indicates that the dissociated ions adopt a solvent-separated ion-pair configuration such that the Cs+ and H3O+ cations flank the NO3- anion in a microhydrated salt bridge. This charge separation motif is evidently assisted by the electrostatic stabilization of the product NO3-/H3O+ ion pair by the proximal metal ion.

Intra- and inter-molecular interactions in choline-based ionic liquids studied by 1D and 2D NMR

Choline acetate [Ch]+[OAc]− and choline trifluoroacetate [Ch]+[TFA]− mixed with 30 wt% H2O have been characterized by variable temperature 1D and 2D 1H/19F NMR. It was found that the cation, anion and water molecules participate in a strong local H-bond network, but there was no evidence for the formation of large clusters of ions. Weaker association effects were observed for the trifluoroacetate than for the acetate, originating from the higher anion acidity of TFA. Rotational correlation times, which were extracted by temperature dependent spin-lattice relaxation measurements, indicated that cations, anions and H2O exhibit collective reorientation even at room temperature. At the same time, no significant long-range correlations in the translational motion of the ions could be observed from PFG diffusion measurements. Moreover, 1D and 2D NOE experiments showed that both ions share close proximity to water molecules for a period on the order of 0.2 s at room temperature. It is deduced that the ions exist in a solvent-shared ion pair (SIP) configuration, with strongly correlated short-range dynamics between the choline cation and the H2O molecules, indicating a similar behavior of these ILs to a deep eutectic solvent with water as Lewis acid. A reason for this observation, which is supported by a high entropy term in the activation process particularly at low temperatures, could be a fairly static configuration with water near the –N(CH3)3 center of the choline, and relaxation caused by pseudorotation due to reordering of the H-bond network.

Optimizing the Energetics of FLP-Type H2 Activation by Modulating the Electronic and Structural Properties of the Lewis Acids: A DFT Study.

The great potential of frustrated Lewis pairs (FLPs) as metal-free catalysts for activation of molecular hydrogen has attracted increasing interest as an alternative to transition-metal catalysts. However, the complexity of FLP systems, involving the simultaneous interaction of three molecules, impedes a detailed understanding of the activation mechanism and the individual roles of the Lewis acid (LA) and Lewis base (LB). In the present work, using density functional theory (DFT) calculations, we examine the reactivity of 75 FLPs for the H2 splitting reaction, including a series of experimentally investigated LAs combined with conventional phosphine-based (tBu3P) and oxygen-based (i.e., ethereal solvent) Lewis bases. We find that the catalytic activity of the FLP is the result of a delicate balance of the LA and LB strengths and their bulkiness. The H2 splitting reaction can be changed from endergonic to exergonic by tuning the electrophilicity of the LA. Also, a more nucleophilic LB results in a more stable ion pair product and a lower barrier for the hydrogen splitting. The bulkiness of the LB leads to an early transition state to reduce steric hindrance and lower the barrier height. The bulkiness of the fragments determines the cavity size in the FLP complex, and a large cavity allows for a larger charge separation in the ion pair configuration. A shorter proton–hydride distance in this product complex correlates with a stronger attraction between the fragments, which forms more reactive ion pairs and facilitates the proton and hydride donations in the subsequent hydrogenation process. These insights may help with rationalizing the experimentally observed reactivities of FLPs and with designing better FLP systems for hydrogenation catalysis and hydrogen storage.

Impact of Size and Electronegativity of Halide Anions on Hydrogen Bonds and Properties of 1-Ethyl-3-methylimidazolium-Based Ionic Liquids.

The effect of the anion size and electronegativity of halide-based anions (Cl-, Br-, I-, and BF4-) on the interionic interaction in 1-ethyl-3-methylimidazolium-based ionic liquids (ILs) C2mim X (X = Cl, Br, I, and BF4) is studied by a combined approach of experiments (Raman, IR, UV-vis spectroscopy) and quantum chemical calculations. The fingerprint region of the Raman spectra of these C2mim X ion-pairs provides evidence of the presence of the conformational isomerism in the alkyl chain of the C2mim+ cation. The Raman and IR bands of the imidazolium C2-H stretch vibration for C2mim X (X = Cl, Br, I, and BF4) were noticeably blue-shifted with the systematic change in size of anions and the electronegativity. The observed blue shift in the C2-H stretch vibration follows the order C2mim BF4 > C2mim I > C2mim Br > C2mim Cl, which essentially indicates the strong hydrogen bonding in the C2mim Cl ion-pair. DFT calculations predict at least four configurations for the cation-anion interaction. On the basis of relative optimized energies and basis-set-superposition-error (BSSE) corrected binding energies for all ion-pair configurations, the most active site for the anion interaction was found at the C2H position of the cation. Besides information about the C2H position, our DFT results give insights into the anion interaction with the ethyl and methyl chain of the cation, which was also confirmed experimentally [ Chem. Commun. 2015 , 51 , 3193 ]. The anion interaction at the C2H site of the cation favors a planar geometry in C2mim X for X = Cl, Br, and I; however, for BF4, the system prefers a nonplanar geometry where the anion is located over the imidazolium ring. TD-DFT results were used to analyze the observed UV-vis absorption spectra in a more adequate way giving insights into the electronic structure of the ILs. Overall, a reasonable correlation between the observed and the DFT-predicted results is established.

Direct Observation of Ion Pairing in Aqueous Nitric Acid Using 2D Infrared Spectroscopy.

Ion-ion interactions and ion pairing play an important role in the properties of concentrated electrolyte solutions, yet remain difficult to study due to the heterogeneous and highly dynamic behavior of these systems. In concentrated acid solutions, these questions take on a further level of complexity because the structure of the aqueous proton itself is uncertain and may be influenced by the counterion. Here, we address these questions by studying the IR spectra of nitric acid as a function of concentration in H2O and comparing these to the spectra of several alkali nitrate salts. We show how the close proximity between cations and NO3- ions in solution at high concentration affect the IR spectra and therefore the molecular structures. Using two-dimensional IR spectroscopy, we demonstrate the formation of contracted ion pair configurations in nitric acid solutions between NO3- ions and H+(aq) via the observation of a distinct anisotropic intermolecular crosspeak between these species. By studying the concentration dependence of this spectral feature, we show that this ion-paired configuration exists in solution at concentrations as low as 2 M and suggests that the structure of H+(aq) solvation complex in these ion pairs differs from the structure in bulk solution.

Effect of Salt Concentration on Properties of Lithium Ion Battery Electrolytes: A Molecular Dynamics Study

Electrolyte solutions of 1 M concentration are typically used in lithium ion batteries (LIB) for optimal performance. However, recently, superconcentrated electrolytes have been proposed to be a promising alternative to 1 M solutions. Despite their improved stability features, application of the concentrated electrolytes is hindered by their poor transport properties. We probe EC-LiPF6 electrolyte system for a range of concentrations: 0.06 to 4 M using molecular dynamics simulations to study the effect of concentration on transport and structural properties. Molecular structure of the solution changes with concentration from a predominantly solvent separated ion pair (SSIP) configuration at the dilute limit to an aggregate rich configuration at high concentration. Depletion of SSIPs and formation of more aggregates at higher concentrations affect the transport properties. The present work provides insights into the relation between molecular structure and performance of the electrolyte solution and sugge...

Rate theory of ion pairing at the water liquid-vapor interface: A case of sodium iodide.

Studies on ion pairing at interfaces have been intensified recently because of their importance in many chemical reactive phenomena, such as ion-ion interactions that are affected by interfaces and their influence on kinetic processes. In this study, we performed simulations to examine the thermodynamics and kinetics of small polarizable sodium iodide ions in the bulk and near the water liquid-vapor interface. Using classical transition state theory, we calculated the dissociation rates and corrected them with transmission coefficients obtained from the reactive flux formalism and Grote-Hynes theory. Our results show that in addition to affecting the free energy of ions in solution, the interfacial environments significantly influence the kinetics of ion pairing. The results on the relaxation time obtained using the reactive flux formalism and Grote-Hynes theory present an unequivocal picture that the interface suppresses ion dissociation. The effects of the use of molecular models on the ion interactions as well as the ion-pair configurations at the interface are also quantified and discussed.

Role of the Counteranion in the Reaction Mechanism of Propylene Polymerization Catalyzed by a (Pyridylamido)hafnium(IV) Complex

We studied the stereoselective propylene polymerization catalyzed by the ion pair (IP) active species of (pyridylamido)hafnium(IV) complex, the catalyst. In particular, we focused on the role of the counteranion (CA) for the reaction mechanism of propylene insertion because it had never been investigated from the molecular point of view. First, we searched for a variety of IP configurations by the molecular dynamics (MD) method and quantum chemical calculation, and located the transition states (TSs) of propylene insertion. By comparing the CA-bound catalysts with the isolated ones, it was revealed that the manner of energetically favorable insertion is interestingly altered due to the location of the CA. Next, we showed that the CA is essential to determine the coordination structures of the Hf atom at transition state. Therefore, it was concluded that the role of the CA is critical to accurately control the polymerization mechanism. Finally, we proposed a synthesis experiment which would be advantageous...

Solvation structures of sodium halides in dimethyl sulfoxide (DMSO)–methanol (MeOH) mixtures

A constrained molecular dynamics technique has been used to study the structures and dynamics of the solvation shells of three sodium halides, namely sodium chloride (Na+–Cl−), sodium bromide (Na+–Br−) and sodium iodide (Na+–I−) in DMSO–MeOH mixtures. In the case of Na+–Cl− and Na+–Br−, Na+ is preferentially solvated by DMSO and Cl− and Br− are preferentially solvated by methanol in the contact ion pair (CIP) state. In the solvent-assisted ion pair (SAIP) configuration, Na+ ions of Na+–Cl− and Na+–Br− are preferentially solvated by methanol and Cl− and Br− also show preferential solvation by methanol over DMSO. In the case of Na+–I−, the only preferential solvation is in the SAIP state for I− ion by methanol. These observations are supported by the calculated excess coordination numbers and spatial density maps. The heights of the transition states barriers for CIPs and SAIPs/solvent-shared ion pairs (SSHIPs) are significantly affected when the mole fraction of methanol (xMeOH) changes from 0.0 to 0.25 because of a significant increase in the methanol density around halides. From the analysis of angular distribution functions of DMSO and methanol around the cations and anions, it is seen that DMSO and methanol molecules are present in parallel dipolar orientations (with respect to cation–solvent vector) in the first coordination shell of these three ion pairs at the CIP and SAIP states. Methanol molecules are nearly in an antiparallel (with respect to ion–solvent vector) orientation around the three halide ions.

Solvent-shared pairs of densely charged ions induce intense but short-range supra-additive slowdown of water rotation.

The question "Can ions exert supra-additive effects on water dynamics?" has had several opposing answers from both simulation and experiment. We address this ongoing controversy by investigating water reorientation in aqueous solutions of two salts with large (magnesium sulfate) and small (cesium chloride) effects on water dynamics using molecular dynamics simulations and classical, polarizable models. The salt models are reparameterized to reproduce properties of both dilute and concentrated solutions. We demonstrate that water rotation in concentrated MgSO4 solutions is unexpectedly slow, in agreement with experiment, and that the slowdown is supra-additive: the observed slowdown is larger than that predicted by assuming that the resultant of the extra forces induced by the ions on the rotating water molecules tilts the free energy landscape associated with water rotation. Supra-additive slow down is very intense but short-range, and is strongly ion-specific: in contrast to the long-range picture initially proposed based on experiment, we find that intense supra-additivity is limited to water molecules directly bridging two ions in solvent-shared ion pair configuration; in contrast to a non-ion-specific origin to supra-additive effects proposed from simulations, we find that the magnitude of supra-additive slowdown strongly depends on the identity of the cations and anions. Supra-additive slowdown of water dynamics requires long-lived solvent-shared ion pairs; long-lived ion pairs should be typical for salts of multivalent ions. We discuss the origin of the apparent disagreement between the various studies on this topic and show that the short-range cooperative slowdown scenario proposed here resolves the existing controversy.

Designing New Conductive Imidazole-Derived Salt for Lithium-Ion Battery Electrolytes: A Computational Approach

Ever since the commercialization of the rechargeable Li-ion cell, LiPF6 is by far the most widely used electrolyte salt. It can produce electrolyte solutions with a superior ion conductivity, good electrochemical oxidation stability, and adequate thermal stability when dissolved in binary or ternary solvent systems of cyclic carbonates. However, LiPF6is also reported to have series problems that may lead to performance and safety issues (1). In order to circumvent these issues, researches have focused on designing new lithium salts which have thermal and electrochemical stability, high ionic conductivity and easy synthesis route (2-6). In this work, we report designing of an alternative salt derived from bis(trifluoroborane)imidazolide (LiIm(BF3)2) by employing density functional theory (DFT) calculations. New anion structures were generated by considering substitution at the C5 position with electron withdrawing groups (-F, -CF3, -SO3H and -NO2) and electron donating groups (-NH2, -OCH3 and –CH3). By using the optimized anion structures, random lithium ion-pair configurations (around 130 structures for each) were generated and optimized to identify the most stable configurations. The effects of different substituents were investigated with respect to ion dissociation energies, anion oxidative stability, and aromaticity of the molecules. Based on our calculations, we have found that, ion dissociation energies and intrinsic oxidative potentials of the anions were largely affected by the type of substituents introduced on the parent structure. Compared to –PF6 and -Im(BF3)2 (standard anions), –CF3 and -NO2 substituted anions showed an increase in anion oxidation stability. Since, BF3 occupies the strong Lewis basic sites of the imidazole ring, the parent structure itself exhibit a charge delocalization over the aromatic ring, which results in a weak coordinating ability with Li+. With the exception of –F, -CF3 and –CH3, further introduction of substituents at C5 position led to a much lower cation–anion interaction. Among the possible anions studied, 5-nitro-bis(trifluoroborane)imidazolide, with oxidation potentials > 5.7 V vs. Li+/Li , can be considered as a good candidates for high voltage Li-ion battery applications. Moreover, we have also found that the ion-pair dissociation energy, and anion stability of the selected salt are much improved compared to experimentally reported LiIm(BF3)2 and LiPF6.

The coupling between stability and ion pair formation in magnesium electrolytes from first-principles quantum mechanics and classical molecular dynamics.

In this work we uncover a novel effect between concentration dependent ion pair formation and anion stability at reducing potentials, e.g., at the metal anode. Through comprehensive calculations using both first-principles as well as well-benchmarked classical molecular dynamics over a matrix of electrolytes, covering solvents and salt anions with a broad range in chemistry, we elucidate systematic correlations between molecular level interactions and composite electrolyte properties, such as electrochemical stability, solvation structure, and dynamics. We find that Mg electrolytes are highly prone to ion pair formation, even at modest concentrations, for a wide range of solvents with different dielectric constants, which have implications for dynamics as well as charge transfer. Specifically, we observe that, at Mg metal potentials, the ion pair undergoes partial reduction at the Mg cation center (Mg(2+) → Mg(+)), which competes with the charge transfer mechanism and can activate the anion to render it susceptible to decomposition. Specifically, TFSI(-) exhibits a significant bond weakening while paired with the transient, partially reduced Mg(+). In contrast, BH4(-) and BF4(-) are shown to be chemically stable in a reduced ion pair configuration. Furthermore, we observe that higher order glymes as well as DMSO improve the solubility of Mg salts, but only the longer glyme chains reduce the dynamics of the ions in solution. This information provides critical design metrics for future electrolytes as it elucidates a close connection between bulk solvation and cathodic stability as well as the dynamics of the salt.

Dissecting ion-specific dielectric spectra of sodium-halide solutions into solvation water and ionic contributions.

Using extensive equilibrium molecular dynamics simulations we determine the dielectric spectra of aqueous solutions of NaF, NaCl, NaBr, and NaI. The ion-specific and concentration-dependent shifts of the static dielectric constants and the dielectric relaxation times match experimental results very well, which serves as a validation of the classical and non-polarizable ionic force fields used. The purely ionic contribution to the dielectric response is negligible, but determines the conductivity of the salt solutions. The ion-water cross correlation contribution is negative and reduces the total dielectric response by about 5%-10% for 1 M solutions. The dominating water dielectric response is decomposed into different water solvation shells and ion-pair configurations, by this the spectral blue shift and the dielectric decrement of salt solutions with increasing salt concentration is demonstrated to be primarily caused by first-solvation shell water. With rising salt concentration the simulated spectra show more pronounced deviations from a single-Debye form and can be well described by a Cole-Cole fit, in quantitative agreement with experiments. Our spectral decomposition into ionic and different water solvation shell contributions does not render the individual contributions more Debye-like, this suggests the non-Debye-like character of the dielectric spectra of salt solutions not to be due to the superposition of different elementary relaxation processes with different relaxation times. Rather, the non-Debye-like character is likely to be an inherent spectral signature of solvation water around ions.