azatriquinacene, a novel aromatic

Aromaticity Steven Bachrach 11 Dec 2017 No Comments

The range of aromatic compounds seems limitless. Mascal and co-workers have prepared the azatriquinacene 1 in a remarkably simple fashion.1 The molecule is a zwitterion, with the carbon atoms forming a 9-center, but 10 π-electron ring, and the quaternary nitrogen sitting above it. The carbon ring satisfies Hückel’s rule (4n+2) and so should be aromatic. The capping nitrogen should help to keep the carbon ring fixed in a shallow bowl.

As expected, the molecule in fact turns out to possess an aromatic 10 π-electron ring. The B3LYP/6-311++G(d,p) geometry is shown in Figure 1. There is little bond alternation among the C-C distances: the mean deviation is only 0.015 Å with the largest difference only 0.024 &Aring. The x-ray crystal structure shows the same trends. The NICS(1) value is -12.31 ppm, larger even than that of benzene (-10.22 ppm).

Figure 1. B3LYP/6-311++G(d,p) geometry of 1.

References

1) Hafezi, N.; Shewa, W. T.; Fettinger, J. C.; Mascal, M., "A Zwitterionic, 10 π Aromatic Hemisphere." Angew. Chem. Int. Ed. 2017, 56, 14141-14144, DOI: 10.1002/anie.201708521.

InChIs

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

Heavy-atom tunneling

Tunneling Steven Bachrach 04 Dec 2017 1 Comment

Though recognized to occur in organic systems, the breadth of involvement of heavy-atom tunneling has not been established. Doubleday, Greer and coworkers have examined 13 simple organic reactions sampling pericyclic reactions, radical rearrangements and SN2 reactions for heavy-atom tunneling.1 A few of these reactions are shown below.

Reaction rates were obtained using the small curvature tunneling approximation (SCT), computed using Gaussrate. Reaction surfaces were computed at B3LYP/6-31G*. The tunneling correction to the rate was also estimated using the model developed by Bell: kBell = (u/2)/sin(u/2) where u = hν/RT and ν is the imaginary frequency associated with the transition state. The temperature was chosen so as to give a common rate constant of 3 x 10-5 s-1. Interestingly, all of the examined reactions exhibited significant tunneling even at temperatures from 270-350 K (See Table 1). The tunneling effect estimated by Bell’s equation is very similar to that of the more computationally demanding SCT computation.

Table 1. Tunneling contribution to the rate constant

Reaction

% tunneling

95

35

17

28

CN + CH3Cl → CH3CN + Cl (aqueous)

45

This study points towards a much broader range of reactions that may be subject to quantum mechanical tunneling than previously considered.

(Note: The original post had swapped some of the values in Table 1. These have now been corrected. My thanks to Drs. Greer and Doubleday for bringing this to my attention.)

References

1. Doubleday, C.; Armas, R.; Walker, D.; Cosgriff, C. V.; Greer, E. M., "Heavy-Atom Tunneling Calculations in Thirteen Organic Reactions: Tunneling Contributions are Substantial, and Bell’s Formula Closely Approximates Multidimensional Tunneling at ≥250 K." Angew. Chem. Int. Ed. 2017, 56, 13099-13102, DOI: 10.1002/anie.201708489.

Perspective on Tunneling Control

Schreiner &Tunneling Steven Bachrach 13 Nov 2017 No Comments

Over the past nine years the Schreiner group, often in collaboration with the Allen group, have produced some remarkable studies demonstrating the role of tunneling control. (I have made quite a number of posts on this topics.) Tunneling control is a third mechanism for dictating product formation, in tandem with kinetic control (the favored product is the one that results from the lowest barrier) and thermodynamic control (the favored product is the one that has the lowest energy). Tunneling control has the favored product resulting from the narrowest mass-considered barrier.

Schreiner has written a very clear perspective on tunneling control. It is framed quite interestingly by some fascinating quotes:

It is probably fair to say that many organic chemists view the concept of tunneling, even of hydrogen atoms, with some skepticism. – Carpenter 19832

