Our simulation scheme focuses on the physically correct timescales for both reactive and diffusive processes. In situations where multiple reactions compete, reaction paths with lower barriers are favoured over more activated paths. While species that can be produced are limited to those included in the reaction network, those that do form are ultimately dictated by activation barrier. When using uncorrected sampling, populations are more an artefact of the number of candidates found and appearances in the reaction network.