An official website of the United States government
Abstract 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.
more »
« less
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.
more »
« less
- Award ID(s):
- 1650942
- PAR ID:
- 10322350
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 15
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Jackson, George; Head-Gordon, Martin; Helgaker, Trygve; Liu, Wenjian; Osterwalder, Adreas (Ed.)We computed vibrational circular dichroism (VCD) and Raman optical activity (ROA) spectra for a test set of six chiral compounds using two standard density-functionals and an array of basis sets. We analysed the performance of property-oriented basis sets using a quadruple-zeta basis as a reference against four key metrics. We find little qualitative difference between the spectra produced by the larger basis sets (ORP, LPolX, aug-cc-pVTZ, and aug-cc-pVQZ), though their quantitative metrics exhibit wide variations. The smaller basis sets (rDPS, augD-3-21G, augT3-3-21G, Sadlej-pVTZ, and aug-cc-pVDZ) performed better for VCD rotatory strengths than for the corresponding ROA circular intensity differences (CIDs). However, this trend diminishes as the basis-set size is increased, lending validity to the conclusion that more robust property-oriented basis sets are required for ROA spectral generation than that of VCD. We observed improved performance in the mid-infrared region compared to the high-frequency regime, as well as overestimation of VCD rotatory strengths in the latter region as compared to the reference. We conclude that the ORP and LPol-ds basis sets are the most efficient and effective choices of basis set for the prediction of VCD and ROA spectra, as they provide both highly accurate results at reduced computational expense.more » « less
-
We present an ab initio method for computing vibro-polariton and phonon-polariton spectra of molecules and solids coupled to the photon modes of optical cavities. We demonstrate that if interactions of cavity photon modes with both nuclear and electronic degrees of freedom are treated on the level of the cavity Born–Oppenheimer approximation, spectra can be expressed in terms of the matter response to electric fields and nuclear displacements, which are readily available in standard density functional perturbation theory implementations. In this framework, results over a range of cavity parameters can be obtained without the need for additional electronic structure calculations, enabling efficient calculations on a wide range of parameters. Furthermore, this approach enables results to be more readily interpreted in terms of the more familiar cavity-independent molecular electric field response properties, such as polarizability and Born effective charges, which enter into the vibro-polariton calculation. Using corresponding electric field response properties of bulk insulating systems, we are also able to obtain the Γ point phonon-polariton spectra of two dimensional (2D) insulators. Results for a selection of cavity-coupled molecular and 2D crystal systems are presented to demonstrate the method.more » « less
-
Abstract This work presents the first steps to modeling synthetic rovibrational spectra for all molecules of astrophysical interest using a new approach implemented in the Prometheus code. The goal is to create a new comprehensive source of first-principles molecular spectra, thus bridging the gap for missing data to help drive future high-resolution studies. Our primary application domain is for molecules identified as signatures of life in planetary atmospheres (biosignatures), but our approach is general and can be applied to other systems. In this work we evaluate the accuracy of our method by studying four diatomic molecules, H 2 , O 2 , N 2 , and CO, all of which have well-known spectra. Prometheus uses the transition-optimised shifted Hermite (TOSH) theory to account for anharmonicity for the fundamental ν = 0 → ν = 1 band, along with thermal-profile modeling for the rotational transitions. To this end, we expand TOSH theory to enable the modeling of rotational constants. We show that our simple model achieves results that are a better approximation of the real spectra than those produced through an harmonic approach. We compare our results with high-resolution HITRAN and ExoMol spectral data. We find that modeling accuracy tends to diminish for rovibrational transition away from the band origin, thus highlighting the need for the theory to be further adapted.more » « less
-
A<sc>bstract</sc> In the decomposition of gauge-theory amplitudes into kinematic and color factors, the color factors (at a given loop orderL) span a proper subspace of the extended trace space (which consists of single and multiple traces of generators of the gauge group, graded by powers ofN). Using an iterative process, we systematically construct theL-loop color space of four-point amplitudes of fields in the adjoint representation of SU(N), SO(N), or Sp(N). We define the null space as the orthogonal complement of the color space. For SU(N), we confirm the existence of four independent null vectors (forL≥ 2) and for SO(N) and Sp(N), we establish the existence of seventeen independent null vectors (forL≥ 5). Each null vector corresponds to a group-theory constraint on the color-ordered amplitudes of the gauge theory.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0088922
