O-Neophyl Rearrangement of 1,1-Diarylalkoxyl Radicals
one 4b (4.6 kcal‚mol-1). A similar trend is observed, although,
as expected, to a lesser degree, for the free-energies of the first
transition-state TSa-b relative to the formation of the bridged
radical intermediate from the 1,1-diarylalkoxyl radical. In
contrast, the second transition-state TSb-c, relative to the
transformation of the radical intermediate into the corresponding
final carbon centered radical, appears to be much more sensitive
to electronic effects, which spans about 5 kcal‚mol-1 upon going
from the OMe to the NO2 group. Interestingly, this large
sensitivity changes the relative importance of the two transition-
states in affecting the observed reaction rate. Indeed, according
to the steady-state approximation, the reaction rate is mostly
affected by the transition-state which has the higher free-energy,
no matter whether it is the first or the second transition-state.10
The second transition-state, TSb-c, has a higher free-energy than
the first one, TSa-b, in the cases of X ) MeO and H and has
comparable free-energy in the case of X ) CF3. Its free-energy
decreases below the level of TSa-b in the case of X ) NO2.
The final carbon centered radical is also largely sensitive to
electronic effects; its free-energy decreases by 5.4 kcal‚mol-1
upon going from the OMe to the NO2 group. The main
contribution to the energetics of the process can be reasonably
ascribed to the cleavage of a CAr-C bond in the starting 1,1-
diarylalkoxyl radical (1a-4a) that is accompanied by the
formation of a CAr-O bond in the rearranged carbon centered
radical (1c-4c). Along this line, relative gas-phase CAr-C and
CAr-O BDEs of para-substituted benzyl alcohols and anisoles,
collected in Table 2, can help to shed light on the relative
importance of the two processes. The data indeed show that,
although all substituents determine a small but comparable
increase in the CAr-C BDE for the benzyl alcohol series,
opposite effects are instead observed in the anisole series. The
CAr-O BDE slightly decreases in the presence of the OMe ring
substituent and progressively increases upon going to the
electron-withdrawing CF3 and NO2 substituents. Very impor-
tantly, the latter trend is analogous to that observed for the
calculated reaction free-energy differences, ∆G°, for the rear-
rangement of radicals 1a-4a (Table 1) discussed above, which
indicates that in the rearrangement the effect of ring substituents
on the strength of the newly formed CAr-O bond plays an
important role. Along this line, the observations that the free-
energy differences calculated for the bridged radical intermediate
and the activation free-energies relative to its conversion into a
carbon centered radical, ∆G°TSb-c, increase in the presence of
electron-releasing substituents and decrease in the presence of
electron-withdrawing ones (see above) can be interpreted
similarly in terms of the substituent effect on the strength of
the CAr-O bond.
reaction path. The present study, moreover, indicates that, among
the factors that govern the rearrangement, the strength of the
CAr-O bond that is formed during the rearrangement appears
to play an important role.
Experimental Section
Computational Details. Hybrid DFT calculations have been
carried out at the UB3LYP/6-31G(d) level of theory using the
Gaussian 03 program package.12 All species have been fully
geometry optimized in the gas-phase and characterized as minima
or transition-states by frequency calculations; the Cartesian coor-
dinates are supplied in Supporting Information. The calculated spin-
squared expectation values (<S2>) were, after spin annihilation,
e0.751 in all cases, which is in good agreement with the
theoretically expected value of 0.75 for a pure doublet state.
Transition-states have been located by the quadratic synchronous
transit method (QST2). The eigenvector relative to the imaginary
frequency of each transition-state has been carefully inspected with
the help of the program GaussView13 to check that it corresponded
to the expected reaction coordinate. IRC calculations have been
also carried out to confirm that the located transition-state actually
connected the reactants and products from which it had been
generated by the QST2 algorithm. Zero-point energies and gas-
phase thermodynamic data at 298.15 K and 1 atm have been
calculated by the standard harmonic approximation without any
correction to vibrational frequencies.
Gas-phase BDEs at 298.15 K and 1 atm relative to the processes
in Scheme 4, with X ) OCH3, H, CF3, and NO2, were also
calculated by geometry optimization and frequency calculations of
the involved species at the UB3LYP/6-31G(d) level of theory [BDE
) H° (p-XC6H4 ) + H° (Y•) - H° (p-XC6H4Y)].
•
298
298
298
Acknowledgment. Financial support from the Ministero
dell’Istruzione dell’Universita` e della Ricerca (MIUR) is grate-
fully acknowledged.
Supporting Information Available: Results and experimental
details of product and time-resolved studies for radical 9a. Cartesian
coordinates, potential energy, and thermodynamic data for all of
the stationary points of the reaction of radicals 1a-6a calculated
at the UB3LYP/6-31G(d) level of theory. This material is available
JO070125N
(12) Gaussian 03, Revision B.01, Frisch, M. J.; Trucks, G. W.; Schlegel,
H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A.,
Jr.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.;
Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.;
Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda,
R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai,
H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken,
V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.;
Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.;
Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski,
V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D.
K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui,
Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.;
Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith,
T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.;
Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; and
Pople, J. A.; Gaussian, Inc.: Wallingford CT, 2004.
In conclusion, the results presented above show that DFT
computational methods based on the hybrid B3LYP functional
and on a basis set of moderate size such as 6-31G(d) can be
conveniently used for the determination of kinetic data relative
to neutral free-radical reactions. This is evidenced by the close
adherence of the computational gas-phase rate constants to the
experimental ones obtained in solution. The results point out
to the reliability of the computed gas-phase mechanism of
reaction. The bridged 1-oxaspiro[2,5]octadienyl radical has
generally been shown to be a minimum along the reaction
coordinate, which corresponds to a reactive intermediate. This
is in contrast with previous hypotheses involving either a long-
lived intermediate or the absence of this intermediate along the
(13) Dennington, R., II; Keith, T.; Millam, J.; Eppinnett, K.; Hovell, W.
L.; and Gilliland, R.; GaussView, Version 3.09, Semichem, Inc.: Shawnee
Mission, KS, 2003.
J. Org. Chem, Vol. 72, No. 12, 2007 4519