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

First Principles Simulations of Optical Rotation of Chiral Molecular Crystals

In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT-PBC). Calculations are compared with experimental measurements and show semi-quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X-ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.

A Perspective on the Simulation of Electronic Circular Dichroism and Circularly Polarized Luminescence Spectra in Chiral Solid Materials

Chiral materials have shown tremendous potential for many technological applications, such as optoelectronics, sensing, magnetism, information technology, and imaging. Characterization of these materials is mostly based on chiroptical spectroscopies, such as electronic circular dichroism (ECD) and circularly polarized luminescence (CPL). These experimental measurements would greatly benefit from theoretical simulations for interpretation of the spectra as well as predictions on new materials. While ECD and CPL simulations are well established for molecular systems, they are not for materials. In this Perspective, we describe the theoretical quantities necessary to simulate ECD and CPL spectra in oriented systems. Then, we discuss the approximate strategies currently used to perform these calculations, what computational machinery is already available to develop more general approaches, and some of the open challenges for the simulation of ECD and CPL spectra in solid materials. When methods that are as reliable and computationally efficient as those for molecules are developed, these simulations will provide invaluable insight and guidance for the rational design of optically active materials.

Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-Accelerated Computer Architectures

We present the development and implementation of relativistic coupled cluster linear response theory (CC-LR), which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin for the magnetic properties) as well as taking into account the finite lifetime of excited states in the framework of damped response theory. We showcase our implementation, which is capable to offload the computationally intensive tensor contractions characteristic of coupled cluster theory onto graphical processing units, in the calculation of (a) frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a particular emphasis on the calculation of valence absorption cross sections for the I2 molecule; (b) indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H2Se-H2O dimer as a prototypical system containing hydrogen bonds; and (c) optical rotations at the sodium D line for hydrogen peroxide analogues (H2Y2, Y = O, S, Se, Te). Thanks to this implementation, we are able to show the similarities in performance, but often the significant discrepancies, between CC-LR and approximate methods such as density functional theory. Comparing standard CC response theory with the flavor based upon the equation of motion formalism, we find that for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region, such as spin-spin couplings, on the other hand, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table, it may become increasingly difficult to measure pure optical rotation at the sodium D line due to the appearance of absorbing states.

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.

A computational approach for modeling electronic circular dichroism of solvated chromophores

The present study consists in a novel computational protocol to model the UV-circular dichroism spectra of solvated species. It makes use of quantum-chemical calculations on a series of conformations of a flexible chromophore or on a series of chromophore/solvent clusters extracted from molecular dynamic simulations. The protocol is described and applied to the aqueous cationic tripeptide GAG⁺ and to the aqueous neutral decapeptide (GVGVP)₂ . The protocol has proven able to: (i) properly consider the conformational motion of solute in the given environment; (ii) give the actual statistical weight of each conformational state; (iii) provide a reliable quantum mechanical method able to reproduce the spectral features. Temperature effects on conformations and spectral properties are properly taken into account. The role of explicit solvent on the conformational analysis and the spectra calculation is discussed. The comparison of the calculated circular dichroism spectra with experimental ones recorded at different temperatures represents a strict validation test of the method.

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.

Probing Molecular Chirality of Ground and Electronically Excited States in the UV-vis and X-ray Regimes: An EOM-CCSD Study

We present several strategies for computing electronic circular dichroism (CD) spectra across different frequency ranges at the equation-of-motion coupled-cluster singles and doubles level of theory. CD spectra of both ground and electronically excited states are discussed. For selected cases, the approach is compared with coupled-cluster linear response results as well as time-dependent density functional theory. The extension of the theory to include the effect of spin-orbit coupling is presented and illustrated by calculations of X-ray CD spectra at the L-edge.

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

We present an origin-invariant approach to compute the full optical rotation tensor (Buckingham/Dunn tensor) in the length dipole gauge without recourse to London atomic orbitals, called LG(OI). The LG(OI) approach is simpler and less computationally demanding than the more common length gauge (LG)-London and modified velocity gauge (MVG) approaches, and it can be used with any approximate wave function or density functional method. We report an implementation at the coupled cluster with single and double excitations level (CCSD), for which we present the first simulations of the origin-invariant Buckingham/Dunn tensor in the LG. We compare LG(OI) and MVG results on a series of 22 organic molecules, showing good linear correlation between the approaches, although for small tensor elements, they provide values of opposite sign. We also attempt to decouple the effects of electron correlation and basis set incompleteness on the choice of gauge for specific rotation calculations on simple test systems. The simulations show a smooth convergence of the LG(OI) and MVG results with the basis set size toward the complete basis set limit. However, these preliminary results indicate that CCSD may not be close to a complete description of the electron correlation effects on this property even for small molecules and that basis set incompleteness may be a less important cause of discrepancy between choices of gauge than electron correlation incompleteness.

Electronic Currents Induced by Optical Fields and Rotatory Power Density in Chiral Molecules

The electric dipole–magnetic dipole polarizability tensor κ′, introduced to interpret the optical activity of chiral molecules, has been expressed in terms of a series of density functions kαβ′, which can be integrated all over the three-dimensional space to evaluate components καβ′ and trace καα′. A computational approach to kαβ′, based on frequency-dependent electronic current densities induced by monochromatic light shining on a probe molecule, has been developed. The dependence of kαβ′ on the origin of the coordinate system has been investigated in connection with the corresponding change of καβ′. It is shown that only the trace kαα′ of the density function defined via dynamic current density evaluated using the continuous translation of the origin of the coordinate system is invariant of the origin. Accordingly, this function is recommended as a tool that is quite useful for determining the molecular domains that determine optical activity to a major extent. A series of computations on the hydrogen peroxide molecule, for a number of different HO–OH dihedral angles, is shown to provide a pictorial documentation of the proposed method.

First-Principles Calculation of the Optical Rotatory Power of Periodic Systems: Application on α-Quartz, Tartaric Acid Crystal, and Chiral (n,m)-Carbon Nanotubes

The self-consistent coupled-perturbed (SC-CP) method in the CRYSTAL program has been adapted to obtain electromagnetic optical rotation properties of chiral periodic systems based on the calculation of the magnetic moment induced by the electric field. Toward that end, an expression for the magnetic transition moment is developed, which involves an appropriate electronic angular momentum operator. This operator is forced to be hermitian so that the chiroptical properties are real. In our formulation, the trace of the optical rotatory power matrix is gauge-origin-invariant as long as the electric dipole transition matrix elements are obtained using the velocity (rather than position) operator. On the other hand, the component along the optic axis is invariant in general for uniaxial and biaxial crystals. Under the same conditions, these properties also do not depend on the so-called missing integers that occur in the treatment of the electric dipole moment of quasi-one-dimensional periodic systems or the analogue of missing integers for the case of higher dimensionality. Tests on a model H₂O₂ polymer confirm the formalism and, as desired, show that the calculated properties are independent of the size and definition of the unit cell. In addition, an empirical relation to a finite oligomer gauge-including atomic orbital (GIAO) calculation is found. Applications, with comparison to experiment, are carried for α-quartz, tartaric acid crystal, and carbon nanotubes. Future developments of this initial approach to chiroptical properties in the solid state are noted.