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.

Finding reliable methodology for optical rotation and correct predictions of (s)-methyloxirane and (1R,5R)-β-pinene

Theoretical calculations of optical rotation (OR), although important to predict absolute configurations (ACs) and corroborate experiments, require efficient methodology able to reproduce enantiomer specificity and real OR values. Also, troublesome molecules are recurring in the literature, such as (S)-methyloxirane and (1R,5R)-β-pinene. This study evaluates DFT functionals B3LYP, CAM-B3LYP, ωB97X-D, M06-2X, and PBE0 considering basis sets aug-cc-pVDZ, aug-cc-pVTZ, 6-311++G(2d,p), and 6-311++G(3df,2p) in OR prediction of 42 rigid organic molecules assessing cases with wrong enantiomeric determination comparing to available experimental data at wavelengths 355, 589, and 633 nm. Functionals CAM-B3LYP and ωB97X-D with aug-cc-pVTZ are indicated here to reproduce experimental values more accurately considering fewer number of wrong AC predictions, normalized RMSD values below 0.70, and a good approximation to experimental values in hierarchical cluster analysis. Methyloxirane AC was reproduced in CAM-B3LYP and PBE0, with [ α ] 355 = 6 . 94 $$ {\left[\alpha \right]}_{355}=6.94 $$ for CAM-B3LYP/aug-cc-pVTZ close to experimental value [ α ] 355 = 7 . 49 ± 0 . 30 $$ {\left[\alpha \right]}_{355}=7.49\pm 0.30 $$ . Good results were found for AC of β-pinene in M06-2X, CAM-B3LYP, and ωB97X-D while the latter in 6-311++G(3df,2p) obtained OR values of [ α ] 589 = 3 . 44 $$ {\left[\alpha \right]}_{589}=3.44 $$ and [ α ] 689 = 4 . 20 $$ {\left[\alpha \right]}_{689}=4.20 $$ close to experimental values [ α ] 589 = 2 . 8 $$ {\left[\alpha \right]}_{589}=2.8 $$ and [ α ] 689 = 4 . 66 ± 0 . 60 $$ {\left[\alpha \right]}_{689}=4.66\pm 0.60 $$ . The two molecules aforementioned are, for the first time, reported to give valid theoretical OR values in such simple methodologies. OR calculations were all performed after geometry optimization at the same level of theory, and analysis of different functional combinations for each step in β-pinene showed it can interfere with AC prediction even in rigid molecules.

Origin invariant electronic circular dichroism in the length dipole gauge without London atomic orbitals

We present a method for obtaining origin-independent electronic circular dichroism (ECD) in the length-gauge representation LG(OI) without the usage of London atomic orbitals. This approach builds upon the work by Caricato [J. Chem. Phys. 153, 151101 (2020)] and is applied to rotatory strengths and ECD spectra from damped response theory. Numerical results are presented for time-dependent Hartree-Fock and density-functional theory, the second-order algebraic diagrammatic construction method, and linear-response coupled-cluster theory with singles and approximate doubles. We can support the finding that the common choice of placing the gauge origin in the center of mass of a molecule in conventional length-gauge calculations involving chiroptical properties might not be optimal and show that LG(OI) is a valuable alternative for the origin-independent calculation of ECD spectra. We show that, for a limited test set, the convergence of the rotatory strengths calculated with the LG(OI) approach toward the basis-set limit tends to be faster than for the established velocity gauge representation. Relationships between the sum-over-states expression of the optical rotation in the LG(OI) framework and its representation in terms of response functions are analyzed.

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.

Gauge Dependence of the S̃ Molecular Orbital Space Decomposition of Optical Rotation

The measurement of optical rotation (OR) is a foundational technique for the detection and characterization of chiral molecules, but it is poorly understood how the observed property relates to the structure of the molecule. Over the years, several schemes have been developed to decompose the OR into more chemically intuitive contributions. In this paper, we introduce two alternative formulations of our previously developed S̃ molecular orbital space decomposition. These new expressions use the modified velocity gauge-magnetic (MVG-M) and -electric (MVG-E) definitions of OR, rather than the length gauge-magnetic (LG-M) definition used in the original paper. Comparing these formulations across a small set of previously studied chiral molecules, we find that these different definitions produce consistent physical interpretations of the OR. These results demonstrate that the S̃ methodology for investigations of structure-property relationships in chiral molecules is insensitive to the choice of gauge.

Full optical rotation tensor at coupled cluster with single and double excitations level in the modified velocity gauge

This work presents the first simulations of the full optical rotation (OR) tensor at coupled cluster with single and double excitations (CCSD) level in the modified velocity gauge (MVG) formalism. The CCSD-MVG OR tensor is origin independent, and each tensor element can in principle be related directly to experimental measurements on oriented systems. We compare the CCSD results with those from two density functionals, B3LYP and CAM-B3LYP, on a test set of 22 chiral molecules. The results show that the functionals consistently overestimate the CCSD results for the individual tensor components and for the trace (which is related to the isotropic OR), by 10%-20% with CAM-B3LYP and 20%-30% with B3LYP. The data show that the contribution of the electric dipole-magnetic dipole polarizability tensor to the OR tensor is on average twice as large as that of the electric dipole-electric quadrupole polarizability tensor. The difficult case of (1S,4S)-(-)-norbornenone also reveals that the evaluation of the former polarizability tensor is more sensitive than the latter. We attribute the better agreement of CAM-B3LYP with CCSD to the ability of this functional to better reproduce electron delocalization compared with B3LYP, consistent with previous reports on isotropic OR. The CCSD-MVG approach allows the computation of reference data of the full OR tensor, which may be used to test more computationally efficient approximate methods that can be employed to study realistic models of optically active materials.

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.

Origin invariant optical rotation in the length dipole gauge without London atomic orbitals

We present an approach to perform origin invariant optical rotation calculations in the length dipole gauge without recourse to London atomic orbitals, called origin invariant length gauge [LG(OI)]. The LG(OI) approach works with any approximate wave function or density functional method, but here we focus on the implementation with the coupled cluster (CC) with single and double excitations method because of the lack of production-level alternatives. Preliminary numerical tests show the efficacy of the LG(OI) procedure and indicate that putting the origin in the center of mass of a molecule may not be an optimal choice for conventional CC-LG calculations.

The ONIOM/PMM Model for Effective Yet Accurate Simulation of Optical and Chiroptical Spectra in Solution: Camphorquinone in Methanol as a Case Study

This paper deals with the development and first validation of a composite approach for the simulation of chiroptical spectra in solution aimed to strongly reduce the number of full QM computations without any significant accuracy loss. The approach starts from the quantum mechanical computation of reference spectra including vibrational averaging effects and taking average solvent effects into account by means of the polarizable continuum model. Next, the snapshots of classical molecular dynamics computations are clusterized and one reference configuration from each cluster is used to compute a reference spectrum. Local fluctuation effects within each cluster are then taken into account by means of the perturbed matrix model. The performance of the proposed approach is tested on the challenging case of the optical and chiroptical spectra of camphorquinone in methanol solution. Although further validations are surely needed, the results of this first study are quite promising also taking into account that agreement with experimental data is reached by just a couple of full quantum mechanical geometry optimizations and frequency computations.

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.