The Kinetics of Joined Action of Triplet-Triplet Annihilation and First-Order Decay of Molecules in the 𝑇₁ State in the Case of Nondominant First-Order Process

Abstract

This paper presents the description of triplet-triplet annihilation in the case of nondominant first-order decay of molecules in the triplet state. The kinetics of the statistical system is influenced by joined action of two processes: the first- and second-order decies. This kinetics can be described with analytical function if the rate parameter of second-order reaction is constant. The approach presented here combines the well-known Smoluchowski formula with the previously published intuitive and non-Fickian models of diffusion-controlled triplet-triplet annihilation. The kinetics of the delayed fluorescence of anthracene is used as a practical example of applicability of the model proposed. The advantages and limits of the proposed model are discussed.

1. Introduction

The reaction where the diffusion plays an important role was studied for many years [1]. Nowadays studied diffusion-controlled reactions like fluorescence quenching [2] very often take place on the picosecond timescale. Although the Smoluchowski-Collins-Kimball (SCK) model shows some difficulties in reproduction of the kinetic parameters, this approach or its variations are used for the interpretation of the kinetic data. An example of the limits of the SCK model in the description of kinetics of fluorescence quenching is the treatment of the solvent as continuous medium. This approximation seems not to be valid on the picosecond timescale [3]. The other inaccuracy arising from the continuous-medium approximation is reflected in the value of the effective annihilation radius. The effect of the screening by the first coordination shell gives this value 20% to 40% larger than the value calculated from dimensions of the molecule [4]. The other example of more complicated kinetics than that predicted by the SCK model is the intermolecular charge transfer (CT) [5, 6]. In this case the important role is played not only by the translational diffusion of the molecules but the rotational one as well. The rotational diffusion influences the mutual orientation of the molecules. And the mutual orientation plays an important role in the efficiency of charge-transfer reaction. The authors of [5, 6] used three exponential phenomenological models for the description of the kinetics of the system with nonnegligible rotational diffusion. The diffusion also plays an important role in the kinetics of intermolecular proton-transfer (PT) reaction [7–15]. For example, in the case of PT from 5-cyano-2 naphthol in methanol/water mixture the nonexponential kinetics is observed [7]. The kinetics of the diffusion-influenced PT is generally described by Debye-Smoluchowski-Einstein (DSE) equation. The most general description is realized with the numerical solution of the DSE equation [11–13]. In some cases analytical or approximate solutions of DSE equation also offer reasonable description of the PT kinetics [12, 13]. In the case of large delay the concentration of the reacting species is so low that the diffusion plays an important role. In such conditions the decay of the fluorescence should be described by power function instead of the exponential one [9, 10], for example in the case of derivatives of 1- and 2-naphthol the asymptotic decay of the reacting species can be described with power law ~𝑡^−3/2 [8]. There are however deviations from the power law ~𝑡^3/2 in the case of some diffusion-controlled reactions, for example, due to finite rate of recombination within the contact pair in acid-base systems [16]. In such systems the internal rate parameter of the proton transfer reaction can be measured for high concentration of reacting species where the diffusion does not play the most important role [16].

In the case of triplet-triplet annihilation (TTA) in viscous solvents where the process takes place on microsecond timescale the deviations from diffusional Smoluchowski model [17] are not observed directly [18]. It means that the course of the intensity of delayed fluorescence as a function of time is well described by the function derived from Smoluchowski approach. The deviations from this model manifest themselves more subtle. For example, the effect of screening by the first coordination shell [4] can be detected as larger value of annihilation radius than that calculated from molecular dimensions [19]. Also partly reflecting boundaries [20] should give no important contribution, because the movement of the molecules in viscous solvent is very slow. The mutual diffusion coefficient of molecules in typical diffusion-controlled TTA experiment is at least three orders of magnitude lower than the typical diffusion coefficient of organic molecules in organic solvents (10⁻⁹ m² s⁻¹) at room temperature. The Smoluchowski model can be applied to describe the kinetics of statistical system with partly reflecting boundaries. However, instead of the geometrical radius of interaction the effective one has to be used [21, pages (105–117)]. The effective radius of interaction is larger in comparison with the parameter used originally by Smoluchowski. For typical organic molecule in organic solvent at room temperature the effective radius in Smoluchowski approach used to describe the system with partly reflecting boundaries should be 10% larger than the geometrical one. This difference between “real” and “effective” radius of interaction is difficult to be detected in the experiment [21, pages (105–117)]. So in the case of TTA the assumptions of the original Smoluchowski approach seems to be a good approximation: the probability of the reaction for molecules within annihilation radius equals 100% (the reaction takes place immediately). With other words: due to high viscosity of the solvent in TTA experiment the molecules which build the contact pair cannot diffuse away, because they stay in the contact pair long enough to neglect the influence of the internal rate parameter on the kinetics of the whole process.

Also the reacting species can be treated as spherically symmetric. There are two reasons that make this approximation valuable. The first one is the very fast rotation of the molecules in comparison with their translational diffusion. This aspect was discussed in earlier papers [19, 22]. In general case the translational and rotational diffusion should be treated simultaneously. However, very often the correction for rotational diffusion results in a constant factor multiplying the time-dependent rate parameter obtained from the second Fick’s law [23, pages (25-26)]. This kind of correction is automatically included in the amplitude of fitted function which describes the kinetics of the process.

The complicated spatial feature of Dexter mechanism of interaction [24] is not taken into account. Due to its very small spatial range of interaction the influence on the kinetics of the process can be treated as negligible. The influence of the spin statistics on the kinetics of TTA was referred to in previous papers [19, 22]. Also the existence of different TTA channels and their contribution to the kinetics of the TTA was summarized there.

To sum up, in the description of TTA kinetics of in viscous solutions the molecules can be treated as hard spheres which react with each other only when they are in the contact.

