Archive for the 'Authors' Category

Mechanism of dimethyldioxirane oxidation

Dimethyldioxirane can oxidize alkanes to alcohols. The mechanism for the oxidation has been controversial, ranging from concerted, to radical intermediates to an H-abstraction—O-rebound mechanism. Yang, Yu, and Houk now offer a molecular dynamics examination of the reaction of dimethyldioxirane with isobutane.1

Gas–phase (U)B3LYP/6-311++G(d,p)//(U)B3LYP/6-31G(d) computations give critical points outlined in the reaction below. The structures of the transition states and the intermediate are shown in Figure 1.

TS1
26.6
(25.9)

INT
8.6
(9.3)

TS2
8.6
(-0.8)

Figure 1. (U)B3LYP/6-31G(d) optimized geometries of TS1, INT, and TS2. Relative free energies (kcal mol-1) in the gas (top) and solution (bottom) phases

The PES indicates a rebound mechanism, though in acetone solution phase, there was no transition state located for the second step; it appears to be barrierless. It should be noted that the size of the barrier is very small even in the gas phase. The energy given in Figure 1 is for the gas phase structure computed in solution.

Trajectories for both gas and solution phase were computed. For the gas phase, about 90% of the trajectories lead to separated radicals, but in an acetone about 90% of the trajectories lead directly to the alcohol, with only 10% leading to radicals. Even so, the acetone trajectories divide into two types, a dynamically concerted path where the time gap between the formation of the new C-O and O-H bonds is less than 60 fs, and a dynamically stepwise path where the time gap is greater than 60 fs, though for the trajectories that lead to product the gap is typically still less than 150 fs.

References

(1)  Yang, Z.; Yu, P.; Houk, K. N. "Molecular Dynamics of Dimethyldioxirane C–H Oxidation," J. Am. Chem. Soc. 2016, 138, 4237-4242, DOI: 10.1021/jacs.6b01028.

InChIs

Dimethyldioxirane: InChI=1S/C3H6O2/c1-3(2)4-5-3/h1-2H3
InChIKey=FFHWGQQFANVOHV-UHFFFAOYSA-N

Isobutane: InChI=1S/C4H10/c1-4(2)3/h4H,1-3H3
InChIKey=NNPPMTNAJDCUHE-UHFFFAOYSA-N

Houk Steven Bachrach 06 Jun 2016 No Comments

Benchmarking Platonic solids and related hydrocarbons

Karton, Schreiner, and Martin have benchmarked the heats of formation of some Platonic Solids and related hydrocarbons.1 The molecules examined are tetrahedrane 1, cubane 2, dodecahedrane 3, trisprismane 4, pentaprismane 5, and octahedrane 6.

The optimized structures (B3LYP-D3BJ/def2-TZVPP) of these compounds are shown in Figure 1.

1

2

3

4

5

6

Figure 1. B3LYP-D3BJ/def2-TZVPP optimized geometries of 1-6.

Using the W1-F12 and W2-F12 composite methods, the estimated the heats of formation of these hydrocarbons are listed in Table 1. Experimental values are available only for 2 and 3; the computed values are off by about 2 kcal mol-1, which the authors argue is just outside the error bars of the computations. They suggest that the experiments might need to be revisited.

Table 1. Heats of formation (kcal mol-1) of 1-6.

compound

ΔHf (comp)

ΔHf (expt)

1

128.57

 

2

144.78

142.7 ± 1.2

3

20.18

22.4 ± 1

4

132.84

 

5

119.42

 

6

63.08

 

They conclude with a comparison of strain energies computed using isogyric, isodesmic, and homodesmotic reactions with a variety of computational methods. Somewhat disappointingly, most DFT methods have appreciable errors compared with the W1-F12 results, and the errors vary depend on the chemical reaction employed. However, the double hybrid method DSD-PBEP86-D3BJ consistently reproduces the W1-F12 results.

References

(1)  Karton, A.; Schreiner, P. R.; Martin, J. M. L. "Heats of formation of platonic hydrocarbon cages by means of high-level thermochemical procedures," J. Comput. Chem. 2016, 37, 49-58, DOI: 10.1002/jcc.23963.

