Thermodynamics for reduced models of chemical reactions by PEA and QSSA

The partial equilibrium approximation (PEA) and the quasi-steady-state approximation (QSSA) are two classical methods for reducing complex macroscopic chemical reactions into simple computable ones. Previous studies have primarily focused on the accuracy of solutions before and after applying model reduction. In this paper, however, we adopt a thermodynamic view and aim to establish a quantitative connection on the essential thermodynamic quantities, such as entropy production rate, free-energy dissipation rate, and entropy flow rate, between the original chemical mass-action equations and the reduced models obtained by PEA or QSSA. Our findings reveal that the reduced models by PEA and QSSA may not preserve the thermodynamic structure of the original full model, particularly when algebraic relations are used instead of differential equations. This is evident in the loss of non-negativity of the free-energy dissipation rate. To validate our results, we apply them to the Michaelis-Menten reactions as a prototype, both analytically and numerically. We anticipate that our study will inspire a reevaluation of the effectiveness of various model reduction or approximation methods from the perspective of nonequilibrium thermodynamics.