A Pragmatic Recipe for the Treatment of Hindered Rotations in the Vibrational Averaging of Molecular Properties

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.

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.

Towards reliable references for electron paramagnetic resonance parameters based on quantum chemistry: the case of verdazyl radicals

We present an efficient and accurate computational procedure to calculate properties measurable by EPR spectroscopy. We simulate a molecular dynamics (MD) trajectory by employing the quantum mechanically derived force field (QMDFF) [S. Grimme, J. Chem. Theory Comput., 2014, 10, 4497] and sample the trajectory at different time steps. For each snapshot EPR properties are calculated with a hybrid density functional theory (DFT) method. EPR spectra are simulated based on the averaged results. We applied the strategy to a number of previously published and novel verdazyl radicals, for which we recorded EPR spectra. The resulting simulated spectra are compatible with experiment already before employing an additional fitting step, in contrast to those from single point electronic-structure calculations. After the refinement, the experimental data are excellently reproduced, and the fitted EPR parameters do not deviate much from the calculated ones. This provides confidence in ascribing a direct physical meaning to the refined data in terms of experimental EPR parameters rather than merely considering them as mathematical fit parameters. We also find that couplings to hydrogen nuclei have a significant influence on the spectra of verdazyl radicals.

Rovibrational and Temperature Effects in Theoretical Studies of NMR Parameters

The demand for high precision calculations of NMR shieldings (or their related values, chemical shifts δ) and spin–spin coupling constants facilitating and supporting detailed interpretations of NMR spectra increases hand in hand with the development of computational techniques and hardware resources. Highly sophisticated calculations including even relativistic effects are nowadays possible for these properties. However, NMR parameters depend not only on molecular structure and environment but also on molecular flexibility and temperature and the apparent success of theoretical predictions for molecular equilibrium geometries creates a demand for zero-point vibrational and temperature corrections. In this chapter we describe briefly the theory behind rovibrational corrections and review then some important contributions to this field.

The electronic circular dichroism of nicotine in aqueous solution: a test case for continuum and mixed explicit-continuum solvation approaches

In this paper, we extend an integrated QM/MM/polarizable continuum model (PCM) method, which combines a fluctuating charge (FQ) approach to the MM polarization with the PCM, to describe electronic circular dichroism (ECD) spectra of systems in aqueous solution. The main features of the approach are presented, and then applications to the UV and ECD spectra of neutral (S)-nicotine in aqueous solution are reported. The performance of the QM/FQ/PCM is compared with that of the PCM against newly measured UV ECD spectra, which are in agreement with previous findings. The inclusion of specific solvation effects via the QM/FQ/PCM method leads to an improvement in the calculated spectra compared to the experimental findings, though the pure PCM results are still qualitatively correct and are a useful tool for the characterization of the states.

NMR J-coupling constants in cisplatin derivatives studied by molecular dynamics and relativistic DFT

Solvent effects on J((195)Pt-(15)N) one-bond nuclear spin-spin coupling constants (J(PtN)) of cisplatin [cis-diamminedichloroplatinum(II)] and three cisplatin derivatives are investigated using a combination of density functional theory (DFT) based ab initio molecular dynamics (aiMD) and all-electron relativistic DFT NMR calculations employing the two-component relativistic zeroth-order regular approximation (ZORA). Good agreement with experiment is obtained when explicit solvent molecules are considered and when the computations are performed with a hybrid functional. Spin-orbit coupling causes only small effects on J(PtN) . Key factors contributing to the magnitude of coupling constants are elucidated, with the most significant being the presence of solvent as well as the quality of the density functional and basis set combination. The solvent effects are of the same magnitude as J(PtN) calculated for gas-phase geometries. However, the trends of J(PtN) among the complexes are already present in the gas phase. Results obtained with a continuum solvent model agree quite well with the aiMD results, provided that the Pt solvent-accessible radius is carefully chosen. The aiMD results support the existence of a partial hydrogen-bond-like inverse-hydration-type interaction affording a weak (1)J(Pt⋅⋅⋅H(w)) coupling between the complexes and the coordinating water molecule.

Time-dependent density functional response theory for electronic chiroptical properties of chiral molecules

Methodology to calculate electronic chiroptical properties from time-dependent density functional theory (TDDFT) is outlined. Applications of TDDFT to computations of electronic circular dichroism, optical rotation, and optical rotatory dispersion are reviewed. Emphasis is put on publications from 2005 to 2010, but much of the older literature is also cited and discussed. The determination of the absolute configuration of chiral molecules by combined measurements and computations is an important application of TDDFT chiroptical methods and discussed in some detail. Raman optical activity (ROA) spectra are obtained from normal-mode derivatives of the optical rotation tensor and other linear response tensors. A few selected (ROA) benchmarks are reviewed.