QC and MD Modelling for Predicting the Electrochemical Stability Window of Electrolytes: New Estimating Algorithm

: The electrolyte is an important component of lithium-ion batteries, especially when it comes to cycling high-voltage cathode materials. In this paper, we propose an algorithm for estimating both the oxidising and reducing potential of electrolytes using molecular dynamics and quantum chemistry techniques. This algorithm can help to determine the composition and structure of the solvate complexes formed when a salt is dissolved in a mixture of solvents. To develop and conﬁrm the efﬁciency of the algorithm, LiBF 4 solutions in binary mixtures of ethylene carbonate (EC)/dimethyl carbonate (DMC) and sulfolane (SL)/dimethyl carbonate (DMC) were studied. The structure and composition of the complexes formed in these systems were determined according to molecular dynamics. Quantum chemical estimation of the thermodynamic and oxidative stability of solvate complexes made it possible to establish which complexes make the most signiﬁcant contribution to the electrochemical stability of the electrolyte system. This method can also be used to determine the additive value of the oxidation and reduction potentials of the electrolyte, along with the contribution of each complex to the overall stability of the electrolyte. Theoretical calculations were conﬁrmed experimentally in the course of studying electrolytes by step-by-step polarisation using inert electrodes. Thus, the main aim of the study is to demonstrate the possibility of using the developed algorithm to select the optimal composition and solvent ratio to achieve predicted redox stability.


