Journal of the American Chemical Society
Article
outset, we recognized that effective computational simulation
of the crystalline lattice would be critical in achieving
meaningful results. While gas-phase and simple solvation
models would be unlikely to accurately simulate the crystalline
environment, full density functional theory (DFT) treatment
of the crystalline lattice would be untenable due to prohibitive
computational cost. Therefore, the ONIOM57 approach within
the hybrid QM/MM58 method was employed, in order to
estimate the intermolecular interactions within the crystal
lattice environment.59 The central part consisted of a single
radical pair computed with the unrestricted open-shell
dispersion-corrected ωB97X-D/6-311G(d,p) level of theory.60
The surrounding shell of molecules, extracted from the
experimentally determined crystal structure, was computed
with the electronically embedded UFF force field.61
Using this methodology, we investigated the recombination
barriers for radical pairs 10 in the solid state (Figure 9).
Radical pairs 10 were generated computationally by removing
the carbonyls from ketones 11 and then performing an energy
minimization calculation. Interestingly, the success of the
ground-state optimization for radical fragments was critically
dependent on embedding them in their respective crystal
cavities. In the gas phase, the radical pair would spontaneously
recombine to form the C−C bond. However, when conforma-
tionally restrained by the crystal cavity, the radical pair was an
energetic minimum and existed as a discrete species.
Transition states for the recombination reaction were then
identified for each of the radical pairs in their respective
lattices. Radical pairs 10a and 10c, which lack ortho
substitution, displayed the lowest barriers of 4.8 and 1.7
kcal/mol, respectively. Recombination barriers for ortho-
functionalized radical pairs 10d−f were uniformly higher,
with values between 6.1 and 6.8 kcal/mol. Given the successful
calculation of transition states for C−C bond formation from
radical pairs 10, we were additionally interested in probing the
competing disproportionation process. The results of these
calculations suggest that hydrogen atom tunneling may be
responsible for the observed disproportionation (for further
After determining the recombination barriers for radical
pairs 10, we turned our attention to radical pair 22d (arising
from decarbonylation of ketone 28d) (Figure 10). For this
calculation, we used the coordinates derived from the reactive
polymorph of 28d. While radical pairs 10 all displayed small
but non-negligible barriers to C−C bond formation, the
corresponding process for radical pair 22d was barrierless. In
fact, all attempts at energy minimization of 22d resulted in
radical recombination to form the C−C bond. This result was
consistent with experimental results described previously,
which demonstrated that photodecarbonylation of ketone
28d (reactive polymorph) led to efficient C−C bond
formation (Figure 6a).
Figure 10. Computational investigation of radical pair 22d.
reaction, realized by the presence of an activation barrier.
Comparatively, the reactive polymorph of ketone 28d likely
undergoes a barrierless recombination due to a favorable
ground state conformation. This is reflected both by the
barrierless recombination observed in computational analysis
and by the high-yielding formation of photoproduct 30d.
Additionally, the magnitude of the barrier to radical
recombination for radical pairs derived from ketones 11 is
related to the presence of an ortho-substituent. The computa-
tional results clearly show a higher activation energy for
recombination when an ortho-substituent is present (6.1−6.8
kcal/mol) and a diminished barrier for radical pairs that lack
ortho substitution (1.7−4.8 kcal/mol).62
Summary of the Optimization and Mechanistic
Analysis of the Solid-State Photodecarbonylation
Reaction. The optimization and analysis en route to the
bis(cyclotryptamine) alkaloids gave insight into the role of
substrate conformation on the success or failure of the solid-
state photodecarbonylation reaction. It was found that the
performance of the transformation was highly sensitive to
remote substituent effects. While the reaction was low-yielding
for N-methyl ketones 11d−f, the solid-state photodecarbony-
lation reaction of N-acyl ketones 28d and 29d,e gave up to
61% yield of the desired vicinal quaternary stereocenter-
containing product. We hypothesize that changes to the solid-
state conformation of the substrate underpin these dramatic
substituent effects on photodecarbonylation. This hypothesis
was further supported computationally through study of the
C−C bond forming events for radical pairs 10a, 10c−f, and
22d. Additionally, a key stereoelectronic relationship between
the conformation of the substrate and the photolability of the
ketone carbonyl was identified. This led to dramatically
different reactivity of two polymorphic forms of 28d, the
failed reaction of 28e, and efficient decarbonylative C−C bond
formation from 29d,e. Ultimately, this optimization employed
crystal engineering of the solid-state conformation to improve
On the basis of the computational results depicted in Figures
9 and 10 and the experimental results previously discussed, we
were able to draw a number of conclusions about the key
photodecarbonylation reaction. Conformational restriction can
sometimes lead to activation barriers for recombination. In the
photodecarbonylation reaction of ketones 11a and 11c−f, the
computationally generated radical pairs 10a and 10c−f were
not energetic minima unless they were confined by the
crystalline lattice. This is consistent with our hypothesis that in
the solid state, unfavorable ground-state substrate conforma-
tions can lead to repulsive interactions in the course of the
4049
J. Am. Chem. Soc. 2021, 143, 4043−4054