Reaction processes have been considered as taking place according to the laws of classical mechanics, quantum mechanical theory being only employed in calculating interatomic forces. – Bell 19333

Schreiner’s article makes it very clear how critical it is to really think about reactions from a truly quantum mechanical perspective. He notes the predominance of potential energy diagrams that focus exclusively on the relative energies and omits any serious consideration of the reaction coordinate metrics, like barrier width. When one also considers the rise in our understanding of the role of reaction dynamics in organic chemistry (see, for example, these many posts), just how long will it take for these critical notions to penetrate into standard organic chemical thinking? As Schreiner puts it:

It should begin by including quantum phenomena in introductory textbooks, where they are, at least in organic chemistry, blatantly absent. To put this oversight in words similar to those used much earlier by Frank Weinhold in a different context: “When will chemistry textbooks begin to serve as aids, rather than barriers, to this enriched quantum-mechanical perspective?”4

References

1) Schreiner, P. R., "Tunneling Control of Chemical Reactions: The Third Reactivity Paradigm." J. Am. Chem. Soc. 2017, 139, 15276-15283, DOI: 10.1021/jacs.7b06035.

2) Carpenter, B. K., "Heavy-atom tunneling as the dominant pathway in a solution-phase reaction? Bond shift in antiaromatic annulenes." J. Am. Chem. Soc. 1983, 105, 1700-1701, DOI: 10.1021/ja00344a073.

3) Bell, R. P., "The Application of Quantum Mechanics to Chemical Kinetics." Proc. R. Soc. London, Ser. A 1933, 139 (838), 466-474, DOI: 10.1098/rspa.1933.0031.

4) Weinhold, F., "Chemistry: A new twist on molecular shape." Nature 2001, 411, 539-541, DOI: 10.1038/35079225.

Heavy atom tunneling in semibullvalene

Borden &Tunneling Steven Bachrach 02 Nov 2017 No Comments

Another prediction made by quantum chemistry has now been confirmed. In 2010, Zhang, Hrovat, and Borden predicted that the degenerate rearrangement of semibullvalene 1 occurs with heavy atom tunneling.1 For example, the computed rate of the rearrangement including tunneling correction is 1.43 x 10-3 s-1 at 40 K, and this rate does not change with decreasing temperature. The predicted half-life of 485 s is 1010 shorter than that predicted by transition state theory.

Now a group led by Sander has examined the rearrangement of deuterated 2.2 The room temperature equilibrium mixture of d42 and d22 was deposited at 3 K. IR observation showed a decrease in signal intensities associated with d42 and concomitant growth of signals associated with d22. The barrier for this interconversion is about 5 kcal mol-1, too large to be crossed at this temperature. Instead, the interconversion is happening by tunneling through the barrier (with a rate about 10-4 s-1), forming the more stable isomer d22 preferentially. This is exactly as predicted by theory!

References

1. Zhang, X.; Hrovat, D. A.; Borden, W. T., "Calculations Predict That Carbon Tunneling Allows the Degenerate Cope Rearrangement of Semibullvalene to Occur Rapidly at Cryogenic Temperatures." Org. Letters 2010, 12, 2798-2801, DOI: 10.1021/ol100879t.

2. Schleif, T.; Mieres-Perez, J.; Henkel, S.; Ertelt, M.; Borden, W. T.; Sander, W., "The Cope Rearrangement of 1,5-Dimethylsemibullvalene-2(4)-d1: Experimental Evidence for Heavy-Atom Tunneling." Angew. Chem. Int. Ed. 2017, 56, 10746-10749, DOI: 10.1002/anie.201704787.

InChIs

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

d42: InChI=1S/C10H12/c1-9-5-3-7-8(4-6-9)10(7,9)2/h3-8H,1-2H3/i5D
InChIKey=WUJOLJNLXLACNA-UICOGKGYSA-N

d22: InChI=1S/C10H12/c1-9-5-3-7-8(4-6-9)10(7,9)2/h3-8H,1-2H3/i7D
InChIKey=WUJOLJNLXLACNA-WHRKIXHSSA-N

Review of the Activation Strain/Distortion-Interaction Model

