Full Paper
Theoretical studies
larization functions, on the previously optimized geometries.
The Orca program[21] was used for this. All the calculations
have been carried out in THF solution, using either PCM[22]
(low level geometry optimizations) or COSMO[23] (energy refine-
ments).
Some of the obtained results seem to be contradictory if we
consider the existence of only one of the proposed mecha-
nisms (Scheme 2). For instance, the unexpected reluctance of
p-methoxy or p-dimethylaminophenyl ethers 1d and 1e, re-
spectively, to undergo Wittig rearrangement could lead us to
think that the anionic pathway is operating. However, the dif-
ferent behavior of the regioisomeric o-methoxyphenyl ether
1g is in disagreement with this assumption. Similarly, the fact
that 1o also undergoes Wittig rearrangement, seems to be
against the anionic proposal that would involve the formation
of a highly strained bridged intermediate. Moreover, the
higher tendency of 4-biphenyl ether 1k, as well as naphthyl
ethers 1m and 1n to undergo Wittig rearrangement, could
not be easily explained through the radical pathway. To shed
more light on the mechanism of the Wittig rearrangement in
aryl benzyl ethers and understand the results obtained for the
different substrates in Table 1, we decided to use DFT to
model two different paths, one of radical nature and the other
of anionic character (Figure 1).
Computational results
The radical pathway would involve the homolitic cleavage of
the ether ArÀO bond in I, leading to the aryl radical II and ben-
zyloxy radical III (Figure 1). Recombination of these two radi-
cals through a barrierless transition state ts(III–IV) provides the
reorganized product IV. An alternative anionic path would in-
volve an attack of the benzylic anion on the aryl group that
would be concerted with the cleavage of the ArÀO bond, di-
rectly leading to the reorganized product IV. For anionic char-
acter, we also proposed a third option where the attack of the
anion on the aryl group would lead to an oxirane intermediate
V that would then undergo ring-opening to produce IV; how-
ever, we were not able to locate such an intermediate, at least
for when Ar is Ph.
In Table 2, we show the relative Gibbs free energies of the
different structures found along these two paths for some rep-
resentative ethers from Table 1. We found that there is not
a generally preferred mechanism, and that in some instances,
the radical and anionic concerted paths can be considered to
be competitive and can coexist. If we compare the activation
energies for ts(I–II) (the later collapse of the two radicals is
found to be a barrierless process) and ts(I–IV), we find that the
concerted path is favored when Ar is Ph (1a), 4-ClC6H4 (1h), 1-
naphthyl (1m), 2-naphthyl (1n), and 6H-benzo[c]chromene
(1o) (entries 1 and 6–9, respectively, underlined in Table 2),
and slightly preferred when Ar is 4-MeC6H4 (1b) (entry 2). In
contrast, ts(I–II) is preferred with more electron-rich aryl
groups, such as when Ar is 4-MeOC6H4 (1d), 4-Me2NC6H4 (1e)
and 2-MeOC6H4 (1g) (entries 3–5, respectively). In these cases
(italicized in Table 2), however, (as with Ar=4-MeC6H4 (1b),
entry 2) both paths are competitive, with differences between
Figure 1. Mechanisms for the Wittig rearrangement considered in the DFT
study.
the two rate-determining steps of less than 1.0 kcalmolÀ1
.
The preference for a radical or anionic concerted mechanism
seems to be electronic in nature, because we found that one
of the largest differences between the activation energy of the
two alternative paths is found in 6H-benzo[c]chromene (1o).
For this ether, the preferred anionic concerted mechanism
could be expected to be prevented by the rigidity of the struc-
ture. The observed trends in reactivity are well replicated in
our calculations: 1o, naphthyl ethers 1m and 1n, and 4-
chlorophenyl ether 1h rearrange faster than the parent phenyl
ether 1a, whereas 4-methoxyphenyl (1d) or 4-aminomethyl
(1e) ethers do not rearrange at all. The fact that we do not
find a rearrangement product for 4-methoxyphenyl ether 1d,
although 2-methoxyphenyl ether 1g rearranges readily, is puz-
zling from a purely electronic point of view, because the differ-
ence in activation of the phenyl ring should not be very large
between the two systems. What is found here is that the ge-
ometry of the substrate is key in determining the reactivity.
The barriers for ts(I–II) and ts(I–IV) are comparable for when
Computational methods
We used DFT to locate the stationary points depicted in the
first two reaction paths in Figure 1. We carried out geometry
optimizations with the B3LYP hybrid exchange-correlation
functional with the double-z 6-31G* basis set, as implemented
in the Gausian 09 code.[18] We confirmed that the wavefunc-
tions are stable and used harmonic analysis to characterize
these structures as either minima or transition states in the po-
tential energy surface. These harmonic frequencies have been
used to calculate the thermal corrections needed to obtain
free energies. In some cases, we ran IRC (intrinsic reaction co-
ordinate) calculations to unambiguously connect transition
states and the corresponding minima. For obtaining more ac-
curate electronic energies, we used the double hybrid B2PLYP
functional with Grimme’s dispersion correction (B2PLYP-D3),[19]
with the triple-z basis def2-tzvpp,[20] which includes extra po-
&
&
Chem. Eur. J. 2016, 22, 1 – 12
4
ꢂ 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim
ÝÝ These are not the final page numbers!