Properties of Ion Complexes and Their Impact on Charge Transport in Organic Solvent-Based Electrolyte Solutions for Lithium Batteries: Insights from a Theoretical Perspective

: Electrolyte formulations in standard lithium ion and lithium metal batteries are complex mixtures of various components. In this article, we review molecular key principles of ion complexes in multicomponent electrolyte solutions in regards of their inﬂuence on charge transport mechanisms. We outline basic concepts for the description of ion–solvent and ion–ion interactions, which can be used to rationalize recent experimental and numerical ﬁndings concerning modern electrolyte formulations. Furthermore, we discuss beneﬁts and drawbacks of empirical concepts in comparison to molecular theories of solution for a more reﬁned understanding of ion behavior in organic solvents. The outcomes of our discussion provide a rational for beneﬁcial properties of ions, solvent, co-solvent and additive molecules, and highlight possible routes for further improvement of novel electrolyte solutions.


Introduction
The use of multicomponent electrolyte solutions in modern electrochemical storage devices such as rechargeable lithium ion and lithium metal batteries (LIBs and LMBs) is of fundamental importance for effective ion shuttling and transport mechanisms [1][2][3][4][5][6][7][8][9][10]. Nowadays well-known fast charging rates as well as the impressive performance of recent LIB and LMB devices [11][12][13][14] are not possible without decades of intense research on advanced electrolyte compositions. Although often not known in detail, standard electrolyte formulations in commercial energy storage devices mostly include various organic solvents, a plethora of distinct lithium ion conducting salts, and a broad range of possible further components in order to fulfill several duties related to the performance and safety of the electrochemical cell [1,3,6]. In addition to ions and the main solvent, components with concentrations at the molar scale are usually called co-solvents, whereas all other species which have total mass ratios lower than 5 wt % are classified as additives [15]. It is worth noting that additives are most often used as flame retardants, for the formation of protective layers at the electrodes, or as ion solvation enhancers [15][16][17]. Prominent examples for standard electrolyte formulations are high dielectric constant solvents of cyclic alkyl carbonates such as propylene carbonate (PC) or ethylene carbonate (EC), which provide an enhanced dissolution of the lithium salts, in combination with linear esters such as dimethyl carbonate (DMC), or ethylmethyl carbonate (EMC) as co-solvents to decrease the viscosity of the solution [1,3,6,18,19]. In combination with one or more lithium ion conducting salts, it becomes clear that the underlying ion correlation and charge transport mechanisms in the solution reveal a high level of complexity. Despite important benefits in terms of relatively low cost, most components of liquid organic electrolyte formulations are highly flammable and their use is restricted to well-defined electrode voltage ranges [5][6][7]16]. Recent research thus focuses on novel electrolyte formulations with lower flammability, higher electrochemical and thermal stability, and the presence of multi-functional molecular groups in regards of improved efficiencies and safeties for next-generation LIBs, LMBs and dual ion batteries (DIBs) [1,3,[6][7][8]16,18,[20][21][22][23][24][25][26][27][28][29][30][31][32].
Whereas most electrolyte solutions are well characterized in terms of their thermophysical and electrochemical properties [1,3], less is known on the underlying molecular interactions between the components. Due to the complexity of the solution, previous theoretical approaches such as scaled particle theory [33] mostly focus on ad hoc molecular approximations, which are often too simple and not applicable for distinct species in non-ideal or non-homogeneous solutions [34]. Consequently, parts of the general research strategy for the development of improved electrolyte solutions rely on empirical concepts and simple rules of thumbs instead of a rigorous molecular understanding of interaction mechanisms. A prominent example is the often cited alchemical law similia similibus solvuntur, which means the same dissolves in the same. In regards of this approach, it is often assumed that polar solutes and specifically lithium salts are highly solvated by polar solvents, whereas they are nearly insoluble in apolar media. It is worth noting that broad discrepancies between the individual saturation concentrations for distinct lithium salts in various organic solvents highlight crucial limitations for this simple approach [8,28,30,[35][36][37][38][39][40].
Hence, the development of novel electrolyte formulations often relies on limited molecular understanding, empirical arguments, theoretical approaches suffering from unspecific approximations, or previous experimental findings instead of a desirable frame-guided design of tailor-made components for distinct purposes. In this context, refined theoretical and numerical approaches often provide useful frameworks to establish a more advanced level of understanding. Although molecular theories of solution usually suffer from an absence of exact analytic results, the benefits of the Kirkwood-Buff (KB) theory [34,41], as such an approach for the study of ion correlation effects, solubilities of salts, development of ion models, or the influence of co-solvents on ion dissociation equilibria, were recently pointed out [30,[42][43][44][45][46][47][48][49][50][51][52].
In this article, we introduce basic theoretical concepts, which allow us to rationalize specific properties of ion complexes and electrolyte solutions in more detail. We discuss benefits and limitations of theoretical models, and review their applicability in comparison to experimental and simulation outcomes. The corresponding findings highlight the importance of specific ion effects, ion-ion correlations, molecular effects of solvation, and the influence of additive and co-solvent molecules on the dissociation behavior of the ions. Our discussion provides important insights into the properties of ions in solution, which allow us to answer relevant questions, and to avoid pitfalls of modern electrolyte research in terms of long-standing, but questionable assumptions.
In particular, we focus on the following questions: 1. Which factors determine ion-ion interactions, and how do the corresponding effects modify charge transport? 2. Which properties of solvents are appropriate discriminators in order to distinguish between good and poor solvents, and which solvents are well-suited to increase the performance of electrochemical cells in terms of high salt solubility and beneficial charge transport behavior? 3. How does the presence of co-solvents or additives affect the solvation of ions, and which conclusions can be drawn for the solubility of salts?
It has to be pointed out that most of the further important effects for electrolyte solutions, for instance chemical reactions of components at electrode interfaces, as well as the influence of strong electric fields in combination with polarization effects, are outside the scope of this review. With regard to these points, we refer the reader to further interesting and recent works besides various other studies [1,3,6,8,9,16,[53][54][55][56][57][58][59].
Our article is organized as follows. With regard to Figure 1, we introduce three main sections, which focus on ion-ion interactions and the resulting consequences for charge transport, the interaction between ions and solvents, and the properties of ions in complex mixtures such as multicomponent solutions. The article starts with a basic introduction into ion correlation effects in terms of simple mean-field approaches. Hereafter, we discuss the properties of ion complex states in organic solvents and their influence on the charge transport. In Section 3, we discuss potential influences on the ion solvation behavior. In more detail, we review beneficial concepts of donor and acceptor numbers in combination with softness and hardness parameters, which provide a rational for the the occurrence of specific ion effects. The influence of co-solvents or additives on ion correlation effects is discussed in Section 4. Herein, we introduce basic concepts from KB theory in order to derive explicit expressions for the ion solubility, and the change of the chemical equilibrium between dissociated and associated ion states in presence of multicomponent solutions. At the end, we combine the concepts in order to rationalize the properties of recent electrolyte formulations, and to discuss possible routes for future research. We briefly conclude and summarize in the last section.

Ions in Solution: Correlation Effects and Their Influence on Charge Transport
In this section, we review basic principles regarding the properties of ions in organic solvent-based electrolyte solutions. In more detail, we discuss general concepts and highlight the crucial role of specific ion effects with regard to the underlying solvation behavior. Moreover, we introduce simple electrostatic mean-field theories and comment on their deficiencies for the description of real electrolyte solutions. At this point, we ignore the molecular details of the solvent, and thus rely on a continuum background model with a global and constant dielectric constant. The corresponding simplifications allow us to point out basic properties of ion complexes in solution. The abundant occurrence of long-range electrostatic correlation effects between the ions and the corresponding consequences for charge transport are discussed at the end of this section.

Electrostatic Interactions and Properties of Ion Complexes
The formation of ion complexes in electrolyte solutions is governed by electrostatic interactions. For two ions in a solvent with dielectric constant r , the electrostatic Coulomb potential Φ(r) between the ions shows a decay on large length scales with [60] where r denotes the distance between the ions. Note that Equation (1) with Boltzmann constant k B and vacuum permittivity 0 . The density of ion species j in bulk phase at Φ(r) = 0 is denoted by ρ j . It is worth noting that, even within the mean-field PB approach, analytic solutions except for simple geometries are hard to derive [60]. In contrast, for simple and dilute ions under the condition of charge neutrality (∑ j z j ρ j = 0), unit valency and moderate maximum electrostatic potential Φ(r s ) = Φ s with Φ s /k B T 1 at the hydrodynamic boundary position r s around the main solute, Equation (2) can be linearized, which yields for the electrostatic potential around a spherical object [60] with the Debye-Hückel length suggesting a faster decay of electrostatic interactions at finite total ion density ρ s when compared to single ion pairs. The corresponding decay of Φ(r) due to the presence of the surrounding ions when compared to the standard Coulomb potential for single ions (Equation (1)) is depicted in Figure 2. As can be seen, even for moderate values of the Debye-Hückel length κ −1 D = σ, where σ denotes the radius of the charged object, a significant reduction of the electrostatic potential and thus a screening effect for r/σ 1 can be observed. However, also for dilute ions in the limit of vanishing values for κ D and thus in presence of non-modified Coulomb potentials, electrostatic interactions become screened beyond a characteristic distance, which is called Bjerrum length and defined by [61] which increases with decreasing temperature and dielectric constant of the solution. In more detail, the corresponding value of the Bjerrum length denotes the distance between the ions, where electrostatic interactions are of comparable magnitude as the thermal energy k B T. Consequently, even in the limit of vanishing salt concentrations, electrostatic interactions between the ions can be ignored for distances r ≥ λ B , as induced by a dielectric screening effect of the surrounding solvent molecules [63].  [1] have rather short Bjerrum lengths of λ PC B ≈ 0.84 nm and λ EC B ≈ 0.61 nm, respectively. As a simple estimate, the value of the Bjerrum length can be used to define the largest distance for electrostatically bound or correlated ions. It has to be noted that other distance criteria for ion pairs and ion complexes have also been introduced [54], but, due to reasons of simplicity, we rely on the Bjerrum length as fundamental length scale for the definition of bound ions or ion complexes without further restrictions. Previous simulation outcomes by some of us highlight the validity of this simple estimate [30,38,51,[64][65][66]. 1/r exp(-κ D r) 1/r Figure 2. Electrostatic potential Φ(r) between two differently charged species at distance r in presence of a finite salt density according to Equation (3) with κ −1 D = σ (red line), and for unscreened Coulomb interactions with Φ(r) ∼ 1/r (blue line). The electrostatic potential at the hydrodynamic boundary of the object r s = σ is arbitrarily chosen with Φ s = −k B T.
In more detail, one can observe contact (CIP), solvent-shared (1SP), and solvent-separated (2SP) ion pairs in accordance with Figure 3. Despite a vast amount of current research efforts, the mechanisms underlying ion complex formation are still not fully understood [53,69,[79][80][81].  In more detail, CIPs usually correspond to closest-distance contact ion pairs, whereas 1SP states include ion pairs, which are separated by one solvation shell, in contrast to 2SP states, which are separated by two or more solvation shells. Consequently, each ion in the 2SP state maintains its individual solvation shell, which is only marginally affected by the presence of the counterion. Consequently, the distance between the ions increases in the order CIP < 1SP < 2SP, which is in agreement with the reciprocal binding energies associated with the individual ion states. Despite the presence of various solvation shells [51,65], the distance between the ions in the 2SP state, as the largest complex in the aforementioned series, is always smaller than the corresponding value of λ B [70]. In terms of a moderate salt concentration of c s 1 mol/L, it has to be noted that most ion complexes in organic solvents are dominated by solvent-separated pairs [8,59]. The underlying reason can be attributed to the molecular details of the solvent molecules, whose coordination around the ions partly compensates the loss of electrostatic energies associated with the formation of CIP states.
At higher salt concentrations and for specific lithium salts such as lithium bis(trifluorosulfonyl)imide (LiTFSI) with bulky anions, the occurrence of ion aggregate (AGG) states in various solvents has also been observed [8,26,32,59,82]. A prominent example for AGG states are super-concentrated electrolytes or solvent-in-salt electrolyte solutions [8,26,40,59,82]. As the name "solvent-in-salt" implies, the number of solvent molecules is significantly reduced in these solutions, such that the individual solvation shells around the ions are only partially filled by solvent molecules [8,40,59]. Consequently, the number of ion complexes including CIP and AGG states increases with increasing salt concentration, although the corresponding detailed mechanisms and structures are yet under debate [82]. Over the last years, super-concentrated electrolytes for various solvents and salts [26,40,59,82] attracted enormous interest with regard to their potential use in LIBs and LMBs [8,59]. First results for acetonitrile, PC, sulfolane, and even water are promising in terms of reaching relatively high electrochemical stabilities and sufficient ionic conductivities [8,26,40,59,82].
Despite many numerical and experimental studies, a consistent theoretical framework to rationalize specific ion effects is still missing [75]. As a first empirical attempt towards this aim for aqueous electrolyte solutions, it has been proposed to distinguish ions in terms of their chaotropic and kosmotropic properties [53,54,67,68,70,73]. Hence, smaller ions such as F − and Li + are kosmotropes, meaning that they are water structure makers, whereas larger ions such as I − and SCN − are chaotropes, which emphasizes the fact that these species have to be considered as water structure breakers [53,54,63,70]. With regard to these definitions, it has to be clearly pointed out that kosmotropic and chaotropic properties are solely defined for ions in aqueous solutions, whereas the occurrence of comparable effects in aprotic solvents is currently under debate [30,76,83]. Chaotropic and kosmotropic effects on the water structure can be observed in terms of viscosity variations associated with the Jones-Dole viscosity B coefficient, as well as with regard to the dynamic behavior of water molecules in terms of dipolar relaxation times [53,63,87]. It is worth noting that recent studies on mixtures of salts with distinct properties reveal that the individual effects compensate each other [88].
As an extension of the chaotrope/kosmotrope approach, the law of matching water affinities (LMWA) was introduced by Collins [70,89], which rationalizes the experimental finding that pairs of kosmotropic or chaotropic ions reveal the highest stability. In this context, the meaning of kosmotropic and chaotropic behavior is no longer understood in terms of the impact on the water structure, rather than the degree of hydration for the corresponding ions [86]. Consequently, it is generally accepted that kosmotropes show a high degree of hydration, whereas the opposite can be assumed for chaotropic ions. The underlying empirical concept of the LMWA relies on calculated data regarding the differences of the individual ion hydration enthalpies W x + and W y − , the so-called heats of hydration ∆∆H hyd = W x + − W y − , whose difference is compared to the experimentally measurable standard heat of solution Q s of a crystalline salt with identical ions in infinite dilution. Corresponding results for distinct salts in terms of so-called volcano plots [70,83,89,90] are depicted in Figure 4A.
Hence, for Q s > 0, the separation of an ion pair is energetically unfavorable (endothermic reaction), whereas the ion pair easily breaks for Q s < 0 (exothermic reaction). With regard to the presented data, it becomes evident that ions with comparable size or influence on the water structure in terms of ∆∆H hyd ≈ 0 mostly reveal values of Q s > 0. In contrast, chaotropic and kosmotropic ions and thus ions of different size significantly differ in their individual heats of hydration according to |∆∆H hyd | 0, which is most often accomplished by unstable ion pairs in terms of Q s < 0. Consequently, ions of comparable size show comparable degrees of hydration, so-called water affinities, such that combinations of small anions and cations or combinations of large anions and cations form the most stable ion pairs ( Figure 4B [70]. For alkali halide salts, it thus can be seen that CsI with large and chaotropic ions forms the most stable ion pair. In contrast, LiI is the least stable ion pair, as it is composed of the kosmotropic lithium ion and the chaotropic iodide ion. The LMWA concept is discussed in more detail in Refs. [63,70,86,90], where it is also highlighted that electrostatic energies between the ions contribute a significant amount to the resulting values of Q s . Although the presence of chaotropic properties with regard to the hydration behavior is often debated, simulation findings and experimental outcomes reveal that the complexation behavior of the ions is strongly influenced by the hydration behavior [69]. Moreover, it has to be noted that the LMWA approach reveals some similarities with the famous hard/soft acids and bases (HSAB) principle, which states that hard Lewis acids as electron acceptors prefer to coordinate with hard Lewis bases as electron donors, and soft acids prefer to coordinate with soft bases [91]. As already briefly mentioned, few studies are devoted to the occurrence of specific ion effects in non-aqueous solutions [30,38,[76][77][78]83]. With regard to the corresponding results, it has been observed that the LMWA approach is also applicable for several organic solvents, such that it can be generalized to a law of matching solvent affinities [83]. Recent experimental results [78,83] thus reveal the applicability of this concept for various organic solvents such as methanol, formamide, PC, and dimethyl sulfoxide (DMSO).
In addition to the importance of these effects for lithium salts in LIBs, comparable conclusions can be drawn for sodium salts in sodium ion batteries (SIBs). In SIBs, the conducting ion is changed from Li + to Na + , which slightly modifies the molecular interactions between the ion species, or between the ions and the solvent molecules. In fact, the differences in the solvation behavior between Li + and Na + in most organic solvents are rather marginal, which can be rationalized by the low charge of both ions. Hence, it has been discussed that identical solvents or anions such as in LIBs can be used in SIBs [10]. Despite recent successful efforts with regard to new materials, it has to be mentioned that SIBs provide lower energy densities when compared to LIBs due to a higher standard electrode potential of Na + in contrast to Li + besides other drawbacks [92]. In contrast, the properties of salts with divalent or trivalent ions change significantly when compared to lithium or sodium salts. It has been reported that the transfer of multivalent ions from water to organic solvents is highly unfavorable, as evidenced by positive transfer free energies [93]. Furthermore, multivalent ions show a high preference being in contact with water molecules, which hinders systematic studies on ion pairing and ion association in hygroscopic organic solvents. In contrast to monovalent ions, the association tendency of multivalent ions is also significantly increased, which can be rationalized by stronger electrostatic interactions between the ion species. As a specific example, it has been shown that large objects reveal a charge reversal behavior in presence of multivalent ions, which means that more ions are bound than needed for charge compensation [94,95].

Ion Correlation Effects and Transport Behavior
With regard to the previous discussion in Section 2.2, we now demonstrate why the occurrence of ion complexes is of fundamental importance for electrolyte research. As prime example, the formation of ion complexes heavily affects the ionic conductivity due to omnipresent correlations between the ions. Strictly speaking, ions which form net neutral complexes, as mentioned before, do not contribute to the ionic conductivity, which negatively affects the power density of electrochemical cells. It is thus beneficial to achieve high values of ion dissociation, which is usually supported by using high dielectric constant solvents [1].
With regard to the influence of ion correlations on charge transport, the ideal ionic conductivity for non-interacting, i.e., ideal, ions is given by the Nernst-Einstein equation, which for a 1:1 salt in the limit of long times reads [96,97] where D s ++ and D s −− denote the individual self-diffusion coefficients for cations and anions. With regard to the presence of ion complexes, it becomes obvious that Equation (6) is only valid in the limit of low salt concentrations due to negligible ion correlation effects. Hence, for higher salt concentrations, the effective, or experimentally observed, ionic conductivity in terms of the Einstein-Helfand equation [51,65,66,96,98] reads with and where D −− and D ++ denote the cross-correlated diffusion coefficients between cations or anions, respectively, and D ± the corresponding symmetric diffusion coefficients including correlation effects between anions and cations. As can be seen, the contributions of D ± diminish the effective ionic conductivity considerably when compared to the ideal ionic conductivity. Consequently, the ideal Nernst-Einstein ionic conductivity suffers from a neglect of ionic correlations, such that Equation (6) provides the maximum possible ionic conductivity in accordance with σ id > σ eff . In addition, the calculation of ideal and effective ionic conductivities allows us to estimate the amount of correlated ion motion. As a specific expression for the dynamic correlation factor, the ratio of both ionic conductivities reads [96,98] with ∆ IC ∈ [0, 1], where a value of ∆ IC = 1 reveals full ion correlation, and vice versa ∆ IC = 0 the absence of any correlation effects. It has to be pointed out that the value of ∆ IC represents the fraction of correlated ion motion, which does not necessarily coincide with the fraction of ion complexes. As discussed in Section 2.2, all CIP, AGG, 1SP and 2SP states are examples of ion complexes and thus contribute to the actual value of the dynamic correlation factor. A further distinction between the individual contributions regarding the value of ∆ IC is hence not obvious without any detailed analysis of the structural arrangement.
In addition to the ionic conductivity, the cation transference number as a further descriptor of charge transport [9] with is also of fundamental importance for LIBs or LMBs, as it mainly determines the charging and discharging rates. For highly effective electrolyte solutions or single ion conductors, one can observe values of t + > 0.5 [9,97], whereas most values in organic liquid electrolyte solutions are t + < 0.5 [6,9]. The corresponding value reveals that anions instead of cations mostly contribute to the ionic conductivity as an unwanted side effect in organic solvent-based LIBs. In the limit of negligible ion correlations or low salt concentrations, the cation transference number coincides with the cation transport number lim ∆ IC →0 t + = t + app , which is defined by [96,99] and hence ignores any correlation effects. Experimental results reveal significant deviations between transference and transport numbers, which highlights the influence of correlation effects for charge transport mechanisms in external electric fields [99][100][101]. Consequently, the formation of ion complexes at high and moderate salt concentrations imposes a significant limitation on ionic conductivities and on transference numbers. As these values are of fundamental importance for the electrochemical performance of LIBs, it is thus desirable to minimize these effects considerably [9].
As an example for charge transport in adiponitrile, the values of ideal and effective ionic conductivities at 300 K by means of atomistic molecular dynamics (MD) simulations in comparison to experimental values for different concentrations of LiBF 4 and LiTFSI are shown in Figure 5 [30]. In agreement with the Walden relation [63], which reads σ id ∼ c s /η, where η is the shear viscosity of the solution following from insertion of the Stokes-Einstein equation with D s ∼ 1/η, one observes a linear increase of the ionic conductivity at low salt concentrations due to ∆ IC ≈ 0. At moderate concentrations, a maximum value of the ionic conductivity can be observed in terms of a balance between ion correlation effects and a high amount of charge carriers. In agreement with the limiting law for the ionic conductivity at very high salt concentrations [102,103] with σ = σ 0 − ξc 3/2 s , the values for the ionic conductivity reveal a nonlinear decay due to a positive Onsager prefactor ξ. With regard to the values for the effective Einstein-Helfand ionic conductivity, as shown in Figure 5, it can be concluded that the smaller BF − 4 anion induces a higher tendency of ion complex formation when compared to the TFSI − anion, which is in agreement with the LMWA and HSAB concept. Further analysis reveals that both anions have comparable values for self-diffusion constants, and thus comparable ideal Nernst-Einstein ionic conductivities, which highlights the fact that the differences in the ionic conductivity can be solely attributed to ion correlation effects [30].
A simple and fast computation of ionic conductivities is provided by the advanced electrolyte model [104]. The corresponding approach strongly relies on the influence of different ion motion contributions , and also considers the complexation behavior of the ions. The validity of the model has recently been demonstrated for a series of ion and solvent combinations [105]. Ion correlation effects are also observed in polymer electrolytes, which are often considered as potential components of solid state batteries. In Reference [106], it is observed that, in mixtures of poly(ethylene oxide) (PEO) and LiI salts, correlated motion between the ions occurs, which differs significantly from previous findings for organic solvents. As the lithium ions are strongly coordinated by the oxygen atoms of the PEO chain, distinct spatiotemporal correlations dominate, which can be separated into several contributions in terms of the total dynamic behavior. However, the correlation between lithium and iodide ions contributes most significantly to the dynamics, which can be brought into agreement with simple analytic models by the introduction of phenomenological parameters.
All findings reveal that ion correlation effects are omnipresent, and dominate the dynamics of distinct systems at high salt concentrations. A deeper understanding of ion complex formation principles is thus of fundamental importance in order to develop improved electrolyte solutions.

Solvent-Ion Interactions at the Molecular Scale
In this section, we discuss the properties of solvents and their influence on the formation of ion complexes in more detail. Ion-solvent interactions are dominated by several factors, which are reviewed in the following subsections in terms of several examples.

Dielectric Decrement Effects
As discussed in Sections 2.2 and 2.4, electrostatic attraction drives the formation of ion complexes, whose occurrence reduces the ionic conductivity due to correlation effects. In terms of a simple solution to this problem, the Bjerrum length for the respective solvents should have low values as achieved by high values for the dielectric constant r . The question arises: Which molecular factors affect the dielectric constant of a solvent? Besides more sophisticated Green-Kubo relations [51,65,98], a simple relation for the dielectric constant reads [107][108][109] which is readily applicable for one-component mixtures by the evaluation of the time-averaged squared net total translational dipole moment M 2 trans of the molecules in the respective volume V. Consequently, well-arranged polar solvents with pronounced molecular dipole moments show a high value of the dielectric constant, and are thus beneficial media to reduce ion correlation effects.
Despite this simple view, it has to be mentioned that the value of the dielectric constant also depends on the chosen position, as it can locally change in close vicinity to ions, uncharged solutes, interfaces, and due to inhomogeneities in the solution [64,[110][111][112][113]. Moreover, it is also affected by solvent-solvent, solvent-ion, and ion-ion interactions, which contribute significantly to the net value of the dielectric constant. However, for dilute salt conditions, a recent publication reports that the dielectric constant is mostly affected by solvent-solvent instead of ion-solvent or ion-ion contributions [51]. Although the corresponding findings significantly change at molar salt concentrations [64,66,81], the use of high dielectric constant solvents still remains beneficial.
To rationalize the change of the dielectric constant around charged objects such as ions, one has to remember that the electrostatic field of the ions perturbs the relaxation of the surrounding solvent molecules in terms of multipole interactions in accordance with [63] where r (E φ ) denotes the modified dielectric constant of the solvent due to a solvent-and concentration-specific coupling constant β < 0 in combination with the electrostatic field E φ . It is worth noting that the corresponding variation of r (E φ ), which is called dielectric decrement effect, thus affects the Coulomb energy between ions, and hence promotes the formation of ion complexes due to larger values of the Bjerrum length. Consequently, the Bjerrum length for polar solvents often increases in presence of higher salt concentrations. Furthermore, polar solvent molecules are strongly bound to the ion, which results in compression effects in terms of electrostriction mechanisms [63], which also diminish the mobility of the ions. Consequently, the dielectric decrement effect can be considered as an important example for crucial solvent-ion interactions at local scales. Whereas these findings are of minor importance for low salt concentrations, their importance increases drastically at higher concentrations, as most solvent molecules are now part of the first or second solvation shell around the ions [64,66,114]. Thus, strongly bound solvent molecules govern the global properties of the solution, and thus also decrease the overall ionic conductivity in combination with stronger electrostatic correlation effects.

Molecular Properties of the Solvent: Donor/Acceptor Numbers and Chemical Hardnesses
The results of mean-field theories with a dielectric continuum solvent approach imply that solvents with low Bjerrum lengths induce a lower tendency of ion complex formation. In contrast, recent experimental and atomistic molecular dynamics (MD) simulation results highlighted a more complex behavior [38,87]. As an illustrative example [87], the counterion condensation behavior for sodium ions around highly charged sulfonated oligosulfonic acids in different solvents water, dimethyl sulfoxide (DMSO) and chloroform is shown in Figure 6. With regard to electrostatic mean-field counterion condensation theories [115][116][117], one would expect that the number of bound counterions increases with the order water < DMSO < chloroform in accordance with decreasing Bjerrum lengths [87].

Water
DMSO Chloroform Despite this reasonable assumption, it can indeed be seen in Figure 6 that more sodium ions are bound to the polyelectrolyte in water when compared to DMSO, which implies the following order for the fraction of bound counterions, DMSO < water < chloroform, such that water and DMSO change their positions. Hence, the corresponding results reveal that some solvents induce a strong counterion condensation behavior, which does not agree with the corresponding value of the dielectric constant nor Bjerrum length.
Hence, it is concluded that the amount and the interaction energy of coordinating solvent molecules in the first solvent shell around the ions reveal a massive influence on the solubility of the ions, and thus also on the formation tendency of ion complexes and the underlying free solvation energy [87]. In contrast to these enthalpic interactions, entropic contributions in terms of solvent conformations only account for around 10-20% [38] of the ion free solvation energy, and are thus of minor importance for most solvent-ion interactions [63]. The overwhelming remainder of the aforementioned enthalpic contributions can be assigned to ion-dipole interactions in presence of polar solvent molecules. Comparable findings were reported for adiponitrile [30]; for the different solvents water, N,N-dimethylacetamide, methanol, PC and EC [30,118]; and for lithium salts in presence of urea [119].
To discuss the results of Figure 6 in more detail, it is well-known that DMSO has two lone pair electrons, which imply a high nucleophilicity of the molecule, and thus very likely interact with positively charged groups or cations. The nucleophilic behavior of DMSO is even more pronounced when compared to water, which is a more favorable anion solvent due to its protic properties. Thus, in addition to the value of the dielectric constant, distinct polar solvent molecules have nucleophilic or electrophilic properties, which significantly influence the solvation behavior of anions and cations.
To categorize solvents with regard to their cation solvation properties, Gutmann et al. introduced the so-called empirical donor number (DN) [63,120,121], which accounts for the electron donating properties (donicity) of a solvent molecule. The value of the donor number DN is defined as the negative reaction enthalpy in units of kcal/mol for antimony pentachloride with the respective solvent in a 1:1 ratio in inert solvent 1,2-dichloroethane. It has to be noted that sometimes a donor number is also assigned to ion species, but, to keep the discussion simple, we merely focus on DN values for solvents [122].
As few examples, the values of donor numbers for common solvents are shown in Table 1. The highest values of DNs can be observed for pyridine and acetonitrile, which are thus beneficial solvents in order to solvate cations. In contrast, chloroform has to be considered as a poor solvent in this regard. In fact, the DN values for DMSO, water, methanol, chloroform, and N,N-dimethylacetamide also exactly follow the order observed for the ion binding behavior around polyelectrolytes, as discussed above and in Refs. [51,87]. Finally, all standard carbonate-based solvents and co-solvents of LIBs reveal sufficient DN values around 15-20 (Table 1).
Moreover, Table 1 also shows the values for the Gutmann-Mayer acceptor number (AN). The acceptor number accounts for the electron accepting properties (electrophilicity) of the solvent molecules. This property is closely related to the protic properties of a solvent, i.e., its tendency to donate protons in order to form hydrogen bonds or salt bridges with the anions, which rationalizes the high values for water (AN = 54.8). The corresponding AN values are obtained by using NMR chemical shift values of δ for the 31 P atom of triethylphosphine oxide in the respective solvent relative to the shift in N-hexane [63]. All AN values are commonly normalized to assess AN = 2.348 (δ/ppm), meaning that the δ-values are corrected for the diamagnetic susceptibility χ s of the solvent. Consequently, a suitable solvent to reduce the formation of ion complexes between cations and anions should have high values of AN and DN, in combination with a low viscosity and a high dielectric constant. Although exact predictions within the donor and acceptor number concept are hard to derive, the DN and AN values provide a reasonable framework to identify suitable solvents for specific salts. Closely related with the donor number is the Kamlet-Taft (KT) β KT scale [63], which relies on mean values obtained from protic solvatochromic probes in comparison to structurally similar aprotic probes. The corresponding values β KT are normalized to the interval defined by 0 (cyclohexane) and 1 (hexamethylphosphoric triamide). Values from KT and DN scales can be transformed via the relation DN = 38.2 β KT + 0.5 (16) as outlined in detail in Ref. [63]. However, it has to be noted that DN, AN, and β KT values usually fail in terms of the solvation behavior of ions in multicomponent solutions [51,65]. Consequently, AN and DN values provide a reasonable starting point to estimate the quality of ion solvation in the respective single solvents at low salt concentration. Furthermore, the interaction strength between the ions and the solvent is intimately linked to the corresponding polarizabilities of the components. With regard to this remark, specifically the chemical hardness/softness of the solvent and the ions are discussed as a potential discriminator for beneficial ion-solvent and ion-ion combinations [63]. Herewith, a connection to the HSAB principle [91] and the closely related LMWA approach [70,89] can be drawn straightforwardly. Strict definition of chemical hardness or softness are provided in the context of conceptual density functional theory (DFT) [128][129][130][131]. As a prerequisite, it can be rigorously shown that the electronegativity of a molecule or ion can be defined as [128] as the derivative of the total electronic energy functional E with the number of electrons n under the constraint of a constant external or nuclear potential V. Hence, the hardness of a species, whose value can be interpreted as resistance against electronic changes, follows from [128,131] which is the inverse of the softness [131] in order to discriminate between polarizable and unpolarizable molecules [128]. As a simple estimate, the hardness and thus also the softness of a species can be estimated by the relation [129,131] η H E LUMO − E HOMO (20) in terms of highest occupied and lowest unoccupied orbital energies E HOMO and E LUMO , respectively. Consequently, species with small energy gaps between E HOMO and E LUMO are soft and thus easily polarizable, whereas the presence of large energy gaps is a property of hard species. Moreover, pronounced values of |E HOMO | 0 coincide with a high oxidative stability in accordance with the Koopmans theorem [132,133], which implies that hard solvents such as adiponitrile can be used in high voltage LIBs without significant degradation effects [30].
In addition to computational studies and with regard to experimental measurements, the softness of ions can also be estimated by the difference of the formation energies in the gas phase (loss or gain of electrons), and the transfer energy of the respective ions from gas phase to aqueous solution [63]. Consequently, the closely related softness parameter S exp for solvents relies on the evaluation of the transfer free energy ∆G 0 tr = ∆G 0 solv − ∆G 0 aq for standardized ions from water as the source solvent (aq) to the solvent of interest (solv). In terms of standard ions from a well-defined experimental reference system, one usually relies on soft silver ions, on the one hand, and averaged values for hard sodium and potassium ions, on the other hand. Consequently, the empirical softness parameter for solvents can be calculated by [63] S exp = ∆G 0 tr (Ag + ) − 1/2 ∆G 0 tr (Na + ) + ∆G 0 tr (K + ) /(100 kJ/mol) whereas other definitions and reference ions are also used [63]. In terms of Equation (21), solvents with S exp ≥ 0.25 are considered as soft solvents, and thus preferably solvate soft ions, whereas solvents with S exp ≤ 0 are regarded as hard solvents, and thus preferably solvate hard ions. As shown by the values in Table 1 and in agreement with Equation (21), pyridine is a soft solvent, whereas PC and ethylene glycol are hard solvents. Furthermore, and despite the high DN value, pyridine is a rather poor solvent for lithium ions as one of the hardest ion species. Moreover, the roughly similar DN and AN values in combination with its pronounced chemical hardness rationalize the abundant use of PC as common component in LIB electrolyte solutions.

Solvation of Ions: Benefits of Computational Approaches
In combination with experimental studies, computational approaches are powerful tools which provide detailed insight into molecular interactions as well as electrochemical stabilities [17,55,56,82,132,133]. Closely connected with the conceptual DFT approach, numerical studies on the solvation behavior of lithium salts in organic solvents mostly rely on quantum chemical calculations in terms of DFT computations [55][56][57]. Here, we briefly discuss some standard protocols for reasons of completeness.
As a common approach to study the strength of ion-solvent interactions, the ions are surrounded by a local solvent cluster in accordance with the assumed first solvent shell coordination number [55].
Most often, the ion-solvent cluster is also embedded into a local continuum background with global dielectric constant to increase the accuracy of the calculations. With regard to this attempt, the lithium ion solvation energy can be computed by [56,134] with E tot (Li + · S x ) as total potential energy between the lithium ion and x solvent molecules (S x ), the total lithium ion energy E tot (Li + ), and the total energy of x = 1, . . . , N s solvent molecules without ions E tot (S x ). The number x of coordinating solvent molecules around the lithium ions can be estimated from Raman spectroscopy measurements, or from previous atomistic molecular dynamics (MD) simulations.
In addition to the solvation energy, the total desolvation energy, meaning the energy that is needed to remove the solvent shell from the lithium ions reads [56,134] which can be computed for an arbitrary value of x solvent molecules. A corresponding example is shown in Figure 7, where different clusters of lithium ion-EC complexes ([Li · S x=1−4 ] + ) are depicted. After geometry-optimization, the corresponding values for the total energy can be calculated with regard to adequate DFT functionals and basis sets. In terms of computational efficiency, B3LYP and PBE functionals often provide reasonable values [135], whereas more sophisticated Gaussian correlated G4MP2 approaches are nowadays also often in use [55,82]. However, it has to be noted that dispersion interactions are ignored in simple DFT calculations [136], and have to be taken additionally into consideration, specifically with regard to large solvent-ion clusters [136][137][138].
In addition to potential energies (enthalpies), the corresponding free energies for the ion-solvent clusters can also be evaluated via the relation [135] where ZPE denotes the zero point energy, R denotes the universal gas constant, pV denotes the pressure-volume relation, and Σ denotes the molecular partition function with zero ground state energy.
The corresponding values can be obtained from an a posteriori vibrational frequency analysis [139]. Additional evaluation of electrochemical stabilities in terms of redox potentials by the use of the Nernst equation highlight the applicability of the approach for distinct solvents in good agreement with experimental results [55,132,133]. Although quantum chemical approaches at different levels of theory are useful to estimate the binding energies and the corresponding oxidation and reduction potentials to a sufficient extent, limited information on the dynamic behavior of the complexes is available with regard to such static approaches. Moreover, and due to large computational costs, only the first or second coordination shells around the ions are taken into consideration. However, from our previous discussion in Section 2.2, it can be concluded that the associated length scales are often too short to properly reflect long-range electrostatic interactions. Thus, although aiming to provide a high degree of accuracy, quantum chemical calculations often suffer from finite volume effects or slight numerical variations due to the use of inappropriate functionals. The results of these calculations are rough estimates, which have to be underpinned and verified by experimental measurements or sophisticated MD techniques.

Ions in Multicomponent Solutions: Influence of Co-Solvent and Additive Molecules on the Properties of Ion Complexes
In Section 3, we discuss the properties of ions in single solvents. As shown, the molecular details of the solvent molecules and further specific properties such as the chemical hardness and the DN/AN values are of fundamental importance. However, the underlying considered binary solution, as composed by the dissolved salt and the solvent, reveals a rather low level of complexity. In this section, we focus on the influence of additional components in terms of co-solvents or additive molecules. Useful theoretical approaches to unravel these complex properties are molecular theories of solutions such as the Kirkwood-Buff (KB) theory or inhomogeneous fluctuation theory [34,41,42,45,51,140], which provide a straightforward definition of basic thermodynamic concepts in regards of molecular distribution functions. This section aims to introduce the basic framework and to highlight the benefits of the considered approaches for the study of multicomponent solutions.
It is worth noting that standard electrolyte formulations in LIBs are indeed often composed of distinct lithium salts, solvents, co-solvents, and several additive components. Besides further tasks, these additional compounds are often introduced to increase the solubility of the lithium ion conducting salt, and to decrease the viscosity of the solution [1,3,4,141]. With regard to both aims, often a high dielectric constant solvent such as EC or PC, which usually suffers from a high viscosity [6], is combined with co-solvent or additive components to decrease the viscosity of the solution, and thus enhance charge transport. Unfortunately, typical co-solvents such as EMC or DMC reveal low dielectric constants, which crucially supports the occurrence of ion correlation effects at higher co-solvent concentrations. In more detail, for low co-solvent concentrations, one usually observes a significant increase of the ionic conductivity due to decreasing viscosities in regards of σ eff ∼ 1/η until a maximum value is reached. For higher co-solvent concentrations, the influence of the main high dielectric constant solvent on the solubility of the salt decreases, such that the higher amount of co-solvent molecules induces an increase of ion correlation effects, which negatively affects the ionic conductivity in terms of the dynamic correlation factor (Equation (10)). The question arises: How can the presence of co-solvent molecules induce a change of ion correlation effects? With regard to our previous discussion in Section 3, it can be assumed that local interactions between the co-solvent molecules and the ions are mainly responsible for the observed variations.

Molecular Theories of Solution: Co-Solvent-Ion Interactions
To clarify this question in more detail, one has to study the composition of the solvent shells around the ions as well as the general distribution of species in the solution. Beneficial approaches in this context are molecular theories of solution, and specifically the Kirkwood-Buff (KB) theory [34,41,42,45,52]. As one would assume, a good solubility of ions in mixtures can be achieved by a preferential coordination of highly polar molecules around the ions, whereas the more apolar species should remain in bulk solution [1,3]. This simple view can be associated with a preferential solvation behavior, which can be studied straightforwardly by means of KB theory. Hence, if both ion species are preferably solvated by either solvent or co-solvent molecules, meaning that only one species forms the vast majority of the first solvation shell around both ions, the corresponding behavior is called homoselective solvation. In contrast, if both ion species are solvated by distinct components, which means that the first solvation shells around the ions differ significantly, a heteroselective solvation behavior can be observed. Furthermore, if the minor component of the solution, meaning the co-solvent or additive molecules, is excluded from the first ion solvation shell, one usually describes this observation as a preferential exclusion behavior, whereas a significant accumulation in close vicinity around the ions is called preferential binding [45,52]. After having defined basic principles of ion solvation in multicomponent solution, we now outline some fundamental concepts of KB theory in order to discuss the influence of co-solvent molecules on the ion dissociation behavior in more detail.

Ion Complexes in Presence of Co-Solvent Molecules: Influence on the Ion Dissociation Equilibrium
To study the influence of co-solvent molecules on the ion dissociation-association chemical equilibrium or the closely related salt solubility, we start this section with a simple ion dissociation-association reaction in terms of where the ion complex, either CIP, 1SP, 2SP or AGG, as defined in Section 2.2, is denoted by (X + Y − ) sol , and the fully dissociated ions by X + sol and Y − sol . In principle, a lower amount of ion complexes can be achieved if the chemical equilibrium is significantly shifted to the right side, such that the solvation of the individual ion species is enhanced. For a theoretical description of this effect, we assume the presence of an infinitesimally small ion concentration, and we do not distinguish between CIP, 1SP, 2SP and AGG states. Thus, all ions within the Bjerrum length are considered as associated ion complexes, whereas all other ions are considered as dissociated ions [51]. In addition, the chemical equilibrium constant K is defined by the relations [142] ∆µ = ∑ j ω j µ j = 0 and with the chemical potentials µ j , the standard chemical potentials µ 0 j , and the stoichiometric coefficients ω j [142]. With regard to Equation (25), the chemical equilibrium constant can also be written as [142] K = ∏ j a ω j j (26) with the chemical activity a j = γ j x j , where γ j denotes the chemical activity coefficient in combination with the mole fraction x j . In terms of nomenclature, the salt ions have index "2", whereas the solvent and co-solvent molecules have indices "1" and "3", respectively. Within standard KB theory, it has to be noted that cations and anions have to be considered as indistinguishable [45,46,143].
If a linear influence of co-solvent molecules on the chemical equilibrium and on the chemical potential of the ions is assumed, it follows [42,45,51,144] ∆µ cs = ∆µ 0 − mRTρ 3 (27) with the co-solvent or additive bulk number density ρ 3 , where ∆µ cs can be regarded as a modified difference between the chemical potentials of dissociated and associated ion states due to the presence of co-solvent molecules. A more detailed discussion of the underlying relations and a verification of the approach can be found in Ref. [51]. As a result, a refined Equation (25) reads [45,52] K cs = exp(−∆µ cs /RT) = exp(−∆µ 0 /RT) · exp(mρ 3 ) = K · K app (28) with where K app stands for an apparent chemical equilibrium constant, which yields the corresponding modified chemical equilibrium constant K cs in presence of a co-solvent.
A connection between Equation (28) and the KB theory including the aforementioned binding concepts as introduced in Section 4.1 can be established by the definition of a preferential binding coefficient, which reads [45,51,145] with excess coordination numbers [45] CN xs where g ij (r) denotes the radial distribution function between the ions (index "2") and solvent (index "1") or co-solvent molecules (index "3") at constant pressure p, constant temperature T, and for vanishing ion densities ρ 2 → 0. In more detail, the value of Equation (30) reveals the excess or the deficit of co-solvent or solvent molecules around the ions in comparison to bulk solution. Hence, a preferential exclusion of co-solvent molecules becomes evident for ν 23 < 0, whereas a preferential binding is dictated by ν 23 > 0. The connection in terms of the preferential binding coefficient with the KB theory can be established by the KB integrals as defined by [34,41] which can be regarded as an expression for the excess volume around the ions associated with the surrounding species [42,45,51,146,147]. Moreover, the relation for KB integrals and excess coordination numbers are defined by G ij = G ji and CN xs ij = CN xs ji [45]. It has to be noted that the validity of the KB theory does not rely on the shape or volume of the considered species, and thus is free from any spurious assumptions [34,42]. A comparison with Equation (32) finally yields for Equation (30) the following expression [42] which corresponds to a difference in the excess volumes of co-solvent and solvent molecules in close vicinity around the respective ion state [42,43,[147][148][149][150][151][152][153][154][155][156][157]. With regard to different chemical potentials in terms of the individual ion states (ion complex, as denoted by superscript A vs. dissociated ions, as denoted by superscript F), one can define ∆µ cs = µ F cs − µ A cs according to Equation (27), and thus the difference between the dissociated and the associated ion state is given by [158] where the corresponding changes in the excess coordination number of co-solvent and solvent molecules, and the differences in the KB integrals around the dissociated and the associated ion state are defined by ∆CN 2j = CN xs,F 2j − CN xs,A 2j and ∆G 2j = G F 2j − G A 2j , respectively. In combination with Equations (28) and (34) and with the definition of the chemical potential µ j = µ 0 j + RT ln a j [142], it follows [45,51] with (∂ ln K/∂ ln a 3 ) p,T,ρ 2 →0 = 0 in terms of Equation (28). In more detail, Equation (35) reveals that a distinct accumulation behavior of co-solvent molecules around the dissociated and the associated ions shifts the chemical equilibrium either towards the dissociated (∆ν 23 > 0) or the associated state (∆ν 23 < 0). Consequently, the corresponding expression (Equation (35)) shows that the mere presence of co-solvent molecules, regardless if charged or uncharged, modifies the chemical equilibrium between distinct ion states, and thus also the solubility of species [45,52].

Chemical Equilibrium Constant and Binding Behavior of Co-Solvent Molecules
To derive a detailed expression for the influence of co-solvent molecules, one can assume the existence of a stable chemical equilibrium between the dissociated and the associated ion complex. In combination with Equation (29), one can rewrite Equation (27) accordingly [45,52] ln K cs = ln K + mρ 3 (36) in terms of a linear influence of co-solvent molecules on the chemical equilibrium constant [51]. This assumption is usually valid for low and moderate mole fractions of co-solvent molecules [51]. With regard to Equations (36) and (28), one can also derive for the associated m-value the following expression [45,51] with the derivative of the chemical activity [34,41,42,44,154,159] and the corresponding KB integrals G ij [34,41] as defined in Equation (32). Often, the derivative of the chemical activity (Equation (38)) is written as [42] with the activity coefficient γ 3 of co-solvent molecules [34,48,49], which can be obtained from experiments or refined computational methods [160]. Large values of m in units of inverse number density thus express a significant influence on co-solvent molecules on the respective amount of ion complexes. The expression for the modified chemical equilibrium constant finally reads [51] K cs = K exp(a 33 ∆ν 23 ), (40) which highlights the close relation with a 33 and ∆ν 23 .
In terms of stability requirements, the value for the derivative of the chemical activity is always a 33 > 0, which reveals that the chemical equilibrium is shifted to the associated ion complex state for ∆ν 23 < 0, meaning that co-solvent molecules are more preferentially excluded from the dissociated ions when compared to the associated state. Moreover, for low co-solvent densities, it becomes clear that the activity coefficient approaches unity with γ 3 ≈ 1, and thus accordingly a 33 ≈ 1 due to nearly ideal distributions, which means that the modified chemical equilibrium reads which highlights the fact that solely the binding properties of the co-solvent molecules changes the chemical equilibrium between associated and dissociated ion states.

Beneficial Properties of Co-Solvent Molecules
With regard to Section 4.3, the question arises if a general expression for ∆ν 23 exists, which provides frame-guided design principles to enhance ion dissociation by a specific choice of co-solvents? Although such an expression is not readily derivable, simple considerations indeed highlight general concepts. As already mentioned and with regard to Equation (32), KB integrals can also be interpreted as excess volumes. With regard to this remark, Equation (33) also reads [45,51] where V xs 21 and V xs 23 denote the excess volumes of solvent and co-solvent molecules [45]. For low co-solvent and salt concentrations, one can assume a nearly ideal behavior, such that the derivative of the chemical activity as an estimate for the partial molar volume [34] approaches unity. Consequently, the excess molecular volumes of the species can be replaced accordingly by the partial molar volumes [34] in terms of V . Hence, the difference in the preferential binding coefficients can be written as [34,51] which can be further simplified by the use of [148,158] ∆V pmv 2 ≈ −∆G 21 . If we further assume that the partial molar volume change of ions is small upon dissociation, which can be easily rationalized by combination of lithium ions with large and bulky anions, it follows that ∆G 21 ≈ 0, such that Equation (34) and equivalently Equation (43) reduce to [148,158] ∆ν 23 ≈ ρ 3 ∆G 23 (45) which highlights the fact that the change of the preferential binding coefficient and the change of the chemical equilibrium solely rely on the binding behavior of the co-solvent molecules. Moreover, it is often reasonable to assume ν 23 < 0 for V pmv 1 V pmv 3 , as discussed in more detail in Refs. [51,149]. This assumption can be mostly related to the excluded volume of species around the ions, such that large co-solvent molecules reveal a more negative excluded volume when compared to smaller solvent molecules. The corresponding implications can also be related to the accumulation behavior around the individual ion states, such that it can be assumed that ∆ν 23 < 0 for ∆V xs 3 /∆V xs 1 1. Hence, if the co-solvent molecules have a smaller partial molar volume than the solvent molecules, the chemical equilibrium is shifted to the dissociated ion state, whereas larger co-solvent molecules imply a shift of the chemical equilibrium towards the associated ion complex state. As a result of this discussion, co-solvent molecules should have smaller sizes when compared to solvent molecules in order to enhance the ion dissociation.

Solubility of Salts in Presence of Co-Solvent Molecules
In addition to co-solvent effects on the chemical equilibrium, the general solubility of salts can also be estimated with the help of the KB theory [159]. To derive the respective expression, one first needs to apply Gibbs phase rule [142] F = C − P + 2, where F denotes the degrees of freedom, C denotes the number of components, and P denotes the number of phases. With regard to a salt at the solubility limit, all C = 3 components (salt, solvent, and co-solvent) fulfill the requirement for the chemical equilibrium between the salt in the crystalline and the liquid phase in terms of P = 2. From the Gibbs phase rule, it follows F = 3, which can be attributed to a constant temperature, a constant pressure, and a constant salt density. To estimate the solubility limit of the solution, the chemical potential of the salt can be written as [159] where µ * 2 is the pseudo-chemical potential and Λ T is the thermal de-Broglie wavelength [34,159]. Chemical equilibrium between the salt in the liquid and the crystalline phase in terms of dµ 2 = 0 in accordance with Equation (46) yields [159] with the maximum saturation density S 2 = ρ max 2 . A straightforward expression for the maximum salt solubility in presence of co-solvent molecules reads [159] where the formal equivalence between µ 2 and µ * 2 is used. Moreover, it follows [159] and by combination of Equations (36) and (37) that θ S 2 = ν 23 a 33 and thus [159] which highlights the importance of the molecular distribution functions of all species in the solution. Thus, a crucial connection among co-solvent-salt, solvent-salt, and co-solvent-solvent distributions imposes a crucial influence on the solubility of the salt. For instance, if the co-solvent molecules are preferentially bound to the ions in terms of G 23 − G 21 0, one can observe a significant increase of the maximum salt solubility. Moreover, large values for co-solvent-co-solvent distributions in terms of local clustering behavior and thus G 33 − G 31 1 significantly reduce the solubility of the salt. It is thus desirable to use solvents and co-solvents with reasonable partitioning behavior in order to reduce non-ideal mixing effects. Furthermore, co-solvent molecules should have smaller sizes when compared to solvent molecules, as outlined in more detail in the last subsection.
The corresponding results reveal that KB theory is a useful approach to provide a rationale for future research on electrolyte compositions. The benefits of the theory are its close connection to basic thermodynamic relations, and the independence of the equations from any size or shape of the molecules. Furthermore, KB theory can be applied straightforwardly, and specifically has shown its benefits in terms of the analysis of computer simulations [42,43,45]. It is worth noting that KB theory can also be applied to experimental data to evaluate the underlying radial distribution functions via an indirect technique [34,161,162]. A recent study on the convergence of experimental and simulation data for water in ionic liquids reveals the benefits of combined experimental and computational approaches [163].

Remarks and Future Perspectives
From this article, it becomes clear that there are numerous benefits in terms of combined experimental, computational, and theoretical approaches for the improvement of electrolyte solutions. In fact, less is known about ion solvation in organic solvents, despite many results in aqueous solutions [30,38,51,83]. Whereas the solvation behavior of lithium ions is often studied by various computational approaches (for a recent overview, see Refs. [55,56,82]), it is often assumed that anions reveal only a weak solvation behavior, which is of minor importance for the electrochemical performance of standard LIB or LMB electrolyte solutions [164,165]. The underlying assumption can be related to the abundant usage of high donor number solvents such as EC or PC in combination with the HSAB principle, which both predict a strong coordination of the solvent molecules around the hard cations, whereas the coordination of hard solvents around soft anions is energetically less favorable [38]. Nevertheless, as the discussion points out, a preferred solvation of anions further improves the dissociation behavior of salts, which positively affects the ionic conductivity. It is thus of fundamental importance, besides further studies on cation-solvent interactions, to unravel the anion-solvent coordination behavior in more detail.
Hence, the development of post LIB technologies with different salts, solvents, co-solvents, additives, and electrode materials requires a basic understanding of theoretical concepts and solvation mechanisms. Modern electrolyte formulations are complex solutions, which are crucially affected by the molecular interactions between the individual components. Consequently, a stepwise understanding and definition of main principles paves the way to a more rational design of electrolyte solutions for specific purposes also in future technologies and sodium ion batteries [10,92,166].
With regard to this aim, beneficial properties of solvents, co-solvents and ions are summarized in Figure 8. Even at the lowest level of mean-field theory, it can be shown that high dielectric constant solvents with low viscosities minimize ion correlation effects and thus enhance charge transport. Ion correlation effects are also less pronounced if moderate salt concentrations are used. Moreover, the solvent should have high DN or AN values for beneficial ion solvation, and should be further hardly polarizable. As discussed above, hard solvents are electrochemically stable and provide a good solvation of lithium ions in terms of the HSAB concept. With regard to the data shown in Table 1, these properties are often fulfilled by electrolyte solvents PC, tetramethylene sulfone and γ-Butyrolactone. With regard to the fact that hard solvents usually reveal high electrochemical stabilities [55,132], the benefits of these solvents with regard to broad redox potentials and enhanced lithium ion solvation become obvious.
The largest room for improvement can be assigned to the specific design of additive and co-solvent molecules. Recent research into this direction fostered the introduction of multi-functional electrolyte components, which fulfill several tasks [6,15,167]. To decrease the ion correlation and to promote the occurrence of free ions, small additives or co-solvent molecules have to be introduced, which are even smaller than the solvent species. However, if the solvent and co-solvent molecules significantly differ in their properties, non-ideal mixing effects might occur, which diminish the salt saturation concentration considerably. Moreover, the co-solvents should also significantly decrease the viscosity of the solution. Hence, the corresponding requirements point to the development of an optimization strategy in order to fulfill all duties within acceptable deviations. The corresponding outcomes of our discussion might help to clarify the most important factors to define general design principles.

Summary and Conclusions
The outcomes of several studies discussed in this article highlight that most processes in organic solvent-based electrolyte solutions are very complex. This conclusion can be mostly related to diverse molecular interactions between the individual components, which heavily affect thermodynamic and electrochemical properties of the solution. An important influence on the electrochemical performance of electrolyte formulations can be assigned to the formation of ion complexes. The corresponding consequences with regard to the ion correlation behavior are most important for charge transport mechanisms in terms of ionic conductivities. Specifically, for high salt concentrations as used in standard LIB electrolyte formulations, the consequences of these correlation effects cannot be ignored.
Thus, ion complexes comprise distinct ion pairing states, which differ in their relative occurrence probability in terms of the considered ionic species or solvents. In principle, one can distinguish between contact ion pairs, ion aggregates and different states of solvent-separated ion pairs, for which the size and the influence of the ions on the solvent behavior plays an important role. Closely related, the molecular properties of the solvent, as reflected by DN/AN numbers and the chemical hardness of the molecules, are also of fundamental importance for the occurrence of these effects. Previous simulation and experimental outcomes highlight the significant influence of molecular and chemical properties in terms of the ion solvation and dissociation behavior, which by far exceeds a simple description in terms of mean-field approaches. Thus, the value of the dielectric constant can only be regarded as a rough estimate for the tendency of ion attraction and the formation of ion complexes. Moreover, dielectric constants and the corresponding Bjerrum length are global constants, which rely explicitly on the presence of homogeneous solutions. Consequently, specific ion and solvent effects impose a massive influence on the properties of the electrolyte solution. To achieve a high ion dissociation behavior, solvents with a pronounced anion and cation affinity have to be used. In more detail and in accordance with the HSAB principle, hard ions should be dissolved in hard solvents and vice versa soft ions in soft solvents. Despite this simple classification, standard salts in LIBs include Li + as rather hard ions, which are most often combined with soft anions such as TFSI − . In consequence and in addition to other requirements such as low flammability, low viscosity, and high electrochemical stability, the search of suitable solvents for LIB electrolytes is a challenging task.
Moreover, standard LIB and LMB electrolyte solutions are composed of distinct components. We show that co-solvent molecules impose a significant influence on the chemical equilibrium between ion complexes and dissociated ions. A general framework for a proper and reliable description of ion solvation in multicomponent solution is thus of urgent need. With regard to this remark, the use of molecular theories of solution such as inhomogeneous fluctuation theory or Kirkwood-Buff theory is highly appreciated to establish a more consistent and fundamental point of view [34,140,168,169]. Although the generalization of most effects is far out of reach, the use of these theories allows us to define main principles and research guidelines for an improvement of electrolyte solutions for future LIB, LMB, and DIB technologies.
The content of this article mostly focused on the diversity and complexity of solution effects in multicomponent electrolyte formulations. We ignored several other important factors such as safety issues [167,170,171], as well as electrode-electrolyte interactions, chemical reactions [16,172,173], effects of solvents on solvent/ion co-intercalation [174][175][176], and active electrode material dissolution [177][178][179][180], which also have to be taken into consideration for the development of novel electrolyte compositions. Thus, distinct requirements have to be fulfilled at the same time, which highlights the difficulty of this challenge. Consequently, more insight into molecular mechanisms has to be gained in order to pave a successful way to future electrochemical energy storage devices.

Conflicts of Interest:
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.