Linear conjugacy of chemical kinetic systems

TL;DR

This paper develops a computational method to determine linear conjugacy in chemical kinetic systems, extending existing criteria to broader classes and enabling dynamic equivalence transformations.

Contribution

It introduces a general computational approach for constructing linear conjugates of RID chemical kinetic systems, extending the Johnston-Siegel Criterion to complex factorizable systems.

Findings

Extended JSC to CF RID systems

NF RID systems can be transformed into CF systems

Linear conjugacy can be generated for all RID systems

Abstract

Two networks are said to be linearly conjugate if the solution of their dynamic equations can be transformed into each other by a positive linear transformation. The study on dynamical equivalence in chemical kinetic systems was initiated by Craciun and Pantea in 2008 and eventually led to the Johnston-Siegel Criterion for linear conjugacy (JSC). Several studies have applied Mixed Integer Linear Programming (MILP) approach to generate linear conjugates of MAK (mass action kinetic) systems, Bio-CRNs (which is a subset of hill-type kinetic systems when the network is restricted to digraphs), and PL-RDK (complex factorizable power law kinetic) systems. In this study, we present a general computational solution to construct linear conjugates of any "rate constant-interaction function decomposable" (RID) chemical kinetic systems, wherein each of its rate function is the product of a rate constant and an interaction function. We generate an extension of the JSC to the complex factorizable (CF) subset of RID kinetic systems and show that any non-complex factorizable (NF) RID kinetic system can be dynamically equivalent to a CF system via transformation. We show that linear conjugacy can be generated for any RID kinetic systems by applying the JSC to any NF kinetic system that are transformed to CF kinetic system.