# Proton-Coupled Electron Transfer in the Reaction of 3,4-Dihydroxyphenylpyruvic Acid with Reactive Species in Various Media

Authors: J. J. Fifen; Z. Dhaouadi; M. Nsangou; O. Holtomo; N. Jaidane
Journal: International Journal of Chemical Physics
 (2015)

---

## Abstract

The distinction of concerted proton-coupled electron transfer (CPCET) from sequential one as well as proton transfer-electron transfer (PT-ET) from electron transfer-proton transfer (ET-PT) in the O–H bond cleavage reactions in various media has always been a difficult task. In this work, the activation barrier of the CPCET mechanism, its rate constants, and reaction free energies related to ET-PT and PT-ET involving coreactive species were presented as good parameters to attempt the problem. DFT calculations were carried out studying the described pathways subsequent to the scavenging ofOH• and OBr- by the 3,4-DHPPA in various media. The solvation was described in a hybrid manner using IEF-PCM model conjointly with a model that takes into account some solute-solvent interactions. As a result, we found that the scavenging of hydroxyl radical by 3,4-DHPPA is thermodynamically governed by a one-step hydrogen atom transfer (CPCET) from the acid to the radical in all media. In kinetic viewpoint, CPCET still dominates in the vacuum and in nonpolar solvents, but in polar solvents it could compete strongly with the ET-PT mechanism so that the latter could slightly dominate.


---

## Body

## 1. Introduction

Understanding proton transfer (PT) or proton-coupled electron transfer (PCET) mechanism is very important in chemical processes. Such mechanisms can be subsequent to the O–H bond cleavage. The O–H bond cleavage subsequent to a transfer of the proton or a hydrogen atom from the concerned molecule to another molecule is governed by two main mechanisms: proton-coupled electron transfer (PCET) and proton transfer process (PT). In contrast to simple electron transfer (ET) or proton transfer (PT) reactions, PCET is more complex as both an electron and a proton must be transferred and their coupling strongly influences the process thermodynamically and kinetically. In PCET, the transfer of an electron and a proton may be sequential (stepwise) or concerted. In sequential transfer, either the electron or proton is transferred first. When the electron is transferred first, the process is termed electron transfer-proton transfer (ET-PT); otherwise, it is termed proton transfer-electron transfer (PT-ET) [1, 2]. In the concerted mechanism, the electron and proton are transferred simultaneously, and this is termed concerted proton-coupled electron transfer (CPCET) [1, 2].Usually, theoretically or experimentally, it is very difficult to distinguish CPCET from PT-ET or ET-PT. It is quite always an open problem. Recently, we studied the 3,4-dihydroxyphenylpyruvic acid (3,4-DHPPA) O–H bond cleavage in various media [3]. The scheme used in such a work does not consider the coreactive species (molecules or atoms which receive the transferred proton or hydrogen atom). It is only based on the calculated bond dissociation enthalpy (BDE), proton affinity (PA), and proton dissociation enthalpy (PDE). The main result derived in that work is that solvent notably influences and changes the nature of hydrogen atom cleavage, favoring by this way the proton transfer-electron transfer (PT-ET) mechanism in polar solvents and the concerted proton-coupled electron transfer (CPCET) mechanism in nonpolar solvents. Using the same scheme, the same results were pointed out by several other authors [4–9], studying the hydrogen atom cleavage on different systems in various media. The question is, do we have a complete insight into such mechanisms without taking into account coreactive species?Aiming to provide a reliable answer to the above question, a theoretical scheme involving reaction free energies (RFEs) related to the described processes was presented and the study is completed by a kinetic investigation. The application case concerns the O–H bond cleavage subsequent to the scavenging ofOH• and OBr- by the 3,4-dihydroxyphenylpyruvic acid (3,4-DHPPA) in the vacuum, cyclohexane, benzene, DMSO, and water. Solvents were chosen in regard to their increasing polarity or dielectric constant.3,4-DHPPA is a phenolic acid (PhA) reported as an alternative substrate for mammalian 4-hydroxyphenylpyruvate dioxygenase, which gives competitive inhibition versus 4-hydroxyphenylpyruvate [10, 11]. PhAs and their derivatives are widely present in plants (vegetables, fruits, grains, and spices), many being metabolites and several functions being attributed to them [12]. PhAs may contribute to the dark colour, bitter taste, and objectionable flavour of some fruits, leaves, and seeds. They have been considered as possible agents that influence toxicological, nutritional, sensory, and antioxidant properties of foods [12]. It has also been reported that phenolic acids may activate or inhibit microbial growth [12–16]. PhAs present antioxidant activity and exert this property through their capacity to donate hydrogen atoms, protons, and/or electrons to inhibit reactive oxygen species. The enol tautomer of 3,4-DHPPA, which was less studied as compared to the keto tautomer, was pointed out by Milane as being the fragment of Quercetin responsible for its higher antioxidant activity [17]. Thus, the results obtained involving 3,4-DHPPA could be generalized to Quercetin. Moreover, metal ions effects on the antioxidant activity of 3,4-DHPPA were studied in the vacuum [18]. A conformational study of this acid was also undertaken [3] and the likeliest conformer is depicted in Figure 1. The OH• radical has been chosen as the most unstable and the most reactive among the available reactive oxygen species (ROS), whereas OBr- is chosen as a nonradical oxygen species involved in several biological processes [19–21]. The thermodynamic study of the reaction of 3,4-DHPPA with the OH• radical was completed by a kinetic study involving tunnelling and recrossing effects.Figure 1
Equilibrium geometry of the likeliest conformer of 3,4-DHPPA in the vacuum, 6-31+G(d).


## 2. Computational Methods

### 2.1. Calculations of Thermodynamic Parameters

When coreactive species are not taken into account, thermodynamic parameters characterizing PCET mechanism can be evaluated through bond dissociation free energy (BDFE), ionization potential (IP), proton dissociation free energy (PDFE), proton affinity (PA), and electron free energy (ETFE) described as follows:(i)
BDFE of a given O–H bond is the free energy change of the dissociation reaction:(1)RO–H⟶RO•+H•(2)BDFE=ΔfG∘RO•+ΔfG∘H•-ΔfG∘RO–H.(ii)
IP of the antioxidant is the free energy change required to remove its electron:(3)RO–H⟶RO–H+•+e-(4)IP=ΔfG∘RO–H+•+ΔfG∘e--ΔfG∘RO–H.(iii)
PDFE is defined as the free energy change of the dissociation reaction:(5)RO–H+•⟶RO•+H+(6)PDFE=ΔfG∘RO•+ΔfG∘H+-ΔfG∘RO–H+•.(iv)
PA of the anionRO- is the free energy change of the dissociation reaction:
(7)RO–H⟶RO-+H+(8)PA=ΔfG∘RO-+ΔfG∘H+-ΔfG∘RO–H.(v)
ETFE represents the free energy change of the reaction:(9)RO-⟶RO•+e-(10)ETFE=ΔfG∘RO•+ΔfG∘e--ΔfG∘RO-.When coreactive species are taken into account in the cleavage process, the mechanism can be presented as summarized in Figure2 and the process is characterized by the reaction free energy (RFE). In such a scheme, RFE is defined as the sum of free energies of the donor and acceptor reactions. For example, RFE related to the mechanism ET(1) in the scheme is denoted by RFE-ET(1) and calculated as
(11)RFE-ET1=ΔfGsolv∘(ROH+)+ΔfGsolv∘(X-)M-ΔfGsolv∘ROH-ΔfGsolv∘X.Figure 2
Thermodynamic cycle showing ET-PT, PT-ET, and CPCET mechanisms involving coreactive species X. The indices 1 and 2 are devoted to ET-PT and PT-ET mechanisms, respectively. ROH is a phenolic acid.According to the above thermodynamic cycle, the possibility of ET-PT and PT-ET mechanisms can be evaluated, respectively, through RFE-ET(1) and RFE-PT(2), since ET(1) and PT(2) are, respectively, the first paths of these mechanisms.


### 2.2. Rate Constants Computations

Rate constants are computed as a function of temperatureT, in the zero-order semiclassical transition state theory [22]. Thus, the rate constant for the canonical transition state theory (TST) is given by
(12)k(T)=κ(T)×kTST(T),
where, for a bimolecular reaction as considered herein,
(13)kTSTT=σβh×Q≠TQATQBT×e-βΔV≠;kTST is the conventional TST rate constant at T temperature, without tunnelling correction. β=1/kBT; QA(T), QB(T), and Q≠(T) are, respectively, the partition functions of the reactants A, B and the transition state at temperature T. h and kB are, respectively, the Plank and Boltzmann constants. ΔV≠ is the classical barrier height while σ is the symmetry factor accounting for the possibility of two or more symmetry-related reaction paths. κ(T) is a ground state transmission coefficient at temperature T, which primarily accounts for tunnelling correction.Usually, tunnelling is considered only in the degree of freedom corresponding to the reaction coordinate. Within this one-dimensional formalism, we have used two different approximations of the transmission coefficient. The first approximation is a truncated parabola for a potential energy through which tunnelling happens. Therefrom, the transmission coefficient is [23]
(14)κT=β∑n=0∞‍-1n1-exp⁡β-n+1αΔV≠n+1α-β+1nα+β,
where β=1/kBT, α=1/ħν≠, and ν≠ is the imaginary frequency associated with the TS structure. For special cases,
(15)κT=1+124hν≠kT2,α≫βπβαsinπβα+ββ-α×eβ-αΔV≠,α≥βββ-α×eβ-αΔV≠-1,α≤β.
The case α≫β is known as the Wigner tunnelling approximation [24] and the associated transmission coefficient can be noted by κW. The two other cases (α≥β and α≤β) are, respectively, known as the first-case and the second-case of Skodje and Truhlar tunnelling approximation [23], and the associated transmission coefficients can be noted, respectively, by κST> and κST<.The second approximation of the transmission coefficient is based on a variational transition state theory (VTST). In this approximation, the transmission coefficient is a zero-order interpolated approximation [25, 26]. By zero-order, the authors mean that noab initio or DFT calculations at points other than reactants, saddle point, and products are available. In addition, tunnelling is assumed to occur along the minimum energy path (MEP) which is interpolated by an Eckart function [27]. Since one is not sure that the highest point in the MEP corresponds to the highest free energy, the variational TST rate constant is defined as
(16)kVTSTT=min⁡sσβh×Q≠T,sQATQBT×e-βΔV≠s,
where s is the distance from the generalized transition state along the MEP. This last approximation of the transmission coefficient is implemented in POLYRATE 2010-A [28] and termed zero-order interpolated variational transition state theory (IVTST-0) which is rigorously equivalent to the traditional transition state with zero-curvature tunnelling (TST/ZCT). The computations of all rate constants have been performed using our homemade Fortran program and POLYRATE 2010-A.


### 2.3. DFT, Basis Sets, and Solvation Methods Description

Due to the solute cavity description, implicit solvation models cannot provide a good description for short range solute-solvent interactions. In this work, such interactions were considered explicitly and the solvation was treated in a hybrid manner considering some explicit interactions between the particle to be transferred (during the scavenging mechanism) and a solvent molecule. The rest of solvent was treated implicitly using integral equation formalism polarized continuum model (IEF-PCM) [29]. Thus, aiming to compute the reaction free energies in solvent, we strongly needed solvation free energies of H•, H+, and e-. All free energy changes were computed at a temperature of 298.15 K and a pressure of 1 atm. Geometry optimizations were carried out using the density functional theory (DFT) methods implemented in the GAUSSIAN03 computational package [30].DFT was chosen because of the excellent compromise between the computational time and the description of the electronic correlation. Throughout our calculations, the B3LYP hybrid functional which consists of Becke’s three parameters exact exchange functional (B3) [31] combined with the nonlocal gradient corrected correlation functional of Lee-Yang-Parr (LYP) [32] was used. For systems with old number of electrons (radicals involved in this work), we used the restricted open shells approach (ROB3LYP). The basis sets used in this work were 6-31+G* and 6-311++G** [33–35]. The last one is used especially for solvation free energies of the proton and electron. Full geometry optimization was carried out without symmetry constraints up to convergence (largest component of nuclear gradient equal to 10-6 a.u./bohr and the change in total energy less than 10-7 a.u.). A subsequent vibrational frequency calculation was undertaken in order to confirm that the resulting optimized geometries correspond to minima (no negative frequencies) or first-order saddle points (one negative frequency) on the PES, as we were looking for equilibria or transition states, respectively. The natural bond orbital (NBO) technique [36, 37] was used for natural population analysis (NPA). All these calculations were undertaken in the vacuum and in polar and nonpolar solvents. Farther in the work, G3B3 method was used for counterpoise computations. G3B3 is the Gaussian-3 (G3) method using B3LYP (B3) structures and frequencies for very accurate energies computations [38, 39].


## 2.1. Calculations of Thermodynamic Parameters

When coreactive species are not taken into account, thermodynamic parameters characterizing PCET mechanism can be evaluated through bond dissociation free energy (BDFE), ionization potential (IP), proton dissociation free energy (PDFE), proton affinity (PA), and electron free energy (ETFE) described as follows:(i)
BDFE of a given O–H bond is the free energy change of the dissociation reaction:(1)RO–H⟶RO•+H•(2)BDFE=ΔfG∘RO•+ΔfG∘H•-ΔfG∘RO–H.(ii)
IP of the antioxidant is the free energy change required to remove its electron:(3)RO–H⟶RO–H+•+e-(4)IP=ΔfG∘RO–H+•+ΔfG∘e--ΔfG∘RO–H.(iii)
PDFE is defined as the free energy change of the dissociation reaction:(5)RO–H+•⟶RO•+H+(6)PDFE=ΔfG∘RO•+ΔfG∘H+-ΔfG∘RO–H+•.(iv)
PA of the anionRO- is the free energy change of the dissociation reaction:
(7)RO–H⟶RO-+H+(8)PA=ΔfG∘RO-+ΔfG∘H+-ΔfG∘RO–H.(v)
ETFE represents the free energy change of the reaction:(9)RO-⟶RO•+e-(10)ETFE=ΔfG∘RO•+ΔfG∘e--ΔfG∘RO-.When coreactive species are taken into account in the cleavage process, the mechanism can be presented as summarized in Figure2 and the process is characterized by the reaction free energy (RFE). In such a scheme, RFE is defined as the sum of free energies of the donor and acceptor reactions. For example, RFE related to the mechanism ET(1) in the scheme is denoted by RFE-ET(1) and calculated as
(11)RFE-ET1=ΔfGsolv∘(ROH+)+ΔfGsolv∘(X-)M-ΔfGsolv∘ROH-ΔfGsolv∘X.Figure 2
Thermodynamic cycle showing ET-PT, PT-ET, and CPCET mechanisms involving coreactive species X. The indices 1 and 2 are devoted to ET-PT and PT-ET mechanisms, respectively. ROH is a phenolic acid.According to the above thermodynamic cycle, the possibility of ET-PT and PT-ET mechanisms can be evaluated, respectively, through RFE-ET(1) and RFE-PT(2), since ET(1) and PT(2) are, respectively, the first paths of these mechanisms.


## 2.2. Rate Constants Computations

Rate constants are computed as a function of temperatureT, in the zero-order semiclassical transition state theory [22]. Thus, the rate constant for the canonical transition state theory (TST) is given by
(12)k(T)=κ(T)×kTST(T),
where, for a bimolecular reaction as considered herein,
(13)kTSTT=σβh×Q≠TQATQBT×e-βΔV≠;kTST is the conventional TST rate constant at T temperature, without tunnelling correction. β=1/kBT; QA(T), QB(T), and Q≠(T) are, respectively, the partition functions of the reactants A, B and the transition state at temperature T. h and kB are, respectively, the Plank and Boltzmann constants. ΔV≠ is the classical barrier height while σ is the symmetry factor accounting for the possibility of two or more symmetry-related reaction paths. κ(T) is a ground state transmission coefficient at temperature T, which primarily accounts for tunnelling correction.Usually, tunnelling is considered only in the degree of freedom corresponding to the reaction coordinate. Within this one-dimensional formalism, we have used two different approximations of the transmission coefficient. The first approximation is a truncated parabola for a potential energy through which tunnelling happens. Therefrom, the transmission coefficient is [23]
(14)κT=β∑n=0∞‍-1n1-exp⁡β-n+1αΔV≠n+1α-β+1nα+β,
where β=1/kBT, α=1/ħν≠, and ν≠ is the imaginary frequency associated with the TS structure. For special cases,
(15)κT=1+124hν≠kT2,α≫βπβαsinπβα+ββ-α×eβ-αΔV≠,α≥βββ-α×eβ-αΔV≠-1,α≤β.
The case α≫β is known as the Wigner tunnelling approximation [24] and the associated transmission coefficient can be noted by κW. The two other cases (α≥β and α≤β) are, respectively, known as the first-case and the second-case of Skodje and Truhlar tunnelling approximation [23], and the associated transmission coefficients can be noted, respectively, by κST> and κST<.The second approximation of the transmission coefficient is based on a variational transition state theory (VTST). In this approximation, the transmission coefficient is a zero-order interpolated approximation [25, 26]. By zero-order, the authors mean that noab initio or DFT calculations at points other than reactants, saddle point, and products are available. In addition, tunnelling is assumed to occur along the minimum energy path (MEP) which is interpolated by an Eckart function [27]. Since one is not sure that the highest point in the MEP corresponds to the highest free energy, the variational TST rate constant is defined as
(16)kVTSTT=min⁡sσβh×Q≠T,sQATQBT×e-βΔV≠s,
where s is the distance from the generalized transition state along the MEP. This last approximation of the transmission coefficient is implemented in POLYRATE 2010-A [28] and termed zero-order interpolated variational transition state theory (IVTST-0) which is rigorously equivalent to the traditional transition state with zero-curvature tunnelling (TST/ZCT). The computations of all rate constants have been performed using our homemade Fortran program and POLYRATE 2010-A.


## 2.3. DFT, Basis Sets, and Solvation Methods Description

Due to the solute cavity description, implicit solvation models cannot provide a good description for short range solute-solvent interactions. In this work, such interactions were considered explicitly and the solvation was treated in a hybrid manner considering some explicit interactions between the particle to be transferred (during the scavenging mechanism) and a solvent molecule. The rest of solvent was treated implicitly using integral equation formalism polarized continuum model (IEF-PCM) [29]. Thus, aiming to compute the reaction free energies in solvent, we strongly needed solvation free energies of H•, H+, and e-. All free energy changes were computed at a temperature of 298.15 K and a pressure of 1 atm. Geometry optimizations were carried out using the density functional theory (DFT) methods implemented in the GAUSSIAN03 computational package [30].DFT was chosen because of the excellent compromise between the computational time and the description of the electronic correlation. Throughout our calculations, the B3LYP hybrid functional which consists of Becke’s three parameters exact exchange functional (B3) [31] combined with the nonlocal gradient corrected correlation functional of Lee-Yang-Parr (LYP) [32] was used. For systems with old number of electrons (radicals involved in this work), we used the restricted open shells approach (ROB3LYP). The basis sets used in this work were 6-31+G* and 6-311++G** [33–35]. The last one is used especially for solvation free energies of the proton and electron. Full geometry optimization was carried out without symmetry constraints up to convergence (largest component of nuclear gradient equal to 10-6 a.u./bohr and the change in total energy less than 10-7 a.u.). A subsequent vibrational frequency calculation was undertaken in order to confirm that the resulting optimized geometries correspond to minima (no negative frequencies) or first-order saddle points (one negative frequency) on the PES, as we were looking for equilibria or transition states, respectively. The natural bond orbital (NBO) technique [36, 37] was used for natural population analysis (NPA). All these calculations were undertaken in the vacuum and in polar and nonpolar solvents. Farther in the work, G3B3 method was used for counterpoise computations. G3B3 is the Gaussian-3 (G3) method using B3LYP (B3) structures and frequencies for very accurate energies computations [38, 39].


