Equation [10] was used to account for the dependency of solid-phase diffusivity on stoichiometry, this equation is illustrated in Figure . The equation is a linear combination of Gaussian functions and captures the features of the profile accurately, meaning non-linear solid-phase diffusivity can be captured by the simulations (the fitting parameters can be found in Table ):