where E a is the self-diffusion activation energy. "Best fit" Arrhenius relationships were determined using nonlinear Bayesian regression, with the full posterior distribution p(D * ) at each simulation temperature used as inputs. Only converged D * values obtained for simulations where the diffusive regime was reached within the simulation timescale (T ≥ 400 K) were used in this part of the analysis.