Density Functional Study of the Relative
Reactivity in the Concerted 1,3-Dipolar
Cycloaddition of Nitrile Ylide to
Disubstituted Ethylenes
HSIN-YI LIAO,
1MING-DER SU,
2WEN-SHENG CHUNG,
3SAN-YAN CHU
11Department of Chemistry, National Tsing Hua University, Hsinchu 30043, Taiwan, ROC 2School of Chemistry, Kaohsiung Medical University, Kaohsiung 80708, Taiwan, ROC
3Department of Applied Chemistry, National Chiao Tung University, Hsinchu 30050, Taiwan, ROC
Received 27 November 2000; revised 21 February 2001; accepted 5 March 2001
ABSTRACT:
Density functional theory was used to perform a theoretical evaluation of (E)-1,2-disubstituted ethylenes as dipolarophiles for the 1,3-dipolar cycloaddition reaction. The reactivities of electron-withdrawing and -donating substituted ethylenes were examined by estimating their activation energies. The calculated activation energies predicted that the most reactive species is (E)-1,2-C2H2(NO)2, whereas the least reactive is(E)-2-butene. Namely, it was demonstrated that 16-electron 1,3-dipole reactants with more electropositive substituents in terminal positions and ethylenes that possess more strongly electron-withdrawing substituents facilitate 1,3-dipolar cycloaddition reactions. All of the theoretical results can be rationalized using the configuration mixing model.
c
2001 John Wiley & Sons, Inc. Int J Quantum Chem 83: 318–323, 2001
Key words:
1,3-dipolar cycloaddition reaction; configuration mixing modelIntroduction
T
he 1,3-dipolar cycloaddition reaction is a powerful tool for synthesizing five-membered heterocyclic compounds. Many 1,3-dipoles are iso-electronic to an allyl anion, which has four πelec-Correspondence to: M.-D. Su.
Contract grant sponsor: National Science Council of Taiwan.
trons and at least one charge-separated resonance structure that has opposite charges in a 1,3 re-lationship. It is this structural feature that leads to the name 1,3-dipolar cycloaddition reaction. This reaction is generally assumed to be a con-certed [π4s + π2s] cycloaddition [1, 2] like the
Diels–Alder reaction [3]. This assumption is sup-ported by their high regio- and stereoselectivity, although a stepwise radical mechanism has been proposed [4].
Perhaps the most outstanding contribution to this field has come from the brilliant work of Hoisgen [5]. His efforts have helped to transform 1,3-dipolar addition from an almost obscure phe-nomenon into a major reaction type. For compre-hensive reviews of 1,3-dipolar cycloaddition, see Ref. [6]. Although there are excellent reviews [5, 6] in this field, several basic questions need to be solved for the sake of synthesis applications. One important question that arouses our interest is how to increase the reactivity of the dipolarophile and/or dipole. Our primary aim in this work is, therefore, to propose a theoretical model that can fa-cilitate the 1,3-dipolar cycloaddition reaction.
It is generally acknowledged that substituents may influence reaction kinetics. In our previous study [7], we found that a 16-electron 1,3-dipole re-actant with more electropositive substituents at the terminal positions may facilitate its cycloaddition with a dipolarophile. Nitrile ylide is thus a good example molecule for the 1,3-dipole. In addition to choosing a highly reactive dipole, increasing the po-tential reactivity of the dipolarophile by a proper choice of substituents is also expected to increase the reactivity. Indeed, it is of theoretical and experimen-tal interest to determine the reactivity of different substituted dipolarophiles toward nitrile ylide. For this purpose, we choose ethylene, a classical dipo-larophile whose reactivity can be changed by either electron-withdrawing or -donating substituents.
In this work, we report a density functional study of the effects of substituents on the 1,3-dipolar cycloaddition reaction. The systems we chose to investigate here are the addition of ni-trile ylide to (E)-1,2-C2H2(NO)2, (E)-1,2-C2H2(CN)2,
(E)-1,2-C2H2Cl2, (E)-1,2-C2H2(OH)2, and
(E)-1,2-C2H2(CH3)2.
There are several ways to determine the reactiv-ity of different compounds in a reaction. One very popular and mostly qualitative method is to de-termine the energy gap between the frontier mole-cular orbitals (FMO) of the reactants. According to this frontier orbital theory [8], the rate of the reaction is determined largely by the degree of high-est occupied molecular orbital (HOMO) to lowhigh-est unoccupied molecular orbital (LUMO) reactant in-teractions. These energy effects were explained by Sustmann [9] and Houk et al. [10] using perturba-tion theory [11]. The smaller the energy gap between the controlling orbitals, the faster the reaction. The reaction is controlled by both electron-donating and -withdrawing substituents on either component [7]. Nevertheless, we believe that a somewhat
differ-ent approach and emphasis on other aspects of the reaction analyzed here may supplement this approach. We will show that the reactivity of the 1,2-disubstituted ethylene is correlated strongly to its singlet–triplet splitting, which is important for un-derstanding the origin of the substituent effects. We mention here that reports concerning the theoretical study of such substituent effects on the 1,3-dipolar cycloadditions are very rare [12], in marked con-trast to the Diels–Alder reactions [5]. In particular, to our knowledge, to date no theoretical work has been devoted to the study of the reactivity of 1,2-disubstituted ethylenes with a 1,3-dipole.
Methodology
Gas phase geometry optimizations without any restraints were carried out with Gaussian 94 [13] at the B3LYP level using the 6-31G∗ basis set [5] (denoted as B3LYP/6-31G∗). The reason for us-ing the B3LYP/6-31G∗ method is that our previ-ous study [7] indicated a good agreement between the B3LYP and CCSD(T) results for the 1,3-dipolar cycloaddition reactions. We have compared the B3LYP/6-31+G∗∗and B3LYP/6-31G∗results for the
systems. The changes in E‡and H are+0.128 and
−1.13 kcal/mol, respectively. All stationary points were identified by performing vibrational analysis with the same basis set used for the optimization. Additionally, we used these frequencies to evalu-ate the corresponding zero-point vibrational energy (ZPE) corrections to the relative energies.
Results and Discussion
The theoretical investigations [7, 16] suggest that the 1,3-dipolar cycloaddition reaction proceeds via initial formation of a precursor complex and sub-sequent rearrangement through a five-membered cyclic transition state to yield the eventual cycload-dition product. Thus, in this reaction the results for four regions on the potential energy surfaces will be presented: the reactants (R; nitrile ylide + 1,2-disubstituted ethylene), the orientation com-plex (C), the transition state (TS), and the cycload-dition product (P). The total energies and rela-tive energies, including the ZPE corrections for the B3LYP/6-31G∗level of theory used in the geometry determinations, are collected in Table I.
Basically, the kinetic stability of the initial orientation complex depends on the depth of this
TABLE I
Total energiesaand relative energiesbcomputed at the B3LYP/6-31G∗levels of theory.
Systems R C TS P C2H2(NO)2 −469.74258 −469.74508 −469.74485 −469.83570 (0.0) (−1.565) (−1.422) (−58.43) [0.0] [−3.078] [−2.810] C2H2(CN)2 −395.64641 −395.65128 −395.64410 −395.73741 (0.0) (−3.057) (+1.448) (−57.10) C2H2Cl2 −1130.37090 −1130.37175 −1130.35485 −1130.47251 (0.0) (−0.5209) (+10.09) (−63.75) C2H2(OH)2 −361.57418 −361.58212 −361.55320 −361.66825 (0.0) (−4.981) (+13.17) (−59.03) C2H2(CH3)2 −289.74496 −289.74293 −289.72214 −289.83862 (0.0) (+1.271) (+14.32) (−58.77)
aValues in atomic units.
bValues in parentheses are the relative energies (in kcal/mol) that corresponding to its reactants. Values in brackets are the relative
energies computed at the CCSD(T)/6-31G∗//B3LYP/6-31G∗levels of theory. local minimum, i.e., on the complexation energy
Ecas well as on the height of the energy barrier for
the transition state (E‡). As can be seen in Table I,
the orientation complexes for (E)-1,2-C2H2Cl2,
(E)-1,2-C2H2(OH)2, and (E)-1,2-C2H2(CH3)2 reside
in a fairly shallow minimum as shown by a complexation energy Ec of −0.5–1.3 kcal/mol
and an activation energy for the addition reaction of 10–14 kcal/mol. On the other hand, in the case of election-withdrawing substituted ethylene, the energy of the transition state relative to its corresponding orientation complex is 4.5 [(E)-1,2-C2H2(CN)2] and 0.14 [(E)-1,2-C2H2(NO)2] kcal/mol,
whereas the complexation energy of the orientation complex is −3.1 [(E)-1,2-C2H2(CN)2] and −1.6
[(E)-1,2-C2H2(NO)2] kcal/mol. It appears that the
potential energy surfaces for the (E)-1,2-C2H2(CN)2
and (E)-1,2-C2H2(NO)2 orientation complexes
are fairly flat. In addition, analytical force con-stant evaluations show that these complexes are indeed local minima. Accordingly, all of these theoretical results strongly indicate that the nitrile ylide forms a loose complex with the substituted ethylene and should not have a sufficient life-time to be a spectroscopically observable species. Nevertheless, orientation complexes, due to van der Waals interactions, have recently have been found for ozone and ethylene in the gas phase and were confirmed by ab initio calculations; see [17].
The geometries of the B3LYP/6-31G∗ transi-tion structures for each reactransi-tion are presented in Figure 1. The most interesting geometrical
parame-ters for these species along with the transition vec-tors are also shown in Figure 1. All transition struc-tures had only one imaginary vibrational frequency, corresponding to the motion for the newly forming C—C bonds in the concerted mechanism. In addi-tion, all transition states represent the concerted, but not necessarily synchronous, mechanism of the reaction. The degree of synchronicity of the reactions can be determined by comparing the bond distances for the two newly forming bonds. As expected, all the transition structures are more or less asynchronous because the reactants are asymmetric. Indeed, the asynchronicity strongly de-pends on the polarity of both the 1,3-dipole and the dipolarophile, as well as the electronic and steric interaction between the constituents in the transition structures. The degree of asynchronic-ity for these transition structures is 0.309 Å [TS-C2H2(NO)2], 0.198 Å [TS-C2H2(CN)2], 0.107 Å
[TS-C2H2Cl2], 0.154 Å [TS-C2H2(OH)2], and 0.025 Å
[TS-C2H2(CH3)2]. Thus, by comparing the
differ-ences in the bond lengths for these cycloaddition reactions, the highest degree of asychronicity is pre-dicted for the transition structure TS-C2H2(NO)2,
whereas the transition structure TS-C2H2(CH3)2
is the most synchronous. In other words, from these theoretical results, it seems possible to conclude that the more electron-withdrawing substituted re-actants should form products through more asyn-chronous transition states, and the more electron-donating substituted reactants will go through more synchronous transition structures as is the case with (E)-1,2-C2H2(CH3)2.
FIGURE 1. B3LYP/6-31G∗optimized geometries (in angstroms) of the transition state of (E)-1,2-C2H2(NO)2,
(E)-1,2-C2H2(CN)2, (E)-1,2-C2H2Cl2, (E)-1,2-C2H2(OH)2,
and (E)-1,2-C2H2(CH3)2. The heavy arrows indicate
the eigenvectors for the main atomic motion in the transition state. The corresponding imaginary frequencies are given.
We can understand the unsymmetrical structures of the transition state by the following consider-ation. We used the quantity R = RCC(right) −
RCC(left) (see Fig. 1) to indicate the difference
be-tween the two newly forming C—C bonds. We found that R for the transition state R(TS) corre-late well with those for the product R(P), as is well expected for cycloaddition reaction with late tran-sition states. For the products, the values of R(P) in decreasing order are NO (0.047) > CH3(0.042) >
CN (0.034) > Cl (0.032) > OH (0.026), which agree qualitatively with R(TS) NO (0.31) > CN (0.20) >
Cl (0.11) > CH3 (0.025) > OH (−0.15) with only
CH3misplaced. We suspect that CH3has an
impor-tant steric effect, but it is short ranged and therefore more sensitive for the product than the transition state.
Furthermore, considering the substituent effect, the density-functional theory (DFT) calculations suggest that the more strongly electron-with-drawing the substituent, the lower the activation energy of the cycloaddition reaction. For instance, because electron-withdrawing ability is ordered NO > CN > Cl ∼ OH > CH3 [18], the barrier
height for the substituted ethylene addition to nitrile ylide increases in the order TS-C2H2(NO)2
(−1.4 kcal/mol) < TS-C2H2(CN)2(1.4 kcal/mol) <
TS-C2H2Cl2 (10 kcal/mol) < TS-C2H2(OH)2
(13 kcal/mol) < TS-C2H2(CH3)2 (14 kcal/mol) at
the B3LYP/6-31G∗ level of theory (see Table I). That is to say, from a kinetic viewpoint, the most reactive dipolarophile is (E)-1,2-C2H2(NO)2,
followed by (E)-1,2-C2H2(CN)2, (E)-1,2-C2H2Cl2,
(E)-1,2-C2H2(OH)2, and (E)-1,2-C2H2(CH3)2.
As mentioned earlier, according to FMO the-ory [8], higher reactivity is predicted for the pair of reactants that have similar HOMO–LUMO en-ergies or small energy gaps between the HOMO of one reactant, and the LUMO of the other [9, 10]. The differences in frontier orbital energies between the reactants are presented in Table II. We have thus examined the relationship between the Hartree–Fock orbital HOMO–LUMO energy gaps and the activation barriers for the afore-mentioned five systems as given in Table II. It is apparent that neither HOMOdipole–LUMOethylene
nor HOMOethylene–LUMOdipole energy gaps
corre-late with the addition barriers. (Compounds with high-lying HOMOs are nucleophilic and those with low-energy LUMOs are electrophilic. Accord-ing to FMO theory, highly nucleophilic 1,3-dipoles and potent electrophilic dipolarophiles are needed.) However, the sum of the two, which leads to a HOMO–LUMO energy gap for various ethylenes, has a better correlation with the barrier. [From Table II, the HOMO–LUMO energy gaps based on the Hartree–Fock orbitals are 0.04372, 0.4555, 0.5100, 0.5469, and 0.5418 au for trans-1,2-C2H2(NO)2,
trans-1,2-C2H2(CN)2, trans-1,2-C2H2Cl2,
trans-1,2-C2H2(OH)2, and trans-1,2-C2H2(CH3)2.] This is not
surprising because FMO theory considers only the separated reactants without incorporating the pos-sible extra interactions of the reactants in the tran-sition state, which can play a crucial role in de-termining the reactivity. An alternative model for
TABLE II
HOMO and LUMO energies (au), singlet–triplet energy gaps (kcal/mol), barriers (kcal/mol), and enthalpies (kcal/mol) for 1,3-dipolar cycloadditions.a
System HOMO LUMO EHO–LUb ELU–HOc Estd E‡e Hf
C2H2(NO)2 −0.42298 +0.01422 0.3303 0.5727 +14.35 −1.4216 −58.43
C2H2(CN)2 −0.41716 +0.03834 0.3544 0.5669 +71.47 1.4480 −57.10
C2H2Cl2 −0.37059 +0.13940 0.4555 0.5204 +91.92 +10.09 −63.75
C2H2(OH)2 −0.33247 +0.21443 0.5305 0.4822 +97.65 +13.17 −59.03
C2H2(CH3)2 −0.34130 +0.20053 0.5166 0.4911 +102.21 +14.32 −58.77
aAll at the B3LYP/6-31G∗level. The HOMO and LUMO are based on the Hartree–Fock orbitals.
bThe energy difference between the HOMO of the 1,3-dipoles and the LUMO of the ethylene. The HOMO of reactant H
2CNCH is
−0.3161 au.
cThe energy difference between the HOMO of ethylene and the LUMO of the 1,3-dipoles. The LUMO of reactant H
2CNCH is
+0.1498 au.
dA positive value indicates a singlet ground state.
eThe activation energy of the transition state, relative to its corresponding reactants. fThe exothermicity of the products, relative to its corresponding reactants (also see Table I).
interpreting the reactivity of cycloaddition reactions is provided by the so-called configuration mixing (CM) model. The details of the CM model have been discussed in the literature [7, 19], and need not be reiterated here. According to this model, the energy barriers that govern processes with similar reaction energies should be proportional to the en-ergy gap Est(= Etriplet−Esinglet) between the singlet
ground state and the first excited triplet state of the disubstituted ethylene. As can be seen in Table II, the trend for the singlet–triplet splitting Estof the
disubstituted ethylene is in general agreement with the computed energy barriers, i.e., the smaller the Est of the disubstituted ethylene, the lower the
barrier height. As a result, our theoretical findings suggest that if a disubstituted ethylene has a singlet ground state with a low-lying triplet state (lead-ing to a smaller Est, also a smaller HOMO–LUMO
gap), it may readily undergo a concerted reaction [7, 19]. The HOMOethylene–LUMOethylenegap, which
is closely related to Est, has a good correlation with
the activation barriers.
Potential energy profiles based on the B3LYP data in Table I are summarized in Figure 2. Moreover, as discussed previously [7], our DFT results indicate that the stronger electron-withdrawing substituents should give a lower barrier for the 1,3-dipolar cy-cloaddition reaction. The reason for this can be simply understood in terms of the singlet–triplet splitting (Est). It is well known that a good electron
withdrawer (or a good π acceptor, such as NO) will stabilize the π∗level of ethylene much more than the
π level of ethylene. This will decrease the HOMO– LUMO gap in the substituted ethylene [19, 20] and, in turn, lower its singlet–triplet splitting. Such an effect will facilitate the 1,3-dipolar cycloaddition re-action, as mentioned earlier.
FIGURE 2. Potential energy surfaces for the cycloadditions of nitrile ylide with (E)-1,2-disubstituted ethylenes. The relative energies are taken from the B3LYP/6-31G∗values as given in Table I. For optimized structures of the transition states, see Figure 1.
In summary, in the present work we have stud-ied disubstituted ethylenes in reaction with ni-trile ylide to determine the influence of electron-withdrawing and -donating substituents on the reactivity. Given the importance of substitution on both the 1,3-dipole [7] and the dipolarophile, it is, therefore, predicted that a 16-electron 1,3-dipole reactant with more electropositive substituents in terminal positions and an ethylene bearing more strongly electron-withdrawing substituents will fa-cilitate 1,3-dipolar cycloaddition reactions.
ACKNOWLEDGMENTS
We are grateful to the National Center for High-Performance Computing of Taiwan and the Com-puting Center at Tsing Hua University for generous amounts of computing time. We also thank the Na-tional Science Council of Taiwan for their financial support. We express our gratitude to the referees for their useful comments.
References
1. Houk, K. N.; Tamaguchi, K. In 1,3-Dipolar Cycloaddition Chemistry; Padwa, A., Ed.; Wiley: New York, 1987; Vol. 2. 2. Suser, J.; Sustmann, R. Angew Chem Int Ed Engl 1980, 19,
779. Houk, K. N.; Li, Y.; Evanseck, J. D. Angew Chem Int Ed Engl 1992, 31, 682. Dewar, M. J. S.; Caoxian, J. Acc Chem Res 1992, 25, 537.
3. Firestone, R. A. J Org Chem 1968, 33, 2285; 1972, 37, 2181; J Am Chem Soc 1970, 1570. Huisgen, R. J Org Chem 1976, 41, 403.
4. Huisgen, R. J.; Weinberger, R. Tetrahedron Lett 1985, 26, 5119. Houk, K. N.; Firestone, R. A.; Munchausen, L. L.; Mueller, P. H.; Arison, B. H.; Garcia, L. A. J Am Chem Soc 1985, 107, 7227. McDouall, J. J. W.; Robb, M. A.; Niazi, U.; Bernardi, F.; Schlegel, H. B. J Am Chem Soc 1987, 109, 4642. Quast, H.; Regnat, D.; Peters, E.-M.; Peters, K.; von Schnering, H. G. Angew Chem Int Ed Engl 1990, 29, 695. Sustmann, R.; Sicking, W.; Quast, H. J. J Comput Chem 1992, 13, 314.
5. Huisgen, R. Angew Chem Int Ed Engl 1963, 2, 562; 1963, 2, 633.
6. Bianch, G.; DeMicheli, C.; Gandolfi, R. In The Chemistry of Double Bonded Functional Groups; Patai, S., Ed.; Wi-ley: New York, 1977; Part I, Suppl. A, p. 369. Padwa, A. Angew Chem Int Ed Engl 1976, 15, 123. Padwa, A., Ed. In
1,3-Dipolar Cycloaddition Chemistry; Wiley: New York, 1984; Vol. 1. Coppola, B. P.; Noe, M. C.; Schwartz, D. J.; Abdon, R. L., II; Trost, B. M. Tetrahedron 1994, 50, 93. Huis-gen, R. In Advances in Cycloaddition; Curan, D. P., Ed.; JAI Press: Greenwich, CT, 1988; Vol. 1, p. 1.
7. Su, M.-D.; Liao, H.-Y.; Chung, W.-S.; Chu, S.-Y. J Org Chem 1999, 64, 6710, and references therein.
8. Fukui, K. In Molecular Orbitals in Chemistry, Physics and Biology; Lowdin, P. O.; Pullman, B., Eds.; Academic Press: New York, 1964. Fleming, I. In Frontier Molecular Orbitals and Organic Chemical Reaction; Wiley: New York, 1976. Fukui, K. Acc Chem Res 1971, 4, 57. Woodward, R. B.; Hoff-mann, R. Angew Chem Int Ed Engl 1969, 8, 781.
9. Sustmann, R. Tetrahedron Lett 1971, 00, 2721.
10. Houk, K. N.; Sims, J.; Duke, R. E., Jr.; Strozier, R. W.; George, J. K. J Am Chem Soc 1973, 95, 7287. Houk, K. N.; Sims, J.; Watts, C. R.; Luskus, T. J. J Am Chem Soc 1973, 95, 7301. 11. Albright, T. A.; Burdett, J. K.; Whangbo, M.-H. In Orbital
Interactions in Chemistry; Wiley: New York, 1985. 12. Klicic, J. J.; Friesner, R. A. J Phys Chem 1999, 103, 1276. 13. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Gill, P. M. W.;
Johnson, B. G.; Robb, M. A.; Cheeseman, J. R.; Keith, T.; Petersson, G. A.; Montgomery, J. A.; Raghavachari, K.; Al-Laham, M. A.; Zakrzewski, V. G.; Ortiz, J. V.; Foresman, J. B.; Cioslowski, J.; Stefanov, B. B.; Nanayakkara, A.; Challa-combe, M.; Peng, C. Y.; Ayala, P. Y.; Chen, W.; Wong, M. W.; Andres, J. L.; Replogle, E. S.; Gomperts, R.; Martin, R. L.; Fox, D. J.; Binkley, J. S.; Defrees, D. J.; Baker, J.; Stewart, J. P.; Head-Gordon, M.; Gonzalez, C.; Pople, J. A. Gaussian 94; Gaussian, Inc.: Pittsburgh, PA, 1995.
14. Becke, A. D. J Chem Phys 1992, 97, 9173; 1993, 98, 5648. Lee, C.; Yang, W.; Parr, R. G. Phys Rev B 1988, 37, 785. 15. Hehre, W. J.; Ditchfield, R.; Pople, J. A. J Chem Phys 1972,
56, 2257.
16. Huisgen, R.; Fisera, L.; Giera, H.; Sustmann, R. J Am Chem Soc 1995, 117, 9671. Sustmann, R.; Sicking, W.; Huisgen, R. J Am Chem Soc 1995, 117, 9679.
17. Gillies, C. W.; Gillies, J. Z.; Suenram, R. D.; Lovas, F. J.; Kraka, E.; Cremer, D. J Am Chem Soc 1991, 113, 2412. Gillies, J. Z.; Gillies, C. W.; Lovas, F. J.; Matsumura, K.; Suenram, R. D.; Kraka, E.; Cremer, D. J Am Chem Soc 1991, 113, 6408. 18. March, J. Advanced Organic Chemistry; Wiley: New York,
1992; pp. 509, 652.
19. Shaik, S.; Schlegel, H. B.; Wolfe, S. In Theoretical Aspects of Physical Organic Chemistry; Wiley: New York, 1992. Pross, A. In Theoretical and Physical Principle of Organic Reactivity; Wiley: New York, 1995. Su, M.-D. Inorg Chem 1995, 34, 3829. Shaik, S.; Shurki, A. Angew Chem Int Ed Engl 1999, 38, 586.
20. Wade, P. A.; Hinney, H. R. J Am Chem Soc 1979, 101, 1319. Brittelli, D. R.; Boswell, G. A. J Org Chem 1981, 46, 312, 316.