Using the default parameter set (k push = 0.05, α = 0.3 and s = 0.8) we find 20 of the 26 reactions. For three of the reactions (R3725, R7207, and R8701) the meta-MD failed to generate the corresponding product structures. Another reaction (R2514) is eliminated due to having a GFN2-xTB reaction energy of 53 kcal/mol, which is significantly higher than the corresponding ωB97X-D3/def2-TZVP-value of -3 kcal/mol. The final two reactions (R4612 and R9011) are eliminated due to high estimated barrier heights of 48 and 47 kcal/mol, respectively -considerably higher than the corresponding ωB97X-D3/def2-TZVPvalues of 26 and 30 kcal/mol. The actual barrier heights at the GFN2-xTB level are 42 and 45 kcal/mol, which indicates that the problem lies primarily with the GFN2-xTB method itself and not the barrierestimation method. If we perform additional meta-MD simulations with a slightly different parameter set (k push =0.03, α = 0.7, s = 0.6) we locate R3725 and R7207, but R3725 is subsequently eliminated due to a high reaction energy (42 kcal/mol) compared to a DFT value of -21 kcal/mol. So, using two set of parameters, the meta-MD simulations identify 25 of the 26 low-barrier reactions reactions found by Grambow et al., but four of these are subsequently eliminated due to overestimated reaction energies and estimated barrier heights by GFN2-xTB. Only R8701 (Figure ) is not found by the meta-MD, presumably due to a combination of high barrier and reaction energy (27 and 19 kcal/mol at the DFT level) plus a relatively small change in structure on going from reactants to products, resulting in a weak biasing potential.