Nonradical Zn−Barbier Synthesis of Amino Alcohols
A R T I C L E S
single-electron transfer (SET) mechanism. In particular the
COOEt and CF3-groups provide an interesting comparison. The
two groups have similar σ and σ+ values, but differ strongly
for σ•; the ester can stabilize a radical through resonance,
whereas the CF3 group cannot. Thus, the two groups should
induce very different rates in a radical reaction, but in the current
study, the rates of the two substrates are similar. In addition,
none of the groups employed here has a negative σ• value, but
still we see that electron-donating substituents retard the reaction.
Going back to the Hammett plot using the regular σ values,
the weakly positive F is consistent with a nucleophilic attack
on the imine and indicates that the electrophilicity of the imine
has a higher influence on the rate than the ability of the imine
nitrogen to coordinate to the metal.
The Hammett study strongly indicates that the allylic bromide
initially reacts with zinc to form an organometallic reagent,
which subsequently reacts with the imine. This opens up the
possibility to study the reaction computationally with a fair
accuracy. The SET mechanism would in principle require high-
level methods able to calculate accurate structures and energies
for open shell species on a metal surface. Such calculations are
possible today, but well beyond our resources. However, with
the current data in hand, we estimate that the reaction path can
be studied starting from the isolated organometallic reagent. We
have previously been involved in several computational studies
on the addition of isolated zinc reagents to carbonyl com-
pounds.20 The goal of the current study is to understand the
observed diastereoselectivity in the reaction.
Computational Study. We have initiated the study using
small model complexes, first to enable us to select a suitable
computational methodology, then to screen all potential paths
for the reaction. Subsequently, the most promising reaction paths
have been studied using the full substrate 3a, to give final
predicted selectivities. All calculations have been run in Jaguar21
using the B3LYP functional22 in combination with the LACVP*
basis set.23 All gas-phase transition states were verified by
normal analysis obtained from a mass-weighted force-constant
matrix built up with analytical energy second derivatives. Each
transition state displayed exactly one mode with a negative
eigenvalue, as required. The corresponding eigenmodes were
inspected visually in Molekel24 and found to correspond to
Zimmerman-Traxler transition states, as indicated in Figure
3. As a further validation, a sample transition state structure of
the small model system was displaced forward and backward
along the negative eigenmode and subjected to energy mini-
mization, yielding the expected reactant and product.
Figure 3. Model systems employed in the computational study.
Calculations in solvent use the PB-SCRF method25 in Jaguar,
employing parameters suitable for THF (dielectric constant 7.43,
probe radius 2.52372 Å). PB-SCRF is a continuum solvation
model, where the molecule is put into a reaction field consisting
of surface charges on a solvent accessible surface constructed
using a hypothetical spherical solvent probe molecule with the
indicated radius.26 The wave function and the reaction field
charges are solved iteratively until self-consistency is reached.
In our experience, the major role of this electrostatic model is
to screen excessive charge separation. For other metal-assisted
reactions, we have shown that solvation models of this type
give agreement with experimental data in cases where gas-phase
calculations deviate strongly.27 In the current case, we found
that for the small model system, the neutral complex (including
the bromide counterion) gave very similar results in the gas
phase and in solvent. The six-membered rings of the Zimmer-
man-Traxler transition states were virtually superimposable.
Furthermore, the energy difference between the boat and chair
forms was found to be the same in the gas phase and in solvent
(16 kJ/mol in favor of the chair). The position of the bromide
did differ between the two calculations, and the loose coordina-
tion of the bromide was found to give severe convergence
problems in the solvent calculations. We also tested the cationic
model system (without the bromide counterion). Gratifyingly,
the structures obtained with the cationic model system in solvent
were very close to those obtained with the bromide ion present,
and the relative energy of the boat and chair forms only shifted
slightly, to 14 kJ/mol. Not unexpectedly, when the cationic
system was reconverged in the gas phase, the structures were
distorted, even though the relative energy of the boat and the
chair was still reasonable, 20 kJ/mol.
From previous projects,27 we have found that energies
calculated with the continuum solvent model gives a fair
correspondence with experimental ratios, but for a quantitative
agreement we frequently also need to account for the vibrational
component of the free energy. With our current resources, this
component must be determined from gas-phase calculations, and
the combination of this contribution with the solvated energies
is only valid if the geometries in the gas phase and solvent are
similar. Thus, in the current approach, we have determined all
transition states both in solvent and in the gas phase, verified
the similarity of the structures by superimposition, and per-
formed a vibrational analysis in the gas phase. The vibrational
analysis was used both to verify that the eigenmode correspond-
ing to the single negative eigenvalue indeed corresponds to the
We first studied the neutral and cationic forms of the small
model system (Figure 3), in the gas phase and in solvent.
(20) (a) Rasmussen, T.; Norrby, P.-O. J. Am. Chem. Soc. 2001, 123, 2464-
2465. (b) Rudolph, J.; Rasmussen, T.; Bolm, C.; Norrby, P.-O. Angew.
Chem., Int. Ed. 2003, 42, 3002-3005. (c) Rudolph, J.; Bolm, C.; Norrby,
P.-O. J. Am. Chem. Soc. 2005, 127, 1548-1552. (d) Cozzi, P. G.; Rudolph,
J.; Bolm, C.; Norrby, P.-O.; Tomasini, C. J. Org. Chem. 2005, 70, 5733-
5736.
(21) Jaguar 4.2; Schro¨dinger, Inc.: Portland, OR, 2000. For the most recent
(22) (a) Becke, A. D. J. Chem. Phys. 1993, 98, 5648-5652. (b) Lee, C.; Yang,
W.; Parr, R. G. Phys. ReV. B 1988, 37, 785-789. (c) Stephens, P. J.; Devlin,
F. J.; Chabalowski, C. F.; Frisch, M. J. J. Phys. Chem. 1994, 98, 8, 11623-
11627.
(25) Marten, B.; Kim, K.; Cortis, C.; Friesner, R. A.; Murphy, R. B.; Ringnalda,
M. N.; Sitkoff, D.; Honig, B. J. Phys. Chem. 1996, 100, 11775-11788.
(26) For an illuminating discussion about implicit solvation models in general,
see: Cramer, C. Essentials of Computational Chemistry: Theories and
Models; Wiley: New York, 2002.
(23) LACVP* uses the 6-31G* basis set for all light elements and the Hay-
Wadt ECP and basis set for Zn and Br; see: Hay, P. J.; Wadt, W. R. J.
Chem. Phys. 1985, 82, 270-283.
(24) Molekel, version 4.3, 2002, Stefan Portmann, CSCS/ETHZ; see http://
(27) (a) Hagelin, H.; A° kermark, B.; Norrby, P.-O. Chem. Eur. J. 1999, 5, 902-
909. (b) Norrby, P.-O.; Mader, M.; Vitale, M.; Prestat, G.; Poli, G.;
Organometallics 2003, 22, 1849-1855. (c) Kieken, E.; Wiest, O.; Helquist,
P.; Cucciolito, M. E.; Flores, G.; Vitagliano, A.; Norrby, P.-O. Organo-
metallics 2005, 24, 3737-3745.
9
J. AM. CHEM. SOC. VOL. 127, NO. 45, 2005 15759