Houk Steven Bachrach 16 Oct 2017 No Comments

Bickelhaupt and Houk present a nice review of their separately developed, but conceptually identical model for assessing reactivity.1 Houk termed this the “distortion/interaction” model,2 while Bickelhaupt named it “activation strain”.3 The concept is that the activation barrier can be dissected in a distortion or stain energy associated with bringing the reactants into the geometry of the transition state, and the interaction energy is the stabilization energy afforded by the molecular orbital interactions of the reactant components with each other in the transition state.

The review discusses a broad range of applications, including SN2 and E2 reactions, pericyclic reactions (including Diels-Alder reactions of enones and the dehdydro Diels-Alder reaction that I have discussed in this blog), a click reaction, a few examples involving catalysis, and the regioselectivity of indolyne (see this post). They also discuss the role of solvent and the relationship of this model to Marcus Theory.

I also want to mention in passing a somewhat related article by Jorgensen and co-authors published in the same issue of Angewandte Chemie as the above review.4 This article discusses the paucity of 10 electron cycloaddition reactions, especially in comparison to the large number of very important cycloaddition reactions involving 6 electrons, such as the Diels-Alder reaction, the Cope rearrangement, and the Claisen rearrangement. While the article does not focus on computational methods, computations have been widely used to discuss 10-electron cycloadditions. The real tie between this paper and the review discussed above is Ken Houk, whose graduate career started with an attempt to perform a [6+4] cycloaddition, and he has revisited the topic multiple times throughout his career.

References

1. Bickelhaupt, F. M.; Houk, K. N., "Analyzing Reaction Rates with the Distortion/Interaction-Activation Strain Model." Angew. Chem. Int. Ed. 2017, 56, 10070-10086, DOI: 10.1002/anie.201701486.

2. Ess, D. H.; Houk, K. N., "Distortion/Interaction Energy Control of 1,3-Dipolar Cycloaddition Reactivity." J. Am. Chem. Soc. 2007, 129, 10646-10647, DOI: 10.1021/ja0734086

3. Bickelhaupt, F. M., "Understanding reactivity with Kohn-Sham molecular orbital theory: E2-SN2 mechanistic spectrum and other concepts." J. Comput. Chem. 1999, 20, 114-128

4. Palazzo, T. A.; Mose, R.; Jørgensen, K. A., "Cycloaddition Reactions: Why Is It So
Challenging To Move from Six to Ten Electrons?" Angew. Chem. Int. Ed. 2017, 56, 10033-10038, DOI: 10.1002/anie.201701085.

More applications of computed NMR spectra

NMR Steven Bachrach 04 Oct 2017 No Comments

In this post I cover two papers discussing application of computed NMR chemical shifts to structure identification and (yet) another review of computational techniques towards NMR structure prediction.

Grimblat, Kaufman, and Sarotti1 take up the structure of rubriflordilactone B 1, which was isolated from Schisandra rubriflora. The compound was then synthesized and its x-ray structure reported, however its NMR did not match with the natural extract. It was suggested that there were actually two compounds in the extract, the minor one was less soluble and is the crystallized 1, and a second compound responsible for the NMR signal.

The authors looked at all stereoisomers of this molecule keeping the three left-most rings intact. The low energy rotamers of these 32 stereoisomers were then optimized at B3LYP/6-31G* and the chemical shifts computed at PCM(pyridine)/mPW1PW91/6-31+G**. To benchmark the method, DP4+ was used to identify which stereoisomer best matches with the observed NMR of authentic 1; the top fit (92.6% probability) was the correct structure.

The 32 stereoisomers were then tested against the experimental NMR of the natural extract. DP4+ with just the proton shifts suggested structure 2 (99.8% probability); however, the 13C chemical shifts predicted a different structure. Re-examination of the reported chemical shifts identifies some mis-assigned signals, which led to a higher C-DP4+ prediction. When all 128 stereoisomers were tested, structure 2 had the highest DP4+ prediction (99.5%), but the C-DP4+ prediction remained problematic (10.8%). Analyzing the geometries of all reasonable alternative for agreement with the NOESY spectrum confirmed 2. These results underscore the importance of using all data sources.

