Abstract

A geminate reaction between 𝐴 and 𝐡 reactants affected by the bulk reaction between 𝐡 and scavengers 𝐢 is discussed. The exact solution of the problem obtained recently for randomly walking 𝐡 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.

1. Introduction

More than fifty years ago Monchick [1] formulated the problem of the influence of scavengers with a uniform concentration [𝐢] 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, 𝐴 and 𝐡 radicals or ion-radicals produced from a β€œparent" molecule form geminate (spatially correlated) pairs [12–17]. The subsequent reaction 𝐴+π΅βŸΆπ‘ƒ(1) 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 distributed 𝐢 acceptors of one (say 𝐡) type of radicals to the solution. Usually 𝐢 reactants are in excess (as compared to 𝐡𝑠), thus their concentration remains almost unchanged 𝐡+𝐢⟢𝐢+𝐷.(2) Note that the excess of 𝐢 reactants is no barrier to taking into account solely pair encounters of 𝐡 with 𝐢 reactants in considering small concentrations of scavengers [𝐢]. The problem consists in calculating the decay kinetics of 𝐴 or 𝐡 reactants (or the formation kinetics of 𝑃 or 𝐷 products of geminate (1) or bulk (2) reactions, resp.), and the dependence of 𝑃 or 𝐷 product yield on 𝐢 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 [𝐢]π‘˜ into the Smoluchowski equation defining the evolution of isolated geminate pair: πœ•π‘‘πΊξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ,0=π·π΄π΅βˆ‡βƒ—π‘Ÿξ‚Έβˆ‡βƒ—π‘ŸπΊπΊ+π‘˜π΅π‘‡βˆ‡βƒ—π‘ŸΞ¦ξ‚Ήβˆ’[𝐢]π‘˜πΊ,(3) where 𝐷𝐴𝐡=𝐷𝐴+𝐷𝐡 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 π‘˜ is a bimolecular rate constant of bulk reaction (2) depending on mutual mobility of 𝐡 and 𝐢 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 πΊξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ[𝐢],0=exp(βˆ’π‘˜π‘‘)πΊπ΄π΅ξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ,0(4) converts (3) to the Smoluchowski equation in the absence of scavengers for 𝐺𝐴𝐡. 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 𝒫𝐷(π‘‘β†’βˆž) (that is defined as the probability that 𝐡 radical react with a scavenger before recombination with 𝐴 radical), can be related to the survival probability of a geminate pair as follows: 𝒫𝐷[𝐢]π‘˜ξ€œ(π‘‘βŸΆβˆž)β‰ƒβˆž0𝑑𝑑Ω(𝑑)π‘’βˆ’[𝐢]π‘˜π‘‘.(5) 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 Ξ©(𝑑). 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 𝐴 and 𝐢 reactants mobility and the Coulomb interaction in a geminate pair [35]. The most essential assumption of the model is that 𝐡 particle of a geminate pair (excess electron) reacting with scavengers 𝐢 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 [𝐢]. 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 Section 3 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 𝐡 and 𝐢 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 𝐡 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, 𝒫𝐡(𝑑) (survival probability of 𝐡 reactant) and 𝒫𝑃(𝑑) (accumulation kinetics of geminate reaction products) are quantities measured experimentally. The survival probability 𝒫𝐴(𝑑) of 𝐴 reactant, the kinetics of scavenging, 𝒫𝐷(𝑑), can be found by using the normalization conditions [34, 35]: 𝒫𝐴(𝑑)+𝒫𝑃(𝑑)=1,𝒫𝐡(𝑑)+𝒫𝐷(𝑑)+𝒫𝑃(𝑑)=1.(6) The last relation is simplified in the stationary limit π‘‘β†’βˆž, thus the experimentally observed ultimate probability of recombination 𝒫𝑃(π‘‘β†’βˆž) and total scavenging yield 𝒫𝐷(π‘‘β†’βˆž) obey the relation: 𝒫𝐷(π‘‘βŸΆβˆž)+𝒫𝑃(π‘‘βŸΆβˆž)=1.(7) 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 ̇𝒫𝑃(𝑑) [48].

Conventional theory [1–10] enters constant decay rate on acceptors [𝐢]π‘˜ 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 [𝐢] and stationary rate constant π‘˜ of 𝐡 particles decay on these acceptors. Thus we have the recipe for the calculation of experimentally measured quantities: 𝒫𝑃(𝑑)≃𝒫𝑃𝑐𝑑(ξ€œπ‘‘)≑𝑑0π‘‘πœπΎπ‘”([𝐢]π’«πœ)exp(βˆ’π‘˜πœ),(8)𝐡(𝑑)≃𝒫𝐡𝑐𝑑[𝐢](𝑑)≑Ω(𝑑)exp(βˆ’π‘˜π‘‘).(9) As usual, geminate reaction rate in the absence of acceptors 𝐾𝑔(𝑑) is expressed in terms of the elementary rate 𝑀𝐴𝐡(βƒ—π‘Ÿ):𝐾𝑔(ξ€œπ‘‘)=π‘‘βƒ—π‘Ÿξ…žξ€œπ‘‘βƒ—π‘Ÿπ‘€π΄π΅ξ€·ξ€ΈπΊβƒ—π‘Ÿπ΄π΅ξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έπ‘ƒξ€·,0βƒ—π‘Ÿξ…žξ€Έ,(10) where 𝑃(βƒ—π‘Ÿξ…ž) is the initial distribution of reactants in geminate pairs, and 𝐺𝐴𝐡(βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…ž,0) is the probability density to find 𝐴 and 𝐡 reactants at the relative distance βƒ—π‘Ÿ at the instant of time 𝑑 if at the initial moment 𝑑=0 they were at the distance βƒ—π‘Ÿξ…ž (the Green function). It obeys the equation πœ•π‘‘πΊπ΄π΅ξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ=ℒ,0π΄π΅βˆ’π‘€π΄π΅ξ€·ξ€Έξ‚„πΊβƒ—π‘Ÿπ΄π΅ξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ,0(11) with the initial condition πΊπ΄π΅ξ€·βƒ—π‘Ÿ,0βˆ£βƒ—π‘Ÿξ…žξ€Έξ€·,0=π›Ώβƒ—π‘Ÿβˆ’βƒ—π‘Ÿξ…žξ€Έ.(12) In (11), ℒ𝐴𝐡 is the integral operator defining Markovian random walks of 𝐴 and 𝐡 in relative coordinates of the pair 𝐴𝐡. The force interaction is included.