Introduction
The almost threefold increase in the energy density of lithium-ion batteries (LIBs) in recent years has been mainly achieved due to increased electrode capacity.The primary factor that limit further increases the energy intensity of cathode materials is the insufficient electrochemical stability of electrolytes.To increase the energy density in the LIB, highvoltage cathode materials with a discharge voltage of ≥4.5 V can be used.These include layered nickel-rich and lithium-and manganese-rich cathode materials [1-4], LiF-MO nanocomposites (M = Mn, Fe, Co) [4,5], LiCoPO 4 [6,7] and others.
The electrolytes used in the production of LIBs typically consist of a 1-1.5 M solution of LiPF 6 in baseline mixtures of linear and cyclic carbonates, to which functional additives are additionally introduced [8][9][10].However, commercial electrolytes, which are stable only to ≤4.3 V, consequently cannot ensure the operability of high-voltage cathode materials [11].In this regard, electrolytes resistant to anodic oxidation at potentials above 5 V are being developed all around the world [12].
Ethylene carbonate (EC), which provides the formation of a stable solid electrolyte layer (SEI) on the surface of the carbon electrode, remains an indispensable component of commercial lithium-conductive electrolytes.EC, like with other components of the electrolyte solution, affects the properties of the cathode-electrolyte interface (CEI), which forms on the surface of cathodes.Under conditions of anodic polarisation of the positive electrode during the charge cycle, molecules of solvent and lithium salt ions undergo electrochemical oxidation at increased potentials.Insoluble products of these poorly studied processes are part of CEI, affecting the stability of cycling and cell performance.Alkyl carbonates still constitute the main components of the electrolyte solution due to their suitable physicochemical properties.Therefore, one of the possible strategies for increasing the electrochemical stability to ≥5 V by modifying the composition of the electrolyte involves the partial replacement of carbonate solvents with more electrochemically stable compounds (whose complete replacement creates problems with the formation of high-quality SEI) [13].To replace carbonate solvents, several electrolytes have been proposed such as fluorinated [14][15][16], nitrile-based [17][18][19], ionic liquids [20] and sulfone-based [21][22][23] electrolytes.Among them, of particular interest are cheap and affordable sulfones representing a group of by-products of oil and chemical industries, whose most widely used member is sulfolane (SL) [23].
LiPF 6 is used as the main lithium salt in commercial electrolytes for LIBs.Despite its many advantages (commercial availability, high conductivity of solutions, the ability to reliably passivate Al foil [10]), the main disadvantage of lithium hexafluorophosphate is its rather low resistance to anodic oxidation [10].In high-voltage electrolytes, more electrochemically stable salts can be used instead of LiPF 6 : LiBF 4 [24], LiBOB [25], LiDFOB [26], LiTFSI [27].In addition, it is known that boron-containing salts improve the electrochemical characteristics of high-voltage cathodes [12].The ancestor of the listed boron-containing salts is lithium tetrafluoroborate LiBF 4 .The anion BF 4 − is characterised by a wider range of electrochemical stability (up to 5.2 V) [28] compared to LiPF 6 − and therefore can serve as an alternative or additional salt in electrolytes for high-voltage LIBs [29].The decomposition product of lithium tetrafluoroborate BF 3 can affect the composition of the protective layers on the electrodes; having the properties of a Lewis acid, it contributes to the dissolution of poorly conductive LiF [30].Electrochemical cells with an electrolyte containing LiBF 4 demonstrate improved electrochemical characteristics [31][32][33][34][35].This makes LiBF 4 promising for use in electrolyte solutions for high-voltage lithium-ion batteries.It should be noted that the electrochemical stability of LiBF solutions depends on the composition of the solvents.It has been shown that a 1 M solution of LiBF 4 in an EC/EMC mixture (3:7) is stable to 6 V [25] and in a mixture of SL/DMC (1:1)-up to 5.5 V [23,[34][35][36][37][38][39][40][41].
Thus, the electrochemical stability of electrolyte solutions is determined not only by the combination of solvents, but also by the nature and concentration of the lithium salt; the need to vary all these parameters in order to obtain the best result makes experimental studies very laborious.In this regard, it seems relevant to develop adequate approaches to the theoretical selection of the optimal composition of an electrolyte solution undergoing redox reactions on the surfaces of both electrodes and forming lithium-conductive passivating layers [42].
The theoretical methods of molecular modelling and quantum chemical calculations widely used to describe electrochemical systems [43][44][45][46][47] are certainly useful for searching for compositions with high electrochemical stability.However, when using computational chemistry methods, the obtained results depend on the construction of an adequate calculation model.Ideally, such a model should take into account all the features of the LIB: the surface of the electrode, the composition and stability of the electrolyte, the introduction and extraction of lithium into/from the surface of the electrodes, the peculiarity of the Batteries 2022, 8, 292 3 of 27 chemical reactions, etc.In reality, the creation of such a model is limited by the possibilities of computational resources and the available level of prediction accuracy.For this reason, the authors of most works use the methods of computational chemistry to solve local problems.For example, classical molecular dynamics methods are used to obtain information about near-order structures in electrolyte solutions, to estimate ionic conductivity, and so on.The obvious advantage of the methods of classical molecular dynamics is the ability to study large atomic-molecular systems, while their main disadvantage is the atomistic level of their description, without taking into account the electronic interaction.On the contrary, quantum chemistry methods can be used to evaluate the electronic parameters of the system, including its thermodynamic and electrochemical stability, as well as the surface profile of the potential energy of chemical reactions.Unfortunately, despite the high accuracy of quantum chemical calculations, this method cannot be used to solve all the necessary problems due to the limited size of the atomic-molecular system under study and the impossibility of assessing the obtained parameters in terms of their dynamics.
For evaluating redox potentials or describing the mechanisms of reactions of electrochemical oxidation of solvents (e.g., sulfolane, propylene carbonate, etc., [48]), quantum chemical calculations are most often used.At the same time, modelling of the studied electrochemical system is typically reduced to a separate consideration of the components that make up this system: isolated solvent molecules, their associates, anionic and cationic solvate complexes and solvated ionic pairs.As a result, a fairly wide range of oxidation (reduction) potential values is obtained, covering all the particles studied.An alternative approach with increased predictive value is to calculate additive values of oxidation (reduction) potentials, taking into account the relative contribution of each particle class [10].Unfortunately, such studies are extremely scarce in the literature.
In quantum chemical calculations, the redox potential assessment is reduced to optimising the geometric parameters of ionic and molecular particles of various grades with the solution of the oscillatory problem.There are two important points here: firstly, the choice of the optimal calculation method and, secondly, the method of accounting for the solvent.In both cases, the adequacy of the choice is tested by correlation with the experiment.Calculations are carried out either using powerful composite methods (G4MP2) [10,[49][50][51][52] or according to approaches based on density functional theory [10,42,[49][50][51][52][53][54][55].As a result, the values of the free Gibbs energies of the initial complex and those of the oxidised (reduced) complex are obtained.Calculations are carried out in the approximation of the gas phase and/or by taking into account the solvent using polarised continuum methods.At the same time, Borodin proposes to use acetone as an implicit solvent, the dielectric constant of which is close to the value of most mixtures of carbonate solvents for commercial electrolytes in LIBs [48].Redox potential values are calculated from the Born-Haber thermodynamic cycle [56].As a rule, adiabatic values are evaluated.If it is impossible to optimise a particular oxidised (reduced) complex, vertical values are evaluated.
Thus, the [30,[57][58][59] solvent molecules (E • ox EC, SL, DMC) and their complexes fall into a range from 5.73 V to 7.42 V, while the experimental data for electrolyte solutions, which appear in a range from 5.3 to 6.7 V [48,57,58], depend on the nature and concentration of the salt, as well as the material of the working electrode.The range of experimental values depends on the type and concentration of salt.Here, the variation in the calculated OP values is obviously related to the calculation and solvent accounting methods.The values of ∆E ox [48] closest to the experimental data are those calculated using composite methods, for example, G4MP2 [59].Calculations in the approximation of the gas phase or in a solvent medium with a low dielectric constant give inflated values of OP [48,60,61].
Quantum chemical calculations of reduction potentials, as well as their experimental measurements, are much less common in the literature.Thus, for EC, ∆E red ; the calculated value is [58,62,63], and experimental values are from 0.11 to 1.9 V [64].
According to experiment [64], the heteromolecular SL-DMC associate is reduced at 0.10-1.38V against 0.67 V-the value estimated by quantum chemistry methods.Anionic solvate complexes are the most "capricious" in terms of optimising the geometric parameters of oxidised and reduced forms.However, for such particles, OP values are often found in the literature.For an unsolvated anion, the BF 4 − value ∆E ox is estimated to be between 8.00 and 8.57 [48], depending on the calculation method; for the anionic complexes of various compositions formed by it, however, the oxidation potential is significantly lower, ranging from 5.79 to 6.39 V.It is important to note that the authors consider anionic complexes only in individual solvents; RP values for anionic complexes with mixed solvate shell have not been found in the scientific literature.
Although features of the structure of cationic complexes have been well described in the scientific literature [27,49,60,61,65,66], these studies typically refer to individual solvents.Reduction potentials are estimated for cationic solvate complexes Li + (EC) 1 and Li + (DMC) 1 [27,49,60,61,66].In the first case, the RP values obtained by different authors fit into the range of 0.45-0.61V, which is quite consistent with the experimental value of 0.54 V.In the second case, the calculated values ∆E red range from 0.22 to 0.60 V.
For the solvated ion pairs of the composition {Li + BF 4 − } (EC) 1 and {Li + BF 4 − } (EC) 2 , calculated values on the acidity potential fall into the range of 6.60-8.74V [48].These neutral particles can undergo both electrochemical oxidation and electrochemical reduction.However, no work has been found in the literature on RP calculations.
As can be seen from the above examples, the redox potential values estimated by the methods of quantum chemistry for various complexes (near-order structures) generally correlate with experimental data.This indicates the adequacy of the quantum chemical methods and methods used to take into account the solvent medium.However, the data available in the scientific literature mainly relate to systems based on a single solvent, whereas, in practice, mixtures of two, three or more solvents are typically used.To predict the electrochemical stability of such solutions, it is necessary to additionally investigate heteromolecular associates and complexes with a mixed solvate shell.In addition, it is noteworthy that the most studied are isolated solvent molecules, their associates and cationic solvate complexes, while solvated ionic pairs are described in less detail, and anionic solvate complexes are practically not considered.At the same time, in order to obtain correctly calculated data on redox potential, it is necessary to consider all local structures existing in the studied electrolyte solution.In addition, in order to increase the predictive value of the calculations, it seems expedient to evaluate the additive values of ORP, taking into account the relative contribution of each of the possible complexes.In this regard, the need to develop a generalised algorithm for theoretical studies of the electrochemical stability of an electrolyte solution becomes obvious, which would allow the analysing of systems of any complexity without restrictions on the component composition, and at the same time reduce the number of particles under study to a reasonable minimum, excluding the unlikely ones from consideration.
In this paper, we investigate a 1 m solution of LiBF 4 in two binary solvent mixtures: EC/DMC (1:1) and SL/DMC (1:1).Earlier [67], we proposed an algorithm for the theoretical evaluation of the electrochemical stability of lithium-conducting electrolytes based on a mixed solvent.In this paper, this algorithm is used to analyse the effect of the nature of solvents on the redox potential of the electrolyte.At the first stage, the dominant types of local structures were determined by the methods of molecular dynamics and their relative contribution was estimated; a molecular model of the electrolyte solution was created.Furthermore, quantum chemical calculations of the structure, thermodynamic parameters and electrochemical stability were performed for all isolated particles; special attention was paid to complexes with a mixed solvate shell.The oxidation and reduction potentials were calculated as additive values.Verification of the adequacy of the constructed models of electrolyte solutions and the correctness of the algorithm of complex theoretical studies was carried out through comparison with experimental data regarding the density and window of the electrochemical stability of electrolyte solutions and chemical shift in NMR spectra.The experimental results obtained for both studied solutions are consistent with the theoretically calculated values, which makes it possible to recommend the approach used in this work for predictive optimisation of the composition of electrolyte solutions.

Materials and Methods
2.1.Theoretical Part 2.1.1.Molecular Dynamic Simulations Molecular dynamic modelling was carried out using the Schrodinger Suite Release 2021-2 software package and the Desmond [68] OPLS3/OPLS4 [69,70].To simulate the selected systems, a model was built where the space is limited by a cubic lattice, taking into account periodic conditions.The number of particles in the design systems corresponded to the composition of the test solution (1 m LiBF4 in EC/DMC (1:1, vol.) and 1 M LiBF4 in SL/DMC (1:1, vol.)).The protocol for preparing the system for simulation included pre-minimisation (conjugate gradient method) and the balancing of components in the NPT ensemble for 2 ns at 298 K and at 400 K, sequentially.The period of the recorded simulation was 50 ns in the ensemble NVT at 298 K. To maintain the constancy of the temperature of the systems, the Nosé-Hoover thermostat was used [71,72].
The study of the structural features of substances, in particular to describe the composition and structure of complexes, was evaluated through the radial distribution function (RDF) and its integral component N (R) [73].Analysis of the results of molecular dynamics made it possible to estimate the frequency of occurrence of complexes of cationic structure and solvated ion pairs.The geometric parameters of the most frequently occurring complexes were used for the subsequent evaluation of their thermodynamic and electrochemical stability by the methods of quantum chemistry.

Quantum Chemical Calculations
All quantum chemical calculations were carried out using GAUSSIAN 16 rev B.01 software [74].Optimisation of the geometric parameters of atomic-molecular systems with the frequency task was carried out at the DFT level using the functionality M052x [75] and the basic set TZVP [76].The empirical dispersion D3 [77] was taken into account.Thermodynamic parameters were calculated in the approximation of the gas phase at 298 K and atmospheric pressure, as well as taking into account the implicit salvation [78].Acetone with a ε = 20.5 was considered as a solvent, which is close to the value of the dielectric constant of most electrolyte mixtures [79].
To assess the thermodynamic stability of the studied complexes of solvated ionic particles, incremental (stepwise) values of the formation energies (∆ f G o inc ) were considered as the Gibbs energy of the reaction of the sequential attachment of the solvent molecule to complexes [67], based on the second consequence of Hess's law.
The adiabatic oxidation potential (AOP) ∆E ox abs and the adiabatic reduction potential (ARP) ∆E red abs for each class of complexes were evaluated according to the Born-Haber cycle [56].
The value of the oxidative potential can be expressed by the formula: The value of the reduction potential can be expressed by the formula: ∆G e -free ionisation energy in the gas phase at T = 298.15K, ∆G • S (M + ) and ∆G • S (M) are the free solvation energy of the oxidised and initial complex, respectively; F is Faraday's constant.The free ionisation energy has the contribution of the ionisation potential, which includes the energy of adiabatic ionisation and ∆G e zero point vibration correction or entropic contribution of −T∆S [48].
The additive value of the oxidative potential characterising the system as a whole was estimated using Maxwell-Boltzmann statistics using the formula E • additive ): where N is the Maxwell-Boltzmann distribution according to the values of the incremental energies of the formation of the i-complex-adiabatic oxidative potential of the E • ads i-complex.
The electrochemical stability window of the electrolyte was determined by step-by-step polarisation according to the procedure described in work [67].For measurements, threeelectrode electrochemical cells with an inert working electrode of platinum, a reversible lithium counter electrode and a lithium comparison electrode were used.During the experiment, a phased polarisation of the cell was carried out with a constant potential difference specified in the range of interest, with stationary current measurement (I ss ).The working electrode was gradually polarised for each given potential value, and the dependence of the current on the time to reach a stationary value was recorded.To determine the value of E ox , the dependences of the stationary current on the potential were constructed.
To determine the oxidation potential, the experiment was conducted in a range from 2.5 V (open circuit voltage) to 5.99 V, relative to Li 0 /Li + .For measurements, potentiostat P-20 × 8 (Elins LLC, Chernogolovka, Russia) was used.
The density of the electrolyte solution was determined using the power of pycnometers with a nominal volume of 5 mL, pre-calibrated with distilled water at a temperature of 25 • C; measurements were performed in three parallels.
The compositions were confirmed using high-resolution NMR.The spectra on the cores 1 H, 7 Li, 11 B, 13 C, 17 O and 19 F were captured on the Bruker Avance III 500 MHz NMR spectrometer.The spectra were recorded at room temperature (24 • C), at frequencies 500, 194, 160, 126, 68 and 471 MHz for 1 H, 7 Li, 11 B, 13 C, 17 O and 19 F, respectively.Liquid samples were placed in standard 5 mm ampoules without the addition of deuteron solvent.To calibrate the chemical shift scale, an additional magnetic field was set from the DMSO-d6 signal as an external standard (2.50 ppm for 1H).

Molecular Dynamics
To determine the density of electrolyte solutions by molecular dynamics, molecular dynamic simulations of the systems under study were carried out within 50 ns.Electrolyte solution density 1 m {Li + BF 4 − } EC/DMC (1:1, vol.) equal to 1.238 ± 0.04 g/cm 3 (Figure S3); for solution 1 m {Li + BF 4 − } SL/DMC (1:1, vol.) equal to 1.213 ± 0.04 g/cm 3 (Figure S4).It has been shown that perturbations in terms of cohesion energy remain within the range of 0.2-0.4kcal/mol (Figures S1 and S2).Insignificant perturbations of the cohesion energy and a slight standard deviation in the values of the solutions' density indicate the stay of the dynamic equilibrium of the system and the adequacy of subsequent measurements of physical values.The density values of the systems under study generally correlate with the values estimated as a result of the experiment.
Analysis of the results of the molecular dynamic simulation of the salt solution in solvents also allows us to estimate the frequency and probability of the formation of various solvated complexes.In the case of the systems studied in this work, when using a mixture of solvents, the formation of complexes with a mixed solvate sphere is probable.In both systems, the lithium cation primarily coordinates the anions BF 4 − around itself via a Li + -F-B 3 − bond at a distance of about 2 Å.This is evidenced by the heights of the first peaks on the curve for Li + -F-B-F 3 − and the function of the radial distribution of particles g(R) in Figures 1a and 2a.
dynamic simulations of the systems under study were carried out within 50 ns.Electrolyte solution density 1 m {Li + BF4 − } EC/DMC (1:1, vol.) equal to 1.238 ± 0.04 g/cm 3 (Figure S3); for solution 1 m {Li + BF4 − } SL/DMC (1:1, vol.) equal to 1.213 ± 0.04 g/cm 3 (Figure S4).It has been shown that perturbations in terms of cohesion energy remain within the range of 0.2-0.4kcal/mol (Figures S1 and S2).Insignificant perturbations of the cohesion energy and a slight standard deviation in the values of the solutions' density indicate the stay of the dynamic equilibrium of the system and the adequacy of subsequent measurements of physical values.The density values of the systems under study generally correlate with the values estimated as a result of the experiment.
Analysis of the results of the molecular dynamic simulation of the salt solution in solvents also allows us to estimate the frequency and probability of the formation of various solvated complexes.In the case of the systems studied in this work, when using a mixture of solvents, the formation of complexes with a mixed solvate sphere is probable.In both systems, the lithium cation primarily coordinates the anions BF4 − around itself via a Li + -F-B 3 − bond at a distance of about 2 Å.This is evidenced by the heights of the first peaks on the curve for Li + -F-B-F3 − and the function of the radial distribution of particles g(R) in Figures 1a and 2a.Analysing the curves in Figure 1a,b corresponding to the {Li + BF4 − } (EC)1(DMC)1 system, it can be noted that the first solvate shell of the lithium cation is likely to include a DMC molecule coordinated to Li + through an oxygen atom of a carbonyl group.EC coordination with respect to the lithium cation seems unlikely, but in this case the cation also tends to be the carbonyl oxygen atom.The total coordination number of the lithium cation tends to be six: up to three anions, two DMC molecules and no more than one EC.Analysing the curves in Figure 1a,b corresponding to the {Li + BF 4 − } (EC) 1 (DMC) 1 system, it can be noted that the first solvate shell of the lithium cation is likely to include a DMC molecule coordinated to Li + through an oxygen atom of a carbonyl group.EC coordination with respect to the lithium cation seems unlikely, but in this case the cation also tends to be the carbonyl oxygen atom.The total coordination number of the lithium cation tends to be six: up to three anions, two DMC molecules and no more than one EC.The formation of anionic solvate complexes, if any, is at a distance of 4-6 Å from the anion (Figure 1c).The oscillation of RDF is weakly expressed; the function N(R) is exponential (Figure 1d).All this leads to a conclusion about the weak interaction of the components.At the same time, the coordination of both solvent molecules with respect to the anion is equally likely.
In the {Li + BF 4 − }(SL) 1 (DMC) 1 system, the coordination of sulfolane with respect to the cation is more pronounced than in the case of the DMC environment.The DMC molecule, just like in a mixture of (EC) 1 (DMC) 1 , is coordinated with the cation via carbonyl oxygen.The total value of the coordination number of the cation tends to be six: the coordination of more than three molecules of the anion, two SL and not more than one DMC is likely (Figure 2a,b).
The behaviour of anionic complexes is generally similar to the {Li + BF 4 − }(EC) 1 (DMC) 1 system: a practical absence of the oscillating nature of g(R) and exponential appearance of N(R) (Figure 2c,d).
Thus, on the basis of data from molecular dynamic simulations, it was shown that in a mixture of EC solvents with DMC, the most likely is the formation of lithium cation complexes, mainly surrounded by DMC molecules (cationic complexes of Li + (DMC) n and solvated pairs {Li + BF 4 − }(EC) n (DMC) n with a predominance of DMC molecules).The most probable complexes in the SL-DMC mixture are cationic complexes of type Li + (SL) n and solvated pairs of type {Li + BF 4 − }(SL) n (DMC) n with a predominance of SL molecules.In addition, anionic complexes are formed in the systems under study.It is logical to consider the composition of possible solvate shells of anions separately from the cation.The near order for the anion BF 4 − is realised at a distance of 4-5 Å from the centre of the mass of the particle.Such an observation seems quite reasonable, since the radius of the anion is much larger than the radius of the lithium cation.
Analysis of the results of the molecular dynamic simulation allows us to estimate the composition and probability of the formation of various solvate complexes in the system under study.Methods of quantum chemistry are used in this work to assess the thermodynamic and electrochemical stability of complexes.

Quantum Chemical Calculations
In this section, we consider the geometric structure of solvents and a number of solvate complexes, and also estimate the incremental energy of the formation of these complexes.
We have divided all the complexes under consideration into: Cationic, in which solvent molecules are coordinated around a lithium cation; Anionic, where BF 4 − is also surrounded by solvent molecules and solvated ionic pairs, and where the LiBF 4 salt is surrounded by EC, DMC or SL molecules.

Structure and Thermodynamic Stability of Complexes
Solvent associates and anionic complexes Solvents are prone to forming associates of type (EC) 2 , (SL) 2 , (EC) 1 (DMC) 1 and (SL) 1 (DMC) 1 .The molecules of many dipolar solvents are prone to self-association and often form dimers or more complex complexes [62].Possible configurations of EC, DMC and SL molecules, as well as their dimers and heteromolecular associates, are given in Section SI (Figure S5).
According to the analysis of molecular dynamic simulations, the formation of negatively charged anionic complexes is possible in the systems under study.Such complexes are interesting from the point of view of the oxidative stability of the entire system as a whole.Such structures are characterised by a low ionisation potential (IP) value and a low electron affinity (Table S1).Theoretically, anionic complexes can make a significant contribution to the electrochemical stability of the system as a whole.
In molecular dynamic simulations, we observe the formation of mixed-type anionic complexes (Figure 3).Such structures are less common in the scientific literature.Geometric parameters of complexes of composition BF 4 − (X) n , where X = EC, DMC, SL, are represented in SM (Figure S6).
(SL)1(DMC)1.The molecules of many dipolar solvents are prone to self-association and often form dimers or more complex complexes [62].Possible configurations of EC, DMC and SL molecules, as well as their dimers and heteromolecular associates, are given in Section SI (Figure S5).
According to the analysis of molecular dynamic simulations, the formation of negatively charged anionic complexes is possible in the systems under study.Such complexes are interesting from the point of view of the oxidative stability of the entire system as a whole.Such structures are characterised by a low ionisation potential (IP) value and a low electron affinity (Table S1).Theoretically, anionic complexes can make a significant contribution to the electrochemical stability of the system as a whole.
In molecular dynamic simulations, we observe the formation of mixed-type anionic complexes (Figure 3).Such structures are less common in the scientific literature.Geometric parameters of complexes of composition BF4 − (X)n, where X = EC, DMC, SL, are represented in SM (Figure S6).In the BF4 − (EC)1(DMC)1 complex, the solvents are coordinated on one side of the anion at a distance of about 2.5 Å, with the DMC molecule slightly closer than the EC In the BF 4 − (EC) 1 (DMC) 1 complex, the solvents are coordinated on one side of the anion at a distance of about 2.5 Å, with the DMC molecule slightly closer than the EC (Figure 3a).A similar description can be given to the complex of type BF 4 − (SL) 1 (DMC) 1 (Figure 3b).
The formation of complexes of type BF 4 − (EC) n , BF 4 − (DMC) n and BF 4 − (SL) n (n = 1.2) (Figure S6) is characterised by positive values of the energy of formation (Table A1 #10-13 and Table A2 #10-13).Interestingly, the addition of a second solvent molecule (Tables A1 and A2 #14) requires less energy than the formation of anionic complexes with pure solvents.
Anionic complexes are unstable structures from the point of view of thermodynamics.Most likely, these are short-lived transition forms of ionic particles.Due to such systems remaining poorly understood, the question of choosing an adequate algorithm for the theoretical study of complexes of this kind becomes an urgent task.
Cation complexes Cationic complexes, in which molecules of only one type of solvent are coordinated around Li + , are well described in the literature [27,49,60,61,65,66].The structures we obtained are similar to those published earlier by other researchers.Typically, the arrangement of molecules is quite symmetrical; thus, with an increase in the size of solvent molecules, the distance from the cation to the solvent atom with which the bond is formed increases proportionally.These patterns are observed for the studied cationic complexes (Figures S7-S9).Most cationic complexes are characterised by negative values of the energy of formation

•
Complexes of the 1 m {Li + BF 4 − } в EC/DMC (1:1) When the DMC and EC molecules are coordinated, the distance from the oxygen atoms to the cation is almost the same (Figure 4).When a mixed complex of Li + (EC) 2 (DMC) 2 is formed, a structure with cross-symmetry of solvents is formed, with the lengths of the bonds differing slightly from each other (Figure 4d).From the point of view of their combinatorics, complexes of mixed composition can be formed along various reaction pathways.Here, the free Gibbs energy of the reaction may vary depending on the sequence of attachment of the electrolyte.In Figure 4, we give the lowest values of the incremental energy of the formation of complexes.For more information, see Tables A1 and A2.

•
Complexes of the 1 m {Li + BF4 − } в SL/DMC (1:1) When SLs are coordinated, the distance from oxygen atoms to the cation exceeds that of DMC and EC molecules when comparing the coordination of a single solvent molecule.In the complex Li + (SL)1(DMC)1, the distance from the cation to the SL is 1.832 Å; with the addition of DMC molecules, this distance increases (Figure 5).From the point of view of their combinatorics, complexes of mixed composition can be formed along various reaction pathways.Here, the free Gibbs energy of the reaction may vary depending on the sequence of attachment of the electrolyte.In Figure 4, we give the lowest values of the incremental energy of the formation of complexes.For more information, see Tables A1 and A2.

•
Complexes of the 1 m {Li + BF 4 − } в SL/DMC (1:1) When SLs are coordinated, the distance from oxygen atoms to the cation exceeds that of DMC and EC molecules when comparing the coordination of a single solvent molecule.In the complex Li + (SL) 1 (DMC) 1 , the distance from the cation to the SL is 1.832 Å; with the addition of DMC molecules, this distance increases (Figure 5).
All mixed complexes of type Li + (SL) n (DMC) n (Figure 5) are characterised by negative values of free energy (∆G f o inc ).Thus, analysing Figures 4 and 5, as well as the data of Tables A1 and A2, it logically follows that the formation of lithium cation complexes with molecules (SL) 1 (DMC) 1 proceeds with a greater energy gain than with molecules (EC) 1 (DMC) 1 .All mixed complexes of type Li + (SL)n(DMC)n (Figure 5) are characterised by negative values of free energy (ΔGf o inc).Thus, analysing Figures 4 and 5, as well as the data of Tables A1 and A2, it logically follows that the formation of lithium cation complexes with molecules (SL)1(DMC)1 proceeds with a greater energy gain than with molecules (EC)1(DMC)1.

Solvated ionic pairs
Negatively charged particles tend to coordinate with a positively charged cation.Given the conformational mobility of anions, the lithium cation coordinates anions around itself in addition to solvent molecules.Here, too, the most interesting complexes are of mixed type (Figure 6).Geometric parameters of the remaining complexes of the type {Li + BF4 − }(X)n, where X = EC, DMC, SL, are represented in SM (Figure S10).

Solvated ionic pairs
Negatively charged particles tend to coordinate with a positively charged cation.Given the conformational mobility of anions, the lithium cation coordinates anions around itself in addition to solvent molecules.Here, too, the most interesting complexes are of mixed type (Figure 6).Geometric parameters of the remaining complexes of the type {Li + BF 4 − }(X) n , where X = EC, DMC, SL, are represented in SM (Figure S10).In the {Li + BF 4 − }(EC) 1 (DMC) 1 complex (Figure 6a), the lithium cation is located closer to the DMC than to the anion or to the EC.The bond length between the oxygen DMC and the cation is 1.874 Å, while the bond length between the boron and Li + is 1.970 Å and between oxygen EC and cation is 1.909 Å.In the case of {Li + BF 4 − } (SL) 1 (DMC) 1 (Figure 6b), the cation is also closer to DMC than to SL, but this difference is very small.In other words, the coordination of solvent molecules in relation to the ionic pair is unchanged.The formation energy of −0.18 eV is the same for both these complexes.In the {Li + BF4 − }(EC)1(DMC)1 complex (Figure 6a), the lithium cation is located closer to the DMC than to the anion or to the EC.The bond length between the oxygen DMC and the cation is 1.874 Å, while the bond length between the boron and Li + is 1.970Å and between oxygen EC and cation is 1.909 Å.In the case of {Li + BF4 − } (SL)1(DMC)1 (Figure 6b), the cation is also closer to DMC than to SL, but this difference is very small.In other words, the coordination of solvent molecules in relation to the ionic pair is unchanged.The formation energy of −0.18 eV is the same for both these complexes.
Unlike with the associations of solvent molecules and anionic solvate complexes, the formation of solvated and unsolvated ionic pairs always proceeds with an energy gain.The values of the incremental energy of the formation of the complex in all cases considered are negative (Tables A1 and A2 # [30][31][32][33][34].
According to the analysis of NMR signals of the systems under consideration, the ratio of proton signal intensity for EC:DMC solvents is 1.47:1, and SL:DMC = 1:1.1 in the mole ratio.In the considered concentration of salt in solution, one ionic pair accounts for up to 6-7 solvent molecules.Given the structure of solvated complexes, it can be assumed that such a number of molecules makes up the sphere of near and far order.In this paper, we considered near-order complexes and estimated chemical shifts using quantum chemistry methods for systems {Li + BF4 − }(EC)1(DMC)1 and {Li + BF4 − }(SL)1(DMC)1.There is a correlation between the values of chemical shifts obtained by NMR and those estimated by QM methods (Figure 7).Since the approximation of QM-calculated chemical shifts relative to the experimentally obtained values on the nuclei 1 H, 7 Li, 11 B, 13 C, 17 O and 19 F is a straight line at an angle of 45°, there is reason to believe in the correctness of the chosen structure of associates for quantum chemical calculation.Thus, the totality of methods of molecular Unlike with the associations of solvent molecules and anionic solvate complexes, the formation of solvated and unsolvated ionic pairs always proceeds with an energy gain.The values of the incremental energy of the formation of the complex in all cases considered are negative (Tables A1 and A2 # [30][31][32][33][34]. According to the analysis of NMR signals of the systems under consideration, the ratio of proton signal intensity for EC:DMC solvents is 1.47:1, and SL:DMC = 1:1.1 in the mole ratio.In the considered concentration of salt in solution, one ionic pair accounts for up to 6-7 solvent molecules.Given the structure of solvated complexes, it can be assumed that such a number of molecules makes up the sphere of near and far order.In this paper, we considered near-order complexes and estimated chemical shifts using quantum chemistry methods for systems {Li + BF 4 − }(EC) 1 (DMC) 1 and {Li + BF 4 − }(SL) 1 (DMC) 1 .There is a correlation between the values of chemical shifts obtained by NMR and those estimated by QM methods (Figure 7).In the {Li + BF4 − }(EC)1(DMC)1 complex (Figure 6a), the lithium cation is located closer to the DMC than to the anion or to the EC.The bond length between the oxygen DMC and the cation is 1.874 Å, while the bond length between the boron and Li + is 1.970Å and between oxygen EC and cation is 1.909 Å.In the case of {Li + BF4 − } (SL)1(DMC)1 (Figure 6b), the cation is also closer to DMC than to SL, but this difference is very small.In other words, the coordination of solvent molecules in relation to the ionic pair is unchanged.The formation energy of −0.18 eV is the same for both these complexes.
Unlike with the associations of solvent molecules and anionic solvate complexes, the formation of solvated and unsolvated ionic pairs always proceeds with an energy gain.The values of the incremental energy of the formation of the complex in all cases considered are negative (Tables A1 and A2 # [30][31][32][33][34]. According to the analysis of NMR signals of the systems under consideration, the ratio of proton signal intensity for EC:DMC solvents is 1.47:1, and SL:DMC = 1:1.1 in the mole ratio.In the considered concentration of salt in solution, one ionic pair accounts for up to 6-7 solvent molecules.Given the structure of solvated complexes, it can be assumed that such a number of molecules makes up the sphere of near and far order.In this paper, we considered near-order complexes and estimated chemical shifts using quantum chemistry methods for systems {Li + BF4 − }(EC)1(DMC)1 and {Li + BF4 − }(SL)1(DMC)1.There is a correlation between the values of chemical shifts obtained by NMR and those estimated by QM methods (Figure 7).Since the approximation of QM-calculated chemical shifts relative to the experimentally obtained values on the nuclei 1 H, 7 Li, 11 B, 13 C, 17 O and 19 F is a straight line at an angle of 45°, there is reason to believe in the correctness of the chosen structure of associates for quantum chemical calculation.Thus, the totality of methods of molecular Since the approximation of QM-calculated chemical shifts relative to the experimentally obtained values on the nuclei 1 H, 7 Li, 11 B, 13 C, 17 O and 19 F is a straight line at an angle of 45 • , there is reason to believe in the correctness of the chosen structure of associates for quantum chemical calculation.Thus, the totality of methods of molecular mechanics, quantum chemistry and experimental NMR data unequivocally confirms the existence of mixed composition complexes.

Electrochemical Stability
The electrochemical stability of the electrolyte solution is determined by the set of formed complexes.The resistance of the electrolyte to oxidation, which can be considered as an additive value, consists of the thermodynamic and electrochemical stability of individual components: isolated solvent molecules, their associates, anionic solvent complexes and unsolvated anions, as well as solvated and unsolvated ionic pairs.Solvents and their associates For solvents and their associates, oxidative stability is estimated in the range of 5.73-7.42V. Thus, evaluated using the methods of quantum chemistry [30,48,62,[79][80][81], values of OP (∆E ox abs ) for solvent molecules (EC, SL, DMC) and their associates fall in a range from 5.58 V to 8.50 V [30,48,62,[79][80][81].The experimental data fall in a range from 4.6 to 6.7 V [57,58,64,[82][83][84].RPs (∆E red abs ) for almost all solvents are negative values except SL (Table A1 #3), (SL) 1 (DMC) 1 -1 (Table A1 #7) and (SL) 1 (DMC) 1 -1 (Table A1 #8).From the literature, it is known that the RP of EC < 1 [30], for SL is 0.4 V [85].There are two different RPs of (SL) 1 (DMC) 1 -1: 0.102 V and 1.386 V [64].According to the incremental Gibbs energy values of these complexes, ∆ f G o inc , (Table A1 #1-8 and Table A2 #1-8), the formation of such associates is unlikely.∆ f G o inc values are positive.Anionic complexes

•
Oxidation potentials Negatively charged complexes are characterised by a low IP value (Table S1), as well as a lower value of the OP (

• Reduction potentials
The RPs for anionic complexes cannot be calculated.It commonly occurs that the optimising procedure of geometric parameters of the N-1 system (electron removal) ends unsuccessfully with a decomposition of the negatively charged anion-solvent complex.Many [27,48] authors bypass this problematic issue by calculating the vertical values of the OP.However, the calculation of vertical values does not imply a complete optimisation of the geometric parameters of the N-1 system, but represents a single (Single Point) calculation, which does not fully take into account the processes occurring in the system.
Cation complexes

• Oxidation
Oxidising potentials for cationic complexes have not been evaluated.Obviously, the removal of the second electron from a positively charged system would require a large quantity of energy and is therefore unlikely.

• Reduction potentials
The RP values range from −1.02 to 0.52 V, depending on the structure of the complex.The highest RP value is characteristic of the cationic complex of mixed type Li + (EC) 1 (DMC) 1 .The RP values of cationic complexes for the system (EC) 1 (DMC) 1 can be arranged in the following series: Li + (DMC The values of Li + (EC) 1 and Li + (DMC) 1 are almost the same.The energy of the Li + (EC) 1 (DMC) 1 complex formation does not depend on the type of reaction (Table A1 #23).The sequence addition of one and more molecules (DMC or EC) to the complex reduces the RP values (Table A2 #24).
Solvated ionic pairs
Figure 8 shows the dependence of the stationary current on the applied potential E in a range from 2.50 to 5.99 V. Based on the shape of the curve, the stability of the solution can be concluded to be based on ethylene carbonate.For a sample of composition 1 m {Li + BF 4 − } in a mixture of EC/DMC (1:1), a sharp increase in current is observed at 5.5 V; for a sample of 1 m {Li + BF 4 − } in a mixture of SL/DMC (1:1), a corresponding increase occurs at 5.75 V.The maximum current for the 1 m {Li + BF 4 − } sample in the EC/DMC mixture (1:1) was 27 µA, while for the 1 m {Li + BF 4 − } sample in the SL/DMC mixture (1:1) it was 17 µA.
in a range from 2.50 to 5.99 V. Based on the shape of the curve, the stability of the solution can be concluded to be based on ethylene carbonate.For a sample of composition 1 m {Li + BF4 − } in a mixture of EC/DMC (1:1), a sharp increase in current is observed at 5.5 V; for a sample of 1 m {Li + BF4 − } in a mixture of SL/DMC (1:1), a corresponding increase occurs at 5.75 V.The maximum current for the 1 m {Li + BF4 − } sample in the EC/DMC mixture (1:1) was 27 µA, while for the 1 m {Li + BF4 − } sample in the SL/DMC mixture (1:1) it was 17 µA.The density of the studied electrolyte solutions was 1.2733 g/cm 3 for 1 m {Li + BF4 − } EC/DMC (1:1, vol.) and 1.2188 g/cm 3 for 1 m {Li + BF4 − } SL/DMC (1:1, vol.).

Discussion
Based on the thermodynamic parameters of each of the possible complexes and considering their electrochemical stability, we estimated the additive values of the redox potentials of 1 m LiBF4 in mixtures of EC/DMC (1:1) and SL/DMC (1:1) solvents.
Each geometric structure that is part of the electrolyte system is an important component of the electrolyte.Due to the fact that the electrolyte is a multicomponent system, we can state that oxidation and reduction in it occurs gradually.So during theoretical study, the contribution of each component to this process was evaluated.
In the case of {Li + BF4 − }(EC)1(DMC)1, the additive redox potentials are 6.19 and −0.02 V, respectively (Figure 9).These values are a window of electrochemical stability.When the AOP value of 6.19 V is reached, all complexes are subject to electrochemical oxidation reactions.As a result, significant rearrangements of their components occur up to the breaking of bonds and the formation of a free HF molecule.Here, the first occurrence at 5.43 V will be oxidised: complex BF4 − (EC)1(DMC)1 (Figure 9A).The ionic pair of mixed composition {Li + BF4 − }(EC)1(DMC)1 is characterised by the highest resistance to oxidation.Analysis of the NMR spectrum with QC calculations of NMR shifts confirms the existence of such a complex in the system under consideration.The density of the studied electrolyte solutions was 1.2733 g/cm 3 for 1 m {Li + BF 4 − } EC/DMC (1:1, vol.) and 1.2188 g/cm 3 for 1 m {Li + BF 4 − } SL/DMC (1:1, vol.).

Discussion
Based on the thermodynamic parameters of each of the possible complexes and considering their electrochemical stability, we estimated the additive values of the redox potentials of 1 m LiBF 4 in mixtures of EC/DMC (1:1) and SL/DMC (1:1) solvents.
Each geometric structure that is part of the electrolyte system is an important component of the electrolyte.Due to the fact that the electrolyte is a multicomponent system, we can state that oxidation and reduction in it occurs gradually.So during theoretical study, the contribution of each component to this process was evaluated.
In the case of {Li + BF 4 − }(EC) 1 (DMC) 1 , the additive redox potentials are 6.19 and −0.02 V, respectively (Figure 9).These values are a window of electrochemical stability.When the AOP value of 6.19 V is reached, all complexes are subject to electrochemical oxidation reactions.As a result, significant rearrangements of their components occur up to the breaking of bonds and the formation of a free HF molecule.Here, the first occurrence at 5.43 V will be oxidised: complex BF 4 − (EC) 1 (DMC) 1 (Figure 9A).The ionic pair of mixed composition {Li + BF 4 − }(EC) 1 (DMC) 1 is characterised by the highest resistance to oxidation.Analysis of the NMR spectrum with QC calculations of NMR shifts confirms the existence of such a complex in the system under consideration.
The ARP is measured at −0.02 V, where the cationic complex of the mixed composition Li + (DMC) 1 (EC) 1 is characterised by the highest reduction potential value of 0.98 V. Here, the tendency to acquire an electron increase depending on the volume of the complex should be noted.At the same time, cationic complexes containing only EC molecules are more prone to reduction than those containing only DMC molecules.Cation complexes of mixed composition are characterised by positive ARP values.Among the solvated ionic pairs, the highest RP values are characteristic of complexes containing pure EC.The lowest RP is characteristic of the mixed system {Li + BF 4 − }(EC) 1 (DMC) 1 .In general, the analysis of thermodynamic and electrochemical values for the system {Li + BF 4 − }(EC) 1 (DMC) 1 indicates that the presence of the DMC molecule in all the complexes under consideration leads to an increase in the thermodynamic stability of the system and an increase in redox potentials.The ARP is measured at −0.02 V, where the cationic complex of the mixed composition Li + (DMC)1(EC)1 is characterised by the highest reduction potential value of 0.98 V. Here, the tendency to acquire an electron increase depending on the volume of the complex should be noted.At the same time, cationic complexes containing only EC molecules are more prone to reduction than those containing only DMC molecules.Cation complexes of mixed composition are characterised by positive ARP values.Among the solvated ionic pairs, the highest RP values are characteristic of complexes containing pure EC.The lowest RP is characteristic of the mixed system {Li + BF4 − }(EC)1(DMC)1.In general, the analysis of thermodynamic and electrochemical values for the system {Li + BF4 − }(EC)1(DMC)1 indicates that the presence of the DMC molecule in all the complexes under consideration leads to an increase in the thermodynamic stability of the system and an increase in redox potentials.
The additive redox potentials of the {Li + BF4 − }(SL)1(DMC)1 system are 6.41 and −0.39 V, respectively (Figure 10).Here, the anionic complex BF4 − (DMC)1 will be oxidised first at 5.60 V. Replacing the EC with SL leads to an increase in AOP for a mixed composition complex.A similar phenomenon is observed for solvated ion pairs; as confirmed experimentally, the addition of SL increases the AOP complexes to make the system more stable to oxidation.The additive redox potentials of the {Li + BF 4 − }(SL) 1 (DMC) 1 system are 6.41 and −0.39 V, respectively (Figure 10).Here, the anionic complex BF 4 − (DMC) 1 will be oxidised first at 5.60 V. Replacing the EC with SL leads to an increase in AOP for a mixed composition complex.A similar phenomenon is observed for solvated ion pairs; as confirmed experimentally, the addition of SL increases the AOP complexes to make the system more stable to oxidation.
For all complexes containing SL surrounded by Li + , the RP is characterised by negative values.The lowest value of the RP is characteristic of the cationic complex Li + (SL) 4 .The tendency to accept an electron in cationic complexes of mixed composition increases with an increase in the number of solvent molecules in the solvate sphere of the cation.For solvated ion pairs, all systems containing SL in the coordinate sphere are characterised by negative values of RP, as also confirmed by literature data.SL does not participate in the formation of SEI because of its low potential of reduction [34].It is one of the problems of using electrolytes based on SL in cells with carbon anodes.Electrolytes containing ethylene carbonate are often reduced at potentials below 0.8 V on the surface of various anode materials [9,27,86] including graphite [87], carbon nanostructures [88], silicon [89] and metal oxides [90], leading to the formation of SEI layers.
Thus, replacing the EC with the SL molecule increases the electrochemical stability window of the system.This is logical, since the stronger electron acceptor group in sulfones compared with carbonates increases resistance to oxidation.When replacing ethylene carbonate with sulfolane in the electrolyte, it was shown experimentally that the beginning of the oxidation processes shifts by 0.25 V from 5.5 V, in the case of a composition of 1 m {Li + BF 4 − } EC/DMC (1:1), to 5.75 V, for a composition of 1 m {Li + BF 4 − } SL/DMC (1:1).For all complexes containing SL surrounded by Li + , the RP is characterised by negative values.The lowest value of the RP is characteristic of the cationic complex Li + (SL)4.The tendency to accept an electron in cationic complexes of mixed composition increases with an increase in the number of solvent molecules in the solvate sphere of the cation.For solvated ion pairs, all systems containing SL in the coordinate sphere are characterised by negative values of RP, as also confirmed by literature data.SL does not participate in the formation of SEI because of its low potential of reduction [34].It is one of the problems of using electrolytes based on SL in cells with carbon anodes.Electrolytes containing ethylene carbonate are often reduced at potentials below 0.8 V on the surface of various anode materials [9,27,86] including graphite [87], carbon nanostructures [88], silicon [89] and metal oxides [90], leading to the formation of SEI layers.
Thus, replacing the EC with the SL molecule increases the electrochemical stability window of the system.This is logical, since the stronger electron acceptor group in sulfones compared with carbonates increases resistance to oxidation.When replacing ethylene carbonate with sulfolane in the electrolyte, it was shown experimentally that the beginning of the oxidation processes shifts by 0.25 V from 5.5 V, in the case of a composition of 1 m {Li + BF4 − } EC/DMC (1:1), to 5.75 V, for a composition of 1 m {Li + BF4 − } SL/DMC (1:1).

Conclusions
The main purpose of the work was to develop and test an algorithm for assessing the stability of electrolytes, using a set of methods derived from quantum chemistry and molecular dynamics.The combination of these methods is the main novelty of our work.Using the approach described in this paper, it is possible to predict the value of the oxidation and reduction potential depending on the solvent ratio.In this work, a complex of theoretical and experimental methods was used to evaluate the window of electrochemical stability of electrolytes of compositions 1 m {Li + BF 4 − } in a mixture of EC/DMC (1:1) and 1 m {Li + BF 4 − } in a mixture of SL/DMC (1:1).The composition and structures of all the complexes that can form in these electrolytes have been described.A number of mixed-composition complexes are characterised for the first time.Additive redox potentials were evaluated as a result of a large-scale theoretical study.This new algorithm can be used for different electrochemical systems including salts, solvents and additives.
It was experimentally shown that the processes of reducing the ethylene carbonate-based electrolyte begin at higher potentials than the sulfolane-based sample.Oxidative processes in a sulfolane-based solution begin at higher potential values than in an ethylene carbonatebased solution.A similar conclusion can be drawn by analysing the results of theoretical calculations.In other words, the obvious correlation observed between experimental and theoretical data confirms the possibility of using computational chemistry methods to select the optimal composition and ratio of a solvent in order to predict redox stability.
Table A1 #10-14 and Table A2 #10-14).The ∆ f G o inc values indicate that the formation of such systems in a salt solution requires energy.However, anionic complexes can contribute to the additive value of the OP salt {Li + BF 4 − } in studied electrolyte solutions.The formation of solvate complexes of the anion BF 4 − with molecules EC and/or DMC, and/or SL (Table A1 #10-14, Table A2 #10-14), is also thermodynamically unlikely, despite which of the solvents is coordinated close to the anion first.The fairly high values (∆E ox abs ) of all anionic solvate complexes correspond to a range from 5.43 to 6.16 V.A high value of ∆E ox abs is also characteristic of the unsolvated anion BF 4 − (8.54 V).All calculated values correlate with the literature data.

Table A1 .
Thermodynamic and oxidative parameters of solvent molecules and various complexes existing in the system 1 m LiBF 4 EC/DMC: ∆ f G o inc -incremental (stepped) Gibbs energy of formation (eV); N-Maxwell-Boltzmann energy distribution; ∆E ox abs -adiabatic oxidation potential estimated in this work (M052X/TZVP); ∆E ox abs-calc -adiabatic oxidation potential, published earlier, estimated on the basis of quantum chemical calculations; ∆E ox abs-exp -experimentally measured value of the oxidation potential; ∆E red abs -adiabatic reduction potential estimated in this work (M052X/TZVP); ∆E red* abs-calc -adiabatic reduction potential published earlier, which was estimated on the basis of quantum chemical calculations; ∆E red* abs-exp -experimentally measured value of the reduction potential.

Table A1 .
Cont.Values obtained on quantum chemical calculation and published in articles; 2 values obtained from experimental data and published in articles; 3 relevant data were not found in the literature. 1

Table A2 .
Thermodynamic and oxidative parameters of solvent molecules and various complexes existing in the system 1 m LiBF 4 SL/DMC: ∆ f G o inc -incremental (stepped) Gibbs energy of formation (eV); N-Maxwell-Boltzmann energy distribution; ∆E ox abs -adiabatic oxidation potential estimated in this work (M052X/TZVP); ∆E ox abs-calc -adiabatic oxidation potential, published earlier, estimated on the basis of quantum chemical calculations; ∆E ox abs-exp -experimentally measured value of the oxidation potential; ∆E red abs -adiabatic reduction potential estimated in this work (M052X/TZVP); ∆E red abs-calc -adiabatic reduction potential published earlier, which was estimated on the basis of quantum chemical calculations; ∆E red abs-exp -experimentally measured value of the reduction potential.