Density Functional Theory Studies on the Addition and Abstraction Reactions of OH Radical with Benzoate Anion

The addition and abstraction reactions of OH radical with benzoate anion are investigated by density functional theory calculations that include solvent effects using UB3LYP, UCAM-B3LYP, UmPW1PW91 and UM06-2X functionals with the 6-311++G(2d,2p) basis set. Geometry optimizations of the reactants, products and transition state species are performed for the possible reaction paths. For the addition reactions, those targeting the ipso-, ortho-, metaand para-carbons are predicted to be exoergic. The H-atom abstraction reactions from ortho, meta and para positions are also predicted to be exoergic. On the basis of the rate constants calculated by means of the transition state theory, the H-atom abstraction reaction from the ortho position is determined to be the favored path followed by the ortho OH addition reaction.


Introduction
Sodium salicylate is known as a good phosphor under the x-ray and vacuum ultraviolet excitation [1].Since its proposal as a sensitive dosimeter [2], the salicylate has been also utilized as a probe to measure OH radical generation in aqueous solutions where salicylate anion is generated by the reaction of benzoate anion with OH radical.This method is widely used in radiochemical, sonochemical and photocatalytic reactions [3] as well as the method using the 2-hydroxyterephthalate dianion [4].Peräkylä et al. [5] have theoretically studied the reactivity of benzoate towards OH addition reaction at the ROHF/6-31G(d) and ROMP2/6-31G(d) levels including solvent effects on the basis of the qualitative frontier molecular orbital theory and compared the stability of the isomeric OH adducts of benzoate.They revealed the correlation between the kinetic reactivity and the thermodynamic stabilities of the addition products.However, no theoretical study of the detailed kinetics and thermochemistry of the reactions has been reported.As concerns the water environmental chemistry, the OH radical initiated reactions of aromatic compounds have been studied [6,7].The rate constants, the branching ratios and product distribution were in good agreement with the experimental values.In the present study, we have focused on the kinetics and thermochemistry of the primary steps; namely, OH radical addition to the benzoate anion and H-atom abstraction from the benzoate anion (Scheme 1).

Computational Methods
The equilibrium geometries of the reactants, transition states, products and complexes were optimized using the density functional theory (DFT) method.Hybrid GGA functionals B3LYP [8] and mPW1PW91 [9], hybrid meta-GGA functional M06-2X [10] and range-separated hybrid GGA functional CAM-B3LYP [11] were employed with the 6-311++G(2d,2p) basis set.As compared with the standard B3LYP functional, the mPW1PW91 and CAM-B3LYP improve long-range behavior and the M06-2X functional has been recommended for thermochemistry and kinetics by the authors.Harmonic vibrational frequencies were calculated to confirm the predicted structures as local minima or transition states (one imaginary frequency) and elucidate zero-point vibra- tional energy corrections (ZPE).The obtained transition states were confirmed as those connecting the investigational species by a calculation of the subsequent intrinsic reaction coordinates.The solvation effects were included using a polarizable continuum method (IEF-PCM).All calculations were performed using Gaussian 09 [12].The spin-squared values were checked and the deviations from the doublet value of <s 2 > = 0.75 were lower than 4.8%.Therefore, the spin contamination is negligible for all the studied radical species.
The rate constants for the addition and abstraction reactions were estimated using the conventional thermodynamic formulation of the transition state theory.The rate constant k is given by where k B is the Boltzmann constant, T is the temperature, h is Planck's constant, ΔG ≠ is the Gibbs free energy of activation and κ is the transmission coefficient.The tunneling effect was approximated using the Wigner correction [13].

OH Addition Reaction
There are five possible sites for OH addition, i.e. carbonyl, ipso, ortho, meta and para carbons of benzoate anion.The optimized structures of benzoate anion and addition products for the addition reaction are shown in Figure 1.For the carbonyl, ipso, meta and para addition products, respectively two minima were found on the CCOH torsion potentials, where the O-H bonds point toward or away from the ring.For the ortho addition product, two conformers are present, where the O-H bond points toward the ring or toward the COO group to form the H-bond.Some optimized geometrical parameters are also given in Figure 1.Table 1 lists the relative electronic energies, including the zero-point vibrational energies (Δ r E), relative enthalpies at 298 K (Δ r H) and relative Gibbs energies at 298 K (Δ r G).The Δ r G for the carbonyl addition reactions are calculated to be positive while others to be negative.The spin density of the carbonyl addition product is localized on the oxygen atoms of the OH added COO group, whereas those of the other products are delocalized to gain the resonance stabilization energies.The ortho 1 addition product was calculated to be the thermodynamiccally most stable one due to an intramolecular hydrogen bond with CO•••HO distances in the range of 1.707 (UmPW1PW91) to 1.826 (UM06-2X) Å.The ipso 1 addition product also forms an intramolecular hydrogen bond with CO•••HO distances in the range of 1.758 (UmPW1PW91) to 1.816 (UM06-2X) Å. Transition states TS carbonyl1 , TS carbonyl2 , TS ipso1 , TS ipso2 , TS ortho1 , TS ortho2 , TS meta and TS para correspond to the formation of addition products carbonyl 1 , carbonyl 2 , ipso 1 , ipso 2 , ortho 1 , ortho 2 , meta 1 and para 1 , respectively.IRC calculations revealed that transition states TS ortho1 and TS ipso1 are connected with complex 1 , and TS ortho2 , TS meta and TS para are connected with complex 2 , complex 3 and complex 4 , respectively.
The imaginary frequencies nsition states and distance ratios of r C•••OH /r C-OH are summarized in Table 2.The C•••OH distances in the transition states for the ipso, ortho, meta and para addition reactions are elongated by over 40% as compared with those of the addition products, indicative of an early Optimized structures of the transition states for the OH  quilibrium disnces in the addition products, indicative of a late transi-for the carbonyl addition reaction, the C•••OH distances are approximately 20% longer than the e ta tion state.The ΔG ≠ value for the addition of OH to carbonyl carbon is found to be much higher than those of others.The ΔG ≠ values for the formation of ortho 1 and para 1 are comparable and the lowest, though it depends on the functionals.Therefore, ortho and para positions should be the kinetically most favored sites for addition of OH radical.In complex 1 , the O-H bond length is elongated by 4.0% (UB3LYP) to 4.6% (UM06-2X) as compared with that of the bare OH radical.The CO bond interacting with the OH radical is also elongated while the other CO bond is contracted as compared with that of the benzoate anion.The Δ r G values for the formation of the complex 1 are calculated to a range from −1.3 (UM06-2X) tios Table 2. Calculated imaginary frequencies, CO distance ra and thermochemical parameters (kcal•mol −1 ) of the transition states and rate constants for addition reactions.

