balance between accuracy and efficiency for all dynamically relevant regions. We finally choose the explicitly correlated coupled cluster single, double, and perturbative triple excitations with the augmented F12 explicit corrections, using the aug-cc-pVTZ basis set of Dunning, 40 CCSD(T)-F12a/AVTZ for short, as the theoretical method for all electronic structure calculations. CCSD(T)-F12a/AVTZ has the ability to efficiently obtain results as accurate as those obtained at the CCSD(T)/AV5Z level, which has been widely used to develop PES for various systems, reactive or non-reactive. The T1 diagnostic values for most of the sampled points are consistently below 0.03, thus, a single-reference-based methodology seems to be appropriate for this reaction. All ab initio points are performed using the MOLPRO 2020 software package.