InChIs

1: InChI=1S/C4H4/c1-2-3(1)4(1)2/h1-4H
InChIKey=FJGIHZCEZAZPSP-UHFFFAOYSA-N

2: InChI=1S/C8H8/c1-2-5-3(1)7-4(1)6(2)8(5)7/h1-8H
InChIKey=TXWRERCHRDBNLG-UHFFFAOYSA-N

3: InChI=1S/C20H20/c1-2-5-7-3(1)9-10-4(1)8-6(2)12-11(5)17-13(7)15(9)19-16(10)14(8)18(12)20(17)19/h1-20H
InChIKey=OOHPORRAEMMMCX-UHFFFAOYSA-N

4: InChI=1S/C6H6/c1-2-3(1)6-4(1)5(2)6/h1-6H
InChIKey=RCJOMOPNGOSMJU-UHFFFAOYSA-N

5: InChI=1S/C10H10/c1-2-5-3(1)7-8-4(1)6(2)10(8)9(5)7/h1-10H
InChIKey=TXSUZIHPJIZTJT-UHFFFAOYSA-N

6: InChI=1S/C12H12/c1-2-4-6-5(11-7(1)10(4)11)3(1)9-8(2)12(6)9/h1-12H
InChIKey=DUDCBTWUEJEMGD-UHFFFAOYSA-N

QM Method &Schreiner Steven Bachrach 10 May 2016 No Comments

Diels-Alder reactions of some arenes

Houk has examined the Diels-Alder reaction involving ethene with benzene 1 and all of its aza-substituted isomers having four or fewer nitrogen atoms 2-11.1 The reactions were computed at M06-2X/6-311+G(d,p).

All of the possible Diels-Alder reactions were examined, and they can be classified in terms of whether two new C-C bonds are formed, one new C-C and one new C-N bond are formed, or two new C-N bonds are formed. Representative transition states of these three reaction types are shown in Figure 1, using the reaction of 7 with ethene.

Figure 1. M06-2X/6-311+G(d,p) optimized transition states for the Diels-Alders reactions of 7 with ethene.

A number of interesting trends are revealed. For a given type of reaction (as defined above), as more nitrogens are introduced into the ring, the activation energy decreases. Forming two C-C bonds has a lower barrier than forming a C-C and a C-N, which has a lower barrier than forming two C-N bonds. The activation barriers are linearly related to the aromaticity of the ring defined by either NICS(0) or aromatic stabilization energy, with the barrier decreasing with decreasing aromaticity. The barrier is also linearly related to the exothermicity of the reaction.

The activation barrier is also linearly related to the distortion energy. With increasing nitrogen substitution, the ring becomes less aromatic, and therefore more readily distorted from planarity to adopt the transition state structure.

References

(1) Yang, Y.-F.; Liang, Y.; Liu, F.; Houk, K. N. "Diels–Alder Reactivities of Benzene, Pyridine, and Di-, Tri-, and Tetrazines: The Roles of Geometrical Distortions and Orbital Interactions," J. Am. Chem. Soc. 2016, 138, 1660-1667, DOI: 10.1021/jacs.5b12054.

Aromaticity &Diels-Alder &Houk Steven Bachrach 26 Apr 2016 No Comments

Cyclization reaction of 1,2-cyclohexadiene

1,2-Cyclohexadiene 1 is a very strained and highly reactive species. Houk, Garg and co-workers report on its use as the ene component in a cyclization with a 1,3-dipole, namely nitrones.1 For example, 1 reacts with nitrone 2 to give the cycloadducts 3a and 3b in a ratio of 8.9:1.

To investigate the mechanism of this reaction, they optimized the structures of all compounds at CPCM(acetonitrile)B3LYP/6-31G(d) and single-point energies were obtained using the B3LYP-D3 functional. The structures of some pertinent critical points are shown in Figure 1. They did locate a concerted transition state (TS1) leading to 3a, with a barrier of 14.5 kcal mol-1, but could not find a concerted TS leading to 3b. (Also, the barriers leading to the other regioisomer are much higher than the ones leading to the observed products.) Rather, they identified a stepwise transition state (TS2) with a barrier of nearly the same energy (14.4 kcal mol-1) that leads to the intermediate (INT), which lies 16.5 kcal mol-1 below reactants. They located two transition states from his intermediate, TS3a and TS3b, leading to the two different products. The barrier to 3a is 1.2 kcal mol-1 lower than the barrier leading to 3b, and this corresponds nicely with the observed diastereoselectivity.