Reddy and Kutateladze point out the importance of using coupling constants along with chemical shifts in structure identification.2 They examined cordycepol A 3, obtained from Cordyceps ophioglossoides. They noted that the computed chemical shifts and coupling constants of originally proposed structure 3a differed dramatically from the experimental values.

They first proposed that the compound has structure 3b. The computed coupling constants using their relativistic force field.3 The experimental coupling constants for the proton H1 are 13.4 and 7.1 Hz. The computed values for 3a are 8.9 and 1.6 Hz, and this structure is clearly incorrect. The coupling constants are improved with 3b, but the 13C chemical shifts are in poor agreement with experiment. So, they proposed structure 3c, the epimer at both C1 and C11 of the original structure.

They optimized four conformations of 3c at B3LYP/6-31G(d) and obtained Boltzmann-weighted chemical shifts at mPW1PW91/6-311+G(d,p). The RMS deviation of the computed 13C chemical shifts relative to the experiment is only 1.54 ppm, and more importantly, the computed coupling constants of 13.54 and 6.90 Hz are in excellent agreement with the experiment values.

Lastly, Grimblat and Sarotti present a review of a number of methods for using computed NMR chemical shifts towards structure prediction.4 These methods include CP3, DP4, DP4+ (all of which I have posted on in the past) and an artificial neural network approach of their own design. They discuss a number of interesting cases where each of these methods has been crucial in identifying the correct chemical structure.

References

1. Grimblat, N.; Kaufman, T. S.; Sarotti, A. M., "Computational Chemistry Driven Solution to Rubriflordilactone B." Org. Letters 2016, 18, 6420-6423, DOI: 10.1021/acs.orglett.6b03318.

2. Reddy, D. S.; Kutateladze, A. G., "Structure Revision of an Acorane Sesquiterpene Cordycepol A." Org. Letters 2016, 18, 4860-4863, DOI: 10.1021/acs.orglett.6b02341.

3. (a) Kutateladze, A. G.; Mukhina, O. A., "Minimalist Relativistic Force Field: Prediction of Proton–Proton Coupling Constants in 1H NMR Spectra Is Perfected with NBO Hybridization Parameters." J. Org. Chem. 2015, 80, 5218-5225, DOI: 10.1021/acs.joc.5b00619; (b) Kutateladze, A. G.; Mukhina, O. A., "Relativistic Force Field: Parametrization of 13C–1H Nuclear Spin–Spin Coupling Constants." J. Org. Chem. 2015, 80, 10838-10848, DOI: 10.1021/acs.joc.5b02001.

4. Grimblat, N.; Sarotti, A. M., "Computational Chemistry to the Rescue: Modern Toolboxes for the Assignment of Complex Molecules by GIAO NMR Calculations." Chem. Eur. J. 2016, 22, 12246-12261, DOI: h10.1002/chem.201601150.

InChIs

1: InChI=1S/C28H30O6/c1-13-9-20(32-26(13)30)25-14(2)24-17-6-5-15-12-28-21(8-7-16(15)18(17)10-19(24)31-25)27(3,4)33-22(28)11-23(29)34-28/h5-9,14,19-22,24-25H,10-12H2,1-4H3/t14-,19+,20-,21-,22+,24-,25-,28+/m0/s1
InChIKey=JGSLSHOXBXVVTQ-NEUKEVNNSA-N

2: InChI=1S/C28H30O6/c1-13-9-20(32-26(13)30)25-14(2)24-17-6-5-15-12-28-21(8-7-16(15)18(17)10-19(24)31-25)27(3,4)33-22(28)11-23(29)34-28/h5-9,14,19-22,24-25H,10-12H2,1-4H3/t14-,19-,20-,21-,22+,24+,25-,28+/m0/s1
InChIKey=JGSLSHOXBXVVTQ-WQIRXNRDSA-N

3c: InChI=1S/C16H28O2/c1-6-11(2)9-14-16(5)12(3)7-8-13(16)15(4,17)10-18-14/h9,12-14,17H,6-8,10H2,1-5H3/b11-9-/t12-,13-,14-,15-,16+/m0/s1
InChIKey=WPQIVUHVYBQTBG-AWEVENECSA-N

