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 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 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: 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 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: 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]: The last relation is simplified in the stationary limit , thus the experimentally observed ultimate probability of recombination and total scavenging yield obey the relation: 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: As usual, geminate reaction rate in the absence of acceptors is expressed in terms of the elementary rate : where is the initial distribution of reactants in geminate pairs, and is the probability density to find and reactants at the relative distance at the instant of time if at the initial moment they were at the distance (the Green function). It obeys the equation with the initial condition 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: 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: while reactant decay kinetics is expressed as the product of bulk and geminate kinetics: 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]: 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 Stationary rate constant of the decay on acceptors appearing in formulae (8) and (9) is determined from (17): The Green function of the pair is defined by the equation 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 (the so-called two-scale migration [54]): Here is a unit operator.
Starting with the initial distance we easily obtain that at geminate reaction rate in the absence of bulk acceptors is defined as where 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 . 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 . So for , The use of (22) in (13) with the subsequent integration gives the expression for the survival probability: where the introduced parameter 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 : Only in the limit it is equal to the probability of the diffusion entry of reactants into the bulk: . In this limiting case of immobile () expression (24) is easily seen to turn into the well-known diffusion result [21]: 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 ββ«. As follows from [55, 56], in this case it will suffice to perform the following substitution: because in highly polar solvents the survival probability of charged geminate pairs is kinetically identical to the corresponding value 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: 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 where the dimensionless parameter (the density parameter) appears. It is equal to the average number of acceptors in the reaction zone of the volume . 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 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 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: 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: where and are bulk reaction rate and time in units, respectively.
Formula (33) involves the function : 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: This gives the expression for the ultimate probability of geminate reaction in the presence of bulk acceptors: 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: 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 the βtarget problemβ kinetics proceeds as the well-known Smoluchowski kinetics [53]: where 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 : The validity condition of formula (39) defines the binary time interval 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 where, in contrast to the case of diffusion mobilities of both reactants, stationary rate constant of the bulk two-scale reaction is and is a function decreasing from unity to zero: At the time it goes to its asymptotic , 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 where the Markovian part of the kinetics corresponding to pure exponential term of (40) is: 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 , where .
The non-Markovian part of the kinetics is defined by the integral 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 (); 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: 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 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: 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. 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): For small values of the parameter the result for the ratio of ultimate probabilities (36) and (51) is easily obtained: 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 and accumulation rate 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 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 and . This gives Here is the exact many-particle kinetics of the bulk reaction representable as 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, , 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 and the many-particle probability of acceptors entering the reaction zone, . 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 (). 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: This is also confirmed by the absence of the dependence on reactant mobility.
As one would expect, the influence of correlations disappears under infinite dilution of bulk scavengers: with expression (55) goes into (24). However, it is seen that the presence of acceptors in (55) results in the acceleration of the kinetics as compared to pure geminate kinetics . Correspondingly, the stationary portion of nonreacting geminate pairs (the escape probability) is less than in (26): 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 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 as the product of the bulk kinetics and some independent effective geminate kinetics is physically unjustified.
Another interpretation of the rate (57) is given by the representation: It shows that the increase in 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 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 and 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): 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): 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: where is the contribution of the Markovian (exponential) term of (39) calculated in the previous section: where the function is defined in (34). The value is determined by diffusion non-Markovian part of kinetics (39) and is equal to 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: where 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 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: 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 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: [35]. Taking into account this requirement, we obtain that where the Markovian part of the kinetics can be represented as 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: 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: where is the escape probability defined by (26), and is the following quadrature: It easy to see that in (70) the factor 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: 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: . 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 , while the smallness of the second type corrections requires (as is seen from (39)) that the condition be fulfilled. Thus the binary time interval where comparison of the theories will be made is given by the criterion: 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 , we have, correspondingly, It is seen that exponential deviation from leads to the fact that even at half-decay times of the kinetics (31) the relative deviation may differ from unity just by several tens percent. In particular, for the hopping model of scavenger motion for , , , respectively. For diffusion model of scavenger motion these distinctions are even more essential. It should be noted that does not depend on initial distance in geminate pair . Thus, (74) are correct for any distribution of initial distance .
The behaviour of the relative deviation 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 (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 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: where the function is given by relation (66).
It is easy to see that at the start at the contact () 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 (). As a result, at , when , an insignificant upper bound is set on the parameter . In the hopping case this restriction is defined by the inequality . 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: For acceptors at rest () condition (76) is realized at 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 is small, correlations become essential under more rigid restrictions: As the motion of reactants becomes more active (), the region of essential deviations is reached under even more strict conditions: For diffusion motion of scavengers the condition of large relative deviations is as follows: where the mean residence time 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 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: . This is equivalent to the condition . 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: Finally, when the motion of scavengers becomes dominant (), condition (80) is met in the following restricted interval: Obviously, large relative deviations can be realized only on condition that the right and the left boundaries do not overlap: . 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 () 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: 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 Examine the behaviour of the function (85). If a geminate pair is formed not at the contact (), then under variation of the parameter from 0 to β this function begins with zero, passes through the minimum in the area , and then smoothly tends to zero again with increasing . The closer the starting point 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 expression (85) for may be simplified: Function (86) has the minimum at the point determined by a transcendental equation: This equation has the only root lying in the interval . Substituting (87) in (86), we have with the accuracy up to the main order in : The most essential deviations arise at . Examining the function in this interval, we easily obtain the analytical expression for the minimum: 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: 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 It follows from the above relation (91) that maximum time correlations are realized when the residence time of diffusing reactant in a thin layer 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): 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: 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 for hopping and 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 βs [43] and typical values ββ« and βcm2/s, we get . 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), , the value of the initial distance in the geminate pair decreases similarly. In real situation when the initial distance in geminate pair is not fixed but is distributed according to 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 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.