In order to make the discussion of the joined action of first-order decay and triplet-triplet annihilation the description is done within basic assumption of the classical Smoluchowki [17] model—mathematically simplest possible treatment of the problem:(i)the Brownian motions of the reacting particles can be described by diffusion in continuous medium; (ii)the molecules at distances larger than annihilation radius do not interact with each other;(iii)the annihilation within contact pair takes place immediately, it means that the molecules in the T₁ state cannot diffuse away if they are within an annihilation radius from each other;(iv)the evolution of the local density 𝜌(r, t) starts from random distribution of the particles in solution 𝜌(𝑟,𝑡=0)=𝜌₀ and tends to 𝜌(𝑟,𝑡→∞).The summarized above assumptions led Smoluchowski to his famous formula for the time-dependent rate parameter 𝑘_2𝐴(𝑡): 𝑘_2𝐴(𝑡)=4𝜋𝑅_𝐴𝐷_𝑇1+√𝜋𝑅𝐷^𝐴_𝑇𝑡,(1) where 𝑅_𝐴 is the annihilation radius and 𝐷_𝑇 is the mutual diffusion coefficient. However, the original Smoluchowski model does not take into account the joined action of two processes: the first-order decay of the molecules in the 𝑇₁ state and the diffusion-controlled TTA. So the correction has to be introduced. In previous papers [19, 22] we concentrated on the case of strongly dominant first-order decay. That condition allowed to introduce the following approximation: the decay of the average density of molecules in the 𝑇₁ state was described by monoexponential function. The Smoluchowski solution of the second Fick’s law was modified in order to introduce the joined action of both processes: the first-order decay and the spatial redistribution, to the expression for time-dependent rate parameter. In the intuitive model the term exp(−𝑘_𝑇𝑡) was introduced to the denominator of so-called “short-time” part of 𝑘_2𝐴(𝑡) [22]. The other approach was presented in [19] where the time-dependent mutual diffusion coefficient in the form 𝐷_𝑇 exp(−2𝑘_𝑇𝑡) was used in the second Fick’s law instead of the constant diffusion coefficient 𝐷_𝑇. The experimental conditions of the measurements that allowed to test the models were very unfavorable. Generally the p-type delayed fluorescence is a few orders of magnitude smaller than the prompt one. In order to obtain good signal-to-noise ratio the measurement is time consuming (the measurement of a single decay took over 10 hours). The investigation of the diffusion-controlled TTA kinetics requires also extremely high stability of the temperature: it must be constant within the range of ±0.01 K. The reason for this extremely high stability is to keep the viscosity of the medium constant within a few percent.

In this paper we extend the modifications of the Smoluchowki time-dependent rate parameter to the case of nondominant first-order decay. The modified models (two variations) are applied to the description of the TTA kinetics of anthracene in the mixture of cis- and trans-dimethylcyclohexane (cis/trans-DMCH).

2. Description of the Model

Let us come back to the analysis of the TTA presented in paper [22]. The reason to use the time-dependent rate parameter 𝑘_2𝐴(𝑡) in the form: 𝑘_2𝐴(𝑡)=4𝜋𝑅_𝐴𝐷_𝑇1+𝑒^{−𝐾𝑇𝑡}𝑅√^𝐴𝜋𝐷_𝑇𝑡,(2) was to introduce the contribution of the first-order decay to the evolution of the local density. Since in that case the first-order reaction was dominant the joined action of two processes, first-order decay and diffusion-controlled TTA, was described by introducing the exponential function in the denominator of component representing so-called “short-time” effect in (1). This led to (2).

In order to extend the intuitive model presented in paper [22] to the case of nondominant first-order decay we carry out similar reasoning as in that paper.

Let us take into account three molecules in the 𝑇₁ state and denote them with letters A, B, and C. Let us look at the kinetics of the system from the point of view of A. The molecule A can react with one of two others (B or C) if the reaction partner is still in the 𝑇₁ state when molecules come to the contact distance. Each of the potential reaction partner for A can decay on the first-order reaction before building the contact pair with A. It is also possible that both B and C decay due to diffusion-controlled TTA: ³𝐵^∗+³𝐶^∗⟶⟶¹₁(𝐵𝐶)^∗∗⟶¹𝐵^∗+¹𝐶or¹𝐵+¹𝐶^∗𝐵+¹𝐶+ℎ𝜈,(3) where ³𝐵^∗ and ³𝐶^∗ denote the molecules in the 𝑇₁ state, ¹(𝐵𝐶)^∗∗ denotes the encounter pair with singlet multiplicity, ¹𝐵^∗ and ¹𝐶^∗ are the molecules in the 𝑆₁ state, and ¹𝐵 and ¹𝐶 are in the ground state; ℎ𝜈 denotes the energy of the emitted photon. Following the same course of reasoning as in previous paper [22] that the kinetics of TTA should be determined by diffusion of molecules in the 𝑇₁ state which is in turn influenced by the decay of the average density of molecules in the triplet state, one can propose the following description of the kinetics. Since the main influence of the joined action can be limited to the “so-called” short-time effect [25] one can approximate the description of the kinetics in the following way:(1)the stationary part of the standard Smoluchowski rate parameter 𝑘_2𝐴(𝑡) remains unchanged;(2)the time-dependent part 𝑅_𝐴/(𝜋𝐷_𝑇𝑡)^1/2 is modified in such a way that takes into account the influence of the decay of the average density on the total amount of diffusing molecules.The denominator in the short-time part of standard Smoluchowski parameter 𝑘_2𝐴(𝑡) is proportional to the total number of diffusing molecules within the period from 𝑡=0 to 𝑡>0 under the assumption, that there are no other processes of decay of molecules. In the case of TTA the amount of diffusing molecules will be decreased due to the decay of the average density. The total amount of molecules which diffuse during the period from 𝑡=0 to 𝑡>0 should be proportional to the (𝜌(𝑡)/𝜌₀)(𝜋𝐷_𝑇𝑡)^1/2. Finally the intuitive expression for time-dependent rate parameter should have the form: 𝑘_2𝐴(𝑡)=4𝜋𝑅_𝐴𝐷_𝑇1+𝜌(𝑡)/𝜌𝑅₀^𝐴√𝜋𝐷_𝑇𝑡.(4)

