The results of the simulation are shown in Figure , and the values of the fitting parameters are presented in Table . As is evident in Figure and, the simulation provides excellent fits of the experimental data at all pH using a self-consistent set of fitting parameters and can explain the observed trends in pH-dependent limiting current and Tafel behavior across 4-orders of magnitude in kinetic current. The formation of the electron polaron in Step 1 of the reaction mechanism is identified as the RLS, which helps to explain the experimental observation of polaron suppression, whereby oxygen rapidly oxidizes any polaron formed through the chemical reaction in Step 2. We also note that, assuming a symmetry coefficient of β = 0.5 -0.6 for the outer-sphere electron transfer in Reaction [8], the observed rate coefficient of 𝑘 𝑓,2 = 12.1 cm s -1 would have a standard electrochemical rate constant of 𝑘 𝑂 2 /𝑂 2 •- 0 = 0.003 -0.011 cm s -1 which agrees very well with the observed rate constant of 𝑘 𝑂 2 /𝑂 2 •- 0 = 0.003 ± 0.005 cm s -1 for oxygen reduction to superoxide reported by McCreery et al. In Figure observed in this work compared to results reported by Bielski and Allen, showing the expected negative linear dependence but a shallower slope, suggesting the disproportionation process in this work is catalyzed.