## 3. Results and Discussion

Throughout this work, the O–H bond cleavage through the reactions of 3,4-DHPPA withOH• or OBr- was considered. In such reactions, the reactive species (RS) considered react with the acid through the phenolic hydrogen atoms (or protons) of the acid, that is, H3 (or H3+), H4 (or H4+), and H8 (or H8+), since phenolic hydrogen atoms (or protons) are more acidic than carboxylic ones [12, 40–42]. Thus, in regard to the diagram of Figure 2, X may be replaced by OH• or OBr- when the considered coreactive is OH• or OBr-, respectively. Figures 3 and 4 depict the optimized geometries of supermolecules resulting from the reactions of 3,4-DHPPA with OH• and OBr- at the 3-OH attacked site in the vacuum, respectively. Aiming to localize the transition states of these reactions, 3D-relaxed potential energy surfaces (3D-RPES) were performed in the vacuum at the B3LYP/6-31+G* level of theory on the supermolecule obtained, varying the distance O3–H3 and the angle C3–O3–H3, Figures 5 and 6.Figure 3
Equilibrium geometry of the supermolecule formed after the reaction of 3,4-DHPPA withOH•, at B3lyp/6-31+G(d).Figure 4
Equilibrium geometry of the supermolecule formed after the reaction of 3,4-DHPPA withOBr-, at B3lyp/6-31+G(d).Figure 5
3D-PES of the reaction of 3,4-DHPPA withOH•, varying the distance O3–H3 and the angle C3–O3–H3, at B3lyp/6-31+G(d).Figure 6
3D-PES of the reaction of 3,4-DHPPA withOBr-, varying the distance O3–H3 and the angle C3–O3–H3, at B3lyp/6-31+G(d).Therefrom, it is noteworthy that the scavenging of hydroxyl radical by 3,4-DHPPA passes through a transition state (TS) while that ofOBr- anion is barrierless. Moreover, it is well-known that the reaction of a phenolic acid (PhA) with neutral ROS leads to PCET mechanism while that with negatively charged ROS undergoes a barrierless PT mechanism [3, 12, 18, 43, 44]. This result may be generalized to PhAs which only have carbon, oxygen, and hydrogen in their skeleton (e.g., Quercetin may be a good candidate) and are involving the same RS or other. Therefore, the scavenging of the hydroxyl radical and the OBr- anion by a PhA (ROH) can be summarized by the following, respectively:
(17)ROH+OH•⟶RO⋯H(TS)⋯OH⟶RO•+H2O(18)ROH+OBr-⟶RO-+HOBrNatural populations analysis (NPA) and spin densities on reactants, TS, and products, conjointly with singly occupied molecular orbitals (SOMOs) distributions, may provide a good insight in the nature of the particle transferred during the scavenging process.


### 3.1. Solvation Free Energies of the Electron, Proton, and Hydrogen Atom

Solvation free energies of the electron, proton, and hydrogen atom are, respectively, the free energies change of reactions below [3]:
(19)solvent(solv)+e(g)-⟶solvent(solv)-(20)solvent(solv)+H(g)+⟶solvent-H(solv)+(21)solvent(solv)+H(g)•⟶solvent-Hsolv•
where solvent(solv) represents a molecule of solvent in its cavity, X(solv) means that X is solvated, and solvent may be one of those used in this work. So, “solv” stands for cyclohexane, benzene, and so forth. Free energies of X(solv) are obtained from geometry optimizations followed by frequency calculations. As a matter of fact, (20) with water as solvent could be written as
(22)H2O(water)+H(g)+⟶H2O-H(water)+
It is worthwhile that applying the proposed scheme in the cases of chloroform and dichloromethane is very difficult, since attaching a single electron or proton to the host molecule (chloroform or dichloromethane molecule) leads to the fragmentation of the formed complex. Thus, these solvents are special cases and could be investigated through a more robust scheme. It is also worth mentioning that the solvation free energies of a single electron, proton, and hydrogen atom are not necessary to derive RFEs but they would be very useful to explain the possibility of partial reactions involved in the studied mechanisms. Thus, such parameters help in the thorough investigation of proton-coupled electron transfer in O–H bond cleavage.Gas phase free energy of the electron used here is that obtained by Bartmess:ΔfG∘e-=-3.72 kJ/mol [45]. Our calculated gas phase free energy of the proton is ΔfG∘H+=-26.26 kJ/mol. Similar results were also proposed recently [46, 47]. Solvation free energies of H atom in studied solvents were taken from published solvation free energies and entropies listed in [48–50] and reported in Table 1. In organic solvents, ΔfHsolv(H•)≃ΔfHsolv(H2)=(5±1) kJ/mol [48–50]. Thus, for organic solvents for which data were unavailable, the value of 5 kJ/mol is assumed. Likewise, ΔfGsolv(H•)≃ΔfGsolv(H2)≃-20kJ/mol. Therefore, this value was retained for the solvents with unavailable experimental values. For water, H-atom hydration free energy reached 27.7 kJ/mol. Besides, since the only available experimental values for solvation free energies of proton and electron were their hydration free energies (ΔfGhydr∘H+=-1104.5±0.3kJ/mol [51], ΔfGhydr∘e-=-156.8kJ/mol [52]), solvation free energies of these particles in different solvents used here were computed at IEF-PCM/B3LYP/6-311++G** level of theory. In comparison to experimental results, the relative errors on calculated free energies changes for hydrated proton and electron are 9.8% and 36.8%, respectively. The first error may be acceptable, while the last one is too high. This could mainly be due to the fact that, in our explicit solvation treatment, the first solvation shell is not completely considered. The more accurate scheme should consider at least the first solvation shell completely. This could be undertaken through a thorough investigation of the clustering of protonated, hydrogenated, or electronized solvent. Such a task is beyond the scope of this work and it is devoted to the forthcoming works. For the moment, we assumed that these relative errors are constant in all media. Thus, all calculated ΔfGsolv∘(H•) and ΔfGsolv∘(e-) were shifted by an amount of the percentage error and presented in Table 1.Table 1
Solvation free energies of hydrogen atom (ΔfGsolv∘(H•)), proton (ΔfGsolv∘(H+)), and electron (ΔfGsolv∘(e-)) in kJ/mol, at 6-311++G(d,p) in various media. ε and μ are, respectively, the dielectric constant and the dipole moment of the considered media.
Solvent
ε
μ
Δ
f
G
solv
∘
(
H
•
)
Δ
f
G
solv
∘
(
H
+
)
Δ
f
G
solv
∘
(
e
-
)
Cyclohexane
2.02
0.00
−20
−879.6
−56.4
Benzene
2.28
0.00
−20
−1076.0
−69.6
Methanol
32.7
1.71
−20
−1119.7
−140.6
DMSO
46.7
3.90
−20
−1199.3
−141.2
Water
78.4
1.80
27.7
−1104.5
−156.8


### 3.2. PCET and PT Mechanisms via SOMOs, NPAs Charges, and Spin Densities Analysis