Triplet cyclobutadiene

Aromaticity &cyclobutadiene Steven Bachrach 11 Sep 2017 1 Comment

Cyclobutadiene has long fascinated organic chemists. It is the 4e analogue of the 6e benzene molecule, yet it could hardly be more different. Despite nearly a century of effort, cyclobutadiene analogues were only first prepared in the 1970s, reflecting its strong antiaromatic character.

Per-trimethylsilylcyclobutadiene 1 offers opportunities to probe the properties of the cyclobutadiene ring as the bulky substituents diminish dimerization and polymerization of the reactive π-bonds. Kostenko and coworkers have now reported on the triplet state of 1.1 They observe three EPR signals of 1 at temperatures above 350 K, and these signals increase in area with increasing temperature. This is strong evidence for the existence of triplet 1 in equilibrium with the lower energy singlet. Using the variable temperature EPR spectra, the singlet triplet gap is 13.9 ± 0.8 kcal mol-1.

The structures of singlet and triplet 1 were optimized at B3LYP-D3/6-311+G(d,p) and shown in Figure 1. The singlet is the expected rectangle, with distinctly different C-C distance around the ring. The triplet is a square, with equivalent C-C distances. Since both the singlet and triplet states are likely to have multireference character, the energies of both states were obtained at RI-MRDDCI2-CASSCF(4,4)/def2-SVP//B3LYPD3/6-311+G(d,p) and give a singlet-triplet gap of 11.8 kcal mol-1, in quite reasonable agreement with experiment.

singlet

triplet

Figure 1. Optimized geometries of singlet and triplet 1.

References

1. Kostenko, A.; Tumanskii, B.; Kobayashi, Y.; Nakamoto, M.; Sekiguchi, A.; Apeloig, Y., "Spectroscopic Observation of the Triplet Diradical State of a Cyclobutadiene." Angew. Chem. Int. Ed. 2017, 56, 10183-10187, DOI: 10.1002/anie.201705228.

InChIs

1: InChI=1S/C16H36Si4/c1-17(2,3)13-14(18(4,5)6)16(20(10,11)12)15(13)19(7,8)9/h1-12H3
InChIkey=AYOHYRSQVCLGKR-UHFFFAOYSA-N

Diatomic densities from DFT

DFT Steven Bachrach 21 Aug 2017 No Comments

I recently blogged about a paper arguing that modern density functional development has strayed from the path of improving density description, in favor of improved energetics. The Medvedev paper1 was met with a number of criticisms. A potential “out” from the conclusions of the work was that perhaps molecular densities do not fare so poorly with more modern functionals, following the argument that better energies might reflect better densities in bonding regions.

The Hammes-Schiffer group have now examined 14 diatomic molecules with the goal of testing just this hypothesis.2 They subjected both homonuclear diatomics, like N2, Cl2, and Li2, and heteronuclear diatomics, like HF, LiF, and SC, to 90 different density functionals using the very large aug-cc-pCVQZ basis set. Using the CCSD density as a reference, they examined the differences in the densities predicted by the various functional both along the internuclear axis and perpendicular to it.

The 20 functionals that do the best job in mimicking the CCSD density are all hybrid GGA functionals, along with the sole double hybrid functional included in the study (B2PLYP). These functionals date from 1993 to 2012. The 20 functionals that do the poorest job include functionals from all rung-types, and date from 1980-2012. A very slight upward trend can be observed in the density error increasing with development year, while the error in the dissociation energy clearly is decreasing over time.

They note that six functionals of the Minnesota-type, those that are highly parameterized and of recent vintage, perform very poorly at predicting atomic densities, but do well with the diatomic densities.

Hammes-Schiffer concludes that their diatomic results support the general trend noted by Medvedev’s atomic results, that density description is lagging in more recently developed functionals. I’d add that this trend is not as dramatic for the diatomics as for atoms.