Now one has to find an expression for 𝜌(𝑡). The start point is the well-known formula for the concentration 𝜌(𝑡) of molecules in the triplet state obtained in the case of time-independent rate constant of the TTA. The expression can be found elsewhere [26]: 𝜌(𝑡)=𝜌₀𝑒^{−𝑘𝑇𝑡}1+𝑘_𝑇𝑇𝜌₀/𝑘_𝑇1−𝑒^{−𝐾𝑇𝑡},(5) where 𝜌₀ is the initial concentration of the molecule in the 𝑇₁ state, 𝑘_𝑇 is the first-order rate constant, and 𝑘_𝑇𝑇 is the second-order rate constant. Since in our case the TTA is diffusion controlled we approximate the 𝜌(𝑡) by replacing 𝑘_𝑇𝑇 with 𝑘_2𝐴(𝑡). The expression for average density has the following form: 𝜌=(𝑡𝜌)₀𝑒^{−𝑘𝑇𝑡}=𝜌₀1+𝑘_2𝐴(𝑡)𝜌₀/𝑘_𝑇1−𝑒^{−𝑘𝑇𝑡}𝑒^{−𝑘𝑇𝑡}1+4𝜋𝑅_𝐴𝐷_𝑇𝜌₀/𝑘_𝑇1+𝑅_𝐴/𝜌(𝑡)/𝜌₀√𝜋𝐷_𝑇𝑡1−𝑒^{−𝑘𝑇𝑡}.(6) Equation (6) can be resolved due to the average density 𝜌(𝑡): 𝜌(𝑡)=𝜌₀𝑒^{−𝑘𝑇𝑡}−4√𝜋𝐷_𝑇𝑅^2𝐴𝜌⁰/𝑘_𝑇√𝑡1−𝑒^{−𝑘𝑇𝑡}1+4𝜋𝐷_𝑇𝑅_𝐴𝜌₀/𝐾_𝑇1−𝑒^{−𝑘𝑇𝑡}.(7) Substituting the ratio 𝜌(𝑡)/𝜌₀ from (7) to (4) one obtains:𝑘_2𝐴(𝑡)=4𝜋𝑅_𝐴𝐷_𝑇×⎛⎜⎜⎝1+𝑅^𝐴𝑒^{−𝑘𝑇𝑡}−4√𝜋𝐷_𝑇𝑅^2𝐴𝜌⁰/𝐾_𝑇√𝑡1−𝑒^{−𝑘𝑇𝑡}/1+4𝜋𝐷_𝑇𝑅_𝐴𝜌₀/𝑘_𝑇1−𝑒^{−𝑘𝑇𝑡}√𝜋𝐷_𝑇𝑡⎞⎟⎟⎠.(8)

The expression for the time-dependent intensity of the delayed fluorescence has the form [27]: 𝐼_𝐷𝐹(𝑡)=𝐶𝑘_2𝐴(𝑡)𝜌(𝑡)²,(9) where 𝐶 is constant and 𝑘_2𝐴(𝑡) is time-dependent rate parameter. Introducing the average density 𝜌(𝑡) is as in (7) and 𝑘_2𝐴(𝑡) from (8) one obtains the following equation for the intensity of the delayed fluorescence: 𝐼=_𝐷𝐹𝐶(4𝑡𝜋)𝑅_𝐴𝐷_𝑇⎛⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜⎝1+𝑅_𝐴𝑒^{−𝑘𝑇𝑡}−4√𝜋𝐷_𝑇𝑅^2𝐴𝜌₀𝑘_𝑇√𝑡1−𝑒^{−𝑘𝑇𝑡}1+4𝜋𝑅^𝐴𝐷_𝑇𝜌₀𝑘_𝑇1−𝑒^{−𝑘𝑇𝑡}√×⎛⎜⎜⎜⎜⎜⎝𝜌₀𝜋𝐷_𝑇𝑡⎞⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟⎠𝑒^{−𝑘𝑇𝑡}−4√𝜋𝐷_𝑇𝑅^2𝐴𝜌₀𝑘_𝑇√𝑡1−𝑒^{−𝑘𝑇𝑡}1+4𝜋𝑅^𝐴𝐷_𝑇𝜌₀𝑘_𝑇1−𝑒^{−𝑘𝑇𝑡}⎞⎟⎟⎟⎟⎟⎠².(10) The above presented formula can be recognized as an extension of the time-dependent rate parameter 𝑘_2𝐴(𝑡) (2). Instead of exp(−𝑘_𝑇𝑡) in the case of the dominant first-order decay the ratio 𝜌(𝑡)/𝜌₀ is introduced in the denominator of the short-time term. If the second term in the nominator in (7) is negligible: 4√𝜋𝐷_𝑇𝑅^2𝐴𝜌₀𝑘_𝑇√𝑡≈0.(11) Equation (10) can be reduced to the following form: 𝐼=_𝐷𝐹𝐶(4𝑡𝜋)𝑅_𝐴𝐷_𝑇⎛⎜⎜⎜⎜⎜⎜⎝1+𝑅_𝐴𝑒^{−𝑘𝑇𝑡}1+4𝜋𝑅^𝐴𝐷_𝑇𝜌₀𝑘_𝑇1−𝑒^{−𝑘𝑇𝑡}√×⎛⎜⎜⎜⎝𝜌₀𝜋𝐷_𝑇𝑡⎞⎟⎟⎟⎟⎟⎟⎠𝑒^{−𝑘𝑇𝑡}1+4𝜋𝑅^𝐴𝐷_𝑇𝜌₀𝑘_𝑇1−𝑒^{−𝑘𝑇𝑡}⎞⎟⎟⎟⎠².(12) If the second term in the denominator in (7) is also negligible: 4⋅𝜋𝐷_𝑇𝑅_𝐴𝜌₀𝑘_𝑇≈0.(13) Equation (10) is reduced to the form discussed in the paper [22]: 𝐼_𝐷𝐹(𝑡)=𝐶4𝜋𝑅_𝐴𝐷_𝑇1+𝑒^{−𝑘𝑇𝑡}𝑅√^𝐴𝜋𝐷_𝑇𝑡𝑒^{−2𝑘𝑇𝑡}.(14)Non-Fickian Effects in the Nonexponential Decay. In order to compare the intuitive model [22] with the non-Fickian one [19] we decided to introduce another type of the modification of 𝐼_𝐷𝐹(𝑡). The starting points are (5) with 𝑘_𝑇𝑇=𝑘_2𝐴(𝑡) and (9). Now we assume that the nonexponential behavior influence the average density 𝜌(𝑡), but in turn the local density can be reasonably described with the exponential contribution of the first-order decay to the redistribution process. It means that TTA can contribute to the behavior of the average concentration of the molecules in the 𝑇₁ state. At the other side the correction to the average density is so small that the time-dependent rate parameter 𝑘_2𝐴(𝑡) can be reasonably described for the dominant first-order decay. This kind of approximation can be treated as a second step after that presented in papers [19] and [22]. In those papers the first-order decay influenced redistribution but the opposite influence was negligible. Now we introduce nonnegligible contribution of the TTA to the average density. This contribution is however so small that it is not recognized in the time-dependent rate parameter 𝑘_2𝐴(𝑡). This lead to the following form of the nonexponential model for the 𝐼_𝐷𝐹(𝑡): 𝐼_𝐷𝐹(𝑡)=𝐶𝑘_2𝐴(𝑡)𝑒^{−𝑘𝑇𝑡}1+𝑘_2𝐴(𝑡)/𝑘_𝑇1−𝑒^{−𝑘𝑇𝑡}²,(15) where 𝑘_2𝐴(𝑡) can be taken from (2)—intuitive version [22] or from paper [19]—non-Fickian version:𝐼_𝐷𝐹(𝑡)=𝐶4𝜋𝐷_𝑇𝑅_𝐴×⎛⎜⎜⎝1+𝑒^{−𝑡𝑇𝑡}√𝑅^𝐴𝜋𝐷_𝑇𝑡𝑒_{−𝑘𝑇𝑡}1+4𝜋𝐷_𝑇𝑅_𝐴𝜌₀/𝑘_𝑇1+𝑅_𝐴/𝑒^{−𝑘𝑇𝑡}√𝜋𝐷_𝑇𝑡1−𝑒^{−2𝑘𝑇𝑡}⎞⎟⎟⎠,²𝐼_𝐷𝐹(𝑡)=𝐶4𝜋𝐷_𝑇𝑅_𝐴⎛⎜⎜⎜⎝1+√2𝑘_𝑇𝑅_𝐴×⎛⎜⎜⎜⎜⎝𝜋𝐷_𝑇1−𝑒^{−2𝑘𝑇𝑡}⎞⎟⎟⎟⎠𝑒_{−𝑘𝑇𝑡}1+4𝜋𝐷_𝑇𝑅_𝐴𝜌₀/𝑘_𝑇1+√2𝑘_𝑇𝑅_𝐴/𝜋𝐷_𝑇1−𝑒^{−2𝑘𝑇𝑡}1−𝑒^{−𝑘𝑇𝑡}⎞⎟⎟⎟⎟⎠².(16)

