An Overview on Transport Phenomena within Solid Electrolyte Interphase and Their Impact on the Performance and Durability of Lithium-Ion Batteries

: The nature of the electrode–electrolyte interface has an impact on the performance and durability of lithium-ion batteries (LIBs). The initial electrolyte’s thermodynamic instability at the anode–electrolyte interface in LIBs results in the formation of a passivation layer, called solid electrolyte interphase (SEI). The initial dense and intact layer allows Li + transport and restricts electron tunneling, thus preventing electrolyte decomposition and ensuring the electrochemical stability of a battery. However, the growth of this layer can reduce the availability of active lithium and electrolyte, and ultimately lead to an irreversible battery capacity fade. Investigating the transport phenomena of lithium ions within SEI is crucial for understanding its formation and growth. Nonetheless, accurately describing all relevant mechanisms is challenging due to its complex and multiscale nature. An overview of current computational efforts to study Li + transport within SEI is given in this article, ranging from electronic/atomistic scale simulations to macroscopic models. The drawbacks and advantages of the proposed numerical approaches are summarized along with the obstacles that need to be overcome to obtain accurate experimental data, identiﬁed on the basis of the most recent literature evidence. We highlight collaboration gaps between modeling and experimental approaches, as well as the urgent need for new multiscale models, to gain a better understanding of such a crucial transport phenomenon.


Introduction
Rechargeable or secondary batteries are devices that can convert chemical energy into electrical energy and vice versa for energy storage purposes. Secondary batteries have been used since the 1800s, but rechargeable lithium-ion batteries (LIBs) were not introduced until a century later [1]. In the 1900s, there was a technological revolution that resulted in the development of portable electronic devices such as video cameras, cell phones, and computers that required batteries with greater capacity or equal capacity with reduced size and weight. As a result, the search for a new rechargeable, small, and lightweight battery began, thus leading to the development of LIBs starting in the 1970s, with significant progress in the 1980s and 1990s [2]. At present, rechargeable batteries also play an important role in enabling the transition to a zero-emission society [3][4][5]. For example, electrochemical storage holds the promise to effectively cope with the intermittent nature of renewable energy technologies and their integration into energy grids [6]. LIBs continue to dominate the rechargeable battery market, but they are approaching their performance limits. Therefore, it is necessary to create the low-cost, high-performance, and long-lasting batteries of the future [7,8]. 2 of 30 In their commercial implementations, lithium-ion batteries are typically composed of several electrochemical cells that are connected in series and/or parallel to provide the necessary voltage and capacity. Each cell consists of a positive electrode, a negative electrode, a separator, and an electrolyte solution [9]. Chemical reactions occur at the two electrodes, which are separated by a separator and immersed in an electrolyte solution. The separator serves as a physical barrier, preventing direct contact between the electrodes and thus reducing the risk of short circuits while allowing the transfer of lithium ions. The electrolyte solution acts as a transport medium for lithium ions between the electrodes. It contains a lithium salt that dissociates into lithium ions and other ions [10][11][12]. Figure 1 shows the operation of LIBs based on the reversible process of lithium-ion intercalation in the two electrodes. A popular option adopted for the positive electrode is lithium cobalt oxide (LiCoO 2 ). For the negative electrode, graphite is often used, with the electrolyte solution being a mixture of organic solvents. The main compound in the electrolyte is ethylene carbonate (EC), in which the lithium hexafluorophosphate (LiPF 6 ) salt is dissolved [13]. It is widely accepted that additional investigation into the chemistry and materials utilized in LIBs are necessary for the development of next-generation energy storage technologies [14,15]. This includes researching new materials for electrodes that can facilitate efficient lithium-ion insertion/extraction and exhibit higher theoretical capacities than conventional materials [16]. Moreover, researchers are exploring the incorporation of novel additives into electrolytes to enhance cycling performance, expand the temperature range of operation, and improve the electrochemical properties of the cells [17,18]. It has been observed that the current commercial electrolytes exhibit poor thermal tolerance and are prone to undesired decomposition and reactions with the electrode, mainly due to the highly reductive nature of EC. As a result, ongoing research endeavors aim to develop EC-free electrolytes to address these issues [19]. Positive electrode (e.g., LiCoO 2 ) Electrolyte (organic solvents) Figure 1. A lithium-ion battery's operation is shown schematically [20]. The negative electrode is composed of a low-potential material, such as graphite, whereas the positive electrode is composed of a transition metal that can release and receive lithium ions in its lattice. The electrolyte can be a liquid with a high ion conductivity and a wide potential window that is created by dissolving a lithium salt in an organic solvent. The performance and lifespan of a battery are impacted by a multitude of factors, including cell design, charge/discharge rates, temperature, and more. To achieve optimal battery performance, it is crucial to consider these factors simultaneously. However, one particular aspect emerges as critically important in the development of high-performing batteries, i.e., the interfacial chemistry between the electrodes and the electrolyte. The characteristics of the electrode-electrolyte interaction play a significant role in determining the battery performance and represent a focal point for advancements in battery design [21][22][23]. The aging mechanisms of a LIB are: 1.
Such phenomena cause reversible or irreversible decreases in current density, capacity, and/or power. The formation of a passivation layer at the electrode-electrolyte interface is one of the most important irreversible phenomena occurring within LIBs [27]. This passivation layer forms as a result of the initial electrolyte's thermodynamics instability, which undergoes unwanted reactions during charge or discharge cycles at the anodic or cathodic interface, respectively [28]. Uncontrolled growth of this layer leads to degradation, whereas its initial formation is desirable for a battery's stability [29]. As a matter of fact, the initial dense and intact layer restricts electron tunneling and prevents the electrolyte from further reduction, which is essential for the chemical and electrochemical stability of a battery [30]. However, its excessive growth can deplete the active lithium and electrolyte, ultimately leading to an irreversible capacity fade of the battery. The choice of battery operation condition is relevant for the formation and growth of the passivation layer. Identifying an optimal charging strategy, for instance, can help avoid overgrowth or degradation of such layer [31]. In order to maximize battery performance and lifetime, greater attention is devoted to understanding the layer formed at the anode. This is due to the characteristics of the current anode materials because their higher reactivity and redox potential make the electrolyte more unstable on the anode side than on the cathode side. In addition, the anode material undergoes larger volume expansions than the cathode, which can induce mechanical stress at the interface [32,33]. Peled et al. [34] first defined this layer in 1979 as an electronically insulating and ionically conducting solid electrolyte interphase (SEI) between the anode and the electrolyte. Currently, there is no unanimous consensus regarding the terminology used to define this layer. Some studies refer to it as the solid electrolyte interface, such as Goodenough et al. [14], whereas others use the more-accurate term solid electrolyte interphase as defined by Peled et al. [35]. This discrepancy arises from the fact that the interphase is essentially a material found at the solid-liquid (or solid-solid in solid-state batteries) interface between the anode and the electrolyte [36]. Therefore, modelers commonly use the term "interface" to refer to the surface that separates two distinct domains, whereas materials scientists prefer the term "interphase" to indicate the actual material. To remove any ambiguity, in this work, we will adopt the term "interface" to denote the surface that separates the two phases and use the term "interphase" to refer to the material itself. The SEI is made up of inorganic and organic compounds. The inorganic component of the SEI is typically composed of lithium salts and their decomposition products, which can include lithium oxides, lithium carbonates, and lithium fluorides [37,38]. The organic component of the SEI is usually composed of polymeric materials, formed by the reduction of the electrolyte organic solvent, such as alkyl carbonates, vinylene carbonate, and various electrolyte solution additives [39,40]. The SEI components can be modified by adjusting the electrolyte content in the solution. For example, Figure 2 illustrates a comparison of SEI components discovered by Yun et al. [41] for two different mixture electrolytes on silicon negative electrode using a computational approach. SEI research has been ongoing for many years, with the goal of shedding light on its complex chemical structure and the mechanisms involved in its formation [32,42]. Nevertheless, SEI is one of the most important aspects of rechargeable batteries that is still poorly understood [43]. The scientific community is constantly looking for new strategies to design batteries capable of developing sufficiently stable SEI layer to reduce irreversible capacity loss and interface resistance. The goal is to create a dense, thin, and intact passivation layer that prevents further electrolyte reductions while also allowing lithium-ion (Li + ) transport [44]. Furthermore, it is essential to minimize the formation time of such layer due to its high cost and time-consuming nature. The formation steps involved in creating the SEI can take up to two or three weeks, resulting in a notable bottleneck in the LIBs production. Therefore, by accelerating the electrolyte wetting and SEI formation processes, considerable advantages in cost and production speed can be achieved [45,46]. Given that the reversible process of lithium-ion intercalation in the electrodes is a key aspect for the long-lasting LIBs, understanding lithium-ion transport phenomena in the SEI is extremely relevant and debated in the battery field [47,48]. The SEI's extreme chemical sensitivity to oxygen and hydrogen, as well as its complex chemical composition, limits experimental observation of this transport phenomenon [49][50][51]. In general, the thermodynamics and kinetic properties of the SEI layer cannot be characterized using current experimental techniques. Therefore, to overcome experimental constraints and study mass transport events in the SEI, a numerical approach based on physical models is required (see Figure 3). Furthermore, multiscale models including chemical reactions must be developed, ranging from the electronic to the continuum level, because the mechanisms of lithium-ion diffusion are characterized by disparate time and space scales [52][53][54]. The purpose of this review is to describe the various modeling approaches that have been suggested in the literature to date to shed light on the phenomenon of lithium-ion transport through SEI in LIBs. It will be stressed that although numerous studies have been conducted over the last forty years to understand the formation of SEI [55][56][57][58], the mechanism of lithium transport in SEI has only recently received more attention. As a result, we believe that it is important at this time to summarize and comment on the numerical methods reported in the literature to investigate the state-of-art research on such diffusive phenomenon. In particular, in the chapter entitled "Modeling of Li-ion diffusion", we will review the various computational studies that have been performed, categorizing them according to spatial and time scales, namely electronic, atomistic, mesoscale, and continuum. Furthermore, we will briefly discuss the experimental challenges associated with the mechanism of interest in the chapter entitled "Experimental measure of Li-ion diffusion". Finally, in the chapter entitled "Challenges and perspectives", we will critically analyze the knowledge gaps within the methods developed so far and indicate desirable future research paths in order to reach an adequate understanding of transport phenomena within SEI.

