In each SEI component, we have investigated native charged point defects, including vacancies (e.g., Li vacancy in charge state 1, V -1 Li ), interstitials (e.g., Li interstitial in charge state +1, Li +1 i ), and antisites (e.g., Li on O site in charge state +1, Li +1 O ). Vacancies and antisites at each Wycko↵ site were treated as distinct defects. For interstitial defects, the plausible sites were determined by a Voronoi tessellation scheme as implemented in pylada-defects. At least five interstitial sites were identified for each of SEI component. The lowest-energy interstitial site was determined from the total energy of the interstitial configurations in the neutral charge state. A broad set of charge states were considered for each defect type, typically q = -3, -2, -1, 0 , +1, +2, and +3. Additional charge states were considered where necessary. The defect formation energy E D,q of a point defect D in charge state q is calculated as: