Observation of laser speckle effects and nonclassical kinetics in an elementary chemical reaction

Upscaling Mixing in Highly Heterogeneous Porous Media via a Spatial Markov Model

In this work, we develop a novel Lagrangian model able to predict solute mixing in heterogeneous porous media. The Spatial Markov model has previously been used to predict effective mean conservative transport in flows through heterogeneous porous media. In predicting effective measures of mixing on larger scales, knowledge of only the mean transport is insufficient. Mixing is a small scale process driven by diffusion and the deformation of a plume by a non-uniform flow. In order to capture these small scale processes that are associated with mixing, the upscaled Spatial Markov model must be extended in such a way that it can adequately represent fluctuations in concentration. To address this problem, we develop downscaling procedures within the upscaled model to predict measures of mixing and dilution of a solute moving through an idealized heterogeneous porous medium. The upscaled model results are compared to measurements from a fully resolved simulation and found to be in good agreement.

Incomplete mixing and reactions in laminar shear flow

Incomplete mixing of reactive solutes is well known to slow down reaction rates relative to what would be expected from assuming perfect mixing. In purely diffusive systems, for example, it is known that small initial fluctuations in reactant concentrations can lead to reactant segregation, which in the long run can reduce global reaction rates due to poor mixing. In contrast, nonuniform flows can enhance mixing between interacting solutes. Thus, a natural question arises: Can nonuniform flows sufficiently enhance mixing to restrain incomplete mixing effects and, if so, under what conditions? We address this question by considering a specific and simple case, namely, a laminar pure shear reactive flow. Two solution approaches are developed: a Lagrangian random walk method and a semianalytical solution. The results consistently highlight that if shear effects in the system are not sufficiently strong, incomplete mixing effects initially similar to purely diffusive systems will occur, slowing down the overall reaction rate. Then, at some later time, dependent on the strength of the shear, the system will return to behaving as if it were well mixed, but represented by a reduced effective reaction rate.

Predicting the enhancement of mixing-driven reactions in nonuniform flows using measures of flow topology

The ability for reactive constituents to mix is often the key limiting factor for the completion of reactions across a huge range of scales in a variety of media. In flowing systems, deformation and shear enhance mixing by bringing constituents into closer proximity, thus increasing reaction potential. Accurately quantifying this enhanced mixing is key to predicting reactions and typically is done by observing or simulating scalar transport. To eliminate this computationally expensive step, we use a Lagrangian stochastic framework to derive the enhancement to reaction potential by calculating the collocation probability of particle pairs in a heterogeneous flow field accounting for deformations. We relate the enhanced reaction potential to three well known flow topology metrics and demonstrate that it is best correlated to (and asymptotically linear with) one: the largest eigenvalue of the (right) Cauchy-Green tensor.

Geometry-controlled kinetics

It has long been appreciated that the transport properties of molecules can control reaction kinetics. This effect can be characterized by the time it takes a diffusing molecule to reach a target-the first-passage time (FPT). Determining the FPT distribution in realistic confined geometries has until now, however, seemed intractable. Here, we calculate this FPT distribution analytically and show that transport processes as varied as regular diffusion, anomalous diffusion, and diffusion in disordered media and fractals, fall into the same universality classes. Beyond the theoretical aspect, this result changes our views on standard reaction kinetics and we introduce the concept of 'geometry-controlled kinetics'. More precisely, we argue that geometry-and in particular the initial distance between reactants in 'compact' systems-can become a key parameter. These findings could help explain the crucial role that the spatial organization of genes has in transcription kinetics, and more generally the impact of geometry on diffusion-limited reactions.

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.

Nonclassical kinetics of an elementary A+B-->C reaction-diffusion system showing effects of a speckled initial reactant distribution and eventual self-segregation: experiments

We demonstrate here the implementation of an experimental system suitable for the study of the diffusion limited A+B-->0, nonclassical reaction behavior. Using a combination of a fluorescent calcium indicator and a calcium ion which is initially "caged," a pulse of near-UV light initiates the reaction which is followed as product formation vs time. Sensitive dependence on the initial reactant distribution is observed through patterns in the uncaging UV light profile. In one case, the reaction progress passes through two nonclassical time regimes, one due to roughness originating from laser speckles, followed by one consistent with the three-dimensional Zeldovich rate of (1/rho(A)-1/rho(A0)) approximately t(3/4), with features matching Monte Carlo simulations on this initial distribution. This behavior is contrasted with reactions initiated by a homogeneous source which induces random initial reactant distributions, though both systems seem to approach the asymptotic limit of self-segregation of reactants.

Dynamics of the depletion zone at a finite-sized imperfect trap in two dimensions: photobleaching experiments and simulations

The kinetics of the growth of depletion zones around a static trap in an effective two-dimensional geometry were studied experimentally with photobleaching of fluorescein dye by a focused laser beam. The phototrap served as an imperfect trap with a finite size. The growth of the depletion zone was monitored by the theta distance, defined as the distance from the trap to the point where the concentration of the reactants reaches a given arbitrary fraction theta (0<theta<1) of its initial value, which could be directly measured experimentally. At the asymptotic limit, the results confirm the theoretical nonuniversal t(theta/2) scaling behavior for the theta distance. We also find an effect of fast expansion at an early time of the depletion zone inside an imperfect trap. Both the imperfect trapping strength and the finite trap size are found to control the early-time behavior, while the trap shape does not much affect the dynamics of the theta distance. A dimensional crossover was found for a perfect trap with a finite radius, when the theta distance was measured from the trap surface. The actual trapping efficiency was determined for different laser powers of the phototrap. Results are supported by analytical equations, exact enumerations, and Monte Carlo simulations.

Reaction-diffusion dynamics: confrontation between theory and experiment in a microfluidic reactor

We confront, quantitatively, the theoretical description of the reaction-diffusion process of a second-order reaction to experiment. The reaction at work is Ca(2+)/CaGreen, a fluorescent tracer for calcium. The reactor is a T-shaped microchannel, 10 microm deep, 200 microm wide, and 2 cm long. The experimental measurements are compared with the two-dimensional numerical simulation of the reaction-diffusion equations. We find good agreement between theory and experiment. From this study, one may propose a method of measurement of various quantities, such as the kinetic rate of the reaction, in conditions yet inaccessible to conventional methods.

Reaction front structure in the diffusion-limited A+B model with initially randomized reactants

Subtle features of the reaction front formation in the A+B-->0 reaction are reported for the initially random and equal A+B reactant distribution. Three nonclassical parameters (initial linewidth, minimum, and maximum), for each interparticle gap and nearest neighbor distance distributions, are derived, as a function of time, using Monte Carlo simulations. These empirical front measures and their temporal scaling exponents are compared with the previously studied ones for the reactant interparticle distributions.