Note that in this work we choose not to screen the found reactions based on barriers but only reaction energies. These choices are possible since the number of intermediates generated by our meta-MD procedure is low enough that we can grow the reaction network quite a bit with screening based only on reaction free energies (in this case three cycles of step 1-4 above) before doing another cycle would become computationally unfeasible (months on our local cluster). Subsequently we can define the parts of the reaction networks found that seem most interesting. In this case the application is to look for potential catalysts, therefore we look for places in the reaction network where the catalyst is regenerated. Meanwhile the reaction energy for the catalyzed reaction should not be very endothermic, the threshold will depend on the accuracy of the DFT method used for calculating the reaction energies.