1
(0.0)

TS1
(14.5)

TS2
(14.4)

 

INT
(-16.5)

 

TS3a
(-7.4)

TS3b
(-6.2)

Figure 1. CPCM(acetonitrile)B3LYP/6-31G(d) optimized geometries and CPCM(acetonitrile)B3LYP-D3/6-31G(d) free energies.

References

(1) Barber, J. S.; Styduhar, E. D.; Pham, H. V.; McMahon, T. C.; Houk, K. N.; Garg, N. K.
"Nitrone Cycloadditions of 1,2-Cyclohexadiene," J. Am. Chem. Soc. 2016, 138, 2512-2515, DOI: 10.1021/jacs.5b13304.

InChIs

1: InChI=1S/C6H8/c1-2-4-6-5-3-1/h1,5H,2,4,6H2
InChIKey=NMGSDTSOSIPXTN-UHFFFAOYSA-N

2: InChI=1S/C11H15NO/c1-11(2,3)12(13)9-10-7-5-4-6-8-10/h4-9H,1-3H3/b12-9-
InChIKey=IYSYLWYGCWTJSG-XFXZXTDPSA-N

3a: InChI=1S/C17H23NO/c1-17(2,3)18-16(13-9-5-4-6-10-13)14-11-7-8-12-15(14)19-18/h4-6,9-11,15-16H,7-8,12H2,1-3H3/t15-,16-/m0/s1
InChIKey=HQIBSHJEOJPFTI-HOTGVXAUSA-N

3b: InChI=1S/C17H23NO/c1-17(2,3)18-16(13-9-5-4-6-10-13)14-11-7-8-12-15(14)19-18/h4-6,9-11,15-16H,7-8,12H2,1-3H3/t15-,16+/m1/s1
InChIkey=HQIBSHJEOJPFTI-CVEARBPZSA-N

cycloadditions &Houk Steven Bachrach 11 Apr 2016 No Comments

Mechanism of organocatalysis by Cinchona alkaloids

Cinchona alkaloids cat catalyze reactions, such as shown in Reaction 1. Wynberg1 proposed a model to explain the reaction, shown in Scheme 1, based on NMR. Grayson and Houk have now used DFT computations to show that the mechanism actually reverses the arrangements of the substrates.2

Reaction 1

Scheme 1.


Wynberg Model


Grayson and Houk Model

M06-2X/def2-TZVPP−IEFPCM(benzene)//M06-2X/6-31G(d)−IEFPCM(benzene) computations show that the precomplex of catalyst 3 with nucleophile 1 and Michael acceptor 2 is consistent with Wynberg’s model. The alternate precomplex is 5.6 kcal mol-1 higher in energy. These precomplexes are shown in Figure 1.

Wynberg precomplex

Grayson/Houk precomplex

Figure 1. Precomplexes structures

However, the lowest energy transition state takes the Grayson/Houk pathway and leads to the major isomer observed in the reaction. The Grayson/Houk TS that leads to the minor product has a barrier that is 3 kcal mol-1 higher in energy. The lowest energy TS following the Wynberg path leads to the minor product, and is 2.2 kcal mol-1 higher than the Grayson/Houk path. These transition states are shown in Figure 2. The upshot is that complex formation is not necessarily indicative of the transition state structure.

Wynberg TS (major)
Rel ΔG = 5.3

Wynberg TS (minor)
Rel ΔG = 2.2

Grayson/Houk TS (major)
Rel ΔG = 0.0

Grayson/Houk TS (minor)
Rel ΔG = 3.0

Figure 2. TS structures and relative free energies (kcal mol-1).

References

