The converged transition states became the starting point for IRC calculations using the Euler Predictor Corrector method version implemented in pysisyphus. For cases where the IRC calculation was used as a verification that a transition state corresponds to the intended reactant and product geometries, default step and Hessian calculation frequencies were used. For production quality reaction path tracing, such as in Figure , a displacement length of 0.01 ao and Hessian recalculation at every step was applied.