The quality of tuning was confirmed by comparing agreement of the different C-rates at 25 °C, see Figure . Since the parameters were tuned based on the 1C case there is excellent agreement here, the voltage profiles for 0.5C and 2C have disagreement for the final voltage during relaxation. However, this could be improved by adjusting the electrode diffusivities explicitly for these cases-tuning is likely needed depending on the operational conditions being used in the simulation. Next, we compare the data for the 1C case at various temperatures, see Figure . This allows us to observe whether the temperature dependencies of the electrochemical parameters have been mapped adequately. The comparison at 10 °C and 25 °C illustrates good agreement for both the voltage and temperature profiles. However, the agreement between experimental and simulated data is worse at 0 °C, this is due to this temperature being outside the range for measured parameter values-there is likely an interplay of different effects meaning that the Arrhenius relationship cannot be applied for the entire operating temperature range. The heat transfer coefficient was adjusted manually to a value of 15 W m -2 K -1 , which is within the expected range of values.