# Geminate Recombination in the Presence of Scavengers: The New Vision of the Old Problem

Authors: S. G. Fedorenko; A. A. Kipriyanov; A. B. Doktorov
Journal: Physics Research International
 (2011)

---

## Abstract

A geminate reaction betweenA and B reactants affected by the bulk reaction between B and scavengers C is discussed. The exact solution of the problem
obtained recently for randomly walking B reactant (excess electron) is compared
with the superposition approximation commonly used to interpret experimental
data. Distinctions related to the important role of time correlations between geminate
and bulk reactions are analyzed. The largest deviations exponentially growing
in time are observed for geminate reaction rate in the presence of scavengers. It
is shown that superposition approximation can decrease essentially the ultimate
probability of geminate recombination. The difference is great enough to lead to
qualitatively incorrect description of the experiment even at small concentration of
scavengers. This, in turn, may give rise to considerable errors in the determination
of geminate pair parameters or, alternatively, to the wrong information about bulk
kinetics of electron scavenging.


---

## Body

## 1. Introduction

More than fifty years ago Monchick [1] formulated the problem of the influence of scavengers with a uniform concentration [C] in the bulk on the evolution of a geminate pair formed by photo  or pulse radiation in solution. This problem subsequently called a “scavenger problem” [2–11] in the most simple case may be described as follows. After thermalization in chemically inert solvent, A and B radicals or ion-radicals produced from a “parent" molecule form geminate (spatially correlated) pairs [12–17]. The subsequent reaction  (1)A+B⟶P
is divided into two stages: geminate and bulk recombinations. Geminate reaction proceeds between spatially correlated radicals formed from one and the same parent molecule. As a rule, this takes several nanoseconds [18, 19]. Bulk reaction is recombination of radicals that left initial geminate pairs to find themselves in the bulk. At small concentration of such pairs this stage differs essentially from a geminate one in the time scale and is commonly observed on the times from micro- to milliseconds. Ignoring this stage enables one to treat reaction (1) as recombination of isolated geminate pairs [20–22].Under the above conditions, the simplest bulk reaction competing with geminate one is possible on addition of uniformly distributedC acceptors of one (say B) type of radicals to the solution. Usually C reactants are in excess (as compared to Bs), thus their concentration remains almost unchanged  (2)B+C⟶C+D.
Note that the excess of C reactants is no barrier to taking into account solely pair encounters of B with C reactants in considering small concentrations of scavengers [C]. The problem consists in calculating the decay kinetics of A or B reactants (or the formation kinetics of P or D products of geminate (1) or bulk (2) reactions, resp.), and the dependence of P or D product yield on C acceptors concentration in such a two-stage reaction.In his first work Monchick [1] allowed for the influence of scavengers by introducing extraterm proportional to constant decay rate on acceptors [C]k into the Smoluchowski equation defining the evolution of isolated geminate pair:  (3)∂tG(r⃗,t∣r⃗′,0)=DAB∇r⃗[∇r⃗G+GkBT∇r⃗Φ]-[C]kG,
where DAB=DA+DB is the mutual diffusion coefficient and Φ is a potential of mean force between two radicals. Geminate recombination is defined by partially absorbing boundary condition at the contact radius, while k is a bimolecular rate constant of bulk reaction (2) depending on mutual mobility of B and C reactants. Physically such a simple description of the scavenging effect has been justified by the assumption that diffusion profile of a geminate pair is gradually reduced in time, but spatial distribution is not affected by the presence of scavengers. Simple transformation  (4)G(r⃗,t∣r⃗′,0)=exp⁡(-[C]kt)GAB(r⃗,t∣r⃗′,0)
converts (3) to the Smoluchowski equation in the absence of scavengers for GAB. In fact, this means that reactions (1) and (2) are considered to be independent of one another. This assumption of the independent course of geminate and bulk reactions is the basic postulate of conventional phenomenological theory developed further in a number of papers [2–10]. Based on this assumption, experimentally important quantity, total scavenging probability 𝒫D(t→∞) (that is defined as the probability that B radical react with a scavenger before recombination with A radical), can be related to the survival probability of a geminate pair as follows:  (5)PD(t⟶∞)≃[C]k∫0∞dtΩ(t)e-[C]kt.
This equation shows that in phenomenological theory the total scavenging probability is proportional to the Laplace transform of the geminate pair survival probability in the absence of scavengers Ω(t). Relation (4) has become a theoretical basis for the development of powerful experimental method for the investigation of geminate recombination kinetics by concentration dependence measurements of total scavenging yield. The devised method, along with the method of external electric field, has become a very important tool for studying the kinetics of ion-radical processes. Its extensive applications is related to the discovery of short-lived particles, in particular, excess electron. This method described in detail in textbooks and reviews [21, 23–26], in fact, compares with the method of investigating quantum yield of luminescence concentration quenching both in ideology and significance. Despite recent progress of the new methods of direct observation of geminate kinetics and excess electron capture scavenger method, as a universal way of affecting geminate recombination by bulk reactants, is still of importance. It has found application in identification of kinetic mechanism of multistage bulk-geminate reactions in various systems [27–30].At the same time, only one attempt has been made in the literature to substantiate the phenomenological theory based on relations (3)–(5). In paper [11] the phenomenological theory (corrected for nonstationary course of bulk reaction) was confirmed by a more consistent many-particle theory based on superposition decoupling of three-particle correlation in the reduced distribution function at small concentrations of acceptors. However, even in the binary approach the use of superposition decoupling in many-particle derivation of kinetic equations for some reactions proved to be incorrect and was criticized in a number of works [31–33]. Recently the exactly solvable many-particle model has been proposed for the “scavenger problem” in order to take into account time correlations between bulk and geminate reactions for the first time [34]. This model allows us to obtain exact kinetic equations for the general type of A and C reactants mobility and the Coulomb interaction in a geminate pair [35]. The most essential assumption of the model is that B particle of a geminate pair (excess electron) reacting with scavengers C moves by infinitely large jumps (the hopping mechanism [36–42]). The hopping mechanism has been used more than once to describe the scavenging of an excess electron moving by random walks in solutions [40–44]. This assumption of the large jumps is justified by the well-known fact that mean free path of excess electrons is fairly large in some solvent of low polarity [23, 43, 45] (≈60–200 Å). In Letters [46, 47] two simple microscopic models of reactions (1) and (2) admitting analytical solution of exactly solvable model equations [34] were examined. It has been shown that the role of time correlations in the course of reactions (1) and (2) may be significant even at small concentration of scavengers [C]. The analysis of relative deviations between the exact solution of the “scavenger problem” and superposition approximation shows that superposition approximation substantially underestimates the rate of geminate recombination and the yield of geminate reaction product in the presence of scavengers. Thus it is very important to investigate the parameters value areas where correlations are so significant that approximate theories become unsuitable for interpretation of experimental evidence.The outline of the paper is as follows: the next section presents conventional theory and superposition approximation. In Section3 two simple microscopic models of the reacting system are described, and calculations of experimentally measured quantities by the recipes of the above approximate theories are performed. In the first model chosen (reacting “black” spheres, the hopping motion of B and C reactants) the bulk reaction kinetics develop exponentially in time, thus the results of both approximate theories become kinetically identical. We ignore deliberately non-stationarity effects responsible for additional correlations in the system, and focus our attention on the simplest case (most favourable for approximate theories) where B reactant decays at a constant rate on acceptors. The second microscopic model refers to a more realistic case of diffusion motion of scavengers in liquid solutions, and, therefore, contains additional correlations related to the presence of diffusion nonstationary part of decay kinetics on scavengers. In Section 4 we define the exactly solvable many-particle problem, and give analytical solution for hopping and diffusion motion of scavengers. Section 5 presents the detailed analysis of relative deviations of approximate theories from the exactly solvable model in the area of small concentrations of scavengers.


## 2. Conventional Theory and Superposition Approximation

Commonly,𝒫B(t) (survival probability of B reactant) and 𝒫P(t) (accumulation kinetics of geminate reaction products) are quantities measured experimentally. The survival probability 𝒫A(t) of A reactant, the kinetics of scavenging, 𝒫D(t), can be found by using the normalization conditions [34, 35]:  (6)PA(t)+PP(t)=1,PB(t)+PD(t)+PP(t)=1.
The last relation is simplified in the stationary limit t→∞, thus the experimentally observed ultimate probability of recombination 𝒫P(t→∞) and total scavenging yield 𝒫D(t→∞) obey the relation:  (7)PD(t⟶∞)+PP(t⟶∞)=1.
Besides, recently, a new method for the determination of the bulk reaction rate constant has been proposed which is based on direct experimental measurements of the formation rate of geminate recombination products in the presence of bulk scavengers 𝒫̇P(t) [48].Conventional theory [1–10] enters constant decay rate on acceptors [C]k into microscopic equations of a geminate pair, therefore, the influence of competing bulk reaction is defined only by this parameter that is the product of the acceptor concentration [C] and stationary rate constant k of B particles decay on these acceptors. Thus we have the recipe for the calculation of experimentally measured quantities:  (8)PP(t)≃PPct(t)≡∫0tdτKg(τ)exp⁡(-[C]kτ),(9)PB(t)≃PBct(t)≡Ω(t)exp⁡(-[C]kt).
As usual, geminate reaction rate in the absence of acceptors Kg(t) is expressed in terms of the elementary rate wAB(r⃗):(10)Kg(t)=∫dr⃗′∫dr⃗wAB(r⃗)GAB(r⃗,t∣r⃗′,0)P(r⃗′),
where P(r⃗′) is the initial distribution of reactants in geminate pairs, and GAB(r⃗,t∣r⃗′,0) is the probability density to find A and B reactants at the relative distance r⃗ at the instant of time t if at the initial moment t=0 they were at the distance r⃗′ (the Green function). It obeys the equation  (11)∂tGAB(r⃗,t∣r⃗′,0)=[L̂AB-wAB(r⃗)]GAB(r⃗,t∣r⃗′,0)
with the initial condition  (12)GAB(r⃗,0∣r⃗′,0)=δ(r⃗-r⃗′).
In (11), ℒ̂AB is the integral operator defining Markovian random walks of A and B in relative coordinates of the pair AB. The force interaction is included.The survival probability of geminate pairs in the absence of competing bulk reaction is defined by the expression:(13)Ω(t)=1-∫0tdτKg(τ).
Superposition approximation [11] adjusts phenomenological theory in view of non-stationarity of the bulk kinetics. As a result, accumulation kinetics of geminate reaction products is obtained by integrating the product of the bulk reaction kinetics and pure geminate reaction rate:  (14)PP(t)≃PPsp(t)≡∫0tdτKg(τ)N(τ),
while B reactant decay kinetics is expressed as the product of bulk and geminate kinetics:  (15)PB(t)≃PBsp(t)≡N(t)Ω(t).
Note that superposition approximation actually substantiates the validity of papers [49, 50] wherein formulae (14), (15) were introduced intuitively.For dilute solution of reactants the bulk reaction kinetics may be expressed as [21, 22]:  (16)N(t)=exp⁡{-[C]∫0tdτK(τ)}.
Note that the expression in exponent (16) is linear in concentration of acceptors [C], since it allows solely for binary collisions of reactants in the bulk. In the general case of arbitrary mobilities of B and C reactants, taking account of many-particle collisions results in nonlinear corrections in concentration corresponding to the competition of acceptors with each other. Only in the case of immobile B particle (“target problem” [51, 52]) nonlinear corrections become zero, and formula (16) gives the exact many-particle description of the bulk reaction.In the absence of initial correlations between reactants the rate constant of the bulk reaction is(17)K(t)=∫dr⃗′∫dr⃗wBC(r⃗)GBC(r⃗,t∣r⃗′,0).
Stationary rate constant of the decay on acceptors k appearing in formulae (8) and (9) is determined from (17):  (18)k=lim⁡t→∞K(t).
The Green function GBC(r⃗,t∣r⃗′,0) of the pair BC is defined by the equation  (19)∂tGBC(r⃗,t∣r⃗′,0)=[L̂BC-wBC(r⃗)]GBC(r⃗,t∣r⃗′,0)
with the initial condition of type (12). As in (11), ℒ̂BC is the integral operator describing translational motion of reactants in relative coordinates of the reaction pair, and bulk reaction between B and C reactants is specified by the elementary rate wBC(r⃗) which depends on their relative position.Thus in both approximate theories the experimentally measured quantities (8), (9), and (14), (15) are expressed in a simple way in terms of pure geminate and bulk reactions kinetics. This indicates that, in fact, both conventional theory and superposition approximation consider these reactions as proceeding independently.