(a) Continuum scale
Graphical illustration of ageing mechanisms in LIBs.

(d) Electronic scale
Energy pathway of Li-ion hopping in inorganic SEI components by quantum mechanical calculations.

(c) Atomistic scale
Li + transport in organic SEI components studied with MD simulations.

(b) Mesoscopic scale
The variation of the diffusion coefficient with Li concentration for layered and spinel Li x TiS 2 calculated with the kMC.

Length
Time Figure 3. Overview of research methodologies using physical models and multiscale simulations.
(a) Pastor-Fernández et al. [25] evaluated the merits of the techniques capable to identify and quantify aging mechanisms non-invasively (Copyright: Elsevier, under a license). (b) Van der Ven et al. [59] reviewed the key factors that govern Li + diffusion in electrode materials (reprinted with permission from Van der Ven et al. [59]. Copyright 2023 American Chemical Society). (c) Bedrov et al. [60] conducted an extensive investigation of structural, dynamic, and mechanical properties of organic SEI components using atomistic simulations (reprinted with the permission from Bedrov et al. [60]. Copyright 2023 American Chemical Society). (d) Ahmad et al. [48] probed the interface of inorganic SEI components to study its effects on electronic transport, Li-ion conduction, and defect formation using first principles (reprinted with permission from Ahmad et al. [48]. Copyright 2023 American Chemical Society).

Modeling of Li-Ion Diffusion
Lithium ions are the primary species that diffuse in a LIB due to their role as charge carriers in the electrochemical reactions during the charging and discharging processes. Other ions and species, such as additives or impurities such as water or metallic ions, may also diffuse to varying degrees depending on the specific battery design and chemistry [61]. Nevertheless, an efficient SEI, which is dense, thin, and intact has the advantage of limiting the transport of additives or impurities to the electrode, while allowing only the lithium ions to diffuse [62]. As a result, the ongoing research attempts to observe only the Li + diffusion mechanisms within SEI on very different spatial and temporal scales, as well as to determine how diffusion on the nano-scale and micro-scale affects the overall behavior of the battery.
Numerical models at the electronic scale can provide information on the various possible mechanism of lithium-ion diffusion, based on the characteristics and defects of the individual components of the SEI [63][64][65]. On the other hand, atomistic simulations can be used to determine the diffusion coefficient at the microscopic level [60,66,67]. Several mechanisms govern lithium mass diffusion, namely:

1.
Interstitial diffusion occurs when lithium ions in interstitial positions in the crystal lattice move through interstitial sites or gaps between atoms.

2.
Substitutional diffusion occurs when lithium ions replace other ions in the crystal lattice, resulting in the migration of Li ions within the lattice typically into vacant lattice sites. It is also known as the "knock-off" mechanism. 3.
Vacancy diffusion can happen when a Li + moves into a vacant site, leaving a vacancy behind, and another lithium ion moves into the empty site, resulting in net lithium diffusion.

