An Efficient Brownian Dynamics Method for Evaluating Inertial Dynamic Effects on Diffusion-Influenced Reactions

Accurate analytical calculation of the rate coefficient for the diffusion-controlled reactions due to hyperbolic diffusion

Using an approach based on the diffusion analog of the Cattaneo-Vernotte differential model, we find the exact analytical solution to the corresponding time-dependent linear hyperbolic initial boundary value problem, describing irreversible diffusion-controlled reactions under Smoluchowski's boundary condition on a spherical sink. By means of this solution, we extend exact analytical calculations for the time-dependent classical Smoluchowski rate coefficient to the case that includes the so-called inertial effects, occurring in the host media with finite relaxation times. We also present a brief survey of Smoluchowski's theory and its various subsequent refinements, including works devoted to the description of the short-time behavior of Brownian particles. In this paper, we managed to show that a known Rice's formula, commonly recognized earlier as an exact reaction rate coefficient for the case of hyperbolic diffusion, turned out to be only its approximation being a uniform upper bound of the exact value. Here, the obtained formula seems to be of great significance for bridging a known gap between an analytically estimated rate coefficient on the one hand and molecular dynamics simulations together with experimentally observed results for the short times regime on the other hand. A particular emphasis has been placed on the rigorous mathematical treatment and important properties of the relevant initial boundary value problems in parabolic and hyperbolic diffusion theories.

Atomistic description of molecular binding processes based on returning probability theory

The efficiency of molecular binding such as host-guest binding is commonly evaluated in terms of kinetics, such as rate coefficients. In general, to compute the coefficient of the overall binding process, we need to consider both the diffusion of reactants and barrier crossing to reach the bound state. Here, we develop a methodology of quantifying the rate coefficient of binding based on molecular dynamics simulation and returning probability (RP) theory proposed by Kim and Lee [J. Chem. Phys. 131, 014503 (2009)]. RP theory provides a tractable formula of the rate coefficient in terms of the thermodynamic stability and kinetics of the intermediate state on a predefined reaction coordinate. In this study, the interaction energy between reactants is utilized as the reaction coordinate, enabling us to effectively describe the reactants' relative position and orientation on one-dimensional space. Application of this method to the host-guest binding systems, which consist of β-cyclodextrin and small guest molecules, yields the rate coefficients consistent with the experimental results.

Dynamics theory for molecular liquids based on an interaction site model

Dynamics theories for molecular liquids based on an interaction site model have been developed over the past few decades and proved to be powerful tools to investigate various dynamical phenomena.

Irreversible bimolecular reactions with inertia: from the trapping to the target setting at finite densities

We investigate numerically pseudo-first-order irreversible bimolecular reactions of the type A + B → B between hard spheres undergoing event-driven Brownian dynamics. We study the encounter rate and the survival probability of A particles as functions of the packing fraction ϕ in the trapping (a single particle diffusing among static non-overlapping traps) and target (many traps diffusing in the presence of a single static target particle) settings, as well as in the case of diffusing traps and particles (full mobility). We show that, since inertial effects are accounted for in our simulation protocol, the standard Smoluchowski theory of coagulation of non-interacting colloids is recovered only at times greater than a characteristic time Δt, marking the transition from the under-damped to the over-damped regime. We show that the survival probability S(t) decays exponentially during this first stage, with a rate 1/τ0 is proportional to φ. Furthermore, we work out a simple analytical expression that is able to capture to an excellent extent the numerical results for t < Δt at low and intermediate densities. Moreover, we demonstrate that the time constant of the asymptotic exponential decay of S(t) for diffusing traps and particles is k(S)(-1), where kS = 4π(DA + DB)Rρ is the Smoluchowski rate. Detailed analyses of the effective decay exponent β = d [log(-logS(t))]/d (logt) and of the steady-state encounter rate reveal that the full mobility and trapping problem are characterized by very similar kinetics, rather different from the target problem. Our results do not allow one to ascertain whether the prediction S(t) is proportional to exp(-at(3/2)) (a = const.) as t → ∞ for the trapping problem in 3D is indeed recovered. In fact, at high density, S(t) is dominated by short encounter times, which makes it exceedingly hard to record the events corresponding to the exploration of large, trap-free regions. As a consequence, at high densities the steady-state rate simply tends to 1/τ0. Finally, we work out an analytical formula for the rate that shows a remarkable agreement with the numerics up φ = 0.4.

Generalization of Wilemski-Fixman-Weiss decoupling approximation to the case involving multiple sinks of different sizes, shapes, and reactivities

We generalize the Wilemski-Fixman-Weiss decoupling approximation to calculate the transient rate of absorption of point particles into multiple sinks of different sizes, shapes, and reactivities. As an application we consider the case involving two spherical sinks. We obtain a Laplace-transform expression for the transient rate that is in excellent agreement with computer simulations. The long-time steady-state rate has a relatively simple expression, which clearly shows the dependence on the diffusion constant of the particles and on the sizes and reactivities of sinks, and its numerical result is in good agreement with the known exact result that is given in terms of recursion relations.

An efficient molecular dynamics simulation method for calculating the diffusion-influenced reaction rates

We present a molecular dynamics (MD) simulation method for calculating the diffusion-influenced reaction rates in the limit of low reactant concentrations. To calculate the reaction rate coefficient, we use MD trajectories of a nonreactive equilibrium system that are initiated with a pair of reactant molecules in reactive configuration. Hence reaction systems involving complicated reactant molecules with geometrically restricted reactivities can be treated with comparable efficiency as the simple hard-sphere reaction system. Compared to the similar MD method proposed by Van Beijeren, Dong, and Bocquet [J. Chem. Phys. 114, 6265 (2001)], the present method has a couple of advantages. First, reactions involving more general sink functions can be treated. Second, more accurate results can be obtained when the reaction probability upon collision is less than unity. As an application, we investigate the effects of nondiffusive dynamics and hydrodynamic interaction of reactants on the reaction rate.