An official website of the United States government
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.
more »
« less
Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures
It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.
more »
« less
- PAR ID:
- 10582894
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 15
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Low temperature plasmas are open driven thermodynamic systems capable of increasing the free energy of the mass that flows through them. An interesting thing about low temperature plasmas is that different species have different temperatures at the same location in space. Since thermal equilibrium cannot be assumed, many of the familiar results of equilibrium thermodynamics cannot be applied in their familiar form to predict, e.g., the direction of a chemical reaction. From the perspective of classical processing governed by thermal equilibrium, examples of highly unexpected gas-phase chemical reactions (CO2 dissociation, NO, N2H4, O3 synthesis) and solid material transformations (surface activation, size-focusing, and hyperdoping) promoted by low temperature plasmas are presented. The lack of a known chemical reaction equilibrium criterion prevents assessment of predictive kinetics models of low temperature plasmas, to ensure that they comply with the laws of thermodynamics. There is a need for a general method to predict chemical reaction equilibrium in low temperature plasmas or an alternative method to establish the thermodynamic admissibility of a proposed kinetics model. Toward those ends, two ideas are explored in this work. The first idea is that chemical reactions in low temperature plasmas proceed toward a thermal equilibrium state at an effective temperature intermediate between the neutral gas temperature and the electron temperature. The effective temperature hypothesis is simple, and surprisingly is adequate for elucidation in some systems, but it lacks generality. The general equation for nonequilibrium reversible–irreversible coupling (GENERIC) is a general beyond equilibrium thermodynamics framework that can be used to rigorously establish the thermodynamic admissibility of a set of dynamic modeling equations, such as a kinetic model, without knowledge of the final state that the system is tending toward. The use of GENERIC is described by way of example using a two-temperature hydrodynamic model from the literature. The conclusion of the GENERIC analysis presented in this work is that the concept of superlocal equilibrium is thermodynamically admissible and may be applied to describe low temperature plasmas, provided that appropriate terms are included for exchange of internal energy and momentum between different species that may have different temperatures and bulk velocities at the same location in space. The concept of superlocal equilibrium is expected to be of utility in future work focused on deriving equilibrium criteria for low temperature plasmas.more » « less
-
Hydrogels cross-linked by dynamic covalent chemistry (DCC) are stiff and remodelable, making them ideal biomimetics for tissue engineering applications. Due to the reversibility of DCC cross-links, the opportunity exists to transiently control hydrogel network formation through the use of small molecule competitors. Specifically, we incorporate low molecular weight competitors that reversibly disrupt the formation of hydrazone cross-links as they diffuse through a recombinant hydrogel. Using complementary experimental, computational, and theoretical polymer physics approaches, we present a family of competitors that predictably alter hydrogel gelation time and mechanics. By changing the competitor chemistry, we connect key reaction parameters (forward and reverse reactions rates and thermodynamic equilibrium constants) to the delayed onset of a percolated network, increased hydrogel gelation time, and transient control of hydrogel stiffness. Using human intestinal organoids as a model system, we demonstrate the ability to tune gelation kinetics of a recombinant hydrogel for uniform encapsulation of individual, patient-derived stem cells and their proliferation into three-dimensional structures. Taken together, our data establish a validated framework to relate molecular-level parameters of transient competitors to predicted macromolecular-network properties. As interest in biomimetic, DCC-cross-linked hydrogels continues to grow, these results will enable the rationale design of bespoke, dynamic biomaterials for tissue engineering.more » « less
-
Many complex fluids can be described by continuum hydrodynamic field equations, to which noise must be added in order to capture thermal fluctuations. In almost all cases, the resulting coarse-grained stochastic partial differential equations carry a short-scale cutoff, which is also reflected in numerical discretisation schemes. We draw together our recent findings concerning the construction of such schemes and the interpretation of their continuum limits, focusing, for simplicity, on models with a purely diffusive scalar field, such as ‘Model B’ which describes phase separation in binary fluid mixtures. We address the requirement that the steady-state entropy production rate (EPR) must vanish for any stochastic hydrodynamic model in a thermal equilibrium. Only if this is achieved can the given discretisation scheme be relied upon to correctly calculate the nonvanishing EPR for ‘active field theories’ in which new terms are deliberately added to the fluctuating hydrodynamic equations that break detailed balance. To compute the correct probabilities of forward and time-reversed paths (whose ratio determines the EPR), we must make a careful treatment of so-called ‘spurious drift’ and other closely related terms that depend on the discretisation scheme. We show that such subtleties can arise not only in the temporal discretisation (as is well documented for stochastic ODEs with multiplicative noise) but also from spatial discretisation, even when noise is additive, as most active field theories assume. We then review how such noise can become multiplicative via off-diagonal couplings to additional fields that thermodynamically encode the underlying chemical processes responsible for activity. In this case, the spurious drift terms need careful accounting, not just to evaluate correctly the EPR but also to numerically implement the Langevin dynamics itself.more » « less
-
We develop a convergent reaction-drift-diffusion master equation (CRDDME) to facilitate the study of reaction processes in which spatial transport is influenced by drift due to one-body potential fields within general domain geometries. The generalized CRDDME is obtained through two steps. We first derive an unstructured grid jump process approximation for reversible diffusions, enabling the simulation of drift-diffusion processes where the drift arises due to a conservative field that biases particle motion. Leveraging the Edge-Averaged Finite Element method, our approach preserves detailed balance of drift-diffusion fluxes at equilibrium, and preserves an equilibrium Gibbs-Boltzmann distribution for particles undergoing drift-diffusion on the unstructured mesh. We next formulate a spatially-continuous volume reactivity particle-based reaction-drift-diffusion model for reversible reactions of the form A + B <–> C. A finite volume discretization is used to generate jump process approximations to reaction terms in this model. The discretization is developed to ensure the combined reaction-drift-diffusion jump process approximation is consistent with detailed balance of reaction fluxes holding at equilibrium, along with supporting a discrete version of the continuous equilibrium state. The new CRDDME model represents a continuous-time discrete-space jump process approximation to the underlying volume reactivity model. We demonstrate the convergence and accuracy of the new CRDDME through a number of numerical examples, and illustrate its use on an idealized model for membrane protein receptor dynamics in T cell signaling.more » « less