Stationary point
Functional  calculated to be negative while Δ r G values to be positive.Therefore, at 298 K the formation of complex 2 , complex 3 and complex 4 are negligible.Complex 3 was unable to be located on the potential energy surfaces at the UM06-2X level.

H-Atom Abstraction Reaction
Optimized structures of the relevant stationary points for H-atom abstraction from the b nzoate an on by the e i OH 11++G(2d, 2p) lists the ble 3. The rate constant for the abstractio radical calculated at the UB3LYP/6-3 level of theory are shown in Figure 3. Table 3 ΔE ≠ , ΔH ≠ and ΔG ≠ values for the H-atom abstraction reactions.For H-atom abstractions from the ortho and meta positions, the transition states were optimized by starting from the three initial geometries where the OH molecular axis was perpendicular to the molecular plane of the benzoate anion or the OH molecular axis was in the plane of the benzoate anion with either the H-atom or O-atom of the OH radical pointing to the carbonyl oxygen atom.For H-atom abstraction from the para position, the calculation was started from the two initial geometries where the OH molecular axis was perpendicular to or in the molecular plane of the benzoate anion.TS ab_ortho1 was calculated to be planar with the OH•••OC distance in the range from 1.650 (mPW1PW91) to 1.704 (B3LYP) Å. TS ab_ortho2 was calculated to be nonplanar with the O-H bond pointing away from the carbonyl group.In the TS structures for the meta and para H-atom abstraction reactions, the OH molecular axis was perpendicular to the molecular plane of the benzoate anion.The breaking C-H bonds are elongated by about 17, 14, 12 and 13% as compared with the equilibrium C-H bond of the benzoate anion while the forming O-H bonds are elongated by about 28, 35, 37 and 35% as compared with the O-H bond length of H 2 O for TS ab_ortho1 , TS ab_ortho2 , TS ab_meta and TS ab_para , respectively for the UB3LYP calculation.The reactant-like geometry is characteristic of an early transition state.The ΔG ≠ values of TS ab_ortho1 are smaller than those of TS ab_ortho2 owing to hydrogen bonding stabilization.
The calculated rate constants for H-atom abstraction are listed in Ta n reaction via TS ab_ortho1 is one to three orders of magnitude larger than those for the other abstraction reactions.Total rate constants for the abstraction reaction were obtained to be 7.0 × 10 −12 , 4.0 × 10 −12 , 1.4 × 10 −12 and 1.3 × 10 −13 cm 3 mol −1 s −1 using UB3LYP, UmPW1-PW91, UCAM-B3LYP and UM06-2X functionals, respectively.Figure 4 gives the potential energy profile for the reaction pathway.The product distribution was obtained for the UB3LYP calculations as follows: carbonyladdition:ipso-addi-tion:ortho-addition:meta-addition:para -addition:ortho-abstraction:meta-abstraction:para-abstrac tion = 0.0:0.4:19.9:8.6:13.7:50.7:4.7:2.1.Comparing the rate constants for the addition reactions versus the abstraction reaction, it can be concluded that the latter is the slightly favored reaction channel for the reaction of OH radical with the benzoate anion.

Conclusion
The kinetics and t abstraction reactions of OH radical with benzoate anion have been investigated using DFT.On the basis of the rate constants calculated by means of the transition state theory, the H-atom abstraction reaction from the ortho position is determined to be the favored path followed by the ortho OH addition reaction.Hydrogen bonds present

Figure 1 .
Figure 1.Optimized structures of the reactants and OH ad

Figure 2 .
Figure 2. Optimized structures of the transition states for OH addition reaction and complexes calculated at the UB3LYP/6-311++G(2d,2p) level of theory.

Figure 3 .
Figure 3. Optimized structures of the transition states for H-atom abstraction reaction and abstraction products calculated at the UB3LYP/6-311++G(2d,2p) level of theory.

Figure 4 .
Figure 4. Potential energy diagram for the addition and H-atom abstraction reactions of OH + benzoate anion calcula at the UB3LYP/6-311++G(2d,2p) level of theory. ted