(1) Hiemstra, H.; Wynberg, H. "Addition of aromatic thiols to conjugated cycloalkenones, catalyzed by chiral .beta.-hydroxy amines. A mechanistic study of homogeneous catalytic asymmetric synthesis," J. Am. Chem. Soc. 1981, 103, 417-430, DOI: 10.1021/ja00392a029.

(2) Grayson, M. N.; Houk, K. N. "Cinchona Alkaloid-Catalyzed Asymmetric Conjugate Additions: The Bifunctional Brønsted Acid–Hydrogen Bonding Model," J. Am. Chem. Soc. 2016, 138, 1170-1173, DOI: 10.1021/jacs.5b13275.

InChIs

1: InChI=1S/C10H14S/c1-10(2,3)8-4-6-9(11)7-5-8/h4-7,11H,1-3H3
InChIKey=GNXBFFHXJDZGEK-UHFFFAOYSA-N

2: InChI=1S/C8H12O/c1-8(2)5-3-4-7(9)6-8/h3-4H,5-6H2,1-2H3
InChIKey=CDDGRARTNILYAB-UHFFFAOYSA-N

3: InChI=1S/C18H22N2O/c1-12-11-20-9-7-13(12)10-17(20)18(21)15-6-8-19-16-5-3-2-4-14(15)16/h2-6,8,12-13,17-18,21H,7,9-11H2,1H3/t12?,13?,17?,18-/m1/s1
InChIKey=ZOZLJWFJLBUKKL-NKHWWFDVSA-N

4: InChI=1S/C18H26OS/c1-17(2,3)13-6-8-15(9-7-13)20-16-10-14(19)11-18(4,5)12-16/h6-9,16H,10-12H2,1-5H3/t16-/m0/s1
InChIKey=XUTYYZOSKLYWLW-INIZCTEOSA-N

Houk &Michael addition &stereoinduction Steven Bachrach 03 Mar 2016 No Comments

QM/MM trajectory of an aqueous Diels-Alder reaction

I discuss the aqueous Diels-Alder reaction in Chapter 7.1 of my book. A key case is the reaction of methyl vinyl ketone with cyclopentadiene, Reaction 1. The reaction is accelerated by a factor of 740 in water over the rate in isooctane.1 Jorgensen argues that this acceleration is due to stronger hydrogen bonding to the ketone than in the transition state than in the reactants.2-4

Rxn 1

Doubleday and Houk5 report a procedure for calculating trajectories including explicit water as the solvent and apply it to Reaction 1. Their process is as follows:

  1. Compute the endo TS at M06-2X/6-31G(d) with a continuum solvent.
  2. Equilibrate water for 200ps, defined by the TIP3P model, in a periodic box, with the transition state frozen.
  3. Continue the equilibration as in Step 2, and save the coordinates of the water molecules after every addition 5 ps, for a total of typically 25 steps.
  4. For each of these solvent configurations, perform an ONIOM computation, keeping the waters fixed and finding a new optimum TS. Call these solvent-perturbed transition states (SPTS).
  5. Generate about 10 initial conditions using quasiclassical TS mode sampling for each SPTS.
  6. Now for each the initial conditions for each of these SPTSs, run the trajectories in the forward and backward directions, typically about 10 of them, using ONIOM to compute energies and gradients.
  7. A few SPTS are also selected and water molecules that are either directly hydrogen bonded to the ketone, or one neighbor away are also included in the QM portion of the ONIOM, and trajectories computed for these select sets.

The trajectory computations confirm the role of hydrogen bonding in stabilizing the TS preferentially over the reactants. Additionally, the trajectories show an increasing asynchronous reactions as the number of explicit water molecules are included in the QM part of the calculation. Despite an increasing time gap between the formation of the first and second C-C bonds, the overwhelming majority of the trajectories indicate a concerted reaction.

References

(1) Breslow, R.; Guo, T. "Diels-Alder reactions in nonaqueous polar solvents. Kinetic
effects of chaotropic and antichaotropic agents and of β-cyclodextrin," J. Am. Chem. Soc. 1988, 110, 5613-5617, DOI: 10.1021/ja00225a003.

