Learning the mechanisms of chemical disequilibria

Molecular Dynamics Simulation of Thermal Nonequilibrium and Chemical Reaction Processes in Hydrogen Combustion

Using reactive force field (ReaxFF) and molecular dynamics simulation, we investigate the combustion process of hydrogen-oxygen systems in initial thermal nonequilibrium states with different translational and rovibrational temperatures for oxygen. The system studied in this work contains 300 oxygen molecules and 700 hydrogen molecules with a density of 7 times the air density. For this system, the characteristic relaxation times of oxygen and hydrogen vibrational energies are 0.173 and 0.249 ns, respectively. 0.6% of hydrogen undergoes a chemical reaction with oxygen during the thermal nonequilibrium relaxation stage. For the distribution of translational energy and vibrational energy of oxygen in the thermal nonequilibrium state, the maximum mean error of the statistical distribution in the simulation and the Boltzmann distribution at temperature calculated from the average kinetic energy of molecules is about 2.25 × 10-5. At the same time, it was observed in the simulation that many-body interactions play a certain role in the combustion process. Furthermore, we compare the ignition time and temperature rise behavior of different combustion mechanisms and molecular dynamics simulations starting from the thermal equilibrium state. These results will provide meaningful references for the construction of thermal nonequilibrium combustion chemical reaction mechanisms.

The Typical Set and Entropy in Stochastic Systems with Arbitrary Phase Space Growth

The existence of the typical set is key for data compression strategies and for the emergence of robust statistical observables in macroscopic physical systems. Standard approaches derive its existence from a restricted set of dynamical constraints. However, given its central role underlying the emergence of stable, almost deterministic statistical patterns, a question arises whether typical sets exist in much more general scenarios. We demonstrate here that the typical set can be defined and characterized from general forms of entropy for a much wider class of stochastic processes than was previously thought. This includes processes showing arbitrary path dependence, long range correlations or dynamic sampling spaces, suggesting that typicality is a generic property of stochastic processes, regardless of their complexity. We argue that the potential emergence of robust properties in complex stochastic systems provided by the existence of typical sets has special relevance to biological systems.

Optimizing dynamical functions for speed with stochastic paths

Living systems are built from microscopic components that function dynamically; they generate work with molecular motors, assemble and disassemble structures such as microtubules, keep time with circadian clocks, and catalyze the replication of DNA. How do we implement these functions in synthetic nanostructured materials to execute them before the onset of dissipative losses? Answering this question requires a quantitative understanding of when we can improve performance and speed while minimizing the dissipative losses associated with operating in a fluctuating environment. Here, we show that there are four modalities for optimizing dynamical functions that can guide the design of nanoscale systems. We analyze Markov models that span the design space: a clock, ratchet, replicator, and self-assembling system. Using stochastic thermodynamics and an exact expression for path probabilities, we classify these models of dynamical functions based on the correlation of speed with dissipation and with the chosen performance metric. We also analyze random networks to identify the model features that affect their classification and the optimization of their functionality. Overall, our results show that the possible nonequilibrium paths can determine our ability to optimize the performance of dynamical functions, despite ever-present dissipation, when there is a need for speed.

Dynamic scaling of stochastic thermodynamic observables for chemical reactions at and away from equilibrium

Physical kinetic roughening processes are well-known to exhibit universal scaling of observables that fluctuate in space and time. Are there analogous dynamic scaling laws that are unique to the chemical reaction mechanisms available synthetically and occurring naturally? Here, we formulate an approach to the dynamic scaling of stochastic fluctuations in thermodynamic observables at and away from equilibrium. Both analytical expressions and numerical simulations confirm our dynamic scaling ansatz with associated scaling exponents, function, and law. A survey of common chemical mechanisms reveals classes that organize according to the molecularity of the reactions involved, the nature of the reaction vessel and external reservoirs, (non)equilibrium conditions, and the extent of autocatalysis in the reaction network. Varying experimental parameters, such as temperature, can cause coupled reactions capable of chemical feedback to transition between these classes. While path observables, such as the dynamical activity, have scaling exponents that are time-independent, the variance in the entropy production and flow can have time-dependent scaling exponents and self-averaging properties as a result of temporal correlations that emerge during thermodynamically irreversible processes. Altogether, these results establish dynamic universality classes in the nonequilibrium fluctuations of thermodynamic observables for well-mixed chemical reactions.

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.

Typical Stochastic Paths in the Transient Assembly of Fibrous Materials

