Despite the importance of lateral interactions in catalytic processes under industrial conditions, the incorporation of coverage-dependent parameters in microkinetic mechanisms remains a persistent challenge for the modeling community. Coverage dependence, if it is addressed at all, is typically included as a linear correction to the activation energy in an Arrhenius expression of an elementary rate constant, E a (θ) = E 0 a + i e i θ i , where E a (θ) is the coverage-dependent activation energy, θ is the array of surface site fractions, E 0 a is the activation energy in the low-coverage limit, and e i is the "correction" to the activation energy due to the coverage of species i, which could be positive or negative . However, arguably the single most important impact of lateral interactions is to the binding energy of adsorbates; consequently, a kinetic parameterization is an awkward solution to a problem that is thermodynamic in nature. More importantly, though, this kinetic parameterization can lead to thermodynamic inconsistency, in which the state properties enthalpy and entropy are no longer conserved . One aspect of this challenge is the absence of a common computational framework. CatMap, for example, includes features for coverage-dependent binding energies (and thus enthalpies of adsorption) . However, it's primary focus is on the steady-state simulations of catalyst surfaces and the evaluation of rate or selectivity maps when combined with linear-scaling relationships. It is not intended to simulate reactor models or even complex networks of reactors, especially for transient systems. At present, there is no agreed-upon standard or formulation for coverage dependence in microkinetic mechanisms.