(2) Blake, J. F.; Lim, D.; Jorgensen, W. L. "Enhanced Hydrogen Bonding of Water to Diels-Alder Transition States. Ab Initio Evidence," J. Org. Chem. 1994, 59, 803-805, DOI: 10.1021/jo00083a021.

(3) Chandrasekhar, J.; Shariffskul, S.; Jorgensen, W. L. "QM/MM Simulations for Diels-Alder
Reactions in Water: Contribution of Enhanced Hydrogen Bonding at the Transition State to the Solvent Effect," J. Phys. Chem. B 2002, 106, 8078-8085, DOI: 10.1021/jp020326p.

(4) Acevedo, O.; Jorgensen, W. L. "Understanding Rate Accelerations for Diels−Alder Reactions in Solution Using Enhanced QM/MM Methodology," J. Chem. Theor. Comput. 2007, 3, 1412-1419, DOI: 10.1021/ct700078b.

(5) Yang, Z.; Doubleday, C.; Houk, K. N. "QM/MM Protocol for Direct Molecular Dynamics of Chemical Reactions in Solution: The Water-Accelerated Diels–Alder Reaction," J. Chem. Theor. Comput. 2015, , 5606-5612, DOI: 10.1021/acs.jctc.5b01029.

Diels-Alder &Houk &Solvation Steven Bachrach 02 Feb 2016 1 Comment

Fractional occupancy density

Assessing when a molecular system might be subject to sizable static (non-dynamic) electron correlation, necessitating a multi-reference quantum mechanical treatment, is perhaps more art than science. In general one suspects that static correlation will be important when the frontier MO energy gap is small, but is there a way to get more guidance?

Grimme reports the use of fractional occupancy density (FOD) as a visualization tool to identify regions within molecules that demonstrate significant static electron correlation.1 The method is based on the use of finite temperature DFT.2,3

The resulting plots of the FOD for a series of test cases follow our notions of static correlation. Molecules, such as alkanes, simple aromatics, and concerted transition states show essentially no fractional orbital density. On the other hand, the FOD plot for ozone shows significant density spread over the entire molecule; the transition state for the cleavage of the terminal C-C bond in octane shows FOD at C1 and C2 but not elsewhere; p-benzyne shows significant FOD at the two radical carbons, while the FOD is much smaller in m-benzyne and is negligible in o-benzyne.

This FOD method looks to be a simple tool for evaluating static correlation and is worth further testing.

References

(1) Grimme, S.; Hansen, A. "A Practicable Real-Space Measure and Visualization of Static Electron-Correlation Effects," Angew. Chem. Int. Ed. 2015, 54, 12308-12313, DOI: 10.1002/anie.201501887.

(2) Mermin, N. D. "Thermal Properties of the Inhomogeneous Electron Gas," Phys. Rev. 1965, 137, A1441-A1443, DOI: 10.1103/PhysRev.137.A1441.

(3) Chai, J.-D. "Density functional theory with fractional orbital occupations," J. Chem. Phys 2012, 136, 154104, DOI: doi: 10.1063/1.3703894.

Grimme Steven Bachrach 20 Jan 2016 2 Comments

Dispersion in organic chemistry – a review and another example

The role of dispersion in organic chemistry has been slowly recognized as being quite critical in a variety of systems. I have blogged on this subject many times, discussing new methods for properly treating dispersion within quantum computations along with a variety of molecular systems where dispersion plays a critical role. Schreiner1 has recently published a very nice review of molecular systems where dispersion is a key component towards understanding structure and/or properties.

In a similar vein, Wegner and coworkers have examined the Z to E transition of azobenzene systems (1a-g2a-g) using both experiment and computation.2 They excited the azobenzenes to the Z conformation and then monitored the rate for conversion to the E conformation. In addition they optimized the geometries of the two conformers and the transition state for their interconversion at both B3LYP/6-311G(d,p) and B3LYP-D3/6-311G(d,p). The optimized structure of the t-butyl-substituted system is shown in Figure 1.


a: R=H; b: R=tBu; c: R=Me; d: R=iPr; e: R=Cyclohexyl; f: R=Adamantyl; g: R=Ph

1b

1b-TS-2b

2b

