Modelling kinetic isotope effects for Swern oxidation using DFT-based transition state theory†

Abstract

We investigate the model reported by Giagou and Meyer in 2010 for comparing deuterium kinetic isotope effects (KIEs) computed using density functional theory (DFT) for the intramolecular hydrogen transfer step in the mechanism of the Swern oxidation of alcohols to aldehydes, with those measured by experiment. Whereas the replication of the original computed values for the gas-phase reaction proved entirely successful, several issues were discovered when a continuum solvent model was used. These included uncertainty regarding the parameters and methods used for the calculations and also the coordinates for the original reactant and transition structures, via their provision as data in the ESI. The original conclusions, in which a numerical mis-match between the magnitude of the computed and experimentally measured KIE was attributed to significant deviations from transition structure theory, are here instead rationalised as a manifestation of basis-set effects in the computation. Transition state theory appears to be operating successfully. We now recommend the use of basis sets of triple- or quadruple-ζ quality, rather than the split-valence level previously employed, that dispersion energy corrections be included and that a continuum solvent model using smoothed reaction cavities is essential for effective geometry optimisation and hence accurate normal coordinate analysis. An outlying experimental KIE obtained for chloroform as solvent is attributed to a small level of an explicit hydrogen bonded interaction with the substrate. A temperature outlier for the measured KIE at 195 K is suggested for further experimental investigation, although it may also be an indication of an unusually abrupt incursion of hydrogen tunnelling which would need non-Born–Oppenheimer methods in which nuclear quantum effects are included to be more accurately modelled. We predict KIEs for new substituents, of which those for R = NMe₂ are significantly larger than those for R = H. This approach could be useful in designing variations of the Swern reagent that could lead to synthesis of aldehydes incorporating much higher levels of deuterium. The use of FAIR data rather than the traditional model of its inclusion in the ESI is discussed, and two data discovery tools exploiting these FAIR attributes are suggested.