Accurate description of thermodynamic, structural, and electronic properties for bulk and surfaces of ceria (CeO₂) necessitates the inclusion of the Hubbard parameter (U) in the density functional theory (DFT) calculations to precisely account for the strongly correlated 4f electrons. Such treatment is a daunting task when attempting to draw a potential energy surface for CeO₂-catalyzed reaction. This is due to the inconsistent change in thermo-kinetics parameters of the reaction in reference to the variation in the U values. As an illustrative example, we investigate herein the discrepancy in activation and reaction energies for steps underlying the partial and full hydrogenation of acetylene over the CeO₂(111) surface. Overall, we find that both activation and reaction energies positively correlate with the increase in the U value. In addition to benchmarking against more accurate theoretical methodologies, we suggest that U values are better optimized against kinetics modelling of experimentally observed profiles of products from the catalytic-assisted system of reactions.