Thermodynamic Circuits: Modeling chemical reaction networks with nonequilibrium conductance matrices

We derive the nonequilibrium conductance matrix for open stationary Chemical Reaction Networks (CRNs) described by a deterministic mass action kinetic equation. As an illustration, we determine the nonequilibrium conductance matrix of a CRN made of two pseudo-linear sub-networks, called chemical modules, in two different ways: First by computing the nonequilibrium conductances of the modules that are then serially connected. Second by computing the nonequilibrium conductance of the CRN directly. The two approaches coincide, as expected from our theory of thermodynamic circuits.