## 3. Microscopic Model of the “Scavenger Problem”

We believe the above independence of bulk and geminate reactions in the “scavenger problem” to be unjustified at least for two physical reasons. First, on intrusion into a geminate pair, an acceptor is to affect essentially the course of geminate reaction, with the most significant changes being determined by the interruption of recontacts betweenA and B reactants. Second, the situation of highly mobile B reactant and relatively slow C reactants is possible. In this case the entry of reactants into the bulk results in arresting geminate reaction with the subsequent initiation of the bulk one. In other words, generally speaking, the course of both reactions should be correlated in time.


### 3.1. Geminate Reaction

For further investigation we choose the simplest microscopic model of the reacting system. We shall describe a geminate reaction by the model of isotropic “black” sphere of the radiusRAB [53]. The “black” sphere model excludes recontacts of reactants and, consequently, related time correlations in the system.Let us take that unchargedA reactant moves by continual diffusion, while the motion of B reactant (excess electron) is described by the model of infinite jumps with the mean frequency τB-1 (the so-called two-scale migration [54]):  (20)L̂AB=DA∇r⃗2-τB-1Î.
Here Î is a unit operator.Starting with the initial distancer0(21)P(r⃗)=P(r)=δ(r-r0)4πrr0,
we easily obtain that at r0>RAB geminate reaction rate in the absence of bulk acceptors is defined as(22)Kg(t,r0)=K̃g(t,r0)e-t/τB,K̃g(t,r0)=RABr0τ̃g4πt3exp⁡(-τ̃g4t),
where τ̃g=(r0-RAB)2/DA is the characteristic decay time of the geminate pair, and K̃g is the geminate reaction rate for the case of immobile B reactant. Further, quantities related to the problem with immobile B particle will be marked by a tilde.By analogy withτ̃g we introduce the characteristic decay time of the geminate pair τg for the case of two-scale migration defining it as the characteristic decay time of the function Kg(t,r0). Obviously, it is proportional to the time moment tm of the maximum of this function τg=λtm. The proportionality coefficient λ is determined from the coincidence condition of τg and τ̃g at τB→∞. It is easily seen that λ=6. So for τg, (23)τg=32τB(9+4τ̃gτB-3)≈3τ̃gτB,τB≪τ̃g.
The use of (22) in (13) with the subsequent integration gives the expression for the survival probability:  (24)Ω(t,r0)=1-RAB2r0[e-ηerfc⁡(τ̃g4t-tτB)+eηerfc⁡(τ̃g4t+tτB)],
where the introduced parameter (25)η=τ̃gτB=(r0-RAB)2DAτB
is defined by the ratio of τ̃g to B reactant residence time in the pair τB. Since geminate reaction is arrested immediately after the first jump of B reactant, the depth of the reaction course depends noticeably on the parameter η:  (26)Ω(∞,r0)≡Ω∞(r0)=1-RABr0exp⁡(-η).
Only in the limit η→0 it is equal to the probability of the diffusion entry of reactants into the bulk: 1-RAB/r0. In this limiting case of immobile B(τB→∞) expression (24) is easily seen to turn into the well-known diffusion result [21]:  (27)lim⁡τB→∞Ω(t,r0)≡Ω̃(t,r0)=1-RABr0erfc⁡(τ̃g4t).
Finally, it is worth noting that our microscopic model neglects the Coulomb interaction in a geminate pair and considers the case of the photodetachment process of excess electron generation in solution. However, our results can readily be applied to the description of charged geminate pairs in water-like solutions for which the Onsager radius |rc|≤10 Å. As follows from [55, 56], in this case it will suffice to perform the following substitution:  (28)τ̃g⟶τ̃gc=rc24DA[coth⁡(rc2r0)-coth⁡(rc2RAB)]2,RABr0⟶reff(RAB)reff(r0),reff(r)=|rc|1-exp⁡(-|rc|/r),
because in highly polar solvents the survival probability of charged geminate pairs Ω̃c(t,r0) is kinetically identical to the corresponding value Ω̃(t,r0) without the Coulomb interaction.


### 3.2. Hopping Motion of Scavengers

When calculating the bulk reaction kinetics, in this subsection we shall take thatC reactants (just as B) move by jumps of infinite length:  (29)L̂BC=-(τB-1+τC-1)Î.
The choice of this model of motion excludes recontacts determined by the kinematics of reactants approach thus reducing the value of time correlations in the system.The model involving jumps of infinite length describes the so-called hopping mechanism of reactions. Physically it is realized when a reactant reaches the reaction zone in one jump. In this case the rate constant is equal to the product of the frequency of jumps and the reaction volume. Being an obvious alternative to diffusion mechanism, the hopping mechanism was first proposed in papers [36, 37] to describe luminescence impurity quenching experiments in solid solutions. For liquid phase reaction it is developed in papers [38–42] and was first found experimentally in excess electron capture reactions in nonpolar liquids [43].The model of the “black” ball of the radiusRBC insight of wich the reaction proceeds instantaneously is used for the calculation of the bulk kinetics N(t). For the hopping mechanism this model is the analog of the commonly used “black” sphere model. In this case the initial nonstationary static stage of the kinetics turns into a point (“instantaneous quenching” [57]), and the bulk reaction proceeds exponentially. Calculations with ℒ̂BC in form (29) yield  (30)N(t)=exp⁡[-ξ-ξ(τB-1+τC-1)t],
where the dimensionless parameter ξ=[C]V (the density parameter) appears. It is equal to the average number of acceptors in the reaction zone of the volume V=(4π/3)RBC3. Here we neglect the intrinsic volume of reactants and refer to long-range reactions.Kinetics (30) consists of two cofactors. The first co-factor, e-ξ, defines the depth of nonstationary stage and is the probability that at the initial instant of time there are no acceptors in the reaction volume. The second co-factor describes the exponential decay at a constant rate equal to the product of the total frequency of jumps τB-1+τC-1 and the average number of acceptors ξ entering the reaction zone. Note that in the limit of immobile B reactant (τB→∞)N(t) tends to the exact many-particle kinetics of the target problem Ñ(t) that for the hopping motion of scavengers is(31)Ñ(t)=exp⁡(-ξ-ξtτC).
Passing to experimentally measured quantities of the “scavenger problem”, note that for both approximate theories the expression for B reactant decay kinetics is obtained immediately after substitution of (24) and (30) in formulae (9) and (15), while calculation of the accumulation kinetics of geminate reaction products calls in addition for integration over time in (8) and (14). Nevertheless, we easily derive for the hopping motion of scavengers:  (32)PPsp(t)=e-ξPPct(t),PBsp(t)=e-ξPBct(t).
The above relations show that in the microscopic model chosen phenomenological theory (8), (9) becomes kinetically identical to the result of superposition approximation (14), (15).The explicit expression for the accumulation kinetics of geminate reaction product in superposition approximation is deduced by substitution of (22), (30) in  (14) and integration:  (33)PPsp(t,r0)=RABr0e-ξφ(η,γ+1,τ),
where γ=ξ(τB/τC+1) and τ=t/τB are bulk reaction rate and time in τB units, respectively.Formula (33) involves the function φ(η,x,τ):  (34)φ(η,x,τ)=12{e-ηxerfc⁡(η4τ-xτ)+eηxerfc⁡(η4τ+xτ)},
which, in accordance with the physical meaning of the product accumulation kinetics, is a monotonically increasing function that attains its steady-state value at rather long times:  (35)φ(η,x,τ⟶∞)=exp⁡(-ηx).
This gives the expression for the ultimate probability of geminate reaction in the presence of bulk acceptors:  (36)PPsp(t⟶∞,r0)=RABr0e-ξexp⁡(-η(γ+1)).
Comparison between (36) and (26) shows that the presence of bulk acceptors leads to decrease in the product yield of geminate reaction. The kinetics of B reactant 𝒫Bsp in accordance with the independence of geminate and bulk reactions in superposition approximation mentioned above is defined by (15), (24), and (30).


### 3.3. Diffusion Motion of Scavengers