4.
Hopping diffusion can occur when lithium ions jump from one lattice site to another through a series of intermediate sites, facilitated by thermal energy.
The specific mechanism or combination of mechanisms that dominate lithium diffusion depends on the material's crystal structure, temperature, and other external factors. Once this information on the nano-scale and micro-scale has been obtained, it must be determined how they can be incorporated into macro-scale models, to ensure that any micrometerscale optimization of the transport processes has a quantitative impact on the battery performance estimated using continuum models.
Fick's laws are then used to calculate the flow of lithium ions in an LIB, where the diffusion coefficient must be known [68]. Hence, an ambitious goal could be to establish a clear relationship between microscopic and macroscopic diffusion rates. Meso-scale models are used to attempt to bridge the micro-and macro-scales. For example, once the possible diffusion paths of Li ions within a crystalline structure at the nano-scale have been determined, the kinetic Monte Carlo (kMC) method can simulate how the macroscopic diffusion rate varies as the concentration of Li ions changes by simultaneously considering all possible diffusion paths [59]. Energy barriers are needed as input parameters in the kMC method in order to calculate the rates of various diffusive events. These energy barriers represent the amount of energy required by diffusing species to move from one configuration to another. They are typically estimated using nano-scale models that calculate energy barriers for relevant diffusive events as the difference in energy between an initial and final configuration along a diffusive pathway. In this respect, the nudged elastic band (NEB) method or transition state theory (TST) can be used to estimate energy barriers for kMC simulations. Current numerical models are based on empirical assumptions that may not fully capture the complex behavior of the real LIB, especially under changing conditions or during different stages of battery operation. For instance, these assumptions include simplifications of the reactions at the electrode-electrolyte interface, such as using the empirical Tafel equation to calculate the reaction rate [52]. Furthermore, the battery materials are assumed to have idealized properties, such as uniform porosity or perfect conductivity, and lithium ions move uniformly through the electrolyte without taking into account factors such as diffusion barriers or interaction with other species [69,70]. As a result, the ability to realize a multiscale model would allow for a reduction in empirical assumptions and the explicit description of the lithium-ion transport phenomenon.

Electronic Scale Models
Quantum chemical models based on electronic theories are often referred to as firstprinciples or ab initio techniques, and they play an important role in understanding the elementary reaction and transport processes within the main SEI components. The potential for accessing aspects of the electronic structure that are essential for analyzing its transport mechanisms is a major advantage of first-principles calculations. These theoretical models can explain the connections between structures and properties and allow for the investigation of mechanisms at the atomic scale and beyond. The ionic conductivity of SEI depends on the ion migration energy and the concentration of structural defects because the existence of point defects, such as vacancies and interstitials, causes ions to diffuse in solids [71]. In order to characterize the crystal structures formed in the SEI, their defects, the activation energy of the atomic hops, and more, first-principles calculations have been applied in numerous research activities [65,72,73].
Iddir et al. [74] conducted one of the first theoretical studies on lithium-ion diffusion in an inorganic SEI component using density functional theory (DFT) [75] calculations in 2010. They focused on the possible diffusion paths of the lithium ion within a lithium carbonate (Li 2 CO 3 ) crystal, as this component is thought to be one of the most stable SEI components on the electrode surface [76]. First, they investigated the equilibrium position of an interstitial Li ion in the Li 2 CO 3 crystal and then obtained the overall migration energy profile by applying the NEB [77] method to different diffusion paths. They were able to conclude from the results that the diffusion of the interstitial Li ion in Li 2 CO 3 is almost two-dimensional, with migration favored along the [0 1 0] direction. Furthermore, the migration barrier calculated for Li-ion diffusion along the favored direction is about 0.28 eV, which is similar to the migration barrier calculated for lithium-ion diffusion in graphite (i.e., 0.3 eV under the lithium excess [78]). As a result, they were able to demonstrate that a thin layer of Li 2 CO 3 on the graphite anode with the proper orientation can ensure good ionic conductivity without slowing down the transport of lithium ion in the graphite, while also improving electrolyte thermodynamics stability at the interface. Iddir et al. [74] focused on interstitial diffusion but, as Chen et al. [79] pointed out, Li ions also diffuse due to the presence of vacancies in the SEI crystal components. Hence, Chen et al. [79] used DFT to continue the investigation of Li-ion diffusion, this time attempting to shed light on the effect of vacancies in the crystal structures. Specifically, they investigated the diffusion of lithium ion into the main components of the inorganic SEI, namely the bulk of Li 2 CO 3 , lithium oxide (Li 2 O), and lithium fluoride (LiF), by removing one lithium atom from the ideal supercell. Their findings demonstrated that the diffusion of Li + is accelerated by the existence of vacancies. In particular, their results showed a slight decrease in the energy barrier for Li + diffusion along the [0 1 0] plane in Li 2 CO 3 , thus reducing the previous value of approximately 0.28 eV for interstitial diffusion [74] to the values of 0.227 eV and 0.273 eV (see Figure 4i energy profiles (c) and (d)). Moreover, for lithium-ion diffusion along the [0 1 0] plane, a substantial decrease in the energy barrier was observed, ranging from the previously reported value of about 0.60 eV for interstitial diffusion [74] to the values of 0.491 eV, 0.267 eV, and 0.324 eV (see Figure 4i energy profiles (a), (b) and (e)). The energy barriers that Li ions must overcome to diffuse in Li 2 CO 3 (see Figure 4i) and Li 2 O (see Figure 4ii) are comparable to or slightly lower than those for Li ion diffusion in graphite, ranging from 0.308 to 0.400 eV [80], although they are greater in LiF (see Figure 4iii). As a consequence, a higher energy barrier for Li-ion diffusion in SEI than in graphite may be attributed to a higher concentration of LiF in the passivation layer.   Nevertheless, as the SEI grows during battery operation, numerous defects and grain boundaries form, making the transport of Li ion to the electrode surface anisotropic and with changing paths. Ramasubramanian et al. [81] used DFT to understand and evaluate lithium-ion diffusion pathways through SEI grain boundaries. They identified grain boundaries between the inorganic SEI components that were mechanically and thermodynamically stable before running NEB calculations for a number of potential Li-ion diffusion paths. Their results revealed that the Li-ion diffusion rate was faster at the heterogeneous LiF/Li 2 O grain boundaries than at the homogeneous LiF/LiF and Li 2 O/Li 2 O grain boundaries, as shown in Table 1. The lithium-ion diffusivity may differ from their estimates because the precise structure and arrangement of the grain boundaries in the SEI during battery operation can be hardly reconstructed. It is also worth mentioning the results found by Zheng et al. [82] using DFT calculations to investigate the transport of lithium ions within the inorganic components of the SEI in the same Table 1. It should be noted that the diffusion coefficients listed in Table 1 were calculated using the Arrhenius equation: where E a is the activation energy, k B is the Boltzmann constant, and T is the temperature. The pre-exponential factor, D 0 , is dependent on the dimensionality of the diffusion (1-3), the Li + migration distance (taken from the NEB path), and the lattice vibration frequency [47]. The values in Table 1 were calculated assuming a standard lattice vibration frequency of 10 13 s −1 . However, this approximation may potentially alter the accuracy of the diffusion coefficient values. We observe some discrepancies in the activation energy values for the various Li + diffusion mechanisms within the SEI components by comparing results from the literature. For example, the E a for lithium-ion diffusion in LiF driven by vacancies was equal to 0.729 eV in the work of Chen et al. [79], whereas Zheng et al. [82] reported a value of 0.66 eV. However, it appears that there is general consensus that the knock-off mechanism, in which an interstitial Li ion moves a crystalline structure atom, and the presence of defects (such as vacancies) are the fundamental aspects underpinning the diffusion of lithium-ion in inorganic SEI components. Several environmental parameters, such as the applied electric field and temperature, have an impact on the transport of lithium ion in SEI as well. To evaluate these effects, various modeling approaches are frequently combined. The ab initio molecular dynamics (AIMD) [83] method is one way to research atomic-level diffusion mechanisms. Atomic interactions are derived from quantum mechanical calculations in this method, and then the atoms are allowed to propagate through the system according to the laws of classical mechanics. AIMD simulations can include systems of a few hundred atoms as well as short-lived phenomena lasting between 10 and 200 ps. For example, Ahmad et al. [48] used AIMD simulations to investigate the conduction of lithium ions at the interface of LiF and Li 2 CO 3 as the temperature changed. Specifically, they studied the effects of this interface on Li-ion diffusion in Li 2 CO 3 . They compared their findings to those obtained for lithium-ion diffusion in the bulk and discovered a significant reduction in the activation barrier from 0.3 eV to 0.1 or 0.22 eV due to the interface (see Figure 5). Their work highlighted the importance of interfaces in the analysis of SEI properties.

