The model for the catalytic partial oxidation of methane on Rh(111) from Mazeau et al. was initially developed with RMG using LSRs to scale the species' binding energies from Pt(111) data. The thermodynamics on Pt(111) used in the model by Mazeau et al. were also obtained using the BEEF-vdW functional in Quantum Espresso, but the DFT parameters such as plane-wave cuto↵ energy, k-point grids etc., were di↵erent from the one used in this work, and the uncertainties were not reported. To make the LSR-estimated thermodynamic data and uncertainties consistent with the DFT calculations on Rh in this work, the species' thermodynamic parameters on platinum were recalculated using the same workflow discussed in Section 2.2.1. The polynomial parameters for each species were then modified to describe the thermodynamics based on the LSR correction from Pt to Rh, using Equation :