The more realistic model, we shall take here for the bulk reaction thatC reactants move by continual diffusion with the diffusion coefficient DC, while B reactant (excess electron), as before, moves by infinite jumps with the mean time τB between them:  (37)L̂BC=DC∇r⃗2-τB-1Î.
During this time interval τB between jumps the bulk reaction develops due to acceptors diffusion to immobile electron (the “target problem” kinetics Ñ(t)), while a regular very large electron jump changes the acceptors surroundings in an uncorrelated way, and the process is repeated until the two-scale bulk kinetics N(t) is formed.The model of a “black” ball of the radiusRBC is used for further calculations. As in the previous subsection, in this case the initial nonstationary static stage of the kinetics turns into a point (“instantaneous quenching” [57]). It is obvious that for τB-1=0 the “target problem” kinetics Ñ(t) proceeds as the well-known Smoluchowski kinetics [53]:  (38)Ñ(t)=exp⁡[-[C](V+k̃t+23πVk̃t)],k̃=4πRBCDC,
where V=(4π/3)RBC3 is the reaction zone volume. The density parameter ξ=[C]V equal to the average number of acceptors in the reaction zone of the volume V can be introduced as in hopping kinetics (30).The first co-factor of (38),  e-ξ, is the same as in (30), the second co-factor describes the exponential decay with the diffusion rate constant k̃, and the last one is the diffusion nonstationary correction. In this subsection we use the binary approximation of the bulk kinetics [58], which allows us to expand the nonstationary part of the kinetics with the accuracy of a linear term over acceptor concentration [C]:  (39)Ñ(t)≈e-[C]k̃t[1-ξ-23πξ[C]k̃t].
The validity condition of formula (39) defines the binary time interval [C]k̃t≪ξ-1 in which time-dependent observable values will be considered.The kinetics of the bulk reactionN(t) (16) can be obtained by solving (19) with operator (37) in the “black” sphere limit. In the binary approximation it takes the form(40)N(t)≈e-[C]kt[1-ξ-2f(tτB)3πξ[C]k̃t],
where, in contrast to the case of diffusion mobilities of both reactants, stationary rate constant of the bulk two-scale reaction is(41)k=k̃+VτB+3k̃VτB,
and f(τ) is a function decreasing from unity to zero:  (42)f(τ)=12[e-τ+(π4τ+πτ)erf⁡(τ)-πτ].
At the time t≫τB it goes to its asymptotic f(τ)≈(π/τ)/4, and the nonstationary part of bulk kinetics becomes time independent.Stationary rate constant (41) is equal to the sum of diffusion rate constant k̃ (38), hopping rate constant V/τB [36–42] and the interference term appearing due to the intricate character of reactants approach. The presence of this term is a distinguishing feature of stationary rate constant (41) of the two-scale bulk reaction. It has been already obtained in [54]. The second feature, as follows from our result (40), consists in the decrease of the nonstationary part of the kinetics as compared to (39) owing to the decreasing function f(t/τB).The expression for the accumulation kinetics of geminate reaction product in the superposition approximation is deduced by substitution of (22), (40) in (14) and integration(43)PPsp(t,r0)=(1-ξ)PP(m)sp(t,r0)-PP(n)sp(t,r0),
where the Markovian part of the kinetics corresponding to pure exponential term of (40) is:  (44)PP(m)sp(t,r0)=RABr0φ(η,γ+1,τ),γ=[C]kτB.
As in the previous subsection, the introduced parameter γ is the bulk reaction rate in τB units. Using (41), it can be rewritten in our case of diffusion motion of acceptor as γ=β+ξ+3ξβ, where β=[C]k̃τB.The non-Markovian part of the kinetics is defined by the integral(45)PP(n)sp(t,r0)=RABr03πξηβ∫0(γ+1)τdxx×exp⁡(-x-η(γ+1)4x)f(xγ+1).
This integral is not taken analytically even in the stationary limit τ→∞.In accordance with (15), the kinetics 𝒫Bsp of B reactant is defined by the product of (24) and (40).


## 3.1. Geminate Reaction

For further investigation we choose the simplest microscopic model of the reacting system. We shall describe a geminate reaction by the model of isotropic “black” sphere of the radiusRAB [53]. The “black” sphere model excludes recontacts of reactants and, consequently, related time correlations in the system.Let us take that unchargedA reactant moves by continual diffusion, while the motion of B reactant (excess electron) is described by the model of infinite jumps with the mean frequency τB-1 (the so-called two-scale migration [54]):  (20)L̂AB=DA∇r⃗2-τB-1Î.
Here Î is a unit operator.Starting with the initial distancer0(21)P(r⃗)=P(r)=δ(r-r0)4πrr0,
we easily obtain that at r0>RAB geminate reaction rate in the absence of bulk acceptors is defined as(22)Kg(t,r0)=K̃g(t,r0)e-t/τB,K̃g(t,r0)=RABr0τ̃g4πt3exp⁡(-τ̃g4t),
where τ̃g=(r0-RAB)2/DA is the characteristic decay time of the geminate pair, and K̃g is the geminate reaction rate for the case of immobile B reactant. Further, quantities related to the problem with immobile B particle will be marked by a tilde.By analogy withτ̃g we introduce the characteristic decay time of the geminate pair τg for the case of two-scale migration defining it as the characteristic decay time of the function Kg(t,r0). Obviously, it is proportional to the time moment tm of the maximum of this function τg=λtm. The proportionality coefficient λ is determined from the coincidence condition of τg and τ̃g at τB→∞. It is easily seen that λ=6. So for τg, (23)τg=32τB(9+4τ̃gτB-3)≈3τ̃gτB,τB≪τ̃g.
The use of (22) in (13) with the subsequent integration gives the expression for the survival probability:  (24)Ω(t,r0)=1-RAB2r0[e-ηerfc⁡(τ̃g4t-tτB)+eηerfc⁡(τ̃g4t+tτB)],
where the introduced parameter (25)η=τ̃gτB=(r0-RAB)2DAτB
is defined by the ratio of τ̃g to B reactant residence time in the pair τB. Since geminate reaction is arrested immediately after the first jump of B reactant, the depth of the reaction course depends noticeably on the parameter η:  (26)Ω(∞,r0)≡Ω∞(r0)=1-RABr0exp⁡(-η).
Only in the limit η→0 it is equal to the probability of the diffusion entry of reactants into the bulk: 1-RAB/r0. In this limiting case of immobile B(τB→∞) expression (24) is easily seen to turn into the well-known diffusion result [21]:  (27)lim⁡τB→∞Ω(t,r0)≡Ω̃(t,r0)=1-RABr0erfc⁡(τ̃g4t).
Finally, it is worth noting that our microscopic model neglects the Coulomb interaction in a geminate pair and considers the case of the photodetachment process of excess electron generation in solution. However, our results can readily be applied to the description of charged geminate pairs in water-like solutions for which the Onsager radius |rc|≤10 Å. As follows from [55, 56], in this case it will suffice to perform the following substitution:  (28)τ̃g⟶τ̃gc=rc24DA[coth⁡(rc2r0)-coth⁡(rc2RAB)]2,RABr0⟶reff(RAB)reff(r0),reff(r)=|rc|1-exp⁡(-|rc|/r),
because in highly polar solvents the survival probability of charged geminate pairs Ω̃c(t,r0) is kinetically identical to the corresponding value Ω̃(t,r0) without the Coulomb interaction.


## 3.2. Hopping Motion of Scavengers

When calculating the bulk reaction kinetics, in this subsection we shall take thatC reactants (just as B) move by jumps of infinite length:  (29)L̂BC=-(τB-1+τC-1)Î.
The choice of this model of motion excludes recontacts determined by the kinematics of reactants approach thus reducing the value of time correlations in the system.The model involving jumps of infinite length describes the so-called hopping mechanism of reactions. Physically it is realized when a reactant reaches the reaction zone in one jump. In this case the rate constant is equal to the product of the frequency of jumps and the reaction volume. Being an obvious alternative to diffusion mechanism, the hopping mechanism was first proposed in papers [36, 37] to describe luminescence impurity quenching experiments in solid solutions. For liquid phase reaction it is developed in papers [38–42] and was first found experimentally in excess electron capture reactions in nonpolar liquids [43].The model of the “black” ball of the radiusRBC insight of wich the reaction proceeds instantaneously is used for the calculation of the bulk kinetics N(t). For the hopping mechanism this model is the analog of the commonly used “black” sphere model. In this case the initial nonstationary static stage of the kinetics turns into a point (“instantaneous quenching” [57]), and the bulk reaction proceeds exponentially. Calculations with ℒ̂BC in form (29) yield  (30)N(t)=exp⁡[-ξ-ξ(τB-1+τC-1)t],
where the dimensionless parameter ξ=[C]V (the density parameter) appears. It is equal to the average number of acceptors in the reaction zone of the volume V=(4π/3)RBC3. Here we neglect the intrinsic volume of reactants and refer to long-range reactions.Kinetics (30) consists of two cofactors. The first co-factor, e-ξ, defines the depth of nonstationary stage and is the probability that at the initial instant of time there are no acceptors in the reaction volume. The second co-factor describes the exponential decay at a constant rate equal to the product of the total frequency of jumps τB-1+τC-1 and the average number of acceptors ξ entering the reaction zone. Note that in the limit of immobile B reactant (τB→∞)N(t) tends to the exact many-particle kinetics of the target problem Ñ(t) that for the hopping motion of scavengers is(31)Ñ(t)=exp⁡(-ξ-ξtτC).
Passing to experimentally measured quantities of the “scavenger problem”, note that for both approximate theories the expression for B reactant decay kinetics is obtained immediately after substitution of (24) and (30) in formulae (9) and (15), while calculation of the accumulation kinetics of geminate reaction products calls in addition for integration over time in (8) and (14). Nevertheless, we easily derive for the hopping motion of scavengers:  (32)PPsp(t)=e-ξPPct(t),PBsp(t)=e-ξPBct(t).
The above relations show that in the microscopic model chosen phenomenological theory (8), (9) becomes kinetically identical to the result of superposition approximation (14), (15).The explicit expression for the accumulation kinetics of geminate reaction product in superposition approximation is deduced by substitution of (22), (30) in  (14) and integration:  (33)PPsp(t,r0)=RABr0e-ξφ(η,γ+1,τ),
where γ=ξ(τB/τC+1) and τ=t/τB are bulk reaction rate and time in τB units, respectively.Formula (33) involves the function φ(η,x,τ):  (34)φ(η,x,τ)=12{e-ηxerfc⁡(η4τ-xτ)+eηxerfc⁡(η4τ+xτ)},
which, in accordance with the physical meaning of the product accumulation kinetics, is a monotonically increasing function that attains its steady-state value at rather long times:  (35)φ(η,x,τ⟶∞)=exp⁡(-ηx).
This gives the expression for the ultimate probability of geminate reaction in the presence of bulk acceptors:  (36)PPsp(t⟶∞,r0)=RABr0e-ξexp⁡(-η(γ+1)).
Comparison between (36) and (26) shows that the presence of bulk acceptors leads to decrease in the product yield of geminate reaction. The kinetics of B reactant 𝒫Bsp in accordance with the independence of geminate and bulk reactions in superposition approximation mentioned above is defined by (15), (24), and (30).


## 3.3. Diffusion Motion of Scavengers

The more realistic model, we shall take here for the bulk reaction thatC reactants move by continual diffusion with the diffusion coefficient DC, while B reactant (excess electron), as before, moves by infinite jumps with the mean time τB between them:  (37)L̂BC=DC∇r⃗2-τB-1Î.
During this time interval τB between jumps the bulk reaction develops due to acceptors diffusion to immobile electron (the “target problem” kinetics Ñ(t)), while a regular very large electron jump changes the acceptors surroundings in an uncorrelated way, and the process is repeated until the two-scale bulk kinetics N(t) is formed.The model of a “black” ball of the radiusRBC is used for further calculations. As in the previous subsection, in this case the initial nonstationary static stage of the kinetics turns into a point (“instantaneous quenching” [57]). It is obvious that for τB-1=0 the “target problem” kinetics Ñ(t) proceeds as the well-known Smoluchowski kinetics [53]:  (38)Ñ(t)=exp⁡[-[C](V+k̃t+23πVk̃t)],k̃=4πRBCDC,
where V=(4π/3)RBC3 is the reaction zone volume. The density parameter ξ=[C]V equal to the average number of acceptors in the reaction zone of the volume V can be introduced as in hopping kinetics (30).The first co-factor of (38),  e-ξ, is the same as in (30), the second co-factor describes the exponential decay with the diffusion rate constant k̃, and the last one is the diffusion nonstationary correction. In this subsection we use the binary approximation of the bulk kinetics [58], which allows us to expand the nonstationary part of the kinetics with the accuracy of a linear term over acceptor concentration [C]:  (39)Ñ(t)≈e-[C]k̃t[1-ξ-23πξ[C]k̃t].
The validity condition of formula (39) defines the binary time interval [C]k̃t≪ξ-1 in which time-dependent observable values will be considered.The kinetics of the bulk reactionN(t) (16) can be obtained by solving (19) with operator (37) in the “black” sphere limit. In the binary approximation it takes the form(40)N(t)≈e-[C]kt[1-ξ-2f(tτB)3πξ[C]k̃t],
where, in contrast to the case of diffusion mobilities of both reactants, stationary rate constant of the bulk two-scale reaction is(41)k=k̃+VτB+3k̃VτB,
and f(τ) is a function decreasing from unity to zero:  (42)f(τ)=12[e-τ+(π4τ+πτ)erf⁡(τ)-πτ].
At the time t≫τB it goes to its asymptotic f(τ)≈(π/τ)/4, and the nonstationary part of bulk kinetics becomes time independent.Stationary rate constant (41) is equal to the sum of diffusion rate constant k̃ (38), hopping rate constant V/τB [36–42] and the interference term appearing due to the intricate character of reactants approach. The presence of this term is a distinguishing feature of stationary rate constant (41) of the two-scale bulk reaction. It has been already obtained in [54]. The second feature, as follows from our result (40), consists in the decrease of the nonstationary part of the kinetics as compared to (39) owing to the decreasing function f(t/τB).The expression for the accumulation kinetics of geminate reaction product in the superposition approximation is deduced by substitution of (22), (40) in (14) and integration(43)PPsp(t,r0)=(1-ξ)PP(m)sp(t,r0)-PP(n)sp(t,r0),
where the Markovian part of the kinetics corresponding to pure exponential term of (40) is:  (44)PP(m)sp(t,r0)=RABr0φ(η,γ+1,τ),γ=[C]kτB.
As in the previous subsection, the introduced parameter γ is the bulk reaction rate in τB units. Using (41), it can be rewritten in our case of diffusion motion of acceptor as γ=β+ξ+3ξβ, where β=[C]k̃τB.The non-Markovian part of the kinetics is defined by the integral(45)PP(n)sp(t,r0)=RABr03πξηβ∫0(γ+1)τdxx×exp⁡(-x-η(γ+1)4x)f(xγ+1).
This integral is not taken analytically even in the stationary limit τ→∞.In accordance with (15), the kinetics 𝒫Bsp of B reactant is defined by the product of (24) and (40).


## 4. Exactly Solvable Model

The choice of the hopping mobility of an excess electron in the above microscopic description of a reacting system makes it possible to use the exactly solvable model of the “scavenger problem” formulated by the authors in [34]. It should be kept in mind that, since the exactly solvable many-particle model takes no account of a hard core in the description of geminate reaction, the microscopic model can be formulated with the “black” sphere with penetration instead of the one without it. This will not affect the results of the preceding section, because we always consider the case where the initial distance in a geminate pair exceeds the “black” sphere radius (r0>RAB); later on repulsive forces cannot manifest themselves due to the instantaneous decay of the pair on the approach of reactants to the distance RAB.Assuming the hopping motion ofB reactant, the exactly solvable many-particle model of the “scavenger problem” gives the recipe for the calculation of the quantities measured:  (46)ṖP(t)=Kg(τ)Ñ(τ),PP(t)=∫0tdτKg(τ)Ñ(τ),(47)PB(t)=Ω̃(t)Ñ(t)e-t/τB+1τB∫0tdτÑ(t-τ)e-(t-τ)/τBPB(τ),
where Kg(t) is the pure geminate reaction rate defined in (10), (22), Ω̃(t) is defined in (27), Ñ(t) is the exact many-particle survival probability of immobile B reactant surrounded by bulk acceptors (see  (31), (38)). The binary kinetics N(t) (16) goes to it with τB→∞.In this limit of immobileB particle recipes (46), (47) are reduced to  (48)PP(t)=∫0tdτK̃g(τ)Ñ(τ);PB(t)=Ñ(t)Ω̃(t).
Comparison of (48) with (14), (15) shows that at τB=∞ the superposition approximation coincides with the exact solution of the problem. Thus for immobile B time correlations between geminate and bulk reactions are absent in the given model at any acceptor concentrations.We begin the examination of time correlation influence with the accumulation kinetics analysis of geminate recombination products. Comparison between (46) and (14) shows that in the general case of mobile B the rate of geminate reaction products formation 𝒫̇P(t)=Kg(t)Ñ(t) cannot be expressed solely in terms of the kinetics of geminate and bulk reactions. It is the product of geminate reaction rate Kg(t) and the above probability Ñ(t) formed only during the residence time of B particle in a geminate pair before it jumps into the bulk. Thus here the use of pure bulk reaction kinetics N(t) is physically unjustified.


### 4.1. Hopping Motion of Scavengers

