Cotton-Mouton Effect Using Algebraic Diagrammatic Construction Schemes in the Intermediate State Representation Formalism.

The Cotton-Mouton effect is theoretically investigated for a selected set of molecules by using a novel computational methodology based on algebraic diagrammatic construction (ADC) schemes in the intermediate state representation (ISR) formulation. Therefore, the electronic contributions to the frequency-dependent polarizabilities and, for the first time, to the magnetizabilities as well as mixed electric and magnetic hypermagnetizabilities have been computed in the ADC/ISR framework. In addition to calculation of the Cotton-Mouton constant and the birefringence, the gauge origin dependence of the computed tensors and the applied methodology are thoroughly investigated. The new ADC/ISR methodology, employing the recently presented responsefun package, is applied to a test set of Ne and small molecules (H2, HF, O2, CO2, and benzene) and compared to data from the experiment as well as other ab initio methods. The presented theoretical ab initio ADC/ISR approach is a substantial extension of the available computational methods for the investigation of complex nonlinear properties, however, with a gauge origin dependence inherent to the method that decreases with increasing perturbation order.

responsefun: Fun with Response Functions in the Algebraic Diagrammatic Construction Framework.

We present the open-source responsefun package, which implements a universally applicable procedure for computing molecular response properties within the algebraic diagrammatic construction (ADC) framework, exploiting the intermediate state representation (ISR) approach. With symbolic mathematics, the user can simply enter textbook sum-over-states (SOS) expressions from time-dependent perturbation theory, which are then automatically translated into the corresponding symbolic ADC/ISR formulations. Using the data structures provided by the hybrid Python/C++ module adcc for calculating excited states with ADC, the specified response property is directly evaluated, and the result is returned to the user. Employing the novel responsefun package, we present the first ADC/ISR calculations of second-order hyperpolarizability tensors and three-photon-absorption matrix elements.

Solving response expressions in the ADC/ISR framework.

We present an implementation for the calculation of molecular response properties using the algebraic-diagrammatic construction (ADC)/intermediate state representation approach. For the second-order ADC model [ADC(2)], a memory-efficient ansatz avoiding the storage of double excitation amplitudes is investigated. We compare the performance of different numerical algorithms for the solution of the underlying response equations for ADC(2) and show that our approach also strongly improves the convergence behavior for the investigated algorithms compared with the standard implementation. All routines are implemented in an open-source Python library.

Magnetic circular dichroism within the algebraic diagrammatic construction scheme of the polarization propagator up to third order.

We present an implementation of the B term of Magnetic Circular Dichroism (MCD) within the Algebraic Diagrammatic Construction (ADC) scheme of the polarization propagator and its Intermediate State Representation. As illustrative results, the MCD spectra of the ADC variants ADC(2), ADC(2)-x, and ADC(3) of the molecular systems uracil, 2-thiouracil, 4-thiouracil, purine, hypoxanthine 1,4-naphthoquinone, 9,10-anthraquinone, and 1-naphthylamine are computed and compared with results obtained by using the Resolution-of-Identity Coupled-Cluster Singles and Approximate Doubles method, with literature Time-Dependent Density Functional Theory results, and with available experimental data.

Exploring the accuracy and usefulness of semi-empirically scaled ADC schemes by blending second and third order terms.

Different approaches to mixed-order algebraic-diagrammatic construction (ADC) schemes are investigated. The performance of two different strategies for scaling third-order contributions to the ADC secular matrix is evaluated. Both considered schemes employ a single tuning parameter and conserve general properties inherent to all ADC methods, such as hermiticity and size-consistency. The first approach, scaled-matrix ADC[(2) + x(3)], scales all contributions first occurring in ADC(3) equally and leads to an improvement of the accuracy of excitation energies compared to ADC(3) for x = 0.4-0.5. However, with respect to excited state dipole moments, this method provides lower accuracy than ADC(3). The second scaling approach, MP[(1) + x(2)] - ISR(3), scales the second order contributions of the ground-state wavefunction and derives a rigorous ADC scheme via the intermediate state representation formalism. Although the error in excitation energies is not improved, this method provides insight into the relevance of the individual terms of the ADC(3) matrix and indicates that the MP(2) wavefunction is, indeed, the optimal reference wavefunction for deriving a third-order single-reference ADC scheme.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Analytical gradients for core-excited states in the algebraic diagrammatic construction (ADC) framework.

Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core-valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn-Teller distorted geometry of the 1s → π* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states.

