Optimal Control of Chemical Reaction Networks with Thermodynamic Cost
Abstract
Thermodynamic control seeks to realize desired functions at minimal thermodynamic cost, yet existing frameworks remain largely confined to slow driving regimes or state-to-state transport. Here we establish a thermodynamic optimal-control theory for Markov jump processes (MJPs) and chemical reaction networks (CRNs) that accommodates objectives over both states and reaction fluxes beyond slow-driving limit. By deriving the thermodynamic Hamilton--Jacobi--Bellman equation, we demonstrate that the optimal dynamics take the form of a perturbed generalized-gradient flow driven by the value function. A convex representation of the entropy-production rate via generalized means identifies its Legendre--Fenchel dual as the dissipation function. The resulting dynamics shares the thermodynamic structure of CRNs and admits an information-geometric characterization. We further derive a thermodynamic Pontryagin equation by recasting the control problem in terms of thermodynamic forces, providing a complementary description of optimal trajectories and clarifying connections to thermodynamic optimal transport. Our results establish a thermodynamic and geometric foundation for controlling MJPs and CRNs beyond mere transport, opening avenues for optimizing biological processes and finite-time thermodynamic operations in far-from-equilibrium regimes.
Related
- Joint Neural Posterior Estimation and Model Selection for Chemical Reaction Networks
- Mutual Linearity of Complexes in Chemical Reaction Networks
- Large-deviations theory for growing chemical reaction networks
- Generative discovery of chemical reaction networks
- Hierarchical models for large chemical reaction networks
BibTeX
@misc{kobayashi2026optimal,
title = {{Optimal Control of Chemical Reaction Networks with Thermodynamic Cost}},
author = {Tetsuya J. Kobayashi and Alix Robin and Shuhei A. Horiguchi},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.01734},
eprint = {2611.01734},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.01734},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Optimal Control of Chemical Reaction Networks with Thermodynamic Cost
AU - Tetsuya J. Kobayashi
AU - Alix Robin
AU - Shuhei A. Horiguchi
PY - 2026
DA - 2026/11/01
PB - arXiv
DO - 10.48550/arxiv.2611.01734
UR - https://arxiv.org/abs/2611.01734
AB - Thermodynamic control seeks to realize desired functions at minimal thermodynamic cost, yet existing frameworks remain largely confined to slow driving regimes or state-to-state transport. Here we establish a thermodynamic optimal-control theory for Markov jump processes (MJPs) and chemical reaction networks (CRNs) that accommodates objectives over both states and reaction fluxes beyond slow-driving limit. By deriving the thermodynamic Hamilton--Jacobi--Bellman equation, we demonstrate that the optimal dynamics take the form of a perturbed generalized-gradient flow driven by the value function. A convex representation of the entropy-production rate via generalized means identifies its Legendre--Fenchel dual as the dissipation function. The resulting dynamics shares the thermodynamic structure of CRNs and admits an information-geometric characterization. We further derive a thermodynamic Pontryagin equation by recasting the control problem in terms of thermodynamic forces, providing a complementary description of optimal trajectories and clarifying connections to thermodynamic optimal transport. Our results establish a thermodynamic and geometric foundation for controlling MJPs and CRNs beyond mere transport, opening avenues for optimizing biological processes and finite-time thermodynamic operations in far-from-equilibrium regimes.
ER -