When chemically fueled, molecular self-assembly can sustain dynamic aggregates of polymeric fibers-hydrogels-with tunable properties. If the fuel supply is finite, the hydrogel is transient, as competing reactions switch molecular subunits between active and inactive states, drive fiber growth and collapse, and dissipate energy. Because the process is away from equilibrium, the structure and mechanical properties can reflect the history of preparation. As a result, the formation of these active materials is not readily susceptible to a statistical treatment in which the configuration and properties of the molecular building blocks specify the resulting material structure. Here, we illustrate a stochastic-thermodynamic and information-theoretic framework for this purpose and apply it to these self-annihilating materials. Among the possible paths, the framework variationally identifies those that are typical-loosely, the minimum number with the majority of the probability. We derive these paths from computer simulations of experimentally-informed stochastic chemical kinetics and a physical kinetics model for the growth of an active hydrogel. The model reproduces features observed by confocal microscopy, including the fiber length, lifetime, and abundance as well as the observation of fast fiber growth and stochastic fiber collapse. The typical mesoscopic paths we extract are less than 0.23% of those possible, but they accurately reproduce material properties such as mean fiber length.

Maximum Entropy Principle in Statistical Inference: Case for Non-Shannonian Entropies

In this Letter, we show that the Shore-Johnson axioms for the maximum entropy principle in statistical estimation theory account for a considerably wider class of entropic functional than previously thought. Apart from a formal side of the proof where a one-parameter class of admissible entropies is identified, we substantiate our point by analyzing the effect of weak correlations and by discussing two pertinent examples: two-qubit quantum system and transverse-momentum behavior of hadrons in high-energy proton-proton collisions.

Entrance and escape dynamics for the typical set

According to the asymptotic equipartition property, sufficiently long sequences of random variables converge to a set that is typical. While the size and probability of this set are central to information theory and statistical mechanics, they can often only be estimated accurately in the asymptotic limit due to the exponential growth in possible sequences. Here we derive a time-inhomogeneous dynamics that constructs the properties of the typical set for all finite length sequences of independent and identically distributed random variables. These dynamics link the finite properties of the typical set to asymptotic results and allow the typical set to be applied to small and transient systems. The main result is a geometric mapping-the triangle map-relating sequences of random variables of length n to those of length n+1. We show that the number of points in this map needed to quantify the properties of the typical set grows linearly with sequence length, despite the exponential growth in the number of typical sequences. We illustrate the framework for the Bernoulli process and the Schlögl model for autocatalytic chemical reactions and demonstrate both the convergence to asymptotic limits and the ability to reproduce exact calculations.

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.

Nonequilibrium phase coexistence and criticality near the second explosion limit of hydrogen combustion

While hydrogen is a promising source of clean energy, the safety and optimization of hydrogen technologies rely on controlling ignition through explosion limits: pressure-temperature boundaries separating explosive behavior from comparatively slow burning. Here, we show that the emergent nonequilibrium chemistry of combustible mixtures can exhibit the quantitative features of a phase transition. With stochastic simulations of the chemical kinetics for a model mechanism of hydrogen combustion, we show that the boundaries marking explosive domains of kinetic behavior are nonequilibrium critical points. Near the pressure of the second explosion limit, these critical points terminate the transient coexistence of dynamical phases-one that autoignites and another that progresses slowly. Below the critical point temperature, the chemistry of these phases is indistinguishable. In the large system limit, the pseudo-critical temperature converges to the temperature of the second explosion limit derived from mass-action kinetics.

On the role of the termolecular reactions 2O2 + H2 → 2HO2 and 2O2 + H2 → H + HO2 + O2 in formation of the first radicals in hydrogen combustion: ab initio predictions of energy barriers

We have investigated the role of termolecular reactions in the early chemistry of hydrogen combustion. We performed molecular chemical dynamics simulations using ReaxFF in LAMMPS to identify potential initial reactions for a 1 : 4 mixture of H₂ : O₂ in the NVT ensemble at density 276.3 kg m^-3 and ∼3000 K (∼4000 atm) and ∼4000 K (∼5000 atm), and then characterized the saddle points for those reactions using ab initio methods: CCSD(T) = FC/cc-pVTZ//MP2/6-31G, CCSD(T) = FULL/aug-cc-pVTZ//CCSD = FC/cc-pVTZ and CASSCF MP2/6-31G//MP2/6-31G. The main initial reaction is H₂ + O₂ → H + HO₂, frequently occurring in the presence of a second O₂ as a third body; that is, 2O₂ + H₂ → H + HO₂ + O₂. The second most frequent reaction is 2O₂ + H₂ → 2HO₂. We found three saddle points on the triplet PES of these termolecular reactions: one for 2O₂ + H₂ → H + HO₂ + O₂ and two for 2O₂ + H₂ → 2HO₂. In the latter case, one has a symmetric structure consistent with simultaneous formation of two HO₂ and the other corresponds to a bimolecular reaction between O₂ and H₂ that is "interrupted" by a second O₂ before going to completion. The classical barrier height of the symmetric saddle point for 2O₂ + H₂ → 2HO₂ is 49.8 kcal mol^-1. The barrier to H₂ + O₂ → H + HO₂ is 58.9 kcal mol^-1. The termolecular reaction will be competitive with H₂ + O₂ → H + HO₂ only at sufficiently high pressures.