Many-particle effects on the relaxation kinetics of fast reversible reactions of the type A+B⇄C

A new method for choosing the computational cell in stochastic reaction-diffusion systems

How to choose the computational compartment or cell size for the stochastic simulation of a reaction-diffusion system is still an open problem, and a number of criteria have been suggested. A generalized measure of the noise for finite-dimensional systems based on the largest eigenvalue of the covariance matrix of the number of molecules of all species has been suggested as a measure of the overall fluctuations in a multivariate system, and we apply it here to a discretized reaction-diffusion system. We show that for a broad class of first-order reaction networks this measure converges to the square root of the reciprocal of the smallest mean species number in a compartment at the steady state. We show that a suitably re-normalized measure stabilizes as the volume of a cell approaches zero, which leads to a criterion for the maximum volume of the compartments in a computational grid. We then derive a new criterion based on the sensitivity of the entire network, not just of the fastest step, that predicts a grid size that assures that the concentrations of all species converge to a spatially-uniform solution. This criterion applies for all orders of reactions and for reaction rate functions derived from singular perturbation or other reduction methods, and encompasses both diffusing and non-diffusing species. We show that this predicts the maximal allowable volume found in a linear problem, and we illustrate our results with an example motivated by anterior-posterior pattern formation in Drosophila, and with several other examples.

Manifestation of macroscopic correlations in elementary reaction kinetics. I. Irreversible reaction A+A-->product

Using an modern many-particle method for the derivation of non-Markovian binary kinetic equations, we have treated theoretically the applicability of the encounter theory (ET) (the prototype of the collision theory) concepts to the widely known diffusion assisted irreversible bulk reaction A+A-->product (for example, radical reaction) in dilute solutions. The method shows that the agreement with the ET is observed when the familiar integral ET is employed which in this method is just a step in the derivation of kinetic equations. It allows for two-particle correlations only, but fails to take account of correlation of reactant simultaneously with the partner of the encounter and the reactant in the bulk. However, the next step leading to the modified ET under transformation of equations to the regular form both extends the time range of the applicability of ET rate equation (as it was for reactions proceeding with one of the reactants in excess), and gives the equation of the generalized ET. In full agreement with physical considerations, this theory reveals macroscopic correlations induced by the encounters in the reservoir of free walks. This means that the encounters of reactants in solution are correlated on a rather large time interval of the reaction. Though any nonstationary (non-Markovian) effects manifest themselves rather weakly in the kinetics of the bimolecular reaction in question, just the existence of the revealed macroscopic correlations in the binary theory is of primary importance. In particular, it means that the well-known phenomena which are generally considered to be associated solely with correlation of particles on the encounter (for example, chemically induced dynamic nuclear polarization) may be induced by correlation in the reservoir of free random walks of radicals in solution.

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.

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

For diffusion-limited reversible A+A<==>B reactions we reexamine two fundamental concepts of classical chemical kinetics-the notion of "chemical equilibrium" and the "law of mass action." We consider a general model with distance-dependent reaction rates, such that any pair of A particles, performing standard random walks on sites of a d-dimensional lattice and being at a distance mu apart of each other at time moment t, may associate forming a B particle at the rate k+(mu). In turn, any randomly moving B particle may spontaneously dissociate at the rate k-(lambda) into a geminate pair of As "born" at a distance lambda apart of each other. Within a formally exact approach based on Gardiner's Poisson representation method we show that the asymptotic t=infinity state attained by such diffusion-limited reactions is generally not a true thermodynamic equilibrium, but rather a nonequilibrium steady state, and that the law of mass action is invalid. The classical concepts hold only in case when the ratio k+(mu)k-(mu) does not depend on mu for any mu.

Elementary Steps in Excited-State Proton Transfer

The absorption of a photon by a hydroxy-aromatic photoacid triggers a cascade of events contributing to the overall phenomenon of intermolecular excited-state proton transfer. The fundamental steps involved were studied over the last 20 years using a combination of theoretical and experimental techniques. They are surveyed in this sequel in sequential order, from fast to slow. The excitation triggers an intramolecular charge transfer to the ring system, which is more prominent for the anionic base than the acid. The charge redistribution, in turn, triggers changes in hydrogen-bond strengths that set the stage for the proton-transfer step itself. This step is strongly influenced by the solvent, resulting in unusual dependence of the dissociation rate coefficient on water content, temperature, and isotopic substitution. The photolyzed proton can diffuse in the aqueous solution in a mechanism that involves collective changes in hydrogen-bonding. On longer times, it may recombine adiabatically with the excited base or quench it. The theory for these diffusion-influenced geminate reactions has been developed, showing nice agreement with experiment. Finally, the effect of inert salts, bases, and acids on these reactions is analyzed.

Mobility of taxol in microtubule bundles

Mobility of taxol inside microtubules was investigated using fluorescence recovery after photobleaching on flow-aligned bundles. Bundles were made of microtubules with either GMPCPP or GTP at the exchangeable site on the tubulin dimer. Recovery times were sensitive to bundle thickness and packing, indicating that taxol molecules are able to move laterally through the bundle. The density of open binding sites along a microtubule was varied by controlling the concentration of taxol in solution for GMPCPP samples. With >63% sites occupied, recovery times were independent of taxol concentration and, therefore, inversely proportional to the microscopic dissociation rate, k(off). It was found that 10k(off)(GMPCPP) approximately equal k(off)(GTP), consistent with, but not fully accounting for, the difference in equilibrium constants for taxol on GMPCPP and GTP microtubules. With <63% sites occupied, recovery times decreased as approximately [Tax](-1/5) for both types of microtubules. We conclude that the diffusion of taxol inside the microtubule bundle is hindered by rebinding events when open sites are within approximately 7 nm of each other.

Rigorous derivation of the long-time asymptotics for reversible binding

Using an iterative solution in Laplace-Fourier space, we supply a rigorous mathematical proof for the long-time asymptotics of reversible binding in one dimension. The asymptotic power law and its concentration dependent prefactor result from diffusional and many-body effects which, unlike for the corresponding irreversible reaction and in classical chemical kinetics, play a dominant role in shaping the approach to equilibrium.