To analyze errors brought about by superposition approximation at arbitrary mobility ofB reactants, we make calculations in the framework of the microscopic model chosen. Comparison of the kinetics Ñ(t) (31) with formula (30) shows that the bulk reaction kinetics N(t) differs from Ñ(t) by the exponential factor:  (49)N(t)=exp⁡(-ξtτB)Ñ(t).
Accordingly, the accumulation rate 𝒫̇Psp(t) of geminate reaction products predicted by superposition approximation differs from the exact result in the same way. So the superposition approximation can decrease essentially the accumulation rate of geminate reaction products.To find the explicit form of the exact many-particle accumulation kinetics𝒫P(t) of geminate reaction products, use formulae (46), (22), and (31). We have.(50)PP(t,r0)=RABr0e-ξφ(η,ν+1,τ),ν=ξτBτC.
The introduced dimensionless parameter ν has the meaning of the bulk reaction rate with immobile B reactant in units of τB. It follows from its definition that the parameter γ from (33) tends to ν with τB tending to ∞. The function φ(η,x,τ) is defined by formula (34). As expression (33) for 𝒫Psp coincides with (50) up to the substitution of ν for γ, the ultimate probability of geminate reaction can be obtained from (36):  (51)PP(t⟶∞,r0)=RABr0e-ξexp⁡(-η(ν+1)).
For small values of the parameter ξ the result for the ratio of ultimate probabilities (36) and (51) is easily obtained:  (52)PPsp(t⟶∞,r0)PP(t⟶∞,r0)=exp⁡(-ξ2ην+1),ξ≪1.
Leaving the detailed analysis of deviations of superposition approximation from the exact solution for Section 5, note that, according to (52), these deviations may be considerable even at small values of ξ, on condition that the value of the parameter η is sufficiently large.The behaviour of the accumulation kinetics𝒫P(t,r0) and accumulation rate 𝒫̇P(t,r0) of geminate reaction products as compared to the results of superposition approximation is given in Figures  1 and 2. The superposition approximation decreases noticeably the ultimate probability of geminate reaction in the presence of scavengers, and, according to (49), overestimates the decay of 𝒫̇Psp(t,r0) as compared to the exact solution denoted by lines with circles.Figure 1
Exact accumulation kinetics𝒫P(t)—(1), and superposition approximation 𝒫Psp(t)—(2) for geminate reaction products in the presence of scavengers. ξ=0.1, η=3,RAB/r0=0.95,τB/τC=0.1.Figure 2
Exact product accumulation rate𝒫̇P(t)—(1), and superposition approximation 𝒫̇Psp(t)—(2). Parameters are the same as in Figure 1.Before passing to the decay kinetics𝒫B(t) of B reactant, note that no radical distinctions should be expected at least at low concentrations of acceptors. The reason is that in the case of rather fast motion of B, correlations are essential, and the contribution of geminate reaction into complete kinetics is insignificant. When B reactant moves slowly, and the contribution of geminate reaction is considerable, correlations in the system are negligible, and the superposition approximation is found to be suitable.To calculate the explicit form of the exact many-particle kinetics ofB reactant, we solve (47) by a conventional method of the Laplace transformation using expressions (27) and (31) for Ω̃(t,r0) and Ñ(t). This gives(53)PB(t,r0)=NB(t)[1-RABr0φ(η,e-ξ,τ)].
Here NB(t) is the exact many-particle kinetics of the bulk reaction representable as(54)NB(t)=exp⁡(-ξ-ξtτC-(1-e-ξ)tτB)≡Ñ(t)N0(t).
Formula (54) is physically clear. The bulk reaction kinetics is defined by two independent processes. The first process described by Ñ(t) (31) is the reaction between mobile acceptors and immobile B particle, while the second one, N0(t), defines the reaction of mobile B particle with the ensemble of immobile acceptors. A regular very large jump of B changes the environment of acceptors in an uncorrelated way, thus the rate of the second process is equal to the product of the frequency of jumps 1/τB and the many-particle probability of acceptors entering the reaction zone, 1-e-ξ. At small parameters ξ expansion (54) gives exactly the formula of binary kinetics (30) depending on a relative mobility of B and C reactants.It is easily seen that the structure of expression (53) coincides with the structure defined by superposition approximation recipe (15). As in the case of immobile B particle, the complete kinetics can be divided into a product of two cofactors. The first co-factor describes the many-particle kinetics of the bulk reaction NB(t) which is the solution of (47) in the absence of geminate recombination (Ω̃(t)≡1). Following the superposition approximation logic, the second co-factor may be considered to be the survival probability of some effective geminate pair making allowance for correlations with the bulk reaction:  (55)Ωeff(t,r0)=1-RABr0φ(η,e-ξ,τ).
This is also confirmed by the absence of the dependence Ωeff(t,r0) on C reactant mobility.As one would expect, the influence of correlations disappears under infinite dilution of bulk scavengers: withξ→0 expression (55) goes into (24). However, it is seen that the presence of acceptors in (55) results in the acceleration of the kinetics Ωeff(t,r0) as compared to pure geminate kinetics Ω(t,r0). Correspondingly, the stationary portion of nonreacting geminate pairs (the escape probability) is less than in (26):  (56)Ω∞eff(r0)=1-RABr0exp⁡(-ηe-ξ).
This means that in a hypothetical effective pair geminate recombination proceeds at a higher rate. Unlike (22), in such a pair the reaction rate may be represented as(57)Kgeff(t,r0)=K̃g(t,r0)exp⁡(-tτBe-ξ).
This is paradoxical: taking account of time correlations with the competing bulk reaction accelerates the reaction in effective geminate pair, though it seems obvious that it must be retarded due to additional decay channel in effective acceptor surroundings. However, the paradox arises from the superposition approximation logic. The attempt to treat the exact kinetics 𝒫B(t,r0) as the product of the bulk kinetics and some independent effective geminate kinetics is physically unjustified.Another interpretation of the rateKgeff(t,r0) (57) is given by the representation:  (58)Kgeff(t,r0)=Kg(t,r0)Ñ(t)NB(t).
It shows that the increase in Kgeff(t,r0) is related to defect of the superposition approximation that underestimates significantly geminate reaction product accumulation rate due to exponential difference between N(t) and Ñ(t) (see (49)). The dependence of Kgeff(t,r0) on the bulk reaction kinetics NB(t) suggests that geminate and bulk reactions are correlated in time. That is why the attempt to represent them as two independent processes leads to physical contradiction. Note that adequate description follows from the derivation of kinetic equations of the “scavenger problem” based on the application of general kinetic theory concepts to the exactly solvable many-particle problem [35].The behaviour of the kinetics𝒫B(t,r0) and 𝒫Bsp(t,r0) is given in Figure 3 for two values of the density parameter ξ. Note that relative deviations at small values of ξ (curves 1, 2) are not so great as for geminate reaction product accumulation. They are proportional to the density parameter. Increased distinctions beyond the limits of the binary interval (curves 3, 4) are due to the many-particle nature of the kinetics NB(t) (54).Figure 3
Exact kinetics𝒫B(t)—(1, 3) and superposition approximation 𝒫Bsp(t)—(2, 4) for different values of ξ: ξ=0.1—(1, 2); ξ=0.3—(3, 4). Other parameters: η=0.1, RAB/r0=0.95, τC/τB=0.1.


### 4.2. Diffusion Motion of Scavengers

For diffusion motion of scavengers comparison between binary formulae (39) and (40) shows that the bulk reaction kinetics N(t) differs from Ñ(t) by the exponential multiplier due to the difference between k and k̃, while the second power factor appears due to nonstationary parts of kinetics (39) and (40):  (59)N(t)Ñ(t)=exp⁡[-(ξ+3ξβ)τ][1-ξ-f(τ)ñ(τ)1-ξ-ñ(τ)],β=[C]k̃τB.
Here β and τ=t/τB are diffusion reaction rate with immobile B reactant and time in τB units, respectively. We introduce the notation for diffusion nonstationary part of kinetics (38), (39): (60)ñ(τ)=23ξβτπ.
Accordingly, the accumulation rate 𝒫̇Psp(t) of geminate reaction products predicted by superposition approximation (14) differs from the exact result in the same way. As is seen, for diffusion motions of scavengers these distinctions become more significant because of the presence of extra root term in the exponent of (59).To find the accumulation kinetics𝒫P(t) of geminate reaction products, use formulae (46), (22), and (39). We have:  (61)PP(t,r0)=(1-ξ)PP(m)(t,r0)-PP(n)(t,r0),
where 𝒫P(m) is the contribution of the Markovian (exponential) term of (39) calculated in the previous section:  (62)PP(m)(t,r0)=RABr0φ(η,β+1,τ),
where the function φ(η,x,τ) is defined in (34). The value 𝒫P(n)(t,r0) is determined by diffusion non-Markovian part of kinetics (39) and is equal to  (63)PP(n)(t,r0)=RABr03πξηβ∫0(β+1)τdxxexp⁡(-x-η(β+1)4x).
In the limit t→∞ the integral in (63) goes to the Bessel function. In this limit, using (35), we can rewrite (61) in an analytical form:  (64)PP(∞,r0)=RABr0[23π(1-ξ)e-η(β+1)-23πξηβK0(η(β+1))],
where K0(x) is the modified Bessel function of a zero order [59].The behaviour of the exact accumulation kinetics𝒫P(t) (61)–(63) and superposition approximation 𝒫Psp(t) (43)–(45) is shown in Figure 4.Figure 4
Exact accumulation kinetics𝒫P(t)—(2), its Markovian part 𝒫P(m)(t)—(4), and superposition approximation 𝒫Psp(t)—(1), its Markovian part 𝒫P(m)sp(t)—(3) for geminate reaction products in the presence of diffusing scavengers. ξ=0.1, η=3, RAB/r0=0.95, τB/τC=3DCτB/RBC2=0.3.The ratio of (43) and (61) at large τ and small ξ is recast as(65)PPsp(∞,r0)PP(∞,r0)=exp⁡[-ξ2ηβ+1(1+3βξ)]Θ,ξ≪1.
The first exponential factor in (65) is defined by the ratio between Markovian contributions (44) and (62), while Θ specifies the relation of non-Markovian parts of geminate product:  (66)Θ=1-ξ-PP(n)sp(∞,r0)/PP(m)sp(∞,r0)1-ξ-PP(n)(∞,r0)/PP(m)(∞,r0).
From (65) as well as from (52) of the previous subsection it follows that even at small ξ the result of superposition approximation can decrease essentially the ultimate probability of geminate recombination in the presence of scavengers in comparison with the exact solution if η is rather large.To calculate the decay kinetics𝒫B(t) of B reactant, we solve (47) by a conventional method of the Laplace transformation using expressions (27) and (39) for Ω̃(t,r0) and Ñ(t). This equation has a many-particle character. The correct binary approximation in the limit of ultrafast geminate reaction has to turn into the product of the escape probability and bulk kinetics: 𝒫B(t,r0)=Ω∞(r0)N(t) [35]. Taking into account this requirement, we obtain that(67)PB(t,r0)=(1-ξ)PB(m)(t,r0)-PB(n)(t,r0),
where the Markovian part of the kinetics can be represented as(68)PB(m)(t,r0)=e-[C]kt[1-RABr0φ(η,α,τ)],α=1-ξ-3ξβ.
This result coincides with the corresponding result of (53) at a small value of ξ except the parameter α. At large t expression (68) goes to the following:  (69)PB(m)(t,r0)≃e-[C]kt[1-RABr0e-αη]=e-[C]ktΩ∞eff(r0),t≫τB.
The non-Markovian part of the kinetics can be represented through quadrature of elementary functions.It is interesting to compare𝒫B and 𝒫Bsp (15), (24), (40) at a long time when B reactant leaves a geminate pair for the bulk. In this case the expressions for these kinetics are simplified, and we obtain the relation:  (70)PBsp(t,r0)PB(t,r0)≃Ω∞Ω∞eff(1-RABr0(e-αη-J(αη))3ξβ/α(1-ξ-3ξβ/2)Ω∞eff)-1,t≫τB,
where Ω∞ is the escape probability defined by (26), and J(x) is the following quadrature:  (71)J(x)=2π∫0∞dyyexp⁡(-y)erfc⁡(x4y).
It easy to see that in (70) the factor Ω∞/Ω∞eff is defined by the ratio of the Markovian parts of 𝒫Bsp and 𝒫B, while the multiplier between the brackets is the relation of the non-Markovian parts of the kinetics.


## 4.1. Hopping Motion of Scavengers