Atomistic Scale Models
Molecular dynamics (MD) [84][85][86] simulations are atomistic simulations with a substantially lower computational cost than AIMD and electronic simulations. In MD, the electronic degrees of freedom are neglected in order to avoid the computational expense of calculating the electronic density. On the contrary, atomic interactions are typically calculated using the gradient of an effective potential energy surface E(r), which is described as a sum of physically inspired interatomic interactions V. Each interaction is represented by a fixed shape function that depends only on the center of mass of the atom nuclei and the system topology V = V(r i , r j , . . . ). This collection of functions is referred to as a force field (FF). Over the past few decades, the need to model diverse atomistic systems has led to the development of several FFs formulations, which can be grouped into three main classes: classical FFs (sometimes also called Lennard-Jones FFs (LJFF)) [86], bond order-based FFs (BOFF) [87,88], and machine learning-based FF (MLFF) [89]. The LJFFs are the most widely used and numerically efficient FF in molecular dynamics simulations. This FF approximates the potential energy using bonded and non-bonded potentials. The bonded potential is typically described by mathematical functions that depend on the bond length, bond angle, and dihedral angle between atoms in a molecule. On the other hand, the non-bonded potential accounts for the van der Waals and Coulomb interactions between atoms and is a function of their relative distance, partial charge, and other electrostatic parameters. Usually, the LJFFs employ harmonic functions to describe the bond stretching, which effectively fixes the system topology during the simulation. However, this approach makes impossible the prediction of reactions within the system because it ignores key information about chemical bonding and reactions. As a result, it is unable to accurately describe reactive phenomena, such as bond breaking and formation. To address this limitation, the concept of bond order (BO) was introduced. Similar to its empirical analog introduced by Pauling [90], it is a continuum function representing the multiplicity of covalent bonds between atoms, which enables the connectivity of an atom to change dynamically during the system trajectory. Several force fields have been developed based on the BO concept, such as the reactive empirical bond order (REBO) [91], intermolecular reactive empirical bond order (AIREBO) [92], and Tersoff potentials [93], which can conveniently describe bond-breaking phenomena for specific systems. However, at present, the most used BOFF is the reactive force field (ReaxFF) proposed by Van Duin [94]. This potential incorporates many-body functions that include classical and quantum mechanical effects, allowing it to predict bond dissociation more accurately. This ability to predict bond dissociation in chemical reactions is the primary reason behind the popularity of ReaxFF, making it one of the most widely used bond order force fields. Although BOFFs enable the prediction of various reactions and phenomena, their accuracy is still far from that of AIMD calculations [87]. In recent years, a more data-driven approach known as MLFF has been proposed. Instead of approximating the energy with a set of empirically inspired functional forms, MLFFs use descriptors (or kernels) and regressors or neural networks (NN) to interpolate the potential energy surface from ab initio calculations in a high-dimensionality space [89]. Although this approach loses the physical meaning of the coefficients calculated for the force fields, the resulting prediction accuracy is comparable to AIMD simulations [89]. For all the FFs mentioned above, the coefficients are usually computed using energy and forces from DFT simulations of a set of system configurations. The number of configurations required depends on the search space dimension for optimization, which is typically less than a hundred for LJFF, less than a thousand for BOFF, and more than a dozen of thousand for MLFF.
Transport properties on greater spatial and temporal scales can be determined by using classical MD instead of AIMD. In classical MD, force fields are further classified into two types: non-polarizable, where charges are fixed at the atomic centers [95,96], and polarizable, in which charges can move relative to the atomic centers [97][98][99]. Depending on the application, one must always evaluate the force field necessary for accurate molecular dynamics. The benefit of using MD simulations is that the potential diffusion paths do not need to be known beforehand, and the impacts of several particles can be taken into account. Benites et al. [67,100,101] used molecular dynamics simulations to investigate the main components of inorganic SEI. The diffusion coefficients of lithium ion were calculated over a temperature range of 250 K to 1800 K. They also investigated the effect of an applied electric field ranging from 0.1 V/Å to 0.85 V/Å. The diffusion coefficient was calculated using the lithium ion's mean square displacement (MSD) at each temperature [102]. Their findings on the dominant mechanisms for Li-ion diffusion in inorganic SEI were consistent with most of the previous theoretical results in the literature, namely knock-off and vacancy in LiF and Li 2 CO 3 , and interstitial in Li 2 O.
In the literature, classical MD simulations have mainly been used to shed light on lithium-ion transport in organic SEI components. Di-lithium ethylene di-carbonate (Li 2 EDC) is one of the main components that has been found experimentally in SEI organic [103,104]. In 2013, Borodin et al. [105] developed and used a many-body polarizable APPLE&P force field to investigate the ionic conductivity of the Li 2 EDC component. They simulated the behavior of both the molten or disordered state and the crystal-like ordered state of Li 2 EDC. Furthermore, they evaluated the impact of the many-body polarizable force field on the Li + diffusion in the molten state of Li 2 EDC. On the one hand, they discovered that ordering slightly speeds up the lithium-ion transport in this SEI component. On the other hand, they found that the removal of polarization significantly slows down the lithium-ion diffusion phenomenon in the disordered Li 2 EDC. Figure 6 displays their findings regarding the influence of ordering and polarizability on ionic transport. They obtained activation energy for lithium-ion transport in Li 2 EDC in the range of about 0.663 eV to 0.871 eV, which is a bit higher than the activation energy of Li + in the graphite anode [106][107][108].
Later, Bedrov et al. [60] performed MD simulations to further investigate on the mechanism of lithium-ion diffusion in the amorphous and ordered states of organic SEI alkyl-carbonate compounds by using the same polarizable force field. They focused on analyzing how the properties of Li 2 EDC and Li 2 BDC changed with temperature. The diffusion coefficients of Li 2 EDC and Li 2 BDC were found to be similar at high temperatures, but as the temperature dropped, they discovered significant differences. For example, at 500 K the diffusion coefficient of Li 2 EDC was approximately 20 times smaller than that of Li 2 BDC. This comparison of the ordered and disordered states revealed that the Li-ion diffusion in the ordered state is larger (see Figure 7a). The temperature range of ordered systems was 393-500 K, and the motion of Li + appeared to be primarily confined in 2D diffusion within the layers of the organic SEI components because they are close enough to allow occasional transfers of Li + , as shown in Figure 7b. Slow SEI growth would give Li 2 EDC enough time to form ordered structures, resulting in higher ionic conductivity. In contrast, despite its higher ionic conductivity, the presence of Li 2 BDC in the SEI should be limited to optimize battery operation because they have a lower shear modulus (1.4 GPa compared to 2.8 GPa for Li 2 EDC at 393 K) and reduce the SEI's ability to prevent dendrite growth.  Li 2 EDC are compared in two scenarios. Firstly, between the ordered and molten phases, both of which are simulated using a polarizable APPLE&P FF [109]. Secondly, between the molten phase simulated using the polarizable FF and the same phase simulated using a two-body FF (TB) with the polarization turned off. Reprinted with permission from Borodin, et al. [105] (Copyright 2023 American Chemical Society).
The following year, Muralidhran et al. [110,111] identified non-polarizable OPLS-AA force field parameters to examine lithium-ion transport properties in Li 2 EDC. The structural information, such as the radial distribution functions, was in agreement with the results of the polarizable FF, indicating that their non-polarizable FF may be used to study the transport properties of Li 2 EDC. Indeed, the diffusion coefficient of Li + in the amorphous Li 2 EDC at a temperature of 333 K was found to be approximately 10 −16 m 2 /s in both MD simulations using polarizable and non-polarizable FFs. The reactive nature of lithium-ion batteries must be considered when studying the Li + transport phenomenon. As a result, reactive dynamics simulations, i.e., using a ReaxFF, should be performed rather than classical molecular dynamics simulations [112]. For instance, the reactive force field developed for Li-C-H-O atoms [113,114] was used to simulate the reduction of EC and dimethyl carbonate (DMC) solvents on a lithium metal surface [115]. The aforementioned works utilizing ReaxFF demonstrate the potential of this many-body force field in analyzing SEI properties. Nonetheless, it is important to note that the applicability of this potential is restricted to the available parametrizations in the literature, which have a limited scope of application [116]. Because the accuracy of ReaxFF depends heavily on the quality of the force field parameters, some of that can be obtained from experimental data when available, such as the equilibrium covalent bond and the dissociation energy, but the majority of them must be fine-tuned. One of the advantages of ReaxFF is that it uses fixed-shape functionals to approximate the potential energy surface. This allows for the use of a database that can be obtained with a limited number of ab initio calculations. For instance, when considering a couple of atom types, ReaxFF typically requires approximately 120-150 coefficients, depending on the number of angles and dihedral interactions to be included in the FF. Therefore, simulating around 300 well-designed configurations can provide enough data to calibrate these coefficients.
However, due to the strong nonlinearity of the equations used in the FF, only gradient-free optimization algorithms can be utilized, which can result in slow convergence. Additionally, the behavior of the resulting ReaxFF strongly depends on the region of phase space explored during the coefficient determination process. Therefore, using ReaxFF to study phenomena beyond the region explored during the ab initio calculations can lead to non-physical predictions [116]. In addition to the technical challenges of using ReaxFF, we must also consider some intrinsic issues in the formulation of the energy terms in the reactive force field developed by Van Duin. Originally developed in 2001 to model organic systems [94], the ReaxFF, over the past two decades, has undergone numerous modifications with the addition of new functional terms to enhance its applicability to a wider range of systems [87]. For example, a "lone-pair" energy term was introduced to improve the accuracy in hydrocarbon combustion [87,117,118], an improved three-body functional term was developed to handle −NO 2 group chemistry [87,119,120], and an energy term was included to account for hydrogen bonds in aqueous systems [117,121]. Therefore, the use of ReaxFF to predict the inorganic component of SEI may require a tailoring process of the energy terms in the FF, which necessitates a comprehensive understanding of the parametrization. This can limit the accessibility of this methodology to researchers outside the ReaxFF community. An alternative solution is to utilize MLFF to upscale AIMD data. Although this approach requires less knowledge about how the energy surface is interpolated, it necessitates an extensive ab initio simulation campaign to incorporate all potential configurations into the training set. Furthermore, reinforcement learning techniques [89] or the use of metadynamics may be necessary to explore the region in the phase space where the MLFF has high errors [122]. This process incurs significantly higher computational costs, even for simple systems, which is the main limit of why this new approach is not yet used in complex systems such as the anode-electrolyte interface. With the continued growth of computational power and the ongoing development of new and advanced machine learning techniques from the computer science community, it is realistic to think that soon it will be possible to describe SEI formation with the accuracy of MLFF, which is indeed is one of the main ambitious objectives of the European project BIG-MAP [8].

