Magnetic circular dichroism spectra from resonant and damped coupled cluster response theory

A computational expression for the Faraday A term of magnetic circular dichroism (MCD) is derived within coupled cluster response theory and alternative computational expressions for the B term are discussed. Moreover, an approach to compute the (temperature-independent) MCD ellipticity in the context of coupled cluster damped response is presented, and its equivalence with the stick-spectrum approach in the limit of infinite lifetimes is demonstrated. The damped response approach has advantages for molecular systems or spectral ranges with a high density of states. Illustrative results are reported at the coupled cluster singles and doubles level and compared to time-dependent density functional theory results.Comment: Submitted to J. Chem. Phys. on May 10, 202

Atomic-scale representation and statistical learning of tensorial properties

This chapter discusses the importance of incorporating three-dimensional symmetries in the context of statistical learning models geared towards the interpolation of the tensorial properties of atomic-scale structures. We focus on Gaussian process regression, and in particular on the construction of structural representations, and the associated kernel functions, that are endowed with the geometric covariance properties compatible with those of the learning targets. We summarize the general formulation of such a symmetry-adapted Gaussian process regression model, and how it can be implemented based on a scheme that generalizes the popular smooth overlap of atomic positions representation. We give examples of the performance of this framework when learning the polarizability and the ground-state electron density of a molecule

Review of biorthogonal coupled cluster representations for electronic excitation

Single reference coupled-cluster (CC) methods for electronic excitation are based on a biorthogonal representation (bCC) of the (shifted) Hamiltonian in terms of excited CC states, also referred to as correlated excited (CE) states, and an associated set of states biorthogonal to the CE states, the latter being essentially configuration interaction (CI) configurations. The bCC representation generates a non-hermitian secular matrix, the eigenvalues representing excitation energies, while the corresponding spectral intensities are to be derived from both the left and right eigenvectors. Using the perspective of the bCC representation, a systematic and comprehensive analysis of the excited-state CC methods is given, extending and generalizing previous such studies. Here, the essential topics are the truncation error characteristics and the separability properties, the latter being crucial for designing size-consistent approximation schemes. Based on the general order relations for the bCC secular matrix and the (left and right) eigenvector matrices, formulas for the perturbation-theoretical (PT) order of the truncation errors (TEO) are derived for energies, transition moments, and property matrix elements of arbitrary excitation classes and truncation levels. In the analysis of the separability properties of the transition moments, the decisive role of the so-called dual ground state is revealed. Due to the use of CE states the bCC approach can be compared to so-called intermediate state representation (ISR) methods based exclusively on suitably orthonormalized CE states. As the present analysis shows, the bCC approach has decisive advantages over the conventional CI treatment, but also distinctly weaker TEO and separability properties in comparison with a full (and hermitian) ISR method

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted

Optoelectronic and Excitonic Properties of Oligoacenes: Substantial Improvements from Range-Separated Time-Dependent Density Functional Theory

The optoelectronic and excitonic properties in a series of linear acenes (naphthalene up to heptacene) are investigated using range-separated methods within time-dependent density functional theory (TDDFT). In these rather simple systems, it is well-known that TDDFT methods using conventional hybrid functionals surprisingly fail in describing the low-lying La and Lb valence states, resulting in large, growing errors for the La state and an incorrect energetic ordering as a function of molecular size. In this work, we demonstrate that the range-separated formalism largely eliminates both of these errors and also provides a consistent description of excitonic properties in these systems. We further demonstrate that re-optimizing the percentage of Hartree-Fock exchange in conventional hybrids to match wavefunction-based benchmark calculations still yields serious errors, and a full 100% Hartree-Fock range separation is essential for simultaneously describing both of the La and Lb transitions. Based on an analysis of electron-hole transition density matrices, we finally show that conventional hybrid functionals overdelocalize excitons and underestimate quasiparticle energy gaps in the acene systems. The results of our present study emphasize the importance of both a range-separated and asymptotically-correct contribution of exchange in TDDFT for investigating optoelectronic and excitonic properties, even for these simple valence excitations.Comment: Accepted by the Journal of Chemical Theory and Computatio