They pose what is really the key question: “Is the purpose to approximate the exact functional or simply to provide chemists with a useful tool for exploring chemical systems?” Since, as they note, the modern highly parameterized functionals have worked so well for predicting energies and geometries, “the observation that many modern functionals produce incorrect densities could be of no great consequence for many studies”. Nonetheless, “the ultimate goal is still to obtain both accurate densities and accurate energies”.

References

1) Medvedev, M. G.; Bushmarinov, I. S.; Sun, J.; Perdew, J. P.; Lyssenko, K. A., "Density functional theory is straying from the path toward the exact functional." Science 2017, 355, 49-52, DOI: 10.1126/science.aah5975.

2) Brorsen, K. R.; Yang, Y.; Pak, M. V.; Hammes-Schiffer, S., "Is the Accuracy of Density Functional Theory for Atomization Energies and Densities in Bonding Regions Correlated?" J. Phys. Chem. Lett. 2017, 8, 2076-2081, DOI: 10.1021/acs.jpclett.7b00774.

Bispericyclic reaction involving two [6+4] cycloadditions

cycloadditions &Dynamics &Houk Steven Bachrach 07 Aug 2017 No Comments

Bispericyclic transition states arise when two pericyclic reactions merge to a common transition state. This leads to a potential energy surface with a bifurcation such that reactions that traverse this type of transition state will head towards two different products. The classic example is the dimerization of cyclopentadiene, involving two [4+2] Diels-Alder reactions. Unusual PESs are discussed in my book and in past blog posts.

Houk and coworkers have now identified a bispericyclic transition state involving two [6+4] cycloadditions.1 Reaching back to work Houk pursued as a graduate student with Woodward for inspiration, these authors examined the reaction of tropone 1 with dimethylfulvene 2. Each moiety can act as the diene or triene component of a [6+4] allowed cycloaddition:

The product with fulvene 2 as the 6 π-e component and tropone as the 4 π-e component [6F+4T] is 3, while reversing their participation in the [6T+4F] cycloaddition leads to 4. A variety of [4+2] reactions are also possible. All of these reactions were investigated at PCM/M06-2X/6-311+G(d,p)//B3LYP-D3/6-31G(d). The reaction leading to 3 is exothermic by 3.0 kcal mol-1, while the reaction to 4 endothermic by 1.3 kcal mol-1.

Interestingly, there is only one transition state that leads to both 3 and 4, the first known bispericyclic transition state for two conjoined [6+4] cycloadditions. The barrier is 27.9 kcal mol-1. The structures of the two products and the transition state leading to them are shown in Figure 1. 3 and 4 can interconvert through a Cope transition state, also shown in Figure 1, with a barrier of 26.3 kcal mol-1 (for 43).

3

4

TS [6+4]

TS Cope

Figure 1. B3LYP-D3/6-31G(d) optimized geometries.

Given that a single transition leads to two products, the product distribution is dependent on the molecular dynamics. A molecular dynamics simulation at B3LYP-D3/6-31G(d) with 117 trajectories indicates that 4 is formed 91% while 3 is formed only 9%. Once again, we are faced with the reality of much more complex reaction mechanisms/processes than simple models would suggest.

References

1) Yu, P.; Chen, T. Q.; Yang, Z.; He, C. Q.; Patel, A.; Lam, Y.-h.; Liu, C.-Y.; Houk, K. N., "Mechanisms and Origins of Periselectivity of the Ambimodal [6 + 4] Cycloadditions of Tropone to Dimethylfulvene." J. Am. Chem. Soc. 2017, 139 (24), 8251-8258, DOI: 10.1021/jacs.7b02966.

InChIs

1: InChI=1S/C7H6O/c8-7-5-3-1-2-4-6-7/h1-6H
InChIKey=QVWDCTQRORVHHT-UHFFFAOYSA-N

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

3:InChI=1S/C15H16O/c1-15(2)10-6-8-12(14(16)9-7-10)11-4-3-5-13(11)15/h3-12H,1-2H3
InChIKey=SEKRUGIZAIQCDA-UHFFFAOYSA-N

