Corrections to the law of mass action and properties of the asymptotic t=∞ state for reversible diffusion-limited reactions

Extinction times in autocatalytic systems

Systems evolving under the influence of autocatalytic processes have been the subject of much study due to their appearance in a wide variety of contexts. Well known examples range from the simplest autocatalytic chemical reactions, biochemical models of ribozyme or prion replication, right up to the duplication of entire organisms in models of population dynamics. While the deterministic approach frequently taken to model such systems is adequate in the limit of large reactant numbers, the intrinsic fluctuations exert an important influence on the dynamics when the supply of reactants is limited. In particular, when combined with spontaneous degradation of the autocatalytic reactant population, such fluctuations can lead to the extinction of that population. In this paper, we study reversible autocatalytic processes of the form X + Y ⇌ 2X in the limit that there exists a surplus of Y and in the presence of a spontaneous degradation process X → Z. Through the use of the Poisson representation, we identify an exact analytical expression for the mean extinction time of the X population. We show that the exact result can be neatly approximated by an Arrhenius-like expression involving an effective activation energy separating a quasi-stationary state from the extinct state.

Distance-dependent diffusion-controlled reaction of •NO and O2•- at chemical equilibrium with ONOO-

The fast reaction of (•)NO and O(2)(•-) to give ONOO(-) has been extensively studied at irreversible conditions, but the reasons for the wide variations in observed forward rate constants (3.8 ≤ k(f) ≤ 20 × 10(9) M(-1) s(-1)) remain unexplained. We characterized the diffusion-dependent aqueous (pH > 12) chemical equilibrium of the form (•)NO + O(2)(•-) = ONOO(-) with respect to its dependence on temperature, viscosity, and [ONOO(-)](eq) by determining [ONOO(-)](eq) and [(•)NO](eq). The equilibrium forward reaction rate constant (k(f)(eq)) has negative activation energy, in contrast to that found under irreversible conditions. In contradiction to the law of mass action, we demonstrate that the equilibrium constant depends on ONOO(-) concentration. Therefore, a wide range of k(f)(eq) values could be derived (7.5-21 × 10(9) M(-1) s(-1)). Of general interest, the variations in k(f) can thus be explained by its dependence on the distance between ONOO(-) particles (sites of generation of (•)NO and O(2)(•-)).

The intranuclear environment

Many of the chapters in this volume are concerned with processes or structures inside the nucleus, and it is relevant to consider the properties of their environment, or rather of the multiple different and specific environments that must exist in local regions of the highly heterogeneous intranuclear space. Relatively little is known about the fundamental physical properties of these environments, and theoretical treatments of phenomena in such concentrated mixtures of charged macromolecules are complex and as yet poorly developed. Some of the phenomena that occur at the molecular level are unexpected and counterintuitive for biologists, although well known to colloid and polymer scientists; for example, the existence of short-range attractive forces between macromolecules or structures with like charges. As a background for the chapters that follow, we consider here some of the particular features of intranuclear environments, how they may influence processes and structures in the nucleus, and their implications for working with nuclei.

Aging processes in reversible reaction-diffusion systems

Reversible reaction-diffusion systems display anomalous dynamics characterized by a power-law relaxation toward stationarity. In this paper we study in the aging regime the nonequilibrium dynamical properties of some model systems with reversible reactions. Starting from the exact Langevin equations describing these models, we derive expressions for two-time correlation and autoresponse functions and obtain a simple aging behavior for these quantities. The autoresponse function is thereby found to depend on the specific nature of the chosen perturbation of the system.

Rate kernel theory for pseudo-first-order kinetics of diffusion-influenced reactions and application to fluorescence quenching kinetics

Theoretical foundation of rate kernel equation approaches for diffusion-influenced chemical reactions is presented and applied to explain the kinetics of fluorescence quenching reactions. A many-body master equation is constructed by introducing stochastic terms, which characterize the rates of chemical reactions, into the many-body Smoluchowski equation. A Langevin-type of memory equation for the density fields of reactants evolving under the influence of time-independent perturbation is derived. This equation should be useful in predicting the time evolution of reactant concentrations approaching the steady state attained by the perturbation as well as the steady-state concentrations. The dynamics of fluctuation occurring in equilibrium state can be predicted by the memory equation by turning the perturbation off and consequently may be useful in obtaining the linear response to a time-dependent perturbation. It is found that unimolecular decay processes including the time-independent perturbation can be incorporated into bimolecular reaction kinetics as a Laplace transform variable. As a result, a theory for bimolecular reactions along with the unimolecular process turned off is sufficient to predict overall reaction kinetics including the effects of unimolecular reactions and perturbation. As the present formulation is applied to steady-state kinetics of fluorescence quenching reactions, the exact relation between fluorophore concentrations and the intensity of excitation light is derived.

Molecular-dynamics simulations for nonclassical kinetics of diffusion-controlled bimolecular reactions

Molecular-dynamics simulations are presented for the diffusion-controlled bimolecular reaction A+B<==>C in two and three dimensions. The reactants and solvent molecules are modeled as spheres interacting via continuous potential-energy functions. The interaction potential between two reactants contains a deep well that results in a reaction. When the solvent concentration is low and the reactant dynamics is essentially ballistic, the system reaches equilibrium rapidly, and the reaction follows classical kinetics with exponential decay to the equilibrium. When the solvent concentration is high the particles enter the normal diffusion regime quickly and nonclassical behavior is observed, i.e., the reactant concentrations approach equilibrium as t(-d/2) where d is the dimensionality of space. When the reaction well depth is large, however, the reaction becomes irreversible within the simulation time. In this case the reactant concentrations decay as t(-d/4). Interestingly this behavior is also observed at intermediate times for reversible reactions.