Bridging Micro-to Macro-Properties
In order to make meaningful comparisons between electronic/atomistic-level information and empirical observations on lithium-ion transport, it is crucial to upscale this information to larger time and space scales. The kinetic Monte Carlo [123] method is a valuable method for bridging the micro-and macro-scales [52,124]. This method is a Monte Carlo variant that allows dynamic simulations of stochastic processes. The time evolution of a system is tracked in the kMC simulations by probabilistically sampling rates of different events, such as particle hops or reactions, and updating the system state accordingly. The transition rates between different configurations of the system are typically calculated using the Boltzmann assumption [125]. According to the Boltzmann assumption, the transition rate (r) from a state A to a state B is proportional to the exponential of the negative energy difference between the two states divided by thermal energy, which is given by the product of the Boltzmann constant (k B ) and temperature (T). In the mathematical expression (Equation (2)), v 0 is the frequency of atomic vibration.
Therefore, kMC can provide insight into the dynamic behavior of complex systems using information from the lower scale (such us ∆E), allowing for the study of diffusion, reactions, and other processes that are inherently random in nature. For example, Esmaeilpour et al. [126] developed a kMC protocol to elucidate the critical parameters governing the chemical, diffusion, and aggregation reactions that influence the formation of an SEI with a thin inorganic component layer and a thick organic component layer (see Figure 8). This protocol leveraged nanoscale information obtained from first-principles calculations to learn about the ethylene carbonate reduction reactions at the graphite interface. One approach to determine ion transport properties in batteries is to perform firstprinciples calculations to extract activation energies for different diffusion paths, followed by kMC simulations to calculate diffusion coefficients numerically using stochastic ion migration simulations. For instance, Van der Ven et al. [59] used this strategy to determine the diffusion coefficient of lithium ions in electrode materials. They analyzed all conceivable pathways for lithium diffusion and used activation energy to determine their respective probabilities, utilizing first-principles calculations to gather the necessary information. Then, by performing a kMC simulation to identify the Li-ion diffusion pathway at each instant, they ensured a comprehensive sampling of the diffusion process and determined the lithium concentration dependence of the diffusion coefficient. In conclusion, their comprehensive first-principles statistical mechanical studies provided valuable insights for designing new electrode materials with a specific concentration dependence of the lithium-ion diffusion coefficient.
Continuum-scale mathematical models describe the transport of lithium ions in the SEI using macroscopic relationships of mass transport with physical parameters obtained from experiments or calculation at lower scales [127,128]. In particular, the model equations to describe the transport of the predominant specie in solids are Fick's laws, where J is the flux of lithium ions, c is the concentration, and D is the diffusion coefficient: The diffusion coefficient should be a known input parameter; however, approximations are often frequently used [129,130]. For example, in the capacity fading model for lithium-ion batteries of Ramadass et al. [131], the diffusion coefficient in the passivation layer was assumed to be the same as that in the graphite layer. Because this assumption is not totally true, it would be necessary to combine data collected at various spatial and temporal scales in order to accurately represent the transport phenomena in the SEI. Furthermore, according to Single et al. [132], the diffusion of interstitial lithium in the SEI could be the main mechanism involved in the long-term growth of this passivation layer. This conclusion would emphasize the importance of better understanding this transport phenomenon. Table 2 summarizes various computational studies, outlining the examined systems, utilized methods, and the computed properties relevant for gaining insights into the lithium-ion transport phenomenon. Ramadass et al. [131] Lithium-ion battery cell model Diffusion coefficient is an input parameter of the mass transport equation Capacity loss due to the side reaction over the negative electrode surface

