The results reveal one reaction that was unexpectedly hard to find, namely R9 , which can be classified as a roaming reaction A possible reason for the particularly rare identification is the extensive traversal of large parts of the PES by roaming reactions due to the separation and reorientation of the radical parts. The roaming nature and the stopping technique at the DS make it improbable for the system to be in the spot where separation and reorientation occur at the same time, i.e. identifying other similar reactions is a more probable outcome. However, PESmapping is capable of shortening the reaction coordinate by biasing the formation of a C -C bond and identifying the reactions after only a short time. It should be noted that triradical products are unlikely to be stable intermediates of any process, and the appearance of R9 in the exploration result may be an artifact of ReaxFF. Nevertheless, these observations are applicable to all roaming reactions. Both reaction exploration methods search for reactions in parallel independent RMD simulations. CTY-TAD simulated for a total of 80.9 core hours in 11 parallel runs, while PESmapping simulated a total of 317.6 core hours of RMD in 999 parallel runs. As CTY-TAD identifies a greater number of reactions, it is more efficient in doing so as compared to PESmapping. The major source of time loss in the PESmapping simulations can be attributed to the runs that are unsuccessful in producing a reaction path. In this study, a total of 111 PESmapping runs simulated until the time limit of two hours without reporting a reaction, which represents almost 70% of the total simulation time. Nevertheless, some PESmapping runs report reactions well after one hour of simulation, which makes it challenging to reduce this time limit by the right amount and improve the efficiency of the exploration method. Fig. shows the average computing times of CTY-TAD for reactions with a TS and without. For reactions with a TS, the computation time correlates roughly exponentially with the activation barrier. For barrierless reactions, the correlation is different. As explained in more detail in the following Section 3.4, dissociation energies need not be fully reached by the RMD simulation for CTY to identify the reactive event.