The survival probability of geminate pairs in the absence of competing bulk reaction is defined by the expression: ξ€œΞ©(𝑑)=1βˆ’π‘‘0π‘‘πœπΎπ‘”(𝜏).(13) 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: 𝒫𝑃(𝑑)≃𝒫𝑃𝑠𝑝(ξ€œπ‘‘)≑𝑑0π‘‘πœπΎπ‘”(𝜏)𝑁(𝜏),(14) while 𝐡 reactant decay kinetics is expressed as the product of bulk and geminate kinetics: 𝒫𝐡(𝑑)≃𝒫𝐡𝑠𝑝(𝑑)≑𝑁(𝑑)Ξ©(𝑑).(15) 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]: ξ‚»βˆ’[𝐢]ξ€œπ‘(𝑑)=exp𝑑0ξ‚Όπ‘‘πœπΎ(𝜏).(16) Note that the expression in exponent (16) is linear in concentration of acceptors [𝐢], since it allows solely for binary collisions of reactants in the bulk. In the general case of arbitrary mobilities of 𝐡 and 𝐢 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 𝐡 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 ξ€œπΎ(𝑑)=π‘‘βƒ—π‘Ÿξ…žξ€œπ‘‘βƒ—π‘Ÿπ‘€π΅πΆξ€·ξ€ΈπΊβƒ—π‘Ÿπ΅πΆξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ,0.(17) Stationary rate constant of the decay on acceptors π‘˜ appearing in formulae (8) and (9) is determined from (17): π‘˜=limπ‘‘β†’βˆžπΎ(𝑑).(18) The Green function 𝐺𝐡𝐢(βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…ž,0) of the pair 𝐡𝐢 is defined by the equation πœ•π‘‘πΊπ΅πΆξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ=ℒ,0π΅πΆβˆ’π‘€π΅πΆξ€·ξ€Έξ‚„πΊβƒ—π‘Ÿπ΅πΆξ€·βƒ—π‘Ÿ,π‘‘βˆ£βƒ—π‘Ÿξ…žξ€Έ,0(19) with the initial condition of type (12). As in (11), ℒ𝐡𝐢 is the integral operator describing translational motion of reactants in relative coordinates of the reaction pair, and bulk reaction between 𝐡 and 𝐢 reactants is specified by the elementary rate 𝑀𝐡𝐢(βƒ—π‘Ÿ) 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 between 𝐴 and 𝐡 reactants. Second, the situation of highly mobile 𝐡 reactant and relatively slow 𝐢 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 radius 𝑅𝐴𝐡 [53]. The β€œblack” sphere model excludes recontacts of reactants and, consequently, related time correlations in the system.

Let us take that uncharged 𝐴 reactant moves by continual diffusion, while the motion of 𝐡 reactant (excess electron) is described by the model of infinite jumps with the mean frequency πœπ΅βˆ’1 (the so-called two-scale migration [54]): ℒ𝐴𝐡=π·π΄βˆ‡2βƒ—π‘Ÿβˆ’πœπ΅βˆ’1𝐼.(20) Here 𝐼 is a unit operator.

Starting with the initial distance π‘Ÿ0π‘ƒξ€·ξ€Έπ›Ώξ€·βƒ—π‘Ÿ=𝑃(π‘Ÿ)=π‘Ÿβˆ’π‘Ÿ0ξ€Έ4πœ‹π‘Ÿπ‘Ÿ0,(21) we easily obtain that at π‘Ÿ0>𝑅𝐴𝐡 geminate reaction rate in the absence of bulk acceptors is defined as𝐾𝑔𝑑,π‘Ÿ0ξ€Έ=𝐾𝑔𝑑,π‘Ÿ0ξ€Έπ‘’βˆ’π‘‘/𝜏𝐡,𝐾𝑔𝑑,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0ξƒŽΜƒπœπ‘”4πœ‹π‘‘3ξ‚΅βˆ’expΜƒπœπ‘”ξ‚Ά,4𝑑(22) where Μƒπœπ‘”=(π‘Ÿ0βˆ’π‘…π΄π΅)2/𝐷𝐴 is the characteristic decay time of the geminate pair, and 𝐾𝑔 is the geminate reaction rate for the case of immobile 𝐡 reactant. Further, quantities related to the problem with immobile 𝐡 particle will be marked by a tilde.

By analogy with Μƒπœπ‘” we introduce the characteristic decay time of the geminate pair πœπ‘” for the case of two-scale migration defining it as the characteristic decay time of the function 𝐾𝑔(𝑑,π‘Ÿ0). Obviously, it is proportional to the time moment π‘‘π‘š of the maximum of this function πœπ‘”=πœ†π‘‘π‘š. The proportionality coefficient πœ† is determined from the coincidence condition of πœπ‘” and Μƒπœπ‘” at πœπ΅β†’βˆž. It is easily seen that πœ†=6. So for πœπ‘”, πœπ‘”=32πœπ΅βŽ›βŽœβŽœβŽξƒŽ9+4Μƒπœπ‘”πœπ΅βŽžβŽŸβŽŸβŽ ξ”βˆ’3β‰ˆ3Μƒπœπ‘”πœπ΅,𝜏𝐡β‰ͺΜƒπœπ‘”.(23) The use of (22) in (13) with the subsequent integration gives the expression for the survival probability: Ω𝑑,π‘Ÿ0𝑅=1βˆ’π΄π΅2π‘Ÿ0ξƒ¬π‘’βˆ’βˆšπœ‚ξƒ©ξ‚™erfcΜƒπœπ‘”βˆ’ξ‚™4π‘‘π‘‘πœπ΅ξƒͺ+π‘’βˆšπœ‚ξƒ©ξ‚™erfcΜƒπœπ‘”+ξ‚™4π‘‘π‘‘πœπ΅,ξƒͺξƒ­(24) where the introduced parameter πœ‚=Μƒπœπ‘”πœπ΅=ξ€·π‘Ÿ0βˆ’π‘…π΄π΅ξ€Έ2𝐷𝐴𝜏𝐡(25) is defined by the ratio of Μƒπœπ‘” to 𝐡 reactant residence time in the pair 𝜏𝐡. Since geminate reaction is arrested immediately after the first jump of 𝐡 reactant, the depth of the reaction course depends noticeably on the parameter πœ‚: Ξ©ξ€·βˆž,π‘Ÿ0ξ€Έβ‰‘Ξ©βˆžξ€·π‘Ÿ0𝑅=1βˆ’π΄π΅π‘Ÿ0ξ‚€βˆ’βˆšexpπœ‚ξ‚.(26) Only in the limit πœ‚β†’0 it is equal to the probability of the diffusion entry of reactants into the bulk: 1βˆ’π‘…π΄π΅/π‘Ÿ0. In this limiting case of immobile 𝐡(πœπ΅β†’βˆž) expression (24) is easily seen to turn into the well-known diffusion result [21]: limπœπ΅β†’βˆžΞ©ξ€·π‘‘,π‘Ÿ0≑Ω𝑑,π‘Ÿ0𝑅=1βˆ’π΄π΅π‘Ÿ0erfcΜƒπœπ‘”ξƒͺ4𝑑.(27) 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 |π‘Ÿπ‘|≀10 Å. As follows from [55, 56], in this case it will suffice to perform the following substitution: Μƒπœπ‘”βŸΆΜƒπœπ‘π‘”=π‘Ÿ2𝑐4π·π΄ξ‚Έξ‚΅π‘Ÿcoth𝑐2π‘Ÿ0ξ‚Άξ‚΅π‘Ÿβˆ’coth𝑐2𝑅𝐴𝐡2,π‘…π΄π΅π‘Ÿ0βŸΆπ‘Ÿeο¬€ξ€·π‘…π΄π΅ξ€Έπ‘Ÿeο¬€ξ€·π‘Ÿ0ξ€Έ,π‘Ÿeff(||π‘Ÿπ‘Ÿ)=𝑐||ξ€·βˆ’||π‘Ÿ1βˆ’exp𝑐||ξ€Έ,/π‘Ÿ(28) because in highly polar solvents the survival probability of charged geminate pairs Ω𝑐(𝑑,π‘Ÿ0) is kinetically identical to the corresponding value Ω(𝑑,π‘Ÿ0) without the Coulomb interaction.