3. Experimental

The purification of anthracene and the components of mixture solvent, as well as the preparation of the sample were described in previous papers [19, 22].

The experimental setup used for these experiments was similar to that described previously [19, 22]. The only one modification was the excitation source. Instead of cw Ar⁺ laser and two slow choppers pulsed laser was used. The excitation wavelength was equal to 363.3 nm. Two configurations of excitation source were used:(1)dye laser pumped N₂ laser;(2)dye laser (Lambda Physik LPD 3000) pumped excimer laser (Lambda Physik LPX 100).The first option of excitation will be called “short excitation pulse” (the temporal half width of the pulse ∼1 ns), the other one “middle excitation pulse” (the temporal half width of the pulse ∼20 ns). The laser pulses were focused on the first disc of the fast chopper as presented in Figure 1. The measurements of the delayed fluorescence and phosphorescence were done as in the case of excitation with the “long pulse” [19, 22].

Figure 1 

Modification of the previously described [19, 22] experimental setup for excitation with pulsed lasers.

The decays were measured in the temperature range from 132 K to 152 K in the case of middle excitation pulse and in the range 132 K–150 K for the short pulse. The increment between the temperature values was equal to 1 K.

3.1. Evaluation Procedure

In order to make our model mathematically as simple as possible we use here similar procedure as in the previous papers [19, 22]. We do not use the modified (Nickel’s) initial conditions. The reasons are the same as reported previously: these initial conditions affect only the initial part of the decay (about 1 ms) leaving the main part of the decay unchanged. We use terms anti-Smoluchowski time range and Smoluchowski time range in the same meaning as previously:(i)anti-Smoluchowski time range: the part of a decay where the kinetics significantly depends on the type of initial conditions;(ii)Smoluchowski time range: the part of a decay where the kinetics is practically independent of the type of initial conditions.The main attention is focused on the period between the end of the anti-Smoluchowski time range and the final instant of the nonstationary part of the 𝑘_2𝐴(𝑡). In this time range the nonstationary part of the kinetics has a significant contribution to the total kinetics of the sample. This contribution offers the possibility to obtain the parameter like annihilation radius.

The delayed fluorescence was analyzed with the following functions: 𝐼_𝐷𝐹(𝑡)=𝑃₁⎛⎜⎜⎜⎜⎜⎜⎜⎜⎝1+𝑃₂𝑒^{−𝑡/𝑃3}−𝑃²𝑃₄1−𝑒^{−𝑡/𝑃3}√𝑡1+𝑃₄1−𝑒^{−𝑡/𝑃3}√×⎛⎜⎜⎜⎜⎜⎝𝑒_{−𝑡/𝑃3}𝑡⎞⎟⎟⎟⎟⎟⎟⎟⎟⎠−𝑃²𝑃₄(1−𝑒^{−𝑡/𝑃3})√𝑡1+𝑃₄1−𝑒^{−𝑡/𝑃3}⎞⎟⎟⎟⎟⎟⎠²+𝑃₀(17)𝐼_𝐷𝐹(𝑡)=𝑃₁⎛⎜⎜⎜⎜⎝1+𝑃₂𝑒^{−𝑡/𝑃3}1+𝑃₄1−𝑒^{−𝑡/𝑃3}√×𝑒_{−𝑡/𝑃3}𝑡⎞⎟⎟⎟⎟⎠1+𝑃₄1−𝑒^{−𝑡/𝑃3}²+𝑃₀,(18) where 𝑃₀ is the dark current, 𝑃₁ is the amplitude, 𝑃₂=𝑅_𝐴/(𝜋𝐷_𝑇)^1/2 is the parameter that describes short-time effect, 𝑃₃=1/𝑘_𝑇 is the decay time of the molecules via first-order reaction, and 𝑃₄=4𝜋𝑅_𝐴𝐷_𝑇𝜌₀/𝑘_𝑇 is the ratio of original Smoluchowski rate parameter in the case of stationary diffusion 4𝜋𝑅_𝐴𝐷_𝑇 multiplied with the initial concentration of molecules in the 𝑇₁ state 𝜌₀ and divided through the rate parameter of the first-order decay 𝑘_𝑇. 𝑃₄ will be called the indicator of nonexponential behavior. The condition 𝑃₄≤0.01 was used by Nickel and others [18, 27] as a limit of dominant first-order decay in TTA kinetics.

The parameters obtained from nonexponential decay were compared with the parameters obtained from model described in [22]: 𝐼_𝐷𝐹(𝑡)=𝑃₁1+𝑒𝑃^{2−𝑡/𝑃}₃√𝑡𝑒^−2t/𝑃3+𝑃₀.(19)

