Although there are significant differences between the experimental and calculated results, we do observe a regime in the perfect crystal where apparent higher order kinetics appear. It is not caused by the holes at the surface, but by the distribution of the holes in the bulk. There is no evidence that the kinetics are governed by steps of higher order dependence on holes involving surface reactive sites. It is clear that one cannot determine a reaction mechanism based on apparent overall trends alone when many processes are involved. The simulations show that different parts of the mechanism dominate in different regimes.