3.2. Hopping Motion of Scavengers

When calculating the bulk reaction kinetics, in this subsection we shall take that 𝐢 reactants (just as 𝐡) move by jumps of infinite length: ξβ„’π΅πΆξ€·πœ=βˆ’π΅βˆ’1+πœπΆβˆ’1𝐼.(29) 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 radius 𝑅𝐡𝐢 insight of wich the reaction proceeds instantaneously is used for the calculation of the bulk kinetics 𝑁(𝑑). 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 ℒ𝐡𝐢 in form (29) yield ξ€Ίξ€·πœπ‘(𝑑)=expβˆ’πœ‰βˆ’πœ‰π΅βˆ’1+πœπΆβˆ’1𝑑,(30) where the dimensionless parameter πœ‰=[𝐢]𝑉 (the density parameter) appears. It is equal to the average number of acceptors in the reaction zone of the volume 𝑉=(4πœ‹/3)𝑅3𝐡𝐢. Here we neglect the intrinsic volume of reactants and refer to long-range reactions.

Kinetics (30) consists of two cofactors. The first co-factor, π‘’βˆ’πœ‰, 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 πœπ΅βˆ’1+πœπΆβˆ’1 and the average number of acceptors πœ‰ entering the reaction zone. Note that in the limit of immobile 𝐡 reactant (πœπ΅β†’βˆž)𝑁(𝑑) tends to the exact many-particle kinetics of the target problem 𝑁(𝑑) that for the hopping motion of scavengers is𝑑𝑁(𝑑)=expβˆ’πœ‰βˆ’πœ‰πœπΆξ‚Ά.(31) Passing to experimentally measured quantities of the β€œscavenger problem”, note that for both approximate theories the expression for 𝐡 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) 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: 𝒫𝑃𝑠𝑝𝑑,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0π‘’βˆ’πœ‰πœ‘(πœ‚,𝛾+1,𝜏),(33) where 𝛾=πœ‰(𝜏𝐡/𝜏𝐢+1) and 𝜏=𝑑/𝜏𝐡 are bulk reaction rate and time in 𝜏𝐡 units, respectively.

Formula (33) involves the function πœ‘(πœ‚,π‘₯,𝜏): 1πœ‘(πœ‚,π‘₯,𝜏)=2ξƒ―π‘’βˆ’βˆšπœ‚π‘₯erfcπœ‚βˆ’βˆš4𝜏ξƒͺπ‘₯𝜏+π‘’βˆšπœ‚π‘₯erfcπœ‚+√4𝜏,π‘₯𝜏ξƒͺξƒ°(34) 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: ξ‚€βˆ’βˆšπœ‘(πœ‚,π‘₯,𝜏⟢∞)=expξ‚πœ‚π‘₯.(35) This gives the expression for the ultimate probability of geminate reaction in the presence of bulk acceptors: π’«π‘ƒπ‘ π‘ξ€·π‘‘βŸΆβˆž,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0π‘’βˆ’πœ‰ξ‚€βˆ’βˆšexpξ‚πœ‚(𝛾+1).(36) 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 𝐡 reactant 𝒫𝐡𝑠𝑝 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 that 𝐢 reactants move by continual diffusion with the diffusion coefficient 𝐷𝐢, while 𝐡 reactant (excess electron), as before, moves by infinite jumps with the mean time 𝜏𝐡 between them: ℒ𝐡𝐢=π·πΆβˆ‡2βƒ—π‘Ÿβˆ’πœπ΅βˆ’1𝐼.(37) During this time interval 𝜏𝐡 between jumps the bulk reaction develops due to acceptors diffusion to immobile electron (the β€œtarget problem” kinetics 𝑁(𝑑)), 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 𝑁(𝑑) is formed.

The model of a β€œblack” ball of the radius 𝑅𝐡𝐢 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 πœπ΅βˆ’1=0 the β€œtarget problem” kinetics 𝑁(𝑑) proceeds as the well-known Smoluchowski kinetics [53]: ξ‚ξƒ¬βˆ’[𝐢]̃𝑁(𝑑)=exp𝑉+π‘˜π‘‘+23πœ‹π‘‰Μƒ,Μƒπ‘˜π‘‘ξƒͺξƒ­π‘˜=4πœ‹π‘…π΅πΆπ·πΆ,(38) where 𝑉=(4πœ‹/3)𝑅3𝐡𝐢 is the reaction zone volume. The density parameter πœ‰=[𝐢]𝑉 equal to the average number of acceptors in the reaction zone of the volume 𝑉 can be introduced as in hopping kinetics (30).

The first co-factor of (38), π‘’βˆ’πœ‰, is the same as in (30), the second co-factor describes the exponential decay with the diffusion rate constant Μƒπ‘˜, 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 [𝐢]: 𝑁(𝑑)β‰ˆπ‘’βˆ’[𝐢]Μƒπ‘˜π‘‘ξƒ¬ξ‚™1βˆ’πœ‰βˆ’23πœ‹πœ‰[𝐢]Μƒξƒ­π‘˜π‘‘.(39) The validity condition of formula (39) defines the binary time interval Μƒ[𝐢]π‘˜π‘‘β‰ͺπœ‰βˆ’1 in which time-dependent observable values will be considered.

The kinetics of the bulk reaction 𝑁(𝑑) (16) can be obtained by solving (19) with operator (37) in the β€œblack” sphere limit. In the binary approximation it takes the form𝑁(𝑑)β‰ˆπ‘’βˆ’[𝐢]π‘˜π‘‘ξƒ¬ξ‚΅π‘‘1βˆ’πœ‰βˆ’2π‘“πœπ΅ξ‚Άξ‚™3πœ‹πœ‰[𝐢]Μƒξƒ­π‘˜π‘‘,(40) where, in contrast to the case of diffusion mobilities of both reactants, stationary rate constant of the bulk two-scale reaction isΜƒπ‘‰π‘˜=π‘˜+𝜏𝐡+ξ‚™3Μƒπ‘˜π‘‰πœπ΅,(41) and 𝑓(𝜏) is a function decreasing from unity to zero: 1𝑓(𝜏)=2ξ‚Έπ‘’βˆ’πœ+ξ‚΅ξ‚™πœ‹+√4πœξ‚Άξ‚€βˆšπœ‹πœerfπœξ‚βˆ’βˆšξ‚Ήπœ‹πœ.(42) At the time π‘‘β‰«πœπ΅ it goes to its asymptotic βˆšπ‘“(𝜏)β‰ˆ(πœ‹/𝜏)/4, and the nonstationary part of bulk kinetics becomes time independent.

