On the importance of vibrational contributions to small-angle optical rotation: Fluoro-oxirane in gas phase and solution

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.

Structure Elucidation and Confirmation of Phloroglucinols from the Roots of Garcinia dauphinensis by Comparison of Experimental and Calculated ECD Spectra and Specific Rotations

Eight phloroglucinols from Garcinia dauphinensis were recently reported to have good to moderate antiplasmodial and anticancer activities, consistent with other phloroglucinol derivatives isolated from natural sources. Chiroptical properties were previously calculated and compared to experimental data for compound 2 as a means to deduce its absolute configuration. Tentative assignments for the remaining compounds were also reported based on these data. In order to arrive at stereochemical assignments for phloroglucinols 1 and 3-8, ECD spectra and specific rotations were computed for all stereoisomers of each compound. Molecular orbital analyses were also carried out for the most energetically favorable conformers of each compound. Absolute configurations are reported for all eight phloroglucinols for the first time.

Modeling Complex Solvent Effects on the Optical Rotation of Chiral Molecules: A Combined Molecular Dynamics and Density Functional Theory Study

The challenge of assigning the absolute stereochemical configuration to a chiral compound can be overcome via accurate ab initio predictions of optical rotation, a sensitive molecular property that is further complicated by solvent effects. The solvent's "chiral imprint"-the transfer of the chirality from the solute to the surrounding achiral solvent-is explored here using conformational averaging and time-dependent density-functional theory. These complex solvent effects are taken into account via simple averaging over a molecular dynamics trajectory together with the explicit quantum mechanical consideration of the solvent molecules within the solute's cybotactic region and implicit modeling of the bulk solvent. We consider several axes along which the system's optical rotation varies, including the sampling of the dynamical trajectory, the quality of the one-electron basis set, and the use of continuum solvent models to account for bulk effects.

A Density Functional Theory Study of Optical Rotation in Some Aziridine and Oxirane Derivatives

We present time-dependent density functional theory (TDDFT) calculations of the electronic optical rotation (ORP) for seven oxirane and two aziridine derivatives in the gas phase and in solution and compare the results with the available experimental values. For seven of the studied molecules it is the first time that their optical rotation was studied theoretically and we have therefore investigated the influence of several settings in the TDDFT calculations on the results. This includes the choice of the one-electron basis set, the exchange-correlation functional or the particular polarizable continuum model (PCM). We can confirm that polarized quadruple zeta basis sets augmented with diffuse functions are necessary for converged results and find that the aug-pc-3 basis set is a viable alternative to the frequently employed aug-cc-pVQZ basis set. Based on our study, we cannot recommend the generalized gradient functional KT3 for calculations of the ORP in these compounds, whereas the hybrid functional PBE0 gives results quite similar to the long-range correct CAM-B3LYP functional. Finally, we observe large differences in the solvent effects predicted by the integral equation formalism of PCM and the SMD variant of PCM. For the majority of solute/solvent combinations in this study, we find that the SMD model in combination with the PBE0 functional and the aug-pc-3 basis set gives the best agreement with the experimental values.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Accurate calculation of optically induced birefringences in chiral systems using efficient polarized basis sets

Using state-of-the-art ab initio methodology, we evaluate universal molecular parameters entering the expressions for various optically induced birefringences in chiral fluids. For this, we use the single and double excitation coupled cluster (CCSD) theory together with Dunning's augmented correlation consistent polarized basis sets of increasing size. As this is the first time these parameters are evaluated for chiral molecules using the CCSD approach, we choose possibly small test systems: a model asymmetric methane molecule, and (R)-fluoro-oxirane. With this choice, the convergence of the molecular parameters with the increase of the basis set size is investigated in detail. The results are compared to those obtained with the LPol-n (n = ds, dl, fs) and the ORP basis sets, and to the corresponding Density Functional Theory (DFT) counterparts. We can conclude that the DFT medium constant values are considerably far from coupled cluster, and therefore, aware of the known excellent performance of the CCSD method in the evaluation of various dynamic responses, we do not recommend the former methodology for accurate evaluation of the present properties. Regarding basis set convergence, the performace of the LPol-ds basis set is more efficient than that of the d-aug-cc-pVDZ set, and therefore, the former basis set can be a good choice when dealing with the study of larger systems.

New Basis Set for the Evaluation of Specific Rotation in Flexible Biological Molecules in Solution

