Microkinetic modeling relies on the numerical steady-state solution of the ordinary differential equations (ODEs) systems defining the material balance for the species in the CRN. The default reactor model in CARE is a zero-conversion differential plug-flow reactor, which enables the extraction of macroscopic reaction rates and selectivity directly comparable to experiments, steady-state surface coverages, and rate-determining steps and species. The kinetic coefficients of the elementary reactions in the network are evaluated via transition state theory for all forward reactions except the adsorption steps, the coefficients of which are given by the Hertz-Knudsen equation . The coefficients of the reverse reactions are obtained by imposing thermodynamic consistency (see Supplementary Note 8 for implementation details).