We denote by y = (y 1 , . . . , y T ) the experimental data trace, that is, the total measured current I data tot, t at each time point t. In our previous work (see Figure and Table of ) we demonstrated that, to a very good approximation, experimental measurements can be assumed to be subject to normally distributed random noise, typically with zero mean and some standard deviation which we will denote by σ (we showed that typical values of the standard deviation of the experimental noise are in the range 1 to 2% of the peak current). Our mathematical model of the system, Eqs 2 to 8 above, assumes that the observed experimental data y is a function of the parameters of interest θ, θ = (E 0 , k 0 , α, C dl , R u ).