All reactant and product structures were reoptimised using GFN2-xTB to verify that the structures have corresponding minima on the GFN2-xTB potential energy surface. This is the case for all reactions but reaction 16, as discussed further below. The DFT geometries for the reactant and product are used as input for the procedure described above.