The functions used to compare the nonexponential versions of the intuitive (20) and the non-Fickian (21) models have the following form: 𝐼_𝐷𝐹=(𝑡)𝑃₁×⎛⎜⎜⎝1+𝑒𝑃^{2−𝑡/𝑃}₃√𝑡𝑒_{−𝑡/𝑃3}1+𝑃₄1+𝑃²/𝑒^{−𝑡/𝑃3}√𝑡1−𝑒^{−𝑡/𝑃3}⎞⎟⎟⎠²+𝑃₀,(20)𝐼=^𝑀𝐷𝐹𝑃(₁𝑡)1+√2/𝜋𝑃₂/√𝑃₃√1−𝑒^−2t/𝑃3×⎛⎜⎜⎝𝑒_{−𝑡/𝑃3}1+𝑃₄1+√2/𝜋𝑃₂/√𝑃₃√1−𝑒^−2t/𝑃31−𝑒^{−𝑡/𝑃3}⎞⎟⎟⎠²+𝑃₀.(21)

As in the previous papers [19, 22] the points from the begining of the decay curve were systematically cut off with the increment of 10 points. The parameters obtained from fitting procedure were plotted as a function of the delay (starting point of evaluation). This procedure was also applied in the case of asymptotic function derived from (17) and (18) by neglecting the short-time term. From asymptotic functions it was obtained that the redistribution term has negligible contribution for the delay longer than the certain value. The time range where the short-time term significant contribution to the kinetics changes with the temperature (viscosity). The reason for such a behavior is the different mobility of the molecules. The value of 10 ms is the upper limit of the time range of particular interest. The results of fitting of asymptotic functions are generally close to the fit of exponential asymptotic function.

Phosphorescence decays were analyzed with the following functions: 𝐼^non-exp._ph(𝑡)=𝑃₁𝑒^{−𝑡/𝑃3}1+𝑃₄1−𝑒^{−𝑡/𝑃3}+𝑃⁰𝐼_exp.ph,(𝑡)=𝑃₁𝑒^{−𝑡/𝑃3}+𝑃₀,(22) where the signs 𝑃₀, 𝑃₁, 𝑃₃, and 𝑃₄ have the same meaning as in the case of the delayed fluorescence. Since the contribution of the short-time effect to the phosphorescence is much smaller in comparison with the delayed fluorescence the simplified version of the expression for the average density 𝜌(𝑡) in the nonexponential case was used. In this formula the term with the product 𝑃₂𝑃₄ is avoided. In order to avoid small inaccuracies coming from short-time effect first few ms were removed from the decay prior to the fit.

4. Results and Discussion

In Figure 2 there is an example of the result of fitting procedure: the decay of the delayed fluorescence of anthracene measured at 146 K, exited with the middle excitation pulse, presented in the plot parameters, obtained for nonexponential fit and are calculated from (17). Since results obtained with application of the simplified function defined by (18) have almost the same values they are omitted. The values of parameters obtained from nonexponential fit, are compared with the result of exponential approximation—(19). The values of kinetic parameters are presented as a function of the delay (𝑡₀). In Figure 2(a) there is presented the first-order rate constant of the molecules in the 𝑇₁ state:(i)squares—calculated from the delayed fluorescence with exponential model (𝑘^df,e_𝑇);(ii)circles—obtained from the delayed fluorescence within nonexponential approach (𝑘^df,ne_𝑇);(iii)solid line—calculated from phosphorescence with exponential model (𝑘^ph,e_𝑇);(iv)dot-dashed line—calculated from phosphorescence with application of nonexponential model (𝑘^ph,ne_𝑇). The dashed lines present the standard deviations of the exponential model: Δ𝑘^ph,e_𝑇 and dash-dot-dotted line marks the standard deviation of the nonexponential model: Δ𝑘^ph,ne_𝑇 in the case of the phosphorescence decay. In the case of exponential fit the calculated standard deviation is so small for the delayed fluorescence that it cannot be reproduced on the plot although the ordinate is expanded: the presented range is between 29 s⁻¹ and 30 s⁻¹. In the case of the exponential fit of the delayed fluorescence and for the delay time larger than about 1 ms the values of the first-order rate constant 𝑘^df,e_𝑇 can be treated as nearly constant. Which the values obtained with nonexponential model for the delay below ~100𝜇s, the 𝑘^df,ne_𝑇 grow up with the increase of 𝑡₀. This course we interpret as anti-Smoluchowski behavior. The values of 𝑘^df,ne_𝑇 are scattered around the constant value for the delay ~100𝜇s<𝑡₀< ~5.5 ms. They are smaller than the average value of 𝑘^df,e_𝑇 calculated for the delay larger than 1 ms. For 𝑡₀ above 5.5 ms the values of the first-order rate constant obtained from nonexponential model tend to the values obtained for exponential fit. The values of the first-order rate constant obtained from evaluation of phosphorescence decays are smaller than that evaluated from the delayed fluorescence. In the range ~100𝜇s–5.5ms the first-order rate constant calculated with different approaches, with exponential and nonexponential models, can be arranged in ascending order 𝑘^ph,ne_𝑇 < 𝑘^ph,e_𝑇 < 𝑘^df,ne_𝑇 < 𝑘^df,e_𝑇. Figure 2(b) presents the value of the indicator of nonexponential behavior 𝑃₄ (circles) as a function of the 𝑡₀. The course of this parameter versus the delay corresponds to the behavior of 𝑘^df,ne_𝑇. For the delay below ~100𝜇s (the first experimental point only) the value decreases with the increase of the 𝑡₀. For the delay from the range ~100𝜇s–5.5ms the values are scattered around the constant value ≈0.025. For larger 𝑡₀ the coefficient 𝑃₄ tends to zero. Two inserts present: the 𝜒² function (right picture) and the parameter 𝑃₂=𝑅_𝐴/(𝜋𝐷_𝑇)^1/2 (left picture) obtained from the temporal behavior of the delayed fluorescence. The values of 𝜒² in the case of nonexponential fit (circles) are equal to 1 with the exception of the first point, whereas the same parameter obtained for exponential fit (squares) decreases with the increase of 𝑡₀ within the range of delay 𝑡₀ below a few ms. For larger delay the values of 𝜒² obtained from both types of model: exponential and nonexponential are equal. The 𝑃₂ parameter calculated from nonexponential model (circles) has the values scattered around constant ≈0.7s^1/2 for 𝑡₀ within the range from ~100 𝜇s to ~5.5 ms. For smaller delay 𝑃₂ increases with the increase of 𝑡₀. For the delay above 5.5 ms the values of 𝑃₂ have large standard deviation Δ𝑃₂ and are spread over large range of values. The values of 𝑃₂ obtained from exponential fit (squares) grow up with the increase of 𝑡₀ up to about 10 ms. The above description of the results can be summarized as follows: at the constant temperature there is the range of 𝑡₀ where 𝑃₂, 𝑘^df,ne_𝑇 and 𝑃₄ obtained from nonexponential fit is scattered around the constant.

