Abstract Random walk simulations based on a molecular trajectory algorithm are performed on critical percolation clusters. The analysis of corrections to scaling is carried out. It has been found that the fractal dimension of the random walk on the incipient infinite cluster is d w = 2.873 ± 0.008 in two dimensions and 3.78 ± 0.02 in three dimensions. If instead the diffusion is averaged over all clusters at the threshold not subject to the infinite restriction, the corresponding critical exponent k is found to be k = 0.3307 ± 0.0014 for two-dimensional space and 0.199 ± 0.002 for three-dimensional space. Moreover, in our simulations the asymptotic behaviors of local critical exponents are reached much earlier than in other numerical methods.