Stationary rate constant (41) is equal to the sum of diffusion rate constant Μƒπ‘˜ (38), hopping rate constant 𝑉/𝜏𝐡 [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 𝑓(𝑑/𝜏𝐡).

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𝒫𝑃𝑠𝑝𝑑,π‘Ÿ0ξ€Έ=(1βˆ’πœ‰)𝒫𝑃(π‘š)𝑠𝑝𝑑,π‘Ÿ0ξ€Έβˆ’π’«π‘ƒ(𝑛)𝑠𝑝𝑑,π‘Ÿ0ξ€Έ,(43) where the Markovian part of the kinetics corresponding to pure exponential term of (40) is: 𝒫𝑃(π‘š)𝑠𝑝𝑑,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0[𝐢]πœ‘(πœ‚,𝛾+1,𝜏),𝛾=π‘˜πœπ΅.(44) As in the previous subsection, the introduced parameter 𝛾 is the bulk reaction rate in 𝜏𝐡 units. Using (41), it can be rewritten in our case of diffusion motion of acceptor as βˆšπ›Ύ=𝛽+πœ‰+3πœ‰π›½, where ̃𝛽=[𝐢]π‘˜πœπ΅.

The non-Markovian part of the kinetics is defined by the integral 𝒫𝑃(𝑛)𝑠𝑝𝑑,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0√3πœ‹βˆšξ€œπœ‰πœ‚π›½0(𝛾+1)πœπ‘‘π‘₯π‘₯ξ‚΅Γ—expβˆ’π‘₯βˆ’πœ‚(𝛾+1)𝑓π‘₯4π‘₯ξ‚Ά.𝛾+1(45) This integral is not taken analytically even in the stationary limit πœβ†’βˆž.

In accordance with (15), the kinetics 𝒫𝐡𝑠𝑝 of 𝐡 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 (π‘Ÿ0>𝑅𝐴𝐡); later on repulsive forces cannot manifest themselves due to the instantaneous decay of the pair on the approach of reactants to the distance 𝑅𝐴𝐡.

Assuming the hopping motion of 𝐡 reactant, the exactly solvable many-particle model of the β€œscavenger problem” gives the recipe for the calculation of the quantities measured: ̇𝒫𝑃(𝑑)=𝐾𝑔(ξ‚πœ)𝑁(𝜏),𝒫𝑃(ξ€œπ‘‘)=𝑑0π‘‘πœπΎπ‘”(ξ‚π’«πœ)𝑁(𝜏),(46)𝐡(𝑑)=Ξ©(𝑑)𝑁(𝑑)π‘’βˆ’π‘‘/𝜏𝐡+1πœπ΅ξ€œπ‘‘0ξ‚π‘‘πœπ‘(π‘‘βˆ’πœ)π‘’βˆ’(π‘‘βˆ’πœ)/πœπ΅π’«π΅(𝜏),(47) where 𝐾𝑔(𝑑) is the pure geminate reaction rate defined in (10), (22), Ω(𝑑) is defined in (27), 𝑁(𝑑) is the exact many-particle survival probability of immobile 𝐡 reactant surrounded by bulk acceptors (see (31), (38)). The binary kinetics 𝑁(𝑑) (16) goes to it with πœπ΅β†’βˆž.

In this limit of immobile 𝐡 particle recipes (46), (47) are reduced to 𝒫𝑃(ξ€œπ‘‘)=𝑑0ξ‚πΎπ‘‘πœπ‘”(ξ‚πœ)𝑁(𝜏);𝒫𝐡(𝑑)=𝑁(𝑑)Ξ©(𝑑).(48) Comparison of (48) with (14), (15) shows that at 𝜏𝐡=∞ the superposition approximation coincides with the exact solution of the problem. Thus for immobile 𝐡 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 𝐡 the rate of geminate reaction products formation ̇𝒫𝑃(𝑑)=𝐾𝑔(𝑑)𝑁(𝑑) cannot be expressed solely in terms of the kinetics of geminate and bulk reactions. It is the product of geminate reaction rate 𝐾𝑔(𝑑) and the above probability 𝑁(𝑑) formed only during the residence time of 𝐡 particle in a geminate pair before it jumps into the bulk. Thus here the use of pure bulk reaction kinetics 𝑁(𝑑) is physically unjustified.

4.1. Hopping Motion of Scavengers

To analyze errors brought about by superposition approximation at arbitrary mobility of 𝐡 reactants, we make calculations in the framework of the microscopic model chosen. Comparison of the kinetics 𝑁(𝑑) (31) with formula (30) shows that the bulk reaction kinetics 𝑁(𝑑) differs from 𝑁(𝑑) by the exponential factor: 𝑑𝑁(𝑑)=expβˆ’πœ‰πœπ΅ξ‚Άξ‚π‘(𝑑).(49) Accordingly, the accumulation rate ̇𝒫𝑃𝑠𝑝(𝑑) 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 𝒫𝑃(𝑑) of geminate reaction products, use formulae (46), (22), and (31). We have.𝒫𝑃𝑑,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0π‘’βˆ’πœ‰πœ‘(πœ‚,𝜈+1,𝜏),𝜈=πœ‰πœπ΅πœπΆ.(50) The introduced dimensionless parameter 𝜈 has the meaning of the bulk reaction rate with immobile 𝐡 reactant in units of 𝜏𝐡. It follows from its definition that the parameter 𝛾 from (33) tends to 𝜈 with 𝜏𝐡 tending to ∞. The function πœ‘(πœ‚,π‘₯,𝜏) is defined by formula (34). As expression (33) for 𝒫𝑃𝑠𝑝 coincides with (50) up to the substitution of 𝜈 for 𝛾, the ultimate probability of geminate reaction can be obtained from (36): π’«π‘ƒξ€·π‘‘βŸΆβˆž,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0π‘’βˆ’πœ‰ξ‚€βˆ’βˆšexpξ‚πœ‚(𝜈+1).(51) For small values of the parameter πœ‰ the result for the ratio of ultimate probabilities (36) and (51) is easily obtained: π’«π‘ƒπ‘ π‘ξ€·π‘‘βŸΆβˆž,π‘Ÿ0ξ€Έπ’«π‘ƒξ€·π‘‘βŸΆβˆž,π‘Ÿ0ξ€Έξ‚΅βˆ’πœ‰=exp2ξ‚™πœ‚ξ‚Άπœˆ+1,πœ‰β‰ͺ1.(52) 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 𝒫𝑃(𝑑,π‘Ÿ0) and accumulation rate ̇𝒫𝑃(𝑑,π‘Ÿ0) 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 ̇𝒫𝑃𝑠𝑝(𝑑,π‘Ÿ0) as compared to the exact solution denoted by lines with circles.

Before passing to the decay kinetics 𝒫𝐡(𝑑) of 𝐡 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 𝐡, correlations are essential, and the contribution of geminate reaction into complete kinetics is insignificant. When 𝐡 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 of 𝐡 reactant, we solve (47) by a conventional method of the Laplace transformation using expressions (27) and (31) for Ω(𝑑,π‘Ÿ0) and 𝑁(𝑑). This gives𝒫𝐡𝑑,π‘Ÿ0ξ€Έ=𝑁𝐡𝑅(𝑑)1βˆ’π΄π΅π‘Ÿ0πœ‘ξ€·πœ‚,π‘’βˆ’πœ‰ξ€Έξ‚Ή,𝜏.(53) Here 𝑁𝐡(𝑑) is the exact many-particle kinetics of the bulk reaction representable as𝑁𝐡𝑑(𝑑)=expβˆ’πœ‰βˆ’πœ‰πœπΆβˆ’ξ€·1βˆ’π‘’βˆ’πœ‰ξ€Έπ‘‘πœπ΅ξ‚Άβ‰‘ξ‚π‘(𝑑)𝑁0(𝑑).(54) Formula (54) is physically clear. The bulk reaction kinetics is defined by two independent processes. The first process described by 𝑁(𝑑) (31) is the reaction between mobile acceptors and immobile 𝐡 particle, while the second one, 𝑁0(𝑑), defines the reaction of mobile 𝐡 particle with the ensemble of immobile acceptors. A regular very large jump of 𝐡 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/𝜏𝐡 and the many-particle probability of acceptors entering the reaction zone, 1βˆ’π‘’βˆ’πœ‰. At small parameters πœ‰ expansion (54) gives exactly the formula of binary kinetics (30) depending on a relative mobility of 𝐡 and 𝐢 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 𝐡 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 𝑁𝐡(𝑑) which is the solution of (47) in the absence of geminate recombination (Ω(𝑑)≑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: Ξ©eff𝑑,π‘Ÿ0𝑅=1βˆ’π΄π΅π‘Ÿ0πœ‘ξ€·πœ‚,π‘’βˆ’πœ‰ξ€Έ,𝜏.(55) This is also confirmed by the absence of the dependence Ξ©eff(𝑑,π‘Ÿ0) on 𝐢 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(𝑑,π‘Ÿ0) as compared to pure geminate kinetics Ξ©(𝑑,π‘Ÿ0). Correspondingly, the stationary portion of nonreacting geminate pairs (the escape probability) is less than in (26): Ξ©eο¬€βˆžξ€·π‘Ÿ0𝑅=1βˆ’π΄π΅π‘Ÿ0ξ‚€βˆ’βˆšexpπœ‚π‘’βˆ’πœ‰ξ‚.(56) 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𝐾eff𝑔𝑑,π‘Ÿ0ξ€Έ=𝐾𝑔𝑑,π‘Ÿ0ξ€Έξ‚΅βˆ’π‘‘expπœπ΅π‘’βˆ’πœ‰ξ‚Ά.(57) 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 𝒫𝐡(𝑑,π‘Ÿ0) as the product of the bulk kinetics and some independent effective geminate kinetics is physically unjustified.

Another interpretation of the rate 𝐾eff𝑔(𝑑,π‘Ÿ0) (57) is given by the representation: 𝐾eff𝑔𝑑,π‘Ÿ0ξ€Έ=𝐾𝑔𝑑,π‘Ÿ0𝑁(𝑑)𝑁𝐡(𝑑).(58) It shows that the increase in 𝐾eff𝑔(𝑑,π‘Ÿ0) is related to defect of the superposition approximation that underestimates significantly geminate reaction product accumulation rate due to exponential difference between 𝑁(𝑑) and 𝑁(𝑑) (see (49)). The dependence of 𝐾eff𝑔(𝑑,π‘Ÿ0) on the bulk reaction kinetics 𝑁𝐡(𝑑) 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 𝒫𝐡(𝑑,π‘Ÿ0) and 𝒫𝐡𝑠𝑝(𝑑,π‘Ÿ0) 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 𝑁𝐡(𝑑) (54).

4.2. Diffusion Motion of Scavengers

For diffusion motion of scavengers comparison between binary formulae (39) and (40) shows that the bulk reaction kinetics 𝑁(𝑑) differs from 𝑁(𝑑) by the exponential multiplier due to the difference between π‘˜ and Μƒπ‘˜, while the second power factor appears due to nonstationary parts of kinetics (39) and (40): 𝑁(𝑑)ξ‚ξ‚ƒβˆ’ξ‚€βˆšπ‘(𝑑)=expπœ‰+ξ‚πœξ‚„ξ‚Έ3πœ‰π›½1βˆ’πœ‰βˆ’π‘“(𝜏)̃𝑛(𝜏)ξ‚Ή,[𝐢]Μƒ1βˆ’πœ‰βˆ’Μƒπ‘›(𝜏)𝛽=π‘˜πœπ΅.(59) Here 𝛽 and 𝜏=𝑑/𝜏𝐡 are diffusion reaction rate with immobile 𝐡 reactant and time in 𝜏𝐡 units, respectively. We introduce the notation for diffusion nonstationary part of kinetics (38), (39): ̃𝑛(𝜏)=23πœ‰π›½πœπœ‹.(60) Accordingly, the accumulation rate ̇𝒫𝑃𝑠𝑝(𝑑) 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 𝒫𝑃(𝑑) of geminate reaction products, use formulae (46), (22), and (39). We have: 𝒫𝑃𝑑,π‘Ÿ0ξ€Έ=(1βˆ’πœ‰)𝒫𝑃(π‘š)𝑑,π‘Ÿ0ξ€Έβˆ’π’«π‘ƒ(𝑛)𝑑,π‘Ÿ0ξ€Έ,(61) where 𝒫𝑃(π‘š) is the contribution of the Markovian (exponential) term of (39) calculated in the previous section: 𝒫𝑃(π‘š)𝑑,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0πœ‘(πœ‚,𝛽+1,𝜏),(62) where the function πœ‘(πœ‚,π‘₯,𝜏) is defined in (34). The value 𝒫𝑃(𝑛)(𝑑,π‘Ÿ0) is determined by diffusion non-Markovian part of kinetics (39) and is equal to 𝒫𝑃(𝑛)𝑑,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0√3πœ‹βˆšξ€œπœ‰πœ‚π›½0(𝛽+1)𝜏dπ‘₯π‘₯ξ‚΅expβˆ’π‘₯βˆ’πœ‚(𝛽+1)ξ‚Ά4π‘₯.(63) In the limit π‘‘β†’βˆž the integral in (63) goes to the Bessel function. In this limit, using (35), we can rewrite (61) in an analytical form: π’«π‘ƒξ€·βˆž,π‘Ÿ0ξ€Έ=π‘…π΄π΅π‘Ÿ0(1βˆ’πœ‰)π‘’βˆ’βˆšπœ‚(𝛽+1)βˆ’2√3πœ‹βˆšπœ‰πœ‚π›½πΎ0ξ‚€βˆšξ‚ξƒ­,πœ‚(𝛽+1)(64) where 𝐾0(π‘₯) is the modified Bessel function of a zero order [59].

The behaviour of the exact accumulation kinetics 𝒫𝑃(𝑑) (61)–(63) and superposition approximation 𝒫𝑃𝑠𝑝(𝑑) (43)–(45) is shown in Figure 4.

The ratio of (43) and (61) at large 𝜏 and small πœ‰ is recast asπ’«π‘ƒπ‘ π‘ξ€·βˆž,π‘Ÿ0ξ€Έπ’«π‘ƒξ€·βˆž,π‘Ÿ0ξ€ΈβŽ‘βŽ’βŽ’βŽ£βˆ’πœ‰=exp2ξ‚™πœ‚βŽ›βŽœβŽœβŽξƒŽπ›½+11+3π›½πœ‰βŽžβŽŸβŽŸβŽ βŽ€βŽ₯βŽ₯⎦Θ,πœ‰β‰ͺ1.(65) 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: Θ=1βˆ’πœ‰βˆ’π’«π‘ƒ(𝑛)π‘ π‘ξ€·βˆž,π‘Ÿ0ξ€Έ/𝒫𝑃(π‘š)π‘ π‘ξ€·βˆž,π‘Ÿ0ξ€Έ1βˆ’πœ‰βˆ’π’«π‘ƒ(𝑛)ξ€·βˆž,π‘Ÿ0ξ€Έ/𝒫𝑃(π‘š)ξ€·βˆž,π‘Ÿ0ξ€Έ.(66) 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 𝒫𝐡(𝑑) of 𝐡 reactant, we solve (47) by a conventional method of the Laplace transformation using expressions (27) and (39) for Ω(𝑑,π‘Ÿ0) and 𝑁(𝑑). 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: 𝒫𝐡(𝑑,π‘Ÿ0)=Ω∞(π‘Ÿ0)𝑁(𝑑) [35]. Taking into account this requirement, we obtain that𝒫𝐡𝑑,π‘Ÿ0ξ€Έ=(1βˆ’πœ‰)𝒫𝐡(π‘š)𝑑,π‘Ÿ0ξ€Έβˆ’π’«π΅(𝑛)𝑑,π‘Ÿ0ξ€Έ,(67) where the Markovian part of the kinetics can be represented as𝒫𝐡(π‘š)𝑑,π‘Ÿ0ξ€Έ=π‘’βˆ’[𝐢]π‘˜π‘‘ξ‚Έπ‘…1βˆ’π΄π΅π‘Ÿ0ξ‚Ή,βˆšπœ‘(πœ‚,𝛼,𝜏)𝛼=1βˆ’πœ‰βˆ’3πœ‰π›½.(68) This result coincides with the corresponding result of (53) at a small value of πœ‰ except the parameter 𝛼. At large 𝑑 expression (68) goes to the following: 𝒫𝐡(π‘š)𝑑,π‘Ÿ0ξ€Έβ‰ƒπ‘’βˆ’[𝐢]π‘˜π‘‘ξ‚Έπ‘…1βˆ’π΄π΅π‘Ÿ0π‘’βˆ’βˆšπ›Όπœ‚ξ‚Ή=π‘’βˆ’[𝐢]π‘˜π‘‘Ξ©eο¬€βˆžξ€·π‘Ÿ0ξ€Έ,π‘‘β‰«πœπ΅.(69) The non-Markovian part of the kinetics can be represented through quadrature of elementary functions.

It is interesting to compare 𝒫𝐡 and 𝒫𝐡𝑠𝑝 (15), (24), (40) at a long time when 𝐡 reactant leaves a geminate pair for the bulk. In this case the expressions for these kinetics are simplified, and we obtain the relation: 𝒫𝐡𝑠𝑝𝑑,π‘Ÿ0𝒫𝐡𝑑,π‘Ÿ0ξ€Έβ‰ƒΞ©βˆžΞ©eο¬€βˆžβŽ›βŽœβŽœβŽπ‘…1βˆ’π΄π΅π‘Ÿ0ξ‚€π‘’βˆ’βˆšπ›Όπœ‚ξ‚βˆšβˆ’π½(π›Όπœ‚)3πœ‰π›½/π›Όβˆš(1βˆ’πœ‰βˆ’3πœ‰π›½/2)Ξ©eο¬€βˆžβŽžβŽŸβŽŸβŽ βˆ’1,π‘‘β‰«πœπ΅,(70) where Ω∞ is the escape probability defined by (26), and 𝐽(π‘₯) is the following quadrature: 2𝐽(π‘₯)=βˆšπœ‹ξ€œβˆž0√d𝑦𝑦exp(βˆ’π‘¦)erfcπ‘₯ξ‚Ά4𝑦.(71) It easy to see that in (70) the factor Ω∞/Ξ©eο¬€βˆž is defined by the ratio of the Markovian parts of 𝒫𝐡𝑠𝑝 and 𝒫𝐡, 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: Δ𝑃dot=ξ€·Μ‡π’«π‘ƒβˆ’Μ‡π’«π‘ƒπ‘ π‘ξ€ΈΜ‡π’«π‘ƒ,Δ𝑃=ξ€·π’«π‘ƒβˆ’π’«π‘ƒπ‘ π‘ξ€Έπ’«π‘ƒ,Δ𝐡=ξ€·π’«π΅βˆ’π’«π΅π‘ π‘ξ€Έπ’«π΅.(72) 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: [𝐢]π‘˜π‘‘β‰ͺπœ‰βˆ’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 Μƒ[𝐢](π‘˜βˆ’π‘˜)𝑑β‰ͺπœ‰βˆ’1, while the smallness of the second type corrections requires (as is seen from (39)) that the condition Μƒ[𝐢]π‘˜π‘‘β‰ͺπœ‰βˆ’1 be fulfilled. Thus the binary time interval where comparison of the theories will be made is given by the criterion: 𝑑β‰ͺ𝑑bin=πœ‰βˆ’1[𝐢]ξ€·Μƒπ‘˜minπ‘˜βˆ’ξ€Έξ€Έβˆ’1,ξ€·[𝐢]Μƒπ‘˜ξ€Έβˆ’1.(73) 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 ̇𝒫𝑃 and ̇𝒫𝑃𝑠𝑝 for the hopping and diffusion models of scavenger motion in Δ𝑃dot, we have, correspondingly, Δ𝑃dot=⎧βŽͺ⎨βŽͺ⎩1βˆ’π‘’βˆ’πœ‰πœ,𝜏β‰ͺπœ‰βˆ’2,1βˆ’π‘’βˆšβˆ’(πœ‰+3πœ‰π›½)πœξ‚Έ1βˆ’πœ‰βˆ’π‘“(𝜏)̃𝑛(𝜏)ξ‚Ή11βˆ’πœ‰βˆ’Μƒπ‘›(𝜏),𝜏β‰ͺπœ‰ξ‚€βˆšπœ‰+.3πœ‰π›½(74) It is seen that exponential deviation ̇𝒫𝑃𝑠𝑝 from ̇𝒫𝑃 leads to the fact that even at half-decay times π‘‘β‰ˆπœπΆ/πœ‰ of the kinetics 𝑁(𝑑) (31) the relative deviation Δ𝑃dot may differ from unity just by several tens percent. In particular, for the hopping model of scavenger motion Δ𝑃dot=0.63;0.86;0.95 for 𝜏𝐡=𝜏𝐢, 𝜏𝐡=𝜏𝐢/2, 𝜏𝐡=𝜏𝐢/3, respectively. For diffusion model of scavenger motion these distinctions are even more essential. It should be noted that Δ𝑃dot does not depend on initial distance in geminate pair π‘Ÿ0. Thus, (74) are correct for any distribution of initial distance 𝑃(π‘Ÿ0).

The behaviour of the relative deviation Δ𝑃dot 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)).