(a)

(b)

Figure 2 

The fitted parameters as dependent on the delay (the starting point of evaluation—𝑡0) obtained for the delay of anthracene in cis/trans-DMCH at 146 K for middle excitation pulse. Part (a): first-order rate constant 𝑘𝑇: squares (𝑘df,e𝑇)—delayed fluorescence, exponential fit (21); circles (𝑘df,ne𝑇)—delayed fluorescence, nonexponential fit (17); solid line (𝑘ph,e𝑇)—phosphorescence, exponential fit; dashed lines—accuracy of 𝑘ph,e𝑇; dash-dotted line (𝑘ph,ne𝑇)—phosphorescence, nonexponential fit; dash-dot-dotted lines—the accuracy of (𝑘ph,ne𝑇). Part (b): main plot indicator of nonexponential (𝑃4) behavior; right inset: comparison of smoluchowski parameter   (𝑃2) for nonexponential fit (lower curve) and exponential fit (upper curve); left inset: autocorrelation function for nonexponential fit (lower curve) and exponential fit (upper curve).

The parameters obtained for the excitation with the “short pulse” have very similar dependence on the delay 𝑡₀ in comparison with data obtained from the middle excitation pulse. The example of experimental data (measured at 150 K) analysis is described below. The values of 𝑘^df,e_𝑇 can be treated as a nearly constant for the delay 𝑡₀ larger than about 1 ms. For 𝑡₀<∼1 ms; 𝑘^df,e_𝑇 decreases with the increase of 𝑡₀. For 𝑡₀<0.5 ms the increase of the value of the first-order rate constant 𝑘^df,ne_𝑇 with the increase of the delay is observed. 𝑘^df,ne_𝑇 can be treated as constant for 𝑡₀ from the range from ~0.5 ms to ~5 ms. In the case of the delay longer than ~5 ms the values of 𝑘^df,ne_𝑇 are much more scattered in comparison with that from the range ~0.5 ms–~5 ms. Also the standard deviation Δ𝑘^df,ne_𝑇 increases with the increase of the delay. This increase is significant for 𝑡₀>∼5 ms. This behavior reflects the decreasing contribution of the short-time effect to the kinetics of TTA with the increase of 𝑡₀. In the range of 𝑡₀ where 𝑘^df,ne_𝑇 can be treated as constant the values of the first-order rate constant obtained from different measurements (delayed fluorescence and phosphorescence) and different models (exponential and nonexponential) can be placed in the similar row as for excitation with the “middle pulse”: 𝑘^ph,ne_𝑇≤𝑘^df,ne_𝑇<𝑘^ph,e_𝑇<𝑘^df,e_𝑇. The course of indicator of nonexponential behavior 𝑃₄ corresponds to that of 𝑘^df,ne_𝑇. For the delay below ~0.5 ms the value of the parameter decreases with the increase of 𝑡₀. In the range from ~0.5 ms to ~4.5 ms–~5 ms the value can be treated as constant within the accuracy of the measurement and the fitting procedure. For larger values of the delay 𝑃₄ is scattered and its accuracy Δ𝑃₄ increases with the increase of the 𝑡₀. The values of 𝑃₂ obtained from exponential model increase with the delay up to 𝑡₀≈8 ms. In the case of the nonexponential fit the value of 𝑃₂ can be treated as constant for the delay from the range ~0.5 ms–~5 ms. Outside this range the parameter 𝑃₂ has different course:(i)it increases with the increase of 𝑡₀, for the 𝑡₀<∼0.5 ms;(ii)in the case of 𝑡₀>∼0.5 ms its scattering and standard deviation grows up with the delay.For 𝑡₀ below 4.5 ms–5 ms the values of 𝜒² calculated from exponential model are larger than the values obtained for nonexponential fit. For the longer delays (𝑡₀>∼0.5 ms) the values of 𝜒² are equal for both models.

Since the kinetics of delayed fluorescence depends on the mobility of molecules in 𝑇₁ state the differences between exponential and nonexponential models of TTA should be much larger in the case of delayed fluorescence than in the case of phosphorescence. Generally the temperature range of the study can be divided into two parts:(i)in the first one the exponential model gives reasonable description of the kinetics of TTA;(ii)in the other one the application of nonexponential model is necessary in order to obtain reasonable values of the kinetic parameters.The temperature where the kinetic model has to be change from exponential to nonexponential will be called transition temperature.

In Figure 3 there is presented the first-order rate constant of the molecules in the 𝑇₁ state as dependent on temperature. Squares represent the data obtained from the fitting procedure of the delayed fluorescence with application of the exponential model—𝑘^df,e_𝑇. Stars represent the result of the fitting of the exponential model to phosphorescence decays—𝑘^ph,e_𝑇. The values of the first-order rate constant obtained from nonexponential model are represented by: circles—delayed fluorescence (𝑘^df,ne_𝑇) and triangles—phosphorescence (𝑘^ph,ne_𝑇). Part (a) presents the data measured when molecules are excited with the middle pulse, part (b) for excitation with the short one. The ordinate of both parts of the figure is expanded: only the range between 28.6 s⁻¹ and 30.8 s⁻¹ is shown. This is done in order to present small changes of the first-order rate constant obtained from different approaches. The behavior of the parameter in both parts of Figure 3 is similar.

(a)

(b)

Figure 3 

First-order rate constant as a function of temperature. Part (a): measurements performed with middle excitation pulse (dye laser pumped excimer laser), part (b) measurements in the case of excitation with short pulse (dye laser pumped nitrogen laser). Results of following measurements and fits are compared: exponential fit of the delayed fluorescence (squares), nonexponential fit of the delayed fluorescence (circles), exponential fit of the phosphorescence (crosses), and nonexponential fit of the phosphorescence (triangles). All fits were performed for intuitively modified model.