A detailed theoretical investigation of specific rotation is carried out in solution for nine flexible molecules of biological importance. Systematic search for the main conformers is followed by time-dependent density functional theory (TD-DFT) calculations of specific rotation employing a wide range of basis sets. Due to conformational flexibility of the compounds under study, the possibility of basis set size reduction without deterioration of the results is investigated. The increasing size (d-)aug-cc-pVXZ (X = D, T, Q) bases of Dunning et al., and the ORP basis set, recently developed to efficiently provide molecular specific rotation, are used for this purpose. The polarizable continuum model is employed at all steps of the investigation. Comparison of the present results with the available data obtained in a vacuum reveals considerable differences, the values in solution being much closer to the experimental specific rotation data available. The ORP basis set proves to be competitive with the d-aug-cc-pVDZ set of Dunning in specific rotation calculations carried out in solution. While having the same number of functions, the former yields, in general, results considerably closer to the reference triple-ζ values. We can thus recommend the ORP basis set to study the optical rotation in conformationally flexible molecules in solution.

Effective Fully Polarizable QM/MM Approach To Model Vibrational Circular Dichroism Spectra of Systems in Aqueous Solution

We propose a methodology, based on the combination of classical Molecular Dynamics (MD) simulations with a fully polarizable Quantum Mechanical (QM)/Molecular Mechanics (MM)/Polarizable Continuum Model (PCM) Hamiltonian, to calculate Vibrational Circular Dichroism (VCD) spectra of chiral systems in aqueous solution. Polarization effects are included in the MM force field by exploiting an approach based on Fluctuating Charges (FQ). By performing the MD, the description of the solvating environment is enriched by taking into account the dynamical aspects of the solute-solvent interactions. On the other hand, the QM/FQ/PCM calculation of the VCD spectrum ensures an accurate description of the electronic density of the solute and a proper account for the specific interactions in solution. The application of our approach to (R)-methyloxirane and (l)-alanine in aqueous solution gives calculated spectra in remarkable agreement with their experimental counterparts and a substantial improvement with respect to the same spectra calculated with the PCM.

Optical Rotation Calculations for a Set of Pyrrole Compounds

Optical rotation of 14 molecules containing the pyrrole group is calculated by employing both time-dependent density functional theory (TDDFT) with the CAM-B3LYP functional and the second-order approximate coupled-cluster singles and doubles (CC2) method. All optical rotations have been provided using the aug-cc-pVDZ basis set at λ = 589 nm. The two methods predict similar results for both sign and magnitude for the optical rotation of all molecules. The obtained signs are consistent with experiments as well, although several conformers for four molecules needed to be studied to reproduce the experimental sign. We have also calculated excitation energies and rotatory strengths for the six lowest lying electronic transitions for several conformers of the two smallest molecules and found that each rotatory strength has various contributions for each conformer which can cause different optical rotations for different conformers of a molecule. Our results illustrate that both methods are able to reproduce the experimental optical rotations, and that the CAM-B3LYP functional, the least computationally expensive method used here, is an applicable and reliable method to predict the optical rotation for these molecules in line with previous studies.

Ro-vibrational averaging of the isotropic hyperfine coupling constant for the methyl radical

We present the first variational calculation of the isotropic hyperfine coupling constant of the carbon-13 atom in the CH3 radical for temperatures T = 0, 96, and 300 K. It is based on a newly calculated high level ab initio potential energy surface and hyperfine coupling constant surface of CH3 in the ground electronic state. The ro-vibrational energy levels, expectation values for the coupling constant, and its temperature dependence were calculated variationally by using the methods implemented in the computer program TROVE. Vibrational energies and vibrational and temperature effects for coupling constant are found to be in very good agreement with the available experimental data. We found, in agreement with previous studies, that the vibrational effects constitute about 44% of the constant's equilibrium value, originating mainly from the large amplitude out-of-plane bending motion and that the temperature effects play a minor role.

Atomic Cholesky decompositions: a route to unbiased auxiliary basis sets for density fitting approximation with tunable accuracy and efficiency

Cholesky decomposition of the atomic two-electron integral matrix has recently been proposed as a procedure for automated generation of auxiliary basis sets for the density fitting approximation [F. Aquilante et al., J. Chem. Phys. 127, 114107 (2007)]. In order to increase computational performance while maintaining accuracy, we propose here to reduce the number of primitive Gaussian functions of the contracted auxiliary basis functions by means of a second Cholesky decomposition. Test calculations show that this procedure is most beneficial in conjunction with highly contracted atomic orbital basis sets such as atomic natural orbitals, and that the error resulting from the second decomposition is negligible. We also demonstrate theoretically as well as computationally that the locality of the fitting coefficients can be controlled by means of the decomposition threshold even with the long-ranged Coulomb metric. Cholesky decomposition-based auxiliary basis sets are thus ideally suited for local density fitting approximations.