Consider the behaviour of the function Δ𝑃 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 Δ𝑃 increases monotonically with time and attains its steady-state value Δ𝑠𝑃 rather rapidly. The behaviour of the relative deviation Δ𝑃 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 Δ𝑃 (line 2) and its Markovian version Δ𝑃(π‘š)=1βˆ’π’«π‘ƒ(π‘š)𝑠𝑝/𝒫𝑃(π‘š) (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).

Using the expansion of the functions 𝒫𝑃𝑠𝑝(∞,π‘Ÿ0)/𝒫𝑃(∞,π‘Ÿ0) 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: Δ𝑠𝑃=⎧βŽͺ⎨βŽͺβŽ©ξ‚΅βˆ’πœ‰1βˆ’exp2ξ‚™πœ‚ξ‚Άξ‚΅βˆ’πœ‰πœˆ+1,hopping,1βˆ’exp2ξ‚™πœ‚ξ‚΅ξ‚™π›½+11+3π›½πœ‰ξ‚Άξ‚ΆΞ˜,diffusion,(75) where the function Θ is given by relation (66).

It is easy to see that at the start at the contact (π‘Ÿ0=𝑅𝐴𝐡,πœ‚=0) the function Δ𝑠𝑃 is equal to zero. Then, as πœ‚ increases, Δ𝑠𝑃 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 Δ𝑠𝑃. Indeed, the use of the binary kinetics in (14) calls for the smallness of the characteristic decay time πœπ‘” (23) of a geminate pair as compared to the binary interval (πœπ‘”β‰ͺ𝑑bin). As a result, at πœ‚β‰«1, when πœπ‘”βˆš=3Μƒπœπ‘”πœπ΅, 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 Δ𝑠𝑃 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: πœ‰2ξ‚™πœ‚=πœ‰πœˆ+16πœπ‘”πœπ΅1√1+πœ‰πœπ΅/𝜏𝐢>1.(76) For acceptors at rest (πœπΆβ†’βˆž) condition (76) is realized at πœ‰πœπ‘”πœπ΅>1.(77) Physically this means that geminate pair lifetime is long enough for the bulk reaction (proceeding at the rate πœ‰/𝜏𝐡) to reach considerable depth. Time correlations in this case are maximum. When 𝐢 reactant starts moving, correlations decrease. In the region πœ‰πœπ΅β‰ͺ𝜏𝐢β‰ͺ𝜏𝐡, where the bulk reaction rate is πœ‰/𝜏𝐢 (since this rate can be calculated neglecting 𝐡 reactant motion), but the bulk reaction depth for the time 𝜏𝐡(𝜈β‰ͺ1) is small, correlations become essential under more rigid restrictions: πœ‰πœπ‘”πœπΆ>πœπ΅πœπΆβ‰«1.(78) As the motion of 𝐢 reactants becomes more active (𝜏𝐢β‰ͺπœ‰πœπ΅), the region of essential deviations is reached under even more strict conditions: πœ‰πœπ‘”πœπΆ>πœπ΅πœπΆξ‚™πœ‰πœπ΅πœπΆβ‰«πœπ΅πœπΆβ‰«1.(79) For diffusion motion of scavengers the condition of large relative deviations is as follows: πœ‰2ξ‚™πœ‚βŽ›βŽœβŽœβŽξƒŽπ›½+11+3π›½πœ‰βŽžβŽŸβŽŸβŽ =πœ‰6πœπ‘”πœπ΅ξ‚€βˆš1+3𝜏𝐡/πœπΆξ‚βˆš1+πœ‰πœπ΅/𝜏𝐢>1,(80) where the mean residence time 𝜏𝐢=𝑅2𝐡𝐢/(3𝐷𝐢) of diffusing scavengers in the reaction zone is introduced. Obviously, with πœπΆβ†’βˆž 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 𝜏𝐢 lies in the interval πœ‰πœπ΅β‰ͺ𝜏𝐢β‰ͺ𝜏𝐡. Then the second term between the brackets in (80) becomes decisive, and this equality may be written as πœ‰πœπ‘”πœπΆ>ξ‚™πœπ΅πœπΆβ‰«1.(81) 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 πœ‰πœπ‘”/𝜏𝐢β‰ͺπœ‰βˆ’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: πœ‰βˆ’1πœβ‰«πœ‰π‘”πœπΆ>ξ‚™πœπ΅πœπΆβ‰«1.(82) Finally, when the motion of scavengers becomes dominant (𝜏𝐢β‰ͺπœ‰πœπ΅), condition (80) is met in the following restricted interval: πœ‰βˆ’1πœβ‰«πœ‰π‘”πœπΆ>βˆšπœ‰πœπ΅πœπΆβ‰«1.(83) Obviously, large relative deviations can be realized only on condition that the right and the left boundaries do not overlap: πœπΆβ‰«πœ‰3/2𝜏𝐡. With progressively decreasing 𝜏𝐢 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 (πœπΆβ†’0) the difference between 𝑁(𝑑) and 𝑁(𝑑) 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 𝐢 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 Δ𝑠𝑃 as a function of parameter πœ‚ in the case of dominate electron motion (πœπΆβ‰«πœπ΅) 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 Δ𝑠𝑃 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.

