Journal of Physical Chemistry A, Vol.105, No.22, 5356-5371, 2001
Calculation of optical rotation using density functional theory
We report calculations of the frequency-dependent electric dipole-magnetic dipole polarizability tenser, beta (alpha beta)(nu), using ab initio density functional theory (DFT). Gauge invariant (including) atomic orbitals (GIAOs) are used to guarantee origin-independent values of beta = (1/3)Tr [beta (alpha beta)] Calculations of beta at the sodium D fine frequency, beta (D), for 30 rigid chiral molecules are used to predict their specific rotations, [alpha](D). Calculations have been carried out using the B3LYP functional and the 6-31G*, DZP, 6-311++G(2d,2g), aug-cc-pVDZ, and aug-cc-pVTZ basis sets. Comparison to experimental [alpha](D) values for 28 of the 30 molecules yields average absolute deviations of calculated and experimental [alpha](D) values in the range 20-250 for the three large basis sets, all of which include diffuse functions. The accuracies of [alpha](D) values calculated using the 6-31G* and DZP basis sets, which do not include diffuse functions, are significantly lower: average deviations from experiment are 33 degrees and 43 degrees, respectively. Hartree-Fock/Self-Consistent Field (HF/SCF) calculations have been carried out in parallel. HF/SCF [alpha](D) values are substantially lower in accuracy than corresponding B3LYP values; at the aug-cc-pVDZ basis set level, the average deviation from experiment is 63 degrees. [alpha](D) values obtained using beta values calculated in the static limit (nu = 0) are also of lower accuracy than values obtained using beta (D) Absolute Configurations of chiral molecules can be assigned by comparison of predicted and experimental optical rotations. Our results demonstrate that DFT provides substantially more accurate rotations than HF/SCF methodologies employed heretofore and therefore constitutes the current method of choice for stereochemical applications.