Ab Initio Excited-State Electronic Circular Dichroism Spectra Exploiting the Third-Order Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator.

Excited-state rotatory strengths are reported for the first time at a correlated ab initio level, here with the algebraic diagrammatic construction scheme of the polarization propagator up to the third order. To demonstrate the capabilities of this computational approach, the gas phase S1 electronic circular dichroism spectra of the bicyclic ketones (1R)-camphor, (1R)-norcamphor, and (1R)-fenchone have been calculated at the ADC(3) level of theory. Furthermore, the solution excited-state spectra of the energetically lowest conformer of R-(+)-1,1'-bi(2-naphthol) have been computed with inclusion of a polarizable continuum model at the ADC(2) level of theory.

Unitary coupled-cluster approach for the calculation of core-excited states and x-ray absorption spectra.

In this work, we present the core-valence separation (CVS) approximation applied to unitary coupled-cluster (UCC) theory for the calculation of core-excited states and the simulation of x-ray absorption spectroscopy (XAS). Excitation energies and oscillator strengths of small- to medium-sized organic molecules have been computed using the second-order and extended second-order UCC schemes (CVS-UCC2 and CVS-UCC2-x) as well as the third-order scheme (CVS-UCC3). All results are compared to the corresponding algebraic-diagrammatic construction methods and experimental data. The agreement between CVS-UCC and experimental data demonstrates its potential as a new approach for the calculation of XAS.

Electronic circular dichroism spectra using the algebraic diagrammatic construction schemes of the polarization propagator up to third order.

Expressions for the calculation of rotatory strengths using the algebraic diagrammatic construction (ADC) scheme of the polarization propagator in both length and velocity gauges have been implemented. This enables the simulation of electronic circular dichroism (ECD) spectra at the ADC level up to third order of perturbation theory. The ADC(n) methods produce rotatory strengths of comparable accuracy to those obtained with coupled cluster methods of corresponding approximation levels, as evaluated for methyloxirane, methylthiirane, dimethyloxirane, dimethylthiirane, hydrogen peroxide, and dihydrogen disulfide. ECD spectra of (1R)-camphor, (1R)-norcamphor, and (1R)-fenchone computed at the third order ADC(3) level of theory are shown to agree very favorably with experimental gas phase spectra, demonstrating the usefulness of ADC for the calculation of chiro-optical properties of organic molecules. ADC(2) in combination with the polarizable continuum model is shown to successfully reproduce the ECD spectrum of the L-epinephrine enantiomer in water, further demonstrating the applicability of this approach.

Deciphering the Electronic Transitions of Thiophene-Based Donor-Acceptor-Donor Pentameric Ligands Utilized for Multimodal Fluorescence Microscopy of Protein Aggregates.

Anionic pentameric thiophene acetates can be used for fluorescence detection and diagnosis of protein amyloid aggregates. Replacing the central thiophene unit by benzothiadiazole (BTD) or quinoxaline (QX) leads to large emission shifts and basic spectral features have been reported [Chem. Eur. J. 2015, 21, 15133-13137]. Here we present new detailed experimental results of solvent effects, time-resolved fluorescence and examples employing multi-photon microscopy and lifetime imaging. Quantum chemical response calculations elucidate how the introduction of the BTD/QX groups changes the electronic states and emissions. The dramatic red-shift follows an increased conjugation and quinoid character of the π-electrons of the thiophene backbone. An efficient charge transfer in the excited states S1 and S2 compared to the all-thiophene analogue makes these more sensitive to the polarity and quenching by the solvent. Taken together, the results guide in the interpretation of images of stained Alzheimer disease brain sections employing advanced fluorescence microscopy and lifetime imaging, and can aid in optimizing future fluorescent ligand development.

Complex excited state polarizabilities in the ADC/ISR framework.

We present the derivation and implementation of complex, frequency-dependent polarizabilities for excited states using the algebraic-diagrammatic construction for the polarization propagator (ADC) and its intermediate state representation. Based on the complex polarizability, we evaluate C6 dispersion coefficients for excited states. The methodology is implemented up to third order in perturbation theory in the Python-driven adcc toolkit for the development and application of ADC methods. We exemplify the approach using illustrative model systems and compare it to results from other ab initio methods and from experiments.

