where = [Li ] + [Li ] and E f is an effective vacancy-formation energy. To fit this model, we use the mean observed defect concentrations from our simulations at each temperature as input, and fit using nonlinear least-squares. For consistency, we restrict ourselves to simulations where the diffusive regime was reached within the simulation timescale (T ≥ 400 K) in this part of the analysis.