Product formation during discharge: a combined modelling and experimental study for Li–O cathodes in LiTFSI/DMSO and LiTFSI/TEGDME electrolytes

Li–air or Li–O 2 batteries are a promising energy storage technology due to the potentially high energy density. However, significant challenges related to reversible charge/discharge of these cells need to be solved. The discharge reaction is generally agreed to proceed via two main routes, which may occur simultaneously. These are the surface mechanism, leading to Li 2 O 2 product formation as surface films, or the solution mechanism, with solid particles formed in the pore structure of the cathode. A detailed understanding of the reaction mechanisms and the dynamic performance of the electrodes is key to further improvements. Here, we present a mathematical model for the discharge process, based on porous electrode theory, including effects of reactant transport and kinetic limitations, as well as the continuous change of properties due to the formation of reaction products via the solution mechanism and the surface mechanism. The model describes the dynamic change in the ratio of the surface and solution mechanism as a function of growth of film thickness, in line with recent findings. The model is able to predict the differences in experimentally obtained discharge curves between dimethyl sulfoxide and tetra ethylene glycol dimethyl ether solvents with 1M LiTFSI, with a minimum of free parameters. The model parameters are based on physical characterization of the materials and the electrodes, or determined by fitting to impedance spectra recorded during the discharge. The developed model and the methodology will provide a powerful tool for optimization of such electrodes.


Introduction
A battery system which has attracted much attention as a possible alternative for energy storage is the lithium-air (Li-air) battery. Due to their high specific energy [1] it is believed that a fully developed Li-air (or in lab-scale Li-O 2 ) battery system can become a serious competitor to fossil fuels in the transportation sector. However, present Li-air batteries still suffer from a poor cycling stability as well as practical specific energies significantly below the theoretical one. One of the main challenges for a successful development of a practical Li-air battery is the search for an electrolyte that is stable in the reactive environment of the cathode with simultaneous support of the electrochemical reactions. Over the last years, good progress has been made in the understanding of the discharge reaction mechanism of these batteries, and the role of the electrolyte; see i.e. Ref. [2]-for a comprehensive review.
Aprotic electrolytes have dominated the Li-air battery development over the past ten years. The reasons for this are the low ionic conductivity of the solid-state electrolyte and the instability of lithium metal against aqueous electrolytes which often causes corrosion and H 2 formation [3,4]. The desired reaction product at the cathode is the lithium peroxide, Li 2 O 2 , and the overall reaction in the cathode is given as It is known that the formation of the discharge product ( Li 2 O 2 ) leads to passivation of the active surface and blockage of the pores [5]. Therefore, the formation and deposition morphology of Li 2 O 2 have been the focus of recent research activities [6,7]. It has been experimentally verified, that the deposition morphology of Li 2 O 2 is strongly dependent on the electrolyte properties [5] and the discharge current [8]. Li 2 O 2 is formed via two main reaction pathways: the surface-based mechanism forming a film on the electrode surface and the solution-based mechanism forming crystalline toroids in the electrolyte, which are attached to the electrode surface [9].
For both the solution and solvation pathway, the formation of Li 2 O 2 follows a multi-step reaction [6]. First, oxygen is reduced by one electron.
The further reaction can either happen on the surface or in the solution. Following the proposed reaction mechanism of Xue et al. [10] it is believed that for the surface mechanism the reduced oxygen reacts with Li + to form LiO 2 as an intermediate product. Subsequently, LiO 2 reacts with one more electron and Li + to form the end product Li 2 O 2 as a film which is covering the electrode surface.
For the solution-based mechanism it is assumed that the reduced oxygen is transported through the electrolyte to a preferred nucleation site to form toroids [10].
It is generally agreed that the dominating mechanisms is determined by the basicity of the solvent, as given by the Gutman donor number (DN) [2], as well as the discharge current density [2]. Typical examples of solvents used for Li-air batteries are TEGDME and DMSO, representing low and high Lewis basicity solvents with DNs of 16.6 and 29.8, respectively [10,12]. The basicity reflects the ability of the solvent to solvate O − 2 /LiO 2 , a governing factor with respect to the dominating mechanism. When high DN electrolytes are used, the solvated LiO 2 intermediates subsequently disproportionate to solid, crystalline, toroidally shaped Li 2 O 2 particles [13][14][15]. With low DN solvents and/or high current densities and overpotentials, the discharge reaction proceeds via the surface growth pathway. For the surface pathway, the discharge capacity is typically limited by the film resistance, which will increase dramatically when the film thickness grows beyond 5-10 nm [2,16,17]. For the solution pathway, the product formation also causes a gradually increasing overpotential and local current density, which again leads to a shift in the reaction mechanism towards the surface pathway, and eventually passivation due to blocking of the active surface. A detailed study of the product formation during charge and discharge (at various states of charge) has been undertaken by Augustin et al. [7], using both SEM surface and cross-sectional analysis.
Mathematical models for the electrochemical performance of porous cathodes for Li-air batteries are usually developed within the framework of the porous electrode theory [18]. A review of modelling work is provided in Ref. [19]. Models have been developed describing either the surface mechanism, see e.g [20], or with a focus on optimization of the pore network [21,22]. In Ref. [22], a continuum model was developed, where the oxygen reduction reaction (ORR) was modelled as a multistep reaction using different rate constants. The rate of solution vs. surface mechanism was controlled simply by the choice of the kinetic parameters. The model was validated based on experimental results. Xue et al. [10] presented an alternative description of the two pathways for the ORR in their model by introducing a so called escape probability as measure of the ratio between the two pathways. The escape probability was a function of the pore diameters and volumes. However, the escape probability was treated as constant and did not take into account the dynamic change of the porosity during discharge.
In this work, we have developed a continuum type model for the Li-air cathode, in order to describe the discharge reaction, incorporating both mechanisms (i.e. the surface and the solution pathway) as simultaneously occurring processes in the electrode. The transport of dissolved oxygen, the kinetics of the ORR, as well as ohmic losses are included in the model. An escape parameter was introduced, taken to be a function of the resulting film thickness, in line with recent findings in the field. This allowed for a dynamic description of the discharge process. The discharge model was coupled to an impedance model to account for the dynamic change of the impedance response during discharge. Model parameters could be determined in a robust manner by combining data from electrochemical impedance spectroscopy and galvanostatic discharge. A good agreement with experimentally obtained discharge curves for LiTFSI/DMSO and LiTFSI/ TEGDME electrolytes could be achieved with a minimum of free parameters.

