Simulating Chemical Kinetics Without Differential Equations: A Quantitative Theory Based on Chemical Pathways

Stochastic paths controlling speed and dissipation

Natural processes occur in a finite amount of time and dissipate energy, entropy, and matter. Near equilibrium, thermodynamic intuition suggests that fast irreversible processes will dissipate more energy and entropy than slow quasistatic processes connecting the same initial and final states. For small systems, recently discovered thermodynamic speed limits suggest that faster processes will dissipate more than slower processes. Here, we test the hypothesis that this relationship between speed and dissipation holds for stochastic paths far from equilibrium. To analyze stochastic paths on finite timescales, we derive an exact expression for the path probabilities of continuous-time Markov chains from the path summation solution to the master equation. We present a minimal model for a driven system in which relative energies of the initial and target states control the speed, and the nonequilibrium currents of a cycle control the dissipation. Although the hypothesis holds near equilibrium, we find that faster processes can dissipate less under far-from-equilibrium conditions because of strong currents. This model serves as a minimal prototype for designing kinetics to sculpt the nonequilibrium path space so that faster paths produce less dissipation.

Sample size dependence of tagged molecule dynamics in steady-state networks with bimolecular reactions: Cycle times of a light-driven pump

Here, steady-state reaction networks are inspected from the viewpoint of individual tagged molecules jumping among their chemical states upon the occurrence of reactive events. Such an agent-based viewpoint is useful for selectively characterizing the behavior of functional molecules, especially in the presence of bimolecular processes. We present the tools for simulating the jump dynamics both in the macroscopic limit and in the small-volume sample where the numbers of reactive molecules are of the order of few units with an inherently stochastic kinetics. The focus is on how an ideal spatial "compartmentalization" may affect the dynamical features of the tagged molecule. Our general approach is applied to a synthetic light-driven supramolecular pump composed of ring-like and axle-like molecules that dynamically assemble and disassemble, originating an average ring-through-axle directed motion under constant irradiation. In such an example, the dynamical feature of interest is the completion time of direct/inverse cycles of tagged rings and axles. We find a surprisingly strong robustness of the average cycle times with respect to the system's size. This is explained in the presence of rate-determining unimolecular processes, which may, therefore, play a crucial role in stabilizing the behavior of small chemical systems against strong fluctuations in the number of molecules.

Individual-Molecule Perspective Analysis of Chemical Reaction Networks: The Case of a Light-Driven Supramolecular Pump

The first study in which stochastic simulations of a two‐component molecular machine are performed in the mass‐action regime is presented. This system is an autonomous molecular pump consisting of a photoactive axle that creates a directed flow of rings through it by exploiting light energy away from equilibrium. The investigation demonstrates that the pump can operate in two regimes, both experimentally accessible, in which light‐driven steps can be rate‐determining or not. The number of photons exploited by an individual molecular pump, as well as the precision of cycling and the overall efficiency, critically rely on the operating regime of the machine. This approach provides useful information not only to guide the chemical design of a self‐assembling molecular device with desired features, but also to elucidate the effect of the environment on its performance, thus facilitating its experimental investigation.

Effects of temperature and mass conservation on the typical chemical sequences of hydrogen oxidation

Macroscopic properties of reacting mixtures are necessary to design synthetic strategies, determine yield, and improve the energy and atom efficiency of many chemical processes. The set of time-ordered sequences of chemical species are one representation of the evolution from reactants to products. However, only a fraction of the possible sequences is typical, having the majority of the joint probability and characterizing the succession of chemical nonequilibrium states. Here, we extend a variational measure of typicality and apply it to atomistic simulations of a model for hydrogen oxidation over a range of temperatures. We demonstrate an information-theoretic methodology to identify typical sequences under the constraints of mass conservation. Including these constraints leads to an improved ability to learn the chemical sequence mechanism from experimentally accessible data. From these typical sequences, we show that two quantities defining the variational typical set of sequences-the joint entropy rate and the topological entropy rate-increase linearly with temperature. These results suggest that, away from explosion limits, data over a narrow range of thermodynamic parameters could be sufficient to extrapolate these typical features of combustion chemistry to other conditions.