Experimental Measure of Li-Ion Diffusion
The transport of lithium ions in the solid electrolyte interphase has an impact on the performance and safety of lithium-ion batteries. The ionic conductivity of the SEI layer can be affected by a variety of factors, some internal (composition and structure) and some external (temperature and voltage). Lithium-ion transport in this layer is primarily accomplished through a combination of diffusion, which is driven by the concentration gradient, and migration, which is driven by an electric field. In the existing LIBs, the SEI layer has lower ionic conductivity than the electrolyte. To improve the ionic transport properties of the passivation layer, the formation of this layer should be controlled to achieve the desired ionic conductivity [133,134]. However, the overview of the computational methods currently available in the literature has revealed gaps in the theoretical understanding of the lithium-ion transport phenomenon in the SEI [63,[135][136][137]. As a result, the design of electrodes and electrolyte giving rise to a SEI with specific ion transport properties remains an unfulfilled challenge.
One of the main limitations of current computational models is that they are created to describe the behavior of a SEI with an ideal composition [138]. For example, electronic scale models have been used to investigate lithium-ion diffusion in ideal inorganic SEI components, which are far from the true complexity of inorganic SEI, in terms of defects, imperfections, grain boundaries, and other aspects [74,79,81]. Furthermore, atomistic scale models have focused on studying Li + diffusion in the main component of organic SEI (Li 2 EDC), but the organic matrix of SEI is much more complex [60,105]. The SEI layer has a very complex chemical structure and composition, that is hard to identify experimentally due to its high sensitivity to oxygen and hydrogen. Another drawback of current numerical approaches is that they are not yet effective at bridging information across various spatial and temporal scales [139,140]. In other words, they are unable to connect data gathered at the atomistic scale with macroscopic observables such as loss capacity, voltage drop, etc. In this section, we provide a brief overview of some of the potential experimental approaches that could reveal useful information on lithium-ion transport in the SEI layer. There are several reasons for the lack of such experimental information. First and foremost, the difficulty of characterizing the structure and composition of SEI via experiment, as well as the challenge of observing its formation and growth in operando. For instance, the interpretation of Fourier-transform infrared spectroscopy (FTIR) and X-ray photo-electron spectroscopy (XPS) spectra is difficult and ambiguous because many SEI compounds are similar to one another and to the electrolyte. Furthermore, the chemistry of the SEI layer is sensitive to air exposure, which occurs frequently during sample transfer [141]. The use of other experimental techniques, such as transmission electron microscopy (TEM), is equally problematic due to the chemical reactivity and sensitivity of the SEI layer to electron beam irradiation [142]. As a result, the majority of experimental methods do not facilitate observation of the distribution and nanostructures of the organic and inorganic components in the SEI film. However, cryo-electron microscopy (cryo-EM) can be used to characterize the structure of SEI with atomic resolution [82]. There are no reactions at cryogenic temperatures, and structural and chemical information is preserved. For example, Figure 9 shows the Yuzhang et al.'s [51] observation of an SEI layer formed at the interface with a lithium metal electrode in (i) a standard electrolyte (i.e., ethylene carbonate-diethyl carbonate) and (ii) the electrolyte mixed with a fluorine-functionalized additive.
(i) (ii) Figure 9. (i) The SEI formed in a standard electrolyte resembles a mosaic structure with a heterogeneous distribution of inorganic and organic components. (ii) The SEI formed in a carbonate-based electrolyte containing 10 % volume fluoroethylene carbonate appears to have a multilayer structure with an inner layer of amorphous polymer matrix and an outer layer of large grains of Li oxide. In both panels, the atomic-resolution images are on the left, and the schematics of the observed structure formed on Li metal are on the right. This figure is from Yuzhang et al. [51] and reprinted with permission from the American Association for the Advancement of Science.
To observe the diffusion of lithium ions in the SEI layer, a technique that investigates the crystalline structure of the SEI during the LIB's operating cycle is required. For example, a non-destructive tool that can be used to study the SEI layer is XPS with sputter depth profiles. By measuring the intensity of the emitted photo-electrons, the composition of the SEI can be determined. The sputter depth profiling technique enables the analysis of the SEI's thickness and its composition from surface to depth. This non-destructive approach allows for the characterization of the SEI's composition and thickness over time, but not its atomic structure [143]. Another non-destructive technique that can be used for determining the structure of the interface at the atomic scale is in situ synchrotron X-ray scattering. Chattapdhysy et al. [49] used this technique to investigate the formation and growth of a SEI layer on a graphene surface. Then, they examined the final structure and composition of SEI by using high-resolution transmission electron microscopy (HRTEM). Their results showed the formation of a LiF crystal that grows in quantity and size during lithiation, whereas the other SEI components remain amorphous because no additional crystals have been found. Amorphous material appeared to be present between the anode surface and the LiF crystal area according to the ex situ HRTEM in Figure 10. As a result, they concluded that crystalline LiF grows randomly in the organic components matrix rather than directly on the anode surface. This result may lead to the development of a novel numerical model for describing the lithium-ion transport in the SEI, specifically representing this layer as a mosaic structure composed of randomly distributed inorganic and organic components. This evidence is in contrast to most current numerical models in which the SEI is defined as a stratified state with an inner layer of inorganic material and an outer layer of organic material. Given the lack of enough experimental analyses, it is unclear what the true final structure of SEI layer, whether mosaic or stratified, will be on the various anodic surfaces. Hence, additional experimental analyses are required in order to shed light on the lithium diffusion phenomenon in the SEI, but more importantly, close coordination and communication between the experimental and computational activity are needed. Synchroton-based operando measurements are costly, time-consuming, and frequently necessitate highly specialized geometries. Merrywaeather et al. [144] proposed the visualization of the dynamics of Li ions in the materials of the LIBs using the optical interferometric scattering microscopy technique (iSCAT). In particular, they focused on analyzing the lithiation and delithiation processes in single crystal particles of the cobalt oxide LiCoO 2 , and discovered different crystal distortions during the two processes. Using phase-field modeling, they explained the distortions observed during lithiation and delithiation in terms of variations in the diffusion velocity of lithium ions in the crystal. This serves as another evidence of the potential for collaborations between cross-disciplinary modeling and experimental efforts. On the other hand, other measurable characteristics should be predicted by theoretical models for studying Li-ion transport in the SEI. Currently, after finding the diffusion coefficient experimentally, there are no relationships that allow it to be correlated with the knowledge attained at the electronic and atomic scales. As a result, we are unable to establish which lithium-ion diffusion mechanism, as well as which SEI structure and composition, influences the measured lithium diffusion coefficient. This last piece of information is critical for guiding the formation of SEI layer and achieving adequate ionic conductivity.

Challenges and Perspectives
Understanding the lithium-ion transport phenomenon in the SEI layer requires the development of new computational approaches based on physical models. To encourage cooperation between all the related research areas, it is also important to develop open, interoperable data infrastructures for LIB materials and interfaces [8,145,146]. In light of the recent increase in interest in investigating the Li-ion transport mechanism in the SEI, additional information from theoretical and experimental activities should be carefully collected and comparable figures of merit defined for model validation. In particular, the efforts of the scientific community should focus on the development of new multiscale models to investigate the lithium-ion transport mechanism in the SEI. Additionally, experiments that provide input parameters, such as the real chemical structure of SEI, and validate output information, such as capacity loss and voltage drop, should be used to enhance these models.
We believe that future computational efforts should focus on overcoming the three main critical issues encountered in this literature review, namely:

1.
Investigating the diffusion of lithium ions in an SEI that is as realistic as possible, rather than just in the components of an ideal SEI that do not fully capture the complexity of the true SEI structure.

2.
Employing a multiscale approach that transfers information from electronic/atomic scale to macro-scale, defined, bottom-up approach. This approach is advantageous for designing battery materials because it allows for a more comprehensive understanding of underlying physics and chemistry at each length scale. In contrast, a top-down approach relies more on empirical or phenomenological relationships and does not capture the full complexity of the system at each scale.

3.
We also suggest the investigation "in vitro" and "in operando" of the SEI components to obtain a real measurement of the diffusion of Li atoms. These measures are crucial to validate all the past and future SEI models.
Concerning the first critical issue, the main challenge is to create an atomistic model of the SEI layer with morphology and composition as realistic as possible. To accomplish this, it would be necessary to run molecular dynamics simulations that faithfully reproduce the formation of the SEI components. Clearly, it is better to use reactive molecular dynamics instead of classical molecular dynamics to accurately simulate these reactive processes [41,147,148]. However, whenever the existing reactive FFs must be modified, the parametrization is computationally time-consuming and ineffective [149]. One possible solution could be to develop machine learning-based force fields [150]. The goal of this is to learn the statistical relationship between chemical structure and potential energy without prior knowledge of fixed chemical bonds or relevant interactions. The strength of these MLFFs lies in their ability to capture complex configurations and behaviors, limited only by the quality and quantity of training data [151][152][153]. Such MLFFs can utilize the knowledge gained from the training data to predict the behavior of unknown molecules more effectively than the reactive FFs that need to be customized for new systems. Although this new approach has yet to be extensively employed in complex systems such as the anode-electrolyte interface, we envision it as a tool that will enable an accurate description of SEI and its formation. Even if an appropriate force field is obtained, the SEI grows at different time and spatial scales than the initial formation of the layer; thus, it is challenging to simulate with MD. The main obstacle encountered when using MD to describe SEI evolution is the duration of the phenomenon, which is of the order of hours. To circumvent this problem, De Angelis et al. [154] developed a code that can build a SEI of a desired thickness by stochastically placing its inorganic components formed during the electrolyte's initial degradation at the electrode interface. A tool such as this enables a subsequent simulation in which the first layer of inorganic SEI is already at the electrode interface and the subsequent growth of SEI is investigated, thereby shortening the overall duration of the phenomenon. Therefore, we believe that this code could enable the study of the lithium-ion transport phenomenon in a passivation layer that is as realistic  [35]), rather than only in the individual ideal components of the SEI.
In regard to the second critical issue, communication between molecular-scale and continuum models is crucial for understanding lithium-ion transport in the heterogeneous SEI structure and allowing comparisons with experimental observables [155,156]. To fill the gap between molecular and continuum scales, mesoscale models, such as the kinetic Monte Carlo model, are frequently used. To date, they have primarily been used to shed light on the growth of SEI, with few multiscale approaches used to investigate the mechanism of lithium diffusion. The reason for this is that the interest in this transport phenomenon is relatively recent, and there are not enough data and information from experiments or from the electron/atom scale. For instance, kinetic Monte Carlo models require data and information from lower scale simulations such as DFT/AIMD. As an example, the kMC models used to study the growth of the SEI for times and lengths greater than the electron/atomic scale required data on the internal structure of the SEI and/or the activation energies of the various reaction from DFT/AIMD simulations [126,157,158]. The phase-field model is another valuable modeling approach capable of accurately describing and analyzing the behavior and evolution of complex interfaces between different phases or materials. This modeling technique has successfully explained distortions observed in experimental studies during the lithiation and delithiation processes of cathode material. Specifically, it has been employed to determine the diffusion velocity of lithium within a cathode material [144]. This modeling approach offers several advantages. Firstly, it may enable the study of diffusion a priori. Indeed, by examining different compounds and morphologies, researchers can gain valuable insights into the diffusion process. However, it is necessary to have knowledge of the diffusion of lithium in bulk conditions for each component of the SEI, which can be obtained from simplified systems using reactive atomistic simulations via well-trained MLFF. Furthermore, the phase-field model can be employed to extract the diffusion coefficient for each component of the SEI a posteriori by replicating experimental observations. In this case, a detailed spatial description of the different phases of the SEI is crucial to capture the diffusion behavior accurately. To achieve this, we suggest creating "in vitro" samples with well-defined and properly characterized SEI components and subjecting them to controlled experimental conditions. The transport of lithium ions within the crystal structures can be observed and analyzed by employing optical instrumentation. Subsequently, the phase-field modeling technique can be used to extract properties, i.e., the diffusion coefficient, which can serve as valuable inputs for continuum models.
This review delves into the transport phenomena of lithium within the SEI of batteries utilizing liquid electrolytes. However, it is crucial to emphasize that the discoveries and insights unveiled here can serve as valuable guidance for achieving a SEI with desired transport properties in future solid-state batteries (SSBs). In this direction, current studies have already been inspired by SEI in LIBs to design an artificial interface (in vitro) in SSBs with high ionic conductivity [134]. SSBs have gained significant relevance due to their notable benefits, including heightened safety features, improved energy density, extended cycle life, compatibility over a broader temperature range, suitability with high-energy materials, reduced internal resistance, and the potential for environmental benefits [159,160]. Promising avenues for SSBs development encompass quasisolid systems such as ionogels and polymer-glass composites, as well as solid oxide electrolytes. Quasisolid systems offer compelling platforms for fabricating SSBs by combining the advantages of both solid and liquid electrolytes [161]. The presence of liquid electrolytes within the solid matrix provides crucial high ionic conductivity, necessary for efficient ion transport in batteries. Solid oxide electrolytes find common application in high-temperature SSBs, enabling the utilization of high-energy-density electrode materials and enhancing overall battery performance [162]. In all of these quasisolid systems and solid oxide electrolytes, the formation and characteristics of the SEI play a pivotal role in determining the performance and stability of SSBs. In SSBs, similar to conventional LIBs, the SEI forms at the electrode-electrolyte interface as a result of side electrolyte reduction reactions. However, the formation mechanisms and the composition of the SEI in SSBs may vary due to the distinct nature of the interface. For instance, in the quasisolid system, the presence of polymer or glass matrices alters the transport properties and ion diffusion kinetics, potentially affecting the reduction reactions of the electrolyte and thus influencing SEI formation [163]. On the other hand, SEI formation in SSBs with solid oxide electrolytes significantly differs from that in batteries with liquid electrolytes, mainly due to the distinct chemical and electrochemical properties of these electrolytes. In this case, the SEI layer typically forms through the interaction between lithium ions in the electrolyte and the oxide surface of the solid electrolyte [164]. Understanding and controlling the SEI formation in these systems is crucial to achieve efficient ion transport, stable electrode-electrolyte interfaces, and long-term cycling stability in SSBs [165,166].

Conclusions
Investigating the transport phenomena of lithium ions within the SEI of LIBs is crucial for understanding the formation and growth of SEI and its impact on the performance and durability of batteries. Although several studies have examined the influence of SEI on LIBs, relatively few have focused on determining the transport properties of SEI. Modeling this mechanism is challenging, as it occurs at different length and time scales and involves various modeling limitations.
In this review, we summarize the current computational efforts to comprehend the Li + diffusion inside the known SEI components of popular LIBs (see Figure 1). Electronic and atomistic scale models provide fairly accurate information on the diffusion mechanisms based on the structural and chemical characteristics of SEI components. Continuum models can predict the macroscopic behavior of lithium-ion transport within SEI and its impact on LIBs performance, whereas mesoscale modeling links small-scale models with large-scale models to properly understand this transport phenomenon. Despite progress, obstacles remain to shed light on lithium-ion diffusion within SEI, which are also related to the difficulty of obtaining experimental data and information.
In conclusion, in order to address these challenges and properly analyze this transport phenomenon, an intense collaboration between modeling and experimental approaches as well as a larger research effort on multiscale modeling are required. Indeed, this literature review has highlighted three critical aspects: the challenge of modeling a realistic SEI, the absence of models capable of effectively coupling different scales, and the difficulty of obtaining experimental observables to validate the models. To bridge these gaps and enhance the theoretical understanding of lithium transport in the SEI, we propose several favorable directions. At the lowest scale, we advocate the use of machine learning force fields, trained by accurate DFT simulations, to simulate the SEI formation and identify electrolyte degradation products. Considering that the SEI growth occurs over longer temporal and spatial scales compared to its initial formation, a statistical approach can be employed to effectively scale up the modeling. For example, we suggest De Angelis et al.'s [154] code that involves statistically positioning the inorganic components obtained from previous electrolyte degradation simulations to build a SEI of desired thickness. Such approach can provide crucial information for mesoscale and continuum models. For instance, once a SEI with composition and morphology as realistic as possible is known, the diffusion process of lithium within it could be sampled through kMC simulations. Furthermore, it is worthy to advance experimental techniques for evaluating the transport properties within the SEI. Therefore, we propose employing optical instrumentation to observe lithium transport in "in vitro" realized SEI components. By combining theoretical prediction with experimental observations, a deeper understanding of the transport properties calculated within the SEI can be achieved. Finally, given the increasing importance of solid-state batteries, we offer a concise overview of the SEI formation in these batteries. It is worth highlighting that the insights uncovered in this study can provide valuable guidance in achieving the desired transport properties of the SEI in future SSBs.
In our vision the proposed approaches, including the machine learning force field, statistical construction of the SEI, kinetic Monte Carlo simulations, and optical instrumentation, offer promising avenues for future research, leading to improved models and a better grasp of lithium-ion transport within SEI.
Author Contributions: R.C. contributed to this work by conducting a literature screening and investigation process and by writing the original draft. P.D.A. reviewed and edited this work and he verified the final considerations and perspectives. M.F., E.C. and P.A. supervised the research. E.C. and P.A. acquired the financial support for the project leading to this publication. All authors have read and agreed to the published version of the manuscript.