The highest energy structure from minimum energy pathway search was used as the starting point for transition state optimization in all cases. Tight convergence settings were also used for all cases except the dinuclear Pd compound. For this system we used very tight settings to avoid problematic premature convergence from using small step sizes. Transition state optimization was conducted using RS-P-RFO approach with a maximum trust radius of 0.1 ao. Like the minimum energy pathway search, the trust region was iteratively decreased by a factor of 2 for convergence issues. These primarily arise because of a balance between the quality of transition state estimate, the flatness of the potential energy surface near the critical point, and the accuracy of AIMNet2-Pd. Hessians were recalculated at either every optimization step or every other optimization step to support the highest quality transition state geometries. In test cases, less frequent Hessian calculations (every 5 or 10 steps) resulted in minimal structural differences.