The free energy and thermal rate coefficients of the title reaction were accurately determined with the PIP-NN PES interfaced to the Caracal program. The reaction coordinate ξ corresponds to the reactant, transition state and product of the system at specific values of 0, 1, 2, respectively. Thus, to ensure comprehensive sampling, a series of tests were conducted. Ultimately, the range of ξ from -0.05 to 1.95 was chosen. Note that the reaction coordinate at ξ = 1.875 corresponds to the product, rather than ξ = 2.0 as discussed in the method section. The umbrella sampling is performed in increments of Δξ = 0.01, resulting in 201 umbrella windows, and with the force constant of the bias potential of 0.05 (atomic units). The time step is selected to be 0. is then computed across a range of partitions (bin = 5000) with the umbrella integration method. A set of tests were conducted to ensure that the results are converged with respect to the parameters, as presented in Supplementary Material. This analysis was carried out at three distinct temperatures (472.55 K, 503.55 K and 532.95 K), for which experimental data at the high-pressure limit are available. To approximate the quantum mechanical results, convergence test with respect to the number of beads is carried out to determine the exact free energy.