Accurate description of thermodynamic, structural, and electronic properties for bulk and surfaces of ceria (CeO2) 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 CeO2-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 CeO2(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.
