Towards an Improved Method for the Prediction of Linear Response Properties of Small Organic Molecules

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.

Name thesis:degree_name

Doctor of Philosophy

Level thesis:degree_level

doctoral

Discipline thesis:degree_discipline

Chemistry

Department dc:contributor.department

Chemistry

Grantor dc:publisher

Virginia Tech

Year dc:date.issued

2021

Author dc:creator

• Dcunha, Ruhee Lancelot

Chair dc:contributor.committeechair

• Crawford, T. Daniel

Committee members dc:contributor.committeemember

• Mayhall, Nicholas
• Troya, Diego
• Embree, Mark P.

• molecular properties
• coupled cluster theory
• reduced scaling

Statement dc:rights

• In Copyright

Licence dc:rights.uri

• http://rightsstatements.org/vocab/InC/1.0/

Dc Identifier Other

vt_gsexam:32231

OAI identifier oai:identifier

oai:vtechworks.lib.vt.edu:10919/104677