Reactive symbol sequences for a model of hydrogen combustion

Transient, macroscopic states of chemical disequilibrium are born out of the microscopic dynamics of molecules. As a reaction mixture evolves, the temporal patterns of chemical species encodes some of this dynamical information, while their statistics are a manifestation of the bulk kinetics. Here, we define a chemically-informed symbolic dynamics as a coarse-grained representation of classical molecular dynamics, and analyze the sequences of chemical species for a model of hydrogen combustion. We use reactive molecular dynamics simulations to generate the sequences and derive probability distributions for sequence observables: the reaction time scales and the chain length - the total number of reactions between initiation of a reactant and termination at products. The time scales and likelihood of the sequences depend strongly on the chain length, temperature, and density. Temperature suppresses the uncertainty in chain length for hydrogen sequences, but enhances the uncertainty in oxygen sequence chain lengths. This method of analyzing a surrogate chemical symbolic dynamics reduces the complexity of the chemistry from the atomistic to the molecular level and has the potential for extension to more complicated reaction systems.