1. Introduction
The viscosity of a solvent (η) is the dominant factor that affects rates of fast bimolecular and some monomolecular reactions. Bimolecular fast reactions include the termination of free radicals, some swift radical addition reactions, and proton/electron/hydrogen atom transfer. The mini-review will focus on the self-termination of reactive free radicals:
For brevity, Reaction (1) is named recombination. (In general, self-termination of radicals may be disproportionation or competition of both recombination and disproportionation).
Some rearrangements and cis–trans isomerization are examples of fast monomolecular reactions.
The viscous drug of a solvent limits the rates of the listed above chemical transformations. When mutual diffusion of reagents controls the reaction, its rate constant is expected to equal the renowned von Smoluchowski equation [
1]:
where ρ is the sum of Van der Waals radii of reagents, and D
12 is a mutual diffusion coefficient D
1 + D
2. The symbol D is used for the diffusion coefficient of an individual molecule or radical. Formula (2) should be modified in the case of diffusion-controlled reactions between free radicals [
2]:
where σ = 1/4 is a spin statistic factor. For a diffusion-controlled reaction between identical radicals (1), 2
k should equal
kdiff. Formula (3) requires further correction for relatively large radicals (molecules) with a localized highly reactive atom. This highly reactive atom is an atom bearing the unpaired electron for radicals. The formula is [
2]:
where f
eff is a steric factor. Such diffusion-controlled reactions with f
eff < 1 are termed
pseudodiffusion. See for details [
2,
3]. The “usual”
diffusion-controlled reactions have f
eff = 1.0, Equation (3). Bimolecular reactions in the liquid phase are usually divided into two classes- activation-controlled and diffusion-controlled [
1]. In this paper, we attract attention to the less common class of biomolecular fast reactions—
pseudodiffusion.
The substitution of D in Equation (2) is expressed by the famous Stokes-Einstein formula
for translational diffusion leads to the well-known Debye formula [
1]:
k
B is the Boltzmann constant, R—universal gas constant, r—Van der Waals radius of a diffusing molecule (radical), η—solvent viscosity in Equations (5) and (6).
Modification of the Equation (4) for the diffusion-controlled reaction between two radicals is:
The main focus of this work is the recombination of reactive organic free radicals and their dependence upon solvent viscosity.
Unless stated otherwise, we present below data obtained at 293–298 K and atmospheric pressure.
2. Recombination of Radicals in the Newtonian Liquids
Newtonian liquids are the main subject of this paper, where the dynamic η is the only characteristic of a liquid at the low and high shear rate [
4]. (η depends upon temperature and external pressure). The values of η for many Newtonian individual liquids and liquid mixtures can be found in reference books and in the literature on chemistry and physics. Many η values have been measured with high accuracy.
It follows from Equations (6) and (7), that the rate constant of a molecular mobility-controlled reaction is inversely proportional to η at the constant T. (Say, water has η = 1 cP whereas glycerol has η ≈ 1000 cP. The fast bimolecular reaction would proceed slower in glycerol compared to water). There are other additional criteria of a (pseudo)diffusion recombination related to activation energy and activation volume of the reaction and the viscous flow [
2,
3].
The authors of [
3] studied the recombination of highly reactive radicals labeled I–III on
Scheme 1 below:
Figure 1 below displays the dependence of 2
k1 for these three radicals upon η:
Changing the viscosity in an experiment
k = F(η) by changing solvents of different chemical nature leads to a large data scatter. (It is a good practice to vary η using alkanes or alcohols of different viscosities. One can start with pentane η = 0.24 cP and end up with heptadecane (C
17), the highest MW alkane, a liquid with η = 4.0 cP. The least viscous alcohol is methanol η = 0.55 cP, and undecan-1-ol (C
11) is the highest MW linear alcohol, a liquid with η = 17 cP.) Binary mixtures of non-viscous and viscous solvents of similar chemical nature are suitable for measuring the viscosity effect. It was found that the binary mixture toluene—dibutyl phthalate 0.6–18.4 cP or toluene—vaseline oil (η = 120 cP) are quality binary mixtures. Measuring η solvents or solvent mixtures with η ≥ 10 cP in η—dependence experiments is recommended). The solvent’s nature may have a pronounced effect on reaction (1), and solvents (S) can form solvation complexes R
.S affecting the rates of recombination of free radicals [
5].
However, in many experiments at constant T for (pseudo)diffusion-controlled reactions, Equation (8) holds:
The latter is a linear plot convenient for visualization; see, e.g.,
Figure 1. In this case, α ≈ 1.0 at η > 2 cP. However, in many experiments α = 0.5–1.2; more often, α ≤ 1.0 [
3,
6].
It would be naïve to expect the precise prediction of the corresponding experimental values by Equation (5) and the consequential Equations (6) and (7).
The derivation of Equation (5) assumes that the molecule/radical is spherical and diffuses in the structureless continuum. A molecule is not necessarily ball-like and diffuses in the structured media among low MW molecules of comparable size. (Probably only fullerene C
60 is the ball-like molecule). These facts lead to the deviation of experimental D from its predicted value by Equation (5). Many attempts have been made to introduce corrections into Equation (5), a “microfriction” coefficient f
micro,
which leads to an agreement between the calculated and experimental values [
7,
8]. There are speculations that a prolate or an oblate form of a molecule, slip or stick boundary condition affects the value of D [
7,
8]. f
micro accounts for radii or volumes of reactive and solvent molecules [
7,
8]. A simple Spernol-Wirtz formula for f
micro has often been used in the correction of D since 1953 [
7]:
where r is a reagent’s Van der Waals radius (Equation (5)), and r
L is such a radius for a solvent molecule. For radicals I–III (
Scheme 1) in toluene f
micro ≈ 0.63 per the estimation of [
3]. Viscosity, which should be placed in Equation (5) leading to the experimental value of D, is often called microviscosity η
micro. If a shear viscosity of the media (a Newtonian or non-Newtonian liquid) is η, then:
As a reminder, shear viscosity is a rheological property that describes a fluid’s resistance to flow under shear stress [
4].
The diffusion of elongated molecules was analyzed in [
9].
3. Experimental Measurement of Diffusion Coefficients or Microviscosity
It would be better to measure the D value experimentally instead of estimating D with Equations (5) and (9). It is difficult for transient species and requires much effort, even for stable molecules. Terazima et al. [
10,
11] developed a transient grating (TG) technique that allows measurements of D of short-live radicals in solution D(R.). In the TG method, two coherent light beams are crossed at a sample to produce the interference pattern of the light intensity. Analysis of the grating pattern eventually allows for D(R.). Several obtained D(R.) were reported [
10]. It was interesting to compare D(R.) with D of the parent hydrogenated molecules D(RH). It was observed that in several cases, D(R.) is two to three times smaller than D(RH). The reasons for such differences need to be better understood [
9]. Unfortunately, this promising work has not been continued yet.
Reversible recombination of 2,6-diphenyl-4-methoxyphenoxyl radical is controlled by diffusion [
12]. It was possible to measure D(R.) in benzene by a capillary method for this system radicals-dimer [
12]. The benzene solution of radicals-dimer was put into a thin capillary of the known dimensions, and the opened capillary with reagents was placed into a thermostatted large vessel filled with the same solvent. The concentration of [R.] in the capillary, which decreased with time (
t), coincides with periodically measured ESR. A mathematical analysis of [R.] vs.
t allowed us to obtain D(R.) [
12,
13]. The experimentally measured 2
k1 = (1.5 ± 0.3) × 10
9 M
−1·s
−1 with the experimental error with a value predicted by Equation (3)
kdiff = (1.1 ± 0.2) × 10
9 M
−1·s
−1. The reaction radius was estimated at ρ = 1.0 ± 0.1 nm [
12].
Despite a few observed differences between D(R.) and D(RH) observed in [
10], D(RH), or even D of another stable molecule of comparable dimensions, many research publications use the obtained D as a reasonable estimation of D of the reactive radical/molecule or η
micro. η
micro can be measured by fluorescence technique with molecular rotors [
14], rotational diffusion correlation times by NMR, diffusion-ordered spectroscopy (DOSY) by NMR [
15], and by ESR with nitroxyl radicals as spin probes [
13]. The reader can find different versions of fluorescence, NMR, and ESR techniques used for the measurements of D in the relevant publications.
Rate constants 2
k1 for recombination of 2,6-di-
tert-butylphenoxyl and 2,6-di-
tert-butyl-4-methylphenoxyl radicals were measured in different solvents by flash photolysis [
16]. Kinetics ESR study reproduced these 2
k1-values [
16]. The abovementioned phenoxyl radicals participate in pseudodifffusion recombination (1). It was concluded that
tert-butyl substituted radicals, contrary to phenyl-substituted phenoxyls, do not form solvate complexes [
5]. Phenyl-substituted radicals are not sterically hindered and have a large π-electron system. They can form solvation complexes R
S. with some solvents.
tert-Butyl-substituted radicals are sterically hindered by
tert-butyl groups around oxygen atoms and are less prone to form R
S [
5]. Rotational correlation times τ for the parent two phenols in the same solvents were measured by NMR [
16]. It was determined that 2
k1 ∝ 1/τ or 2
k1·τ = const for these two
tert-butyl phenoxyls. Thus, the molecular mobility of the parent phenols nicely reflects the mobility of the corresponding phenoxyls [
16].
The kinetic ESR study in different solvents studied recombination (1) of benzyl radicals and substituted benzyls [
17]. The reactions are diffusion-controlled [
17]. D of toluene and substituted toluene were measured in the same solvents. It was found that 2
k1 coincides with the experimental error with corresponding values calculated by Equation (3) [
17]. Thus, D(RH) is usually a good estimation of D(R.).
4. Cage Effect in Different Solvents
The cage effect (φ) plays a vital role in the photoinitiation of chemical reactions and free-radical chemistry [
18,
19,
20,
21].
φ is the fraction of radicals that decayed near each other after photolysis or thermolysis of molecules in solution; 0 ≤
φ ≤ 1. It is known that
φ increases with the solvent viscosity η increase [
18,
19,
20,
21]. Tyler and coworkers [
18,
19] studied the dependence of experimentally measured
φ vs. η under photolysis of an organometallic dimer. Solvents of different chemical nature, including a solution of a polymer, were used. Dependences
φ vs. η demonstrated much scatter. However, the authors measured the microviscosity η
micro of a stable molecule of similar structure to the studied organometallic radical by NMR [
18,
19]. Interestingly, all dependencies of
φ vs. η
micro lay on the one smooth and increased with viscosity as expected [
19]:
Figure 2 is an awe-inspiring result confirming the expectation related to the
φ and its microviscosity dependence [
15,
18,
19]. The work [
19] aimed to find a method to quantitatively correlate
φ values to the properties of the solvent system. Although bulk viscosity qualitatively correlates with
φ cage, the shortcomings of using solvent η for quantitative work are well-known. The new method developed in [
19] employs η
micro, a parameter straightforwardly measured using an NMR probe. As demonstrated herein, microviscosity provides a reliable correlation between the cage recombination efficiency in all categories of solvents, thus lending itself to quantitative predictions. The correlation between
φ and microviscosity holds for a wide range of solvent systems, spanning nonpolar, polar, aromatic, and hydrogen-bonding solvents [
19]. In addition, selective solvation from mixed-type solvent systems was shown to insignificantly affect the predictive power of the method [
19]. It is important to emphasize the straightforwardness of the method: to predict
φ of a particular couple of radicals in a particular solvent, all that is required is to measure the diffusion coefficient of an appropriate NMR probe [
19].
It is expected that 1/
φ = 1 + const/microviscosity is to be a linear plot [
20,
21].
The kinetics of free radicals recombination in the solvent cage and recombination in the solvent bulk are analyzed in [
22].
5. Microviscosity of Complex Liquids
A solvent or a reaction media is not necessarily a Newtonian liquid. Much effort was devoted to studying diffusion in micelles, biological objects (vesicles, membranes), and liquid polymers or polymer solutions. These liquids are usually viscous, their shear η does not reflect a viscous drag imposed on low MW molecules or radicals.
Such media viscosity is measured with one or another viscometer. In particular, the polymer and coatings industry measures η with Brookfield viscometer under specified conditions. That measurement can be replicated at another laboratory with a standard Brookfield.
In non-Newtonian complex liquids, a notion of microviscosity (ηmicro) is often used. ηmicro is measured with a molecular or a radical probe; see above.
The local viscosity or η
micro supposedly plays an important role in biology across several length scales [
23]. There is growing evidence that the fluidity of the cell’s lipid plasma membrane is regulated to enable optimal cellular function, and η
micro within the cell’s nucleus may influence the formation of non-canonical DNA assemblies [
23]. However, unlike in the cases of bulk industrial applications (see above), including additive manufacturing, biological viscosity is challenging to measure and requires the development of methodologies that can probe viscoelastic properties of microscopic volumes of individual biological cells and organelles [
23]. While a myriad of different dyes have been developed to study the microenvironment of lipid membranes, Molecular fluorescent rotors give valuable η
micro of their surroundings in a biological object [
23]. At the same time, many developed dyes-probes of η
micro often respond to multiple biophysical features and/or are not fully calibrated against the various membrane parameters [
23]. This is a problem for biologists and biophysicists. Microviscosity can indicate the condition of living cells [
24]. It is closely linked to numerous diseases [
23,
24]. It is necessary to design tools to effectively monitor η
micro changes, which could provide promising avenues for treating diseases. A novel mitochondria-targeting fluorescent probe
P was suggested in [
25] for the detection of η
micro changes in vivo and in vitro.
P has advantages such as long emission wavelength (650 nm), large Stokes shift (105 nm), significant fluorescence enhancement (59-fold), high sensitivity, good biocompatibility, and so on. Biological experiments showed that
P could target mitochondria and detect viscosity alterations in HeLa cells. (HeLa cells are the first immortal human cell line [
24]). Moreover,
P allows the exploration of fluctuations in η
micro within living organisms [
24].
The authors [
25] aimed to measure microviscosity in dimyristoyllecithin vesicles with a fluorescent probe. The system turns out to be rather complex [
25]. At the same time, nanosecond time-dependent depolarization measurements of fluorescence are capable of reporting in detail the complex rotational motion of fluorescent probes in biological membranes thereby providing information on both their η
micro and structure [
25].
The viscosity of cell membranes is a crucial parameter that affects the diffusion of small molecules both across and within the lipidic membrane [
26]. It was mentioned above that such diffusion can be related to several diseases [
26]. The possibility to measure quantitatively membrane viscosity η
micro on the nanoscale is of great interest [
26]. The paper [
26] reports a complete investigation of the photophysics of an amphiphilic membrane-targeted azobenzene, and the results validate its use as a viscosity probe for some other cell membranes.
Trans-cis photoisomerization azobenzene was a method to develop a molecular viscometer and to assess the viscosity of
E. coli bacteria membranes employing time-resolved fluorescence spectroscopy. Lifetime measurements of the probe in
E. coli bacteria suspensions correctly indicate membrane η
micro value in live
E. coli cells [
26]. η
micro changes from 10 to 5 cP when increasing the temperature from 295 K up to 313 K [
26].
Microviscositry of the studied mitochondria measured with a fluorescent probe varies from 63 to 130 cP [
27]. Supposedly, it gives valuable information on the subcellular microenvironment [
27].
According to [
28], “how intracellular contents affect the mobility of single molecules is a key aspect of cell biology”. A molecular probe gives such info for solutions of calmodulin labeled with green fluorescent protein as a diffusing probe and containing polymer dextran with MW 10 to 500 kg/mol [
28]. Not surprisingly, microviscosity essentially differs from the shear η and depends upon the concentration and MW of a polymer [
28].
Near-infrared emissive indolizine squaraine fluorophores turn out to be efficient molecular viscosity sensors of intracellular microenvironments [
29]. Moreover, the quantum yield of fluorescence of these fluorophores strongly depends upon viscosity [
29].
Below, we list a few examples not related to biology.
Typical micelles sodium dodecyl sulfate (SDS), cetyltrimethylammonium bromide (CTAB), and some others) have η
micro in the range of 12–45 cP measured with an ESR probe [
30]. Microviscosities of the several studied micelles obtained by the ESR method are in concordance with those obtained by other techniques [
30].
The processing of initially liquid materials often involves a drying step, during which a solvent evaporates from a multicomponent mixture. Such materials are usually confined in a thin film or a droplet [
31]. Examples of such materials may be paints or coatings. Evaporation leads to the concentration of non-volatile material in the product. The observed η of the material measured one or another way cannot reflect the occurring changes in different parts of the drying material. The authors of [
31] found a way to measure microviscosity for local and time-dependent mapping of concentration fields in such confined systems. The use of fluorescence lifetime imaging microscopy and fluorescent molecular rotors is suggested as a versatile quantitative method to map η
micro. Measured microviscosity and concentration fields allow the characterization of drying time dependence in agreement with simple modeling of the evaporation [
31].
High MW commercial poly(dimethylsiloxane) (PDMS) with added reagents has η ~ 10
3 P [
32]. However, reactive radicals recombine in such media very fast like it had η ~ 1 cP. Thus, microviscosity is ~10
5 times lower than the shear viscosity [
32].
Knowledge of microviscosity is important in polymer-based lubrication systems [
33]. A relationship between the change in microviscosity and the friction coefficient is established in [
33]. η
micro was measured with molecular rotors [
33]. In contrast to viscosity η, η
micro enables a precise characterization of each component’s diffusion and adsorption rates within the contact region at micro and nano scales. Therefore, η
micro exhibits a more rapid response, which is very important for studying polymer aqueous lubrication friction regulation systems and for revealing the friction coefficient’s response to viscosity changes within the contact region and the underlying mutation mechanism. However, the existing mechanical viscosity measurement methods pose challenges when it comes to detecting η
micro at polymer-modified interfaces [
33].
Thus, microviscosity indicates the diffusivity of low MW molecules/radicals in a complex media. At the same time, ηmicro is a characteristic of a specific media and cannot be considered a parameter that could be used for another similar media.
6. Fast Monomolecular Reactions
Almost any rearrangement in the condensed phase experiences the viscous drag of the media to one or another extent. Much effort was devoted to
cis–trans (photo)isomerization in the liquid phase and solvent viscosity effect on this process. The movement of relatively large groups (phenyl and bulkier) behaving like blades is retarded by η [
34,
35,
36,
37,
38,
39]. At the same time, an excellent approximation to
cis–trans isomerization is a movement along
one reaction coordinate, usually the angle between the moving groups. Such a reaction is expected to be described by Kramers’ theory [
38]. There are many problems related to applying Kramers’ theory to isomerization and other monomolecular reactions [
39,
40,
41]. They are beyond the scope of this manuscript. Some authors consider volumes of the rotating group and a solvent and f
micro [
36,
37]. One of the appealing to the experimenter applications of Kramers’ theory is presented by Formula (12) [
36]:
where
kiso/s
−1 is a (photo)isomerization rate constant, 0 < α < 1.0 is a constant parameter for an experiment, and E
o (kJ/mol) is the height of the isoviscosity barrier. E
o is measured at different temperatures, preserving a certain viscosity. It is an intrinsic activation energy for overcoming a barrier separating the reagent and its isomer.
It is improbable to observe in the experiment α = 1.0, Equation (12), especially in a system with a high and narrow barrier, the solvent molecules will not have enough time to slow the reaction [
30].
One can suggest that in the following Equation (12), the experimental activation energy of isomerization E
obs should be equal to a sum:
where B (kJ/mol) is an activation energy of viscous flow obtained from Andrade Equation (14):
We will use experimental values obtained for the isomerization of the cyanine dye DODCI (3,3′-diethyloxadicarbocynanine iodide) in the ground (gr) and the excited (ex) singlet state in two solvents—ethanol (EtOH) and decan-1-ol (DecOH) [
36], see
Table 1.
There is good agreement between the observed calculated values of E
obs generally; see the two last columns in
Table 1. In the case of EtOH, “ex” demonstrates a worse agreement.
Equation (12) may be applied for the analysis of different cis–trans isomerizations.
Photochemical and thermal spiropyran–merocyanine interconversion in photoswitches depends upon viscosity. However, the process is too complex to be analyzed by Kramers’ theory [
42].
It was demonstrated that in some enzyme catalysis (allosteric regulation), free energy surfaces along the reaction coordinate and the diffusion coefficient of the reaction dynamics can be analyzed within the framework of Kramers’ theory [
43]. That approach allows the authors [
43] to explore various factors underlying the reaction kinetics and their allosteric regulation. The results provide new insights into the detailed mechanism of the allosteric fine-tuning of enzyme catalysis, and the tools developed can be valuable for studying other allosteric systems [
43].