Quasi-classical
molecular dynamics (MD) simulations were carried
out to study the mechanism of iron porphyrin-catalyzed hydroxylation
of ethylbenzene. The hydrogen atom abstraction from ethylbenzene by
iron-oxo species is the rate-determining step, which generates the
radical pair of iron-hydroxo species and the benzylic radical. In
the subsequent radical rebound step, the iron-hydroxo species and
benzylic radical recombine to form the hydroxylated product, which
is barrierless on the doublet energy surface. In the gas-phase quasi-classical
MD study on the doublet energy surface, 45% of the reactive trajectories
lead directly to the hydroxylated product, and this increases to 56%
in implicit solvent model simulations. The percentage of reactive
trajectories leading to the separated radical pair is 98–100%
on high-spin (quartet/sextet) energy surfaces. The low-spin state
reactivity dominates in the hydroxylation of ethylbenzene, which is
dynamically both concerted and stepwise, since the time gap between
C–H bond cleavage and C–O bond formation ranges from
41 to 619 fs. By contrast, the high-spin state catalysis is an energetically
stepwise process, which has a negligible contribution to the formation
of hydroxylation products.