Third-Order Unitary Coupled Cluster (UCC3) for Excited Electronic States: Efficient Implementation and Benchmarking.

The efficient implementation of the third-order unitary coupled-cluster scheme (UCC3) for the calculation of excited electronic states is reported. The UCC3 scheme and its second-order UCC2 variant have been benchmarked and compared to Jacquemin's recently introduced, as well as Thiel's well-established, benchmark sets for excitation energies and oscillator strengths. For the latter, the calculation of 134 excited singlet and 71 excited triplet states of 28 small- to medium-sized organic molecules has revealed that UCC2 exhibits a mean error and standard deviation of 0.36 ± 0.41 eV for singlet states and 0.22 ± 0.21 eV for triplet states, whereas UCC3 revealed an accuracy of 0.06 ± 0.27 eV for singlet and -0.22 ± 0.15 eV for triplet states. In addition, the oscillator strengths obtained with effective transition moments correct through second order in perturbation theory are in very good agreement with literature data.

Hermitian second-order methods for excited electronic states: Unitary coupled cluster in comparison with algebraic-diagrammatic construction schemes.

Employing an intermediate state representation (ISR) approach, Hermitian second-order methods for the calculation of electronic excitation energies are presented and compared in detail. These comprise the algebraic-diagrammatic construction scheme for the polarization propagator, a hybrid second-order ISR scheme based on traditional coupled-cluster theory as well as two similar approaches based on a unitary coupled-cluster (UCC) ansatz. Although in a strict perturbation-theoretical framework all prove to be identical, differences emerge when the corresponding converged cluster amplitudes are used and depending on how the similarity-transformed UCC Hamiltonian is evaluated. The resulting excitation energies, however, do not significantly differ for systems well described by means of perturbation theory.

Nuclear dynamics in resonant inelastic X-ray scattering and X-ray absorption of methanol.

We report on a combined theoretical and experimental study of core-excitation spectra of gas and liquid phase methanol as obtained with the use of X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS). The electronic transitions are studied with computational methods that include strict and extended second-order algebraic diagrammatic construction [ADC(2) and ADC(2)-x], restricted active space second-order perturbation theory, and time-dependent density functional theory-providing a complete assignment of the near oxygen K-edge XAS. We show that multimode nuclear dynamics is of crucial importance for explaining the available experimental XAS and RIXS spectra. The multimode nuclear motion was considered in a recently developed "mixed representation" where dissociative states and highly excited vibrational modes are accurately treated with a time-dependent wave packet technique, while the remaining active vibrational modes are described using Franck-Condon amplitudes. Particular attention is paid to the polarization dependence of RIXS and the effects of the isotopic substitution on the RIXS profile in the case of dissociative core-excited states. Our approach predicts the splitting of the 2a″ RIXS peak to be due to an interplay between molecular and pseudo-atomic features arising in the course of transitions between dissociative core- and valence-excited states. The dynamical nature of the splitting of the 2a″ peak in RIXS of liquid methanol near pre-edge core excitation is shown. The theoretical results are in good agreement with our liquid phase measurements and gas phase experimental data available from the literature.

Similarities and differences of the Lagrange formalism and the intermediate state representation in the treatment of molecular properties.

When dealing with approximate wave functions, molecular properties can be computed either as expectation values or as derivatives of the energy with respect to a corresponding perturbation. In this work, the intermediate state representation (ISR) formalism for the computation of expectation values is compared to the Lagrange formalism following a derivative ansatz, which are two alternative approaches of which neither one can be considered superior in general. Within the ISR formalism, terms are included up to a given order of perturbation theory only, while in the Lagrange formalism, all terms are accounted for arising through the differentiation. Similarities and differences of the Lagrange and ISR formalism are illustrated using explicit working equations for selected methods and analyzing numerical results for a range of coupled-cluster as well as algebraic-diagrammatic construction (ADC) methods for excited states. The analysis explains why the ADC(3/2) method is able to yield a large amount of the orbital-relaxation effects for p-h states in contrast to ADC(2) although the same second-order ISR is used to represent the corresponding operator.

Analytic nuclear gradients of the algebraic-diagrammatic construction scheme for the polarization propagator up to third order of perturbation theory.

