Candidates are identified and reactions executed using a pre-defined list of reactions. Separation of the reactive and dynamical steps in each cycle means that classical non-reactive MD can be used. Control over the sampling of candidates and length of MD being run is maintained, and reactions are guaranteed to be chemically viable. Similar methods have also been used to model the curing of composites and could, in principle, be applicable to any reactive set of molecules.