To analyze errors brought about by superposition approximation at arbitrary mobility ofB reactants, we make calculations in the framework of the microscopic model chosen. Comparison of the kinetics Ñ(t) (31) with formula (30) shows that the bulk reaction kinetics N(t) differs from Ñ(t) by the exponential factor:  (49)N(t)=exp⁡(-ξtτB)Ñ(t).
Accordingly, the accumulation rate 𝒫̇Psp(t) of geminate reaction products predicted by superposition approximation differs from the exact result in the same way. So the superposition approximation can decrease essentially the accumulation rate of geminate reaction products.To find the explicit form of the exact many-particle accumulation kinetics𝒫P(t) of geminate reaction products, use formulae (46), (22), and (31). We have.(50)PP(t,r0)=RABr0e-ξφ(η,ν+1,τ),ν=ξτBτC.
The introduced dimensionless parameter ν has the meaning of the bulk reaction rate with immobile B reactant in units of τB. It follows from its definition that the parameter γ from (33) tends to ν with τB tending to ∞. The function φ(η,x,τ) is defined by formula (34). As expression (33) for 𝒫Psp coincides with (50) up to the substitution of ν for γ, the ultimate probability of geminate reaction can be obtained from (36):  (51)PP(t⟶∞,r0)=RABr0e-ξexp⁡(-η(ν+1)).
For small values of the parameter ξ the result for the ratio of ultimate probabilities (36) and (51) is easily obtained:  (52)PPsp(t⟶∞,r0)PP(t⟶∞,r0)=exp⁡(-ξ2ην+1),ξ≪1.
Leaving the detailed analysis of deviations of superposition approximation from the exact solution for Section 5, note that, according to (52), these deviations may be considerable even at small values of ξ, on condition that the value of the parameter η is sufficiently large.The behaviour of the accumulation kinetics𝒫P(t,r0) and accumulation rate 𝒫̇P(t,r0) of geminate reaction products as compared to the results of superposition approximation is given in Figures  1 and 2. The superposition approximation decreases noticeably the ultimate probability of geminate reaction in the presence of scavengers, and, according to (49), overestimates the decay of 𝒫̇Psp(t,r0) as compared to the exact solution denoted by lines with circles.Figure 1
Exact accumulation kinetics𝒫P(t)—(1), and superposition approximation 𝒫Psp(t)—(2) for geminate reaction products in the presence of scavengers. ξ=0.1, η=3,RAB/r0=0.95,τB/τC=0.1.Figure 2
Exact product accumulation rate𝒫̇P(t)—(1), and superposition approximation 𝒫̇Psp(t)—(2). Parameters are the same as in Figure 1.Before passing to the decay kinetics𝒫B(t) of B reactant, note that no radical distinctions should be expected at least at low concentrations of acceptors. The reason is that in the case of rather fast motion of B, correlations are essential, and the contribution of geminate reaction into complete kinetics is insignificant. When B reactant moves slowly, and the contribution of geminate reaction is considerable, correlations in the system are negligible, and the superposition approximation is found to be suitable.To calculate the explicit form of the exact many-particle kinetics ofB reactant, we solve (47) by a conventional method of the Laplace transformation using expressions (27) and (31) for Ω̃(t,r0) and Ñ(t). This gives(53)PB(t,r0)=NB(t)[1-RABr0φ(η,e-ξ,τ)].
Here NB(t) is the exact many-particle kinetics of the bulk reaction representable as(54)NB(t)=exp⁡(-ξ-ξtτC-(1-e-ξ)tτB)≡Ñ(t)N0(t).
Formula (54) is physically clear. The bulk reaction kinetics is defined by two independent processes. The first process described by Ñ(t) (31) is the reaction between mobile acceptors and immobile B particle, while the second one, N0(t), defines the reaction of mobile B particle with the ensemble of immobile acceptors. A regular very large jump of B changes the environment of acceptors in an uncorrelated way, thus the rate of the second process is equal to the product of the frequency of jumps 1/τB and the many-particle probability of acceptors entering the reaction zone, 1-e-ξ. At small parameters ξ expansion (54) gives exactly the formula of binary kinetics (30) depending on a relative mobility of B and C reactants.It is easily seen that the structure of expression (53) coincides with the structure defined by superposition approximation recipe (15). As in the case of immobile B particle, the complete kinetics can be divided into a product of two cofactors. The first co-factor describes the many-particle kinetics of the bulk reaction NB(t) which is the solution of (47) in the absence of geminate recombination (Ω̃(t)≡1). Following the superposition approximation logic, the second co-factor may be considered to be the survival probability of some effective geminate pair making allowance for correlations with the bulk reaction:  (55)Ωeff(t,r0)=1-RABr0φ(η,e-ξ,τ).
This is also confirmed by the absence of the dependence Ωeff(t,r0) on C reactant mobility.As one would expect, the influence of correlations disappears under infinite dilution of bulk scavengers: withξ→0 expression (55) goes into (24). However, it is seen that the presence of acceptors in (55) results in the acceleration of the kinetics Ωeff(t,r0) as compared to pure geminate kinetics Ω(t,r0). Correspondingly, the stationary portion of nonreacting geminate pairs (the escape probability) is less than in (26):  (56)Ω∞eff(r0)=1-RABr0exp⁡(-ηe-ξ).
This means that in a hypothetical effective pair geminate recombination proceeds at a higher rate. Unlike (22), in such a pair the reaction rate may be represented as(57)Kgeff(t,r0)=K̃g(t,r0)exp⁡(-tτBe-ξ).
This is paradoxical: taking account of time correlations with the competing bulk reaction accelerates the reaction in effective geminate pair, though it seems obvious that it must be retarded due to additional decay channel in effective acceptor surroundings. However, the paradox arises from the superposition approximation logic. The attempt to treat the exact kinetics 𝒫B(t,r0) as the product of the bulk kinetics and some independent effective geminate kinetics is physically unjustified.Another interpretation of the rateKgeff(t,r0) (57) is given by the representation:  (58)Kgeff(t,r0)=Kg(t,r0)Ñ(t)NB(t).
It shows that the increase in Kgeff(t,r0) is related to defect of the superposition approximation that underestimates significantly geminate reaction product accumulation rate due to exponential difference between N(t) and Ñ(t) (see (49)). The dependence of Kgeff(t,r0) on the bulk reaction kinetics NB(t) suggests that geminate and bulk reactions are correlated in time. That is why the attempt to represent them as two independent processes leads to physical contradiction. Note that adequate description follows from the derivation of kinetic equations of the “scavenger problem” based on the application of general kinetic theory concepts to the exactly solvable many-particle problem [35].The behaviour of the kinetics𝒫B(t,r0) and 𝒫Bsp(t,r0) is given in Figure 3 for two values of the density parameter ξ. Note that relative deviations at small values of ξ (curves 1, 2) are not so great as for geminate reaction product accumulation. They are proportional to the density parameter. Increased distinctions beyond the limits of the binary interval (curves 3, 4) are due to the many-particle nature of the kinetics NB(t) (54).Figure 3
Exact kinetics𝒫B(t)—(1, 3) and superposition approximation 𝒫Bsp(t)—(2, 4) for different values of ξ: ξ=0.1—(1, 2); ξ=0.3—(3, 4). Other parameters: η=0.1, RAB/r0=0.95, τC/τB=0.1.


## 4.2. Diffusion Motion of Scavengers

For diffusion motion of scavengers comparison between binary formulae (39) and (40) shows that the bulk reaction kinetics N(t) differs from Ñ(t) by the exponential multiplier due to the difference between k and k̃, while the second power factor appears due to nonstationary parts of kinetics (39) and (40):  (59)N(t)Ñ(t)=exp⁡[-(ξ+3ξβ)τ][1-ξ-f(τ)ñ(τ)1-ξ-ñ(τ)],β=[C]k̃τB.
Here β and τ=t/τB are diffusion reaction rate with immobile B reactant and time in τB units, respectively. We introduce the notation for diffusion nonstationary part of kinetics (38), (39): (60)ñ(τ)=23ξβτπ.
Accordingly, the accumulation rate 𝒫̇Psp(t) of geminate reaction products predicted by superposition approximation (14) differs from the exact result in the same way. As is seen, for diffusion motions of scavengers these distinctions become more significant because of the presence of extra root term in the exponent of (59).To find the accumulation kinetics𝒫P(t) of geminate reaction products, use formulae (46), (22), and (39). We have:  (61)PP(t,r0)=(1-ξ)PP(m)(t,r0)-PP(n)(t,r0),
where 𝒫P(m) is the contribution of the Markovian (exponential) term of (39) calculated in the previous section:  (62)PP(m)(t,r0)=RABr0φ(η,β+1,τ),
where the function φ(η,x,τ) is defined in (34). The value 𝒫P(n)(t,r0) is determined by diffusion non-Markovian part of kinetics (39) and is equal to  (63)PP(n)(t,r0)=RABr03πξηβ∫0(β+1)τdxxexp⁡(-x-η(β+1)4x).
In the limit t→∞ the integral in (63) goes to the Bessel function. In this limit, using (35), we can rewrite (61) in an analytical form:  (64)PP(∞,r0)=RABr0[23π(1-ξ)e-η(β+1)-23πξηβK0(η(β+1))],
where K0(x) is the modified Bessel function of a zero order [59].The behaviour of the exact accumulation kinetics𝒫P(t) (61)–(63) and superposition approximation 𝒫Psp(t) (43)–(45) is shown in Figure 4.Figure 4
Exact accumulation kinetics𝒫P(t)—(2), its Markovian part 𝒫P(m)(t)—(4), and superposition approximation 𝒫Psp(t)—(1), its Markovian part 𝒫P(m)sp(t)—(3) for geminate reaction products in the presence of diffusing scavengers. ξ=0.1, η=3, RAB/r0=0.95, τB/τC=3DCτB/RBC2=0.3.The ratio of (43) and (61) at large τ and small ξ is recast as(65)PPsp(∞,r0)PP(∞,r0)=exp⁡[-ξ2ηβ+1(1+3βξ)]Θ,ξ≪1.
The first exponential factor in (65) is defined by the ratio between Markovian contributions (44) and (62), while Θ specifies the relation of non-Markovian parts of geminate product:  (66)Θ=1-ξ-PP(n)sp(∞,r0)/PP(m)sp(∞,r0)1-ξ-PP(n)(∞,r0)/PP(m)(∞,r0).
From (65) as well as from (52) of the previous subsection it follows that even at small ξ the result of superposition approximation can decrease essentially the ultimate probability of geminate recombination in the presence of scavengers in comparison with the exact solution if η is rather large.To calculate the decay kinetics𝒫B(t) of B reactant, we solve (47) by a conventional method of the Laplace transformation using expressions (27) and (39) for Ω̃(t,r0) and Ñ(t). This equation has a many-particle character. The correct binary approximation in the limit of ultrafast geminate reaction has to turn into the product of the escape probability and bulk kinetics: 𝒫B(t,r0)=Ω∞(r0)N(t) [35]. Taking into account this requirement, we obtain that(67)PB(t,r0)=(1-ξ)PB(m)(t,r0)-PB(n)(t,r0),
where the Markovian part of the kinetics can be represented as(68)PB(m)(t,r0)=e-[C]kt[1-RABr0φ(η,α,τ)],α=1-ξ-3ξβ.
This result coincides with the corresponding result of (53) at a small value of ξ except the parameter α. At large t expression (68) goes to the following:  (69)PB(m)(t,r0)≃e-[C]kt[1-RABr0e-αη]=e-[C]ktΩ∞eff(r0),t≫τB.
The non-Markovian part of the kinetics can be represented through quadrature of elementary functions.It is interesting to compare𝒫B and 𝒫Bsp (15), (24), (40) at a long time when B reactant leaves a geminate pair for the bulk. In this case the expressions for these kinetics are simplified, and we obtain the relation:  (70)PBsp(t,r0)PB(t,r0)≃Ω∞Ω∞eff(1-RABr0(e-αη-J(αη))3ξβ/α(1-ξ-3ξβ/2)Ω∞eff)-1,t≫τB,
where Ω∞ is the escape probability defined by (26), and J(x) is the following quadrature:  (71)J(x)=2π∫0∞dyyexp⁡(-y)erfc⁡(x4y).
It easy to see that in (70) the factor Ω∞/Ω∞eff is defined by the ratio of the Markovian parts of 𝒫Bsp and 𝒫B, while the multiplier between the brackets is the relation of the non-Markovian parts of the kinetics.