Figure 1. B3LYP-D3/6-311G(d,p) optimized geometries of 1a, 2a, and the TS connecting them.

The experiment finds that the largest activation barriers are for the adamantly 1f and t-butyl 1b azobenzenes, while the lowest barriers are for the parent 1a and methylated 1c azobenzenes.

The trends in these barriers are not reproduced at B3LYP but are reproduced at B3LYP-D3. This suggests that dispersion is playing a role. In the Z conformations, the two phenyl groups are close together, and if appropriately substituted with bulky substituents, contrary to what might be traditionally thought, the steric bulk does not destabilize the Z form but actually serves to increase the dispersion stabilization between these groups. This leads to a higher barrier for conversion from the Z conformer to the E conformer with increasing steric bulk.

References

(1) Wagner, J. P.; Schreiner, P. R. "London Dispersion in Molecular Chemistry—Reconsidering Steric Effects," Angew. Chem. Int. Ed. 2015, 54, 12274-12296, DOI: 10.1002/anie.201503476.

(2) Schweighauser, L.; Strauss, M. A.; Bellotto, S.; Wegner, H. A. "Attraction or Repulsion? London Dispersion Forces Control Azobenzene Switches," Angew. Chem. Int. Ed. 2015, 54, 13436-13439, DOI: 10.1002/anie.201506126.

InChIs

1b: InChI=1S/C28H42N2/c1-25(2,3)19-13-20(26(4,5)6)16-23(15-19)29-30-24-17-21(27(7,8)9)14-22(18-24)28(10,11)12/h13-18H,1-12H3/b30-29-
InChIKey=SOCNVTNVHBWFKC-FLWNBWAVSA-N

2b: InChI=1S/C28H42N2/c1-25(2,3)19-13-20(26(4,5)6)16-23(15-19)29-30-24-17-21(27(7,8)9)14-22(18-24)28(10,11)12/h13-18H,1-12H3/b30-29+
InChIKey=SOCNVTNVHBWFKC-QVIHXGFCSA-N

DFT &Schreiner Steven Bachrach 04 Jan 2016 No Comments

Dynamics in the reaction of tetrazine with cyclopropene

Houk and Doubleday report yet another example of dynamic effects in reactions that appear to be simple, ordinary organic reactions.1 Here they look at the Diels-Alder reaction of tetrazine 1 with cyclopropene 2. The reaction proceeds by first crossing the Diels-Alder transition state 3 to form the intermediate 4. This intermediate can then lose the anti or syn N2, through 5a or 5s, to form the product 6. The structures and relative energies, computed at M06-2X/6-31G(d), of these species are shown in Figure 1.

3
17.4

4
-33.2

5a
-28.9

5s
-20.0

6
-86.2

Figure 1. M06-2X/6-31G(d) optimized geometries and energies (relative to 1 + 2) of the critical points along the reaction of tetrazine with cyclopropene.

The large difference in the activation barriers between crossing 5a and 5s (nearly 9 kcal mol-1) suggests, by transition state theory, a preference of more than a million for loss of the anti N2 over the syn N2. However, quasiclassical trajectory studies, using B3LYP/6-31G(d), finds a different situation. The anti pathway is preferred, but only by a 4:1 ratio! This dynamic effect arises from a coupling of the v3 mode which involves a rocking of the cyclopropane ring that brings a proton near the syn N2 functionality, promoting its ejection. In addition, the trajectory studies find short residence times within the intermediate neighborhood for the trajectories that lead to the anti product and longer residence times for the trajectories that lead to the syn product. All together, a very nice example of dynamic effects playing a significant role in a seemingly straightforward organic reaction.

References

(1) Törk, L.; Jiménez-Osés, G.; Doubleday, C.; Liu, F.; Houk, K. N. "Molecular Dynamics of the Diels–Alder Reactions of Tetrazines with Alkenes and N2 Extrusions from Adducts," J. Am. Chem. Soc. 2015, 137, 4749-4758, DOI: 10.1021/jacs.5b00014.

InChIs

1: InChI=1S/C2H2N4/c1-3-5-2-6-4-1/h1-2H
InChIKey=HTJMXYRLEDBSLT-UHFFFAOYSA-N

