For a simple kinetic model with a small set of surface species, the direct computation of thermophysical properties via statistical mechanics methods at each time step and coverage would be tenable. However for a larger microkinetic model with multiple surface species that are relevant for technical systems, such direct thermochemistry calculations can easily accumulate computational burden. Instead, analytical forms of thermodynamic properties to reduce computational load are desired. Thus, we describe our schemes to parameterize thermodynamic coverage dependence. We have implemented a coverage-dependent thermochemistry module in Cantera. In Cantera, the user provides an input file in the YAML format (an easy, human-readable syntax for representing data) that contains the essential information: the names of the species, their thermophysical properties (typically represented as NASA polynomials, but other representations are possible), and elementary rate constants in the forward direction (typically in Arrhenius format). The reverse rate constant is computed internally assuming microscopic reversibility, k rev = k fwd /K, by first computing the free energy of reaction and thence the equilibrium constant K. This approach ensures that the mechanism is always thermodynamically consistent, even when coverage effects are included. Coverage-dependent thermophysical properties are calculated according to Equations 4-11. The coefficients in Equation -11 can now be included as input parameters in the mechanism YAML file, thereby enabling the evaluation of coverage-dependent thermophysical properties via Cantera for reactor simulations. Full details of the implementation are documented on the Cantera website ().