In this work, we characterize the temperature dependence of kinetic properties in heavy atom tunneling reactions by means of molecular dynamics simulations, including nuclear quantum effects (NQEs) via Path Integral theory. To this end, we consider the prototypical Cope rearrangement of semibullvalene. The reaction was studied in the 25-300 K temperature range observing that the inclusion of NQEs modifies the temperature behavior of both free energy barriers and dynamical recrossing factors with respect to classical dynamics. Notably, while in classical simulations the activation free energy shows a very weak temperature dependence, it becomes strongly dependent on temperature when NQEs are included. This temperature behavior shows a transition from a regime where the quantum effects are limited and can mainly be traced back to zero point energy, to a low temperature regime where tunneling plays a dominant role. In this regime, the free energy curve tunnels below the potential energy barrier along the reaction coordinate, allowing much faster reaction rates. Finally, the temperature dependence of the rate constants obtained from molecular dynamics simulations was compared with available experimental data and with semi-classical transition state theory calculations, showing comparable behaviors and similar transition temperatures from thermal to (deep) tunneling regime.When heavy atom tunneling regime is attained, the free energy reaction profile tunnels below the potential energy profile due to the inclusion of nuclear quantum effects. Ring Polymer Molecular Dynamics simulations are able to provide free energy and rate constants describing this phenomenon. Here we use the Cope's rearrangement of semibullvalene as prototypical reaction. imag
