Quantum chemical methods to predict experimental chiroptical properties by solving the time-dependent Schrödinger equation are useful in the assignment of absolute configurations. Chiroptical properties, being very sensitive to the electronic structure of the system, require highly-accurate methods on the one hand and on the other, need to be able to be computed with limited computational resources.

The calculation of the optical rotation in the solution phase is complicated by solvent effects. In order to capture those solvent effects, we present a study that uses conformational averaging and time-dependent density functional theory calculations that incorporate solvent molecules explicitly in the quantum mechanical region. While considering several controllable parameters along which the system's optical rotation varies, we find that the sampling of the dynamical trajectory and the density functional chosen have the largest impact on the value of the rotation.

In order to eliminate the arbitrariness of the choice of density functional, we would prefer to use coupled cluster theory, a robust and systematically improvable method. However, the high-order polynomial scaling of coupled cluster theory makes it intractable for numerous large calculations, including the conformational averaging required for optical rotation calculations in solution. We therefore attempt to reduce the scaling of a linear response coupled cluster singles and doubles (LR-CCSD) calculation via a perturbed pair natural orbital (PNO++) local correlation approach which uses an orbital space created using a perturbed density matrix. We find that by creating a "combined PNO++" space, incorporating a set of orbitals from the unperturbed pair natural orbital (PNO) space into the PNO++ space, we can obtain well-behaved convergence behavior for both CCSD correlation energies and linear response properties, including dynamic polarizabilities and optical rotations, for the small systems considered.

The PNO++ and combined PNO++ methods require aggressive truncation to keep the computational cost low, due to an expensive two-electron integral transformation at the beginning of the calculation. We apply the methods to larger systems than previously studied and refine them for more aggressive truncation by exploring an alternative form of the perturbed density and a perturbation-including weak pair approximation.