Arlo method [23]. The waiting time between two consecutive (reaction) events follows an is the sum exponential distribution with a mean waiting time = 1/tot, where of the rule rates. The probability for a rule to be chosen is proportional to its rate. The rulebased kinetic Monte Carlo algorithm requires several additional operations per reaction event following the selections of and e: (1) sele