Abstract
In this work, the isotope effect in optical rotation (OR) is examined by exploring structure- property relationships for H → D substitutions in chiral molecules. While electronic effects serve as the dominant source of optical activity, there is a non-negligible contribution from nuclear vibrations, which changes with isotopic substitution. We employ a test set of 50 small organic molecules: three-membered rings with varying heteroatoms (PCl, PH, S, NCl, NH, O, NBr) and functional groups (Me, F), and simulations were run at the B3LYP/aug- cc-pVDZ level of theory. The objectives of this work are to determine locations of isotopic substitution that result in significant changes in the vibrational correction to the OR and to evaluate which vibrational modes and electronic response are the major contributors to the isotope effect. Molecules with more polarizable heteroatoms in the ring (e.g., S and P) have the largest change in the vibrational correction compared to the unsubstituted parent molecules. In many cases, isotopic substitution made to the hydrogens on the opposite side of the ring from the functional group provides the largest change in the OR. H/D wagging modes and C vibrations (for D–C centers) are the largest contributors to the isotope effect. This is explained with a molecular orbital decomposition analysis of the OR. The relevant vibrational modes affect the orbital transitions that are already significant at the equilibrium geometry. However, this effect is only large when polarizable heteroatoms are involved, because the electron density surrounding them is diffuse enough to feel the subtle effect of change in mass due to isotopic substitution on the relevant vibrational modes.
Supplementary materials
Title
Supporting Information for: Theoretical Study of the Isotope Effect in Optical Rotation
Description
The Supporting Information includes: the calculated values of [α]eq, [α]V1, and [α]V2 for all test molecules and isotopic substitutions, Tables S1-S7; the corresponding individual normal mode contributions, Tables S8-S23; and localized mode contributions, Tables S24- S31; the [α]V1 and [α]V2 normal mode plots for each isotopologue of molecules S−1F , PCl−1F−Diast.2, S−1Me and PCl−1Me−Diast.2, Figures S1-S2; the [α]V1 and [α]V2 localized mode plots for each isotopologue of the same four molecules, Figures S3-S4; the ∆Sia heatmap plots for selected normal modes of PCl−1F−Diast.2 along with visualizations of displacement vectors, Figure S5; the ∆Sia heatmap plots for selected localized modes of PCl−1F−Diast.2, Figure S6; the ∆S ̃ia heatmap plots for selected normal modes of PCl−1Me−Diast.2 along with visualizations of displacement vectors, Figure S7; the ∆Sia heatmap plots for selected localized modes of PCl−1Me−Diast.2, Figure S8.
Actions
Title
Optimized geometries
Description
This file contains the optimized geometries in xyz file format of the 50 molecules used in this work.
Actions