## 5. Analysis of Deviations

To describe the distinctions between the exact solution and superposition approximation, we introduce relative deviations of the quantities measured:(72)ΔPdot=(ṖP-ṖPsp)ṖP,ΔP=(PP-PPsp)PP,ΔB=(PB-PBsp)PB.
Taking into consideration that comparison is worthwhile in the region of small concentrations, in calculating relative deviations the exact many-particle quantities should be expanded in terms of a small parameter ξ. This imposes restrictions on the time interval of the comparison. In particular, the binary interval where the approximation linear in ξ is defined by the quadratic term smallness condition: [C]kt≪ξ-1. This criterion can be defined more precisely for our exactly solvable model. Note that, in fact, quadratic corrections are determined both by the expansion of the many-particle rate constant nonlinearly dependent on acceptor concentration and by the expansion of the nonstationary part of the kinetics. The smallness condition for corrections of the first type is determined by the inequality [C](k-k̃)t≪ξ-1, while the smallness of the second type corrections requires (as is seen from (39)) that the condition [C]k̃t≪ξ-1 be fulfilled. Thus the binary time interval where comparison of the theories will be made is given by the criterion:  (73)t≪tbin=ξ-1min⁡{([C](k-k̃))-1,([C]k̃)-1}.
The most considerable deviation exponentially increasing with time takes place for the first quantity in (72) whatever the relation between geminate and bulk reactions. Substituting the values of the rates 𝒫̇P and 𝒫̇Psp for the hopping and diffusion models of scavenger motion in ΔPdot, we have,  correspondingly,  (74)ΔPdot={1-e-ξτ,τ≪ξ-2,1-e-(ξ+3ξβ)τ[1-ξ-f(τ)ñ(τ)1-ξ-ñ(τ)],τ≪1ξ(ξ+3ξβ).
It is seen that exponential deviation 𝒫̇Psp from 𝒫̇P leads to the fact that even at half-decay times t≈τC/ξ of the kinetics Ñ(t) (31) the relative deviation ΔPdot may differ from unity just by several tens percent. In particular, for the hopping model of scavenger motion ΔPdot=0.63;0.86;0.95 for τB=τC, τB=τC/2, τB=τC/3, respectively. For diffusion model of scavenger motion these distinctions are even more essential. It should be noted that ΔPdot does not depend on initial distance in geminate pair r0. Thus, (74) are correct for any distribution of initial distance P(r0).The behaviour of the relative deviationΔPdot is given in Figure 5. It is evident that deviations are already very large even at the very beginning of a binary time interval both for diffusion (line 2) and for hopping (line 3) motion of scavengers. The second power factor of (49) slightly decreases the deviation in comparison with the Markovian exponential part of the kinetics designated in Figure 5 by line 1. Nevertheless the deviation for diffusion remains larger than that for the hopping motion of acceptors (the first line in (74)).Figure 5
Relative deviationΔPdot for diffusion (2) and for hopping motion of scavengers (3). Line (1) is the Markovian part ΔPdot(m) of deviation for the diffusion case. Parameters: ξ=0.1, τB/τC=3τBDC/RBC2=0.3.Consider the behaviour of the functionΔP that describes correlations arising in the course of geminate recombination product accumulation. Since time correlations accumulate with time, and superposition approximation takes no account of this process, the quantity ΔP increases monotonically with time and attains its steady-state value ΔPs rather rapidly. The behaviour of the relative deviation ΔP of accumulation kinetics (50) and (61) as compared to the results of superposition approximation (33) and (43) is given in Figure 6. As we can see, for diffusion model relative deviations ΔP (line 2) and its Markovian version ΔP(m)=1-𝒫P(m)sp/𝒫P(m) (line 1) grow up with τ and rapidly go to stationary values. The non-Markovian part of the kinetics decreases the relative deviation but nevertheless, as in Figure 5, it remain larger than for the model of the hopping motion of scavengers (line 3).Figure 6
Relative deviationΔP (2) and its Markovian part ΔP(m) (1) for diffusing scavengers in comparison with ΔP for hopping scavengers (3). Parameters: ξ=0.1, η=10, RAB/r0=0.95, τB/τC=3τBDC/RBC2=0.3.Using the expansion of the functions𝒫Psp(∞,r0)/𝒫P(∞,r0) from (52) and (65) for hopping and diffusion models of scavenger motion at small values of the parameter ξ, we have the relation for stationary relative deviations:  (75)ΔPs={1-exp⁡(-ξ2ην+1),hopping,1-exp⁡(-ξ2ηβ+1(1+3βξ))Θ,diffusion,
where the function Θ is given by relation (66).It is easy to see that at the start at the contact (r0=RAB,η=0) the function ΔPs is equal to zero. Then, as η increases, ΔPs rises monotonically gradually approaching its asymptotic value equal to unity. However, it should be taken into account that the use of (75) is justified only in the case where the attainment of steady-state values occurs in the binary interval. This imposes definite restrictions on the values of the parameter η essential in the analysis of ΔPs. Indeed, the use of the binary kinetics in (14) calls for the smallness of the characteristic decay time τg (23) of a geminate pair as compared to the binary interval (τg≪tbin). As a result, at η≫1, when τg=3τ̃gτB, an insignificant upper bound is set on the parameter η. In the hopping case this restriction is defined by the inequality η≪ξ-2. However, this restriction is no barrier to the realization of considerable relative deviations ΔPs possible on condition that the exponent index in (75) is greater than unity. For the hopping model of scavengers motion this condition is defined by the inequality:  (76)ξ2ην+1=ξ6τgτB11+ξτB/τC>1.
For acceptors at rest (τC→∞) condition (76) is realized at  (77)ξτgτB>1.
Physically this means that geminate pair lifetime is long enough for the bulk reaction (proceeding at the rate ξ/τB) to reach considerable depth. Time correlations in this case are maximum. When C reactant starts moving, correlations decrease. In the region ξτB≪τC≪τB, where the bulk reaction rate is ξ/τC (since this rate can be calculated neglecting B reactant motion), but the bulk reaction depth for the time τB(ν≪1) is small, correlations become essential under more rigid restrictions:  (78)ξτgτC>τBτC≫1.
As the motion of C reactants becomes more active (τC≪ξτB), the region of essential deviations is reached under even more strict conditions:  (79)ξτgτC>τBτCξτBτC≫τBτC≫1.
For diffusion motion of scavengers the condition of large relative deviations is as follows:  (80)ξ2ηβ+1(1+3βξ)=ξ6τgτB(1+3τB/τC)1+ξτB/τC>1,
where the mean residence time τC=RBC2/(3DC) of diffusing scavengers in the reaction zone is introduced. Obviously, with τC→∞ inequality (80), just as (76), holds subject to condition (77). However, the situation reverses when diffusion of scavengers becomes so rapid that the mean residence time τC lies in the interval ξτB≪τC≪τB. Then the second term between the brackets in (80) becomes decisive, and this equality may be written as  (81)ξτgτC>τBτC≫1.
It is seen that restriction (81) is much less rigid than (77). However, in this range of parameters the binary interval boundary for diffusion motion of scavengers is already determined by the smallness condition of non-Markovian corrections (the second term in (73)). Thus a more rigid upper bound on the value of the parameter η is set: η≪(ξβ)-1. This is equivalent to the condition ξτg/τC≪ξ-1. This inequality does not prevent the fulfillment of the realization condition of large relative deviations (80), though it sets the upper bound on this value. Thus condition (81) takes the form:  (82)ξ-1≫ξτgτC>τBτC≫1.
Finally, when the motion of scavengers becomes dominant (τC≪ξτB), condition (80) is met in the following restricted interval:  (83)ξ-1≫ξτgτC>ξτBτC≫1.
Obviously, large relative deviations can be realized only on condition that the right and the left boundaries do not overlap: τC≫ξ3/2τB. With progressively decreasing τC the upper and the lower boundaries in inequality (83) approach one another thus eliminating the possibility of the realization of considerable relative deviations in the limit of dominating motion of scavengers. In the limit of infinitely high mobility of acceptors (τC→0) the difference between N(t) and Ñ(t) vanishes, so the superposition approximation coincides with the exact solution, and the relative deviation tends to zero.Thus conditions (77)–(79) and (81)–(83) show that correlations decrease with increasing intensity of C reactants motion. It should be noted that in the region where deviations are small (reversal of inequalities (76) and (80)) they are still linear in ξ, that is, measurable in the framework of the binary approach.  The relative deviation ΔPs as a function of parameter η in the case of dominate electron motion (τC≫τB) is shown in Figure 7. It is readily seen that taking into account the non-Markovian part of deviation defined by the factor Θ in  (52) decreases noticeably the values of ΔPs for diffusion model (line 2) in comparison with pure Markovian part (line 1) including only the first exponential part of relation (75). The corresponding deviation for the hopping acceptors motion (line 3) lies everywhere below.Figure 7
Relative deviationΔPs defined by (75) for hopping (3) and for diffusing (2) scavengers versus η=τ̃g/τB. Line (1) is the Markovian part of (2). Parameters: ξ=0.1, τB/τC=3τBDC/RBC2=0.3.Now we pass to the analysis of the functionΔB. Since the superposition approximation retards the decay kinetics of B reactant as compared to the exact solution, this function is everywhere negative, except the initial instant of time. For the model of the hopping motion of scavengers in the binary interval it is of a simplest form:  (84)ΔB=1-[1-(RAB/r0)φ(η,1,τ)][1-(RAB/r0)φ(η,1-ξ,τ)].
We will not give similar expression for the model of diffusion motion of scavengers due to its cumbersome form. As shown below it behaves similarly to (84). One can readily see that function (84), though always negative, attains rapidly its steady-state value ΔBs defined by asymptotic (35) of the function φ(η,x,τ). In the range of small values of ξ the expression for ΔBs can be brought to the form  (85)ΔBs=-RABr0e-η(eξη/2-1)1-(RAB/r0)e-η(1-ξ/2).
Examine the behaviour of the function ΔBs (85). If a geminate pair is formed not at the contact (r0>RAB), then under variation of the parameter η from 0 to ∞ this function begins with zero, passes through the minimum in the area η<1, and then smoothly tends to zero again with increasing  η. The closer the starting point r0 to the reaction zone, the sharper and deeper the minimum, and the greater its shift to the area of small values of  η. In the region of the extremum η<1 expression (85) for ΔBs may be simplified:  (86)ΔBs=-ξ2RABr0ηeη-(RAB/r0)(1+ξη/2).
Function (86) has the minimum at the point ηm determined by a transcendental equation:  (87)RABr0e-ηm=1-ηm.
This equation has the only root lying in the interval 0≤ηm<1. Substituting (87) in (86), we have with the accuracy up to the main order in ξ:  (88)min⁡ΔBs≃-ξ2(1-ηm).
The most essential deviations |ΔBs|≃ξ/2 arise at ηm≪1. Examining the function ΔBs in this interval, we easily obtain the analytical expression for the minimum:  (89)ηm≈2(1-RABr0)≪1.
As follows from (88), (89), the greatest deviation is realized in the case where a geminate pair is formed immediately near the contact. The definition of the parameter η makes it possible to estimate the thickness of the layer the initial start from which leads to maximum time correlations:  (90)r0-RAB≃2DAτBRAB≪RAB.
Result (90) may be rewritten as the ratio of times in view of the fact that the residence time of a diffusing particle in a thin spherical layer of the thickness δ of the radius RAB is δRAB/DA [22]. Really, expression (90) can be recast as  (91)(r0-RAB)RAB2DA≃τB.
It follows from the above relation (91) that maximum time correlations are realized when the residence time of diffusing A reactant in a thin layer (r0-RAB)/2 is equal to the time τB between B reactant jumps. Note that in this case relative deviations (88) are not as large as they are in the analysis of geminate reaction product accumulation, however, they are proportional to the density parameter, and, therefore, can be measured in the frame of the binary approximation.Figure8 shows the dependence of the relative deviation ΔBs versus parameter η. Underestimating the rate of geminate reaction, the superposition approximation thus overestimates the probability of B reactant to survive and escape into the bulk. Apparently, ΔBs is always negative. The largest relative deviations are observed for small η when the contribution of geminate reaction to total kinetics 𝒫B(t) is most substantial.Figure 8
Relative deviationΔBs for the diffusing scavengers (2) and its Markovian part (1) versus η=τ̃g/τB. Line 3 is the ΔBs for the model of hopping scavengers. Parameters: ξ=0.1,τB/τC=3τBDC/RBC2=0.3.


## 6. Conclusion

In this brief review the “scavenger problem” (1), (2) is analyzed with two examples of many-particle model of reacting particles that admits exact analytical solution for any concentrations of bulk acceptors. The derived exact solutions are compared with the results of approximate theories [1–11] widely used to describe experimental data. They are shown to coincide only in the limiting case where the mobility of excess electron (geminate B partner) involved in the competing bulk reaction is equal to zero. In this case time correlations between bulk and geminate reactions are absent at any concentrations of acceptors. Thus the accumulation rate 𝒫̇P of geminate reaction products and the decay kinetics 𝒫B of B reactant are expressed as the product of bulk kinetics (31), (38) and, respectively, geminate reaction rate (22) and geminate kinetics (27):  (92)ṖP(t)=Ñ(t)K̃g(t),PB(t)=Ñ(t)Ω̃(t),1τB=0.
As soon as B reactant starts moving, time correlations appear in the system. As a consequence, the rate in a geminate channel ceases to be proportional to the bulk reaction kinetics, as it was in the case of independent course of reactions (92). When B reactant becomes mobile, the bulk channel is modified accordingly to formulae (49), (59), but the rate in the geminate channel is equal to the product of pure geminate reaction rate Kg(t) (22) and the kinetics Ñ(t) formed solely during the residence time of B particle in a geminate pair: 𝒫̇P(t)=Ñ(t)Kg(t). As a result, the reaction rate in the geminate channel differs strongly from the prediction following from the superposition approximation: 𝒫̇Psp(t)=Kg(t)N(t). This, in turn, leads to the fact that the superposition approximation substantially underestimates the yield of geminate reaction products in the cases where the bulk reaction proceeds deeply during geminate pair lifetime. These distinctions are more pronounced in the case of diffusion motion of acceptors rather than in that of the hopping one.Another consequence of time correlations between bulk and geminate channels of reactions is that the survival probability𝒫B(t) of B reactant cannot be represented as a product of the bulk kinetics and the kinetics of geminate reaction proceeding independently in some effective pair. If considered in such a way, formula (53) gives rise to a physical contradiction: taking account of additional decay channel in effective acceptor surroundings accelerates the reaction in effective geminate pair, though obviously it can do nothing but slow it down.Evidently, large distinctions in the theoretical treatment of the measured quantities are to affect experimental data interpretation. In a classical variant [21] the experiment is run as follows: the dependence of the scavenging yield F on the concentration [C] of bulk acceptors is measured with the subsequent determination of the geminate pair kinetics Ω(t) from the measured concentration dependence F(α[C]). Here α is the experimentally measured constant depending on the initial distance distribution in geminate pair and acceptor and solvent properties. Commonly, conventional theory formula (8) is employed assuming that the scavenging yield is the Laplace transform of geminate reaction kinetics:  (93)PDct(t⟶∞)=[C]k∫0∞dtΩ(t)e-[C]kt≈F(α[C]).
The inverse Laplace transformation of (93) with [C]k as transformation variable gives the geminate kinetics Ω(t). Its characteristic decay time is determined from the ratio τg~α/k, with the rate constant value k of the bulk reaction taken from independent experiments on pure bulk reaction. However, it follows from the exact solution (46) that even if “instantaneous quenching” effect is neglected, formula (93) can be reproduced but with an essential distinction: instead of the bulk reaction rate [C]k (equal to ξ(1/τC+1/τB) for hopping and ξ(1/τC+1/τB+3/(τCτB)) for diffusing scavengers), the role of the Laplace variable will be played by the reaction rate with immobile B reactant: [C]k̃=ξ/τC that can be essentially less than [C]k. Thus the value τg found by formula (93) is k/k̃ times less than it actually is. For instance, in liquid hydrocarbons where the mobility of an excess electron substantially exceeds that of acceptors, the constant k is determined mainly by electron mobility. Indeed, using both the lifetime of the excess electron in the localized state in cyclohexane τB~10-10 s [43] and typical values RBC=10 Å and DC=10-5 cm2/s, we get k/k̃≈7.5. Note that this value is almost twice greater than the corresponding ratio k/k̃ for the hopping motion of acceptors. Since for our model, according to (23), τg≈3τ̃gτB=3(r0-RAB)τB/DA, the value of the initial distance in the geminate pair decreases similarly. In real situation when the initial distance in geminate pair r0 is not fixed but is distributed according to P(r0) the characteristic time of geminate reaction is proportional to the width L on the distribution: τg≈LτB/DA. So one can considerably underestimate the width of required distribution.Another part of experimental works is concerned with the rate constant of the bulk reaction of a geminate partner decay on acceptors. In particular, paper [48] suggests an ingenious method of excluding pure geminate kinetics from geminate reaction channel based on formula (14) of the superposition decoupling: 𝒫̇Psp(t)=Kg(t)N(t). It is seen from this formula that the bulk kinetics N(t) may be obtained by dividing the experimentally measured fluorescence rate 𝒫̇Psp(t) of geminate reaction products in the presence of acceptors by the fluorescence rate Kg(t) of pure geminate reaction products. However, as follows from exact solution (46), the kinetics Ñ(t) thus obtained is not the bulk kinetics and can differ considerably from it. Since the distinction of the exact model and superposition approximation is determined by the difference between the bulk kinetics N(t) and kinetics Ñ(t) the relative deviations ΔPdot does not depend of the distribution of initial distances in geminate pair.In conclusion it should be noted that the exactly solvable model considered above refers to the hopping motion of excess electrons with the mean length of a jump large enough for an electron to leave geminate pair in one jump and change its surroundings in the process of random walks in the bulk in an uncorrelated way. Such a model of excess electron mobility appears to be realizable in such hydrocarbons as isooctane, neopentane, where the mean free path of an excess electron is about 60–200 Å [23, 43, 45]. On the other hand, due to a large value of Onsager radius in such systems (~300 Å) it is necessary for an electron to execute several jumps to leave geminate pair and escape into the bulk, so an addition correlation caused by Coulomb attraction should be taken into account. That is why, strictly speaking, our many-particle exactly solvable model can first of all be used to describe photodetachment processes of a high mobility excess electron generation in solution.Also it should be noted that the exactly solvable model was used to establish the form of kinetic equations [60, 61]. It is consistent with the concepts of general kinetic theory [62], according to which in the general case these equations are inhomogeneous integrodifferential ones. This allows one to elucidate the physical meaning of the involved kinetic coefficients, and thus to use the exactly solvable model for testing the kinetic equations obtained on the basis of consistent many-particle approach for a more general character of reactant mobility.



---

*Source: 451670-2011-04-10.xml*