As a proton and an electron transfer from a donor to an acceptor site in a molecular system, the charge, spin density, and molecular orbitals (MOs) localized on the donor and acceptors sites change simultaneously. Thus, the calculation of the charge, spin density, and plots of MOs along a reaction coordinate (transferring hydrogen or proton) provides useful insight in understanding the PCET mechanism.SOMOs distributions on the equilibrium structures resulting from the reaction of 3,4-DHPPA withOH• in all media show that the radical is entirely delocalised on the phenoxyl moiety (see Figure 7) and could lead to the conclusion that the scavenging of OH• by 3,4-DHPPA is governed by a hydrogen atom transfer (PCET mechanism) whatever the media.Figure 7
SOMO distribution on the equilibrium structure of the supermolecule resulting from the reaction of 3,4-DHPPA withOH• radical: a B3lyp/6-31+G(d). SOMO is entirely delocalized on the phenoxyl moiety, suggesting hydrogen atom transform from the acid to the radical.NPAs and spin densities computations were performed at the same level of computations as that of geometry optimization. Computed NPA charges for products of the reaction illustrated by (17) are neutral, while spin densities are, respectively, 1.00 and 0.00 for the phenoxyl moiety and the water molecule formed, whatever the medium. In addition, the computation of such charges on reactants, products, and transition state of this reaction shows that the transferred particle has a positive charge with no spin associated, Table 2. Moreover, it is worth mentioning that the spin density associated with the phenoxyl moiety at the TS is increasing while that of the hydroxyl radical is decreasing, going from reactants to products. This indicates that the radical transfer is occurring. Moreover, at the TS, the negative NPA charges on OH increase with the solvent polarity while its spin density decreases with it, until spin density is no longer associated with OH in polar solvents. Note that, in such solvents, the electron spin density on OH is quite nil, while it is a unit on the phenoxyl moiety (RO), Table 2. Moreover, the electron charge on OH is more than the half of elementary electron charge. Thus, while the reaction of (17) is proceeding in polar solvents, the radical character is completely transferred from the free radical (OH•) to the phenoxyl moiety before the proton is transferred from the acid to the free radical. Therefore, in polar solvents, the scavenging of the hydroxyl radical by DHPPA is stepwise and the electron is transferred first. In the vacuum and nonpolar solvents, the electron spin densities on OH and RO are about 0.2 and 0.8, respectively. In addition, the electron charge on OH is less than the unit. Thus, the located TS for the reaction of (17) is a mixture of ET and PT characters. Therefore, a clear distinction of ET-PT, PT-ET, and CPCET mechanisms for such a reaction is not ensured using NPA charges and spin densities at the TS. Parameters based on reaction free energies and activation energies could be good indicators for such a target.Table 2
NPA charges (e) and spin densities carried by elements of the reaction at the transition state. The superscripts 3 and 4 on DHPPA indicate that the hydrogen atom concerned by the abstraction is H3 and H4, respectively. The symbol X denotes the unknown transferred particle.
Media
NPA charges
Spin densities
RO
X
OH
RO
X
OH
DHPPA3
Vacuum
−0.02
0.52
−0.50
0.73
0.00
0.27
Cyclohexane
0.11
0.52
−0.63
0.82
0.00
0.18
Benzene
0.10
0.53
−0.63
0.83
0.00
0.17
Methanol
0.26
0.52
−0.78
0.96
0.00
0.04
Water
0.25
0.52
−0.79
0.98
0.00
0.02
DHPPA4
Vacuum
−0.02
0.52
−0.50
0.76
0.00
0.24
Cyclohexane
—
—
—
—
—
—
Benzene
0.12
0.53
−0.65
0.89
0.00
0.11
Methanol
0.28
0.55
−0.83
0.99
0.00
0.01
Water
0.29
0.55
−0.84
0.99
0.00
0.01Moreover, we plotted the SOMOs at the TS of this reaction (17) in different solvents. These SOMOs are delocalised entirely on the supermolecule, Figures 8 and 9. These SOMOs show that the electron and proton transfer from different orbitals on the donor to different orbitals on the acceptor. The same spin density distribution was noted in other solvents (see Figure S1 in the Supplementary Material available online at http://dx.doi.org/10.1155/2015/835707).Figure 8
SOMO distribution on a transition state structure of the reaction of D-DHPPA withOH• radical in the vacuum.Figure 9
SOMO distribution on a transition state structure of the reaction of D-DHPPA withOH• radical in water.Thus, based on the works of Mayer et al. [53] and Hammes-Schiffer and Iordanova [54], these SOMOs distributions allow us to conclude that the scavenging of hydroxyl radical by D-DHPPA is governed by a CPCET mechanism, whatever the medium. Hence, in contrast to the results published before (the scavenging of hydroxyl radical by DHPPA is governed by CPCET in nonpolar solvents and PT-ET mechanisms in polar solvents [3]), PT-ET mechanism may no longer be favored in polar solvents if coreactive species are taken into account.Regarding the reaction illustrated by (18), for which there is neither radical nor transition state involved, NPA calculations show that the resulting phenoxyl moiety bears a negative charge while the HOBr formed is neutral. This reveals a proton transfer and confirms the well-known result concerning the nature of the scavenging of a negatively charged species by a PhA.


### 3.3. Thermodynamic Investigation of PCET and PT Mechanisms

Reaction free energies relevant to different pathways related to PT-ET, ET-PT, CPCET, and PT mechanisms are based on coreaction free energies (Table3) and abstraction free energies of H•, H+, and e- (Table 4). Coreaction free energies are assumed to be the free energies of the reactions involving coreactive species and abstracted or transferred particle (H•, H+, and e-). It is noteworthy that RFE-CPCET may not be used in concert with RFE-ET(1) and RFE-PT(2) to derive the preferred mechanism (PCET? PT-ET? or ET-PT?) in a given medium. In fact, the former characterizes the one-step H-atom transfer, while the other characterizes one step in a stepwise H-atom transfer. Thus, a reliable thermodynamic way to derive the preferred mechanism in a given medium is to compare RFE-ET(1), RFE-PT(2), and the activation barrier (ΔG≠) of the CPCET mechanism.Table 3
Free energies of reactions involving coreactive species and abstracted or transferred particle (H•, H+, and e-) in kJ/mol at B3LYP/6-31+G(d).
Media
•OH + e-→ OH−
•OH + •H → H2O
•OH + H+→H2O+
OBr− + H+→ HOBr
Vacuum
−160.3
−436.0
−544.1
−1438.1
Cyclohexane
−329.0
−418.8
35.4
−536.3
Benzene
−332.2
−419.1
212.6
−327.4
Methanol
−397.7
−422.0
90.1
−180.3
DMSO
−400.1
−422.1
165.8
−98.6
Water
−387.3
−469.8
67.3
−191.3Table 4
Abstraction free energies of electron, proton, and hydrogen atom from a phenolic acid ROH in kJ/mol.
Media
ROH→ ROH+ + e-
ROH→ RO− + H+
ROH→ RO• + •H
DHPPA3
Vacuum
699.4
1339.2
270.0
Cyclohexane
604.8
488.3
254.4
Benzene
581.9
283.5
254.8
Methanol
428.7
172.9
263.6
DMSO
428.0
91.9
263.7
Water
408.1
185.8
310.8
DHPPA4
Vacuum
699.4
1355.3
288.4
Cyclohexane
604.8
500.3
264.0
Benzene
581.9
295.4
263.2
Methanol
428.7
174.4
256.3
DMSO
428.0
92.9
256.1
Water
408.1
186.5
304.6
DHPPA8
Vacuum
699.40
1363.5
280.5
Cyclohexane
604.78
506.9
257.3
Benzene
581.87
301.6
256.8
Methanol
428.71
174.9
252.7
DMSO
428.04
94.9
253.0
Water
408.1
187.4
298.5


#### 3.3.1. CPCET Mechanism

Previously, studies on the CPCET mechanism at 3-OH, 4-OH, and 8-OH positions in 3,4-DHPPA have shown that BDE3 < BDE8 < BDE4 in the vacuum and nonpolar solvents [3]. This order changes in polar solvents and becomes BDE8 < BDE4 < BDE3. These behaviors were also obtained when coreactive species were taken into account, Figure 10 (or Table S1). The justification of such a result was based on the relative strengths of hydrogen bonds O4H4⋯O8, O8H8⋯O9 and C2H2⋯O8 in 3,4-DHPPA. Details were discussed in our previous manuscript [3]. Thus, it may be sufficient to study the CPCET mechanism using BDFEs or BDEs rather than RFEs, in order to derive the relative ability to abstract the H-atom through the abstraction sites of the PhA. Note that the negative values obtained for RFE-CPCET in all media are the confirmation that the transfer of the H-atom from the acid to the hydroxyl radical is always possible in such media.Figure 10
Reaction free energies for the reactionROH+OH•→RO•+H2O, in different media. ROH stands for 3,4-DHPPA, while RO• stands for its phenoxyl moiety.


#### 3.3.2. PT-ET Mechanisms

Assume that the reaction of 3,4-DHPPA withOH• radical processes through a PT-ET mechanism leads to the following reaction:
(23)ROH+OH•⟶RO-+H2O+⟶RO•+H2O
All the RFE-PT(2)(OH) values obtained are positive, Table S1. Thus, such a reaction is always thermodynamically unfavorable, especially in the vacuum. This may be explained by the fact that the fragmentation of a neutral molecule into a charged species is more difficult in the vacuum than in solvent, since solvent will make the molecule more polar and thereafter easy to be fragmented in charged species. Moreover, the reaction involves the reaction of OH• with a proton H+. This reaction is usually less favorable, especially in solvent. Indeed, the free radical OH• is neutral and polar, so it cannot interact easily with a charged and nonpolar molecular system. This physical insight is confirmed through RFEs of the reactions OH•+H+→H2O+, Table 3.Furthermore, the scavenging of theOH• radical by 3,4-DHPPA through the PT-ET mechanism is done in all abstraction sites (3-OH, 4-OH, and 8-OH) with the same ability in polar solvents, Figure 11. In the vacuum and nonpolar solvents, the 3-OH site is slightly preferred to the 4-OH one, which in turn is slightly preferred to the 8-OH site. In addition, this mechanism (PT-ET) is more favorable in solvent than in the vacuum, and the more the solvent polarity, the more favorable the mechanism, Figure 12. This may be explained through its dipole moment by the fact that solvent eases charge separation and, the more the solvent polarity, the easier the charge separation. These results are consistent with those found without taking into account the coreactive species [3].Figure 11
Reaction free energies evaluating the possibility of the mechanism PT-ET in the reaction of 3,4-DHPPA withOH• in different media. These RFEs are calculated as the sum of partial reaction free energies of the following reactions: ROH→RO-+H+, OH+H+→H2O+. ROH stands for 3,4-DHPPA, while RO- stands for its deprotonated structure. 3-OH, 4-OH, and 8-OH denote the cleavage sites related, respectively, to the protons H3+, H4+, and H8+ of 3,4-DHPPA.Figure 12
Reaction free energies and activation barriers evaluating, respectively, the possibility of the mechanisms PT-ET, ET-PT, and CPCET, in the reaction of 3,4-DHPPA withOH• in different media. For CPCET, the reaction involved is ROH+OH•→RO•+H2O. For PT-ET, the partial reactions involved are ROH→RO-+H+, OH+H+→H2O+. For ET-PT, the partial reactions involved are ROH→ROH++e-, OH•+e-→OH-. ROH stands for 3,4-DHPPA.By comparing the activation barrierΔG≠ of the CPCET mechanism with RFE-PT(2), one can predict that, even kinetically, the CPCET mechanism would be preferred to the PT-ET mechanism in all media, Figure 12. The justification of this statement could also be performed by a kinetic investigation of the PT-ET mechanism which is beyond the scope of this work.


#### 3.3.3. ET-PT Mechanism

Assume that the reaction of 3,4-DHPPA withOH• radical processes through a ET-PT mechanism leads to the following reaction:
(24)ROH+OH•⟶ROH++OH-⟶RO•+H2O
This reaction is thermodynamically more favorable in solvent than in the vacuum, and the more the solvent polarity, the more favorable the mechanism, Figure 12. These results are consistent with the earlier explanation relevant to the easiness of the solvent to separate charged species. They are also consistent with those found without taking into account the coreactive species [3].By comparing the activation barrierΔG≠ of the CPCET mechanism with RFE-ET(1), one can predict that, even kinetically, the CPCET mechanism would be preferred in the vacuum and nonpolar solvents compared to ET-PT mechanism which in turn would be preferred in polar solvents, Figure 12. Furthermore, in polar solvents, ET-PT mechanism competes strongly with CPCET mechanism, Figure 12. The justification of these kinetic viewpoints could also be performed by a kinetic investigation of the ET-PT mechanism which is beyond the scope of this work.Therefore, during the scavenging of the hydroxyl radical by DHPPA, the transfer of the hydrogen atom is thermodynamically concerted in all media and favoring by this way the CPCET mechanism in all media. However, in kinetic viewpoint, this mechanism could compete strongly with the ET-PT mechanism so that the latter could slightly dominate. These results which corroborate the results obtained through the SOMOs distribution on the TS structures are also in agreement with those derived in polar solvents using NPA charges and electron spin densities at TS. They also bring a clear decision on the preferred mechanism in the vacuum and nonpolar solvents. It is noteworthy that the obtained results are consistent in the vacuum and nonpolar solvents with previous results [3, 5, 6, 9] which do not involve coreactive species. The differences were noted in polar solvents. This difference could be explained by the fact that the approach free of coreactives involvement is unable to give a kinetic insight of such a reaction.


#### 3.3.4. PT Mechanism

PT mechanism is involved in the scavenging of the anionOBr- by 3,4-DHPPA. This mechanism is less favorable in solvents than in the vacuum, and the more the solvent polarity increases, the less the mechanism is favored, Figure 13 (or Table S1). This result may also be explained by the easiness to separate charged species in the vacuum compared to in solvents, especially in polar solvents.Figure 13
Reaction free energies related to proton transfer (PT) involving the scavenging of the anionOBr- by 3,4-DHPPA in different media. RFE-PT is calculated as the sum of partial RFEs of the reactions ROH→RO-+H+ and OBr-+H+→HOBr.


### 3.4. Kinetic Study of the PCET Mechanism

To complete the above thermodynamic description, kinetic investigation of the PCET mechanism (17) was undertaken in various media. The media retained were vacuum, cyclohexane, benzene, methanol, and water. The basis set superposition error (BSSE) was computed for the transition state structure at the G3B3 level of computation in the vacuum. For the TS structures related to the H3-atom, H4-atom, and H8-atom abstractions in the vacuum, the obtained values were 12.6, 12.4, and 12.3kJ/mol, respectively. The same values were used in solvents, since up to now there are no possible counterpoise computations in solvent. Rate constants were computed at T=298.15 K, using TST/ST and TST/ZCT methods. The TST/ST method is the traditional TST method with κST< for the transmission coefficient and is implemented in our homemade Fortran program. The computed activation free energies (ΔG≠) were achieved including zero-point energies.The rate constants obtained from both described methods are quite consistent (see Figure14) and lead to the same information that rate constants decrease as the solvent polarity increases. The result is also consistent with the increasing distance d(On⋯Hn) with the solvent polarity, n=3,4. It is also in good agreement with the decreasing frequency ν(On⋯Hn). Indeed, the longer distance matches a lower frequency, since a longer distance corresponds to a weaker bond. The discrepancy observed in the vacuum between both methods of rate constants calculation could be due to the large Wigner tunnelling effects in such a medium. This shows that the tunnelling effects on hydrogen atom transfer could be more important in the vacuum than in a solvent. Indeed, the TS vibrational frequency in the vacuum is -2429cm-1, while in a solvent the value is about -1500cm-1, which is lower than the former.Figure 14
Rate constants for the scavenging of the free radicalOH• by the 3,4-DHPPA, computed using TST/ST and TST/ZCT methods in different media. The attacked site considered is O3–H3.This result is consistent with the increasing activation barrier with solvent polarity (Table5), suggesting that the higher the solvent polarity, the slower the CPCET mechanism. This result is not surprising, since it is shown early that solvent makes the one-step H-atom transfer harder.Table 5
Parameters characterizing the O–H bond cleavage in selected media.d(On⋯Hn) in Å is the distance from the oxygen On of the phenoxyl moiety to the moved hydrogen atom Hn. ν(On⋯Hn) in cm-1 is the corresponding frequency associated with the movement of this hydrogen atom, while kTST/ST and kTST/ZCT in cm3⋅molecule-1⋅s-1 are the associated rate constants using, respectively, TST/ST and TST/ZCT methods. ΔG≠ in kJ/mol is the activation free energy, and n = 3, 4. NA stands for not applicable.
Media
d
ν
Δ
G
≠
k
TST/ST
k
TST/ZCT
Related to theO3–H3 attacked site
Vacuum
0.969
3769
31.2
6.8× 10-14
5.0× 10-13
Cyclohexane
0.975
3645
47.0
1.6× 10-14
3.5× 10-14
Benzene
0.976
3625
53.8
5.1× 10-15
7.1× 10-15
Methanol
0.990
3331
71.8
3.6× 10-15
2.5× 10-16
Water
0.991
3326
78.3
4.8× 10-16
1.6× 10-16
Related to theO4–H4 attacked site
Vacuum
0.964
3710
68.1
1.9× 10-17
6.2× 10-16
Cyclohexane
0.975
3661
—
—
—
Benzene
0.976
3652
78.6
2.9× 10-18
NA
Methanol
0.983
3493
93.3
1.7× 10-19
6.1× 10-18
Water
0.985
3457
90.9
—
—Moreover, whatever the medium, the hydrogen atomH3 is transferred more rapidly than the hydrogen atom H4, since the related rate constant of H3 transfer is higher than that related to H4 transfer in a given medium, Table 5. In addition, the activation free energy and rate constants using TST/ST and TST/ZCT methods for the O8–H8 bond cleavage in the vacuum are 70.1 kJ/mol, 1.3×10-17 cm3·molecule−1·s−1, and 4.2×10-16 cm3·molecule−1·s−1, respectively. Thus, the rate constant related to H4 transfer is higher than that related to H8 in the vacuum.


## 3.1. Solvation Free Energies of the Electron, Proton, and Hydrogen Atom

Solvation free energies of the electron, proton, and hydrogen atom are, respectively, the free energies change of reactions below [3]:
(19)solvent(solv)+e(g)-⟶solvent(solv)-(20)solvent(solv)+H(g)+⟶solvent-H(solv)+(21)solvent(solv)+H(g)•⟶solvent-Hsolv•
where solvent(solv) represents a molecule of solvent in its cavity, X(solv) means that X is solvated, and solvent may be one of those used in this work. So, “solv” stands for cyclohexane, benzene, and so forth. Free energies of X(solv) are obtained from geometry optimizations followed by frequency calculations. As a matter of fact, (20) with water as solvent could be written as
(22)H2O(water)+H(g)+⟶H2O-H(water)+
It is worthwhile that applying the proposed scheme in the cases of chloroform and dichloromethane is very difficult, since attaching a single electron or proton to the host molecule (chloroform or dichloromethane molecule) leads to the fragmentation of the formed complex. Thus, these solvents are special cases and could be investigated through a more robust scheme. It is also worth mentioning that the solvation free energies of a single electron, proton, and hydrogen atom are not necessary to derive RFEs but they would be very useful to explain the possibility of partial reactions involved in the studied mechanisms. Thus, such parameters help in the thorough investigation of proton-coupled electron transfer in O–H bond cleavage.Gas phase free energy of the electron used here is that obtained by Bartmess:ΔfG∘e-=-3.72 kJ/mol [45]. Our calculated gas phase free energy of the proton is ΔfG∘H+=-26.26 kJ/mol. Similar results were also proposed recently [46, 47]. Solvation free energies of H atom in studied solvents were taken from published solvation free energies and entropies listed in [48–50] and reported in Table 1. In organic solvents, ΔfHsolv(H•)≃ΔfHsolv(H2)=(5±1) kJ/mol [48–50]. Thus, for organic solvents for which data were unavailable, the value of 5 kJ/mol is assumed. Likewise, ΔfGsolv(H•)≃ΔfGsolv(H2)≃-20kJ/mol. Therefore, this value was retained for the solvents with unavailable experimental values. For water, H-atom hydration free energy reached 27.7 kJ/mol. Besides, since the only available experimental values for solvation free energies of proton and electron were their hydration free energies (ΔfGhydr∘H+=-1104.5±0.3kJ/mol [51], ΔfGhydr∘e-=-156.8kJ/mol [52]), solvation free energies of these particles in different solvents used here were computed at IEF-PCM/B3LYP/6-311++G** level of theory. In comparison to experimental results, the relative errors on calculated free energies changes for hydrated proton and electron are 9.8% and 36.8%, respectively. The first error may be acceptable, while the last one is too high. This could mainly be due to the fact that, in our explicit solvation treatment, the first solvation shell is not completely considered. The more accurate scheme should consider at least the first solvation shell completely. This could be undertaken through a thorough investigation of the clustering of protonated, hydrogenated, or electronized solvent. Such a task is beyond the scope of this work and it is devoted to the forthcoming works. For the moment, we assumed that these relative errors are constant in all media. Thus, all calculated ΔfGsolv∘(H•) and ΔfGsolv∘(e-) were shifted by an amount of the percentage error and presented in Table 1.Table 1
Solvation free energies of hydrogen atom (ΔfGsolv∘(H•)), proton (ΔfGsolv∘(H+)), and electron (ΔfGsolv∘(e-)) in kJ/mol, at 6-311++G(d,p) in various media. ε and μ are, respectively, the dielectric constant and the dipole moment of the considered media.
Solvent
ε
μ
Δ
f
G
solv
∘
(
H
•
)
Δ
f
G
solv
∘
(
H
+
)
Δ
f
G
solv
∘
(
e
-
)
Cyclohexane
2.02
0.00
−20
−879.6
−56.4
Benzene
2.28
0.00
−20
−1076.0
−69.6
Methanol
32.7
1.71
−20
−1119.7
−140.6
DMSO
46.7
3.90
−20
−1199.3
−141.2
Water
78.4
1.80
27.7
−1104.5
−156.8


## 3.2. PCET and PT Mechanisms via SOMOs, NPAs Charges, and Spin Densities Analysis

As a proton and an electron transfer from a donor to an acceptor site in a molecular system, the charge, spin density, and molecular orbitals (MOs) localized on the donor and acceptors sites change simultaneously. Thus, the calculation of the charge, spin density, and plots of MOs along a reaction coordinate (transferring hydrogen or proton) provides useful insight in understanding the PCET mechanism.SOMOs distributions on the equilibrium structures resulting from the reaction of 3,4-DHPPA withOH• in all media show that the radical is entirely delocalised on the phenoxyl moiety (see Figure 7) and could lead to the conclusion that the scavenging of OH• by 3,4-DHPPA is governed by a hydrogen atom transfer (PCET mechanism) whatever the media.Figure 7
SOMO distribution on the equilibrium structure of the supermolecule resulting from the reaction of 3,4-DHPPA withOH• radical: a B3lyp/6-31+G(d). SOMO is entirely delocalized on the phenoxyl moiety, suggesting hydrogen atom transform from the acid to the radical.NPAs and spin densities computations were performed at the same level of computations as that of geometry optimization. Computed NPA charges for products of the reaction illustrated by (17) are neutral, while spin densities are, respectively, 1.00 and 0.00 for the phenoxyl moiety and the water molecule formed, whatever the medium. In addition, the computation of such charges on reactants, products, and transition state of this reaction shows that the transferred particle has a positive charge with no spin associated, Table 2. Moreover, it is worth mentioning that the spin density associated with the phenoxyl moiety at the TS is increasing while that of the hydroxyl radical is decreasing, going from reactants to products. This indicates that the radical transfer is occurring. Moreover, at the TS, the negative NPA charges on OH increase with the solvent polarity while its spin density decreases with it, until spin density is no longer associated with OH in polar solvents. Note that, in such solvents, the electron spin density on OH is quite nil, while it is a unit on the phenoxyl moiety (RO), Table 2. Moreover, the electron charge on OH is more than the half of elementary electron charge. Thus, while the reaction of (17) is proceeding in polar solvents, the radical character is completely transferred from the free radical (OH•) to the phenoxyl moiety before the proton is transferred from the acid to the free radical. Therefore, in polar solvents, the scavenging of the hydroxyl radical by DHPPA is stepwise and the electron is transferred first. In the vacuum and nonpolar solvents, the electron spin densities on OH and RO are about 0.2 and 0.8, respectively. In addition, the electron charge on OH is less than the unit. Thus, the located TS for the reaction of (17) is a mixture of ET and PT characters. Therefore, a clear distinction of ET-PT, PT-ET, and CPCET mechanisms for such a reaction is not ensured using NPA charges and spin densities at the TS. Parameters based on reaction free energies and activation energies could be good indicators for such a target.Table 2
NPA charges (e) and spin densities carried by elements of the reaction at the transition state. The superscripts 3 and 4 on DHPPA indicate that the hydrogen atom concerned by the abstraction is H3 and H4, respectively. The symbol X denotes the unknown transferred particle.
Media
NPA charges
Spin densities
RO
X
OH
RO
X
OH
DHPPA3
Vacuum
−0.02
0.52
−0.50
0.73
0.00
0.27
Cyclohexane
0.11
0.52
−0.63
0.82
0.00
0.18
Benzene
0.10
0.53
−0.63
0.83
0.00
0.17
Methanol
0.26
0.52
−0.78
0.96
0.00
0.04
Water
0.25
0.52
−0.79
0.98
0.00
0.02
DHPPA4
Vacuum
−0.02
0.52
−0.50
0.76
0.00
0.24
Cyclohexane
—
—
—
—
—
—
Benzene
0.12
0.53
−0.65
0.89
0.00
0.11
Methanol
0.28
0.55
−0.83
0.99
0.00
0.01
Water
0.29
0.55
−0.84
0.99
0.00
0.01Moreover, we plotted the SOMOs at the TS of this reaction (17) in different solvents. These SOMOs are delocalised entirely on the supermolecule, Figures 8 and 9. These SOMOs show that the electron and proton transfer from different orbitals on the donor to different orbitals on the acceptor. The same spin density distribution was noted in other solvents (see Figure S1 in the Supplementary Material available online at http://dx.doi.org/10.1155/2015/835707).Figure 8
SOMO distribution on a transition state structure of the reaction of D-DHPPA withOH• radical in the vacuum.Figure 9
SOMO distribution on a transition state structure of the reaction of D-DHPPA withOH• radical in water.Thus, based on the works of Mayer et al. [53] and Hammes-Schiffer and Iordanova [54], these SOMOs distributions allow us to conclude that the scavenging of hydroxyl radical by D-DHPPA is governed by a CPCET mechanism, whatever the medium. Hence, in contrast to the results published before (the scavenging of hydroxyl radical by DHPPA is governed by CPCET in nonpolar solvents and PT-ET mechanisms in polar solvents [3]), PT-ET mechanism may no longer be favored in polar solvents if coreactive species are taken into account.Regarding the reaction illustrated by (18), for which there is neither radical nor transition state involved, NPA calculations show that the resulting phenoxyl moiety bears a negative charge while the HOBr formed is neutral. This reveals a proton transfer and confirms the well-known result concerning the nature of the scavenging of a negatively charged species by a PhA.


## 3.3. Thermodynamic Investigation of PCET and PT Mechanisms

Reaction free energies relevant to different pathways related to PT-ET, ET-PT, CPCET, and PT mechanisms are based on coreaction free energies (Table3) and abstraction free energies of H•, H+, and e- (Table 4). Coreaction free energies are assumed to be the free energies of the reactions involving coreactive species and abstracted or transferred particle (H•, H+, and e-). It is noteworthy that RFE-CPCET may not be used in concert with RFE-ET(1) and RFE-PT(2) to derive the preferred mechanism (PCET? PT-ET? or ET-PT?) in a given medium. In fact, the former characterizes the one-step H-atom transfer, while the other characterizes one step in a stepwise H-atom transfer. Thus, a reliable thermodynamic way to derive the preferred mechanism in a given medium is to compare RFE-ET(1), RFE-PT(2), and the activation barrier (ΔG≠) of the CPCET mechanism.Table 3
Free energies of reactions involving coreactive species and abstracted or transferred particle (H•, H+, and e-) in kJ/mol at B3LYP/6-31+G(d).
Media
•OH + e-→ OH−
•OH + •H → H2O
•OH + H+→H2O+
OBr− + H+→ HOBr
Vacuum
−160.3
−436.0
−544.1
−1438.1
Cyclohexane
−329.0
−418.8
35.4
−536.3
Benzene
−332.2
−419.1
212.6
−327.4
Methanol
−397.7
−422.0
90.1
−180.3
DMSO
−400.1
−422.1
165.8
−98.6
Water
−387.3
−469.8
67.3
−191.3Table 4
Abstraction free energies of electron, proton, and hydrogen atom from a phenolic acid ROH in kJ/mol.
Media
ROH→ ROH+ + e-
ROH→ RO− + H+
ROH→ RO• + •H
DHPPA3
Vacuum
699.4
1339.2
270.0
Cyclohexane
604.8
488.3
254.4
Benzene
581.9
283.5
254.8
Methanol
428.7
172.9
263.6
DMSO
428.0
91.9
263.7
Water
408.1
185.8
310.8
DHPPA4
Vacuum
699.4
1355.3
288.4
Cyclohexane
604.8
500.3
264.0
Benzene
581.9
295.4
263.2
Methanol
428.7
174.4
256.3
DMSO
428.0
92.9
256.1
Water
408.1
186.5
304.6
DHPPA8
Vacuum
699.40
1363.5
280.5
Cyclohexane
604.78
506.9
257.3
Benzene
581.87
301.6
256.8
Methanol
428.71
174.9
252.7
DMSO
428.04
94.9
253.0
Water
408.1
187.4
298.5


### 3.3.1. CPCET Mechanism

Previously, studies on the CPCET mechanism at 3-OH, 4-OH, and 8-OH positions in 3,4-DHPPA have shown that BDE3 < BDE8 < BDE4 in the vacuum and nonpolar solvents [3]. This order changes in polar solvents and becomes BDE8 < BDE4 < BDE3. These behaviors were also obtained when coreactive species were taken into account, Figure 10 (or Table S1). The justification of such a result was based on the relative strengths of hydrogen bonds O4H4⋯O8, O8H8⋯O9 and C2H2⋯O8 in 3,4-DHPPA. Details were discussed in our previous manuscript [3]. Thus, it may be sufficient to study the CPCET mechanism using BDFEs or BDEs rather than RFEs, in order to derive the relative ability to abstract the H-atom through the abstraction sites of the PhA. Note that the negative values obtained for RFE-CPCET in all media are the confirmation that the transfer of the H-atom from the acid to the hydroxyl radical is always possible in such media.Figure 10
Reaction free energies for the reactionROH+OH•→RO•+H2O, in different media. ROH stands for 3,4-DHPPA, while RO• stands for its phenoxyl moiety.


### 3.3.2. PT-ET Mechanisms

Assume that the reaction of 3,4-DHPPA withOH• radical processes through a PT-ET mechanism leads to the following reaction:
(23)ROH+OH•⟶RO-+H2O+⟶RO•+H2O
All the RFE-PT(2)(OH) values obtained are positive, Table S1. Thus, such a reaction is always thermodynamically unfavorable, especially in the vacuum. This may be explained by the fact that the fragmentation of a neutral molecule into a charged species is more difficult in the vacuum than in solvent, since solvent will make the molecule more polar and thereafter easy to be fragmented in charged species. Moreover, the reaction involves the reaction of OH• with a proton H+. This reaction is usually less favorable, especially in solvent. Indeed, the free radical OH• is neutral and polar, so it cannot interact easily with a charged and nonpolar molecular system. This physical insight is confirmed through RFEs of the reactions OH•+H+→H2O+, Table 3.Furthermore, the scavenging of theOH• radical by 3,4-DHPPA through the PT-ET mechanism is done in all abstraction sites (3-OH, 4-OH, and 8-OH) with the same ability in polar solvents, Figure 11. In the vacuum and nonpolar solvents, the 3-OH site is slightly preferred to the 4-OH one, which in turn is slightly preferred to the 8-OH site. In addition, this mechanism (PT-ET) is more favorable in solvent than in the vacuum, and the more the solvent polarity, the more favorable the mechanism, Figure 12. This may be explained through its dipole moment by the fact that solvent eases charge separation and, the more the solvent polarity, the easier the charge separation. These results are consistent with those found without taking into account the coreactive species [3].Figure 11
Reaction free energies evaluating the possibility of the mechanism PT-ET in the reaction of 3,4-DHPPA withOH• in different media. These RFEs are calculated as the sum of partial reaction free energies of the following reactions: ROH→RO-+H+, OH+H+→H2O+. ROH stands for 3,4-DHPPA, while RO- stands for its deprotonated structure. 3-OH, 4-OH, and 8-OH denote the cleavage sites related, respectively, to the protons H3+, H4+, and H8+ of 3,4-DHPPA.Figure 12
Reaction free energies and activation barriers evaluating, respectively, the possibility of the mechanisms PT-ET, ET-PT, and CPCET, in the reaction of 3,4-DHPPA withOH• in different media. For CPCET, the reaction involved is ROH+OH•→RO•+H2O. For PT-ET, the partial reactions involved are ROH→RO-+H+, OH+H+→H2O+. For ET-PT, the partial reactions involved are ROH→ROH++e-, OH•+e-→OH-. ROH stands for 3,4-DHPPA.By comparing the activation barrierΔG≠ of the CPCET mechanism with RFE-PT(2), one can predict that, even kinetically, the CPCET mechanism would be preferred to the PT-ET mechanism in all media, Figure 12. The justification of this statement could also be performed by a kinetic investigation of the PT-ET mechanism which is beyond the scope of this work.


### 3.3.3. ET-PT Mechanism

Assume that the reaction of 3,4-DHPPA withOH• radical processes through a ET-PT mechanism leads to the following reaction:
(24)ROH+OH•⟶ROH++OH-⟶RO•+H2O
This reaction is thermodynamically more favorable in solvent than in the vacuum, and the more the solvent polarity, the more favorable the mechanism, Figure 12. These results are consistent with the earlier explanation relevant to the easiness of the solvent to separate charged species. They are also consistent with those found without taking into account the coreactive species [3].By comparing the activation barrierΔG≠ of the CPCET mechanism with RFE-ET(1), one can predict that, even kinetically, the CPCET mechanism would be preferred in the vacuum and nonpolar solvents compared to ET-PT mechanism which in turn would be preferred in polar solvents, Figure 12. Furthermore, in polar solvents, ET-PT mechanism competes strongly with CPCET mechanism, Figure 12. The justification of these kinetic viewpoints could also be performed by a kinetic investigation of the ET-PT mechanism which is beyond the scope of this work.Therefore, during the scavenging of the hydroxyl radical by DHPPA, the transfer of the hydrogen atom is thermodynamically concerted in all media and favoring by this way the CPCET mechanism in all media. However, in kinetic viewpoint, this mechanism could compete strongly with the ET-PT mechanism so that the latter could slightly dominate. These results which corroborate the results obtained through the SOMOs distribution on the TS structures are also in agreement with those derived in polar solvents using NPA charges and electron spin densities at TS. They also bring a clear decision on the preferred mechanism in the vacuum and nonpolar solvents. It is noteworthy that the obtained results are consistent in the vacuum and nonpolar solvents with previous results [3, 5, 6, 9] which do not involve coreactive species. The differences were noted in polar solvents. This difference could be explained by the fact that the approach free of coreactives involvement is unable to give a kinetic insight of such a reaction.


### 3.3.4. PT Mechanism

PT mechanism is involved in the scavenging of the anionOBr- by 3,4-DHPPA. This mechanism is less favorable in solvents than in the vacuum, and the more the solvent polarity increases, the less the mechanism is favored, Figure 13 (or Table S1). This result may also be explained by the easiness to separate charged species in the vacuum compared to in solvents, especially in polar solvents.Figure 13
Reaction free energies related to proton transfer (PT) involving the scavenging of the anionOBr- by 3,4-DHPPA in different media. RFE-PT is calculated as the sum of partial RFEs of the reactions ROH→RO-+H+ and OBr-+H+→HOBr.


## 3.3.1. CPCET Mechanism

Previously, studies on the CPCET mechanism at 3-OH, 4-OH, and 8-OH positions in 3,4-DHPPA have shown that BDE3 < BDE8 < BDE4 in the vacuum and nonpolar solvents [3]. This order changes in polar solvents and becomes BDE8 < BDE4 < BDE3. These behaviors were also obtained when coreactive species were taken into account, Figure 10 (or Table S1). The justification of such a result was based on the relative strengths of hydrogen bonds O4H4⋯O8, O8H8⋯O9 and C2H2⋯O8 in 3,4-DHPPA. Details were discussed in our previous manuscript [3]. Thus, it may be sufficient to study the CPCET mechanism using BDFEs or BDEs rather than RFEs, in order to derive the relative ability to abstract the H-atom through the abstraction sites of the PhA. Note that the negative values obtained for RFE-CPCET in all media are the confirmation that the transfer of the H-atom from the acid to the hydroxyl radical is always possible in such media.Figure 10
Reaction free energies for the reactionROH+OH•→RO•+H2O, in different media. ROH stands for 3,4-DHPPA, while RO• stands for its phenoxyl moiety.


## 3.3.2. PT-ET Mechanisms

Assume that the reaction of 3,4-DHPPA withOH• radical processes through a PT-ET mechanism leads to the following reaction:
(23)ROH+OH•⟶RO-+H2O+⟶RO•+H2O
All the RFE-PT(2)(OH) values obtained are positive, Table S1. Thus, such a reaction is always thermodynamically unfavorable, especially in the vacuum. This may be explained by the fact that the fragmentation of a neutral molecule into a charged species is more difficult in the vacuum than in solvent, since solvent will make the molecule more polar and thereafter easy to be fragmented in charged species. Moreover, the reaction involves the reaction of OH• with a proton H+. This reaction is usually less favorable, especially in solvent. Indeed, the free radical OH• is neutral and polar, so it cannot interact easily with a charged and nonpolar molecular system. This physical insight is confirmed through RFEs of the reactions OH•+H+→H2O+, Table 3.Furthermore, the scavenging of theOH• radical by 3,4-DHPPA through the PT-ET mechanism is done in all abstraction sites (3-OH, 4-OH, and 8-OH) with the same ability in polar solvents, Figure 11. In the vacuum and nonpolar solvents, the 3-OH site is slightly preferred to the 4-OH one, which in turn is slightly preferred to the 8-OH site. In addition, this mechanism (PT-ET) is more favorable in solvent than in the vacuum, and the more the solvent polarity, the more favorable the mechanism, Figure 12. This may be explained through its dipole moment by the fact that solvent eases charge separation and, the more the solvent polarity, the easier the charge separation. These results are consistent with those found without taking into account the coreactive species [3].Figure 11
Reaction free energies evaluating the possibility of the mechanism PT-ET in the reaction of 3,4-DHPPA withOH• in different media. These RFEs are calculated as the sum of partial reaction free energies of the following reactions: ROH→RO-+H+, OH+H+→H2O+. ROH stands for 3,4-DHPPA, while RO- stands for its deprotonated structure. 3-OH, 4-OH, and 8-OH denote the cleavage sites related, respectively, to the protons H3+, H4+, and H8+ of 3,4-DHPPA.Figure 12
Reaction free energies and activation barriers evaluating, respectively, the possibility of the mechanisms PT-ET, ET-PT, and CPCET, in the reaction of 3,4-DHPPA withOH• in different media. For CPCET, the reaction involved is ROH+OH•→RO•+H2O. For PT-ET, the partial reactions involved are ROH→RO-+H+, OH+H+→H2O+. For ET-PT, the partial reactions involved are ROH→ROH++e-, OH•+e-→OH-. ROH stands for 3,4-DHPPA.By comparing the activation barrierΔG≠ of the CPCET mechanism with RFE-PT(2), one can predict that, even kinetically, the CPCET mechanism would be preferred to the PT-ET mechanism in all media, Figure 12. The justification of this statement could also be performed by a kinetic investigation of the PT-ET mechanism which is beyond the scope of this work.


## 3.3.3. ET-PT Mechanism

Assume that the reaction of 3,4-DHPPA withOH• radical processes through a ET-PT mechanism leads to the following reaction:
(24)ROH+OH•⟶ROH++OH-⟶RO•+H2O
This reaction is thermodynamically more favorable in solvent than in the vacuum, and the more the solvent polarity, the more favorable the mechanism, Figure 12. These results are consistent with the earlier explanation relevant to the easiness of the solvent to separate charged species. They are also consistent with those found without taking into account the coreactive species [3].By comparing the activation barrierΔG≠ of the CPCET mechanism with RFE-ET(1), one can predict that, even kinetically, the CPCET mechanism would be preferred in the vacuum and nonpolar solvents compared to ET-PT mechanism which in turn would be preferred in polar solvents, Figure 12. Furthermore, in polar solvents, ET-PT mechanism competes strongly with CPCET mechanism, Figure 12. The justification of these kinetic viewpoints could also be performed by a kinetic investigation of the ET-PT mechanism which is beyond the scope of this work.Therefore, during the scavenging of the hydroxyl radical by DHPPA, the transfer of the hydrogen atom is thermodynamically concerted in all media and favoring by this way the CPCET mechanism in all media. However, in kinetic viewpoint, this mechanism could compete strongly with the ET-PT mechanism so that the latter could slightly dominate. These results which corroborate the results obtained through the SOMOs distribution on the TS structures are also in agreement with those derived in polar solvents using NPA charges and electron spin densities at TS. They also bring a clear decision on the preferred mechanism in the vacuum and nonpolar solvents. It is noteworthy that the obtained results are consistent in the vacuum and nonpolar solvents with previous results [3, 5, 6, 9] which do not involve coreactive species. The differences were noted in polar solvents. This difference could be explained by the fact that the approach free of coreactives involvement is unable to give a kinetic insight of such a reaction.


## 3.3.4. PT Mechanism

PT mechanism is involved in the scavenging of the anionOBr- by 3,4-DHPPA. This mechanism is less favorable in solvents than in the vacuum, and the more the solvent polarity increases, the less the mechanism is favored, Figure 13 (or Table S1). This result may also be explained by the easiness to separate charged species in the vacuum compared to in solvents, especially in polar solvents.Figure 13
Reaction free energies related to proton transfer (PT) involving the scavenging of the anionOBr- by 3,4-DHPPA in different media. RFE-PT is calculated as the sum of partial RFEs of the reactions ROH→RO-+H+ and OBr-+H+→HOBr.


## 3.4. Kinetic Study of the PCET Mechanism

To complete the above thermodynamic description, kinetic investigation of the PCET mechanism (17) was undertaken in various media. The media retained were vacuum, cyclohexane, benzene, methanol, and water. The basis set superposition error (BSSE) was computed for the transition state structure at the G3B3 level of computation in the vacuum. For the TS structures related to the H3-atom, H4-atom, and H8-atom abstractions in the vacuum, the obtained values were 12.6, 12.4, and 12.3kJ/mol, respectively. The same values were used in solvents, since up to now there are no possible counterpoise computations in solvent. Rate constants were computed at T=298.15 K, using TST/ST and TST/ZCT methods. The TST/ST method is the traditional TST method with κST< for the transmission coefficient and is implemented in our homemade Fortran program. The computed activation free energies (ΔG≠) were achieved including zero-point energies.The rate constants obtained from both described methods are quite consistent (see Figure14) and lead to the same information that rate constants decrease as the solvent polarity increases. The result is also consistent with the increasing distance d(On⋯Hn) with the solvent polarity, n=3,4. It is also in good agreement with the decreasing frequency ν(On⋯Hn). Indeed, the longer distance matches a lower frequency, since a longer distance corresponds to a weaker bond. The discrepancy observed in the vacuum between both methods of rate constants calculation could be due to the large Wigner tunnelling effects in such a medium. This shows that the tunnelling effects on hydrogen atom transfer could be more important in the vacuum than in a solvent. Indeed, the TS vibrational frequency in the vacuum is -2429cm-1, while in a solvent the value is about -1500cm-1, which is lower than the former.Figure 14
Rate constants for the scavenging of the free radicalOH• by the 3,4-DHPPA, computed using TST/ST and TST/ZCT methods in different media. The attacked site considered is O3–H3.This result is consistent with the increasing activation barrier with solvent polarity (Table5), suggesting that the higher the solvent polarity, the slower the CPCET mechanism. This result is not surprising, since it is shown early that solvent makes the one-step H-atom transfer harder.Table 5
Parameters characterizing the O–H bond cleavage in selected media.d(On⋯Hn) in Å is the distance from the oxygen On of the phenoxyl moiety to the moved hydrogen atom Hn. ν(On⋯Hn) in cm-1 is the corresponding frequency associated with the movement of this hydrogen atom, while kTST/ST and kTST/ZCT in cm3⋅molecule-1⋅s-1 are the associated rate constants using, respectively, TST/ST and TST/ZCT methods. ΔG≠ in kJ/mol is the activation free energy, and n = 3, 4. NA stands for not applicable.
Media
d
ν
Δ
G
≠
k
TST/ST
k
TST/ZCT
Related to theO3–H3 attacked site
Vacuum
0.969
3769
31.2
6.8× 10-14
5.0× 10-13
Cyclohexane
0.975
3645
47.0
1.6× 10-14
3.5× 10-14
Benzene
0.976
3625
53.8
5.1× 10-15
7.1× 10-15
Methanol
0.990
3331
71.8
3.6× 10-15
2.5× 10-16
Water
0.991
3326
78.3
4.8× 10-16
1.6× 10-16
Related to theO4–H4 attacked site
Vacuum
0.964
3710
68.1
1.9× 10-17
6.2× 10-16
Cyclohexane
0.975
3661
—
—
—
Benzene
0.976
3652
78.6
2.9× 10-18
NA
Methanol
0.983
3493
93.3
1.7× 10-19
6.1× 10-18
Water
0.985
3457
90.9
—
—Moreover, whatever the medium, the hydrogen atomH3 is transferred more rapidly than the hydrogen atom H4, since the related rate constant of H3 transfer is higher than that related to H4 transfer in a given medium, Table 5. In addition, the activation free energy and rate constants using TST/ST and TST/ZCT methods for the O8–H8 bond cleavage in the vacuum are 70.1 kJ/mol, 1.3×10-17 cm3·molecule−1·s−1, and 4.2×10-16 cm3·molecule−1·s−1, respectively. Thus, the rate constant related to H4 transfer is higher than that related to H8 in the vacuum.


## 4. Conclusion

In this work, we investigated the reaction of 3,4-dihydroxyphenylpyruvic acid (3,4-DHPPA) with hydroxyl radical (OH•) and bromite ion (OBr-) in various media (vacuum, benzene, cyclohexane, methanol, DMSO, and water). The goal of this work was to distinguish the concerted proton-coupled electron transfer (CPCET) from sequential one, as well as proton transfer-electron transfer (PT-ET) from electron transfer-proton transfer (ET-PT) in the O–H bond cleavage reactions in various media.As a result, we found that the scavenging of hydroxyl radical by 3,4-DHPPA is thermodynamically governed by a one-step hydrogen atom transfer (CPCET) from the acid to the radical in all media. In kinetic viewpoint, CPCET still dominates in the vacuum and in nonpolar solvents, but in polar solvents it could compete strongly with the ET-PT mechanism so that the latter could slightly dominate. In kinetic viewpoint, this mechanism could compete strongly with the ET-PT mechanism so that the latter could slightly dominate. In addition, the rate constant of the CPCET mechanism decreases going from the vacuum to the solvent, and the higher the solvent polarity is, the lower the rate constant is. Thus, in the vacuum and nonpolar solvents, results that take into account the coreactive species are similar to those that do not. The differences were noted in polar solvents. For the reactionROH+OBr-→RO-+HOBr, one can retain that, being more favored in the vacuum than in solvent, it proceeds through a proton transfer from the acid to the reactive oxygen species, and the more the solvent polarity increases, the less the mechanism is favored.Therefore, studying the O–H bond cleavage subsequent to the scavenging of reactive oxygen species by a phenolic acid through simple thermodynamic parameters such as bond dissociation free energy, proton affinity, or ionization potential is only satisfactory in the thermodynamic viewpoint. Reaction free energy (RFE) parameters could be reliable to have both thermodynamic and kinetic viewpoints in such a reaction. Moreover, although the RFE may not be necessary for the concerted PCET mechanism, it may be very useful and may be a deciding factor for the investigations of sequential proton-coupled electron transfer and proton transfer mechanisms.



---

*Source: 835707-2015-01-19.xml*