The values of the first-order rate constant obtained from phosphorescence with application of both models are not very different from each other. Generally, the values of 𝑘^ph,ne_𝑇 are only slightly smaller in comparison with those of 𝑘^ph,e_𝑇. There is one exception. The value of 𝑘^ph,ne_𝑇 obtained from fit of the data measured at 146 K and with short excitation pulse is placed significantly below the course of the first-order rate constant obtained for other temperatures (Figure 3(b)). Also the accuracy of 𝑘^ph,ne_𝑇 at 146 K is significantly larger in comparison with the values obtained for measurements performed at other temperatures. It suggests that in these conditions the fit of the nonexponential model to the phosphorescence decay gives larger uncertainty in comparison with other temperatures. Since this temperature can be treated as a transition temperature between the ranges of applicability of both models (exponential and nonexponential) the uncertainty may result from the large number of parameters in the nonexponential model. Parameters in the nonexponential model like indicator of nonexponential behavior has small values and standard deviation almost equal to 100% of the value (see Figure 4 and its discussion). Small values and large standard deviation of this parameter influence the values of other parameters.

(a)

(b)

Figure 4 

Indicator of nonexponential behavior as dependent on the temperature measured for excitation with the middle pulse—part (a) and short pulse—part (b). Comparison of the data obtained from the fit of the delayed fluorescence (circles) and phosphorescence (triangles).

Since the middle excitation pulse has significantly larger energy than the short one the transition temperature should be different for different types of excitation pulses. The temperature range of the applicability of exponential model is defined as a temperature region where the values of the first-order rate constant calculated from delayed fluorescence and phosphorescence build two almost parallel lines [22]. This range is as follows:(i)138 K–142 K for middle excitation pulse;(ii)138 K–146 K in the case of the short excitation pulse.

The temperature dependence of the 𝑘^df,e_𝑇 below the transition temperature is the same as that obtained for intuitively modified model [22]. Below 138 K the differences between 𝑘^df,e_𝑇 and 𝑘^ph_𝑇 increase with the decreasing of the temperature. The reasons for this behavior were suggested in the previous papers [19, 22]. For the temperature values above the transition temperature the difference between 𝑘^df,e_𝑇 and 𝑘^ph_𝑇 increases with the increase of the temperature. The temperature dependence of the 𝑘^df,ne_𝑇 shows the advantage of the nonexponential model in the temperature range above the transition temperature. In this temperature range the values of the 𝑘^df,ne_𝑇 are much closer to 𝑘^ph_𝑇 than the values of 𝑘^df,e_𝑇. For the middle excitation pulse the values of 𝑘^df,ne_𝑇 are placed between the values of 𝑘^df,e_𝑇 and 𝑘^ph_𝑇 in the temperature range 143 K–145 K. For the values of the temperature higher than 145 K the values of 𝑘^df,ne_𝑇 are placed very close to the 𝑘^ph,ne_𝑇.

The values of 𝑘^df,ne_𝑇 obtained for the temperature below the transition temperature have much larger standard deviation (Δ𝑘^df,ne_𝑇) in comparison with the values obtained from decays measured above the transition temperature. Also the increase of Δ𝑘^df,ne_𝑇 with the decrease of the temperature is observed within the range below the transition temperature. This suggests that the nonexponential model may produce some artificial effects due to too large number of parameters if it is applied for the experimental conditions where the nonexponential effects are negligible. In this case the parameter describing nonexponential behavior has small value and the uncertainty over 100% influences the values of other parameters and their standard deviations.

In Figure 4, there is presented the indicator of nonexponential behavior—𝑃₄. The data for middle excitation pulse are shown in part (a), for short excitation pulse—in part (b). The circles represent the data obtained from fitting procedure applied to the decays of delayed fluorescence, triangles-phosphorescence. The temperature dependence of the parameter 𝑃₄ confirms the conclusion that the nonexponential model offers better description of the kinetics for the temperatures higher than the transition temperature. Starting from 146 K the value of 𝑃^df₄ increases with the increase of the temperature. In the case of the middle excitation pulse (Figure 4(a)) the values of the indicator of nonexponential behavior calculated from delayed fluorescence (𝑃^df₄) increase with the increase of the temperature within the temperature range 142 K–149 K and reaches the value of ~0.06 at 149 K. In the temperature range 149 K–152 K the values of 𝑃^df₄ are scattered around the constant. Below the transition temperature the values of 𝑃^df₄ are scattered around zero and the standard deviation Δ𝑃^df₄ is generally larger than |𝑃^df₄|. The data obtained from the excitation with the short pulse show very similar behavior. Since the energy of short pulses was much smaller the transition temperature is higher and equals 146 K. Also the range of changes of 𝑃^df₄ is smaller in comparison with the data obtained for middle excitation pulse. The values of 𝑃^df₄ increase with the increase of the temperature within the temperature range from 146 K to 150 K. The maximum of 𝑃^df₄ equals to ~0.02 at 150 K. Below the transition temperature the behavior of the indicator of nonexponential behavior is the same as for the middle excitation pulse. The values of the indicator of nonexponential behavior calculated from the phosphorescence (𝑃^ph₄) are generally slightly larger than zero in the whole temperature range of investigation. The tendency that the values of 𝑃₄ grow up with the increase of the temperature in the range above the transition temperature is observed only for the delayed fluorescence. The reason for this behavior is that the contribution of the diffusion-controlled process to the phosphorescence decay is much smaller than to the delayed fluorescence. In turn: the contribution of nonexponential effect should be here much smaller in comparison with the delayed fluorescence.

Comparison of the kinetic parameters obtained from intuitive and non-Fickian models is presented in Figures 5 and 6. In each figure part (a) presents the data obtained for intuitive treatment, part (b)—for non-Fickian model. In Figure 5 the first-order rate constant is shown as a function of temperature. The data were obtained for middle excitation pulse. In the case of application of the intuitive model to the experimental data measured for middle excitation pulse the kinetics of diffusion-controlled TTA is well described within the exponential approximation for the temperature below 142 K. For the temperature between142 K (transition temperature) and 149 K the values of the first-order rate constant obtained from (20) have the same course as presented in Figure 3(a). Above 149 K the approximation used in (20) do not give reasonable fit of the model to the experimental decay. In Figure 5(a) there are also presented the values of the first-order rate constant obtained from the fit of phosphorescence. The comparison of 𝑘^df,ne_𝑇 obtained from (20) with 𝑘^df,e_𝑇 and 𝑘^ph_𝑇, shows that simplified nonexponential model defined by (20) gives the same results as model defined by (17).

(a)

(b)

Figure 5 