Now we pass to the analysis of the function Δ𝐡. Since the superposition approximation retards the decay kinetics of 𝐡 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: Δ𝐡𝑅=1βˆ’1βˆ’π΄π΅/π‘Ÿ0ξ€Έπœ‘ξ€»(πœ‚,1,𝜏)𝑅1βˆ’π΄π΅/π‘Ÿ0ξ€Έξ€»πœ‘(πœ‚,1βˆ’πœ‰,𝜏).(84) 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 Δ𝑠𝐡 defined by asymptotic (35) of the function πœ‘(πœ‚,π‘₯,𝜏). In the range of small values of πœ‰ the expression for Δ𝑠𝐡 can be brought to the form Δ𝑠𝐡𝑅=βˆ’π΄π΅π‘Ÿ0π‘’βˆ’βˆšπœ‚ξ‚€π‘’πœ‰βˆšπœ‚/2ξ‚βˆ’1𝑅1βˆ’π΄π΅/π‘Ÿ0ξ€Έπ‘’βˆ’βˆšπœ‚(1βˆ’πœ‰/2).(85) Examine the behaviour of the function Δ𝑠𝐡 (85). If a geminate pair is formed not at the contact (π‘Ÿ0>𝑅𝐴𝐡), 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 π‘Ÿ0 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 Δ𝑠𝐡 may be simplified: Ξ”π‘ π΅πœ‰=βˆ’2π‘…π΄π΅π‘Ÿ0βˆšπœ‚π‘’βˆšπœ‚βˆ’ξ€·π‘…π΄π΅/π‘Ÿ0βˆšξ€Έξ€·1+πœ‰ξ€Έπœ‚/2.(86) Function (86) has the minimum at the point πœ‚π‘š determined by a transcendental equation: π‘…π΄π΅π‘Ÿ0π‘’βˆ’βˆšπœ‚π‘šβˆš=1βˆ’πœ‚π‘š.(87) This equation has the only root lying in the interval 0β‰€πœ‚π‘š<1. Substituting (87) in (86), we have with the accuracy up to the main order in πœ‰: minΞ”π‘ π΅πœ‰β‰ƒβˆ’2ξ€·βˆš1βˆ’πœ‚π‘šξ€Έ.(88) The most essential deviations |Δ𝑠𝐡|β‰ƒπœ‰/2 arise at πœ‚π‘šβ‰ͺ1. Examining the function Δ𝑠𝐡 in this interval, we easily obtain the analytical expression for the minimum: πœ‚π‘šξ‚΅π‘…β‰ˆ21βˆ’π΄π΅π‘Ÿ0ξ‚Άβ‰ͺ1.(89) 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: π‘Ÿ0βˆ’π‘…π΄π΅π·β‰ƒ2π΄πœπ΅π‘…π΄π΅β‰ͺ𝑅𝐴𝐡.(90) 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 𝑅𝐴𝐡 is 𝛿𝑅𝐴𝐡/𝐷𝐴 [22]. Really, expression (90) can be recast as ξ€·π‘Ÿ0βˆ’π‘…π΄π΅ξ€Έπ‘…π΄π΅2π·π΄β‰ƒπœπ΅.(91) It follows from the above relation (91) that maximum time correlations are realized when the residence time of diffusing 𝐴 reactant in a thin layer (π‘Ÿ0βˆ’π‘…π΄π΅)/2 is equal to the time 𝜏𝐡 between 𝐡 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.