Experimental
The cathode was manufactured by casting a paste, consisting of Vulcan Ⓡ XC72 carbon black (from Cabot Corp.) and a binder-solution (Polyvinylidene difluoride, PVDF, ≥ 99%, from AME Energy Co., Ltd.), onto carbon paper (PTFE treated TGP-H-60 Toray from Alfa Aesar). For the binder solution 0.6 g of PVDF was mixed with 6 g N-methyl-2-pyrrolidine (NMP, anhydrous, 99.5% from Sigma Aldrich) and stirred for 45 min. NMP was added 3 times each time followed by 25 min shaking in ball milling jars in order to obtain a homogeneous paste. The total carbon:PVDF ratio was 9:1. The paste was applied with a wet casting thickness of 120 μm onto the carbon paper. The wet electrode cast was pre-dried at 60 • C for 45 min and subsequently heated overnight at 85 • C under vacuum in order to remove residual NMP. The dry cast was transferred to a glovebox with argon atmosphere, where discs with diameters of 16 mm were cut from the cast. The gas diffusion electrodes (GDE) were weighed and stored in argon atmosphere until used. The carbon black loading of the electrodes was 3.28 mg ± 0.15 mg. The two electrolyte solvents used were tetraethylene glycol dimethyl ether (TEGDME, 99%, from and Acros Organics), and dimethyl sulfoxide (DMSO, anhydrous, ≥ 99.9%, from Sigma Aldrich), which were mixed with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI, 99.95% trace metals basis, from Sigma Aldrich). Electrolytes were prepared in a glove box (MBRAUN) in an Ar atmosphere with < 0.1 ppm H 2 O and < 0.1 ppm O 2 . A molar concentration of 1M LiTFSI in DMSO and TEGDME was used. All chemicals were used as purchased without further purification. The assembly of the test cells was carried out in a glovebox system with argon atmosphere, using the ECC-Air setup from EL-CELLⓇ GmbH The measurements were carried out at ambient temperature with a VMP-300 multipotentiostat (Bio-Logic Science Instruments). The cells were constantly supplied with an oxygen flow of 4 mL/min (gas purity: 5.0, AGA AS). The cells were discharged and charged with a current density of i = 75 mA/g C , where g C is the mass of Vulcan Ⓡ XC72 carbon black. The potential range was E = 2.4-4.15 V vs. Li/Li + . Before the first discharge the battery system was oxygenated for two hours at open-circuit voltage (OCV).
Impedance measurements were performed during the first discharge at regular intervals, in the frequency range between 20 kHz and 2 mHz. At every point impedance was recorded with 5 measurement points per frequency decade. All impedance measurements were done with the use of a Li metal reference electrode. The capacity of the cells was calculated based on the mass of the carbon (Vulcan Ⓡ XC72). This was done in order to account for the reaction surface area which scales proportionally to the carbon mass, and is common practice for Li-O 2 porous cathodes [2,[23][24][25].

Discharge model
The active layer of the cathode is treated as a network of ideal spherical particles with a diameter of 50 nm to account for the spherical shape and the size of the Vulcan Ⓡ XC72 carbon black particles. In the model the active layer with a measured total thickness of 35 μm is treated as one dimensional and divided into n = 20 segments. The change of the system due to product formation, is updated with a time step size of dt = 10 s. Concentration variations of Li + in the active layer are neglected, as the Li + transport is assumed to not be a limiting factor. This assumption is made on the basis of the high mobility of Li + in the chosen electrolytes [12], and the much higher Li + concentration compared to the oxygen concentration. Also, Li depletion is avoided, as the anode is metallic lithium.
The mass conservation equation for oxygen (i.e. noncharged species in a stagnant electrolyte without convection), can be written as: where c O 2 ,j is the concentration of oxygen, eff j is the effective porosity and i vol,j the volumetric discharge current in each segment j. D O 2 is the diffusion coefficient of oxygen in the electrolyte, = 2 is the number of electrons involved in the overall reaction and F is the Faraday constant. Equation 7 was solved with the following boundary conditions The tortuosity of a bed of spherical particles is described by the well known Bruggeman correlation for spherical particles.
where b j is a factor that describes the blockage of pores as a function of the initial porosity, 0 , and the volume fraction prod,j of the discharge product in each segment Following Faraday's law, the total amount of the formed product ( Li 2 O 2 ) for each segment can be described by integrating the total current I j in each segment over time: Here, the total current I j in each segment is directly related to the volumetric discharge current i vol,j in Eq. 7 by where V j is the volume of each segment. The total current I j in each segment is calculated using the Butler-Volmer equation as described later. The total amount of formed product n j (t) is then used to calculate the volume fraction of the product at a specific time t in each segment. Due to the two different pathways for product formation, Xue et al. [10] proposed a model that is considering both, the Li 2 O 2 growth as a film on the surface of the cathode and the growth as toroids in the electrolyte. The group introduced an escape function based on the pore size distribution in order to determine the proportion of product which contributes to the film formation. The idea of an escape function is adopted in this work, but instead of a simple dependence on the initial pore size distribution, the escape probability is described as a function of film thickness. It is assumed that the probability, escape , of product being formed in the electrolyte decreases with an increase of film thickness and corresponding increase of the resistance. This assumption is supported by a work of Adams et al. [8] who stated that LiO 2 has a stronger adsorption to the Li 2 O 2 film than to the carbon surface. Thus, it is possible that an increase of film thickness leads to a stronger adsorption towards LiO 2 and therefore decreases the probability of the solution-based mechanism to occur. The following escape function was introduced: sol was fitted to the discharge curve and corresponds to the percentage of product which is formed by the solution-based mechanism at the beginning of the discharge. The decrease of the escape probability with a larger film thickness is described by the use of an error function, where is the film thickness averaged over all segments and crit is the critical film thickness. In accordance with Viswanathan et al. [26] the critical film thickness was set to crit = 5 nm. The steepness of the error function was adjusted by division with a factor of 10.
The escape probability is used to calculate the volume fraction of the product resulting from the film formation on the active surface of the cathode, and the following equation is assumed: For the calculation of the film thickness, the film is assumed to grow homogeneously on the surface of the spherical carbon particles. With the given assumptions the film thickness can be derived from the volume fraction of the particle part and the discharge product film that is formed by the surface mechanism: During discharge, an increasing amount of Li 2 O 2 deposits on the surface of the carbon particles, thus causing a decrease of the active surface area and hindering the supply with Li + and O 2 . Consequently, the reduction of the carbon surface by coverage of Li 2 O 2 plays a key role in the discharge behavior of the battery.
Following the empirical equation of Sahapatsombut et al. [27], the reduction of active surface area is calculated by: where A 0,j is the surface area. Exponent p is a geometrical factor that describes the morphology of Li 2 O 2 covering the active surface area. In accordance with the literature, where values between 0.4 and 0.5 are used [27,28], p was fixed to a value of 0.45.
The passivation of the surface area during discharge is known to be caused by the insulating Li 2 O 2 film, which inhibits the charge transfer to the Li 2 O 2 -electrolyte interface, where the electrochemical reactions occur and thus causing an ohmic voltage loss [5,29]. Similar to other approaches in the (12) [30][31][32][33], the film resistance is described in this work by a linear dependence of the film thickness . However, when the film approaches a critical film thickness crit , representing the critical thickness for electron tunnelling [29], electrochemical reactions are kinetically hindered, and we expect a strong increase of the resistance. This strong resistance increase is accounted for by the use of the error function with the same critical film thickness crit = 5 nm as in Eq. 12: where s is the resistivity and is used as a fitting parameter within the range of the literature values [30][31][32]. Note that the unit of R film,dis is [ m 2 ]. It is assumed that the faradaic current can be described by a Butler-Volmer equation. Following the proposed equation of Bard et al. [34] the Butler-Volmer equation can be expressed as: where k is the kinetic rate, the charge transfer coefficient and the activation overpotential. The oxygen reduction rate k is fitted to the first impedance measurement. The Butler-Volmer overpotential is treated as homogeneous across the active layer. For simplification a value of 0.5 was chosen for .
With the given assumptions it is possible to express as a function of the applied discharge current density i a and the averaged oxygen concentration of all segments c O 2 : with where I is the overall total discharge current and A ac the total active surface area. The overall cell voltage can be calculated by the following equation: where R s is the series resistance, i a,cell the discharge current regarding the cross sectional cell area and R film,dis the film resistance. Thus, i a is the discharge current density with respect to the active surface of the carbon, while i a,cell is the discharge current density with respect to the cross section of the battery cell.
The parameters used in the discharge model are given in Table 1.

Impedance model
The elements for the equivalent circuit model were chosen analogous to Eq. 20 of the discharge model. Accordingly, the contributions to the total impedance consist of a series resistance ( R s ), a charge transfer resistance ( R ct ) and the resistance caused by the film ( R film,imp ). To account for the capacitance of the electrode and the film, R ct and R film,imp were connected in parallel to constant phase elements. In electrochemical systems, the impedance response of the charge transfer resistance often follows a distributed reactivity which is generally represented as a constant phase element (CPE) [37]. As a reference electrode was used (lithium), the equivalent circuit models are representing the mechanisms within the cathode. Using the same approach for the description of the film resistance in combination with the series resistance, the overall expression of the equivalent circuit is: The equation for the charge transfer resistance can be directly derived from Eq. 17 of the discharge model: The oxygen reduction rate k was determined from the value of R ct obtained by fitting to the first impedance measurement. This value of the oxygen reduction rate was used in Eq. 18 of the discharge model. Thus, the charge transfer resistance, R ct , is changing with discharge time, depending on the oxygen concentration c O 2 and the faradaic current i far (both computed by the discharge model). The capacitance of the electrode, Q ct was calculated to be 35 F/g in relation to the mass of the carbon black. Also the contribution of the film to the total impedance is changing with discharge time and is calculated by division of the total film resistance (as computed in Eq. 16) with the surface of the film.
where the film surface is calculated by: where r part is the radius of the carbon black particle and s is the film thickness as calculated in Eq. 14 in the discharge model. A value of 50 μF/cm 2 is used for the calculation of the double layer capacitance of the Li 2 O 2 -electrolyte interface, Q film , which is typically in the range of 10-50 μF/cm 2 as given by Huang et al. [38]. As A film is increasing with time the capacitance of the film is increasing accordingly.
To account for heterogeneity's or non-uniform current distributions, ct and film were used as fitting parameters, with reasonable values in the range of 1 to 0.7 [39].
The coupling between the discharge and the impedance model of the first discharge is based on the direct relation between Eqs. 18 and 22 as well as Eqs. 16 and 23. The time-dependent parameters needed for the simulation of the impedance spectra are calculated from the discharge model.

Results and discussion
The discharge performance was simulated based on the parameters given in Table 1, providing the parameters determined from the experimental discharge curves and impedance spectra as shown in Table 2.
The simulated discharge curves are compared to experimental data in Fig. 1 for (a) LiTFSI/DMSO and (b) LiTFSI/ TEGDME. The points at which impedance measurements were made during discharge are marked with crosses and the simulated time of reaching the critical film thickness with crit . The capacity of the cells was calculated based on the mass of the carbon (Vulcan Ⓡ XC72). A discussion on true electrode capacities for porous electrodes is given in [2].
A comparison of the first discharge of LiTFSI/DMSO and LiTFSI/TEGDME in Fig. 1 shows better initial electrochemical properties of the DMSO-based electrolyte than the one based on TEGDME. Three advantages of the system with LiTFSI/DMSO are particularly striking: The higher discharge capacity, the potential at which the discharge takes place and the oxygen reduction rate k, which was found to (24) A film = 4 (r part + s ) 2 be almost 2.5 times higher in LiTFSI/DMSO than in LiTFSI/ TEGDME (see Table 2). The impedance spectra recorded during the first discharge are shown in Fig. 2. Good agreements are obtained by fitting of the simple model to both the discharge curve and the impedance spectra with a minimum of free parameters (see Table 2). Thus, incorporating passivation of electrodes by the growth of the surface film, with a critical thickness corresponding to 5 nm, provides realistic model description.
The first impedance measurements in Fig. 2 of (a) LiTFSI/DMSO shows a smaller semi-circle in the system Fig. 1 Comparison between experimental and simulated data of the voltage-capacity curve for a LiO 2 battery with a LiTFSI/ DMSO and b LiTFSI/TEG-DME as electrolyte. The points for impedance measurements during discharge are marked with circles and the simulated time of reaching the critical film thickness with crit . The cells were discharged with a current density of i = 75 mA/g C within a potential range of E = 2.4-4.15 V vs. Li/Li + Fig. 2 Experimental Impedance measurements performed during the first discharge as marked in Fig. 1 for LiTFSI/TEGDME and LiTFSI/DMSO. The impedance illustration are divided in two sections (I and II) corresponding to the discharge curves. The graphs are a section I of LiTFSI/DMSO b section I of LiTFSI/TEGDME c section II of LiTFSI/DMSO and d section II of LiTFSI/TEGDME compared to the one in (b) LiTFSI/TEGDME. Since the first impedance measurement was made in the very beginning with a negligible amount of product, the size of the semi-circle is dominated by the charge transfer resistance R ct . The difference between the impedance spectra of the DMSO and TEGDME electrolytes reflects the higher oxygen reduction rate k in LiTFSI/DMSO (compare Eq. 22). This is particularly remarkable since the oxygen concentration is lower in LiTFSI/DMSO than in LiTFSI/TEGDME (see Fig. 3b) which has a negative impact on the charge transfer resistance. Laoire et al. [12] were able to show a dependence of the oxygen reduction rate on the Gutmann donor number (DN) of the solvent. Here, the oxygen reduction rate was higher in solvents with a high DN compared to the ones with a lower one. This finding of Laoire et al. is in good agreement with the oxygen reduction rates obtained in this work and is furthermore supported by a more recent work of Amin et al. [40] who also showed the influence of solvents on the reaction kinetics. The results here confirm the previous findings that a high DN does not only result in a higher discharge capacity but also increases the oxygen reduction rate resulting in a lower charge transfer resistance and a higher discharge voltage.
Section II of the discharge curve of LiTFSI/TEGDME shows a stronger voltage decrease and a double plateau formation, which is not in agreement with the curve predicted by the model. Similarly, when comparing the experimental impedance measurements to the simulated ones (see Fig. 2d), the model cannot account for the increase in resistance as observed in section II. The fact that the discharge curve exhibit two plateaus, as demonstrated also by the presence of two peaks in the differential capacity plots of Augustin et al. [7], indicates that a second electrochemical process occurs. In accordance with Ref. [7], the most likely explanation is the Li 2 O 2 formation via the surface growth pathways (Eqs. [3][4][5], known to dominate for the TEGDME solvent. Another process that might lead to this, is the formation of Li 2 CO 3 through side-reactions, as shown by Ottakam et al. [41], who identified side reactions with Li 2 CO 3 as main product both for DMSO and TEGDME electrolyte solvents. DMSO was found to be more stable than TEGDME which results in smaller amounts of Li 2 CO 3 in the DMSO-based electrolyte than in the one based on TEGDME. This finding is supported by the work of Augustin et al. [7] who found Li 2 CO 3 to be present after the first discharge in both electrolytes (LiTFSI/DMSO and LiTFSI/TEGDME). Albertus et al. [32] compared the bandgaps of Li 2 O 2 and Li 2 CO 3 and stated that the electronic resistivity of Li 2 O 2 is orders of magnitude smaller than the one of Li 2 CO 3 . This indicates that already a small layer of Li 2 CO 3 would significantly increase the resistivity of the film, thus contributing to a voltage decay in the discharge curve and increasing the semi-circle diameters of the impedance measurements. Figure 3 shows the simulated distribution of (a) product saturation, i.e. the free pore volume filled with product, and (b) oxygen concentration within the active layer at the end of the first discharge. Both, the product saturation and the oxygen concentration are plotted against the location within the active layer for LiTFSI/DMSO and LiTFSI/TEGDME. The oxygen concentration at 35 μm was fixed to the oxygen saturation concentration of the corresponding electrolyte (see Table 1). The product saturation is on a higher level for LiTFSI/DMSO compared to the one of LiTFSI/TEG-DME. Both product saturation distributions show a gradual decrease with a higher saturation at the oxygen inlet and a lower one at the side facing the separator. Here, the gradient is significantly larger in the system with LiTFSI/TEGDME than in the one with LiTFSI/DMSO. The gradient of the oxygen concentration in behaves in a similar way with a significantly larger gradient for LiTFSI/TEGDME than for LiTFSI/DMSO. For both electrolytes, with a strong oxygen gradient in LiTFSI/TEGDME and almost constant distribution in LiTFSI/DMSO, there is no oxygen shortage during the first discharge.
The simulated film thickness and the ratio between the two product formation mechanisms are shown in Fig. 4a, b,   Fig. 3 Distribution of a product saturation and b oxygen concentration across the active layer for LiTFSI/DMSO and LiTFSI/ TEGDME. The oxygen inlet is on the right side of the graphs at 35 μm and the separator on the left side at 0 μm respectively. The main increase of the electrode resistance as observed from the impedance spectra are attributed to an increase in the film resistance. Figure 4b shows that the percentage of the current consumed by the solution-based mechanism in LiTFSI/TEGDME is 19% at the beginning and linearly decreases towards 10%. On the opposite, the dominant mechanism at the beginning of the discharge of the system with LiTFSI/DMSO as electrolyte is, according to the model, Li 2 O 2 formation via the solution pathway. However, the ratio between the two mechanisms for LiTFSI/ DMSO changes during discharge, as the percentage of 65% for the solution mechanism gradually decreases to 35% from start to end of discharge. The decrease of the escape probability with progressing discharge in the model is caused by the assumed dependency of the film thickness in Eq. 12 and induces a further increase of the film formation rate.
There are several reports in the literature which suggest that electrolyte solvents with a high Gutmann donor number (DN) enhance Li 2 O 2 formation via the solution mechanism rather than a formation via the surface pathway [12,42,43]. Another explanation of the solution-based mechanism was given by Adams et al. [8] who proposed that the solution pathway is caused by solvation of the intermediate product LiO 2 , which is similar to the solvation of O − 2 , dependent on the stabilization within the electrolyte. Disregarding the uncertainty which intermediate product is responsible for the solution mechanism, the product formation takes place via both mechanisms in both electrolytes, as has been reported [2,7]. In the work of Augustin et al. [7], Li 2 O 2 was found as toroids and as film covering the electrode surfaces in both electrolytes, LiTFSI/DMSO and LiTFSI/TEGDME.
In the model, the film thickness is not only influencing the ratio between surface to solution mechanism, but is also responsible for the voltage decrease in the discharge curve of LiTFSI/TEGDME in section I. This finding is based on the good correlation between the discharge curves and the impedance measurements, as a comparison of R ct and R film,imp showed that it is particularly the resistance of the film which is contributing to the increase of the impedance response. The linear increase of the corresponding impedance measurements can therefore be directly related to an increase of R film,dis (Eqs. 16 and 23) in the model. This is in contrast to the discharge model of Lau et al. [44] who assumed the voltage decrease to be caused by the reduction of the active surface, but corresponds to other studies who proposed the film to be responsible for the voltage decline [26,45]. Furthermore, the good agreement of the discharge and impedance model in section I for both electrolytes points out that Li 2 O 2 is most likely the main discharge product during the first part of discharge for both electrolytes. Here it is notable that the resistivity of the film was found to be higher in LiTFSI/TEGDME than in LiTFSI/DMSO. This is in agreement with the studies of Xue et al. [10] who investigated the dependence of the film resistivity on the discharge current and found the resistivity of the film in an electrolyte with DMSO as solvent to be lower than the one in an electrolyte with TEGDME as solvent. In section II of the discharge the model fits to the experimental discharge and impedance behaviour of LiTFSI/DMSO but shows deviations compared to experimental data of LiTFSI/TEGDME. The gradual decrease of the discharge curve with LiTFSI/DMSO as electrolyte begins when approaching the calculated critical film thickness and shows a sharp voltage drop after reaching it (see Fig. 1a). In contrast, the discharge curve of LiTFSI/TEG-DME decreases earlier and independent of the calculated critical film thickness but shows the same sharp voltage drop after reaching it (see Fig. 1b). The calculated film thicknesses at the end of discharge of 6.6 nm for LiTFSI/ DMSO and 6.1 nm for LiTFSI/TEGDME are within the range of 5-10 nm reported as the critical film thickness for electron tunneling [26].
The given results of the discharge behavior and the presence of enough oxygen in both electrolytes, indicate that Fig. 4 Simulated a film thickness and b escape probability escape plotted over discharge time for LiTFSI/DMSO and LiTFSI/TEGDME the main limitation of the first discharge for both electrolytes is the film covering the electrode surface. Since the voltage of both systems drops significantly when reaching the calculated critical film thickness, this suggests electron tunneling through the film which stops when the film is too thick to support electrochemical reactions. This is also in good agreement with the higher capacity (i.e. longer discharge time) of the system with LiTFSI/DMSO as electrolyte, as the growth of the film is slower compared to the one in LiTFSI/TEGDME (see Fig. 4a)). This is assigned here to the enhanced product formation via the solution mechanism in LiTFSI/DMSO.
Besides the analogous behavior of the discharge curve and the impedance measurements, it is also noticeable that the series resistance in Fig. 5 is behaving in accordance to the discharge curve too. Here, the series resistance of the system with LiTFSI/TEGDME is increasing almost from the start, whereas the one of LiTFSI/DMSO remains nearly constant, followed by a gradual increase when approaching the end of discharge. Since the shape of the discharge curve seems to be influenced by the effective thickness of the film, this suggests that the series resistance is also a function of the film thickness. As the series resistance of batteries is dependent on the electrolyte and electrode conductivities as well as the contact resistance [46], three scenarios are possible here. First, the electrolytes are degrading during discharge, which leads to a decay of the electrolyte conductivity. The second and the third scenario are related to the film formation on the surface of the electrodes. It is conceivable that the film growth leads to a contact loss between (a) the carbon particles themselves and (b) the carbon particles and current collector. This loss of contact would result in a decrease of conductivity between the carbon particles and increase of contact resistance between the carbon particles and the current collector. Since significant electrolyte degradation in electrolytes based on TEGDME was reported to occur only during charge [47,48], scenario two and three are more likely to influence the serial resistance in both electrolytes during the first discharge.

Conclusions
A mathematical model has been developed describing the discharge performance of Li-air batteries, including the transport of dissolved oxygen in the porous carbon electrode, the kinetics of the oxygen reduction reaction, as well as the growth of a surface film and formation of dissolved, solid products inside the pores, by the so-called surface and solution mechanisms, respectively. The model accounts for the recent findings in the field by introduction of an escape parameter describing the ratio of the surface and solution mechanism, as a function of the thickness of the discharge product surface film. The model is able to predict the differences in the observed discharge curves between LiTFSI/ DMSO and the LiTFSI/TEGDME electrolytes based on prior assumptions, with a minimum of free parameters. The model parameters are all either based on physical characterization of the materials and the electrodes, or determined by fitting to impedance spectra recorded during the discharge. The ratio between the solution and surface mechanism was found to shift towards the latter with progressing discharge and hence increasing film thickness. The model and the methodology will provide a powerful tool for the optimization of electrodes for future Li-air batteries.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.

Fig. 5
Series resistance plotted over discharge time for LiTFSI/ DMSO and LiTFSI/TEGDME