Analytic gradient expressions for the algebraic diagrammatic construction (ADC) scheme for the polarization propagator up to third order are derived using a Lagrangian approach. An implementation within the Q-CHEM electronic structure package for excited-state nuclear gradients of the ADC(2), ADC(2)-x, and ADC(3) models based on restricted and unrestricted Hartree-Fock references is presented. Details of the implementation and the applicability of the newly derived gradients for geometry optimizations and the quality of the resulting structures are discussed.

Algebraic-diagrammatic construction scheme for the polarization propagator including ground-state coupled-cluster amplitudes. I. Excitation energies.

An ad hoc modification of the algebraic-diagrammatic construction (ADC) scheme for the polarization propagator is presented. Within this approach, all first-order Møller-Plesset correlation coefficients occurring in the second-order ADC secular matrix are replaced by amplitudes obtained from a coupled cluster doubles (CCD) calculation. This new hybrid method, denoted CCD-ADC(2), has been tested on a series of small diatomic and triatomic molecules and benchmarked with respect to Thiel's benchmark set of medium-sized organic molecules. For the latter, the calculation of 134 singlet and 71 triplet states has shown that CCD-ADC(2) exhibits a mean error and standard deviation of 0.15 ± 0.34 eV for singlet states and 0.0 ± 0.17 eV for triplet states with respect to the provided theoretical best estimates, whereas standard ADC(2) has a mean error and standard deviation of 0.22 ± 0.30 eV for singlet and 0.12 ± 0.16 eV for triplet states. The corresponding extended second-order schemes ADC(2)-x and CCD-ADC(2)-x revealed accuracies of -0.70 ± 0.32 eV and -0.76 ± 0.33 eV for singlet states and -0.55 ± 0.20 eV and -0.67 ± 0.22 eV for triplet states, respectively. Furthermore, the investigation of excited-state potential energy curves along the dissociation of the N2 molecule has shown that the higher reliability of the ground-state CCD method as compared to MP2 is also inherent to the excited states. While the curves obtained at the ADC(2) level break down at around 2 Å, the ones obtained at CCD-ADC(2) remain reasonable up to about 3.5 Å.

Algebraic-diagrammatic construction scheme for the polarization propagator including ground-state coupled-cluster amplitudes. II. Static polarizabilities.

The modification of the algebraic-diagrammatic construction (ADC) scheme for the polarization propagator using ground-state coupled-cluster (CC) instead of Møller-Plesset (MP) amplitudes, referred to as CC-ADC, is extended to the calculation of molecular properties, in particular, dipole polarizabilities. Furthermore, in addition to CC with double excitations (CCD), CC with single and double excitations (CCSD) amplitudes can be used, also in the second-order transition moments of the ADC(3/2) method. In the second-order CC-ADC(2) variants, the MP correlation coefficients occurring in ADC are replaced by either CCD or CCSD amplitudes, while in the F/CC-ADC(2) and F/CC-ADC(3/2) variants, they are replaced only in the second-order modified transition moments. These newly implemented variants are used to calculate the static dipole polarizability of several small- to medium-sized molecules, and the results are compared to the ones obtained by full configuration interaction or experiment. It is shown that the results are consistently improved by the use of CC amplitudes, in particular, for aromatic systems such as benzene or pyridine, which have proven to be difficult cases for standard ADC approaches. In this case, the second-order CC-ADC(2) and F/CC-ADC(2) variants yield significantly better results than the standard third-order ADC(3/2) method, at a computational cost amounting to only about 1% of the latter.

Characterizing Bonding Patterns in Diradicals and Triradicals by Density-Based Wave Function Analysis: A Uniform Approach.

Density-based wave function analysis enables unambiguous comparisons of the electronic structure computed by different methods and removes ambiguity of orbital choices. We use this tool to investigate the performance of different spin-flip methods for several prototypical diradicals and triradicals. In contrast to previous calibration studies that focused on energy gaps between high- and low spin-states, we focus on the properties of the underlying wave functions, such as the number of effectively unpaired electrons. Comparison of different density functional and wave function theory results provides insight into the performance of the different methods when applied to strongly correlated systems such as polyradicals. We show that canonical molecular orbitals for species like large copper-containing diradicals fail to correctly represent the underlying electronic structure due to highly non-Koopmans character, while density-based analysis of the same wave function delivers a clear picture of the bonding pattern.