Figure 8 shows the dependence of the relative deviation Δ𝑠𝐡 versus parameter πœ‚. Underestimating the rate of geminate reaction, the superposition approximation thus overestimates the probability of 𝐡 reactant to survive and escape into the bulk. Apparently, Δ𝑠𝐡 is always negative. The largest relative deviations are observed for small πœ‚ when the contribution of geminate reaction to total kinetics 𝒫𝐡(𝑑) is most substantial.

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 𝐡 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 ̇𝒫𝑃 of geminate reaction products and the decay kinetics 𝒫𝐡 of 𝐡 reactant are expressed as the product of bulk kinetics (31), (38) and, respectively, geminate reaction rate (22) and geminate kinetics (27): ̇𝒫𝑃𝐾(𝑑)=𝑁(𝑑)𝑔(𝑑),𝒫𝐡1(𝑑)=𝑁(𝑑)Ξ©(𝑑),𝜏𝐡=0.(92) As soon as 𝐡 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 𝐡 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 𝐾𝑔(𝑑) (22) and the kinetics 𝑁(𝑑) formed solely during the residence time of 𝐡 particle in a geminate pair: ̇𝒫𝑃(𝑑)=𝑁(𝑑)𝐾𝑔(𝑑). As a result, the reaction rate in the geminate channel differs strongly from the prediction following from the superposition approximation: ̇𝒫𝑃𝑠𝑝(𝑑)=𝐾𝑔(𝑑)𝑁(𝑑). 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 𝒫𝐡(𝑑) of 𝐡 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 𝐹 on the concentration [𝐢] of bulk acceptors is measured with the subsequent determination of the geminate pair kinetics Ξ©(𝑑) from the measured concentration dependence 𝐹(𝛼[𝐢]). 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: 𝒫𝐷𝑐𝑑[𝐢]π‘˜ξ€œ(π‘‘βŸΆβˆž)=∞0𝑑𝑑Ω(𝑑)π‘’βˆ’[𝐢]π‘˜π‘‘[𝐢]β‰ˆπΉ(𝛼).(93) The inverse Laplace transformation of (93) with [𝐢]π‘˜ as transformation variable gives the geminate kinetics Ξ©(𝑑). Its characteristic decay time is determined from the ratio πœπ‘”βˆΌπ›Ό/π‘˜, with the rate constant value π‘˜ 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 [𝐢]π‘˜ (equal to πœ‰(1/𝜏𝐢+1/𝜏𝐡) for hopping and πœ‰(1/𝜏𝐢+1/𝜏𝐡+√3/(𝜏𝐢𝜏𝐡)) for diffusing scavengers), the role of the Laplace variable will be played by the reaction rate with immobile 𝐡 reactant: Μƒ[𝐢]π‘˜=πœ‰/𝜏𝐢 that can be essentially less than [𝐢]π‘˜. Thus the value πœπ‘” found by formula (93) is Μƒπ‘˜π‘˜/ 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 π‘˜ is determined mainly by electron mobility. Indeed, using both the lifetime of the excess electron in the localized state in cyclohexane 𝜏𝐡∼10βˆ’10 s [43] and typical values 𝑅𝐡𝐢=10 Å and 𝐷𝐢=10βˆ’5 cm2/s, we get Μƒπ‘˜/π‘˜β‰ˆ7.5. Note that this value is almost twice greater than the corresponding ratio Μƒπ‘˜π‘˜/ for the hopping motion of acceptors. Since for our model, according to (23), πœπ‘”βˆšβ‰ˆ3Μƒπœπ‘”πœπ΅=3(π‘Ÿ0βˆ’π‘…π΄π΅)√𝜏𝐡/𝐷𝐴, the value of the initial distance in the geminate pair decreases similarly. In real situation when the initial distance in geminate pair π‘Ÿ0 is not fixed but is distributed according to 𝑃(π‘Ÿ0) the characteristic time of geminate reaction is proportional to the width 𝐿 on the distribution: πœπ‘”βˆšβ‰ˆπΏπœπ΅/𝐷𝐴. 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: ̇𝒫𝑃𝑠𝑝(𝑑)=𝐾𝑔(𝑑)𝑁(𝑑). It is seen from this formula that the bulk kinetics 𝑁(𝑑) may be obtained by dividing the experimentally measured fluorescence rate ̇𝒫𝑃𝑠𝑝(𝑑) of geminate reaction products in the presence of acceptors by the fluorescence rate 𝐾𝑔(𝑑) of pure geminate reaction products. However, as follows from exact solution (46), the kinetics 𝑁(𝑑) 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 𝑁(𝑑) and kinetics 𝑁(𝑑) the relative deviations Δ𝑃dot 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.

Acknowledgments

The authors are grateful to the Russian Foundation of Fundamental Research (projects 05-03-32651, 09-03-00456) and to the Program of Fundamental Research of Chemical Division and Material Science RAS (project 5.1.5) for the support of this work.