4: InChI=1S/C15H16O/c1-9(2)14-10-7-8-11(14)13-6-4-3-5-12(10)15(13)16/h3-8,10-13H,1-2H3
InChIKey=AQQAMUGJSGJKLC-UHFFFAOYSA-N

A few review articles

Dynamics &Houk Steven Bachrach 25 Jul 2017 No Comments

A few nice review/opinion pieces have been piling up in my folder of papers of interest for this Blog. So, this post provides a short summary of a number of review articles that computationally-oriented chemists may find of interest.

Holy Grails in computational chemistry

Houk and Liu present a short list of “Holy Grails” in computationally chemistry.1 They begin by pointing out a few technical innovations that must occur for the Grails to be found: development of a universal density functional; an accurate, generic force field; improved sampling for MD; and dealing with the combinatorial explosion with regards to conformations and configurations. Their list of Grails includes predicting crystal structures and structure of amorphous materials, catalyst design, reaction design, and device design. These Grails overlap with the challenges I laid out in my similarly-themed article in 2014.2

Post-transition state bifurcations and dynamics

Hare and Tantillo review the current understanding of post-transition state bifurcations (PTSB).3 This type of potential energy surface has been the subject of much of Chapter 8 of my book and many of my blog posts. What is becoming clear is the possibility of a transition state followed by a valley-ridge inflection leads to reaction dynamics where trajectories cross a single transition state but lead to two different products. This new review updates the state-of-the-art from Houk’s review4 of 2008 (see this post). Mentioned are a number of studies that I have included in this Blog, along with reactions involving metals, and biochemical systems (many of these examples come from the Tantillo lab). They close with the hope that their review might “inspire future studies aimed at controlling selectivity for reactions with PTSBs” (italics theirs). I might offer that controlling selectivity in these types of dynamical systems is another chemical Grail!

The Hase group has a long review of direct dynamics simulations.5 They describe a number of important dynamics studies that provide important new insight to reaction mechanism, such as bimolecular SN2 reactions (including the roundabout mechanism) and unimolecular dissociation. They write a long section on post-transition state bifurcations, and other dynamic effects that cannot be interpreted using transition state theory or RRKM. This section is a nice complement to the Tantillo review.

Benchmarking quantum chemical methods

Mata and Suhm look at our process of benchmarking computational methods.6 They point out the growing use of high-level quantum computations as the reference for benchmarking new methods, often with no mention of any comparison to experiment. In defense of theoreticians, they do note the paucity of useful experimental data that may exist for making suitable comparisons. They detail a long list of better practices that both experimentalists and theoreticians can take to bolster both efforts, leading to stronger computational tools that are more robust at helping to understand and discriminate difficult experimental findings.

References

1) Houk, K. N.; Liu, F., "Holy Grails for Computational Organic Chemistry and Biochemistry." Acc. Chem. Res. 2017, 50 (3), 539-543, DOI: 10.1021/acs.accounts.6b00532.

2) Bachrach, S. M., "Challenges in computational organic chemistry." WIRES: Comput. Mol. Sci. 2014, 4, 482-487, DOI: 10.1002/wcms.1185.

3) Hare, S. R.; Tantillo, D. J., "Post-transition state bifurcations gain momentum – current state of the field." Pure Appl. Chem. 2017, 89, 679-698, DOI: 0.1515/pac-2017-0104.

4) Ess, D. H.; Wheeler, S. E.; Iafe, R. G.; Xu, L.; Çelebi-Ölçüm, N.; Houk, K. N., "Bifurcations on Potential Energy Surfaces of Organic Reactions." Angew. Chem. Int. Ed. 2008, 47, 7592-7601, DOI: 10.1002/anie.200800918

5) Pratihar, S.; Ma, X.; Homayoon, Z.; Barnes, G. L.; Hase, W. L., "Direct Chemical Dynamics Simulations." J. Am. Chem. Soc. 2017, 139, 3570-3590, DOI: 10.1021/jacs.6b12017.

6) Mata, R. A.; Suhm, M. A., "Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction?" Angew. Chem. Int. Ed. 2017, ASAP, DOI: 10.1002/anie.201611308.

Next Page »