2: InChI=1S/C3H4/c1-2-3-1/h1-2H,3H2
InChIKey=OOXWYYGXTJLWHA-UHFFFAOYSA-N

4: InChI=1S/C5H6N4/c1-2-3(1)5-8-6-4(2)7-9-5/h2-5H,1H2
InChIKey=JGSMBFYJCNPYDM-UHFFFAOYSA-N

6: InChI=1S/C5H6N2/c1-4-2-6-7-3-5(1)4/h2-5H,1H2
InChIKey=RYJFHKGQZKUXEH-UHFFFAOYSA-N

Diels-Alder &Dynamics &Houk Steven Bachrach 09 Nov 2015 No Comments

Diels-Alder of yne-diyne species

Cramer, Hoye, Kuwata and coworkers have examined the intramolecular cyclization of an alkyne with a diyne.1 Their model system is 1, which can cyclize through a concerted transition state TSC togive the benzyne product 2, or it can proceed through a stepwise pathway, first going through TS1 to form the intermediate INT¸ before traversing through a second transition state TS2 and on to product 2. Using both computations and experiments, they examined a series of compounds with
differing substituents at the ends of the two yne moieties.

The experiments show almost the exact same rate of reaction regardless of the terminal substituents. This includes the parent case where the terminal substituents are hydrogens and also the case where the terminal substituents (which end up on adjacent centers on the benzyne ring) are bulky TMS groups. And though there is no real rate effect due to changes in solvent or the presence of light or triplet oxygen, which suggest a concerted reaction, these substituent effects argue for a step wise reaction.

SMD(o-dichlorobenzene)/B3LYP-D3BJ/6-311+G-(d,p)//M06-2X/6-311+G(d,p)
computations help explain these observations. Shown in Figure 1 are the optimized geometries and relative energies of the critical points on the reaction surface for the conversion of 1 into 2, and these results are similar with the other substituents as well.

1
(0.0)

2
(-56.9)

TSC
(31.5)

 

TS1
(25.5)

INT
(18.8)

TS2
(18.1)

 

Figure 1. SMD(o-dichlorobenzene)/B3LYP-D3BJ/6-311+G-(d,p)//M06-2X/6-311+G(d,p) optimized geometries and relative energies (kcal mol-1).

The first thing to note is that the concerted TSC is higher in energy than the stepwise TS1. The wavefunction for TSC unstable towards moving from a restricted to unrestricted formalism. Reoptimization of some of these concerted TSs actually led to the stepwise TS.

The next item of note is that TS2 for this case is actually lower in energy than the intermediate (this is a true TS on the energy surface, but when zero-point energy and other thermal corrections are included, it becomes lower in energy than INT). For some of the cases the second TS lies above the intermediate, but typically by a small amount.

Therefore, the mechanism of this reaction is stepwise, but the second step might have such a small barrier (or even no barrier) that one might consider this to be concerted, though highly asymmetric and really bearing little resemblance to more traditional concerted pericyclic reactions.

The authors obliquely hinted at some potential interesting dynamics. I suspect that molecular dynamics calculations will show no effect of that second TS, and one might observe some interesting dynamics, which could be seen in kinetic isotope experiments.

References

(1)  Marell, D. J.; Furan, L. R.; Woods, B. P.; Lei, X.; Bendelsmith, A. J.; Cramer, C. J.; Hoye, T. R.; Kuwata, K. T. "Mechanism of the Intramolecular Hexadehydro-Diels–Alder Reaction," J. Org. Chem. 2015 ASAP, DOI: 10.1021/acs.joc.5b01356.

InChIs

1: InChI=1S/C8H4O2/c1-3-5-6-7-10-8(9)4-2/h1-2H,7H2
InChIKey=MGXDIFXPYGGQLF-UHFFFAOYSA-N

2: InChI=1S/C8H4O2/c9-8-7-4-2-1-3-6(7)5-10-8/h2,4H,5H2
InChIKey=MYFORDRJCVOBTH-UHFFFAOYSA-N

Cramer &Diels-Alder &diradicals Steven Bachrach 05 Oct 2015 2 Comments

Next Page »