Theoretical Study of the Isotope Effect in Optical Rotation.

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, and 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.

First Principles Simulations of Optical Rotation of Chiral Molecular Crystals.

In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT-PBC). Calculations are compared with experimental measurements and show semi-quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X-ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.

On the choice of coordinate origin in length gauge optical rotation calculations.

In this work, we explore the issue of origin dependence in optical rotation (OR) calculations in the length dipole gauge (LG) using standard approximate methods belonging to density functional theory (DFT) and coupled cluster (CC) theory. We use the origin-invariant LG approach, LG(OI), that we recently proposed as reference for the calculations, and we study whether a proper choice of coordinate origin and molecular orientation can be made such that diagonal elements of the LG-OR tensor match those of the LG(OI) tensor. Using a numerical search algorithm, we show that multiple spatial orientations can be found where the LG and LG(OI) results match. However, a simple analytical procedure provides a spatial orientation where the origin of the coordinate system is close to the center of mass of the molecule. At the same time, we also show that putting the origin at the center of mass is not an ideal choice for every molecule (relative errors in the OR up to 70% can be obtained in out test set). Finally, we show that the choice of coordinate origin based on the analytical procedure is transferable across different methods and it is superior to putting the origin in the center of mass or center of nuclear charge. This is important because the LG(OI) approach is trivial to implement for DFT, but not necessarily for nonvariational methods in the CC family. Therefore, one can determine an optimal coordinate origin at DFT level and use it for standard LG-CC response calculations.

Basis Set Dependence of Optical Rotation Calculations with Different Choices of Gauge.

In this work, the basis set dependence of optical rotation (OR) calculations is examined for various choices of gauge/level of theory. The OR is calculated for a set of 50 molecules using B3LYP and CAM-B3LYP and 17 molecules using coupled cluster with single and double excitations (CCSD). The calculations employ the correlation-consistent basis sets, aug-cc-pVζZ with ζ = D, T, Q. An inverse-power extrapolation formula is then utilized to obtain OR values at the complete basis set (CBS) limit. We investigate the basis set convergence for these methods and three choices of gauge: length gauge (with gauge-including atomic orbitals, LG(GIAOs), for DFT), the origin-invariant length gauge [LG(OI)], and the modified velocity gauge (MVG). The results show that all methods converge smoothly to the CBS limit and that the LG(OI) approach has a slightly faster convergence rate than the other choices of gauge. While the DFT methods reach gauge invariance at the CBS limit, CCSD does not. The significant difference between the MVG and LG(OI) results at the CBS limit, 26%, indicates that CCSD is not quite at convergence in the description of electron correlation for this property. On the other hand, gauge invariance at the CBS limit for DFT does not lead to the same OR values for the two density functionals, which is also due to electron correlation incompleteness. A limited comparison to gas-phase experimental OR values for the DFT methods shows that CAM-B3LYP seems more accurate than B3LYP. Overall, this study shows that the LG(OI) approach with the aug-cc-pVTZ basis set for DFT, and with the CBS(DT) extrapolation for CCSD, provides a good cost/accuracy balance.