Comparison of the values of the first-order rate constant calculated with intuitive (a) and non-Fickian (b) models. The values of the parameter are presented as dependent on temperature. The measurements were performed for so-called middle excitation pulse.

(a)

(b)

Figure 6 

Comparison of the values of the indicator of nonexponential behavior calculated with intuitive (a) and non-Fickian (b) models. The values of the parameter are presented as dependent on temperature. The measurements were performed for socalled middle excitation pulse.

In Figure 5(b) the temperature dependence of the first-order rate constant calculated from the nonexponential version of non-Fickian model is compared with the data obtained for this model within exponential approximation. The nonexponential modification introduces no significant correction to the values of the first-order rate constant calculated from phosphorescence decays below 146 K. For higher temperatures the values of 𝑘^ph,ne_𝑇 have the tendency to be smaller than 𝑘^ph,e_𝑇. The values of the first-order rate constant calculated from delayed fluorescence within non-Fickian exponential model are almost equal to that calculated from phosphorescence in the temperature range below 138 K. For higher temperatures the difference between 𝑘^df,e_𝑇 and 𝑘^ph,e_𝑇 is observed: 𝑘^df,e_𝑇 > 𝑘^ph,e_𝑇. This difference increases with the increase of the temperature. Application of nonexponential model extends the range of coincidence of values of 𝑘^df,ne_𝑇 and 𝑘^ph,e_𝑇 up to about 146 K. For higher temperatures the values of 𝑘^df,ne_𝑇 are smaller than 𝑘^ph,e_𝑇. This difference increases with the increase of the temperature.

The behavior of the first-order rate constant for short excitation pulse is similar to that presented for middle excitation pulse. The differences between exponential and nonexponential models are smaller because the energy of excitation pulse was smaller in comparison with the energy of middle excitation pulse. The lower energy of excitation pulse results in lower concentration of molecules in the triplet state.

The temperature dependence of the indicator of nonexponential behavior (𝑃₄) for middle excitation pulse is presented in Figure 6. As in the case of the first-order rate constant part (a) of this figure concerns the data obtained from intuitive model, part (b)—from non-Fickian model. The values of the parameter 𝑃₄ calculated within intuitive model from the delayed fluorescence are equal to zero within the accuracy of the measurement and the evaluation procedure for the temperature below 142 K. For the temperature above 142 K the values of the indicator of nonexponential behavior increase with the increase of the temperature up to about 0.06 for 149 K. In the case of the temperature values larger than 149 K the fit of (21) to the decays of the delayed fluorescence do not give stable result. The values calculated from phosphorescence grow up with the temperature in the range from 145 K to 149 K. In the temperature below 145 K the values of 𝑃₄ (calculated from phosphorescence) are scattered between 0 and 0.01. In the case of non-Fickian model (Figure 6(b)) the values of 𝑃₄ calculated from delayed fluorescence are scattered around 0.01 for the temperature below 142 K. Above this temperature the values of the indicator of nonexponential behavior increase with the increase of temperature. The largest value of 𝑃₄ equals to about 0.06 for 149 K. The behavior of the indicator of nonexponential behavior observed for short excitation pulse is qualitatively the same as for middle pulse. Here however the transition temperature equals to 145 K instead of 142 K. Also the largest value of 𝑃₄ reached at 149 K equals to 0.02 instead of 0.06 reached at the same temperature for middle excitation pulse. The reason for this difference is, as mentioned above, the lower energy of the exciting pulse. Pulses having lower excitation energy generate lower concentration of the molecules in the 𝑇₁ state. This implies that nonexponential behavior gives measurable component to the decay in the system with higher mobility of the molecules, that means for the conditions where the diffusion coefficient is larger.

To sum up, the introducing of the nonexponential model offers better description of the kinetics within this temperature range in which the combination of the energy of excitation pulse and the mobility of molecules results in significant deviation from the exponential model. Both models: the intuitive and non-Fickian has their limits. The intuitive model has the tendency to overestimate the first-order rate constant calculated from the delayed fluorescence in comparison with the value obtained from the phosphorescence. The course of the first-order rate constant calculated from the delayed fluorescence with non-Fickian model does not show this systematic difference. At the other side the non-Fickian model gives a good fit to the experimental decay within the smaller range of the temperature.

The indicator of nonexponential behavior obtained from intuitive model equals to zero within the temperature range where the exponential model offers good description. In the case of non-Fickian model this parameter tends to 0.01 instead of zero with the standard deviation significantly smaller than 100% in the temperature range where the exponential approach offers good description of the kinetics.

The most important conclusion coming out from the paper is that the nonexponential approach significantly enhanced the range of experimental conditions where the kinetics of diffusion-controlled TTA has reasonable description. The indicator of nonexponential behavior—parameter 𝑃₄ should be <0.01 if the exponential approach provides reasonable description of the kinetics. The nonexponential model allowed the use of the model up to values of 𝑃₄ about six times larger (𝑃^4max≈0.06).

The annihilation radius can be calculated from equation for 𝑃₂=𝑅_A/(𝜋𝐷_𝑇)^1/2 parameter. This evaluation requires the value of mutual diffusion coefficient 𝐷_𝑇. The values of diffusion coefficient measured in other solvents and extrapolated to experimental conditions used here give reasonable values of annihilation radius 𝑅_𝐴 at least in the temperature 136 K–149 K [22]. However, the values of 𝐷_𝑇 extrapolated from other solvents do not take into account the specific features of the solvent mixture used in this work. Sophisticated proof of the model requires the values diffusion coefficient measured at the same experimental conditions as presented here for kinetic study. Values of 𝐷_𝑇 can be calculated from decays of delayed fluorescence obtained from measurements with spatially periodic excitation. This kind of excitation requires in turn further extension of the nonexponential approach. The introducing of the concept of joined action of the diffusion and the first-order decay to the triplet-triplet annihilation kinetics in the case of spatially periodic excitation will be the subject of subsequent work.

Acknowledgments

The authors thank Professor Jürgen Troe (Max-Planck Institute, Göttingen, Germany) for generous support of the research project. The discussion with Professor K. Rotkiewicz (Institute of Physical Chemistry PAS, Warsaw, Poland) was very fruitful for preparation of this manuscript. The support of stay in Göttingen (Germany) of one of the authors (P. B.) by Max-Planck-Gesellschaft (from 1.10.2000 to 30.09.2001) and by Alexander von Humboldt-Foundation (from 1.10.2001 to 31.12.2002) is gratefully acknowledged. Dr B. Nickel passed away 27/01/2002.