Search papers, labs, and topics across Lattice.
This paper introduces a computational tool for analyzing large chemical reaction networks using a scale-splitting algorithm inspired by Wilson's renormalization group, focusing on dilute regimes. The method simplifies complex networks into effective graphs and provides hierarchical formulas for dynamics, enabling the identification of autocatalytic behavior and kinetic phases based on parameter thresholds. Key results demonstrate the tool's ability to infer kinetic rates from concentration time series efficiently, offering a robust multiscale description of chemical dynamics.
Autocatalytic behavior can be quantitatively analyzed in complex chemical networks using a novel scale-splitting algorithm that simplifies reaction pathways and dynamics.
The quest for the origin of life, especially in the metabolism-first scenario inspired by the celebrated Miller-Urey experiment, has triggered a research program dedicated to studying the emergence of complex dynamical behaviors in large chemical mixtures. Though autocatalysis, understood as the capacity of a reaction network to grow exponentially, has been recognized as a potential driver of instability and multistability, no quantitative theory has yet emerged, partly because of the lack of available kinetic data. We introduce a computational tool for large chemical reaction networks based on a scale-splitting algorithm inspired by Wilson's renormalization group. We focus on dilute regimes, where species of interest have low concentration, non-unimolecular reactions may be neglected, and the dynamics is close to linear. Depending on parameter thresholds, such networks can exhibit autocatalytic behavior. Our algorithm takes as input a network structure and outputs (1) a simplified effective graph containing the dominant reaction pathways, obtained through recursive coarse-graining; and (2) analytical formulas for the dynamics in terms of kinetic rates, called hierarchical formulas. These formulas are approximate but interpretable, accurate when scale separation is effective, and provide a reliable multiscale description of the dynamics. Their domains of validity define kinetic phases, each typically associated with a distinct pattern of chemical composition. We show on a simple example that this approach enables fast and reliable inference of kinetic rates from concentration time series. Hierarchical formulas have been implemented as a Python package and are illustrated on a simplified model of the formose reaction.
