Understanding Solid Electrolyte Interphase Nucleation and Growth on Lithium Metal Surfaces

: Experiments and theory are needed to decode the exact structure and distribution of components of a passivation layer formed at the anode surface of Li metal batteries, known as the Solid Electrolyte Interphase (SEI). Due to the inherent dynamic behavior as well as the lithium reactivity, the SEI structure and its growth mechanisms are still unclear. This study uses molecular simulation and computational chemistry tools to investigate the initial nucleation and growth dynamics of LiOH and Li 2 O that provide us with thermodynamics and structural information about the nucleating clusters of each species. Following the most favorable pathways for the addition of each of the components to a given nascent SEI cluster reveals their preferential nucleation mechanisms and illustrates different degrees of crystallinity and electron density distribution that are useful to understand ionic transport through SEI blocks.


Introduction
As modern society shifts apart from the use of fossil fuels and towards renewable sources, the improvement and development of energy storage devices become of great importance. Lithium metal is one of the best anode options due to its high theoretical capacity and negative electrochemical potential [1,2]. For next generation batteries, the use of Lithium metal anodes is crucial, especially for energy storage applications such as electric vehicles [3][4][5][6]. Nonetheless, issues related to poor cyclability and safety concerns related to dendrite growth and to unstable SEI formation have kept lithium metal anodes away from massive implementation [7].
Several strategies implemented to control dendrite growth such as electrolyte additives [8][9][10], artificial SEI layers [11,12], high salt concentration electrolytes, and others have been studied to create a dendrite-free battery [13][14][15][16]. An important conclusion from these analyses is that ion diffusion in the electrolyte and through the SEI nucleating on the metal may be key factors on regulating Li nucleation and growth [17].
The formation of a passivating layer is a hot research topic; the SEI model was first introduced by Peled in 1979 [18,19]. This SEI layer forms instantaneously as soon as the Li metal is exposed to the environment or in contact with the electrolyte due to Li metal reactivity and electrolyte instability. An SEI layer may consist of a solid or mixed solid-polymer-liquid product from the reaction between the Li metal slab and electrolytes. An ideal SEI should have Li ionic conductivity but electron blocking capabilities, thus preventing further electrolyte decomposition and parasitic reactions [20]. In addition, other properties such as high elastic strength are useful to mechanically suppress dendrites, along with proper thickness while being a compact structure [21].
Many electrolytes have been tested with Li metal anodes to improve both the efficiency and stability or to promote the formation of the desired SEI [13,14,[22][23][24]. The presence of multicomponent electrolytes and the high reactivity of lithium metal anodes lead to a wide Batteries 2021, 7, 73 2 of 14 variety of SEI components including Li 2 O, LiOH, LiF, and Li 2 CO 3 , along with organic and inorganic fragments and polymers. The SEI composition and structure strongly depend on the electrolyte components: solvents, salts, and additives [25][26][27][28][29][30]. Various aspects of the SEI have been studied such as its formation mechanisms [21,[31][32][33][34][35][36][37] and the passivation effects that SEI components exert on the Li metal [20,38,39]. Additionally, the formation of stable SEI with fluorinated layers was suggested to prolong the stability of Li metal batteries by regulating Li deposition [22,[40][41][42][43]. Several experimental and theoretical attempts have been performed to decode SEI structure and composition, but due to its strongly evolving behavior, as well as to the lithium's reactivity, such decoding is not yet accomplished [35,[44][45][46][47][48].
The SEI's chemistry and morphology may affect its properties and the Li growth mechanisms [49,50]. Although recent work with improved experimental techniques can provide more information about the SEI [51], understanding of its nucleation remains unclear. Crystalline and amorphous SEI components [52] have been found in different electrolytes studied using pair distribution function (PDF) and X-ray diffraction (XRD) techniques [53] as well as cryoelectron microscopy (cryo-EM) [54] and cryogenic transmission electron microscopy (cryo-TEM) [50], emphasizing the importance of understanding the SEI structure and its correlation with electrolyte composition.
Because SEI nucleation mechanisms can determine SEI evolution and properties, here, we investigate the initial nucleation and growth dynamics of LiOH and Li 2 O clusters and analyze the associated energetics and structural SEI nucleation information for these components. Two different mechanisms are studied: molecular addition and fragment reduction. The molecular addition aims to understand the formation energies and structures of SEI nanoclusters that were already formed SEI molecules at or near the interface are added to grow larger SEI clusters. Fragment reduction simulates the reduction of fragments dissociated from solvent or anion molecules over a Li metal cluster. The Li metal cluster is representative of an inhomogeneous plating on the electrode surface. A molecule is surrounded by the molecular electrostatic potential (ESP) that is an electrical cloud created by the molecule's electrons and nuclei [55]. The ESP at a point p(x,y,z) in the vicinity of a molecule can be defined as a force acting on a positive test charge located at point p. The calculated ESP surface of the Li cluster displays attractive and repulsive regions which are used to determine the possible sites where radicals or radical anions may attach during SEI nucleation. The obtained SEI nuclei geometries and morphologies are compared to those obtained in ab initio molecular dynamics (AIMD) simulations of precursor decomposition on a Li surface.

Growth of Lithium Hydroxide SEI Nuclei by Molecular Addition
The formation of LiOH molecules close to the surface has been observed in electrolytes with trace concentrations of water. Water-in-salt electrolytes were recently studied to identify the role of water in electrolyte decomposition. In these studies, the strong solvation formed around Li + by the water molecules causes the [Li(H 2 O) 2x )] + cation to be a 'soft' cation that, approaching the metal surface, can lead to the formation of LiOH by the cleavage of the O-H bond in the water molecules [56]. LiOH has also been detected as a component of native SEIs as a result of exposure of Li metal to the environment [20]. Figure 1 shows the growth of LiOH by the molecular addition mechanism. Assuming that the first LiOH molecule is formed in the electrolyte near the Li surface, the location of ESP surfaces are used to determine the most preferred route/location where the next LiOH molecules will be added to the SEI LiOH nuclei. The electrostatic potential is calculated by the integration of the electronic density of a given cluster or molecule. Orange on the ESP plots indicates negative values on the isosurface, an attractive behavior for the Li-ion in the LiOH molecule and a repulsive one on the OH anion. This characteristic on the ESP plots helps to narrow down the possible growth structures and predicts the next deposition site of one LiOH molecule. low an energetically favorable trend. Each change in the Free Energy of formation shown in Figure 3 illustrates the sequentially favorable addition of each LiOH molecule from 1 to 10. Free energy changes of formation ∆ shown in Figure 3 were calculated using Equation (1), where ∆ is the free energy of the SEI cluster of size and ∆ is the free energy of one LiOH molecule:  The ESP shown on the right side of Figure 1 highlights that each sequential addition of LiOH produces a symmetric cluster with positive and negative regions characteristic of the interactions between ion and dipole forces. Although not perfectly planar, the evolution of the structures in Figure 1 resembles the formation of a monolayer of LiOH which may lay almost flat on top of a surface by the fourth addition ( Figure 1d). Figure 2 showcases the beginning of nanocluster formation where the structure in Figure 1d evolves to a 3D structure from the fifth molecule addition (Figure 2a). By the time it reaches 10 LiOH molecules, the nanocluster resembles a bilayer of 5 LiOH molecules each; these layers can be identified on the side view in Figure 2f. All the structures obtained and optimized follow an energetically favorable trend. Each change in the Free Energy of formation shown in Figure 3 illustrates the sequentially favorable addition of each LiOH molecule from 1 to 10.
Free energy changes of formation ∆G f (n) shown in Figure 3 were calculated using Equation (1), where ∆G(n) is the free energy of the SEI cluster of size n and ∆G LiOH is the free energy of one LiOH molecule:  The ESP surfaces also reveal how a Li-ion may travel through such nascent structures. The orange regions would tend to attract Li-ions and could act as channels for Li diffusion, whereas the yellow regions may promote a knock-off transport mechanism. Moreover, note that electric field gradients may change the locations of these regions, as shown in previous work [57], thus defining alternative diffusion pathways. The ESP surfaces also reveal how a Li-ion may travel through such nascent structures. The orange regions would tend to attract Li-ions and could act as channels for Li diffusion, whereas the yellow regions may promote a knock-off transport mechanism. Moreover, note that electric field gradients may change the locations of these regions, as shown in previous work [57], thus defining alternative diffusion pathways.   Figure 4 shows an alternative nucleation route, via the addition of an OH radical that oxidizes the neutral Li6 nanocluster and forms LiOH. The main finding in this fragmental growth mechanism is the interaction that the OH − radical has with Li metal atoms. This growth mechanism would be consistent with the addition of OH radicals formed by electrolyte side reactions to the Li metal surface. Every OH fragment interacts with three lithium metal atoms that are initially equidistant between them when OH has added the three Li atoms and OH interact, forming a triangular pyramid ( Figure 4a). As more OH fragments are added, the initial Li cluster becomes separated. However, the ion-dipole LiOH molecules formed to connect via the initial Li atoms (Figure 4b Figure 4 shows an alternative nucleation route, via the addition of an OH radical that oxidizes the neutral Li 6 nanocluster and forms LiOH. The main finding in this fragmental growth mechanism is the interaction that the OH − radical has with Li metal atoms. This growth mechanism would be consistent with the addition of OH radicals formed by electrolyte side reactions to the Li metal surface. Every OH fragment interacts with three lithium metal atoms that are initially equidistant between them when OH has added the three Li atoms and OH interact, forming a triangular pyramid ( Figure 4a). As more OH fragments are added, the initial Li cluster becomes separated. However, the ion-dipole LiOH molecules formed to connect via the initial Li atoms (Figure 4b-d). The ESP shown in Figure 4 illustrates the electrostatic interaction behavior of the OH fragment with the Li metal cluster. An example is Figure 4b,d: When the OH is interacting with only two Li atoms, there is a large orange area in the ESP, indicating negative values on the isosurface and attractive behavior for incorporation of an additional Li-ion. The free energy profile is shown in Figure S1a. Figure S1b shows the relative Gibbs Free Energy of formation for the addition of OH to the Li metal cluster and highlights the favorable addition of the OH to the metal cluster even in the presence of extra negative charge on the metal cluster surface. Each OH added to the metal cluster interacts with three Li metal atoms, and the addition is favorable.

Lithium Hydroxide
To test the calculated structures, we carried out AIMD simulations where LiOH was generated by water dissociation and reaction on a Li surface. Figure 5 shows snapshots of an AIMD study in which H 2 O molecules are placed on top of a lithium metal surface. After 80 ps of reaction at 330 K, the system was progressively cooled by lowering 10 K every 5 ps. The cooling rate is shown in Figure 5d. Water molecules react with the lithium metal and form LiOH and H 2 as shown in Figure 5b. Figure 5c illustrates similarities of the OH-Li interaction between the AIMD water-induced growth of LiOH on Li metal surface and the OH addition to Li clusters obtained with quantum chemistry calculations in small clusters shown in Figure 4. After the cooling process, the structure becomes more organized as shown in the increasing intensity of the main peak (at~1.97 Å) in the g LiO (r) radial distribution function in Figure 5e. This distance is in the range of the Li-O distances in the cluster shown in Figure 4 (1.89 to 1.98 Å). Figure S2 shows the similarities between the LiOH crystalline structure and the AIMD structure after the cooling process, the most common distances between Li-O, Li-H, and O-H are similar to the ones found in the crystalline structure. ments are added, the initial Li cluster becomes separated. However, the ion-dipole LiOH molecules formed to connect via the initial Li atoms (Figure 4b-d). The ESP shown in Figure 4 illustrates the electrostatic interaction behavior of the OH fragment with the Li metal cluster. An example is Figure 4b,d: When the OH is interacting with only two Li atoms, there is a large orange area in the ESP, indicating negative values on the isosurface and attractive behavior for incorporation of an additional Li-ion. The free energy profile is shown in Figure S1a.   Figure S1b shows the relative Gibbs Free Energy of formation for the addition of OH to the Li metal cluster and highlights the favorable addition of the OH to the metal cluster even in the presence of extra negative charge on the metal cluster surface. Each OH added to the metal cluster interacts with three Li metal atoms, and the addition is favorable.
To test the calculated structures, we carried out AIMD simulations where LiOH was generated by water dissociation and reaction on a Li surface. Figure 5 shows snapshots of an AIMD study in which H2O molecules are placed on top of a lithium metal surface. After 80 ps of reaction at 330 K, the system was progressively cooled by lowering 10 K every 5 ps. The cooling rate is shown in Figure 5d. Water molecules react with the lithium metal and form LiOH and H2 as shown in Figure 5b. Figure 5c illustrates similarities of the OH-Li interaction between the AIMD water-induced growth of LiOH on Li metal surface and the OH addition to Li clusters obtained with quantum chemistry calculations in small clusters shown in Figure 4. After the cooling process, the structure becomes more organized as shown in the increasing intensity of the main peak (at ~1.97 Å) in the gLiO(r) radial distribution function in Figure 5e. This distance is in the range of the Li-O distances in the cluster shown in Figure 4 (1.89 to 1.98 Å). Figure S2 shows the similarities between the LiOH crystalline structure and the AIMD structure after the cooling process, the most common distances between Li-O, Li-H, and O-H are similar to the ones found in the crystalline structure.

Lithium Oxide
Oxygen atoms are generated by the decomposition of products of the reduction of organic carbonates such as ethylene carbonate (EC), propylene carbonate (PC), fluoroethylene carbonate (FEC), but also from anions or other additives, and Li oxides are also very often found as native layers on Li surfaces. The sequential addition of oxygen atoms to the Li 6 metal cluster, along with its ESP, is shown in Figure S3. The structure of a Li 2 O cluster obtained after the sequential addition of oxygen atoms to Li 6 reveals interesting similarities with the crystalline bulk structure. Figure 6 shows Li 2 O cluster when added oxygen atoms reach the stoichiometric ratio.

Lithium Oxide
Oxygen atoms are generated by the decomposition of products of the reduction of organic carbonates such as ethylene carbonate (EC), propylene carbonate (PC), fluoroethylene carbonate (FEC), but also from anions or other additives, and Li oxides are also very often found as native layers on Li surfaces. The sequential addition of oxygen atoms to the Li6 metal cluster, along with its ESP, is shown in Figure S3. The structure of a Li2O cluster obtained after the sequential addition of oxygen atoms to Li6 reveals interesting similarities with the crystalline bulk structure. Figure 6 shows Li2O cluster when added oxygen atoms reach the stoichiometric ratio. The Li6 metal cluster structure becomes completely modified throughout this process, and the cluster formed has similar atomic distances and angles to the crystalline structure of the Li2O (110) facet. Additional studies performed in a larger Li metal cluster (15 atoms) shown in Figure 7 highlight the preferential formation site of Li2O. Figure 7ad show the optimized structures for the initial structure until the addition of the third oxygen atom. For each new oxygen addition, four to five different sites were tested (not necessarily closer to the previously deposited atoms); however, the first three lowest energy oxygen atoms were found located closer to each other (Figure 7d). In the next addition, shown in Figure 7e, the top of the cluster rearranges creating a Li2O cluster separated that the Li metal region, and the fourth O addition takes place in that Li metal region. In the subsequent additions (Figure 7f,g), a new Li2O cluster starts to be formed following new O incorporations, and in Figure 7h, we observe that the top and bottom Li2O clusters try to rearrange in a unified structure. Note also the ESP diagrams with alternating positive and negative regions. As discussed with LiOH, as the SEI is being formed by the fusion of nascent clusters (in this case Li2O), the ESP provides an idea about the possible Liion migration patterns. Looking at Figure 7b-h, these alternating positive and negative regions may offer a favorable migration channel via hopping or knock-off mechanisms. From Figure 7, the presence of large yellow regions suggests that the Li-ion may migrate easily or become reduced at the periphery of these clusters. However, as the oxide forms, the orange alternate with yellow regions, and that describes a more complex surface for For each new oxygen addition, four to five different sites were tested (not necessarily closer to the previously deposited atoms); however, the first three lowest energy oxygen atoms were found located closer to each other (Figure 7d). In the next addition, shown in Figure 7e, the top of the cluster rearranges creating a Li 2 O cluster separated that the Li metal region, and the fourth O addition takes place in that Li metal region. In the subsequent additions (Figure 7f,g), a new Li 2 O cluster starts to be formed following new O incorporations, and in Figure 7h, we observe that the top and bottom Li 2 O clusters try to rearrange in a unified structure. Note also the ESP diagrams with alternating positive and negative regions. As discussed with LiOH, as the SEI is being formed by the fusion of nascent clusters (in this case Li 2 O), the ESP provides an idea about the possible Li-ion migration patterns. Looking at Figure 7b-h, these alternating positive and negative regions may offer a favorable migration channel via hopping or knock-off mechanisms. From Figure 7, the presence of large yellow regions suggests that the Li-ion may migrate easily or become reduced at the periphery of these clusters. However, as the oxide forms, the orange alternate with yellow regions, and that describes a more complex surface for ion migration. Given that the yellow regions are still present in the inner part of the nascent Li 2 O clusters (see for example Figure 7f), it is difficult to think that Li ions will migrate through vacancies inside the clusters. Instead, it is more plausible that ionic migration takes place through the surface of the nascent clusters. This also suggests that the morphology of an optimum SEI layer should be one where small clusters are dispersed in a matrix, allowing ionic migration through their surfaces. ion migration. Given that the yellow regions are still present in the inner part of the nascent Li2O clusters (see for example Figure 7f), it is difficult to think that Li ions will migrate through vacancies inside the clusters. Instead, it is more plausible that ionic migration takes place through the surface of the nascent clusters. This also suggests that the morphology of an optimum SEI layer should be one where small clusters are dispersed in a matrix, allowing ionic migration through their surfaces.  As with LiOH, we performed AIMD simulations to investigate the structure that would be formed by an SEI reaction. Figure 8 depicts the geometries obtained when O 2 is in contact with a Li metal surface. The distances and angles found after Li oxidation and formation of a Li 2 O surface film are similar to Li 2 O clusters obtained with quantum chemistry calculations shown in Figure 6.
Batteries 2021, 7, x FOR PEER REVIEW 9 of 14 As with LiOH, we performed AIMD simulations to investigate the structure that would be formed by an SEI reaction. Figure 8 depicts the geometries obtained when O2 is in contact with a Li metal surface. The distances and angles found after Li oxidation and formation of a Li2O surface film are similar to Li2O clusters obtained with quantum chemistry calculations shown in Figure 6. The AIMD simulations were carried out for 35 ps at 330 K. After this time, the system was progressively cooled at the same cooling rate used for the LiOH system, as shown in Figure 8d. After the cooling process, the structure becomes more organized. Although the structural change is small compared to that of LiOH, the red line in Figure 8e and the increasing intensity of the main peak in the radial distribution function can be observed. Comparison with the crystalline structure is highlighted in Figure S4. The main interactions between Li and O and between O and O are reasonably reproduced by the AIMD The AIMD simulations were carried out for 35 ps at 330 K. After this time, the system was progressively cooled at the same cooling rate used for the LiOH system, as shown in Figure 8d. After the cooling process, the structure becomes more organized. Although the structural change is small compared to that of LiOH, the red line in Figure 8e and the increasing intensity of the main peak in the radial distribution function can be observed. Comparison with the crystalline structure is highlighted in Figure S4. The main interactions between Li and O and between O and O are reasonably reproduced by the AIMD structure, which follows similar patterns as the clusters in Figures 6 and 7. The comparison of Figures 5e and 8e shows a different extent of rearrangement after the cooling process. This is credited to the differences in crystalline structure and bonding. LiOH has been defined as being at the intersection between ionic bonding and hydrogen bonding [58], where the H-bonding characteristics are temperature-dependent, in contrast to the case of a pure ionic bonding crystal such as Li 2 O. Thus, the T-dependent behavior of H-bonding explains the substantial change in the radial distribution function shown in Figure 5e, in contrast to the almost unchanged structure of Li 2 O (Figure 8e) under the same cooling process. When the cell cools down the OH fragment and residual water molecules become less mobile and more organized resulting in the peak of the RDF increasing. In Li 2 O the Li-O bonds are more rigid and the change after cooling is much less significant.
The advantage of analyzing in parallel the cluster calculations of the growing structures and those of dynamic (amorphous) structures generated in AIMD by decomposition of precursor molecules, is that allows a cross-check of the structure of the nascent clusters and simultaneously provide insights into possible amorphous structures that could be present in certain SEI layers.

Computational and System Details
The methodology used in this study includes quantum chemistry models and calculations performed with Gaussian09 [59] and GaussView6 for visualization. To obtain the ground state structures geometry optimization and frequency calculations were carried out and the thermodynamic properties were based on the structure obtained after energy minimization. The electrostatic potential obtained from the optimized structures is used to determine the most preferred growth location for a molecule or atom adding into the nucleating SEI. The B3PW91 hybrid exchange-correlation functional that was successfully applied to many chemistries [25,[60][61][62] was used along with the 6-311++G(d,p) basis set for geometry optimization and evaluation of the molecular orbitals for LiOH and Li 2 O models. For the cluster calculations, we use an implicit solvent Tetrahydrofuran (THF) with a value of 7.2 for the dielectric constant to emulate 1,2-dimethoxyethane (DME). This implicit solvent model approximates the surrounding electrolyte environment. The implicit solvation is implemented via the Polarizable Continuum Model and the IEFPCM integral equation formalism [63,64] that follows a self-consistent reactive field technique.
The Vienna ab initio Simulation Package (VASP) [65][66][67] was employed to study the behavior of SEI nucleation on a slightly larger scale, taking advantage of the periodic boundary conditions to represent an infinite surface. Projector augmented wave (PAW) pseudopotentials from the VASP database were used to describe electron-ion interactions [68,69]. As exchange-correlation functional, we employed the Perdew-Burke-Ernzerhof generalized gradient approximation (GGA-PBE) [70]. The optimized energy cut-off for the plane-wave basis expansion was 400 eV. The ions were relaxed to their ground state via a conjugate-gradient technique. A width of 0.1 eV was used for Gaussian smearing. A 2 × 2 × 1 k-point Monkhorst-Pack [71] mesh was employed for dynamics calculations. For electronic self-consistent iteration and forces ionic relaxation, the convergence criteria were set to 10 −4 and 0.02 eV, respectively. To study the nucleation of LiOH and Li 2 O by reaction of precursors (water and O 2 ) on a Li surface, AIMD simulations were used. A canonical NVT ensemble was used with a Nose thermostat [72,73], the initial temperature of the simulation is 330 K, and a progressive cooling is performed by lowering the temperature 10 K every 5 ps until a final temperature of 270 K is achieved using a time step of 1 femtosecond with a damping parameter of 0.5.

Conclusions
In this paper, we studied the possible nucleation and initial growth mechanisms for the SEI components Li 2 O and LiOH. Our approach assumes that a Li metal cluster represents a portion of the surface where inhomogeneous plating is taking place. This allows the study of SEI structures in lithium-rich environments simulating the initial stages of electrolyte decomposition and SEI formation. An additional model is studied in which already formed fragments of SEI in the electrolyte-electrode interface are joined together to form larger SEI clusters. We use ab initio methods in cluster calculations to determine the preferential addition sites for OH or O radicals, and the structures are compared to those induced by precursors on Li metal slabs, as obtained in the AIMD simulations. Studying the nucleation mechanism of different SEI components allows for an improved fundamental understanding of the SEI properties and a prediction of the possible SEI structures that the SEI components may have, as these reactions continue during battery cycling. The ability to characterize SEI structures and growing morphologies, along with their chemical characteristics, is crucial to understand ionic and electronic transport through the SEI. Moreover, the research presented here provides thermodynamics and kinetics of nucleation and growth that are useful for the development of coarse-grained models (for example, Kinetic Monte Carlo) and also for force field development for classical molecular dynamics simulations. In addition, a crucial need in battery research is the understanding of ion transport through SEI layers, using methods such as constrained-molecular dynamics [74,75] that require knowledge of the detailed structure and morphology of the growing phase as shown in this work.
In summary, some of the emerging patterns of the SEI growth include the highly symmetrical and crystalline bilayer of LiOH if a molecular addition mechanism happens. Differently, OH fragments interact with the Li metal atoms located in a triangular pyramidal geometry with the OH at the top being equally distant of all three Li atoms, forming a triangular structure. The calculated Li 2 O SEI nanocluster reveals similarities with the crystalline structure of the Li 2 O (110) facet, and the initial Li metal cluster becomes completely modified in the sequential addition of Oxygen atoms. However, if partial oxidation of a Li metal cluster happens, the SEI will nucleate in one site of the Li metal creating pure Li 2 O before oxidizing any other part of the metal cluster. In a larger system such as the one shown in Figure 8c, the Li 2 O starts to form an antifluorite structure characterized by Lithium located at the tetrahedral interstitial sites surrounded by Oxygen. However, in smaller clusters, the O atoms are located in the tetrahedral site (Figure 7) or in the middle of three Lithium atoms resembling a tri-capped trigonal prismatic coordination, a structure usually found at higher pressure [76].
The study of the possible mechanisms of SEI nucleation and growth is of great importance towards the elucidation of the SEI structure and the Li cation deposition and reduction sites and to provide the basic understanding of the structure of the nucleating SEI. The analysis of the electrostatic potential surfaces allows a preliminary discussion on possible cationic migration patterns through the SEI.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.