The ORP basis set designed for optical rotation calculations

Reduced Scaling Real-Time Coupled Cluster Theory

Real-time coupled cluster (CC) methods have several advantages over their frequency-domain counterparts, namely, response and equation of motion CC theories. Broadband spectra, strong fields, and pulse manipulation allow for the simulation of complex spectroscopies that are unreachable using frequency-domain approaches. Due to the high-order polynomial scaling, the required numerical time propagation of the CC residual expressions is a computationally demanding process. This scaling may be reduced by local correlation schemes, which aim to reduce the size of the (virtual) orbital space by truncation according to user-defined parameters. We present the first application of local correlation to real-time CC. As in previous studies of locally correlated frequency-domain CC, traditional local correlation schemes are of limited utility for field-dependent properties; however, a perturbation-aware scheme proves promising. A detailed analysis of the amplitude dynamics suggests that the main challenge is a strong time dependence of the wave function sparsity.

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.

Comparative Study of Vibrational Raman Optical Activity with Different Time-Dependent Density Functional Approximations: The VROA36 Database

A new database, VROA36, is introduced to investigate the performance of computational approaches for vibrational Raman optical activity (VROA) calculations. The database is composed of 36 molecules with known experimental VROA spectra. It includes 93 conformers. Normal modes calculated with B3LYP-D3(BJ)/def2-TZVP are used to compute the VROA spectra with four functionals, B3LYP-D3(BJ), ωB97X-D, M11, and optimally tuned LC-PBE, as well as several basis sets. SimROA indices and frequency scaling factors are used to compare calculated spectra with each other and with experimental data. The four functionals perform equally well independently of the basis set and usually achieve good agreement with the experimental data. For molecules in near- or at-resonance conditions, the inclusion of a complex (damped) linear response approach is important to obtain physically meaningful VROA intensities. The use of any of the tested functional approximations with the def2-SVPD Gaussian-type basis set, or a basis of similar flexibility, can be recommended for efficient and reliable theoretical VROA studies.

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.

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.

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.

New basis set for the prediction of the specific rotation in flexible biological molecules

A novel method based on increasingly accurate calculations is used to obtain the main conformers of a set of nine flexible molecules. Their specific rotation is evaluated using the recently developed ORP basis set.

Comparison of Property-Oriented Basis Sets for the Computation of Electronic and Nuclear Relaxation Hyperpolarizabilities

In the present work, we perform an assessment of several property-oriented atomic basis sets in computing (hyper)polarizabilities with a focus on the vibrational contributions. Our analysis encompasses the Pol and LPol-ds basis sets of Sadlej and co-workers, the def2-SVPD and def2-TZVPD basis sets of Rappoport and Furche, and the ORP basis set of Baranowska-Łączkowska and Łączkowski. Additionally, we use the d-aug-cc-pVQZ and aug-cc-pVTZ basis sets of Dunning and co-workers to determine the reference estimates of the investigated electric properties for small- and medium-sized molecules, respectively. We combine these basis sets with ab initio post-Hartree-Fock quantum-chemistry approaches (including the coupled cluster method) to calculate electronic and nuclear relaxation (hyper)polarizabilities of carbon dioxide, formaldehyde, cis-diazene, and a medium-sized Schiff base. The primary finding of our study is that, among all studied property-oriented basis sets, only the def2-TZVPD and ORP basis sets yield nuclear relaxation (hyper)polarizabilities of small molecules with average absolute errors less than 5.5%. A similar accuracy for the nuclear relaxation (hyper)polarizabilites of the studied systems can also be reached using the aug-cc-pVDZ basis set (5.3%), although for more accurate calculations of vibrational contributions, i.e., average absolute errors less than 1%, the aug-cc-pVTZ basis set is recommended. It was also demonstrated that anharmonic contributions to first and second hyperpolarizabilities of a medium-sized Schiff base are particularly difficult to accurately predict at the correlated level using property-oriented basis sets. For instance, the value of the nuclear relaxation first hyperpolarizability computed at the MP2/def2-TZVPD level of theory is roughly 3 times larger than that determined using the aug-cc-pVTZ basis set. We link the failure of the def2-TZVPD basis set with the difficulties in predicting the first-order field-induced coordinates. On the other hand, the aug-cc-pVDZ and ORP basis sets overestimate the property in question only by roughly 30%. In this study, we also propose a low-cost composite treatment of anharmonicity that relies on the combination of two basis sets, i.e., a large-sized basis set is employed to determine lowest-order derivatives with respect to the field-induced coordinates, and a medium-sized basis set is used to compute the higher-order derivatives. The results of calculations performed at the MP2 level of theory demonstrate that this approximate scheme is very successful at predicting nuclear relaxation hyperpolarizabilities.