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

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.

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.

Comparison of Length, Velocity, and Symmetric Gauges for the Calculation of Absorption and Electric Circular Dichroism Spectra with Real-Time Time-Dependent Density Functional Theory

A velocity and symmetric gauge implementation for real-time time-dependent density functional theory (RT-TDDFT) in the CP2K package using a Gaussian and plane wave approach is presented, including the explicit gauge-transformed contributions due to the nonlocal part of pseudopotentials. Absorption spectra of gas-phase α-pinene are calculated in length and velocity gauges in the long-wavelength approximation for the application of a δ pulse in linear and full order. The velocity gauge implementation is also applied to a solvated uracil molecule to showcase its use within periodic boundary conditions (PBC). For the calculation of the expectation value of the electric dipole moment in PBC, both the velocity representation and the modern theory of polarization give equivalent absorption spectra if a distributed reference point is used for the nonlocal term of the velocity operator. The discussion of linear response theory takes place in a unified framework in terms of linear response functions in propagator notation, distinguishing the parts of the linear response functions associated with perturbation and response. To further investigate gauge dependence, electric circular dichroism (ECD) spectra of α-pinene were calculated either as magnetic response to an electric field perturbation, in length or velocity gauge, or as electric response to a magnetic field perturbation in the symmetric gauge. Both approaches, electric and magnetic perturbations, have been found to yield equivalent ECD spectra.

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.

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.

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.