Abstract
For ethylene carbonate on the (100) surface of lithium, we calculate the adsorption energy in two binding motifs and the barrier height for a ring-opening decomposition reaction. We validate a scheme for producing results in the thermodynamic limit by correcting results obtained on finite lithium clusters containing only 40–100 atoms, which enables the use of hybrid density functionals, random-phase approximation, and correlated wavefunction theories, such as coupled-cluster theory and auxiliary-field quantum Monte Carlo. We find that the high-level theories agree to within 2–5 kcal/mol and can, therefore, serve as benchmarks for more affordable methods. Using our reference data, we demonstrate that generalized gradient approximation functionals, such as PBE, are not sufficiently accurate for reaction barrier heights, and we identify ωB97X-V as an especially promising functional for the interfacial chemistry of electrolyte solvents at lithium metal anodes.